1
|
Toutounji M. Mixed Quantum-Classical Liouville Equation Treatment of Electronic Spectroscopy of Condensed Systems: Harmonic and Anharmonic Electron-Phonon Coupling. J Chem Theory Comput 2023. [PMID: 37365487 DOI: 10.1021/acs.jctc.2c01225] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/28/2023]
Abstract
This Review integrates the use of electronic optical response function theory and the mixed quantum-classical (MQC) Liouville equation (MQCLE), thereby leading to electronic spectroscopy in MQC media. It further sheds light on the applicability, utility, and efficiency of the mixed quantum-classical dynamics (MQCD) formalism, which starts off with the MQCLE, in probing spectroscopy and dynamics of condensed systems, whereby quantum mechanics and classical mechanics are combined systematically. The author has been exploring and implementing MQCD to investigate electron-phonon coupling effects on electronic dephasing in harmonic and anharmonic systems by calculating linear and nonlinear optical transition analytically and numerically dipole moment time correlation functions in an MQC environment, thereby presenting an in depth spectral profile analysis and their shape and symmetry. The distinctive capability of the MQC time correlation functions is that ergodicity and stationarity properties are inherently satisfied as part of the mixed quantum-classical dynamics (MQCD) framework, unlike classical correlation functions. While some research groups have applied MQCLE to calculate vibrational spectra to study hydrogen-bonded complexes in a MQC environment and other groups calculated Optical response function to probe electron transfer dynamics using the basis mapping technique, the approach, purpose, rigor, applications, and path to the end results reported herein are different. Finally, the same framework is employed to study dissipative systems in the MQC limit, whereby the zero-phonon line adopts the correct width and eliminates its asymmetry. While the full quantum mechanical model, like the multimode Brownian oscillator (MBO) model, yields the correct width and inaccurate shape in the low-temperature limit, the MQCD formalism seems to produce an accurate zero-phonon profile. Nonlinear optical signals are also reviewed in MQC media to show the applicability and utility of this approach. The vibronic optical response functions developed here will account for geometry change, frequency change, and anharmonicity upon electronic excitation to accurately probe electronic dephasing, electron-phonon coupling, shape, and symmetry of profiles and present differences and similarities to the MBO model on pure electronic dephasing. Frequency change and anharmonicity are vitally crucial for accurately assessing electron-phonon coupling upon electronic excitation. This is an additional unique result obtained by the author to further demonstrate the applicability and utility of this approach over other approximation schemes in probing electronic dephasing including that of the MBO model.
Collapse
Affiliation(s)
- Mohamad Toutounji
- College of Science, Department of Chemistry, P.O. Box 15551, UAE University, Al-Ain, 33000, UAE
| |
Collapse
|
2
|
Gelin MF, Chen L, Domcke W. Equation-of-Motion Methods for the Calculation of Femtosecond Time-Resolved 4-Wave-Mixing and N-Wave-Mixing Signals. Chem Rev 2022; 122:17339-17396. [PMID: 36278801 DOI: 10.1021/acs.chemrev.2c00329] [Citation(s) in RCA: 20] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023]
Abstract
Femtosecond nonlinear spectroscopy is the main tool for the time-resolved detection of photophysical and photochemical processes. Since most systems of chemical interest are rather complex, theoretical support is indispensable for the extraction of the intrinsic system dynamics from the detected spectroscopic responses. There exist two alternative theoretical formalisms for the calculation of spectroscopic signals, the nonlinear response-function (NRF) approach and the spectroscopic equation-of-motion (EOM) approach. In the NRF formalism, the system-field interaction is assumed to be sufficiently weak and is treated in lowest-order perturbation theory for each laser pulse interacting with the sample. The conceptual alternative to the NRF method is the extraction of the spectroscopic signals from the solutions of quantum mechanical, semiclassical, or quasiclassical EOMs which govern the time evolution of the material system interacting with the radiation field of the laser pulses. The NRF formalism and its applications to a broad range of material systems and spectroscopic signals have been comprehensively reviewed in the literature. This article provides a detailed review of the suite of EOM methods, including applications to 4-wave-mixing and N-wave-mixing signals detected with weak or strong fields. Under certain circumstances, the spectroscopic EOM methods may be more efficient than the NRF method for the computation of various nonlinear spectroscopic signals.
Collapse
Affiliation(s)
- Maxim F Gelin
- School of Science, Hangzhou Dianzi University, Hangzhou 310018, China
| | - Lipeng Chen
- Max-Planck-Institut für Physik komplexer Systeme, Nöthnitzer Strasse 38, D-01187 Dresden, Germany
| | - Wolfgang Domcke
- Department of Chemistry, Technical University of Munich, D-85747 Garching,Germany
| |
Collapse
|
3
|
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
|
4
|
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
|
5
|
Provazza J, Segatta F, Coker DF. Modeling Nonperturbative Field-Driven Vibronic Dynamics: Selective State Preparation and Nonlinear Spectroscopy. J Chem Theory Comput 2020; 17:29-39. [DOI: 10.1021/acs.jctc.0c01035] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Justin Provazza
- Department of Chemistry, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States
| | - Francesco Segatta
- Department of Chemistry, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States
- Dipartimento di Chimica Industriale “Toso Montanari”, University of Bologna, Viale del Risorgimento, 4, 40136 Bologna, Italy
| | - David F. Coker
- Department of Chemistry, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States
| |
Collapse
|
6
|
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
|
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
|
Quasi-Lie Brackets and the Breaking of Time-Translation Symmetry for Quantum Systems Embedded in Classical Baths. Symmetry (Basel) 2018. [DOI: 10.3390/sym10100518] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022] Open
Abstract
Many open quantum systems encountered in both natural and synthetic situations are embedded in classical-like baths. Often, the bath degrees of freedom may be represented in terms of canonically conjugate coordinates, but in some cases they may require a non-canonical or non-Hamiltonian representation. Herein, we review an approach to the dynamics and statistical mechanics of quantum subsystems embedded in either non-canonical or non-Hamiltonian classical-like baths which is based on operator-valued quasi-probability functions. These functions typically evolve through the action of quasi-Lie brackets and their associated Quantum-Classical Liouville Equations, or through quasi-Lie brackets augmented by dissipative terms. Quasi-Lie brackets possess the unique feature that, while conserving the energy (which the Noether theorem links to time-translation symmetry), they violate the time-translation symmetry of their algebra. This fact can be heuristically understood in terms of the dynamics of the open quantum subsystem. We then describe an example in which a quantum subsystem is embedded in a bath of classical spins, which are described by non-canonical coordinates. In this case, it has been shown that an off-diagonal open-bath geometric phase enters into the propagation of the quantum-classical dynamics. Next, we discuss how non-Hamiltonian dynamics may be employed to generate the constant-temperature evolution of phase space degrees of freedom coupled to the quantum subsystem. Constant-temperature dynamics may be generated by either a classical Langevin stochastic process or a Nosé–Hoover deterministic thermostat. These two approaches are not equivalent but have different advantages and drawbacks. In all cases, the calculation of the operator-valued quasi-probability function allows one to compute time-dependent statistical averages of observables. This may be accomplished in practice using a hybrid Molecular Dynamics/Monte Carlo algorithms, which we outline herein.
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
|
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
|
11
|
Dell’Angelo D, Hanna G. Self-Consistent Filtering Scheme for Efficient Calculations of Observables via the Mixed Quantum-Classical Liouville Approach. J Chem Theory Comput 2016; 12:477-85. [DOI: 10.1021/acs.jctc.5b00991] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- David Dell’Angelo
- 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
|
12
|
Bajo JJ, Granucci G, Persico M. Interplay of radiative and nonradiative transitions in surface hopping with radiation-molecule interactions. J Chem Phys 2015; 140:044113. [PMID: 25669511 DOI: 10.1063/1.4862738] [Citation(s) in RCA: 24] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We implemented a method for the treatment of field induced transitions in trajectory surface hopping simulations, in the framework of the local diabatization scheme, especially suited for on-the-fly dynamics. The method is applied to a simple one-dimensional model with an avoided crossing and compared with quantum wavepacket dynamics. The results show the importance of introducing a proper decoherence correction to surface hopping, in order to obtain meaningful results. Also the energy conservation policy of standard surface hopping must be revised: in fact, the quantum wavepacket energetics is well reproduced if energy absorption/emission is allowed for in the hops determined by radiation-molecule coupling. To our knowledge, this is the first time the issues of decoherence and energy conservation have been analyzed in depth to devise a mixed quantum-classical method for dynamics with molecule-field interactions.
Collapse
Affiliation(s)
- Juan José Bajo
- Departamento de Química-Física I, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Giovanni Granucci
- Università di Pisa, Dipartimento di Chimica e Chimica Industriale, via Risorgimento 35, 56126 Pisa, Italy
| | - Maurizio Persico
- Università di Pisa, Dipartimento di Chimica e Chimica Industriale, via Risorgimento 35, 56126 Pisa, Italy
| |
Collapse
|
13
|
Martinez F, Hanna G. Assessment of approximate solutions of the quantum–classical Liouville equation for dynamics simulations of quantum subsystems embedded in classical environments. MOLECULAR SIMULATION 2014. [DOI: 10.1080/08927022.2014.923573] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
|