1
|
Reiter S, Gordiy I, Kollmannsberger KL, Liu F, Thyrhaug E, Leister D, Warnan J, Hauer J, de Vivie-Riedle R. Molecular interactions of photosystem I and ZIF-8 in bio-nanohybrid materials. Phys Chem Chem Phys 2024; 26:23228-23239. [PMID: 39192757 DOI: 10.1039/d4cp03021d] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/29/2024]
Abstract
Bio-nanohybrid devices featuring natural photocatalysts bound to a nanostructure hold great promise in the search for sustainable energy conversion. One of the major challenges of integrating biological systems is protecting them against harsh environmental conditions while retaining, or ideally enhancing their photophysical properties. In this mainly computational work we investigate an assembly of cyanobacterial photosystem I (PS I) embedded in a metal-organic framework (MOF), namely the zeolitic imidazolate framework ZIF-8. This complex has been reported experimentally [Bennett et al., Nanoscale Adv., 2019, 1, 94] but so far the molecular interactions between PS I and the MOF remained elusive. We show via absorption spectroscopy that PS I remains intact throughout the encapsulation-release cycle. Molecular dynamics (MD) simulations further confirm that the encapsulation has no noticeable structural impact on the photosystem. However, the MOF building blocks frequently coordinate to the Mg2+ ions of chlorophylls in the periphery of the antenna complex. High-level quantum mechanical calculations reveal charge-transfer interactions, which affect the excitonic network and thereby may reversibly change the fluorescence properties of PS I. Nevertheless, our results highlight the stability of PS I in the MOF, as the reaction center remains unimpeded by the heterogeneous environment, paving the way for applications in the foreseeable future.
Collapse
Affiliation(s)
- Sebastian Reiter
- Department of Chemistry, Ludwig-Maximilians-Universität München, Butenandtstr. 11, 81377 Munich, Germany.
| | - Igor Gordiy
- Department of Chemistry, Ludwig-Maximilians-Universität München, Butenandtstr. 11, 81377 Munich, Germany.
| | - Kathrin L Kollmannsberger
- Chair of Inorganic and Metal-Organic Chemistry, Department of Chemistry and Catalysis Research Center (CRC), TUM School of Natural Sciences, Technical University of Munich, Lichtenbergstr. 4, 85748 Garching, Germany
| | - Feng Liu
- Faculty of Biology, Ludwig-Maximilians-Universität München, Großhaderner Str. 2-4, 82152 Planegg-Martinsried, Germany
| | - Erling Thyrhaug
- Professorship of Dynamic Spectroscopy, Department of Chemistry and Catalysis Research Center (CRC), TUM School of Natural Sciences, Technical University of Munich, Lichtenbergstr. 4, 85748 Garching, Germany.
| | - Dario Leister
- Faculty of Biology, Ludwig-Maximilians-Universität München, Großhaderner Str. 2-4, 82152 Planegg-Martinsried, Germany
| | - Julien Warnan
- Chair of Inorganic and Metal-Organic Chemistry, Department of Chemistry and Catalysis Research Center (CRC), TUM School of Natural Sciences, Technical University of Munich, Lichtenbergstr. 4, 85748 Garching, Germany
| | - Jürgen Hauer
- Professorship of Dynamic Spectroscopy, Department of Chemistry and Catalysis Research Center (CRC), TUM School of Natural Sciences, Technical University of Munich, Lichtenbergstr. 4, 85748 Garching, Germany.
| | - Regina de Vivie-Riedle
- Department of Chemistry, Ludwig-Maximilians-Universität München, Butenandtstr. 11, 81377 Munich, Germany.
| |
Collapse
|
2
|
Sheng Z, Jiang T, Li W, Shuai Z. TD-DMRG Study of Exciton Dynamics with both Thermal and Static Disorders for Fenna-Matthews-Olson Complex. J Chem Theory Comput 2024. [PMID: 39087905 DOI: 10.1021/acs.jctc.4c00493] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/02/2024]
Abstract
Photosynthesis is a fundamental process that converts solar energy into chemical energy. Understanding the microscopic mechanisms of energy transfer in photosynthetic systems is crucial for the development of novel optoelectronic materials. Simulating these processes poses significant challenges due to the intricate interactions between electrons and phonons, compounded by static disorder. In this work, we present a numerically nearly exact study using the time-dependent density matrix renormalization group (TD-DMRG) method to simulate the quantum dynamics of the Fenna-Matthews-Olson (FMO) complex considering an eight-site model with both thermal and static disorders. We employ the thermo-field dynamics formalism for temperature effects. We merge all electronic interactions into one large matrix product state (MPS) site, boosting accuracy efficiently without increasing complexity. Previous combined experimental and computational studies indicated that the static disorders range from 30 to 90 cm-1 for different FMO sites. We employ a Gaussian distribution and the auxiliary bosonic operator approach to consider the static disorder in our TD-DMRG algorithm. We investigate the impact of different initial excitation sites, temperatures, and degrees of static disorder on the exciton dynamics and temporal coherence. It is found that under the influence of the experimentally determined static disorder strength, the exciton population evolution shows a non-negligible difference at zero temperature, while it is hardly affected at room temperature.
Collapse
Affiliation(s)
- Zirui Sheng
- School of Science and Engineering, The Chinese University of Hong Kong, Shenzhen 518172, P. R. China
| | - Tong Jiang
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, 100084 Beijing, P. R. China
| | - Weitang Li
- School of Science and Engineering, The Chinese University of Hong Kong, Shenzhen 518172, P. R. China
| | - Zhigang Shuai
- School of Science and Engineering, The Chinese University of Hong Kong, Shenzhen 518172, P. R. China
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, 100084 Beijing, P. R. China
| |
Collapse
|
3
|
Khanna A, Shedge SV, Zuehlsdorff TJ, Isborn CM. Calculating absorption and fluorescence spectra for chromophores in solution with ensemble Franck-Condon methods. J Chem Phys 2024; 161:044121. [PMID: 39077907 DOI: 10.1063/5.0217080] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/02/2024] [Accepted: 07/03/2024] [Indexed: 07/31/2024] Open
Abstract
Accurately modeling absorption and fluorescence spectra for molecules in solution poses a challenge due to the need to incorporate both vibronic and environmental effects, as well as the necessity of accurate excited state electronic structure calculations. Nuclear ensemble approaches capture explicit environmental effects, Franck-Condon methods capture vibronic effects, and recently introduced ensemble-Franck-Condon approaches combine the advantages of both methods. In this study, we present and analyze simulated absorption and fluorescence spectra generated with combined ensemble-Franck-Condon approaches for three chromophore-solvent systems and compare them to standard ensemble and Franck-Condon spectra, as well as to the experiment. Employing configurations obtained from ground and excited state ab initio molecular dynamics, three combined ensemble-Franck-Condon approaches are directly compared to each other to assess the accuracy and relative computational time. We find that the approach employing an average finite-temperature Franck-Condon line shape generates spectra nearly identical to the direct summation of an ensemble of Franck-Condon spectra at one-fourth of the computational cost. We analyze how the spectral simulation method, as well as the level of electronic structure theory, affects spectral line shapes and associated Stokes shifts for 7-nitrobenz-2-oxa-1,3-diazol-4-yl and Nile red in dimethyl sulfoxide and 7-methoxy coumarin-4-acetic acid in methanol. For the first time, our studies show the capability of combined ensemble-Franck-Condon methods for both absorption and fluorescence spectroscopy and provide a powerful tool for simulating linear optical spectra.
Collapse
Affiliation(s)
- Ajay Khanna
- Chemistry and Biochemistry, University of California Merced, Merced, California 95343, USA
| | - Sapana V Shedge
- Chemistry and Biochemistry, University of California Merced, Merced, California 95343, USA
| | - Tim J Zuehlsdorff
- Department of Chemistry, Oregon State University, Corvallis, Oregon 97331, USA
| | - Christine M Isborn
- Chemistry and Biochemistry, University of California Merced, Merced, California 95343, USA
| |
Collapse
|
4
|
Myers CA, Lu SY, Shedge S, Pyuskulyan A, Donahoe K, Khanna A, Shi L, Isborn CM. Axial H-Bonding Solvent Controls Inhomogeneous Spectral Broadening, While Peripheral H-Bonding Solvent Controls Vibronic Broadening: Cresyl Violet in Methanol. J Phys Chem B 2024; 128:5685-5699. [PMID: 38832562 DOI: 10.1021/acs.jpcb.4c01401] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/05/2024]
Abstract
The dynamics of the nuclei of both a chromophore and its condensed-phase environment control many spectral features, including the vibronic and inhomogeneous broadening present in spectral line shapes. For the cresyl violet chromophore in methanol, we here analyze and isolate the effect of specific chromophore-solvent interactions on simulated spectral densities, reorganization energies, and linear absorption spectra. Employing both chromophore and its condensed-phase environment control many spectral features, including the vibronic and inhomogeneous broadening present in spectral line shapes. For the cresyl violet chromophore in methanol, we here analyze and isolate the effect of specific chromophore-solvent interactions on simulated spectral densities, reorganization energies, and linear absorption spectra. Employing both force field and ab initio molecular dynamics trajectories along with the inclusion of only certain solvent molecules in the excited-state calculations, we determine that the methanol molecules axial to the chromophore are responsible for the majority of inhomogeneous broadening, with a single methanol molecule that forms an axial hydrogen bond dominating the response. The strong peripheral hydrogen bonds do not contribute to spectral broadening, as they are very stable throughout the dynamics and do not lead to increased energy-gap fluctuations. We also find that treating the strong peripheral hydrogen bonds as molecular mechanical point charges during the molecular dynamics simulation underestimates the vibronic coupling. Including these peripheral hydrogen bonding methanol molecules in the quantum-mechanical region in a geometry optimization increases the vibronic coupling, suggesting that a more advanced treatment of these strongly interacting solvent molecules during the molecular dynamics trajectory may be necessary to capture the full vibronic spectral broadening.
Collapse
Affiliation(s)
- Christopher A Myers
- Department of Chemistry and Biochemistry, University of California Merced, Merced, California 95343, United States
| | - Shao-Yu Lu
- Department of Chemistry and Biochemistry, University of California Merced, Merced, California 95343, United States
| | - Sapana Shedge
- Department of Chemistry and Biochemistry, University of California Merced, Merced, California 95343, United States
| | - Arthur Pyuskulyan
- Department of Chemistry and Biochemistry, University of California Merced, Merced, California 95343, United States
| | - Katherine Donahoe
- Department of Chemistry and Biochemistry, University of California Merced, Merced, California 95343, United States
| | - Ajay Khanna
- Department of Chemistry and Biochemistry, University of California Merced, Merced, California 95343, United States
| | - Liang Shi
- Department of Chemistry and Biochemistry, University of California Merced, Merced, California 95343, United States
| | - Christine M Isborn
- Department of Chemistry and Biochemistry, University of California Merced, Merced, California 95343, United States
| |
Collapse
|
5
|
Lyratzakis A, Daskalakis V, Xie H, Tsiotis G. The synergy between the PscC subunits for electron transfer to the P 840 special pair in Chlorobaculum tepidum. PHOTOSYNTHESIS RESEARCH 2024; 160:87-96. [PMID: 38625595 PMCID: PMC11108878 DOI: 10.1007/s11120-024-01093-7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/13/2023] [Accepted: 03/08/2024] [Indexed: 04/17/2024]
Abstract
The primary photochemical reaction of photosynthesis in green sulfur bacteria occurs in the homodimer PscA core proteins by a special chlorophyll pair. The light induced excited state of the special pair producing P840+ is rapidly reduced by electron transfer from one of the two PscC subunits. Molecular dynamics (MD) simulations are combined with bioinformatic tools herein to provide structural and dynamic insight into the complex between the two PscA core proteins and the two PscC subunits. The microscopic dynamic model involves extensive sampling at atomic resolution and at a cumulative time-scale of 22µs and reveals well defined protein-protein interactions. The membrane complex is composed of the two PscA and the two PscC subunits and macroscopic connections are revealed within a putative electron transfer pathway from the PscC subunit to the special pair P840 located within the PscA subunits. Our results provide a structural basis for understanding the electron transport to the homodimer RC of the green sulfur bacteria. The MD based approach can provide the basis to further probe the PscA-PscC complex dynamics and observe electron transfer therein at the quantum level. Furthermore, the transmembrane helices of the different PscC subunits exert distinct dynamics in the complex.
Collapse
Affiliation(s)
- Alexandros Lyratzakis
- Department of Chemistry, School of Science and Engineering, University of Crete, Heraklion, 70013, Greece
| | - Vangelis Daskalakis
- Department of Chemical Engineering, School of Engineering, University of Patras, Rion, Patras, 26504, Greece
| | - Hao Xie
- Max Planck Institute of Biophysics, 60438, Frankfurt am Main, Germany
| | - Georgios Tsiotis
- Department of Chemistry, School of Science and Engineering, University of Crete, Heraklion, 70013, Greece.
| |
Collapse
|
6
|
Sarngadharan P, Holtkamp Y, Kleinekathöfer U. Protein Effects on the Excitation Energies and Exciton Dynamics of the CP24 Antenna Complex. J Phys Chem B 2024; 128:5201-5217. [PMID: 38756003 PMCID: PMC11145653 DOI: 10.1021/acs.jpcb.4c01637] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/12/2024] [Revised: 04/29/2024] [Accepted: 05/06/2024] [Indexed: 05/18/2024]
Abstract
In this study, the site energy fluctuations, energy transfer dynamics, and some spectroscopic properties of the minor light-harvesting complex CP24 in a membrane environment were determined. For this purpose, a 3 μs-long classical molecular dynamics simulation was performed for the CP24 complex. Furthermore, using the density functional tight binding/molecular mechanics molecular dynamics (DFTB/MM MD) approach, we performed excited state calculations for the chlorophyll a and chlorophyll b molecules in the complex starting from five different positions of the MD trajectory. During the extended simulations, we observed variations in the site energies of the different sets as a result of the fluctuating protein environment. In particular, a water coordination to Chl-b 608 occurred only after about 1 μs in the simulations, demonstrating dynamic changes in the environment of this pigment. From the classical and the DFTB/MM MD simulations, spectral densities and the (time-dependent) Hamiltonian of the complex were determined. Based on these results, three independent strongly coupled chlorophyll clusters were revealed within the complex. In addition, absorption and fluorescence spectra were determined together with the exciton relaxation dynamics, which reasonably well agrees with experimental time scales.
Collapse
Affiliation(s)
- Pooja Sarngadharan
- School of Science, Constructor
University, Campus Ring
1, 28759 Bremen, Germany
| | - Yannick Holtkamp
- School of Science, Constructor
University, Campus Ring
1, 28759 Bremen, Germany
| | | |
Collapse
|
7
|
Ahad S, Lin C, Reppert M. PigmentHunter: A point-and-click application for automated chlorophyll-protein simulations. J Chem Phys 2024; 160:154111. [PMID: 38639311 DOI: 10.1063/5.0198443] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/18/2024] [Accepted: 03/21/2024] [Indexed: 04/20/2024] Open
Abstract
Chlorophyll proteins (CPs) are the workhorses of biological photosynthesis, working together to absorb solar energy, transfer it to chemically active reaction centers, and control the charge-separation process that drives its storage as chemical energy. Yet predicting CP optical and electronic properties remains a serious challenge, driven by the computational difficulty of treating large, electronically coupled molecular pigments embedded in a dynamically structured protein environment. To address this challenge, we introduce here an analysis tool called PigmentHunter, which automates the process of preparing CP structures for molecular dynamics (MD), running short MD simulations on the nanoHUB.org science gateway, and then using electrostatic and steric analysis routines to predict optical absorption, fluorescence, and circular dichroism spectra within a Frenkel exciton model. Inter-pigment couplings are evaluated using point-dipole or transition-charge coupling models, while site energies can be estimated using both electrostatic and ring-deformation approaches. The package is built in a Jupyter Notebook environment, with a point-and-click interface that can be used either to manually prepare individual structures or to batch-process many structures at once. We illustrate PigmentHunter's capabilities with example simulations on spectral line shapes in the light harvesting 2 complex, site energies in the Fenna-Matthews-Olson protein, and ring deformation in photosystems I and II.
Collapse
Affiliation(s)
- S Ahad
- Department of Chemistry, Purdue University, West Lafayette, Indiana 47907, USA
| | - C Lin
- Department of Chemistry, Purdue University, West Lafayette, Indiana 47907, USA
| | - M Reppert
- Department of Chemistry, Purdue University, West Lafayette, Indiana 47907, USA
| |
Collapse
|
8
|
Cignoni E, Suman D, Nigam J, Cupellini L, Mennucci B, Ceriotti M. Electronic Excited States from Physically Constrained Machine Learning. ACS CENTRAL SCIENCE 2024; 10:637-648. [PMID: 38559300 PMCID: PMC10979507 DOI: 10.1021/acscentsci.3c01480] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 11/30/2023] [Revised: 01/16/2024] [Accepted: 01/30/2024] [Indexed: 04/04/2024]
Abstract
Data-driven techniques are increasingly used to replace electronic-structure calculations of matter. In this context, a relevant question is whether machine learning (ML) should be applied directly to predict the desired properties or combined explicitly with physically grounded operations. We present an example of an integrated modeling approach in which a symmetry-adapted ML model of an effective Hamiltonian is trained to reproduce electronic excitations from a quantum-mechanical calculation. The resulting model can make predictions for molecules that are much larger and more complex than those on which it is trained and allows for dramatic computational savings by indirectly targeting the outputs of well-converged calculations while using a parametrization corresponding to a minimal atom-centered basis. These results emphasize the merits of intertwining data-driven techniques with physical approximations, improving the transferability and interpretability of ML models without affecting their accuracy and computational efficiency and providing a blueprint for developing ML-augmented electronic-structure methods.
Collapse
Affiliation(s)
- Edoardo Cignoni
- Dipartimento
di Chimica e Chimica Industriale, Università
di Pisa, 56126 Pisa, Italy
| | - Divya Suman
- Laboratory
of Computational Science and Modeling, Institut des Matériaux, École Polytechnique Fédérale
de Lausanne, 1015 Lausanne, Switzerland
| | - Jigyasa Nigam
- Laboratory
of Computational Science and Modeling, Institut des Matériaux, École Polytechnique Fédérale
de Lausanne, 1015 Lausanne, Switzerland
| | - Lorenzo Cupellini
- Dipartimento
di Chimica e Chimica Industriale, Università
di Pisa, 56126 Pisa, Italy
| | - Benedetta Mennucci
- Dipartimento
di Chimica e Chimica Industriale, Università
di Pisa, 56126 Pisa, Italy
| | - Michele Ceriotti
- Laboratory
of Computational Science and Modeling, Institut des Matériaux, École Polytechnique Fédérale
de Lausanne, 1015 Lausanne, Switzerland
- Division
of Chemistry and Chemical Engineering, California
Institute of Technology, Pasadena, California 91125, United States
| |
Collapse
|
9
|
Allan L, Zuehlsdorff TJ. Taming the third order cumulant approximation to linear optical spectroscopy. J Chem Phys 2024; 160:074108. [PMID: 38380749 DOI: 10.1063/5.0182745] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/20/2023] [Accepted: 01/23/2024] [Indexed: 02/22/2024] Open
Abstract
The second order cumulant method offers a promising pathway to predicting optical properties in condensed phase systems. It allows for the computation of linear absorption spectra from excitation energy fluctuations sampled along molecular dynamics (MD) trajectories, fully accounting for vibronic effects, direct solute-solvent interactions, and environmental polarization effects. However, the second order cumulant approximation only guarantees accurate line shapes for energy gap fluctuations obeying Gaussian statistics. A third order correction has recently been derived but often yields unphysical spectra or divergent line shapes for moderately non-Gaussian fluctuations due to the neglect of higher order terms in the cumulant expansion. In this work, we develop a corrected cumulant approach, where the collective effect of neglected higher order contributions is approximately accounted for through a dampening factor applied to the third order cumulant term. We show that this dampening factor can be expressed as a function of the skewness and kurtosis of energy gap fluctuations and can be parameterized from a large set of randomly sampled model Hamiltonians for which exact spectral line shapes are known. This approach is shown to systematically remove unphysical contributions in the form of negative absorbances from cumulant spectra in both model Hamiltonians and condensed phase systems sampled from MD and dramatically improves over the second order cumulant method in describing systems exhibiting Duschinsky mode mixing effects. We successfully apply the approach to the coumarin-153 dye in toluene, obtaining excellent agreement with experiment.
Collapse
Affiliation(s)
- Lucas Allan
- Department of Chemistry, Oregon State University, Corvallis, Oregon 97331, USA
| | - Tim J Zuehlsdorff
- Department of Chemistry, Oregon State University, Corvallis, Oregon 97331, USA
| |
Collapse
|
10
|
Bose A, Walters PL. Impact of Spatial Inhomogeneity on Excitation Energy Transport in the Fenna-Matthews-Olson Complex. J Phys Chem B 2023; 127:7663-7673. [PMID: 37647510 DOI: 10.1021/acs.jpcb.3c03062] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/01/2023]
Abstract
The dynamics of the excitation energy transfer (EET) in photosynthetic complexes is an interesting question both from the perspective of fundamental understanding and the research in artificial photosynthesis. Over the past decade, very accurate spectral densities have been developed to capture spatial inhomogeneities in the Fenna-Matthews-Olson (FMO) complex. However, challenges persist in numerically simulating these systems, both in terms of parameterizing them and following their dynamics over long periods of time because of long non-Markovian memories. We investigate the dynamics of FMO with the exact treatment of various theoretical spectral densities using the new tensor network path integral-based methods, which are uniquely capable of addressing the difficulty of long memory length and incoherent Förster theory. It is also important to be able to analyze the pathway of EET flow, which can be difficult to identify given the non-trivial structure of connections between bacteriochlorophyll molecules in FMO. We use the recently introduced ideas of relating coherence to population derivatives to analyze the transport process and reveal some new routes of transport. The combination of exact and approximate methods sheds light on the role of coherences in affecting the fine details of the transport and promises to be a powerful toolbox for future exploration of other open systems with quantum transport.
Collapse
Affiliation(s)
- Amartya Bose
- Department of Chemical Sciences, Tata Institute of Fundamental Research, Mumbai 400005, India
| | - Peter L Walters
- Department of Chemistry, University of California, Berkeley, California 94720, United States
- Miller Institute for Basic Research in Science, University of California Berkeley, Berkeley, California 94720, United States
| |
Collapse
|
11
|
Kim Y, Mitchell Z, Lawrence J, Morozov D, Savikhin S, Slipchenko LV. Predicting Mutation-Induced Changes in the Electronic Properties of Photosynthetic Proteins from First Principles: The Fenna-Matthews-Olson Complex Example. J Phys Chem Lett 2023; 14:7038-7044. [PMID: 37524046 DOI: 10.1021/acs.jpclett.3c01461] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/02/2023]
Abstract
Multiscale molecular modeling is utilized to predict optical absorption and circular dichroism spectra of two single-point mutants of the Fenna-Matthews-Olson photosynthetic pigment-protein complex. The modeling approach combines classical molecular dynamics simulations with structural refinement of photosynthetic pigments and calculations of their excited states in a polarizable protein environment. The only experimental input to the modeling protocol is the X-ray structure of the wild-type protein. The first-principles modeling reproduces changes in the experimental optical spectra of the considered mutants, Y16F and Q198V. Interestingly, the Q198V mutation has a negligible effect on the electronic properties of the targeted bacteriochlorophyll a pigment. Instead, the electronic properties of several other pigments respond to this mutation. The molecular modeling demonstrates that a single-point mutation can induce long-range effects on the protein structure, while extensive structural changes near a pigment do not necessarily lead to significant changes in the electronic properties of that pigment.
Collapse
Affiliation(s)
- Yongbin Kim
- Department of Chemistry, Purdue University, 560 Oval Drive, West Lafayette, Indiana 47907, United States
| | - Zach Mitchell
- Department of Physics and Astronomy, Purdue University, 525 Northwestern Avenue, West Lafayette, Indiana 47907, United States
| | - Jack Lawrence
- Department of Chemistry, Purdue University, 560 Oval Drive, West Lafayette, Indiana 47907, United States
| | - Dmitry Morozov
- Nanoscience Center and Department of Chemistry, University of Jyväskylä, P.O. Box 35, 40014 Jyväskylä, Finland
| | - Sergei Savikhin
- Department of Physics and Astronomy, Purdue University, 525 Northwestern Avenue, West Lafayette, Indiana 47907, United States
| | - Lyudmila V Slipchenko
- Department of Chemistry, Purdue University, 560 Oval Drive, West Lafayette, Indiana 47907, United States
| |
Collapse
|
12
|
Sarngadharan P, Maity S, Kleinekathöfer U. Spectral densities and absorption spectra of the core antenna complex CP43 from photosystem II. J Chem Phys 2022; 156:215101. [DOI: 10.1063/5.0091005] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Besides absorbing light, the core antenna complex CP43 of photosystem II is of great importance in transferring excitation energy from the antenna complexes to the reaction center. Excitation energies, spectral densities, and linear absorption spectra of the complex have been evaluated by a multiscale approach. In this scheme, quantum mechanics/molecular mechanics molecular dynamics simulations are performed employing the parameterized density functional tight binding (DFTB) while the time-dependent long-range-corrected DFTB scheme is applied for the excited state calculations. The obtained average spectral density of the CP43 complex shows a very good agreement with experimental results. Moreover, the excitonic Hamiltonian of the system along with the computed site-dependent spectral densities was used to determine the linear absorption. While a Redfield-like approximation has severe shortcomings in dealing with the CP43 complex due to quasi-degenerate states, the non-Markovian full second-order cumulant expansion formalism is able to overcome the drawbacks. Linear absorption spectra were obtained, which show a good agreement with the experimental counterparts at different temperatures. This study once more emphasizes that by combining diverse techniques from the areas of molecular dynamics simulations, quantum chemistry, and open quantum systems, it is possible to obtain first-principle results for photosynthetic complexes, which are in accord with experimental findings.
Collapse
Affiliation(s)
- Pooja Sarngadharan
- Department of Physics and Earth Sciences, Jacobs University Bremen, Campus Ring 1, 28759 Bremen, Germany
| | - Sayan Maity
- Department of Physics and Earth Sciences, Jacobs University Bremen, Campus Ring 1, 28759 Bremen, Germany
| | - Ulrich Kleinekathöfer
- Department of Physics and Earth Sciences, Jacobs University Bremen, Campus Ring 1, 28759 Bremen, Germany
| |
Collapse
|
13
|
Cignoni E, Slama V, Cupellini L, Mennucci B. The atomistic modeling of light-harvesting complexes from the physical models to the computational protocol. J Chem Phys 2022; 156:120901. [DOI: 10.1063/5.0086275] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022] Open
Abstract
The function of light-harvesting complexes is determined by a complex network of dynamic interactions among all the different components: the aggregate of pigments, the protein, and the surrounding environment. Complete and reliable predictions on these types of composite systems can be only achieved with an atomistic description. In the last few decades, there have been important advances in the atomistic modeling of light-harvesting complexes. These advances have involved both the completeness of the physical models and the accuracy and effectiveness of the computational protocols. In this Perspective, we present an overview of the main theoretical and computational breakthroughs attained so far in the field, with particular focus on the important role played by the protein and its dynamics. We then discuss the open problems in their accurate modeling that still need to be addressed. To illustrate an effective computational workflow for the modeling of light harvesting complexes, we take as an example the plant antenna complex CP29 and its H111N mutant.
Collapse
Affiliation(s)
- Edoardo Cignoni
- Dipartimento di Chimica e Chimica Industriale, University of Pisa, via G. Moruzzi 13, 56124 Pisa, Italy
| | - Vladislav Slama
- Dipartimento di Chimica e Chimica Industriale, University of Pisa, via G. Moruzzi 13, 56124 Pisa, Italy
| | - Lorenzo Cupellini
- Dipartimento di Chimica e Chimica Industriale, University of Pisa, via G. Moruzzi 13, 56124 Pisa, Italy
| | - Benedetta Mennucci
- Dipartimento di Chimica e Chimica Industriale, University of Pisa, via G. Moruzzi 13, 56124 Pisa, Italy
| |
Collapse
|
14
|
Cho KH, Rhee YM. Computational elucidations on the role of vibrations in energy transfer processes of photosynthetic complexes. Phys Chem Chem Phys 2021; 23:26623-26639. [PMID: 34842245 DOI: 10.1039/d1cp04615b] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Coupling between pigment excitations and nuclear movements in photosynthetic complexes is known to modulate the excitation energy transfer (EET) efficiencies. Toward providing microscopic information, researchers often apply simulation techniques and investigate how vibrations are involved in EET processes. Here, reports on such roles of nuclear movements are discussed from a theory perspective. While vibrations naturally present random thermal fluctuations that can affect energy transferring characteristics, they can also be intertwined with exciton structures and create more specific non-adiabatic energy transfer pathways. For reliable simulations, a bath model that accurately mimics a given molecular system is required. Methods for obtaining such a model in combination with quantum chemical electronic structure calculations and molecular dynamics trajectory simulations are discussed. Various quantum dynamics simulation tools that can handle pigment-to-pigment energy transfers together with their vibrational characters are also touched on. Behaviors of molecular vibrations often deviate from ideality, especially when all-atom details are included, which practically forces us to treat them classically. We conclude this perspective by considering some recent reports that suggest that classical descriptions of bath effects with all-atom details may still produce valuable information for analyzing sophisticated contributions by vibrations to EET processes.
Collapse
Affiliation(s)
- Kwang Hyun Cho
- Department of Chemistry, Korea Advanced Institute of Science and Technology (KAIST), Daejeon 34141, Korea.
| | - Young Min Rhee
- Department of Chemistry, Korea Advanced Institute of Science and Technology (KAIST), Daejeon 34141, Korea.
| |
Collapse
|
15
|
Lu SY, Zuehlsdorff TJ, Hong H, Aguirre VP, Isborn CM, Shi L. The Influence of Electronic Polarization on Nonlinear Optical Spectroscopy. J Phys Chem B 2021; 125:12214-12227. [PMID: 34726915 DOI: 10.1021/acs.jpcb.1c05914] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The environment surrounding a chromophore can dramatically affect the energy absorption and relaxation process, as manifested in optical spectra. Simulations of nonlinear optical spectroscopy, such as two-dimensional electronic spectroscopy (2DES) and transient absorption (TA), will be influenced by the computational model of the environment. We here compare a fixed point charge molecular mechanics model and a quantum mechanical (QM) model of the environment in computed 2DES and TA spectra of Nile red in water and the chromophore of photoactive yellow protein (PYP) in water and protein environments. In addition to simulating these nonlinear optical spectra, we directly juxtapose the computed excitation energy correlation function to the dynamic Stokes shift function often used to analyze environment dynamics. Overall, we find that for the three systems studied here the mutual electronic polarization provided by the QM environment manifests in broader 2DES signals, as well as a larger reorganization energy and a larger static Stokes shift due to stronger coupling between the chromophore and the environment.
Collapse
Affiliation(s)
- Shao-Yu Lu
- Department of Chemistry and Biochemistry, University of California Merced, Merced, California 95343, United States
| | - Tim J Zuehlsdorff
- Department of Chemistry, Oregon State University, Corvallis, Oregon 97331, United States
| | - Hanbo Hong
- Department of Chemistry and Biochemistry, University of California Merced, Merced, California 95343, United States
| | - Vincent P Aguirre
- Department of Chemistry and Biochemistry, University of California Merced, Merced, California 95343, United States
| | - Christine M Isborn
- Department of Chemistry and Biochemistry, University of California Merced, Merced, California 95343, United States
| | - Liang Shi
- Department of Chemistry and Biochemistry, University of California Merced, Merced, California 95343, United States
| |
Collapse
|
16
|
Abstract
Numerous linear and non-linear spectroscopic techniques have been developed to elucidate structural and functional information of complex systems ranging from natural systems, such as proteins and light-harvesting systems, to synthetic systems, such as solar cell materials and light-emitting diodes. The obtained experimental data can be challenging to interpret due to the complexity and potential overlapping spectral signatures. Therefore, computational spectroscopy plays a crucial role in the interpretation and understanding of spectral observables of complex systems. Computational modeling of various spectroscopic techniques has seen significant developments in the past decade, when it comes to the systems that can be addressed, the size and complexity of the sample types, the accuracy of the methods, and the spectroscopic techniques that can be addressed. In this Perspective, I will review the computational spectroscopy methods that have been developed and applied for infrared and visible spectroscopies in the condensed phase. I will discuss some of the questions that this has allowed answering. Finally, I will discuss current and future challenges and how these may be addressed.
Collapse
Affiliation(s)
- Thomas L C Jansen
- Zernike Institute for Advanced Materials, University of Groningen, Nijenborgh 4, 9747 AG Groningen, The Netherlands
| |
Collapse
|
17
|
Srivastava A, Ahad S, Wat JH, Reppert M. Accurate prediction of mutation-induced frequency shifts in chlorophyll proteins with a simple electrostatic model. J Chem Phys 2021; 155:151102. [PMID: 34686046 DOI: 10.1063/5.0064567] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Photosynthetic pigment-protein complexes control local chlorophyll (Chl) transition frequencies through a variety of electrostatic and steric forces. Site-directed mutations can modify this local spectroscopic tuning, providing critical insight into native photosynthetic functions and offering the tantalizing prospect of creating rationally designed Chl proteins with customized optical properties. Unfortunately, at present, no proven methods exist for reliably predicting mutation-induced frequency shifts in advance, limiting the method's utility for quantitative applications. Here, we address this challenge by constructing a series of point mutants in the water-soluble chlorophyll protein of Lepidium virginicum and using them to test the reliability of a simple computational protocol for mutation-induced site energy shifts. The protocol uses molecular dynamics to prepare mutant protein structures and the charge density coupling model of Adolphs et al. [Photosynth. Res. 95, 197-209 (2008)] for site energy prediction; a graphical interface that implements the protocol automatically is published online at http://nanohub.org/tools/pigmenthunter. With the exception of a single outlier (presumably due to unexpected structural changes), we find that the calculated frequency shifts match the experiment remarkably well, with an average error of 1.6 nm over a 9 nm spread in wavelengths. We anticipate that the accuracy of the method can be improved in the future with more advanced sampling of mutant protein structures.
Collapse
Affiliation(s)
- Amit Srivastava
- Department of Chemistry, Purdue University, West Lafayette, Indiana 47907, USA
| | - Safa Ahad
- Department of Chemistry, Purdue University, West Lafayette, Indiana 47907, USA
| | - Jacob H Wat
- Department of Chemistry, Purdue University, West Lafayette, Indiana 47907, USA
| | - Mike Reppert
- Department of Chemistry, Purdue University, West Lafayette, Indiana 47907, USA
| |
Collapse
|
18
|
Dunnett AJ, Gowland D, Isborn CM, Chin AW, Zuehlsdorff TJ. Influence of non-adiabatic effects on linear absorption spectra in the condensed phase: Methylene blue. J Chem Phys 2021; 155:144112. [PMID: 34654312 DOI: 10.1063/5.0062950] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023] Open
Abstract
Modeling linear absorption spectra of solvated chromophores is highly challenging as contributions are present both from coupling of the electronic states to nuclear vibrations and from solute-solvent interactions. In systems where excited states intersect in the Condon region, significant non-adiabatic contributions to absorption line shapes can also be observed. Here, we introduce a robust approach to model linear absorption spectra accounting for both environmental and non-adiabatic effects from first principles. This model parameterizes a linear vibronic coupling (LVC) Hamiltonian directly from energy gap fluctuations calculated along molecular dynamics (MD) trajectories of the chromophore in solution, accounting for both anharmonicity in the potential and direct solute-solvent interactions. The resulting system dynamics described by the LVC Hamiltonian are solved exactly using the thermalized time-evolving density operator with orthogonal polynomials algorithm (T-TEDOPA). The approach is applied to the linear absorption spectrum of methylene blue in water. We show that the strong shoulder in the experimental spectrum is caused by vibrationally driven population transfer between the bright S1 and the dark S2 states. The treatment of the solvent environment is one of many factors that strongly influence the population transfer and line shape; accurate modeling can only be achieved through the use of explicit quantum mechanical solvation. The efficiency of T-TEDOPA, combined with LVC Hamiltonian parameterizations from MD, leads to an attractive method for describing a large variety of systems in complex environments from first principles.
Collapse
Affiliation(s)
- Angus J Dunnett
- Sorbonne Université, CNRS, Institut des NanoSciences de Paris, 4 place Jussieu, 75005 Paris, France
| | - Duncan Gowland
- Department of Physics, King's College London, London WC2R 2LS, United Kingdom
| | - Christine M Isborn
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, USA
| | - Alex W Chin
- Sorbonne Université, CNRS, Institut des NanoSciences de Paris, 4 place Jussieu, 75005 Paris, France
| | - Tim J Zuehlsdorff
- Department of Chemistry, Oregon State University, Corvallis, Oregon 97331, USA
| |
Collapse
|
19
|
Maity S, Sarngadharan P, Daskalakis V, Kleinekathöfer U. Time-dependent atomistic simulations of the CP29 light-harvesting complex. J Chem Phys 2021; 155:055103. [PMID: 34364345 DOI: 10.1063/5.0053259] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023] Open
Abstract
Light harvesting as the first step in photosynthesis is of prime importance for life on earth. For a theoretical description of photochemical processes during light harvesting, spectral densities are key quantities. They serve as input functions for modeling the excitation energy transfer dynamics and spectroscopic properties. Herein, a recently developed procedure is applied to determine the spectral densities of the pigments in the minor antenna complex CP29 of photosystem II, which has recently gained attention because of its active role in non-photochemical quenching processes in higher plants. To this end, the density functional-based tight binding (DFTB) method has been employed to enable simulation of the ground state dynamics in a quantum-mechanics/molecular mechanics (QM/MM) scheme for each chlorophyll pigment. Subsequently, the time-dependent extension of the long-range corrected DFTB approach has been used to obtain the excitation energy fluctuations along the ground-state trajectories also in a QM/MM setting. From these results, the spectral densities have been determined and compared for different force fields and to spectral densities from other light-harvesting complexes. In addition, time-dependent and time-independent excitonic Hamiltonians of the system have been constructed and applied to the determination of absorption spectra as well as exciton dynamics.
Collapse
Affiliation(s)
- Sayan Maity
- Department of Physics and Earth Sciences, Jacobs University Bremen, Campus Ring 1, 28759 Bremen, Germany
| | - Pooja Sarngadharan
- Department of Physics and Earth Sciences, Jacobs University Bremen, Campus Ring 1, 28759 Bremen, Germany
| | - Vangelis Daskalakis
- Department of Chemical Engineering, Cyprus University of Technology, 30 Archbishop Kyprianou Str. 3603, Limassol, Cyprus
| | - Ulrich Kleinekathöfer
- Department of Physics and Earth Sciences, Jacobs University Bremen, Campus Ring 1, 28759 Bremen, Germany
| |
Collapse
|
20
|
Zuehlsdorff TJ, Shedge SV, Lu SY, Hong H, Aguirre VP, Shi L, Isborn CM. Vibronic and Environmental Effects in Simulations of Optical Spectroscopy. Annu Rev Phys Chem 2021; 72:165-188. [DOI: 10.1146/annurev-physchem-090419-051350] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
Including both environmental and vibronic effects is important for accurate simulation of optical spectra, but combining these effects remains computationally challenging. We outline two approaches that consider both the explicit atomistic environment and the vibronic transitions. Both phenomena are responsible for spectral shapes in linear spectroscopy and the electronic evolution measured in nonlinear spectroscopy. The first approach utilizes snapshots of chromophore-environment configurations for which chromophore normal modes are determined. We outline various approximations for this static approach that assumes harmonic potentials and ignores dynamic system-environment coupling. The second approach obtains excitation energies for a series of time-correlated snapshots. This dynamic approach relies on the accurate truncation of the cumulant expansion but treats the dynamics of the chromophore and the environment on equal footing. Both approaches show significant potential for making strides toward more accurate optical spectroscopy simulations of complex condensed phase systems.
Collapse
Affiliation(s)
- Tim J. Zuehlsdorff
- Department of Chemistry, Oregon State University, Corvallis, Oregon 97331, USA
| | - Sapana V. Shedge
- Department of Chemistry and Chemical Biology, University of California, Merced, California 95343, USA
| | - Shao-Yu Lu
- Department of Chemistry and Chemical Biology, University of California, Merced, California 95343, USA
| | - Hanbo Hong
- Department of Chemistry and Chemical Biology, University of California, Merced, California 95343, USA
| | - Vincent P. Aguirre
- Department of Chemistry and Chemical Biology, University of California, Merced, California 95343, USA
| | - Liang Shi
- Department of Chemistry and Chemical Biology, University of California, Merced, California 95343, USA
| | - Christine M. Isborn
- Department of Chemistry and Chemical Biology, University of California, Merced, California 95343, USA
| |
Collapse
|
21
|
Maity S, Daskalakis V, Elstner M, Kleinekathöfer U. Multiscale QM/MM molecular dynamics simulations of the trimeric major light-harvesting complex II. Phys Chem Chem Phys 2021; 23:7407-7417. [PMID: 33876100 DOI: 10.1039/d1cp01011e] [Citation(s) in RCA: 28] [Impact Index Per Article: 9.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
Photosynthetic processes are driven by sunlight. Too little of it and the photosynthetic machinery cannot produce the reductive power to drive the anabolic pathways. Too much sunlight and the machinery can get damaged. In higher plants, the major Light-Harvesting Complex (LHCII) efficiently absorbs the light energy, but can also dissipate it when in excess (quenching). In order to study the dynamics related to the quenching process but also the exciton dynamics in general, one needs to accurately determine the so-called spectral density which describes the coupling between the relevant pigment modes and the environmental degrees of freedom. To this end, Born-Oppenheimer molecular dynamics simulations in a quantum mechanics/molecular mechanics (QM/MM) fashion utilizing the density functional based tight binding (DFTB) method have been performed for the ground state dynamics. Subsequently, the time-dependent extension of the long-range-corrected DFTB scheme has been employed for the excited state calculations of the individual chlorophyll-a molecules in the LHCII complex. The analysis of this data resulted in spectral densities showing an astonishing agreement with the experimental counterpart in this rather large system. This consistency with an experimental observable also supports the accuracy, robustness, and reliability of the present multi-scale scheme. To the best of our knowledge, this is the first theoretical attempt on this large complex system is ever made to accurately simulate the spectral density. In addition, the resulting spectral densities and site energies were used to determine the exciton transfer rate within a special pigment pair consisting of a chlorophyll-a and a carotenoid molecule which is assumed to play a role in the balance between the light harvesting and quenching modes.
Collapse
Affiliation(s)
- Sayan Maity
- Department of Physics and Earth Sciences, Jacobs University Bremen, Campus Ring 1, 28759 Bremen, Germany.
| | | | | | | |
Collapse
|
22
|
Chaillet M, Lengauer F, Adolphs J, Müh F, Fokas AS, Cole DJ, Chin AW, Renger T. Static Disorder in Excitation Energies of the Fenna-Matthews-Olson Protein: Structure-Based Theory Meets Experiment. J Phys Chem Lett 2020; 11:10306-10314. [PMID: 33227205 PMCID: PMC7751012 DOI: 10.1021/acs.jpclett.0c03123] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/14/2020] [Accepted: 11/17/2020] [Indexed: 06/11/2023]
Abstract
Inhomogeneous broadening of optical lines of the Fenna-Matthews-Olson (FMO) light-harvesting protein is investigated by combining a Monte Carlo sampling of low-energy conformational substates of the protein with a quantum chemical/electrostatic calculation of local transition energies (site energies) of the pigments. The good agreement between the optical spectra calculated for the inhomogeneous ensemble and the experimental data demonstrates that electrostatics is the dominant contributor to static disorder in site energies. Rotamers of polar amino acid side chains are found to cause bimodal distribution functions of site energy shifts, which can be probed by hole burning and single-molecule spectroscopy. When summing over the large number of contributions, the resulting distribution functions of the site energies become Gaussians, and the correlations in site energy fluctuations at different sites practically average to zero. These results demonstrate that static disorder in the FMO protein is in the realm of the central limit theorem of statistics.
Collapse
Affiliation(s)
- Marten
L. Chaillet
- Bijvoet
Centre for Biomolecular Research, University
of Utrecht, Heidelberglaan 8, 3584 CS Utrecht, The Netherlands
| | - Florian Lengauer
- Institute
of Theoretical Physics, Johannes Kepler
University Linz, Altenberger Str. 69, 4040 Linz, Austria
| | - Julian Adolphs
- Leibniz
Institute for Agricultural Engineering and Bioeconomy, Max-Eyth-Allee 100, 14469 Potsdam, Germany
| | - Frank Müh
- Institute
of Theoretical Physics, Johannes Kepler
University Linz, Altenberger Str. 69, 4040 Linz, Austria
| | - Alexander S. Fokas
- TCM
Group, Cavendish Laboratory, 19 J. J. Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| | - Daniel J. Cole
- School
of Natural and Environmental Sciences, Newcastle
University, Newcastle
upon Tyne NE1 7RU, United
Kingdom
| | - Alex W. Chin
- Centre
National de la Recherce Scientifique, Institute des Nanosciences de
Paris, Sorbonne Université, Paris, France
| | - Thomas Renger
- Institute
of Theoretical Physics, Johannes Kepler
University Linz, Altenberger Str. 69, 4040 Linz, Austria
| |
Collapse
|
23
|
Klinger A, Lindorfer D, Müh F, Renger T. Normal mode analysis of spectral density of FMO trimers: Intra- and intermonomer energy transfer. J Chem Phys 2020; 153:215103. [PMID: 33291900 DOI: 10.1063/5.0027994] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022] Open
Abstract
The intermolecular contribution to the spectral density of the exciton-vibrational coupling of the homotrimeric Fenna-Matthews-Olson (FMO) light-harvesting protein of green sulfur bacteria P. aestuarii is analyzed by combining a normal mode analysis of the protein with the charge density coupling method for the calculation of local transition energies of the pigments. Correlations in site energy fluctuations across the whole FMO trimer are found at low vibrational frequencies. Including, additionally, the high-frequency intrapigment part of the spectral density, extracted from line-narrowing spectra, we study intra- and intermonomer exciton transfer. Whereas the intrapigment part of the spectral density is important for fast intramonomer exciton relaxation, the intermolecular contributions (due to pigment-environment coupling) determine the intermonomer exciton transfer. Neither the variations of the local Huang-Rhys factors nor the correlations in site energy fluctuations have a critical influence on energy transfer. At room temperature, the intermonomer transfer in the FMO protein occurs on a 10 ps time scale, whereas intramonomer exciton equilibration is roughly two orders of magnitude faster. At cryogenic temperatures, intermonomer transfer limits the lifetimes of the lowest exciton band. The lifetimes are found to increase between 20 ps in the center of this band up to 100 ps toward lower energies, which is in very good agreement with the estimates from hole burning data. Interestingly, exciton delocalization in the FMO monomers is found to slow down intermonomer energy transfer, at both physiological and cryogenic temperatures.
Collapse
Affiliation(s)
- Alexander Klinger
- Institut für Theoretische Physik, Johannes Kepler Universität Linz, Altenberger Str. 69, 4040 Linz, Austria
| | - Dominik Lindorfer
- Institut für Theoretische Physik, Johannes Kepler Universität Linz, Altenberger Str. 69, 4040 Linz, Austria
| | - Frank Müh
- Institut für Theoretische Physik, Johannes Kepler Universität Linz, Altenberger Str. 69, 4040 Linz, Austria
| | - Thomas Renger
- Institut für Theoretische Physik, Johannes Kepler Universität Linz, Altenberger Str. 69, 4040 Linz, Austria
| |
Collapse
|
24
|
Maity S, Bold BM, Prajapati JD, Sokolov M, Kubař T, Elstner M, Kleinekathöfer U. DFTB/MM Molecular Dynamics Simulations of the FMO Light-Harvesting Complex. J Phys Chem Lett 2020; 11:8660-8667. [PMID: 32991176 DOI: 10.1021/acs.jpclett.0c02526] [Citation(s) in RCA: 39] [Impact Index Per Article: 9.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
Because of the size of light-harvesting complexes and the involvement of electronic degrees of freedom, computationally these systems need to be treated with a combined quantum-classical description. To this end, Born-Oppenheimer molecular dynamics simulations have been employed in a quantum mechanics/molecular mechanics (QM/MM) fashion for the ground state followed by excitation energy calculations again in a QM/MM scheme for the Fenna-Matthews-Olson (FMO) complex. The self-consistent-charge density functional tight-binding (DFTB) method electrostatically coupled to a classical description of the environment was applied to perform the ground-state dynamics. Subsequently, long-range-corrected time-dependent DFTB calculations were performed to determine the excitation energy fluctuations of the individual bacteriochlorophyll a molecules. The spectral densities obtained using this approach show an excellent agreement with experimental findings. In addition, the fluctuating site energies and couplings were used to estimate the exciton transfer dynamics.
Collapse
Affiliation(s)
- Sayan Maity
- Department of Physics and Earth Sciences, Jacobs University Bremen, Campus Ring 1, 28759 Bremen, Germany
| | - Beatrix M Bold
- Institute of Physical Chemistry, Karlsruhe Institute of Technology, Kaiserstrasse 12, 76131 Karlsruhe, Germany
| | | | - Monja Sokolov
- Institute of Physical Chemistry, Karlsruhe Institute of Technology, Kaiserstrasse 12, 76131 Karlsruhe, Germany
| | - Tomáš Kubař
- Institute of Physical Chemistry, Karlsruhe Institute of Technology, Kaiserstrasse 12, 76131 Karlsruhe, Germany
| | - Marcus Elstner
- Institute of Physical Chemistry, Karlsruhe Institute of Technology, Kaiserstrasse 12, 76131 Karlsruhe, Germany
- Institute of Biological Interfaces, Karlsruhe Institute of Technology, Kaiserstrasse 12, 76131 Karlsruhe, Germany
| | - Ulrich Kleinekathöfer
- Department of Physics and Earth Sciences, Jacobs University Bremen, Campus Ring 1, 28759 Bremen, Germany
| |
Collapse
|
25
|
Analysis of Photosynthetic Systems and Their Applications with Mathematical and Computational Models. APPLIED SCIENCES-BASEL 2020. [DOI: 10.3390/app10196821] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
In biological and life science applications, photosynthesis is an important process that involves the absorption and transformation of sunlight into chemical energy. During the photosynthesis process, the light photons are captured by the green chlorophyll pigments in their photosynthetic antennae and further funneled to the reaction center. One of the most important light harvesting complexes that are highly important in the study of photosynthesis is the membrane-attached Fenna–Matthews–Olson (FMO) complex found in the green sulfur bacteria. In this review, we discuss the mathematical formulations and computational modeling of some of the light harvesting complexes including FMO. The most recent research developments in the photosynthetic light harvesting complexes are thoroughly discussed. The theoretical background related to the spectral density, quantum coherence and density functional theory has been elaborated. Furthermore, details about the transfer and excitation of energy in different sites of the FMO complex along with other vital photosynthetic light harvesting complexes have also been provided. Finally, we conclude this review by providing the current and potential applications in environmental science, energy, health and medicine, where such mathematical and computational studies of the photosynthesis and the light harvesting complexes can be readily integrated.
Collapse
|
26
|
Gao X, Lai Y, Geva E. Simulating Absorption Spectra of Multiexcitonic Systems via Quasiclassical Mapping Hamiltonian Methods. J Chem Theory Comput 2020; 16:6465-6480. [DOI: 10.1021/acs.jctc.0c00709] [Citation(s) in RCA: 23] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Xing Gao
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| | - Yifan Lai
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| | - Eitan Geva
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| |
Collapse
|
27
|
Gao X, Geva E. A Nonperturbative Methodology for Simulating Multidimensional Spectra of Multiexcitonic Molecular Systems via Quasiclassical Mapping Hamiltonian Methods. J Chem Theory Comput 2020; 16:6491-6502. [DOI: 10.1021/acs.jctc.0c00843] [Citation(s) in RCA: 23] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Xing Gao
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| | - Eitan Geva
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| |
Collapse
|
28
|
Zuehlsdorff TJ, Hong H, Shi L, Isborn CM. Nonlinear spectroscopy in the condensed phase: The role of Duschinsky rotations and third order cumulant contributions. J Chem Phys 2020; 153:044127. [PMID: 32752702 DOI: 10.1063/5.0013739] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022] Open
Abstract
First-principles modeling of nonlinear optical spectra in the condensed phase is highly challenging because both environment and vibronic interactions can play a large role in determining spectral shapes and excited state dynamics. Here, we compute two dimensional electronic spectroscopy (2DES) signals based on a cumulant expansion of the energy gap fluctuation operator, with specific focus on analyzing mode mixing effects introduced by the Duschinsky rotation and the role of the third order term in the cumulant expansion for both model and realistic condensed phase systems. We show that for a harmonic model system, the third order cumulant correction captures effects introduced by a mismatch in curvatures of ground and excited state potential energy surfaces, as well as effects of mode mixing. We also demonstrate that 2DES signals can be accurately reconstructed from purely classical correlation functions using quantum correction factors. We then compute nonlinear optical spectra for the Nile red and methylene blue chromophores in solution, assessing the third order cumulant contribution for realistic systems. We show that the third order cumulant correction is strongly dependent on the treatment of the solvent environment, revealing the interplay between environmental polarization and the electronic-vibrational coupling.
Collapse
Affiliation(s)
- Tim J Zuehlsdorff
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, USA
| | - Hanbo Hong
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, USA
| | - Liang Shi
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, USA
| | - Christine M Isborn
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, USA
| |
Collapse
|
29
|
Jiang L, Rogers DM, Hirst JD, Do H. Force Fields for Macromolecular Assemblies Containing Diketopyrrolopyrrole and Thiophene. J Chem Theory Comput 2020; 16:5150-5162. [DOI: 10.1021/acs.jctc.0c00399] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Affiliation(s)
- Ling Jiang
- Department of Chemical and Environmental Engineering, University of Nottingham—Ningbo China, Ningbo 315100, China
| | - David M. Rogers
- School of Chemistry, University of Nottingham, University Park, Nottingham NG7 2RD, United Kingdom
| | - Jonathan D. Hirst
- School of Chemistry, University of Nottingham, University Park, Nottingham NG7 2RD, United Kingdom
| | - Hainam Do
- Department of Chemical and Environmental Engineering, University of Nottingham—Ningbo China, Ningbo 315100, China
- New Materials Institute, University of Nottingham—Ningbo China, Ningbo 315042, China
| |
Collapse
|
30
|
Cao J, Cogdell RJ, Coker DF, Duan HG, Hauer J, Kleinekathöfer U, Jansen TLC, Mančal T, Miller RJD, Ogilvie JP, Prokhorenko VI, Renger T, Tan HS, Tempelaar R, Thorwart M, Thyrhaug E, Westenhoff S, Zigmantas D. Quantum biology revisited. SCIENCE ADVANCES 2020; 6:eaaz4888. [PMID: 32284982 PMCID: PMC7124948 DOI: 10.1126/sciadv.aaz4888] [Citation(s) in RCA: 184] [Impact Index Per Article: 46.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/12/2019] [Accepted: 01/06/2020] [Indexed: 05/21/2023]
Abstract
Photosynthesis is a highly optimized process from which valuable lessons can be learned about the operating principles in nature. Its primary steps involve energy transport operating near theoretical quantum limits in efficiency. Recently, extensive research was motivated by the hypothesis that nature used quantum coherences to direct energy transfer. This body of work, a cornerstone for the field of quantum biology, rests on the interpretation of small-amplitude oscillations in two-dimensional electronic spectra of photosynthetic complexes. This Review discusses recent work reexamining these claims and demonstrates that interexciton coherences are too short lived to have any functional significance in photosynthetic energy transfer. Instead, the observed long-lived coherences originate from impulsively excited vibrations, generally observed in femtosecond spectroscopy. These efforts, collectively, lead to a more detailed understanding of the quantum aspects of dissipation. Nature, rather than trying to avoid dissipation, exploits it via engineering of exciton-bath interaction to create efficient energy flow.
Collapse
Affiliation(s)
- Jianshu Cao
- Department of Chemistry, Massachusetts Institute of Technology, 77 Massachusetts Avenue, Cambridge, MA 02139, USA
| | - Richard J. Cogdell
- Institute of Molecular, Cell and Systems Biology, College of Medical, Veterinary and Life Science, University of Glasgow, Glasgow G12 8QQ, UK
| | - David F. Coker
- Department of Chemistry, Boston University, 590 Commonwealth Avenue, Boston, MA 02215, USA
| | - Hong-Guang Duan
- Atomically Resolved Dynamics Department, Max Planck Institute for the Structure and Dynamics of Matter, 22761 Hamburg, Germany
- I. Institut für Theoretische Physik, Universität Hamburg, Jungiusstrasse 9, 20355 Hamburg, Germany
- The Hamburg Centre for Ultrafast Imaging, Universität Hamburg, 22761 Hamburg, Germany
| | - Jürgen Hauer
- Technische Universität München, Dynamische Spektroskopien, Fakultät für Chemie, Lichtenbergstr. 4, 85748 Garching, Germany, and Photonics Institute, TU Wien, 1040 Vienna, Austria
| | - Ulrich Kleinekathöfer
- Department of Physics and Earth Science, Jacobs University Bremen, Campus Ring 1, 28759 Bremen, Germany
| | - Thomas L. C. Jansen
- Zernike Institute for Advanced Materials, University of Groningen, Nijenborgh 4, 9747 AG Groningen, Netherlands
| | - Tomáš Mančal
- Faculty of Mathematics and Physics, Charles University, Ke Karlovu 5, CZ-12116 Prague 2, Czech Republic
| | - R. J. Dwayne Miller
- Atomically Resolved Dynamics Department, Max Planck Institute for the Structure and Dynamics of Matter, 22761 Hamburg, Germany
- The Hamburg Centre for Ultrafast Imaging, Universität Hamburg, 22761 Hamburg, Germany
- Departments of Chemistry and Physics, University of Toronto, Toronto, ON M5S 3H6, Canada
| | | | - Valentyn I. Prokhorenko
- Atomically Resolved Dynamics Department, Max Planck Institute for the Structure and Dynamics of Matter, 22761 Hamburg, Germany
| | - Thomas Renger
- Institute of Theoretical Physics, Department of Theoretical Biophysics, Johannes Kepler University Linz, Altenberger Str. 69, 4040 Linz, Austria
| | - Howe-Siang Tan
- Division of Chemistry and Biological Chemistry, School of Physical and Mathematical Sciences, Nanyang Technological University, 21 Nanyang Link, Singapore 637371, Singapore
| | - Roel Tempelaar
- Department of Chemistry, Columbia University, 3000 Broadway, New York, NY 10027, USA
| | - Michael Thorwart
- I. Institut für Theoretische Physik, Universität Hamburg, Jungiusstrasse 9, 20355 Hamburg, Germany
- The Hamburg Centre for Ultrafast Imaging, Universität Hamburg, 22761 Hamburg, Germany
| | - Erling Thyrhaug
- Technische Universität München, Dynamische Spektroskopien, Fakultät für Chemie, Lichtenbergstr. 4, 85748 Garching, Germany, and Photonics Institute, TU Wien, 1040 Vienna, Austria
| | - Sebastian Westenhoff
- Department of Chemistry and Molecular Biology, University of Gothenburg, Gothenburg 40530, Sweden
| | | |
Collapse
|
31
|
Cui X, Yan Y, Wei J. Theoretical Study on the Effect of Environment on Excitation Energy Transfer in Photosynthetic Light-Harvesting Systems. J Phys Chem B 2020; 124:2354-2362. [PMID: 32130013 DOI: 10.1021/acs.jpcb.0c00266] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
In recent years, two-dimensional (2D) electronic spectroscopy experiments prove that the excitation energy transfer (EET) in photosynthetic light-harvesting systems presents long-lived electronic quantum beating signals. After being discovered in the light-harvesting system, the quantum coherence effect has aroused widespread discussion. To illustrate the EET process in the Fenna-Matthews-Olson (FMO) and phycocyanin 645 (PC645) complex, the local protein environment is often thought to be the same; however, this is ambivalent to the practical structural analysis of the light-harvesting complex. By adopting the dissipaton equation of motion theory, we present the effect of a heterogeneous protein environment on the energy transfer process with accurate numerical results. We demonstrate that the energy transfer process relies on the local heterogeneous environment for the FMO complex. A similar good agreement is found for the PC645 complex. Furthermore, we discuss the optimal value of different chromophores in the excitation energy transfer process by controlling the environmental characteristics.
Collapse
Affiliation(s)
- XueYan Cui
- Department of Physics & Beijing Key Laboratory of Optoelectronic Functional Materials and Micro-nano Devices, Renmin University of China, Beijing 100872, China
| | - YiJing Yan
- Department of Chemical Physics & Hefei National Laboratory for Physical Sciences at the Microscale, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - JianHua Wei
- Department of Physics & Beijing Key Laboratory of Optoelectronic Functional Materials and Micro-nano Devices, Renmin University of China, Beijing 100872, China
| |
Collapse
|
32
|
Kim Y, Morozov D, Stadnytskyi V, Savikhin S, Slipchenko LV. Predictive First-Principles Modeling of a Photosynthetic Antenna Protein: The Fenna-Matthews-Olson Complex. J Phys Chem Lett 2020; 11:1636-1643. [PMID: 32013435 DOI: 10.1021/acs.jpclett.9b03486] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
High efficiency of light harvesting in photosynthetic pigment-protein complexes is governed by evolutionary-perfected protein-assisted tuning of individual pigment properties and interpigment interactions. Due to the large number of spectrally overlapping pigments in a typical photosynthetic complex, experimental methods often fail to unambiguously identify individual chromophore properties. Here, we report a first-principles-based modeling protocol capable of predicting properties of pigments in protein environment to a high precision. The technique was applied to successfully uncover electronic properties of the Fenna-Matthews-Olson (FMO) pigment-protein complex. Each of the three subunits of the FMO complex contains eight strongly coupled bacteriochlorophyll a (BChl a) pigments. The excitonic structure of FMO can be described by an electronic Hamiltonian containing excitation (site) energies of BChl a pigments and electronic couplings between them. Several such Hamiltonians have been developed in the past based on the information from various spectroscopic measurements of FMO; however, fine details of the excitonic structure and energy transfer in FMO, especially assignments of short-lived high-energy sites, remain elusive. Utilizing polarizable embedding quantum mechanics/molecular mechanics with the effective fragment potentials, we computed the electronic Hamiltonian of FMO that is in general agreement with previously reported empirical Hamiltonians and quantitatively reproduces experimental absorption and circular dichroism spectra of the FMO protein. The developed computational protocol is sufficiently simple and can be utilized for predictive modeling of other wild-type and mutated photosynthetic pigment-protein complexes.
Collapse
Affiliation(s)
- Yongbin Kim
- Department of Chemistry, Purdue University, 560 Oval Drive, West Lafayette, Indiana 47907, United States
| | - Dmitry Morozov
- Nanoscience Center and Department of Chemistry, University of Jyväskylä, P.O. Box 35, Jyväskylä 40014, Finland
| | - Valentyn Stadnytskyi
- Department of Physics and Astronomy, Purdue University, 525 Northwestern Avenue, West Lafayette, Indiana 47907, United States
- Laboratory of Chemical Physics, National Institute of Diabetes, Digestion and Kidney Diseases, National Institutes of Health, 5 Memorial Drive, Bethesda, Maryland 20892, United States
| | - Sergei Savikhin
- Department of Physics and Astronomy, Purdue University, 525 Northwestern Avenue, West Lafayette, Indiana 47907, United States
| | - Lyudmila V Slipchenko
- Department of Chemistry, Purdue University, 560 Oval Drive, West Lafayette, Indiana 47907, United States
| |
Collapse
|
33
|
Kaliakin DS, Nakata H, Kim Y, Chen Q, Fedorov DG, Slipchenko LV. FMOxFMO: Elucidating Excitonic Interactions in the Fenna-Matthews-Olson Complex with the Fragment Molecular Orbital Method. J Chem Theory Comput 2020; 16:1175-1187. [PMID: 31841349 DOI: 10.1021/acs.jctc.9b00621] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
In order to study Förster resonance energy transfer (FRET), the fragment molecular orbital (FMO) method is extended to compute electronic couplings between local excitations via the excited state transition density model, enabling efficient calculations of nonlocal excitations in a large molecular system and overcoming the previous limitation of being able to compute only local excitations. The results of these simple but accurate models are validated against full quantum calculations without fragmentation. The developed method is applied to a very important photosynthetic pigment-protein complex, the Fenna-Matthews-Olson complex (FMOc), that is responsible for the energy transfer from a chlorosome to the reaction center in the green sulfur bacteria. Absorption and circular dichroism spectra of FMOc are simulated, and the role of the molecular environment on the excitations is revealed.
Collapse
Affiliation(s)
- Danil S Kaliakin
- Department of Chemistry , Purdue University , 560 Oval Drive , West Lafayette , Indiana 47907 , United States
| | - Hiroya Nakata
- Research Institute for Advanced Materials and Devices , Kyocera , 5-3 Hikaridai-3 , Seika-cho Soraku-gun, Kyoto 619-0237 , Japan
| | - Yongbin Kim
- Department of Chemistry , Purdue University , 560 Oval Drive , West Lafayette , Indiana 47907 , United States
| | - Qifeng Chen
- Department of Chemistry , Purdue University , 560 Oval Drive , West Lafayette , Indiana 47907 , United States
| | - Dmitri G Fedorov
- Research Center for Computational Design of Advanced Functional Materials (CD-FMat) , National Institute of Advanced Industrial Science and Technology (AIST) , Central 2, Umezono 1-1-1 , Tsukuba 305-8568 , Japan
| | - Lyudmila V Slipchenko
- Department of Chemistry , Purdue University , 560 Oval Drive , West Lafayette , Indiana 47907 , United States
| |
Collapse
|
34
|
Bold BM, Sokolov M, Maity S, Wanko M, Dohmen PM, Kranz JJ, Kleinekathöfer U, Höfener S, Elstner M. Benchmark and performance of long-range corrected time-dependent density functional tight binding (LC-TD-DFTB) on rhodopsins and light-harvesting complexes. Phys Chem Chem Phys 2020; 22:10500-10518. [PMID: 31950960 DOI: 10.1039/c9cp05753f] [Citation(s) in RCA: 40] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
The chromophores of rhodopsins (Rh) and light-harvesting (LH) complexes still represent a major challenge for a quantum chemical description due to their size and complex electronic structure. Since gradient corrected and hybrid density functional approaches have been shown to fail for these systems, only range-separated functionals seem to be a promising alternative to the more time consuming post-Hartree-Fock approaches. For extended sampling of optical properties, however, even more approximate approaches are required. Recently, a long-range corrected (LC) functional has been implemented into the efficient density functional tight binding (DFTB) method, allowing to sample the excited states properties of chromophores embedded into proteins using quantum mechanical/molecular mechanical (QM/MM) with the time-dependent (TD) DFTB approach. In the present study, we assess the accuracy of LC-TD-DFT and LC-TD-DFTB for rhodopsins (bacteriorhodopsin (bR) and pharaonis phoborhodopsin (ppR)) and LH complexes (light-harvesting complex II (LH2) and Fenna-Matthews-Olson (FMO) complex). This benchmark study shows the improved description of the color tuning parameters compared to standard DFT functionals. In general, LC-TD-DFTB can exhibit a similar performance as the corresponding LC functionals, allowing a reliable description of excited states properties at significantly reduced cost. The two chromophores investigated here pose complementary challenges: while huge sensitivity to external field perturbation (color tuning) and charge transfer excitations are characteristic for the retinal chromophore, the multi-chromophoric character of the LH complexes emphasizes a correct description of inter-chromophore couplings, giving less importance to color tuning. None of the investigated functionals masters both systems simultaneously with satisfactory accuracy. LC-TD-DFTB, at the current stage, although showing a systematic improvement compared to TD-DFTB cannot be recommended for studying color tuning in retinal proteins, similar to some of the LC-DFT functionals, because the response to external fields is still too weak. For sampling of LH-spectra, however, LC-TD-DFTB is a viable tool, allowing to efficiently sample absorption energies, as shown for three different LH complexes. As the calculations indicate, geometry optimization may overestimate the importance of local minima, which may be averaged over when using trajectories. Fast quantum chemical approaches therefore may allow for a direct sampling of spectra in the near future.
Collapse
Affiliation(s)
- Beatrix M Bold
- Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT), Kaiserstrasse 12, 76131 Karlsruhe, Germany.
| | | | | | | | | | | | | | | | | |
Collapse
|
35
|
Cerezo J, Aranda D, Avila Ferrer FJ, Prampolini G, Santoro F. Adiabatic-Molecular Dynamics Generalized Vertical Hessian Approach: A Mixed Quantum Classical Method To Compute Electronic Spectra of Flexible Molecules in the Condensed Phase. J Chem Theory Comput 2020; 16:1215-1231. [PMID: 31855424 DOI: 10.1021/acs.jctc.9b01009] [Citation(s) in RCA: 45] [Impact Index Per Article: 11.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
We present a general mixed quantum classical method that couples classical molecular dynamics (MD) and vibronic models to compute the shape of electronic spectra of flexible molecules in the condensed phase without, in principle, any phenomenological broadening. It is based on a partition of the nuclear motions of the solute + solvent system in "soft" and "stiff" vibrational modes and an adiabatic hypothesis that assumes that stiff modes are much faster than soft ones. In this framework, the spectrum is rigorously expressed as a conformational integral of quantum vibronic spectra along the stiff coordinates only. Soft modes enter at the classical level through the conformational distribution that is sampled with classical MD runs. In each configuration, reduced-dimensionality quadratic Hamiltonians are built in the space of the stiff coordinates only, thanks to a generalization of the Vertical Hessian harmonic model and an iterative application of projectors in internal coordinates to remove soft modes. Quantum vibronic spectra, specific for each sampled configuration of the soft coordinates, are then computed at the desired temperature with efficient time-dependent techniques, and the global spectrum simply arises from their average. For consistency of the whole procedure, classical MD runs are performed with quantum-mechanically derived force fields, parameterized at the same level of theory selected for generating the quadratic Hamiltonians along the stiff coordinates. Application to N-methyl-6-oxyquinolinium betaine in water, dithiophene in ethanol, and cyanidine in water is presented to show the performance of the method.
Collapse
Affiliation(s)
| | - Daniel Aranda
- CNR-Consiglio Nazionale delle Ricerche, Istituto di Chimica dei Composti Organo Metallici (ICCOM-CNR), SS di Pisa , Area della Ricerca, via G. Moruzzi 1 , I-56124 Pisa , Italy.,Departamento de Química Física, Facultad de Ciencias, Andalucía Tech , Universidad de Málaga , E-29071 Málaga , Spain
| | - Francisco José Avila Ferrer
- Departamento de Química Física, Facultad de Ciencias, Andalucía Tech , Universidad de Málaga , E-29071 Málaga , Spain
| | - Giacomo Prampolini
- CNR-Consiglio Nazionale delle Ricerche, Istituto di Chimica dei Composti Organo Metallici (ICCOM-CNR), SS di Pisa , Area della Ricerca, via G. Moruzzi 1 , I-56124 Pisa , Italy
| | - Fabrizio Santoro
- CNR-Consiglio Nazionale delle Ricerche, Istituto di Chimica dei Composti Organo Metallici (ICCOM-CNR), SS di Pisa , Area della Ricerca, via G. Moruzzi 1 , I-56124 Pisa , Italy
| |
Collapse
|
36
|
Zuehlsdorff TJ, Hong H, Shi L, Isborn CM. Influence of Electronic Polarization on the Spectral Density. J Phys Chem B 2019; 124:531-543. [DOI: 10.1021/acs.jpcb.9b10250] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023]
Affiliation(s)
- Tim J. Zuehlsdorff
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, United States
| | - Hanbo Hong
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, United States
| | - Liang Shi
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, United States
| | - Christine M. Isborn
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, United States
| |
Collapse
|
37
|
Saito S, Higashi M, Fleming GR. Site-Dependent Fluctuations Optimize Electronic Energy Transfer in the Fenna–Matthews–Olson Protein. J Phys Chem B 2019; 123:9762-9772. [DOI: 10.1021/acs.jpcb.9b07456] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Affiliation(s)
- Shinji Saito
- Institute for Molecular Science, Myodaiji, Okazaki, Aichi 444-8585, Japan
- The Graduate University for Advanced Studies, Myodaiji, Okazaki, Aichi 444-8585, Japan
| | - Masahiro Higashi
- Department of Molecular Engineering, Graduate School of Engineering, Kyoto University, Kyoto daigaku-katsura, Nishikyo-ku, Kyoto 615-8510, Kyoto, Japan
- JST, PRESTO, 4-1-8 Honcho, Kawaguchi, Saitama 332-0012, Japan
| | - Graham R. Fleming
- Department of Chemistry, University of California, Berkeley, California 94720, United States
- Molecular Biophysics and Integrated Bioimaging Division, Lawrence Berkeley National Lab, Berkeley, California 94720, United States
- Kavli Energy Nanosciences Institute at Berkeley, Berkeley, California 94720, United States
| |
Collapse
|
38
|
Daskalakis V, Maity S, Hart CL, Stergiannakos T, Duffy CDP, Kleinekathöfer U. Structural Basis for Allosteric Regulation in the Major Antenna Trimer of Photosystem II. J Phys Chem B 2019; 123:9609-9615. [DOI: 10.1021/acs.jpcb.9b09767] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Affiliation(s)
- Vangelis Daskalakis
- Department of Environmental Science and Technology, Cyprus University of Technology, 30 Archbishop Kyprianou Street, 3603, Limassol, Cyprus
| | - Sayan Maity
- Department of Physics & Earth Sciences, Jacobs University Bremen, Campus Ring 1, 28759 Bremen, Germany
| | - Cameron Lewis Hart
- School of Biological and Chemical Sciences, Queen Mary University of London, Mile End Road, London E1 4NS, United Kingdom
| | - Taxiarchis Stergiannakos
- Department of Environmental Science and Technology, Cyprus University of Technology, 30 Archbishop Kyprianou Street, 3603, Limassol, Cyprus
| | - Christopher D. P. Duffy
- School of Biological and Chemical Sciences, Queen Mary University of London, Mile End Road, London E1 4NS, United Kingdom
| | - Ulrich Kleinekathöfer
- Department of Physics & Earth Sciences, Jacobs University Bremen, Campus Ring 1, 28759 Bremen, Germany
| |
Collapse
|
39
|
Duan HG, Thorwart M, Miller RJD. Does electronic coherence enhance anticorrelated pigment vibrations under realistic conditions? J Chem Phys 2019; 151:114115. [PMID: 31542003 DOI: 10.1063/1.5119248] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The light-harvesting efficiency of a photoactive molecular complex is largely determined by the properties of its electronic quantum states. Those, in turn, are influenced by molecular vibrational states of the nuclear degrees of freedom. Here, we reexamine two recently formulated concepts that a coherent vibronic coupling between molecular states would either extend the electronic coherence lifetime or enhance the amplitude of the anticorrelated vibrational mode at longer times. For this, we study a vibronically coupled dimer and calculate the nonlinear two-dimensional (2D) electronic spectra that directly reveal electronic coherence. The time scale of electronic coherence is initially extracted by measuring the antidiagonal bandwidth of the central peak in the 2D spectrum at zero waiting time. Based on the residual analysis, we identify small-amplitude long-lived oscillations in the cross-peaks, which, however, are solely due to groundstate vibrational coherence, regardless of having resonant or off-resonant conditions. Our studies neither show an enhancement of the electronic quantum coherence nor an enhancement of the anticorrelated vibrational mode by the vibronic coupling under ambient conditions.
Collapse
Affiliation(s)
- Hong-Guang Duan
- Max Planck Institute for the Structure and Dynamics of Matter, Luruper Chaussee 149, 22761 Hamburg, Germany
| | - Michael Thorwart
- I. Institut für Theoretische Physik, Universität Hamburg, Jungiusstraße 9, 20355 Hamburg, Germany
| | - R J Dwayne Miller
- Max Planck Institute for the Structure and Dynamics of Matter, Luruper Chaussee 149, 22761 Hamburg, Germany
| |
Collapse
|
40
|
Tamascelli D, Smirne A, Lim J, Huelga SF, Plenio MB. Efficient Simulation of Finite-Temperature Open Quantum Systems. PHYSICAL REVIEW LETTERS 2019; 123:090402. [PMID: 31524443 DOI: 10.1103/physrevlett.123.090402] [Citation(s) in RCA: 39] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/19/2018] [Indexed: 06/10/2023]
Abstract
Chain-mapping techniques in combination with the time-dependent density matrix renormalization group are a powerful tool for the simulation of open-system quantum dynamics. For finite-temperature environments, however, this approach suffers from an unfavorable algorithmic scaling with increasing temperature. We prove that the system dynamics under thermal environments can be nonperturbatively described by temperature-dependent system-environmental couplings with the initial environment state being in its pure vacuum state, instead of a mixed thermal state. As a consequence, as long as the initial system state is pure, the global system-environment state remains pure at all times. The resulting speed-up and relaxed memory requirements of this approach enable the efficient simulation of open quantum systems interacting with highly structured environments in any temperature range, with applications extending from quantum thermodynamics to quantum effects in mesoscopic systems.
Collapse
Affiliation(s)
- D Tamascelli
- Institut für Theoretische Physik, Albert-Einstein-Allee 11, Universität Ulm, 89069 Ulm, Germany
- Dipartimento di Fisica "Aldo Pontremoli", Università degli Studi di Milano, via Celoria 16, 20133 Milano, Italy
| | - A Smirne
- Institut für Theoretische Physik, Albert-Einstein-Allee 11, Universität Ulm, 89069 Ulm, Germany
| | - J Lim
- Institut für Theoretische Physik, Albert-Einstein-Allee 11, Universität Ulm, 89069 Ulm, Germany
| | - S F Huelga
- Institut für Theoretische Physik, Albert-Einstein-Allee 11, Universität Ulm, 89069 Ulm, Germany
| | - M B Plenio
- Institut für Theoretische Physik, Albert-Einstein-Allee 11, Universität Ulm, 89069 Ulm, Germany
| |
Collapse
|
41
|
Zuehlsdorff TJ, Montoya-Castillo A, Napoli JA, Markland TE, Isborn CM. Optical spectra in the condensed phase: Capturing anharmonic and vibronic features using dynamic and static approaches. J Chem Phys 2019; 151:074111. [PMID: 31438704 DOI: 10.1063/1.5114818] [Citation(s) in RCA: 54] [Impact Index Per Article: 10.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023] Open
Abstract
Simulating optical spectra in the condensed phase remains a challenge for theory due to the need to capture spectral signatures arising from anharmonicity and dynamical effects, such as vibronic progressions and asymmetry. As such, numerous simulation methods have been developed that invoke different approximations and vary in their ability to capture different physical regimes. Here, we use several models of chromophores in the condensed phase and ab initio molecular dynamics simulations to rigorously assess the applicability of methods to simulate optical absorption spectra. Specifically, we focus on the ensemble scheme, which can address anharmonic potential energy surfaces but relies on the applicability of extreme nuclear-electronic time scale separation; the Franck-Condon method, which includes dynamical effects but generally only at the harmonic level; and the recently introduced ensemble zero-temperature Franck-Condon approach, which straddles these limits. We also devote particular attention to the performance of methods derived from a cumulant expansion of the energy gap fluctuations and test the ability to approximate the requisite time correlation functions using classical dynamics with quantum correction factors. These results provide insights as to when these methods are applicable and able to capture the features of condensed phase spectra qualitatively and, in some cases, quantitatively across a range of regimes.
Collapse
Affiliation(s)
- Tim J Zuehlsdorff
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, USA
| | | | - Joseph A Napoli
- Department of Chemistry, Stanford University, Stanford, California 94305, USA
| | - Thomas E Markland
- Department of Chemistry, Stanford University, Stanford, California 94305, USA
| | - Christine M Isborn
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, USA
| |
Collapse
|
42
|
Cupellini L, Bondanza M, Nottoli M, Mennucci B. Successes & challenges in the atomistic modeling of light-harvesting and its photoregulation. BIOCHIMICA ET BIOPHYSICA ACTA-BIOENERGETICS 2019; 1861:148049. [PMID: 31386831 DOI: 10.1016/j.bbabio.2019.07.004] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/10/2019] [Revised: 07/10/2019] [Accepted: 07/15/2019] [Indexed: 12/21/2022]
Abstract
Light-harvesting is a crucial step of photosynthesis. Its mechanisms and related energetics have been revealed by a combination of experimental investigations and theoretical modeling. The success of theoretical modeling is largely due to the application of atomistic descriptions combining quantum chemistry, classical models and molecular dynamics techniques. Besides the important achievements obtained so far, a complete and quantitative understanding of how the many different light-harvesting complexes exploit their structural specificity is still missing. Moreover, many questions remain unanswered regarding the mechanisms through which light-harvesting is regulated in response to variable light conditions. Here we show that, in both fields, a major role will be played once more by atomistic descriptions, possibly generalized to tackle the numerous time and space scales on which the regulation takes place: going from the ultrafast electronic excitation of the multichromophoric aggregate, through the subsequent conformational changes in the embedding protein, up to the interaction between proteins.
Collapse
Affiliation(s)
- Lorenzo Cupellini
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via G. Moruzzi 13, Pisa 56124, Italy
| | - Mattia Bondanza
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via G. Moruzzi 13, Pisa 56124, Italy
| | - Michele Nottoli
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via G. Moruzzi 13, Pisa 56124, Italy
| | - Benedetta Mennucci
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via G. Moruzzi 13, Pisa 56124, Italy.
| |
Collapse
|
43
|
Rahman H, Kleinekathöfer U. Chebyshev hierarchical equations of motion for systems with arbitrary spectral densities and temperatures. J Chem Phys 2019; 150:244104. [PMID: 31255062 DOI: 10.1063/1.5100102] [Citation(s) in RCA: 21] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023] Open
Abstract
The time evolution in open quantum systems, such as a molecular aggregate in contact with a thermal bath, still poses a complex and challenging problem. The influence of the thermal noise can be treated using a plethora of schemes, several of which decompose the corresponding correlation functions in terms of weighted sums of exponential functions. One such scheme is based on the hierarchical equations of motion (HEOM), which is built using only certain forms of bath correlation functions. In the case where the environment is described by a complex spectral density or is at a very low temperature, approaches utilizing the exponential decomposition become very inefficient. Here, we utilize an alternative decomposition scheme for the bath correlation function based on Chebyshev polynomials and Bessel functions to derive a HEOM approach up to an arbitrary order in the environmental coupling. These hierarchical equations are similar in structure to the popular exponential HEOM scheme, but are formulated using the derivatives of the Bessel functions. The proposed scheme is tested up to the fourth order in perturbation theory for a two-level system and compared to benchmark calculations for the case of zero-temperature quantum Ohmic and super-Ohmic noise. Furthermore, the benefits and shortcomings of the present Chebyshev-based hierarchical equations are discussed.
Collapse
Affiliation(s)
- Hasan Rahman
- Department of Physics and Earth Sciences, Jacobs University Bremen, Campus Ring 1, 28759 Bremen, Germany
| | - Ulrich Kleinekathöfer
- Department of Physics and Earth Sciences, Jacobs University Bremen, Campus Ring 1, 28759 Bremen, Germany
| |
Collapse
|
44
|
Teh HH, Jin BY, Cheng YC. Frozen-mode small polaron quantum master equation with variational bound for excitation energy transfer in molecular aggregates. J Chem Phys 2019; 150:224110. [DOI: 10.1063/1.5096287] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Hung-Hsuan Teh
- Department of Chemistry, National Taiwan University, Taipei City 10617, Taiwan
| | - Bih-Yaw Jin
- Department of Chemistry, National Taiwan University, Taipei City 10617, Taiwan
| | - Yuan-Chung Cheng
- Department of Chemistry, National Taiwan University, Taipei City 10617, Taiwan
| |
Collapse
|
45
|
|
46
|
Zajac G, Machalska E, Kaczor A, Kessler J, Bouř P, Baranska M. Structure of supramolecular astaxanthin aggregates revealed by molecular dynamics and electronic circular dichroism spectroscopy. Phys Chem Chem Phys 2019; 20:18038-18046. [PMID: 29932184 DOI: 10.1039/c8cp01742e] [Citation(s) in RCA: 23] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
Biomolecular aggregation is omnipresent in nature and important for metabolic processes or in medical treatment; however, the phenomenon is rather difficult to predict or understand on the basis of computational models. Recently, we found that electronic circular dichroism (ECD) spectroscopy and closely related resonance Raman optical activity (RROA) are extremely sensitive to the aggregation mechanism and structure of the astaxanthin dye. In the present study, molecular dynamics (MD) and quantum chemical (QC) computations (ZIndo/S, TDDFT) are used to link the aggregate structure with ECD spectral shapes. Realistic absorption and ECD intensities were obtained and the simulations reproduced many trends observed experimentally, such as the prevalent sign pattern and dependence of the aggregate structure on the solvent type. The computationally cheaper ZIndo/S method provided results very similar to those obtained by TDDFT. In the future, the accuracy of the combined MD/QC methodology of spectra interpretation should be improved to provide more detailed information on astaxanthin aggregates and similar macromolecular systems.
Collapse
Affiliation(s)
- Grzegorz Zajac
- Faculty of Chemistry, Jagiellonian University, Gronostajowa 2, Krakow 30-387, Poland.
| | - Ewa Machalska
- Faculty of Chemistry, Jagiellonian University, Gronostajowa 2, Krakow 30-387, Poland.
| | - Agnieszka Kaczor
- Faculty of Chemistry, Jagiellonian University, Gronostajowa 2, Krakow 30-387, Poland. and Jagiellonian Centre for Experimental Therapeutics (JCET), Jagiellonian University, Bobrzynskiego 14, Krakow 30-348, Poland
| | - Jiří Kessler
- Institute of Organic Chemistry and Biochemistry, Academy of Sciences, Flemingovo náměstí 2, Prague, 16610, Czech Republic.
| | - Petr Bouř
- Institute of Organic Chemistry and Biochemistry, Academy of Sciences, Flemingovo náměstí 2, Prague, 16610, Czech Republic.
| | - Malgorzata Baranska
- Faculty of Chemistry, Jagiellonian University, Gronostajowa 2, Krakow 30-387, Poland. and Jagiellonian Centre for Experimental Therapeutics (JCET), Jagiellonian University, Bobrzynskiego 14, Krakow 30-348, Poland
| |
Collapse
|
47
|
Claridge K, Troisi A. Developing Consistent Molecular Dynamics Force Fields for Biological Chromophores via Force Matching. J Phys Chem B 2018; 123:428-438. [DOI: 10.1021/acs.jpcb.8b10746] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Kirsten Claridge
- Department of Chemistry, University of Liverpool, Liverpool L69 7ZD, U.K
| | - Alessandro Troisi
- Department of Chemistry, University of Liverpool, Liverpool L69 7ZD, U.K
| |
Collapse
|
48
|
Mallus MI, Shakya Y, Prajapati JD, Kleinekathöfer U. Environmental effects on the dynamics in the light-harvesting complexes LH2 and LH3 based on molecular simulations. Chem Phys 2018. [DOI: 10.1016/j.chemphys.2018.08.013] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
|
49
|
Cupellini L, Corbella M, Mennucci B, Curutchet C. Electronic energy transfer in biomacromolecules. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2018. [DOI: 10.1002/wcms.1392] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Affiliation(s)
- Lorenzo Cupellini
- Dipartimento di Chimica e Chimica Industriale University of Pisa Pisa Italy
| | - Marina Corbella
- Departament de Farmàcia i Tecnologia Farmacèutica i Fisicoquímica and Institut de Química Teòrica i Computacional (IQTC‐UB), Facultat de Farmàcia i Ciències de l'Alimentació Universitat de Barcelona Barcelona Spain
| | - Benedetta Mennucci
- Dipartimento di Chimica e Chimica Industriale University of Pisa Pisa Italy
| | - Carles Curutchet
- Departament de Farmàcia i Tecnologia Farmacèutica i Fisicoquímica and Institut de Química Teòrica i Computacional (IQTC‐UB), Facultat de Farmàcia i Ciències de l'Alimentació Universitat de Barcelona Barcelona Spain
| |
Collapse
|
50
|
Zhang PP, Bentley CDB, Eisfeld A. Flexible scheme to truncate the hierarchy of pure states. J Chem Phys 2018; 148:134103. [DOI: 10.1063/1.5022225] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- P.-P. Zhang
- Max Planck Institute for the Physics of Complex Systems, Nöthnitzer Strasse 38, D-01187 Dresden, Germany
| | - C. D. B. Bentley
- Max Planck Institute for the Physics of Complex Systems, Nöthnitzer Strasse 38, D-01187 Dresden, Germany
| | - A. Eisfeld
- Max Planck Institute for the Physics of Complex Systems, Nöthnitzer Strasse 38, D-01187 Dresden, Germany
| |
Collapse
|