1
|
Kim HW, Rhee YM. Two‐dimensional electronic spectrum simulation of simple photosynthetic complex models with semi‐classical Poisson bracket mapping equation. B KOREAN CHEM SOC 2022. [DOI: 10.1002/bkcs.12469] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Affiliation(s)
- Hyun Woo Kim
- Center for Molecular Modeling and Simulation, Korea Research Institute of Chemical Technology (KRICT) Daejeon South Korea
| | - Young Min Rhee
- Department of Chemistry Korea Advanced Institute of Science and Technology (KAIST) Daejeon South Korea
| |
Collapse
|
2
|
Braver Y, Valkunas L, Gelzinis A. Stark absorption and Stark fluorescence spectroscopies: Theory and simulations. J Chem Phys 2021; 155:244101. [PMID: 34972359 DOI: 10.1063/5.0073962] [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
Stark spectroscopy experiments are widely used to study the properties of molecular systems, particularly those containing charge-transfer (CT) states. However, due to the small transition dipole moments and large static dipole moments of the CT states, the standard interpretation of the Stark absorption and Stark fluorescence spectra in terms of the Liptay model may be inadequate. In this work, we provide a theoretical framework for calculations of Stark absorption and Stark fluorescence spectra and propose new methods of simulations that are based on the quantum-classical theory. In particular, we use the forward-backward trajectory solution and a variant of the Poisson bracket mapping equation, which have been recently adapted for the calculation of conventional (field-free) absorption and fluorescence spectra. For comparison, we also apply the recently proposed complex time-dependent Redfield theory, while exact results are obtained using the hierarchical equations of motion approach. We show that the quantum-classical methods produce accurate results for a wide range of systems, including those containing CT states. The CT states contribute significantly to the Stark spectra, and the standard Liptay formalism is shown to be inapplicable for the analysis of spectroscopic data in those cases. We demonstrate that states with large static dipole moments may cause a pronounced change in the total fluorescence yield of the system in the presence of an external electric field. This effect is correctly captured by the quantum-classical methods, which should therefore prove useful for further studies of Stark spectra of real molecular systems. As an example, we calculate the Stark spectra for the Fenna-Matthews-Olson complex of green sulfur bacteria.
Collapse
Affiliation(s)
- Yakov Braver
- Faculty of Physics, Institute of Chemical Physics, Vilnius University, Saulėtekio Ave. 9-III, LT-10222 Vilnius, Lithuania and Department of Molecular Compound Physics, Center for Physical Sciences and Technology, Saulėtekio Ave. 3, LT-10257 Vilnius, Lithuania
| | - Leonas Valkunas
- Faculty of Physics, Institute of Chemical Physics, Vilnius University, Saulėtekio Ave. 9-III, LT-10222 Vilnius, Lithuania and Department of Molecular Compound Physics, Center for Physical Sciences and Technology, Saulėtekio Ave. 3, LT-10257 Vilnius, Lithuania
| | - Andrius Gelzinis
- Faculty of Physics, Institute of Chemical Physics, Vilnius University, Saulėtekio Ave. 9-III, LT-10222 Vilnius, Lithuania and Department of Molecular Compound Physics, Center for Physical Sciences and Technology, Saulėtekio Ave. 3, LT-10257 Vilnius, Lithuania
| |
Collapse
|
3
|
Braver Y, Valkunas L, Gelzinis A. Quantum-Classical Approach for Calculations of Absorption and Fluorescence: Principles and Applications. J Chem Theory Comput 2021; 17:7157-7168. [PMID: 34618457 PMCID: PMC8719324 DOI: 10.1021/acs.jctc.1c00777] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/03/2021] [Indexed: 01/20/2023]
Abstract
Absorption and fluorescence spectroscopy techniques provide a wealth of information on molecular systems. The simulations of such experiments remain challenging, however, despite the efforts put into developing the underlying theory. An attractive method of simulating the behavior of molecular systems is provided by the quantum-classical theory─it enables one to keep track of the state of the bath explicitly, which is needed for accurate calculations of fluorescence spectra. Unfortunately, until now there have been relatively few works that apply quantum-classical methods for modeling spectroscopic data. In this work, we seek to provide a framework for the calculations of absorption and fluorescence lineshapes of molecular systems using the methods based on the quantum-classical Liouville equation, namely, the forward-backward trajectory solution (FBTS) and the non-Hamiltonian variant of the Poisson bracket mapping equation (PBME-nH). We perform calculations on a molecular dimer and the photosynthetic Fenna-Matthews-Olson complex. We find that in the case of absorption, the FBTS outperforms PBME-nH, consistently yielding highly accurate results. We next demonstrate that for fluorescence calculations, the method of choice is a hybrid approach, which we call PBME-nH-Jeff, that utilizes the effective coupling theory [Gelzinis, A.; J. Chem. Phys. 2020, 152, 051103] to estimate the excited state equilibrium density operator. Thus, we find that FBTS and PBME-nH-Jeff are excellent candidates for simulating, respectively, absorption and fluorescence spectra of real molecular systems.
Collapse
Affiliation(s)
- Yakov Braver
- Institute
of Chemical Physics, Faculty of Physics, Vilnius University, Saulėtekio Avenue 9-III, LT-10222 Vilnius, Lithuania
- Department
of Molecular Compound Physics, Center for
Physical Sciences and Technology, Saulėtekio Avenue 3, LT-10257 Vilnius, Lithuania
| | - Leonas Valkunas
- Institute
of Chemical Physics, Faculty of Physics, Vilnius University, Saulėtekio Avenue 9-III, LT-10222 Vilnius, Lithuania
- Department
of Molecular Compound Physics, Center for
Physical Sciences and Technology, Saulėtekio Avenue 3, LT-10257 Vilnius, Lithuania
| | - Andrius Gelzinis
- Institute
of Chemical Physics, Faculty of Physics, Vilnius University, Saulėtekio Avenue 9-III, LT-10222 Vilnius, Lithuania
- Department
of Molecular Compound Physics, Center for
Physical Sciences and Technology, Saulėtekio Avenue 3, LT-10257 Vilnius, Lithuania
| |
Collapse
|
4
|
Abstract
We introduce a simple and effective method to decompose the energy dissipation in the dynamics of open quantum systems into contributions due to individual bath components. The method is based on a vibronic extension of the Förster resonance energy transfer theory that enables quantifying the energy dissipated by specific bath degrees of freedom. Its accuracy is determined by benchmarking against mixed quantum-classical simulations that reveal that the method provides a semi-quantitative frequency-dependent decomposition of the overall dissipation. The utility of the method is illustrated by using a model donor-acceptor pair interacting to a thermal harmonic bath with different coupling strengths. The method can be used to identify the key features of a bath that leads to energy dissipation as required to develop a deep understanding of the dynamics of open quantum systems and to engineer environments with desired dissipative features.
Collapse
Affiliation(s)
- Chang Woo Kim
- Department of Chemistry, University of Rochester, Rochester, New York 14627, USA
| | - Ignacio Franco
- Department of Chemistry, University of Rochester, Rochester, New York 14627, USA
| |
Collapse
|
5
|
Kim HW, Rhee YM. Two-oscillator mapping modification of the Poisson bracket mapping equation formulation of the quantum-classical Liouville equation. J Chem Phys 2020; 153:214103. [PMID: 33291890 DOI: 10.1063/5.0027799] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022] Open
Abstract
Mapping basis solutions provide efficient ways for simulating mixed quantum-classical (MQC) dynamics in complex systems by matching multiple quantum states of interest to some fictitious physical states. Recently, various MQC methods were devised such that two harmonic oscillators are employed to represent each electronic state, showing improvements over one-oscillator-based methods. Here, we introduce and analyze newly modified mapping approximations of the quantum-classical Liouville equation (QCLE) using two oscillators for each electronic state. We design two separate mapping relations that we can adopt toward simulating dynamics and computing expectation values. Through the process, two MQC methods can be constructed, one of which actually reproduces the population dynamics of the forward and backward trajectory solution of QCLE. By applying the methods to spin-boson systems with a range of parameters, we find out that the choice of mapping relations greatly affects the simulation results. We also show that further improvement is possible through using modified identity operator formulations. Our findings may be helpful in constructing improved MQC methods in the future.
Collapse
Affiliation(s)
- Hyun Woo Kim
- Center for Molecular Modeling and Simulation, Korea Research Institute of Chemical Technology (KRICT), Daejeon 34114, South Korea
| | - Young Min Rhee
- Department of Chemistry, Korea Advanced Institute of Science and Technology (KAIST), Daejeon 34141, South Korea
| |
Collapse
|
6
|
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
|
7
|
Braver Y, Valkunas L, Gelzinis A. Benchmarking the forward–backward trajectory solution of the quantum-classical Liouville equation. J Chem Phys 2020; 152:214116. [DOI: 10.1063/5.0006538] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Yakov Braver
- Institute of Chemical Physics, Faculty of Physics, Vilnius University, Saulėtekio Ave. 9-III, LT-10222 Vilnius, Lithuania
- Department of Molecular Compound Physics, Center for Physical Sciences and Technology, Saulėtekio Ave. 3, LT-10257 Vilnius, Lithuania
| | - Leonas Valkunas
- Institute of Chemical Physics, Faculty of Physics, Vilnius University, Saulėtekio Ave. 9-III, LT-10222 Vilnius, Lithuania
- Department of Molecular Compound Physics, Center for Physical Sciences and Technology, Saulėtekio Ave. 3, LT-10257 Vilnius, Lithuania
| | - Andrius Gelzinis
- Institute of Chemical Physics, Faculty of Physics, Vilnius University, Saulėtekio Ave. 9-III, LT-10222 Vilnius, Lithuania
- Department of Molecular Compound Physics, Center for Physical Sciences and Technology, Saulėtekio Ave. 3, LT-10257 Vilnius, Lithuania
| |
Collapse
|
8
|
Liu J, Hanna G. Efficient and Deterministic Propagation of Mixed Quantum-Classical Liouville Dynamics. J Phys Chem Lett 2018; 9:3928-3933. [PMID: 29947224 DOI: 10.1021/acs.jpclett.8b01619] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
We propose a highly efficient mixed quantum-classical molecular dynamics scheme based on a solution of the quantum-classical Liouville equation (QCLE). By casting the equations of motion for the quantum subsystem and classical bath degrees of freedom onto an approximate set of coupled first-order differential equations for c-numbers, this scheme propagates the composite system in time deterministically in terms of independent classical-like trajectories. To demonstrate its performance, we apply the method to the spin-boson model, a photoinduced electron transfer model, and a Fenna-Matthews-Olsen complex model, and find excellent agreement out to long times with the numerically exact results, using several orders of magnitude fewer trajectories than surface-hopping solutions of the QCLE. Owing to its accuracy and efficiency, this method promises to be very useful for studying the dynamics of mixed quantum-classical systems.
Collapse
Affiliation(s)
- Junjie Liu
- Department of Chemistry , University of Alberta , Edmonton , Alberta T6G 2G2 , Canada
| | - Gabriel Hanna
- Department of Chemistry , University of Alberta , Edmonton , Alberta T6G 2G2 , Canada
| |
Collapse
|
9
|
Freedman H, Hanna G. Mixed quantum–classical Liouville simulation of vibrational energy transfer in a model α-helix at 300 K. Chem Phys 2016. [DOI: 10.1016/j.chemphys.2016.08.015] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
10
|
Lee MK, Huo P, Coker DF. Semiclassical Path Integral Dynamics: Photosynthetic Energy Transfer with Realistic Environment Interactions. Annu Rev Phys Chem 2016; 67:639-68. [DOI: 10.1146/annurev-physchem-040215-112252] [Citation(s) in RCA: 96] [Impact Index Per Article: 12.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Affiliation(s)
- Mi Kyung Lee
- Department of Chemistry, Boston University, Boston, Massachusetts 02215;
| | - Pengfei Huo
- Department of Chemistry, University of Rochester, Rochester, New York 14627;
| | - David F. Coker
- Department of Chemistry, Boston University, Boston, Massachusetts 02215;
| |
Collapse
|
11
|
Martinez F, Hanna G. Mixed Quantum-Classical Simulations of Transient Absorption Pump–Probe Signals for a Photo-Induced Electron Transfer Reaction Coupled to an Inner-Sphere Vibrational Mode. J Phys Chem A 2016; 120:3196-205. [DOI: 10.1021/acs.jpca.5b11727] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Franz Martinez
- Department
of Chemistry, University of Alberta, Edmonton, Alberta T6G 2R3, Canada
| | - Gabriel Hanna
- Department
of Chemistry, University of Alberta, Edmonton, Alberta T6G 2R3, Canada
| |
Collapse
|