1
|
Jaiswal VK, Aranda Ruiz D, Petropoulos V, Kabaciński P, Montorsi F, Uboldi L, Ugolini S, Mukamel S, Cerullo G, Garavelli M, Santoro F, Nenov A. Sub-100-fs energy transfer in coenzyme NADH is a coherent process assisted by a charge-transfer state. Nat Commun 2024; 15:4900. [PMID: 38851775 PMCID: PMC11162464 DOI: 10.1038/s41467-024-48871-4] [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: 09/18/2023] [Accepted: 05/15/2024] [Indexed: 06/10/2024] Open
Abstract
Excitation energy transfer (EET) is a key photoinduced process in biological chromophoric assemblies. Here we investigate the factors which can drive EET into efficient ultrafast sub-ps regimes. We demonstrate how a coherent transport of electronic population could facilitate this in water solvated NADH coenzyme and uncover the role of an intermediate dark charge-transfer state. High temporal resolution ultrafast optical spectroscopy gives a 54±11 fs time constant for the EET process. Nonadiabatic quantum dynamical simulations computed through the time-evolution of multidimensional wavepackets suggest that the population transfer is mediated by photoexcited molecular vibrations due to strong coupling between the electronic states. The polar aqueous solvent environment leads to the active participation of a dark charge transfer state, accelerating the vibronically coherent EET process in favorably stacked conformers and solvent cavities. Our work demonstrates how the interplay of structural and environmental factors leads to diverse pathways for the EET process in flexible heterodimers and provides general insights relevant for coherent EET processes in stacked multichromophoric aggregates like DNA strands.
Collapse
Affiliation(s)
- Vishal Kumar Jaiswal
- Dipartimento di Chimica industriale "Toso Montanari", Università di Bologna, Viale del Risorgimento 4, 40136, Bologna, Italy
| | - Daniel Aranda Ruiz
- ICMol, Universidad de Valencia, Catedrático José Beltrán Martínez, 2, 46980, Paterna, Spain
| | - Vasilis Petropoulos
- Dipartimento di Fisica, Politecnico di Milano, Piazza Leonardo da Vinci 32, 20133, Milano, Italy
| | - Piotr Kabaciński
- Dipartimento di Fisica, Politecnico di Milano, Piazza Leonardo da Vinci 32, 20133, Milano, Italy
| | - Francesco Montorsi
- Dipartimento di Chimica industriale "Toso Montanari", Università di Bologna, Viale del Risorgimento 4, 40136, Bologna, Italy
| | - Lorenzo Uboldi
- Dipartimento di Fisica, Politecnico di Milano, Piazza Leonardo da Vinci 32, 20133, Milano, Italy
| | - Simone Ugolini
- Dipartimento di Chimica industriale "Toso Montanari", Università di Bologna, Viale del Risorgimento 4, 40136, Bologna, Italy
| | - Shaul Mukamel
- Department of Chemistry and Department of Physics and Astronomy, University of California, Irvine, CA, 92697, USA
| | - Giulio Cerullo
- Dipartimento di Fisica, Politecnico di Milano, Piazza Leonardo da Vinci 32, 20133, Milano, Italy.
| | - Marco Garavelli
- Dipartimento di Chimica industriale "Toso Montanari", Università di Bologna, Viale del Risorgimento 4, 40136, Bologna, Italy.
| | - Fabrizio Santoro
- Istituto di Chimica dei Composti Organometallici (ICCOM-CNR), Area della Ricerca del CNR, Via Moruzzi 1, I-56124, Pisa, Italy
| | - Artur Nenov
- Dipartimento di Chimica industriale "Toso Montanari", Università di Bologna, Viale del Risorgimento 4, 40136, Bologna, Italy.
| |
Collapse
|
2
|
Liu Z, Hu H, Sun X. Multistate Reaction Coordinate Model for Charge and Energy Transfer Dynamics in the Condensed Phase. J Chem Theory Comput 2023; 19:7151-7170. [PMID: 37815937 PMCID: PMC10601487 DOI: 10.1021/acs.jctc.3c00770] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/13/2023] [Indexed: 10/12/2023]
Abstract
Constructing multistate model Hamiltonians from all-atom electronic structure calculations and molecular dynamics simulations is crucial for understanding charge and energy transfer dynamics in complex condensed phases. The most popular two-level system model is the spin-boson Hamiltonian, where the nuclear degrees of freedom are represented as shifted normal modes. Recently, we proposed the general multistate nontrivial extension of the spin-boson model, i.e., the multistate harmonic (MSH) model, which is constructed by extending the spatial dimensions of each nuclear mode so as to satisfy the all-atom reorganization energy restrictions for all pairs of electronic states. In this work, we propose the multistate reaction coordinate (MRC) model with a primary reaction coordinate and secondary bath modes as in the Caldeira-Leggett form but in extended spatial dimensions. The MRC model is proven to be equivalent to the MSH model and offers an intuitive physical picture of the nuclear-electronic feedback in nonadiabatic processes such as the inherent trajectory of the reaction coordinate. The reaction coordinate is represented in extended dimensions, carrying the entire reorganization energies and bilinearly coupled to the secondary bath modes. We demonstrate the MRC model construction for photoinduced charge transfer in an organic photovoltaic caroteniod-porphyrin-C60 molecular triad dissolved in tetrahydrofuran as well as excitation energy transfer in a photosynthetic light-harvesting Fenna-Matthews-Olson complex. The MRC model provides an effective and robust platform for investigating quantum dissipative dynamics in complex condensed-phase systems since it allows a consistent description of realistic spectral density, state-dependent system-bath couplings, and heterogeneous environments due to static disorder in reorganization energies.
Collapse
Affiliation(s)
- Zengkui Liu
- Division
of Arts and Sciences, NYU Shanghai, 567 West Yangsi Road, Shanghai, 200124, China
- NYU-ECNU
Center for Computational Chemistry at NYU Shanghai, 3663 Zhongshan Road North, Shanghai, 200062, China
- Department
of Chemistry, New York University, New York, New York, 10003, United States
| | - Haorui Hu
- Division
of Arts and Sciences, NYU Shanghai, 567 West Yangsi Road, Shanghai, 200124, China
| | - Xiang Sun
- Division
of Arts and Sciences, NYU Shanghai, 567 West Yangsi Road, Shanghai, 200124, China
- NYU-ECNU
Center for Computational Chemistry at NYU Shanghai, 3663 Zhongshan Road North, Shanghai, 200062, China
- Department
of Chemistry, New York University, New York, New York, 10003, United States
- Shanghai
Frontiers Science Center of Artificial Intelligence and Deep Learning, NYU Shanghai, 567 West Yangsi Road, Shanghai, 200124, China
| |
Collapse
|
3
|
Lyu N, Mulvihill E, Soley MB, Geva E, Batista VS. Tensor-Train Thermo-Field Memory Kernels for Generalized Quantum Master Equations. J Chem Theory Comput 2023; 19:1111-1129. [PMID: 36719350 DOI: 10.1021/acs.jctc.2c00892] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2023]
Abstract
The generalized quantum master equation (GQME) approach provides a rigorous framework for deriving the exact equation of motion for any subset of electronic reduced density matrix elements (e.g., the diagonal elements). In the context of electronic dynamics, the memory kernel and inhomogeneous term of the GQME introduce the implicit coupling to nuclear motion and dynamics of electronic density matrix elements that are projected out (e.g., the off-diagonal elements), allowing for efficient quantum dynamics simulations. Here, we focus on benchmark quantum simulations of electronic dynamics in a spin-boson model system described by various types of GQMEs. Exact memory kernels and inhomogeneous terms are obtained from short-time quantum-mechanically exact tensor-train thermo-field dynamics (TT-TFD) simulations and are compared with those obtained from an approximate linearized semiclassical method, allowing for assessment of the accuracy of these approximate memory kernels and inhomogeneous terms. Moreover, we have analyzed the computational cost of the full and reduced-dimensionality GQMEs. The scaling of the computational cost is dependent on several factors, sometimes with opposite scaling trends. The TT-TFD memory kernels can provide insights on the main sources of inaccuracies of GQME approaches when combined with approximate input methods and pave the road for the development of quantum circuits that implement GQMEs on digital quantum computers.
Collapse
Affiliation(s)
- Ningyi Lyu
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States
| | - Ellen Mulvihill
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States
| | - Micheline B Soley
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States.,Yale Quantum Institute, Yale University, New Haven, Connecticut 06511, United States.,Department of Chemistry, University of Wisconsin-Madison, Madison, Wisconsin 53706, United States
| | - Eitan Geva
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| | - Victor S Batista
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States.,Yale Quantum Institute, Yale University, New Haven, Connecticut 06511, United States
| |
Collapse
|
4
|
Amati G, Runeson JE, Richardson JO. On detailed balance in nonadiabatic dynamics: From spin spheres to equilibrium ellipsoids. J Chem Phys 2023; 158:064113. [PMID: 36792511 DOI: 10.1063/5.0137828] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023] Open
Abstract
Trajectory-based methods that propagate classical nuclei on multiple quantum electronic states are often used to simulate nonadiabatic processes in the condensed phase. A long-standing problem of these methods is their lack of detailed balance, meaning that they do not conserve the equilibrium distribution. In this article, we investigate ideas for restoring detailed balance in mixed quantum-classical systems by tailoring the previously proposed spin-mapping approach to thermal equilibrium. We find that adapting the spin magnitude can recover the correct long-time populations but is insufficient to conserve the full equilibrium distribution. The latter can however be achieved by a more flexible mapping of the spin onto an ellipsoid, which is constructed to fulfill detailed balance for arbitrary potentials. This ellipsoid approach solves the problem of negative populations that has plagued previous mapping approaches and can therefore be applied also to strongly asymmetric and anharmonic systems. Because it conserves the thermal distribution, the method can also exploit efficient sampling schemes used in standard molecular dynamics, which drastically reduces the number of trajectories needed for convergence. The dynamics does however still have mean-field character, as is observed most clearly by evaluating reaction rates in the golden-rule limit. This implies that although the ellipsoid mapping provides a rigorous framework, further work is required to find an accurate classical-trajectory approximation that captures more properties of the true quantum dynamics.
Collapse
Affiliation(s)
- Graziano Amati
- Laboratory of Physical Chemistry, ETH Zürich, 8093 Zürich, Switzerland
| | - Johan E Runeson
- Laboratory of Physical Chemistry, ETH Zürich, 8093 Zürich, Switzerland
| | | |
Collapse
|
5
|
Burger A, Kwek LC, Poletti D. Digital Quantum Simulation of the Spin-Boson Model under Markovian Open-System Dynamics. ENTROPY (BASEL, SWITZERLAND) 2022; 24:1766. [PMID: 36554171 PMCID: PMC9777563 DOI: 10.3390/e24121766] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 11/07/2022] [Revised: 11/28/2022] [Accepted: 11/29/2022] [Indexed: 06/17/2023]
Abstract
Digital quantum computers have the potential to simulate complex quantum systems. The spin-boson model is one of such systems, used in disparate physical domains. Importantly, in a number of setups, the spin-boson model is open, i.e., the system is in contact with an external environment which can, for instance, cause the decay of the spin state. Here, we study how to simulate such open quantum dynamics in a digital quantum computer, for which we use an IBM hardware. We consider in particular how accurate different implementations of the evolution result as a function of the level of noise in the hardware and of the parameters of the open dynamics. For the regimes studied, we show that the key aspect is to simulate the unitary portion of the dynamics, while the dissipative part can lead to a more noise-resistant simulation. We consider both a single spin coupled to a harmonic oscillator, and also two spins coupled to the oscillator. In the latter case, we show that it is possible to simulate the emergence of correlations between the spins via the oscillator.
Collapse
Affiliation(s)
- Andreas Burger
- Faculty of Physics, Ludwig-Maximilians-Universität Munich, Geschwister-Scholl-Platz 1, 80539 Munich, Germany
- Science, Mathematics and Technology Cluster, Singapore University of Technology and Design, Singapore 487372, Singapore
- Centre for Quantum Technologies, National University of Singapore, Singapore 117543, Singapore
| | - Leong Chuan Kwek
- Centre for Quantum Technologies, National University of Singapore, Singapore 117543, Singapore
- National Institute of Education and Institute of Advanced Studies, Nanyang Technological University, Singapore 637616, Singapore
- MajuLab, CNRS-UNS-NUS-NTU International Joint Research Unit, Singapore 117543, Singapore
| | - Dario Poletti
- Science, Mathematics and Technology Cluster, Singapore University of Technology and Design, Singapore 487372, Singapore
- Centre for Quantum Technologies, National University of Singapore, Singapore 117543, Singapore
- National Institute of Education and Institute of Advanced Studies, Nanyang Technological University, Singapore 637616, Singapore
- EPD Pillar, Singapore University of Technology and Design, 8 Somapah Road, Singapore 487372, Singapore
- The Abdus Salam International Centre for Theoretical Physics, Strada Costiera 11, 34151 Trieste, Italy
| |
Collapse
|
6
|
Liu Z, Xu W, Tuckerman ME, Sun X. Imaginary-Time Open-Chain Path-Integral Approach for Two-State Time Correlation Functions and Applications in Charge Transfer. J Chem Phys 2022; 157:114111. [DOI: 10.1063/5.0098162] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Quantum time correlation functions (TCFs) involving two states are important for describing nonadiabatic dynamical processes such as charge transfer. Based on a previous single-state method, we propose an imaginary-time open-chain path-integral (OCPI) approach for evaluating the two-state symmetrized TCFs. Expressing the forward and backward propagation on different electronic potential energy surfaces as a complex-time path integral, we then transform the path variables to average and difference variables such that the integration over the difference variables up to the second order can be performed analytically. The resulting expression for the symmetrized TCF is equivalent to sampling the open-chain configurations in an effective potential that corresponds to the average surface. Using importance sampling over the extended OCPI space via open path integral molecular dynamics, we tested the resulting path-integral approximation by calculating the Fermi's golden rule charge transfer rate constant within a widely-used spin-boson model. Comparing with the real-time linearized semiclassical method and analytical result, we show that the imaginary-time OCPI provides an accurate two-state symmetrized TCF and rate constant in the typical turnover region. It is shown that the first bead of the open chain corresponds to physical zero-time, and the endpoint bead corresponds to final time t; oscillations of the end-to-end distance perfectly match the nuclear mode frequency. The two-state OCPI scheme is seen to capture the tested model's electronic quantum coherence and nuclear quantum effects accurately.
Collapse
Affiliation(s)
- Zengkui Liu
- Division of Arts and Sciences, New York University Shanghai, China
| | - Wen Xu
- New York University Shanghai, China
| | - Mark E. Tuckerman
- Department of Chemistry and Courant Institute of Mathematical Sciences, New York University, United States of America
| | - Xiang Sun
- Division of Arts and Sciences, New York University Shanghai, China
| |
Collapse
|
7
|
Lai Y, Geva E. Electronic Absorption Spectra from Off-Diagonal Quantum Master Equations. J Chem Phys 2022; 157:104115. [DOI: 10.1063/5.0106888] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Quantum master equations (QMEs) provide a general framework for describing electronic dynamics within a complex molecular system. Off-diagonal QMEs (OD-QMEs) correspond to a family of QMEs that describe the electronic dynamics in the interaction picture based on treating the off-diagonal coupling terms between electronic states as a small perturbation within the framework of second-order perturbation theory. The fact that OD-QMEs are given in terms of the interaction picture makes it non-trivial to obtain Schrodinger picture electronic coherences from them. A key experimental quantity that relies on the ability to obtain accurate Schrodinger picture electronic coherences is the absorption spectrum. In this paper, we propose using a recently introduced procedure for extracting Schrodinger picture electronic coherences from interaction picture inputs to calculate electronic absorption spectra from electronic dynamics generated by OD-QMEs. The accuracy of the absorption spectra obtained in this way is studied in the context of a biexciton benchmark model, by comparing spectra calculated based on time-local and time-nonlocal OD-QMEs to spectra calculated based on a Redfield-type QME and the non-perturbative and quantum-mechanically exact hierarchical equations of motion (HEOM) method.
Collapse
Affiliation(s)
- Yifan Lai
- Chemistry, University of Michigan, United States of America
| | - Eitan Geva
- Department of Chemistry, University of Michigan Department of Chemistry, United States of America
| |
Collapse
|
8
|
Mulvihill E, Geva E. Simulating the dynamics of electronic observables via reduced-dimensionality generalized quantum master equations. J Chem Phys 2022; 156:044119. [DOI: 10.1063/5.0078040] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023] Open
Affiliation(s)
- Ellen Mulvihill
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
| | - Eitan Geva
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
| |
Collapse
|
9
|
Wu D, Hu Z, Li J, Sun X. Forecasting nonadiabatic dynamics using hybrid convolutional neural network/long short-term memory network. J Chem Phys 2021; 155:224104. [PMID: 34911307 DOI: 10.1063/5.0073689] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022] Open
Abstract
Modeling nonadiabatic dynamics in complex molecular or condensed-phase systems has been challenging, especially for the long-time dynamics. In this work, we propose a time series machine learning scheme based on the hybrid convolutional neural network/long short-term memory (CNN-LSTM) framework for predicting the long-time quantum behavior, given only the short-time dynamics. This scheme takes advantage of both the powerful local feature extraction ability of CNN and the long-term global sequential pattern recognition ability of LSTM. With feature fusion of individually trained CNN-LSTM models for the quantum population and coherence dynamics, the proposed scheme is shown to have high accuracy and robustness in predicting the linearized semiclassical and symmetrical quasiclassical mapping dynamics as well as the mixed quantum-classical Liouville dynamics of various spin-boson models with learning time up to 0.3 ps. Furthermore, if the hybrid network has learned the dynamics of a system, this knowledge is transferable that could significantly enhance the accuracy in predicting the dynamics of a similar system. The hybrid CNN-LSTM network is thus believed to have high predictive power in forecasting the nonadiabatic dynamics in realistic charge and energy transfer processes in photoinduced energy conversion.
Collapse
Affiliation(s)
- Daxin Wu
- Division of Arts and Sciences, NYU Shanghai, 1555 Century Avenue, Shanghai 200122, China
| | - Zhubin Hu
- Division of Arts and Sciences, NYU Shanghai, 1555 Century Avenue, Shanghai 200122, China
| | - Jiebo Li
- Beijing Advanced Innovation Center for Biomedical Engineering, School of Biological Science and Medical Engineering, Institute of Medical Photonics, Beihang University, Beijing 100191, China
| | - Xiang Sun
- Division of Arts and Sciences, NYU Shanghai, 1555 Century Avenue, Shanghai 200122, China
| |
Collapse
|
10
|
Lai Y, Geva E. On simulating the dynamics of electronic populations and coherences via quantum master equations based on treating off-diagonal electronic coupling terms as a small perturbation. J Chem Phys 2021; 155:204101. [PMID: 34852488 DOI: 10.1063/5.0069313] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Quantum master equations provide a general framework for describing the dynamics of electronic observables within a complex molecular system. One particular family of such equations is based on treating the off-diagonal coupling terms between electronic states as a small perturbation within the framework of second-order perturbation theory. In this paper, we show how different choices of projection operators, as well as whether one starts out with the time-convolution or the time-convolutionless forms of the generalized quantum master equation, give rise to four different types of such off-diagonal quantum master equations (OD-QMEs), namely, time-convolution and time-convolutionless versions of a Pauli-type OD-QME for only the electronic populations and an OD-QME for the full electronic density matrix (including both electronic populations and coherences). The fact that those OD-QMEs are given in terms of the interaction picture makes it non-trivial to obtain Schrödinger picture electronic coherences from them. To address this, we also extend a procedure for extracting Schrödinger picture electronic coherences from interaction picture populations recently introduced by Trushechkin in the context of time-convolutionless Pauli-type OD-QME to the other three types of OD-QMEs. The performance of the aforementioned four types of OD-QMEs is explored in the context of the Garg-Onuchic-Ambegaokar benchmark model for charge transfer in the condensed phase across a relatively wide parameter range. The results show that time-convolution OD-QMEs can be significantly more accurate than their time-convolutionless counterparts, particularly in the case of Pauli-type OD-QMEs, and that rather accurate Schrödinger picture coherences can be obtained from interaction picture electronic inputs.
Collapse
Affiliation(s)
- Yifan Lai
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
| | - Eitan Geva
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
| |
Collapse
|
11
|
Brian D, Sun X. Generalized quantum master equation: A tutorial review and recent advances. CHINESE J CHEM PHYS 2021. [DOI: 10.1063/1674-0068/cjcp2109157] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Affiliation(s)
- Dominikus Brian
- Division of Arts and Sciences, NYU Shanghai, Shanghai 200122, China
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China
- Department of Chemistry, New York University, New York 10003, USA
| | - Xiang Sun
- Division of Arts and Sciences, NYU Shanghai, Shanghai 200122, China
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China
- Department of Chemistry, New York University, New York 10003, USA
- State Key Laboratory of Precision Spectroscopy, East China Normal University, Shanghai 200241, China
| |
Collapse
|
12
|
Hu Z, Brian D, Sun X. Multi-state harmonic models with globally shared bath for nonadiabatic dynamics in the condensed phase. J Chem Phys 2021; 155:124105. [PMID: 34598571 DOI: 10.1063/5.0064763] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023] Open
Abstract
Model Hamiltonians constructed from quantum chemistry calculations and molecular dynamics simulations are widely used for simulating nonadiabatic dynamics in the condensed phase. The most popular two-state spin-boson model could be built by mapping the all-atom anharmonic Hamiltonian onto a two-level system bilinearly coupled to a harmonic bath using the energy gap time correlation function. However, for more than two states, there lacks a general strategy to construct multi-state harmonic (MSH) models since the energy gaps between different pairs of electronic states are not entirely independent and need to be considered consistently. In this paper, we extend the previously proposed approach for building three-state harmonic models for photoinduced charge transfer to the arbitrary number of electronic states with a globally shared bath and the system-bath couplings are scaled differently according to the reorganization energies between each pair of states. We demonstrate the MSH model construction for an organic photovoltaic carotenoid-porphyrin-C60 molecular triad dissolved in explicit tetrahydrofuran solvent. Nonadiabatic dynamics was simulated using mixed quantum-classical techniques, including the linearized semiclassical and symmetrical quasiclassical dynamics with the mapping Hamiltonians, mean-field Ehrenfest, and mixed quantum-classical Liouville dynamics in two-state, three-state, and four-state harmonic models of the triad system. The MSH models are shown to provide a general and flexible framework for simulating nonadiabatic dynamics in complex systems.
Collapse
Affiliation(s)
- Zhubin Hu
- Division of Arts and Sciences, NYU Shanghai, 1555 Century Avenue, Shanghai 200122, China
| | - Dominikus Brian
- Division of Arts and Sciences, NYU Shanghai, 1555 Century Avenue, Shanghai 200122, China
| | - Xiang Sun
- Division of Arts and Sciences, NYU Shanghai, 1555 Century Avenue, Shanghai 200122, China
| |
Collapse
|
13
|
Mulvihill E, Geva E. A Road Map to Various Pathways for Calculating the Memory Kernel of the Generalized Quantum Master Equation. J Phys Chem B 2021; 125:9834-9852. [PMID: 34424700 DOI: 10.1021/acs.jpcb.1c05719] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The generalized quantum master equation (GQME) provides a powerful framework for simulating electronic energy, charge, and coherence transfer dynamics in molecular systems. Within this framework, the effect of the nuclear degrees of freedom on the time evolution of the electronic reduced density matrix is fully captured by a memory kernel superoperator. However, the actual memory kernel depends on the choice of projection operator and is therefore not unique. Furthermore, calculating the memory kernel can be done in multiple ways that use different forms of projection-free inputs. Although the electronic dynamics is invariant to those choices when quantum-mechanically exact projection-free inputs are used, this is not the case when they are obtained via more feasible semiclassical or mixed quantum-classical approximate methods. Furthermore, the accuracy and numerical stability of the resulting electronic dynamics has been observed to be sensitive to the above-mentioned choices when approximate methods are used to calculate the projection-free inputs. In this article, we provide a systematic road map to 30 possible pathways for calculating the memory kernel and highlight how they are related as well as the ways in which they differ. We also compare the performance of different pathways in the context of the spin-boson benchmark model, with the projection-free inputs obtained via a mapping Hamiltonian linearized semiclassical method. In this case, we find that expressing the memory kernel with an exponential operator where the projection operator precedes the Liouvillian yields the most accurate and most numerically stable results.
Collapse
Affiliation(s)
- Ellen Mulvihill
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan48109, United States
| | - Eitan Geva
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan48109, United States
| |
Collapse
|
14
|
Osswald M, Fingerhut BP. Electron Transfer-Induced Active Site Structural Relaxation in 64-Photolyase of Drosophila melanogaster. J Phys Chem B 2021; 125:8690-8702. [PMID: 34323497 DOI: 10.1021/acs.jpcb.1c02951] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
While catalytic electron flow and photoreactivation of CPD-photolyases are increasingly understood, the microscopic details of the 64-photolyase repair mechanism are perpetually debated. Here, we investigate in long-time (μs) molecular dynamics simulations combined with extensive quantum mechanical/molecular mechanical (QM/MM) simulations the primary electron transfer (ET) reactions in 64-photolyase of Drosophila melanogaster (D. melanogaster). The characterization of the relative energetics of locally excited and charge separated states in the (6-4) photoproduct enzyme repair complex reveals a charge-separated state involving the adenine moiety of the FADH- cofactor that facilitates reduction of the photoproduct. Microscopic details of the collective reaction coordinate of ET reactions are identified that involve the reorganization of the hydrogen bond network and structural relaxation of the active site. The simulations reveal complex active site relaxation dynamics involving distinguished amino acids (Lys246, His365, and His369), conformational reorganization of the hydroxyl group of the (6-4) photoproduct, and a strengthening of hydrogen bonds with immobilized water molecules. In particular, rotation of the Lys246 side chain is found to impose a double-well character along the reaction coordinate of the ET reaction. Our findings suggest that the primary ET reactions in the (6-4) photoproduct enzyme repair complex of D. melanogaster are governed by a complex multi-minima active site relaxation dynamics and potentially precede the equilibration of the protein. ET pathways mediated by the adenine moiety and the 5' side of the photoproduct are proposed to be relevant for triggering the catalytic (6-4) photoproduct reactivation.
Collapse
Affiliation(s)
- Mara Osswald
- Max-Born-Institut für Nichtlineare Optik und Kurzzeitspektroskopie, D-12489 Berlin, Germany
| | - Benjamin P Fingerhut
- Max-Born-Institut für Nichtlineare Optik und Kurzzeitspektroskopie, D-12489 Berlin, Germany
| |
Collapse
|
15
|
Mulvihill E, Lenn KM, Gao X, Schubert A, Dunietz BD, Geva E. Simulating energy transfer dynamics in the Fenna-Matthews-Olson complex via the modified generalized quantum master equation. J Chem Phys 2021; 154:204109. [PMID: 34241158 DOI: 10.1063/5.0051101] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
The generalized quantum master equation (GQME) provides a general and formally exact framework for simulating the reduced dynamics of open quantum systems. The recently introduced modified approach to the GQME (M-GQME) corresponds to a specific implementation of the GQME that is geared toward simulating the dynamics of the electronic reduced density matrix in systems governed by an excitonic Hamiltonian. Such a Hamiltonian, which is often used for describing energy and charge transfer dynamics in complex molecular systems, is given in terms of diabatic electronic states that are coupled to each other and correspond to different nuclear Hamiltonians. Within the M-GQME approach, the effect of the nuclear degrees of freedom on the time evolution of the electronic density matrix is fully captured by a memory kernel superoperator, which can be obtained from short-lived (compared to the time scale of energy/charge transfer) projection-free inputs. In this paper, we test the ability of the M-GQME to predict the energy transfer dynamics within a seven-state benchmark model of the Fenna-Matthews-Olson (FMO) complex, with the short-lived projection-free inputs obtained via the Ehrenfest method. The M-GQME with Ehrenfest-based inputs is shown to yield accurate results across a wide parameter range. It is also found to dramatically outperform the direct application of the Ehrenfest method and to provide better-behaved convergence with respect to memory time in comparison to an alternative implementation of the GQME approach previously applied to the same FMO model.
Collapse
Affiliation(s)
- Ellen Mulvihill
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
| | - Kristina M Lenn
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
| | - Xing Gao
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
| | - Alexander Schubert
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
| | - Barry D Dunietz
- Department of Chemistry and Biochemistry, Kent State University, Kent, Ohio 44242, USA
| | - Eitan Geva
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
| |
Collapse
|
16
|
Gao X, Geva E. Improving the Accuracy of Quasiclassical Mapping Hamiltonian Methods by Treating the Window Function Width as an Adjustable Parameter. J Phys Chem A 2020; 124:11006-11016. [DOI: 10.1021/acs.jpca.0c09750] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Xing Gao
- School of Materials, Sun Yat-sen University, Shenzhen, Guangdong 518100, China
| | - Eitan Geva
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| |
Collapse
|
17
|
Gao X, Saller MAC, Liu Y, Kelly A, Richardson JO, Geva E. Benchmarking Quasiclassical Mapping Hamiltonian Methods for Simulating Electronically Nonadiabatic Molecular Dynamics. J Chem Theory Comput 2020; 16:2883-2895. [DOI: 10.1021/acs.jctc.9b01267] [Citation(s) in RCA: 37] [Impact Index Per Article: 9.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Xing Gao
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| | | | - Yudan Liu
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| | - Aaron Kelly
- Department of Chemistry, Dalhousie University, 15000 Halifax, Nova Scotia, Canada
| | | | - Eitan Geva
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| |
Collapse
|
18
|
Tacchino F, Succurro A, Ebenhöh O, Gerace D. Optimal efficiency of the Q-cycle mechanism around physiological temperatures from an open quantum systems approach. Sci Rep 2019; 9:16657. [PMID: 31723177 PMCID: PMC6853958 DOI: 10.1038/s41598-019-52842-x] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/27/2019] [Accepted: 10/14/2019] [Indexed: 01/17/2023] Open
Abstract
The Q-cycle mechanism entering the electron and proton transport chain in oxygenic photosynthesis is an example of how biological processes can be efficiently investigated with elementary microscopic models. Here we address the problem of energy transport across the cellular membrane from an open quantum system theoretical perspective. We model the cytochrome \documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${b}_{6}\,f$$\end{document}b6f protein complex under cyclic electron flow conditions starting from a simplified kinetic model, which is hereby revisited in terms of a Markovian quantum master equation formulation and spin-boson Hamiltonian treatment. We apply this model to theoretically demonstrate an optimal thermodynamic efficiency of the Q-cycle around ambient and physiologically relevant temperature conditions. Furthermore, we determine the quantum yield of this complex biochemical process after setting the electrochemical potentials to values well established in the literature. The present work suggests that the theory of quantum open systems can successfully push forward our theoretical understanding of complex biological systems working close to the quantum/classical boundary.
Collapse
Affiliation(s)
| | - Antonella Succurro
- Institute of Quantitative and Theoretical Biology, Heinrich Heine University, 40225, Düsseldorf, Germany.,Cluster of Excellence on Plant Sciences (CEPLAS), Heinrich Heine University, 40225, Düsseldorf, Germany.,Life and Medical Sciences (LIMES) Institute and West German Genome Center, University of Bonn, Bonn, Germany
| | - Oliver Ebenhöh
- Institute of Quantitative and Theoretical Biology, Heinrich Heine University, 40225, Düsseldorf, Germany.,Cluster of Excellence on Plant Sciences (CEPLAS), Heinrich Heine University, 40225, Düsseldorf, Germany
| | - Dario Gerace
- Department of Physics, University of Pavia, I-27100, Pavia, Italy.
| |
Collapse
|
19
|
Mulvihill E, Gao X, Liu Y, Schubert A, Dunietz BD, Geva E. Combining the mapping Hamiltonian linearized semiclassical approach with the generalized quantum master equation to simulate electronically nonadiabatic molecular dynamics. J Chem Phys 2019; 151:074103. [DOI: 10.1063/1.5110891] [Citation(s) in RCA: 22] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Ellen Mulvihill
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
| | - Xing Gao
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
| | - Yudan Liu
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
| | - Alexander Schubert
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
- Department of Chemistry and Biochemistry, Kent State University, Kent, Ohio 44242, USA
| | - Barry D. Dunietz
- Department of Chemistry and Biochemistry, Kent State University, Kent, Ohio 44242, USA
| | - Eitan Geva
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
| |
Collapse
|
20
|
Mulvihill E, Schubert A, Sun X, Dunietz BD, Geva E. A modified approach for simulating electronically nonadiabatic dynamics via the generalized quantum master equation. J Chem Phys 2019; 150:034101. [DOI: 10.1063/1.5055756] [Citation(s) in RCA: 26] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Ellen Mulvihill
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
| | - Alexander Schubert
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
- Department of Chemistry and Biochemistry, Kent State University, Kent, Ohio 44242, USA
| | - Xiang Sun
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
| | - Barry D. Dunietz
- Department of Chemistry and Biochemistry, Kent State University, Kent, Ohio 44242, USA
| | - Eitan Geva
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
| |
Collapse
|
21
|
Kananenka AA, Sun X, Schubert A, Dunietz BD, Geva E. A comparative study of different methods for calculating electronic transition rates. J Chem Phys 2018; 148:102304. [DOI: 10.1063/1.4989509] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Alexei A. Kananenka
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
| | - Xiang Sun
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
| | - Alexander Schubert
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
- Department of Chemistry and Biochemistry, Kent State University, Kent, Ohio 44242, USA
| | - Barry D. Dunietz
- Department of Chemistry and Biochemistry, Kent State University, Kent, Ohio 44242, USA
| | - Eitan Geva
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
| |
Collapse
|
22
|
Kananenka AA, Hsieh CY, Cao J, Geva E. Nonadiabatic Dynamics via the Symmetrical Quasi-Classical Method in the Presence of Anharmonicity. J Phys Chem Lett 2018; 9:319-326. [PMID: 29239614 DOI: 10.1021/acs.jpclett.7b03002] [Citation(s) in RCA: 28] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
The symmetrical quasi-classical (SQC) method recently proposed by Miller and Cotton allows one to simulate nonadiabatic dynamics based on an algorithm with classical-like scaling with respect to system size. This is made possible by casting the electronic degrees of freedom in terms of mapping variables that can be propagated in a classical-like manner. While SQC was shown to be rather accurate when applied to benchmark models with harmonic electronic potential energy surfaces, it was also found to become inaccurate and to suffer numerical instabilities when applied to anharmonic systems. In this paper, we propose an extended SQC (E-SQC) methodology for overcoming those discrepancies by describing the anharmonic nuclear modes, which are coupled to the electronic degrees of freedom, in terms of classical-like mapping variables. The accuracy of E-SQC relative to standard SQC is demonstrated on benchmark models with quartic and Morse potential energy surfaces.
Collapse
Affiliation(s)
- Alexei A Kananenka
- Department of Chemistry, University of Michigan , Ann Arbor, Michigan 48109, United States
| | - Chang-Yu Hsieh
- Department of Chemistry, Massachusetts Institute of Technology , Cambridge, Massachusetts 02319, United States
- Singapore-MIT Alliance for Research and Technology (SMART) Center , Singapore 487373
| | - Jianshu Cao
- Department of Chemistry, Massachusetts Institute of Technology , Cambridge, Massachusetts 02319, United States
- Singapore-MIT Alliance for Research and Technology (SMART) Center , Singapore 487373
| | | |
Collapse
|
23
|
Mendive-Tapia D, Mangaud E, Firmino T, de la Lande A, Desouter-Lecomte M, Meyer HD, Gatti F. Multidimensional Quantum Mechanical Modeling of Electron Transfer and Electronic Coherence in Plant Cryptochromes: The Role of Initial Bath Conditions. J Phys Chem B 2017; 122:126-136. [DOI: 10.1021/acs.jpcb.7b10412] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- David Mendive-Tapia
- Institut
Charles Gerhardt Montpellier, UMR 5253, CNRS-UM-ENSCM, CTMM, Université Montpellier, CC 15001, Place Eugène Bataillon, 34095 Montpellier, France
- Theoretische
Chemie, Physikalisch-Chemisches Institut, Universität Heidelberg, INF 229, D-69120 Heidelberg, Germany
| | - Etienne Mangaud
- Laboratoire
Collisions Agrégats Réactivité, UMR 5589, IRSAMC, Université Toulouse III Paul Sabatier, F-31062 Toulouse, France
| | - Thiago Firmino
- Laboratoire
de Chimie Physique, CNRS, Université Paris-Sud, Université Paris Saclay, Orsay F-91405, France
| | - Aurélien de la Lande
- Laboratoire
de Chimie Physique, CNRS, Université Paris-Sud, Université Paris Saclay, Orsay F-91405, France
| | - Michèle Desouter-Lecomte
- Laboratoire
de Chimie Physique, CNRS, Université Paris-Sud, Université Paris Saclay, Orsay F-91405, France
| | - Hans-Dieter Meyer
- Theoretische
Chemie, Physikalisch-Chemisches Institut, Universität Heidelberg, INF 229, D-69120 Heidelberg, Germany
| | - Fabien Gatti
- Institut
des Sciences Moléculaires d’Orsay, UMR-CNRS 8214, Université Paris-Sud, Université Paris Saclay, Orsay F-91405, France
| |
Collapse
|
24
|
Mangaud E, Meier C, Desouter-Lecomte M. Analysis of the non-Markovianity for electron transfer reactions in an oligothiophene-fullerene heterojunction. Chem Phys 2017. [DOI: 10.1016/j.chemphys.2017.07.011] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
|
25
|
Hait D, Mavros MG, Van Voorhis T. A hybrid memory kernel approach for condensed phase non-adiabatic dynamics. J Chem Phys 2017; 147:014108. [PMID: 28688393 DOI: 10.1063/1.4990739] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
The spin-boson model is a simplified Hamiltonian often used to study non-adiabatic dynamics in large condensed phase systems, even though it has not been solved in a fully analytic fashion. Herein, we present an exact analytic expression for the dynamics of the spin-boson model in the infinitely slow-bath limit and generalize it to approximate dynamics for faster baths. We achieve the latter by developing a hybrid approach that combines the exact slow-bath result with the popular non-interacting blip approximation (NIBA) method to generate a memory kernel that is formally exact to second-order in the diabatic coupling but also contains higher-order contributions approximated from the second-order term alone. This kernel has the same computational complexity as the NIBA, but is found to yield dramatically superior dynamics in regimes where the NIBA breaks down-such as systems with large diabatic coupling or energy bias. This indicates that this hybrid approach could be used to cheaply incorporate higher-order effects into second-order methods and could potentially be generalized to develop alternate kernel resummation schemes.
Collapse
Affiliation(s)
- Diptarka Hait
- Department of Chemistry, Massachusetts Institute of Technology, 77 Massachusetts Avenue, Cambridge, Massachusetts 02139, USA
| | - Michael G Mavros
- Department of Chemistry, Massachusetts Institute of Technology, 77 Massachusetts Avenue, Cambridge, Massachusetts 02139, USA
| | - Troy Van Voorhis
- Department of Chemistry, Massachusetts Institute of Technology, 77 Massachusetts Avenue, Cambridge, Massachusetts 02139, USA
| |
Collapse
|
26
|
Sun X, Geva E. Non-Condon nonequilibrium Fermi’s golden rule rates from the linearized semiclassical method. J Chem Phys 2016. [DOI: 10.1063/1.4960337] [Citation(s) in RCA: 26] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Xiang Sun
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109-1055, USA
| | - Eitan Geva
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109-1055, USA
| |
Collapse
|
27
|
Sun X, Geva E. Nonequilibrium Fermi’s Golden Rule Charge Transfer Rates via the Linearized Semiclassical Method. J Chem Theory Comput 2016; 12:2926-41. [DOI: 10.1021/acs.jctc.6b00236] [Citation(s) in RCA: 37] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Xiang Sun
- Department
of Chemistry, University of Michigan, Ann Arbor, Michigan 48109-1055, United States
| | - Eitan Geva
- Department
of Chemistry, University of Michigan, Ann Arbor, Michigan 48109-1055, United States
| |
Collapse
|
28
|
Sun X, Geva E. Exact vs. asymptotic spectral densities in the Garg-Onuchic-Ambegaokar charge transfer model and its effect on Fermi's golden rule rate constants. J Chem Phys 2016; 144:044106. [PMID: 26827201 DOI: 10.1063/1.4940308] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The Garg-Onuchic-Ambegaokar model [J. Chem. Phys. 83, 4491 (1985)] has been used extensively for benchmarking methods aimed at calculating charge transfer rates. Within this model, the donor and acceptor diabats are described as shifted parabolas along a single primary mode, which is bilinearly coupled to a harmonic bath consisting of secondary modes, characterized by an Ohmic spectral density with exponential cutoff. Rate calculations for this model are often performed in the normal mode representation, with the corresponding effective spectral density given by an asymptotic expression derived at the limit where the Ohmic bath cutoff frequency is much larger than the primary mode frequency. We compare Fermi's golden rule rate constants obtained with the asymptotic and exact effective spectral densities. We find significant deviations between rate constants obtained from the asymptotic spectral density and those obtained from the exact one in the deep inverted region. Within the range of primary mode frequencies commonly employed, we find that the discrepancies increase with decreasing temperature and with decreasing primary mode frequency.
Collapse
Affiliation(s)
- Xiang Sun
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109-1055, USA
| | - Eitan Geva
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109-1055, USA
| |
Collapse
|
29
|
Mangaud E, de la Lande A, Meier C, Desouter-Lecomte M. Electron transfer within a reaction path model calibrated by constrained DFT calculations: application to mixed-valence organic compounds. Phys Chem Chem Phys 2015; 17:30889-903. [PMID: 26041466 DOI: 10.1039/c5cp01194a] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The quantum dynamics of electron transfer in mixed-valence organic compounds is investigated using a reaction path model calibrated by constrained density functional theory (cDFT). Constrained DFT is used to define diabatic states relevant for describing the electron transfer, to obtain equilibrium structures for each of these states and to estimate the electronic coupling between them. The harmonic analysis at the diabatic minima yields normal modes forming the dissipative bath coupled to the electronic states. In order to decrease the system-bath coupling, an effective one dimensional vibronic Hamiltonian is constructed by partitioning the modes into a linear reaction path which connects both equilibrium positions and a set of secondary vibrational modes, coupled to this reaction coordinate. Using this vibronic model Hamiltonian, dissipative quantum dynamics is carried out using Redfield theory, based on a spectral density which is determined from the cDFT results. In a first benchmark case, the model is applied to a series of mixed-valence organic compounds formed by two 1,4-dimethoxy-3-methylphenylene fragments linked by an increasing number of phenylene bridges. This allows us to examine the coherent electron transfer in extreme situations leading to a ground adiabatic state with or without a barrier and therefore to the trapping of the charge or to an easy delocalization.
Collapse
Affiliation(s)
- E Mangaud
- Laboratoire Collisions Agrégats Réactivité, UMR 5589, IRSAMC, Université Toulouse III Paul Sabatier, Bât. 3R1b4, 118 route de Narbonne, F-31062, Toulouse, France. and Laboratoire de Chimie Physique, UMR 8000, Université Paris-Sud, Bât. 349, 15 avenue Jean Perrin, F-91405 Orsay, France
| | - A de la Lande
- Laboratoire de Chimie Physique, UMR 8000, Université Paris-Sud, Bât. 349, 15 avenue Jean Perrin, F-91405 Orsay, France
| | - C Meier
- Laboratoire Collisions Agrégats Réactivité, UMR 5589, IRSAMC, Université Toulouse III Paul Sabatier, Bât. 3R1b4, 118 route de Narbonne, F-31062, Toulouse, France.
| | - M Desouter-Lecomte
- Laboratoire de Chimie Physique, UMR 8000, Université Paris-Sud, Bât. 349, 15 avenue Jean Perrin, F-91405 Orsay, France and Département de Chimie, Université de Liège, Sart Tilman, B6, B-4000 Liège, Belgium
| |
Collapse
|
30
|
Chęcińska A, Pollock FA, Heaney L, Nazir A. Dissipation enhanced vibrational sensing in an olfactory molecular switch. J Chem Phys 2015; 142:025102. [PMID: 25591386 DOI: 10.1063/1.4905377] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022] Open
Abstract
Motivated by a proposed olfactory mechanism based on a vibrationally activated molecular switch, we study electron transport within a donor-acceptor pair that is coupled to a vibrational mode and embedded in a surrounding environment. We derive a polaron master equation with which we study the dynamics of both the electronic and vibrational degrees of freedom beyond previously employed semiclassical (Marcus-Jortner) rate analyses. We show (i) that in the absence of explicit dissipation of the vibrational mode, the semiclassical approach is generally unable to capture the dynamics predicted by our master equation due to both its assumption of one-way (exponential) electron transfer from donor to acceptor and its neglect of the spectral details of the environment; (ii) that by additionally allowing strong dissipation to act on the odorant vibrational mode, we can recover exponential electron transfer, though typically at a rate that differs from that given by the Marcus-Jortner expression; (iii) that the ability of the molecular switch to discriminate between the presence and absence of the odorant, and its sensitivity to the odorant vibrational frequency, is enhanced significantly in this strong dissipation regime, when compared to the case without mode dissipation; and (iv) that details of the environment absent from previous Marcus-Jortner analyses can also dramatically alter the sensitivity of the molecular switch, in particular, allowing its frequency resolution to be improved. Our results thus demonstrate the constructive role dissipation can play in facilitating sensitive and selective operation in molecular switch devices, as well as the inadequacy of semiclassical rate equations in analysing such behaviour over a wide range of parameters.
Collapse
Affiliation(s)
- Agata Chęcińska
- Centre for Quantum Technologies, National University of Singapore, Singapore 117543, Singapore
| | - Felix A Pollock
- Atomic and Laser Physics, Clarendon Laboratory, University of Oxford, Parks Road, Oxford OX1 3PU, United Kingdom
| | - Libby Heaney
- Centre for Quantum Technologies, National University of Singapore, Singapore 117543, Singapore
| | - Ahsan Nazir
- Photon Science Institute and School of Physics and Astronomy, University of Manchester, Oxford Road, Manchester M13 9PL, United Kingdom
| |
Collapse
|
31
|
Nesterov AI, Berman GP. Decoherence and spin echo in biological systems. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2015; 91:052702. [PMID: 26066194 DOI: 10.1103/physreve.91.052702] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/08/2015] [Indexed: 06/04/2023]
Abstract
The spin-echo approach is extended to include biocomplexes for which the interaction with dynamical noise, produced by the protein environment, is strong. Significant restoration of the free induction decay signal due to homogeneous (decoherence) and inhomogeneous (dephasing) broadening is demonstrated analytically and numerically for both an individual dimer of interacting chlorophylls and for an ensemble of dimers. Our approach does not require the use of small interaction constants between the electron states and the protein fluctuations. It is based on an exact and closed system of ordinary differential equations that can be easily solved for a wide range of parameters that are relevant for bioapplications.
Collapse
Affiliation(s)
- Alexander I Nesterov
- Departamento de Física, CUCEI, Universidad de Guadalajara, Av. Revolución 1500, Guadalajara, CP 44420, Jalisco, México
| | - Gennady P Berman
- Theoretical Division, T-4, Los Alamos National Laboratory, and the New Mexico Consortium, Los Alamos, New Mexico 87544, USA
| |
Collapse
|
32
|
Nesterov AI, Berman GP. Role of protein fluctuation correlations in electron transfer in photosynthetic complexes. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2015; 91:042702. [PMID: 25974522 DOI: 10.1103/physreve.91.042702] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/25/2014] [Indexed: 06/04/2023]
Abstract
We consider the dependence of the electron transfer in photosynthetic complexes on correlation properties of random fluctuations of the protein environment. The electron subsystem is modeled by a finite network of connected electron (exciton) sites. The fluctuations of the protein environment are modeled by random telegraph processes, which act either collectively (correlated) or independently (uncorrelated) on the electron sites. We derived an exact closed system of first-order linear differential equations with constant coefficients, for the average density matrix elements and for their first moments. Under some conditions, we obtained analytic expressions for the electron transfer rates and found the range of parameters for their applicability by comparing with the exact numerical simulations. We also compared the correlated and uncorrelated regimes and demonstrated numerically that the uncorrelated fluctuations of the protein environment can, under some conditions, either increase or decrease the electron transfer rates.
Collapse
Affiliation(s)
- Alexander I Nesterov
- Departamento de Física, CUCEI, Universidad de Guadalajara, Avenida Revolución 1500, Guadalajara, Codigo Postal 44420, Jalisco, Mexico
| | - Gennady P Berman
- Theoretical Division, T-4, Los Alamos National Laboratory and the New Mexico Consortium, Los Alamos, New Mexico 87544, USA
| |
Collapse
|
33
|
|
34
|
Borrelli R, Peluso A. Elementary electron transfer reactions: from basic concepts to recent computational advances. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2013. [DOI: 10.1002/wcms.1147] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
|
35
|
Ritz T, Schulten K. Physik der Photosynthese: Wie Bakterien die Quantenphysik ausnutzen, um effizient Photosynthese zu betreiben. ACTA ACUST UNITED AC 2013. [DOI: 10.1002/phbl.20010570212] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
|
36
|
Skourtis SS. Reviewprobing protein electron transfer mechanisms from the molecular to the cellular length scales. Biopolymers 2013; 100:82-92. [DOI: 10.1002/bip.22169] [Citation(s) in RCA: 31] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/30/2012] [Revised: 09/14/2012] [Accepted: 09/23/2012] [Indexed: 11/10/2022]
|
37
|
Solov'yov IA, Chang PY, Schulten K. Vibrationally assisted electron transfer mechanism of olfaction: myth or reality? Phys Chem Chem Phys 2012; 14:13861-71. [PMID: 22899100 DOI: 10.1039/c2cp41436h] [Citation(s) in RCA: 48] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Smell is a vital sense for animals. The mainstream explanation of smell is based on recognition of the odorant molecules through characteristics of their surface, e.g., shape, but certain experiments suggest that such recognition is complemented by recognition of vibrational modes. According to this suggestion an olfactory receptor is activated by electron transfer assisted through odorant vibrational excitation. The hundreds to thousands of different olfactory receptors in an animal recognize odorants over a discriminant landscape with surface properties and vibrational frequencies as the two major dimensions. In the present paper we introduce the vibrationally assisted mechanism of olfaction and demonstrate for several odorants that, indeed, a strong enhancement of an electron tunneling rate due to odorant vibrations can arise. We discuss in this regard the influence of odorant deuteration and explain, thereby, recent experiments performed on Drosophila melanogaster. Our demonstration is based on known physical properties of biological electron transfer and on ab initio calculations on odorants carried out for the purpose of the present study. We identify a range of physical characteristics which olfactory receptors and odorants must obey for the vibrationally assisted electron transfer mechanism to function. We argue that the stated characteristics are feasible for realistic olfactory receptors, noting, though, that the receptor structure presently is still unknown, but can be studied through homology modeling.
Collapse
Affiliation(s)
- Ilia A Solov'yov
- Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, 405 N. Mathews Ave., Urbana, Illinois 61801, USA.
| | | | | |
Collapse
|
38
|
Strümpfer J, Schulten K. Open Quantum Dynamics Calculations with the Hierarchy Equations of Motion on Parallel Computers. J Chem Theory Comput 2012; 8:2808-2816. [PMID: 23105920 PMCID: PMC3480185 DOI: 10.1021/ct3003833] [Citation(s) in RCA: 118] [Impact Index Per Article: 9.8] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Abstract
Calculating the evolution of an open quantum system, i.e., a system in contact with a thermal environment, has presented a theoretical and computational challenge for many years. With the advent of supercomputers containing large amounts of memory and many processors, the computational challenge posed by the previously intractable theoretical models can now be addressed. The hierarchy equations of motion present one such model and offer a powerful method that remained under-utilized so far due to its considerable computational expense. By exploiting concurrent processing on parallel computers the hierarchy equations of motion can be applied to biological-scale systems. Herein we introduce the quantum dynamics software PHI, that solves the hierarchical equations of motion. We describe the integrator employed by PHI and demonstrate PHI's scaling and efficiency running on large parallel computers by applying the software to the calculation of inter-complex excitation transfer between the light harvesting complexes 1 and 2 of purple photosynthetic bacteria, a 50 pigment system.
Collapse
Affiliation(s)
- Johan Strümpfer
- Center for Biophysics and Computational Biology, University of Illinois at Urbana-Champaign
| | - Klaus Schulten
- Center for Biophysics and Computational Biology, University of Illinois at Urbana-Champaign
- Department of Physics and Beckman Institute, University of Illinois at Urbana-Champaign
| |
Collapse
|
39
|
Li Y, Lagowski JB. Charge carrier mobility in conjugated organic polymers – Case studies using multi-step computational approach. POLYMER 2011. [DOI: 10.1016/j.polymer.2011.08.011] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/17/2022]
|
40
|
Şener M, Strümpfer J, Hsin J, Chandler D, Scheuring S, Hunter CN, Schulten K. Förster energy transfer theory as reflected in the structures of photosynthetic light-harvesting systems. Chemphyschem 2011; 12:518-31. [PMID: 21344591 DOI: 10.1002/cphc.201000944] [Citation(s) in RCA: 101] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
Abstract
Förster’s theory of resonant energy transfer underlies a fundamental process in nature, namely the harvesting of sunlight by photosynthetic life forms. The theoretical framework developed by Förster and others describes how electronic excitation migrates in the photosynthetic apparatus of plants, algae, and bacteria from light absorbing pigments to reaction centers where light energy is utilized for the eventual conversion into chemical energy. The demand for highest possible efficiency of light harvesting appears to have shaped the evolution of photosynthetic species from bacteria to plants which, despite a great variation in architecture, display common structural themes founded on the quantum physics of energy transfer as described first by Förster. Herein, Förster’s theory of excitation transfer is summarized, including recent extensions, and the relevance of the theory to photosynthetic systems as evolved in purple bacteria, cyanobacteria, and plants is demonstrated. Förster’s energy transfer formula, as used widely today in many fields of science, is also derived.
Collapse
Affiliation(s)
- Melih Şener
- Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
| | | | | | | | | | | | | |
Collapse
|
41
|
Strümpfer J, Schulten K. The effect of correlated bath fluctuations on exciton transfer. J Chem Phys 2011; 134:095102. [PMID: 21385000 PMCID: PMC3064689 DOI: 10.1063/1.3557042] [Citation(s) in RCA: 86] [Impact Index Per Article: 6.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/01/2010] [Accepted: 02/01/2011] [Indexed: 01/27/2023] Open
Abstract
Excitation dynamics of various light harvesting systems have been investigated with many theoretical methods including various non-Markovian descriptions of dissipative quantum dynamics. It is typically assumed that each excited state is coupled to an independent thermal environment, i.e., that fluctuations in different environments are uncorrelated. Here the assumption is dropped and the effect of correlated bath fluctuations on excitation transfer is investigated. Using the hierarchy equations of motion for dissipative quantum dynamics it is shown for models of the B850 bacteriochlorophylls of LH2 that correlated bath fluctuations have a significant effect on the LH2→LH2 excitation transfer rate. It is also demonstrated that inclusion of static disorder is crucial for an accurate description of transfer dynamics.
Collapse
Affiliation(s)
- Johan Strümpfer
- Center for Biophysics and Computational Biology, University of Illinois at Urbana-Champaign, Urbana, USA
| | | |
Collapse
|
42
|
Borrelli R, Peluso A. The temperature dependence of radiationless transition rates from ab initio computations. Phys Chem Chem Phys 2011; 13:4420-6. [DOI: 10.1039/c0cp02307h] [Citation(s) in RCA: 50] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
43
|
Paulo PMR, Lopes JNC, Costa SMB. Molecular dynamics simulations of porphyrin-dendrimer systems: toward modeling electron transfer in solution. J Phys Chem B 2009; 112:14779-92. [PMID: 18954105 DOI: 10.1021/jp806849y] [Citation(s) in RCA: 14] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
We have performed computational simulations of porphyrin-dendrimer systems--a cationic porphyrin electrostatically associated to a negatively charged dendrimer--using the method of classical molecular dynamics (MD) with an atomistic force field. Previous experimental studies have shown a strong quenching effect of the porphyrin fluorescence that was assigned to electron transfer (ET) from the dendrimer's tertiary amines (Paulo, P. M. R.; Costa, S. M. B. J. Phys. Chem. B 2005, 109, 13928). In the present contribution, we evaluate computationally the role of the porphyrin-dendrimer conformation in the development of a statistical distribution of ET rates through its dependence on the donor-acceptor distance. We started from simulations without explicit solvent to obtain trajectories of the donor-acceptor distance and the respective time-averaged distributions for two dendrimer sizes and different initial configurations of the porphyrin-dendrimer pair. By introducing explicit solvent (water) in our simulations, we were able to estimate the reorganization energy of the medium for the systems with the dendrimer of smaller size. The values obtained are in the range 0.6-1.5 eV and show a linear dependence with the inverse of the donor-acceptor distance, which can be explained by a two-phase dielectric continuum model taking into account the medium heterogeneity provided by the dendrimer organic core. Dielectric relaxation accompanying ET was evaluated from the simulations with explicit solvent showing fast decay times of some tens of femtoseconds and slow decay times in the range of hundreds of femtoseconds to a few picoseconds. The variations of the slow relaxation times reflect the heterogeneity of the dendrimer donor sites which add to the complexity of ET kinetics as inferred from the experimental fluorescence decays.
Collapse
Affiliation(s)
- Pedro M R Paulo
- Centro de Química Estrutural-Complexo I, Instituto Superior Técnico, Lisboa, Portugal.
| | | | | |
Collapse
|
44
|
Sener M, Hsin J, Trabuco LG, Villa E, Qian P, Hunter CN, Schulten K. Structural model and excitonic properties of the dimeric RC-LH1-PufX complex from Rhodobacter sphaeroides. Chem Phys 2009; 357:188-197. [PMID: 20161332 DOI: 10.1016/j.chemphys.2009.01.003] [Citation(s) in RCA: 42] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
The light-harvesting apparatus of the purple bacterial photosynthetic unit consists of a pool of peripheral light-harvesting complexes that transfer excitation energy to a reaction center (RC) via the surrounding pigment-protein complex LH1. Recent electron microscopy and atomic force microscopy studies have revealed that RC-LH1 units of Rhodobacter sphaeroides form membrane-bending dimeric complexes together with the polypeptide PufX. We present a structural model for these RC-LH1-PufX dimeric complexes constructed using the molecular dynamics flexible fitting method based on an EM density map. The arrangement of the LH1 BChls displays a distortion near the proposed location of the PufX polypeptide. The resulting atomic model for BChl arrays is used to compute the excitonic properties of the dimeric RC-LH1 complex. A comparison is presented between the structural and excitonic features of the S-shaped dimeric BChl array of Rhodobacter sphaeroides and the circular BChl arrangement found in other purple bacteria.
Collapse
Affiliation(s)
- Melih Sener
- Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801
| | | | | | | | | | | | | |
Collapse
|
45
|
From Atomic-Level Structure to Supramolecular Organization in the Photosynthetic Unit of Purple Bacteria. THE PURPLE PHOTOTROPHIC BACTERIA 2009. [DOI: 10.1007/978-1-4020-8815-5_15] [Citation(s) in RCA: 14] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/01/2022]
|
46
|
Gilmore J, McKenzie RH. Quantum Dynamics of Electronic Excitations in Biomolecular Chromophores: Role of the Protein Environment and Solvent. J Phys Chem A 2008; 112:2162-76. [DOI: 10.1021/jp710243t] [Citation(s) in RCA: 63] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2023]
Affiliation(s)
- Joel Gilmore
- Department of Physics, University of Queensland, Brisbane 4072 Australia
| | - Ross H. McKenzie
- Department of Physics, University of Queensland, Brisbane 4072 Australia
| |
Collapse
|
47
|
Kosztin I, Schulten K. Molecular Dynamics Methods for Bioelectronic Systems in Photosynthesis. ACTA ACUST UNITED AC 2008. [DOI: 10.1007/978-1-4020-8250-4_22] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/07/2023]
|
48
|
Solov'yov IA, Chandler DE, Schulten K. Magnetic field effects in Arabidopsis thaliana cryptochrome-1. Biophys J 2007; 92:2711-26. [PMID: 17259272 PMCID: PMC1831705 DOI: 10.1529/biophysj.106.097139] [Citation(s) in RCA: 142] [Impact Index Per Article: 8.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/28/2006] [Accepted: 01/04/2007] [Indexed: 11/18/2022] Open
Abstract
The ability of some animals, most notably migratory birds, to sense magnetic fields is still poorly understood. It has been suggested that this "magnetic sense" may be mediated by the blue light receptor protein cryptochrome, which is known to be localized in the retinas of migratory birds. Cryptochromes are a class of photoreceptor signaling proteins that are found in a wide variety of organisms and that primarily perform regulatory functions, such as the entrainment of circadian rhythm in mammals and the inhibition of hypocotyl growth in plants. Recent experiments have shown that the activity of cryptochrome-1 in Arabidopsis thaliana is enhanced by the presence of a weak external magnetic field, confirming the ability of cryptochrome to mediate magnetic field responses. Cryptochrome's signaling is tied to the photoreduction of an internally bound chromophore, flavin adenine dinucleotide. The spin chemistry of this photoreduction process, which involves electron transfer from a chain of three tryptophans, can be modulated by the presence of a magnetic field in an effect known as the radical-pair mechanism. Here we present and analyze a model of the flavin-adenine-dinucleotide-tryptophan chain system that incorporates realistic hyperfine coupling constants and reaction rate constants. Our calculations show that the radical-pair mechanism in cryptochrome can produce an increase in the protein's signaling activity of approximately 10% for magnetic fields on the order of 5 G, which is consistent with experimental results. These calculations, in view of the similarity between bird and plant cryptochromes, provide further support for a cryptochrome-based model of avian magnetoreception.
Collapse
Affiliation(s)
- Ilia A Solov'yov
- Frankfurt Institute for Advanced Studies, Johann Wolfgang Goethe University, Frankfurt am Main, Germany
| | | | | |
Collapse
|
49
|
Pollard WT, Felts AK, Friesner RA. The Redfield Equation in Condensed-Phase Quantum Dynamics. ADVANCES IN CHEMICAL PHYSICS 2007. [DOI: 10.1002/9780470141526.ch3] [Citation(s) in RCA: 124] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
|
50
|
Tornow S, Tong NH, Bulla R. Dissipative exciton transfer in donor-bridge-acceptor systems: numerical renormalization group calculation of equilibrium properties. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2006; 18:5985-6000. [PMID: 21690813 DOI: 10.1088/0953-8984/18/26/017] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/30/2023]
Abstract
We present a detailed model study of exciton transfer processes in donor-bridge-acceptor (DBA) systems. Using a model which includes the intermolecular Coulomb interaction and the coupling to a dissipative environment we calculate the phase diagram, the absorption spectrum as well as dynamic equilibrium properties with the numerical renormalization group. This method is non-perturbative and therefore allows one to cover the full parameter space, especially the case when the intermolecular Coulomb interaction is of the same order as the coupling to the environment and perturbation theory cannot be applied. For DBA systems with up to six sites we found a transition to the localized phase (self-trapping) depending on the coupling to the dissipative environment. We discuss various criteria which favour delocalized exciton transfer.
Collapse
Affiliation(s)
- Sabine Tornow
- Theoretische Physik III, Elektronische Korrelationen und Magnetismus, Universität Augsburg, 86135 Augsburg, Germany
| | | | | |
Collapse
|