1
|
Cupellini L, Gwizdala M, Krüger TPJ. Energetic Landscape and Terminal Emitters of Phycobilisome Cores from Quantum Chemical Modeling. J Phys Chem Lett 2024; 15:9746-9756. [PMID: 39288324 DOI: 10.1021/acs.jpclett.4c02005] [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: 09/19/2024]
Abstract
Phycobilisomes (PBs) are giant antenna supercomplexes of cyanobacteria that use phycobilin pigments to capture sunlight and transfer the collected energy to membrane-bound photosystems. In the PB core, phycobilins are bound to particular allophycocyanin (APC) proteins. Some phycobilins are thought to be terminal emitters (TEs) with red-shifted fluorescence. However, the precise identification of TEs is still under debate. In this work, we employ multiscale quantum-mechanical calculations to disentangle the excitation energy landscape of PB cores. Using the recent atomistic PB structures from Synechoccoccus PCC 7002 and Synechocystis PCC 6803, we compute the spectral properties of different APC trimers and assign the low-energy pigments. We show that the excitation energy of APC phycobilins is determined by geometric and electrostatic factors and is tuned by the specific protein-protein interactions within the core. Our findings challenge the simple picture of a few red-shifted bilins in the PB core and instead suggest that the red-shifts are established by the entire TE-containing APC trimers. Our work provides a theoretical microscopic basis for the interpretation of energy migration and time-resolved spectroscopy in phycobilisomes.
Collapse
Affiliation(s)
- Lorenzo Cupellini
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via G. Moruzzi 13, 56124 Pisa, Italy
| | - Michal Gwizdala
- Department of Physics, University of Pretoria, Lynnwood Road, Pretoria 0002, South Africa
- Forestry and Agricultural Biotechnology Institute (FABI), University of Pretoria, Lynnwood Road, Pretoria 0002, South Africa
| | - Tjaart P J Krüger
- Department of Physics, University of Pretoria, Lynnwood Road, Pretoria 0002, South Africa
- Forestry and Agricultural Biotechnology Institute (FABI), University of Pretoria, Lynnwood Road, Pretoria 0002, South Africa
- National Institute of Theoretical and Computational Sciences (NITheCS), https://nithecs.ac.za/
| |
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
|
Kim CW, Franco I. General framework for quantifying dissipation pathways in open quantum systems. II. Numerical validation and the role of non-Markovianity. J Chem Phys 2024; 160:214112. [PMID: 38833365 DOI: 10.1063/5.0202862] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/07/2024] [Accepted: 05/13/2024] [Indexed: 06/06/2024] Open
Abstract
In the previous paper [C. W. Kim and I. Franco, J. Chem. Phys. 160, 214111-1-214111-13 (2024)], we developed a theory called MQME-D, which allows us to decompose the overall energy dissipation process in open quantum system dynamics into contributions by individual components of the bath when the subsystem dynamics is governed by a Markovian quantum master equation (MQME). Here, we contrast the predictions of MQME-D against the numerically exact results obtained by combining hierarchical equations of motion (HEOM) with a recently reported protocol for monitoring the statistics of the bath. Overall, MQME-D accurately captures the contributions of specific bath components to the overall dissipation while greatly reducing the computational cost compared to exact computations using HEOM. The computations show that MQME-D exhibits errors originating from its inherent Markov approximation. We demonstrate that its accuracy can be significantly increased by incorporating non-Markovianity by exploiting time scale separations (TSS) in different components of the bath. Our work demonstrates that MQME-D combined with TSS can be reliably used to understand how energy is dissipated in realistic open quantum system dynamics.
Collapse
Affiliation(s)
- Chang Woo Kim
- Department of Chemistry, Chonnam National University, Gwangju 61186, South Korea
| | - Ignacio Franco
- Department of Chemistry, University of Rochester, Rochester, New York 14627, USA
- Department of Physics, University of Rochester, Rochester, New York 14627, USA
| |
Collapse
|
4
|
Lorenzoni N, Cho N, Lim J, Tamascelli D, Huelga SF, Plenio MB. Systematic Coarse Graining of Environments for the Nonperturbative Simulation of Open Quantum Systems. PHYSICAL REVIEW LETTERS 2024; 132:100403. [PMID: 38518302 DOI: 10.1103/physrevlett.132.100403] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/23/2023] [Revised: 01/22/2024] [Accepted: 02/13/2024] [Indexed: 03/24/2024]
Abstract
Conducting precise electronic-vibrational dynamics simulations of molecular systems poses significant challenges when dealing with realistic environments composed of numerous vibrational modes. Here, we introduce a technique for the construction of effective phonon spectral densities that capture accurately open-system dynamics over a finite time interval of interest. When combined with existing nonperturbative simulation tools, our approach can reduce significantly the computational costs associated with many-body open-system dynamics.
Collapse
Affiliation(s)
- Nicola Lorenzoni
- Institut für Theoretische Physik und IQST, Albert-Einstein-Allee 11, Universität Ulm, D-89081 Ulm, Germany
| | - Namgee Cho
- Institut für Theoretische Physik und IQST, Albert-Einstein-Allee 11, Universität Ulm, D-89081 Ulm, Germany
| | - James Lim
- Institut für Theoretische Physik und IQST, Albert-Einstein-Allee 11, Universität Ulm, D-89081 Ulm, Germany
| | - Dario Tamascelli
- Institut für Theoretische Physik und IQST, Albert-Einstein-Allee 11, Universität Ulm, D-89081 Ulm, Germany
- Dipartimento di Fisica "Aldo Pontremoli," Università degli Studi di Milano, Via Celoria 16, 20133 Milano, Italy
| | - Susana F Huelga
- Institut für Theoretische Physik und IQST, Albert-Einstein-Allee 11, Universität Ulm, D-89081 Ulm, Germany
| | - Martin B Plenio
- Institut für Theoretische Physik und IQST, Albert-Einstein-Allee 11, Universität Ulm, D-89081 Ulm, Germany
| |
Collapse
|
5
|
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
|
6
|
Liu R, Zhen ZH, Li W, Ge B, Qin S. How can Phycobilisome, the unique light harvesting system in certain algae working highly efficiently: The connection in between structures and functions. PROGRESS IN BIOPHYSICS AND MOLECULAR BIOLOGY 2024; 186:39-52. [PMID: 38030044 DOI: 10.1016/j.pbiomolbio.2023.11.005] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/11/2023] [Revised: 11/02/2023] [Accepted: 11/23/2023] [Indexed: 12/01/2023]
Abstract
Algae, which are ubiquitous in ecosystems, have evolved a variety of light-harvesting complexes to better adapt to diverse habitats. Phycobilisomes/phycobiliproteins, unique to cyanobacteria, red algae, and certain cryptomonads, compensate for the lack of chlorophyll absorption, allowing algae to capture and efficiently transfer light energy in aquatic environments. With the advancement of microscopy and spectroscopy, the structure and energy transfer processes of increasingly complex phycobilisomes have been elucidated, providing us with a vivid portrait of the dynamic adaptation of their structures to the light environment in which algae thrive: 1) Cyanobacteria living on the surface of the water use short, small phycobilisomes to absorb red-orange light and reduce the damage from blue-violet light via multiple methods; 2) Large red algae inhabiting the depths of the ocean have evolved long and dense phycobilisomes containing phycoerythrin to capture the feeble blue-green light; 3) In far-red light environments such as caves, algae use special allophycocyanin cores to optimally utilize the far-red light; 4) When the environment shifts, algae can adjust the length, composition and density of their rods to better adapt; 5) By carefully designing the position of the pigments, phycobilisomes can transfer light energy to the reaction center with nearly 100% efficiency via three energy transfer processes.
Collapse
Affiliation(s)
- Runze Liu
- Yantai Institute of Coastal Zone Research, Chinese Academy of Sciences, Yantai, Shandong, 264003, China; University of Chinese Academy of Sciences, Beijing, 100000, China
| | - Zhang-He Zhen
- Institute of Physics, Chinese Academy of Sciences, Beijing, 100190, China
| | - Wenjun Li
- Yantai Institute of Coastal Zone Research, Chinese Academy of Sciences, Yantai, Shandong, 264003, China
| | - Baosheng Ge
- China University of Petroleum (HUADONG), Qingdao, Shandong, 266580, China
| | - Song Qin
- Yantai Institute of Coastal Zone Research, Chinese Academy of Sciences, Yantai, Shandong, 264003, China.
| |
Collapse
|
7
|
Li Y, Li Q, Gillilan RE, Abbaspourrad A. Reversible disassembly-reassembly of C-phycocyanin in pressurization-depressurization cycles of high hydrostatic pressure. Int J Biol Macromol 2023; 253:127623. [PMID: 37879586 PMCID: PMC10842036 DOI: 10.1016/j.ijbiomac.2023.127623] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/12/2023] [Revised: 10/20/2023] [Accepted: 10/21/2023] [Indexed: 10/27/2023]
Abstract
Hydrostatic pressure can reversibly modulate protein-protein and protein-chromophore interactions of C-phycocyanin (C-PC) from Spirulina platensis. Small-angle X-ray scattering combined with UV-Vis spectrophotometry and protein modeling was used to explore the color and structural changes of C-PC under high pressure conditions at different pH levels. It was revealed that pressures up to 350 MPa were enough to fully disassemble C-PC from trimers to monomers at pH 7.0, or from monomers to detached subunits at pH 9.0. These disassemblies were accompanied by protein unfolding that caused these high-pressure induced structures to be more extended. These changes were reversible following depressurization. The trimer-to-monomer transition proceeded through a collection of previously unrecognized, L-shaped intermediates resembling C-PC dimers. Additionally, pressurized C-PC showed decayed Q-band absorption and fortified Soret-band absorption. This was evidence that the folded tetrapyrroles, which had folded at ambient pressure, formed semicyclic unfolded conformations at a high pressure. Upon depressurization, the peak intensity and shift all recovered stepwise, showing pressure can precisely manipulate C-PC's structure as well as its color. Overall, a protein-chromophore regulatory theory of C-PC was unveiled. The pressure-tunability could be harnessed to modify and stabilize C-PC's structure and photochemical properties for designing new delivery and optical materials.
Collapse
Affiliation(s)
- Ying Li
- Department of Food Science, College of Agriculture and Life Sciences, Cornell University, Ithaca, NY, USA
| | - Qike Li
- Department of Food Science, College of Agriculture and Life Sciences, Cornell University, Ithaca, NY, USA
| | - Richard E Gillilan
- Cornell High Energy Synchrotron Source (MacCHESS), Cornell University, Ithaca, NY, USA
| | - Alireza Abbaspourrad
- Department of Food Science, College of Agriculture and Life Sciences, Cornell University, Ithaca, NY, USA.
| |
Collapse
|
8
|
Gustin I, Kim CW, McCamant DW, Franco I. Mapping electronic decoherence pathways in molecules. Proc Natl Acad Sci U S A 2023; 120:e2309987120. [PMID: 38015846 PMCID: PMC10710033 DOI: 10.1073/pnas.2309987120] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/13/2023] [Accepted: 10/25/2023] [Indexed: 11/30/2023] Open
Abstract
Establishing the fundamental chemical principles that govern molecular electronic quantum decoherence has remained an outstanding challenge. Fundamental questions such as how solvent and intramolecular vibrations or chemical functionalization contribute to the decoherence remain unanswered and are beyond the reach of state-of-the-art theoretical and experimental approaches. Here we address this challenge by developing a strategy to isolate electronic decoherence pathways for molecular chromophores immersed in condensed phase environments that enables elucidating how electronic quantum coherence is lost. For this, we first identify resonance Raman spectroscopy as a general experimental method to reconstruct molecular spectral densities with full chemical complexity at room temperature, in solvent, and for fluorescent and non-fluorescent molecules. We then show how to quantitatively capture the decoherence dynamics from the spectral density and identify decoherence pathways by decomposing the overall coherence loss into contributions due to individual molecular vibrations and solvent modes. We illustrate the utility of the strategy by analyzing the electronic decoherence pathways of the DNA base thymine in water. Its electronic coherences decay in [Formula: see text]30 fs. The early-time decoherence is determined by intramolecular vibrations while the overall decay by solvent. Chemical substitution of thymine modulates the decoherence with hydrogen-bond interactions of the thymine ring with water leading to the fastest decoherence. Increasing temperature leads to faster decoherence as it enhances the importance of solvent contributions but leaves the early-time decoherence dynamics intact. The developed strategy opens key opportunities to establish the connection between molecular structure and quantum decoherence as needed to develop chemical strategies to rationally modulate it.
Collapse
Affiliation(s)
- Ignacio Gustin
- Department of Chemistry, University of Rochester, Rochester, NY14627
| | - Chang Woo Kim
- Department of Chemistry, Chonnam National University, Gwangju61186, South Korea
| | - David W. McCamant
- Department of Chemistry, University of Rochester, Rochester, NY14627
| | - Ignacio Franco
- Department of Chemistry, University of Rochester, Rochester, NY14627
- Department of Physics, University of Rochester, Rochester, NY14627
| |
Collapse
|
9
|
Liu W, Chen ZH, Su Y, Wang Y, Dou W. Predicting rate kernels via dynamic mode decomposition. J Chem Phys 2023; 159:144110. [PMID: 37823462 DOI: 10.1063/5.0170512] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/03/2023] [Accepted: 09/25/2023] [Indexed: 10/13/2023] Open
Abstract
Simulating dynamics of open quantum systems is sometimes a significant challenge, despite the availability of various exact or approximate methods. Particularly when dealing with complex systems, the huge computational cost will largely limit the applicability of these methods. In this work, we investigate the usage of dynamic mode decomposition (DMD) to evaluate the rate kernels in quantum rate processes. DMD is a data-driven model reduction technique that characterizes the rate kernels using snapshots collected from a small time window, allowing us to predict the long-term behaviors with only a limited number of samples. Our investigations show that whether the external field is involved or not, the DMD can give accurate prediction of the result compared with the traditional propagations, and simultaneously reduce the required computational cost.
Collapse
Affiliation(s)
- Wei Liu
- Department of Chemistry, School of Science, Westlake University, Hangzhou 310024 Zhejiang, China
- Institute of Natural Sciences, Westlake Institute for Advanced Study, Hangzhou 310024 Zhejiang, China
| | - Zi-Hao Chen
- Department of Chemical Physics, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Yu Su
- Department of Chemical Physics, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Yao Wang
- Department of Chemical Physics, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Wenjie Dou
- Department of Chemistry, School of Science, Westlake University, Hangzhou 310024 Zhejiang, China
- Institute of Natural Sciences, Westlake Institute for Advanced Study, Hangzhou 310024 Zhejiang, China
- Department of Physics, School of Science, Westlake University, Hangzhou 310024, Zhejiang, China
| |
Collapse
|
10
|
Novoderezhkin VI. Resonant vibrations produce quantum bridge over high-energy states in heterogeneous antenna. PHOTOSYNTHESIS RESEARCH 2023; 158:13-21. [PMID: 37584896 DOI: 10.1007/s11120-023-01042-w] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/17/2023] [Accepted: 07/29/2023] [Indexed: 08/17/2023]
Abstract
Photosynthetic light-harvesting complexes usually contain several pools of molecules with a big difference in transition energies, for example, chlorophylls a and b in plant antennas. Some pathways of the excitation energy transfer may include pigments from the low-energy pool separated by a site occupied by a high-energy molecule. We demonstrate that such pathways may be functional if high-frequency intramolecular vibrations fall in resonance with the energy gap between the neighboring molecules belonging to different pools. In this case, a vibration-assisted mixing of the excited states can produce delocalized vibronic states playing a role of 'quantum bridge' that facilitates a passage over high-energy barrier. We perform calculations of the excitation dynamics in the model three-state system with the parameters emerging from our previous studies of real antennas. Simulation of the dynamics in an explicit electron-vibrational basis demonstrates that the rate of transfer between the two chlorophylls a through the chlorophyll b intermediate is increased by a factor of 1.7-2 in the presence of resonant vibration. A possible influence of energetic disorder and other (non-resonant) vibrations on this effect is discussed.
Collapse
Affiliation(s)
- Vladimir I Novoderezhkin
- A. N. Belozersky Institute of Physico-Chemical Biology, Lomonosov Moscow State University, Leninskie Gory, 119992, Moscow, Russia.
| |
Collapse
|
11
|
Ozaydin B, Curutchet C. Unraveling the role of thermal fluctuations on the exciton structure of the cryptophyte PC612 and PC645 photosynthetic antenna complexes. Front Mol Biosci 2023; 10:1268278. [PMID: 37790875 PMCID: PMC10544999 DOI: 10.3389/fmolb.2023.1268278] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/28/2023] [Accepted: 08/30/2023] [Indexed: 10/05/2023] Open
Abstract
Protein scaffolds play a crucial role in tuning the light harvesting properties of photosynthetic pigment-protein complexes, influencing pigment-protein and pigment-pigment excitonic interactions. Here, we investigate the influence of thermal dynamic effects on the protein tuning mechanisms of phycocyanin PC645 and PC612 antenna complexes of cryptophyte algae, featuring closed or open quaternary structures. We employ a dual molecular dynamics (MD) strategy that combines extensive classical MD simulations with multiple short Born-Oppenheimer quantum/molecular mechanical (QM/MM) simulations to accurately account for both static and dynamic disorder effects. Additionally, we compare the results with an alternative protocol based on multiple QM/MM geometry optimizations of the pigments. Subsequently, we employ polarizable QM/MM calculations using time-dependent density functional theory (TD-DFT) to compute the excited states, and we adopt the full cumulant expansion (FCE) formalism to describe the absorption and circular dichroism spectra. Our findings indicate that thermal effects have only minor impacts on the energy ladder in PC612, despite its remarkable flexibility owing to an open quaternary structure. In striking contrast, thermal effects significantly influence the properties of PC645 due to the absence of a hydrogen bond controlling the twist of ring D in PCB β82 bilins, as well as the larger impact of fluctuations on the excited states of MBV pigments, which possess a higher conjugation length compared to other bilin types. Overall, the dual MD protocol combined with the FCE formalism yields excellent spectral properties for PC612 and PC645, and the resultant excitonic Hamiltonians pave the way for future investigations concerning the implications of open and closed quaternary structures on phycocyanin light harvesting properties.
Collapse
Affiliation(s)
- Beste Ozaydin
- Departament de Farmàcia i Tecnologia Farmacèutica, i Fisicoquímica, Facultat de Farmàcia i Ciències de l’Alimentació, Universitat de Barcelona (UB), Barcelona, Spain
- Institut de Química Teòrica i Computacional (IQTCUB), Universitat de Barcelona (UB), Barcelona, Spain
| | - Carles Curutchet
- Departament de Farmàcia i Tecnologia Farmacèutica, i Fisicoquímica, Facultat de Farmàcia i Ciències de l’Alimentació, Universitat de Barcelona (UB), Barcelona, Spain
- Institut de Química Teòrica i Computacional (IQTCUB), Universitat de Barcelona (UB), Barcelona, Spain
| |
Collapse
|
12
|
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
|
13
|
Roy PP, Leonardo C, Orcutt K, Oberg C, Scholes GD, Fleming GR. Infrared Signatures of Phycobilins within the Phycocyanin 645 Complex. J Phys Chem B 2023; 127:4460-4469. [PMID: 37192324 DOI: 10.1021/acs.jpcb.3c01352] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/18/2023]
Abstract
Aquatic photosynthetic organisms evolved to use a variety of light frequencies to perform photosynthesis. Phycobiliprotein phycocyanin 645 (PC645) is a light-harvesting complex in cryptophyte algae able to transfer the absorbed green solar light to other antennas with over 99% efficiency. The infrared signatures of the phycobilin pigments embedded in PC645 are difficult to access and could provide useful information to understand the mechanism behind the high efficiency of energy transfer in PC645. We use visible-pump IR-probe and two-dimensional electronic vibrational spectroscopy to study the dynamical evolution and assign the fingerprint mid-infrared signatures to each pigment in PC645. Here, we report the pigment-specific vibrational markers that enable us to track the spatial flow of excitation energy between the phycobilin pigment pairs. We speculate that two high-frequency modes (1588 and 1596 cm-1) are involved in the vibronic coupling leading to fast (<ps) and direct energy transfer from the highest to lowest exciton, bypassing the intermediate excitons.
Collapse
Affiliation(s)
- Partha Pratim Roy
- Department of Chemistry, University of California, Berkeley, California 94720, United States
- Molecular Biophysics and Integrated Bioimaging Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Cristina Leonardo
- Department of Chemistry, University of California, Berkeley, California 94720, United States
- Molecular Biophysics and Integrated Bioimaging Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Kaydren Orcutt
- Department of Chemistry, University of California, Berkeley, California 94720, United States
- Molecular Biophysics and Integrated Bioimaging Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Catrina Oberg
- Department of Chemistry, Princeton University, Washington Road, Princeton, New Jersey 08540, United States
| | - Gregory D Scholes
- Department of Chemistry, Princeton University, Washington Road, Princeton, New Jersey 08540, United States
| | - Graham R Fleming
- Department of Chemistry, University of California, Berkeley, California 94720, United States
- Molecular Biophysics and Integrated Bioimaging Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
- Kavli Energy Nanoscience Institute at Berkeley, Berkeley, California 94720, United States
| |
Collapse
|
14
|
Brütting M, Foerster JM, Kümmel S. Understanding Primary Charge Separation in the Heliobacterial Reaction Center. J Phys Chem Lett 2023; 14:3092-3102. [PMID: 36951395 DOI: 10.1021/acs.jpclett.3c00377] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/18/2023]
Abstract
The homodimeric reaction center of heliobacteria retains features of the ancestral reaction center and can thus provide insights into the evolution of photosynthesis. Primary charge separation is expected to proceed in a two-step mechanism along either of the two reaction center branches. We reveal the first charge-separation step from first-principles calculations based on time-dependent density functional theory with an optimally tuned range-separated hybrid and ab initio Born-Oppenheimer molecular dynamics: the electron is most likely localized on the electron transfer cofactor 3 (EC3, OH-chlorophyll a), and the hole on the adjacent EC2. Including substantial parts of the surrounding protein environment into the calculations shows that a distinct structural mechanism is decisive for the relative energetic positioning of the electronic excitations: specific charged amino acids in the vicinity of EC3 lower the energy of charge-transfer excitations and thus facilitate efficient charge separation. These results are discussed considering recent experimental insights.
Collapse
|
15
|
Maity S, Kleinekathöfer U. Recent progress in atomistic modeling of light-harvesting complexes: a mini review. PHOTOSYNTHESIS RESEARCH 2023; 156:147-162. [PMID: 36207489 PMCID: PMC10070314 DOI: 10.1007/s11120-022-00969-w] [Citation(s) in RCA: 11] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 06/30/2022] [Accepted: 09/22/2022] [Indexed: 06/16/2023]
Abstract
In this mini review, we focus on recent advances in the atomistic modeling of biological light-harvesting (LH) complexes. Because of their size and sophisticated electronic structures, multiscale methods are required to investigate the dynamical and spectroscopic properties of such complexes. The excitation energies, in this context also known as site energies, excitonic couplings, and spectral densities are key quantities which usually need to be extracted to be able to determine the exciton dynamics and spectroscopic properties. The recently developed multiscale approach based on the numerically efficient density functional tight-binding framework followed by excited state calculations has been shown to be superior to the scheme based on pure classical molecular dynamics simulations. The enhanced approach, which improves the description of the internal vibrational dynamics of the pigment molecules, yields spectral densities in good agreement with the experimental counterparts for various bacterial and plant LH systems. Here, we provide a brief overview of those results and described the theoretical foundation of the multiscale protocol.
Collapse
Affiliation(s)
- 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
|
16
|
Feighan O, Manby FR, Bourne-Worster S. An efficient protocol for excited states of large biochromophores. J Chem Phys 2023; 158:024107. [PMID: 36641400 DOI: 10.1063/5.0132417] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023] Open
Abstract
Efficient energy transport in photosynthetic antenna is a long-standing source of inspiration for artificial light harvesting materials. However, characterizing the excited states of the constituent chromophores poses a considerable challenge to mainstream quantum chemical and semiempirical excited state methods due to their size and complexity and the accuracy required to describe small but functionally important changes in their properties. In this paper, we explore an alternative approach to calculating the excited states of large biochromophores, exemplified by a specific method for calculating the Qy transition of bacteriochlorophyll a, which we name Chl-xTB. Using a diagonally dominant approximation to the Casida equation and a bespoke parameterization scheme, Chl-xTB can match time-dependent density functional theory's accuracy and semiempirical speed for calculating the potential energy surfaces and absorption spectra of chlorophylls. We demonstrate that Chl-xTB (and other prospective realizations of our protocol) can be integrated into multiscale models, including concurrent excitonic and point-charge embedding frameworks, enabling the analysis of biochromophore networks in a native environment. We exploit this capability to probe the low-frequency spectral densities of excitonic energies and interchromophore interactions in the light harvesting antenna protein LH2 (light harvesting complex 2). The impact of low-frequency protein motion on interchromophore coupling and exciton transport has routinely been ignored due to the prohibitive costs of including it in simulations. Our results provide a more rigorous basis for continued use of this approximation by demonstrating that exciton transition energies are unaffected by low-frequency vibrational coupling to exciton interaction energies.
Collapse
Affiliation(s)
- Oliver Feighan
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom
| | - Frederick R Manby
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom
| | - Susannah Bourne-Worster
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom
| |
Collapse
|
17
|
Exact simulation of pigment-protein complexes unveils vibronic renormalization of electronic parameters in ultrafast spectroscopy. Nat Commun 2022; 13:2912. [PMID: 35614049 PMCID: PMC9133012 DOI: 10.1038/s41467-022-30565-4] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/12/2022] [Accepted: 05/04/2022] [Indexed: 12/04/2022] Open
Abstract
The primary steps of photosynthesis rely on the generation, transport, and trapping of excitons in pigment-protein complexes (PPCs). Generically, PPCs possess highly structured vibrational spectra, combining many discrete intra-pigment modes and a quasi-continuous of protein modes, with vibrational and electronic couplings of comparable strength. The intricacy of the resulting vibronic dynamics poses significant challenges in establishing a quantitative connection between spectroscopic data and underlying microscopic models. Here we show how to address this challenge using numerically exact simulation methods by considering two model systems, namely the water-soluble chlorophyll-binding protein of cauliflower and the special pair of bacterial reaction centers. We demonstrate that the inclusion of the full multi-mode vibronic dynamics in numerical calculations of linear spectra leads to systematic and quantitatively significant corrections to electronic parameter estimation. These multi-mode vibronic effects are shown to be relevant in the longstanding discussion regarding the origin of long-lived oscillations in multidimensional nonlinear spectra. Multimode vibronic mixing in model photosynthetic systems revealed by numerically exact simulations is shown to strongly modify linear and non-linear optical responses and facilitate the persistence of coherent dynamics.
Collapse
|
18
|
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
|
19
|
Absorption and stationary fluorescent spectra of molecular sensors in solution: A computational study. COMPUT THEOR CHEM 2022. [DOI: 10.1016/j.comptc.2022.113596] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
|
20
|
Photosynthesis tunes quantum-mechanical mixing of electronic and vibrational states to steer exciton energy transfer. Proc Natl Acad Sci U S A 2021; 118:2018240118. [PMID: 33688046 DOI: 10.1073/pnas.2018240118] [Citation(s) in RCA: 23] [Impact Index Per Article: 7.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022] Open
Abstract
Photosynthetic species evolved to protect their light-harvesting apparatus from photoxidative damage driven by intracellular redox conditions or environmental conditions. The Fenna-Matthews-Olson (FMO) pigment-protein complex from green sulfur bacteria exhibits redox-dependent quenching behavior partially due to two internal cysteine residues. Here, we show evidence that a photosynthetic complex exploits the quantum mechanics of vibronic mixing to activate an oxidative photoprotective mechanism. We use two-dimensional electronic spectroscopy (2DES) to capture energy transfer dynamics in wild-type and cysteine-deficient FMO mutant proteins under both reducing and oxidizing conditions. Under reducing conditions, we find equal energy transfer through the exciton 4-1 and 4-2-1 pathways because the exciton 4-1 energy gap is vibronically coupled with a bacteriochlorophyll-a vibrational mode. Under oxidizing conditions, however, the resonance of the exciton 4-1 energy gap is detuned from the vibrational mode, causing excitons to preferentially steer through the indirect 4-2-1 pathway to increase the likelihood of exciton quenching. We use a Redfield model to show that the complex achieves this effect by tuning the site III energy via the redox state of its internal cysteine residues. This result shows how pigment-protein complexes exploit the quantum mechanics of vibronic coupling to steer energy transfer.
Collapse
|
21
|
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
|
22
|
Varvelo L, Lynd JK, Bennett DIG. Formally exact simulations of mesoscale exciton dynamics in molecular materials. Chem Sci 2021; 12:9704-9711. [PMID: 34349941 PMCID: PMC8293828 DOI: 10.1039/d1sc01448j] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/11/2021] [Accepted: 05/31/2021] [Indexed: 02/04/2023] Open
Abstract
Excited state carriers, such as excitons, can diffuse on the 100 nm to micron length scale in molecular materials but only delocalize over short length scales due to coupling between electronic and vibrational degrees-of-freedom. Here, we leverage the locality of excitons to adaptively solve the hierarchy of pure states equations (HOPS). We demonstrate that our adaptive HOPS (adHOPS) methodology provides a formally exact and size-invariant (i.e., ) scaling algorithm for simulating mesoscale quantum dynamics. Finally, we provide proof-of-principle calculations for exciton diffusion on linear chains containing up to 1000 molecules. The adaptive hierarchy of pure states (adHOPS) algorithm leverages the locality of excitons in molecular materials to perform formally-exact simulations with size-invariant (i.e., ) scaling, enabling efficient simulations of mesoscale exciton dynamics.![]()
Collapse
Affiliation(s)
- Leonel Varvelo
- Department of Chemistry, Southern Methodist University PO Box 750314 Dallas TX USA
| | - Jacob K Lynd
- Department of Chemistry, Southern Methodist University PO Box 750314 Dallas TX USA
| | - Doran I G Bennett
- Department of Chemistry, Southern Methodist University PO Box 750314 Dallas TX USA
| |
Collapse
|
23
|
Renger T. Semiclassical Modified Redfield and Generalized Förster Theories of Exciton Relaxation/Transfer in Light-Harvesting Complexes: The Quest for the Principle of Detailed Balance. J Phys Chem B 2021; 125:6406-6416. [PMID: 34126008 PMCID: PMC8237266 DOI: 10.1021/acs.jpcb.1c01479] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
![]()
A conceptual problem
of transfer theories that use a semiclassical
description of the electron-vibrational coupling is the neglect of
the correlation between momenta and coordinates of nuclei. In the
Redfield theory of exciton relaxation, this neglect leads to a violation
of the principle of detailed balance; equal “uphill”
and “downhill” transfer rate constants are obtained.
Here, we investigate how this result depends on nuclear reorganization
effects, neglected in Redfield but taken into account in the modified
Redfield theory. These reorganization effects, resulting from a partial
localization of excited states, are found to promote a preferential
“downhill” relaxation of excitation energy. However,
for realistic spectral densities of light-harvesting antennae in photosynthesis,
the reorganization effects are too small to compensate for the missing
coordinate–momentum uncertainty. For weaker excitonic couplings
as they occur between domains of strongly coupled pigments, we find
the principle of detailed balance to be fulfilled in a semiclassical
variant of the generalized Förster theory. A qualitatively
correct description of the transfer is obtained with this theory at
a significantly lower computational cost as with the quantum generalized
Förster theory. Larger deviations between the two theories
are expected for large energy gaps as they occur in complexes with
chemically different pigments.
Collapse
Affiliation(s)
- Thomas Renger
- Institut für Theoretische Physik, Johannes Kepler Universität Linz, Altenberger Str. 69, 4040 Linz, Austria
| |
Collapse
|
24
|
Friederich P, Häse F, Proppe J, Aspuru-Guzik A. Machine-learned potentials for next-generation matter simulations. NATURE MATERIALS 2021; 20:750-761. [PMID: 34045696 DOI: 10.1038/s41563-020-0777-6] [Citation(s) in RCA: 123] [Impact Index Per Article: 41.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/21/2019] [Accepted: 07/17/2020] [Indexed: 05/18/2023]
Abstract
The choice of simulation methods in computational materials science is driven by a fundamental trade-off: bridging large time- and length-scales with highly accurate simulations at an affordable computational cost. Venturing the investigation of complex phenomena on large scales requires fast yet accurate computational methods. We review the emerging field of machine-learned potentials, which promises to reach the accuracy of quantum mechanical computations at a substantially reduced computational cost. This Review will summarize the basic principles of the underlying machine learning methods, the data acquisition process and active learning procedures. We highlight multiple recent applications of machine-learned potentials in various fields, ranging from organic chemistry and biomolecules to inorganic crystal structure predictions and surface science. We furthermore discuss the developments required to promote a broader use of ML potentials, and the possibility of using them to help solve open questions in materials science and facilitate fully computational materials design.
Collapse
Affiliation(s)
- Pascal Friederich
- Chemical Physics Theory Group, Department of Chemistry, University of Toronto, Toronto, Ontario, Canada
- Department of Computer Science, University of Toronto, Toronto, Ontario, Canada
- Institute of Nanotechnology, Karlsruhe Institute of Technology, Eggenstein-Leopoldshafen, Germany
- Institute of Theoretical Informatics, Karlsruhe Institute of Technology, Karlsruhe, Germany
| | - Florian Häse
- Chemical Physics Theory Group, Department of Chemistry, University of Toronto, Toronto, Ontario, Canada
- Department of Computer Science, University of Toronto, Toronto, Ontario, Canada
- Vector Institute for Artificial Intelligence, Toronto, Ontario, Canada
- Department of Chemistry and Chemical Biology, Harvard University, Cambridge, MA, USA
| | - Jonny Proppe
- Chemical Physics Theory Group, Department of Chemistry, University of Toronto, Toronto, Ontario, Canada
- Department of Computer Science, University of Toronto, Toronto, Ontario, Canada
- Institute of Physical Chemistry, Georg-August University, Göttingen, Germany
| | - Alán Aspuru-Guzik
- Chemical Physics Theory Group, Department of Chemistry, University of Toronto, Toronto, Ontario, Canada.
- Department of Computer Science, University of Toronto, Toronto, Ontario, Canada.
- Vector Institute for Artificial Intelligence, Toronto, Ontario, Canada.
- Department of Chemistry and Chemical Biology, Harvard University, Cambridge, MA, USA.
- Lebovic Fellow, Canadian Institute for Advanced Research (CIFAR), Toronto, Ontario, Canada.
| |
Collapse
|
25
|
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
|
26
|
Ishizaki A, Fleming GR. Insights into Photosynthetic Energy Transfer Gained from Free-Energy Structure: Coherent Transport, Incoherent Hopping, and Vibrational Assistance Revisited. J Phys Chem B 2021; 125:3286-3295. [PMID: 33724833 DOI: 10.1021/acs.jpcb.0c09847] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Giant strides in ultrashort laser pulse technology have enabled real-time observation of dynamical processes in complex molecular systems. Specifically, the discovery of oscillatory transients in the two-dimensional electronic spectra of photosynthetic systems stimulated a number of theoretical investigations exploring the possible physical mechanisms of the remarkable quantum efficiency of light harvesting processes. In this work, we revisit the elementary aspects of environment-induced fluctuations in the involved electronic energies and present a simple way to understand energy flow with the intuitive picture of relaxation in a funnel-type free-energy landscape. The presented free-energy description of energy transfer reveals that typical photosynthetic systems operate in an almost barrierless regime. The approach also provides insights into the distinction between coherent and incoherent energy transfer and the criteria by which the necessity of the vibrational assistance is considered.
Collapse
Affiliation(s)
- Akihito Ishizaki
- Institute for Molecular Science, National Institutes of Natural Sciences, Okazaki 444-8585, Japan.,School of Physical Sciences, Graduate University for Advanced Studies, Okazaki 444-8585, Japan
| | - Graham R Fleming
- Department of Chemistry, University of California, Berkeley, California 94720, United States.,Molecular Biophysics and Integrated Bioimaging Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States.,Kavli Energy NanoSciences Institute at Berkeley, Berkeley, California 94720, United States
| |
Collapse
|
27
|
Abstract
The fundamental principle governing the robust excitation energy transfer (EET) in the light-harvesting complexes remained elusive. Recent spectroscopic observations of the EET attract attention toward the role of the quantum effects in the biological systems. In this study, along with the time evolution, the distance between the quantum states of the pigments is also considered for a detailed illustration of the effect of the quantum coherence on the speed of the excitation transfer. The comparative analysis of the coherently delocalized EET with the incoherent discrete hopping reveals that the coherent superposition speedup the excitation transfer process, for instance, with the initially localized excitation, coherence provides a ∼4-fold enhancement to the excitation transfer rate.
Collapse
Affiliation(s)
- Davinder Singh
- Korea Institute for Advanced Study, Seoul 02455, South Korea
| |
Collapse
|
28
|
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
|
29
|
Wang YC, Zhao Y. The hierarchical stochastic schrödinger equations: Theory and applications. CHINESE J CHEM PHYS 2020. [DOI: 10.1063/1674-0068/cjcp2009165] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Affiliation(s)
- Yu-Chen Wang
- State Key Laboratory of Physical Chemistry of Solid Surfaces, iChEM, Fujian Provincial Key Lab of Theoretical and Computational Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, China
| | - Yi Zhao
- State Key Laboratory of Physical Chemistry of Solid Surfaces, iChEM, Fujian Provincial Key Lab of Theoretical and Computational Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, China
| |
Collapse
|
30
|
Ng K, Webster M, Carbery WP, Visaveliya N, Gaikwad P, Jang SJ, Kretzschmar I, Eisele DM. Frenkel excitons in heat-stressed supramolecular nanocomposites enabled by tunable cage-like scaffolding. Nat Chem 2020; 12:1157-1164. [PMID: 33199886 DOI: 10.1038/s41557-020-00563-4] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/09/2019] [Accepted: 09/16/2020] [Indexed: 11/09/2022]
Abstract
Delocalized Frenkel excitons-coherently shared excitations among chromophores-are responsible for the remarkable efficiency of supramolecular light-harvesting assemblies within photosynthetic organisms. The translation of nature's design principles to applications in optoelectronic devices has been limited by the fragility of the supramolecular structures used and the delicate nature of Frenkel excitons, particularly under mildly changing solvent conditions and elevated temperatures and upon deposition onto solid substrates. Here, we overcome those functionalization barriers through composition of stable supramolecular light-harvesting nanotubes enabled by tunable (~4.3-4.9 nm), uniform (±0.3 nm) cage-like scaffolds. High-resolution cryogenic electron microscopy, combined with scanning electron microscopy, broadband femtosecond transient absorption spectroscopy and near-field scanning optical microscopy revealed that excitons within the cage-like scaffolds are robust, even under extreme heat stress, and control over nanocomposite dimensions is maintained on solid substrates. Our bio-inspired nanocomposites provide a general framework for the development of next-generation organic devices made from stable supramolecular materials.
Collapse
Affiliation(s)
- Kara Ng
- PhD Program in Chemistry, Graduate Center, The City University of New York, New York, NY, USA.,Department of Chemistry and Biochemistry, The City College of New York at The City University of New York, New York, NY, USA
| | - Megan Webster
- Department of Chemical Engineering, The City College of New York at The City University of New York, New York, NY, USA
| | - William P Carbery
- Department of Chemistry and Biochemistry, The City College of New York at The City University of New York, New York, NY, USA.,Department of Chemistry, New York University, New York, NY, USA
| | - Nikunjkumar Visaveliya
- Department of Chemistry and Biochemistry, The City College of New York at The City University of New York, New York, NY, USA
| | - Pooja Gaikwad
- PhD Program in Chemistry, Graduate Center, The City University of New York, New York, NY, USA.,Department of Chemistry and Biochemistry, The City College of New York at The City University of New York, New York, NY, USA
| | - Seogjoo J Jang
- Department of Chemistry and Biochemistry, Queens College at The City University of New York, New York, NY, USA
| | - Ilona Kretzschmar
- PhD Program in Chemistry, Graduate Center, The City University of New York, New York, NY, USA.,Department of Chemical Engineering, The City College of New York at The City University of New York, New York, NY, USA
| | - Dorthe M Eisele
- PhD Program in Chemistry, Graduate Center, The City University of New York, New York, NY, USA. .,Department of Chemistry and Biochemistry, The City College of New York at The City University of New York, New York, NY, USA.
| |
Collapse
|
31
|
Qian Y, Li X, Harutyunyan AR, Chen G, Rao Y, Chen H. Herzberg-Teller Effect on the Vibrationally Resolved Absorption Spectra of Single-Crystalline Pentacene at Finite Temperatures. J Phys Chem A 2020; 124:9156-9165. [PMID: 33103890 DOI: 10.1021/acs.jpca.0c07896] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The line shape of an electronic spectrum conveys the coupling between electronic and vibrational degrees of freedom. In the present study, the light absorption spectra of single-crystalline pentacene were measured by polarized UV-vis microscopy at 77, 185, and 293 K. The vibronic coupling encoded in each spectrum was resolved by the Herzberg-Teller theory that considers the contributions from the Franck-Condon (FC) factor, Franck-Condo/Herzberg-Teller (FC/HT) interference, and Herzberg-Teller (HT) coupling. Specifically, excitation energies, electronic transition dipole moments, and their nuclear gradients were evaluated by the GW method to ensure numerical accuracy, while the computationally efficient density function theory was employed to determine the optimized structures and vibrational normal modes. For every pair of electronic transition and normal mode that gives rise to a strong vibronic transition intensity, we examined their spatial characteristics by projecting them onto the three crystal axes. It was found that all normal modes strongly coupled to the lowest-lying a-polarized electronic transitions oscillate along axis a, whereas none of their counterparts for the lowest-lying b-polarized electronic transitions is predominantly along axis b. This notable difference on the alignment between the electronic transition and molecular vibration could help the directional control of charge dissociation and/or spin separation. Moreover, a significant variance of the destructive FC/HT interference was discovered with increasing temperatures that can well explain the a-polarized fading tableland near 650 nm. Finally, the importance of HT coupling was corroborated by comparing its intensity with those of FC factor and FC/HT interference. Taken all together, the vibrational dependence of the electronic wave function is critical to resolve the light absorption spectra of single-crystalline pentacene and its temperature effects, facilitating the systematic design of functional optical materials based on pentacene and its derivatives.
Collapse
Affiliation(s)
- Yuqin Qian
- Department of Chemistry and Biochemistry, Utah State University, Logan, Utah 84322, United States
| | - Xia Li
- Department of Chemistry and Biochemistry, Utah State University, Logan, Utah 84322, United States
| | | | - Gugang Chen
- Honda Research Institute USA, Inc., San Jose, California 95134, United States
| | - Yi Rao
- Department of Chemistry and Biochemistry, Utah State University, Logan, Utah 84322, United States
| | - Hanning Chen
- Department of Chemistry, American University, Washington, D.C. 20016, United States
| |
Collapse
|
32
|
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
|
33
|
Chen MS, Zuehlsdorff TJ, Morawietz T, Isborn CM, Markland TE. Exploiting Machine Learning to Efficiently Predict Multidimensional Optical Spectra in Complex Environments. J Phys Chem Lett 2020; 11:7559-7568. [PMID: 32808797 DOI: 10.1021/acs.jpclett.0c02168] [Citation(s) in RCA: 31] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/26/2023]
Abstract
The excited-state dynamics of chromophores in complex environments determine a range of vital biological and energy capture processes. Time-resolved, multidimensional optical spectroscopies provide a key tool to investigate these processes. Although theory has the potential to decode these spectra in terms of the electronic and atomistic dynamics, the need for large numbers of excited-state electronic structure calculations severely limits first-principles predictions of multidimensional optical spectra for chromophores in the condensed phase. Here, we leverage the locality of chromophore excitations to develop machine learning models to predict the excited-state energy gap of chromophores in complex environments for efficiently constructing linear and multidimensional optical spectra. By analyzing the performance of these models, which span a hierarchy of physical approximations, across a range of chromophore-environment interaction strengths, we provide strategies for the construction of machine learning models that greatly accelerate the calculation of multidimensional optical spectra from first principles.
Collapse
Affiliation(s)
- Michael S Chen
- Department of Chemistry, Stanford University, Stanford, California 94305, United States
| | - Tim J Zuehlsdorff
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, United States
| | - Tobias Morawietz
- Department of Chemistry, Stanford University, Stanford, California 94305, United States
| | - Christine M Isborn
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, United States
| | - Thomas E Markland
- Department of Chemistry, Stanford University, Stanford, California 94305, United States
| |
Collapse
|
34
|
Häse F, Roch LM, Friederich P, Aspuru-Guzik A. Designing and understanding light-harvesting devices with machine learning. Nat Commun 2020; 11:4587. [PMID: 32917886 PMCID: PMC7486390 DOI: 10.1038/s41467-020-17995-8] [Citation(s) in RCA: 31] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/20/2019] [Accepted: 07/16/2020] [Indexed: 01/27/2023] Open
Abstract
Understanding the fundamental processes of light-harvesting is crucial to the development of clean energy materials and devices. Biological organisms have evolved complex metabolic mechanisms to efficiently convert sunlight into chemical energy. Unraveling the secrets of this conversion has inspired the design of clean energy technologies, including solar cells and photocatalytic water splitting. Describing the emergence of macroscopic properties from microscopic processes poses the challenge to bridge length and time scales of several orders of magnitude. Machine learning experiences increased popularity as a tool to bridge the gap between multi-level theoretical models and Edisonian trial-and-error approaches. Machine learning offers opportunities to gain detailed scientific insights into the underlying principles governing light-harvesting phenomena and can accelerate the fabrication of light-harvesting devices.
Collapse
Affiliation(s)
- Florian Häse
- Department of Chemistry and Chemical Biology, Harvard University, 12 Oxford Street, Cambridge, 02138, MA, USA
- CIFAR AI Chair, Vector Institute for Artificial Intelligence, 661 University Avenue, Toronto, ON, M5S 1M1, Canada
- Department of Computer Science, University of Toronto, 214 College Street, Toronto, ON, M5S 3H6, Canada
- Chemical Physics Theory Group, Department of Chemistry, University of Toronto, 80 St. George Street, Toronto, ON, M5S 3H6, Canada
| | - Loïc M Roch
- CIFAR AI Chair, Vector Institute for Artificial Intelligence, 661 University Avenue, Toronto, ON, M5S 1M1, Canada
- Department of Computer Science, University of Toronto, 214 College Street, Toronto, ON, M5S 3H6, Canada
- Chemical Physics Theory Group, Department of Chemistry, University of Toronto, 80 St. George Street, Toronto, ON, M5S 3H6, Canada
- ChemOS Sàrl, Lausanne, VD, 1006, Switzerland
| | - Pascal Friederich
- Department of Computer Science, University of Toronto, 214 College Street, Toronto, ON, M5S 3H6, Canada
- Chemical Physics Theory Group, Department of Chemistry, University of Toronto, 80 St. George Street, Toronto, ON, M5S 3H6, Canada
- Institute of Nanotechnology, Karlsruhe Insititute of Technology, Hermann-von-Helmholtz-Platz 1, 76344, Eggenstein-Leopoldshafen, Germany
| | - Alán Aspuru-Guzik
- CIFAR AI Chair, Vector Institute for Artificial Intelligence, 661 University Avenue, Toronto, ON, M5S 1M1, Canada.
- Department of Computer Science, University of Toronto, 214 College Street, Toronto, ON, M5S 3H6, Canada.
- Chemical Physics Theory Group, Department of Chemistry, University of Toronto, 80 St. George Street, Toronto, ON, M5S 3H6, Canada.
- Lebovic Fellow, Canadian Institute for Advanced Research (CIFAR), 661 University Avenue, Toronto, ON, M5S 1M1, Canada.
| |
Collapse
|
35
|
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
|
36
|
Kim CW, Rhee YM. Toward monitoring the dissipative vibrational energy flows in open quantum systems by mixed quantum-classical simulations. J Chem Phys 2020; 152:244109. [PMID: 32610983 DOI: 10.1063/5.0009867] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023] Open
Abstract
In open quantum system dynamics, rich information about the major energy relaxation channels and corresponding relaxation rates can be elucidated by monitoring the vibrational energy flow among individual bath modes. However, such calculations often become tremendously difficult as the complexity of the subsystem-bath coupling increases. In this paper, we attempt to make this task feasible by using a mixed quantum-classical method, the Poisson-bracket mapping equation with non-Hamiltonian modification (PBME-nH) [H. W. Kim and Y. M. Rhee, J. Chem. Phys. 140, 184106 (2014)]. For a quantum subsystem bilinearly coupled to harmonic bath modes, we derive an expression for the mode energy in terms of the classical positions and momenta of the nuclei, while keeping consistency with the energy of the quantum subsystem. The accuracy of the resulting expression is then benchmarked against a numerically exact method by using relatively simple models. Although our expression predicts a qualitatively correct dissipation rate for a range of situations, cases involving a strong vibronic resonance are quite challenging. This is attributed to the inherent lack of quantum back reaction in PBME-nH, which becomes significant when the subsystem strongly interacts with a small number of bath modes. A rigorous treatment of such an effect will be crucial for developing quantitative simulation methods that can handle generic subsystem-bath coupling.
Collapse
Affiliation(s)
- Chang Woo Kim
- Department of Chemistry, Korea Advanced Institute of Science and Technology (KAIST), Daejeon 34141, South Korea
| | - Young Min Rhee
- Department of Chemistry, Korea Advanced Institute of Science and Technology (KAIST), Daejeon 34141, South Korea
| |
Collapse
|
37
|
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
|
38
|
Li Y, Zhao L, Yao Y, Guo X. Single-Molecule Nanotechnologies: An Evolution in Biological Dynamics Detection. ACS APPLIED BIO MATERIALS 2019; 3:68-85. [DOI: 10.1021/acsabm.9b00840] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Yu Li
- Beijing National Laboratory for Molecular Sciences, State Key Laboratory for Structural Chemistry of Unstable and Stable Species, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, P. R. China
| | - Lihua Zhao
- Beijing National Laboratory for Molecular Sciences, State Key Laboratory for Structural Chemistry of Unstable and Stable Species, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, P. R. China
| | - Yuan Yao
- Beijing National Laboratory for Molecular Sciences, State Key Laboratory for Structural Chemistry of Unstable and Stable Species, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, P. R. China
| | - Xuefeng Guo
- Beijing National Laboratory for Molecular Sciences, State Key Laboratory for Structural Chemistry of Unstable and Stable Species, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, P. R. China
| |
Collapse
|
39
|
Oh SA, Coker DF, Hutchinson DAW. Variety, the spice of life and essential for robustness in excitation energy transfer in light-harvesting complexes. Faraday Discuss 2019; 221:59-76. [PMID: 31552998 DOI: 10.1039/c9fd00081j] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
For over a decade there has been some significant excitement and speculation that quantum effects may be important in the excitation energy transport process in the light harvesting complexes of certain bacteria and algae, in particular via the Fenna-Matthews-Olsen (FMO) complex. Whilst the excitement may have waned somewhat with the realisation that the observed long-lived oscillations in two-dimensional electronic spectra of FMO are probably due to vibronic coherences, it remains a question whether these coherences may play any important role. We review our recent work showing how important the site-to-site variation in coupling between chloroplasts in FMO and their protein scaffold environment is for energy transport in FMO and investigate the role of vibronic modes in this transport. Whilst the effects of vibronic excitations seem modest for FMO, we show that for bilin-based pigment-protein complexes of marine algae, in particular PC645, the site-dependent vibronic excitations seem essential for robust excitation energy transport, which may again open the door for important quantum effects to be important in these photosynthetic complexes.
Collapse
Affiliation(s)
- Sue Ann Oh
- Dodd-Walls Centre for Photonic and Quantum Technologies, Department of Physics, University of Otago, Dunedin, New Zealand.
| | | | | |
Collapse
|
40
|
Claridge K, Padula D, Troisi A. On the arrangement of chromophores in light harvesting complexes: chance versus design. Faraday Discuss 2019; 221:133-149. [PMID: 31544201 DOI: 10.1039/c9fd00045c] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We used a homogeneous computational approach to derive the excitonic Hamiltonian for five light harvesting complexes containing only one type of chromophore and compare them in terms of statistical descriptors. We then studied the approximate exciton dynamics for the five complexes introducing a measure, the (averaged and time-dependent) inverse participation ratio, that enables the comparison between very diverse complexes on the same ground. We find that the global dynamics are very similar across the set of systems despite the variety of geometric structures of the complexes. In particular, the dynamics of four out of five light harvesting complexes are barely distinguishable with a small variation from the norm seen only for the Fenna-Matthews-Olson complex. We use the information from the realistic Hamiltonians to build a reduced model system that shows how the global dynamics are ultimately dominated by a single parameter, the degree of localization of the excitonic Hamiltonian eigenstates. Considering the physically plausible range of system parameters, the reduced model explains why the dynamics are so similar across most light harvesting complexes containing a single type of chromophore regardless of the detailed pattern of the inter-chromophore excitonic coupling.
Collapse
Affiliation(s)
- Kirsten Claridge
- Department of Chemistry, University of Liverpool, Liverpool L69 7ZD, UK.
| | | | | |
Collapse
|
41
|
Lim J, Bösen CM, Somoza AD, Koch CP, Plenio MB, Huelga SF. Multicolor Quantum Control for Suppressing Ground State Coherences in Two-Dimensional Electronic Spectroscopy. PHYSICAL REVIEW LETTERS 2019; 123:233201. [PMID: 31868446 DOI: 10.1103/physrevlett.123.233201] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/24/2019] [Indexed: 06/10/2023]
Abstract
The measured multidimensional spectral response of different light harvesting complexes exhibits oscillatory features which suggest an underlying coherent energy transfer. However, making this inference rigorous is challenging due to the difficulty of isolating excited state coherences in highly congested spectra. In this work, we provide a coherent control scheme that suppresses ground state coherences, thus making rephasing spectra dominated by excited state coherences. We provide a benchmark for the scheme using a model dimeric system and numerically exact methods to analyze the spectral response. We argue that combining temporal and spectral control methods can facilitate a second generation of experiments that are tailored to extract desired information and thus significantly advance our understanding of complex open many-body structure and dynamics.
Collapse
Affiliation(s)
- J Lim
- Institut für Theoretische Physik and IQST, Albert-Einstein-Allee 11, Universität Ulm, 89081 Ulm, Germany
| | - C M Bösen
- Institut für Theoretische Physik and IQST, Albert-Einstein-Allee 11, Universität Ulm, 89081 Ulm, Germany
| | - A D Somoza
- Institut für Theoretische Physik and IQST, Albert-Einstein-Allee 11, Universität Ulm, 89081 Ulm, Germany
| | - C P Koch
- Theoretische Physik, Universität Kassel, Heinrich-Plett-Strasse 40, 34132 Kassel, Germany
| | - M B Plenio
- Institut für Theoretische Physik and IQST, Albert-Einstein-Allee 11, Universität Ulm, 89081 Ulm, Germany
| | - S F Huelga
- Institut für Theoretische Physik and IQST, Albert-Einstein-Allee 11, Universität Ulm, 89081 Ulm, Germany
| |
Collapse
|
42
|
Das S, Singh AK, Biswas DS, Datta A. Dynamics of Preferential Solvation of 5-Aminoquinoline in Hexane–Alcohol Solvent Mixtures. J Phys Chem B 2019; 123:10267-10274. [DOI: 10.1021/acs.jpcb.9b09143] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Sharmistha Das
- Department of Chemistry, Indian Institute of Technology Bombay, Mumbai 400076, India
| | - Avinash Kumar Singh
- Department of Chemistry, Indian Institute of Technology Bombay, Mumbai 400076, India
| | - Deep Sekhar Biswas
- Department of Chemistry, Indian Institute of Technology Bombay, Mumbai 400076, India
| | - Anindya Datta
- Department of Chemistry, Indian Institute of Technology Bombay, Mumbai 400076, India
| |
Collapse
|
43
|
Provazza J, Coker DF. Multi-level description of the vibronic dynamics of open quantum systems. J Chem Phys 2019; 151:154114. [DOI: 10.1063/1.5120253] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022] Open
Affiliation(s)
- Justin Provazza
- Department of Chemistry, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, USA
| | - David F. Coker
- Department of Chemistry, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, USA
| |
Collapse
|
44
|
Kosenkov YK, Kosenkov D. Quantum dynamics of vibration-assisted excitation energy transfer in phycobiliprotein light-harvesting complex. J Chem Phys 2019; 151:144101. [PMID: 31615220 DOI: 10.1063/1.5116555] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Phycobiliprotein is a light-harvesting complex containing linear tetrapyrrole bilin pigments that are responsible for absorption and funneling the sun's energy in cryptophytes algae. In particular, the protein structure determines relative positions and orientations of the pigments and thus controls energy transfer pathways. The present research reveals the impact of molecular vibrations (in the 850-2700 cm-1 region) on excitation energy transfer in phycobiliprotein. The analysis of the excitation energy transfer pathways indicates a possibility of the coherent mechanism of energy transfer (delocalization) in central dihydrobiliverdin pigments and incoherent vibration-assisted energy transfer to peripheral phycocyanobilin pigments at a sub-picosecond time scale. A computational approach that enables modeling the dynamics of the excitation energy transfer with the quantum master equation formalism employing Huang-Rhys factors to describe electronic-vibrational coupling has been developed. The computational methodology has been implemented in PyFREC software.
Collapse
Affiliation(s)
- Yana Kholod Kosenkov
- Department of Chemistry and Physics, Monmouth University, 400 Cedar Avenue, West Long Branch, New Jersey 07764, USA
| | - Dmitri Kosenkov
- Department of Chemistry and Physics, Monmouth University, 400 Cedar Avenue, West Long Branch, New Jersey 07764, USA
| |
Collapse
|
45
|
Somoza AD, Marty O, Lim J, Huelga SF, Plenio MB. Dissipation-Assisted Matrix Product Factorization. PHYSICAL REVIEW LETTERS 2019; 123:100502. [PMID: 31573298 DOI: 10.1103/physrevlett.123.100502] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/03/2019] [Revised: 06/24/2019] [Indexed: 06/10/2023]
Abstract
Charge and energy transfer in biological and synthetic organic materials are strongly influenced by the coupling of electronic states to a highly structured dissipative environment. Nonperturbative simulations of these systems require a substantial computational effort, and current methods can only be applied to large systems if environmental structures are severely coarse grained. Time evolution methods based on tensor networks are fundamentally limited by the times that can be reached due to the buildup of entanglement in time, which quickly increases the size of the tensor representation, i.e., the bond dimension. In this Letter, we introduce a dissipation-assisted matrix product factorization (DAMPF) method that combines a tensor network representation of the vibronic state within a pseudomode description of the environment where a continuous bosonic environment is mapped into a few harmonic oscillators under Lindblad damping. This framework is particularly suitable for a tensor network representation, since damping suppresses the entanglement growth among oscillators and significantly reduces the bond dimension required to achieve a desired accuracy. We show that dissipation removes the "time-wall" limitation of existing methods, enabling the long-time simulation of large vibronic systems consisting of 10-50 sites coupled to 100-1000 underdamped modes in total and for a wide range of parameter regimes. For these reasons, we believe that our formalism will facilitate the investigation of spatially extended systems with applications to quantum biology, organic photovoltaics, and quantum thermodynamics.
Collapse
Affiliation(s)
- Alejandro D Somoza
- Institut für Theoretische Physik and IQST, Albert-Einstein-Allee 11, Universität Ulm, D-89069 Ulm, Germany
| | - Oliver Marty
- Institut für Theoretische Physik and IQST, Albert-Einstein-Allee 11, Universität Ulm, D-89069 Ulm, Germany
| | - James Lim
- Institut für Theoretische Physik and IQST, Albert-Einstein-Allee 11, Universität Ulm, D-89069 Ulm, Germany
| | - Susana F Huelga
- Institut für Theoretische Physik and IQST, Albert-Einstein-Allee 11, Universität Ulm, D-89069 Ulm, Germany
| | - Martin B Plenio
- Institut für Theoretische Physik and IQST, Albert-Einstein-Allee 11, Universität Ulm, D-89069 Ulm, Germany
| |
Collapse
|
46
|
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
|
47
|
Adir N, Bar-Zvi S, Harris D. The amazing phycobilisome. BIOCHIMICA ET BIOPHYSICA ACTA-BIOENERGETICS 2019; 1861:148047. [PMID: 31306623 DOI: 10.1016/j.bbabio.2019.07.002] [Citation(s) in RCA: 74] [Impact Index Per Article: 14.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/20/2019] [Revised: 06/19/2019] [Accepted: 07/09/2019] [Indexed: 10/26/2022]
Abstract
Cyanobacteria and red-algae share a common light-harvesting complex which is different than all other complexes that serve as photosynthetic antennas - the Phycobilisome (PBS). The PBS is found attached to the stromal side of thylakoid membranes, filling up most of the gap between individual thylakoids. The PBS self assembles from similar homologous protein units that are soluble and contain conserved cysteine residues that covalently bind the light absorbing chromophores, linear tetra-pyrroles. Using similar construction principles, the PBS can be as large as 16.8 MDa (68×45×39nm), as small as 1.2 MDa (24 × 11.5 × 11.5 nm), and in some unique cases smaller still. The PBS can absorb light between 450 nm to 650 nm and in some cases beyond 700 nm, depending on the species, its composition and assembly. In this review, we will present new observations and structures that expand our understanding of the distinctive properties that make the PBS an amazing light harvesting system. At the end we will suggest why the PBS, for all of its excellent properties, was discarded by photosynthetic organisms that arose later in evolution such as green algae and higher plants.
Collapse
Affiliation(s)
- Noam Adir
- Schulich Faculty of Chemistry, Technion-Israel Institute of Technology, Haifa 32000, Israel.
| | - Shira Bar-Zvi
- Schulich Faculty of Chemistry, Technion-Israel Institute of Technology, Haifa 32000, Israel
| | - Dvir Harris
- Schulich Faculty of Chemistry, Technion-Israel Institute of Technology, Haifa 32000, Israel
| |
Collapse
|
48
|
Corbella M, Cupellini L, Lipparini F, Scholes GD, Curutchet C. Spectral Variability in Phycocyanin Cryptophyte Antenna Complexes is Controlled by Changes in the α‐Polypeptide Chains. CHEMPHOTOCHEM 2019. [DOI: 10.1002/cptc.201900045] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Affiliation(s)
- Marina Corbella
- Department of Pharmacy and Pharmaceutical Technology and Physical Chemistry and Institute of Theoretical and Computational Chemistry (IQTC-UB), Faculty of Pharmacy and Food SciencesUniversity of Barcelona Av. Joan XXIII s/n 08028 Barcelona Spain
- Department of ChemistryUppsala University BMC Box 576 Uppsala S-751 23 Sweden
| | - Lorenzo Cupellini
- Dipartimento di Chimica e Chimica IndustrialeUniversity of Pisa Via Risorgimento 35 56126 Pisa Italy
- Institute for Research in Biomedicine (IRB Barcelona)The Barcelona Institute of Science and Technology Baldiri Reixac 10 08028 Barcelona Spain
| | - Filippo Lipparini
- Dipartimento di Chimica e Chimica IndustrialeUniversity of Pisa Via Risorgimento 35 56126 Pisa Italy
| | - Gregory D. Scholes
- Department of ChemistryPrinceton University Washington Road, Princeton New Jersey 08544 United States
| | - Carles Curutchet
- Department of Pharmacy and Pharmaceutical Technology and Physical Chemistry and Institute of Theoretical and Computational Chemistry (IQTC-UB), Faculty of Pharmacy and Food SciencesUniversity of Barcelona Av. Joan XXIII s/n 08028 Barcelona Spain
| |
Collapse
|
49
|
Wang YC, Zhao Y. Effect of an underdamped vibration with both diagonal and off-diagonal exciton-phonon interactions on excitation energy transfer. J Comput Chem 2019; 40:1097-1104. [PMID: 30549065 DOI: 10.1002/jcc.25611] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/29/2018] [Revised: 08/29/2018] [Accepted: 09/06/2018] [Indexed: 11/08/2022]
Abstract
A numerically exact approach, named as the hierarchical stochastic Schrödinger equation, is employed to investigate the resonant vibration-assisted excitation energy transfer in a dimer system, where an underdamped vibration with both diagonal and off-diagonal exciton-phonon interactions is incorporated. From a large parameter space over the site-energy difference, excitonic coupling, and reorganization energy, it is found that the promotion effect of the underdamped vibration is significant only when the excitonic coupling is smaller than the site-energy difference. Under the circumstance, there is an optimal strength ratio between diagonal and off-diagonal exciton-phonon interactions for the resonant vibration-assisted excitation energy transfer as the site-energy difference is greater than the reorganization energy, whereas in the opposite situation the most efficient energy transfer occurs as the exciton-phonon interaction is totally off-diagonal. © 2018 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Yu-Chen Wang
- State Key Laboratory for Physical Chemistry of Solid Surfaces, Collaborative Innovation Center of Chemistry for Energy Materials, Fujian Provincial Key Lab of Theoretical and Computational Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, China
| | - Yi Zhao
- State Key Laboratory for Physical Chemistry of Solid Surfaces, Collaborative Innovation Center of Chemistry for Energy Materials, Fujian Provincial Key Lab of Theoretical and Computational Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, China
| |
Collapse
|
50
|
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
|