1
|
He X, Cheng X, Wu B, Liu J. Nonadiabatic Field with Triangle Window Functions on Quantum Phase Space. J Phys Chem Lett 2024; 15:5452-5466. [PMID: 38747729 PMCID: PMC11129318 DOI: 10.1021/acs.jpclett.4c00793] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/14/2024] [Revised: 04/29/2024] [Accepted: 05/01/2024] [Indexed: 05/24/2024]
Abstract
Recent progress on the constraint coordinate-momentum phase space (CPS) formulation of finite-state quantum systems has revealed that the triangle window function approach is an isomorphic representation of the exact population-population correlation function of the two-state system. We use the triangle window (TW) function and the CPS mapping kernel element to formulate a novel useful representation of discrete electronic degrees of freedom (DOFs). When it is employed with nonadiabatic field (NaF) dynamics, a new variant of the NaF approach (i.e., NaF-TW) is proposed. The NaF-TW expression of the population of any adiabatic state is always positive semidefinite. Extensive benchmark tests of model systems in both the condensed phase and gas phase demonstrate that the NaF-TW approach is able to faithfully capture the dynamical interplay between electronic and nuclear DOFs in a broad region, including where the states remain coupled all the time, as well as where the bifurcation characteristic of nuclear motion is important.
Collapse
Affiliation(s)
- Xin He
- Beijing National Laboratory
for Molecular Sciences, Institute of Theoretical and Computational
Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| | - Xiangsong Cheng
- Beijing National Laboratory
for Molecular Sciences, Institute of Theoretical and Computational
Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| | - Baihua Wu
- Beijing National Laboratory
for Molecular Sciences, Institute of Theoretical and Computational
Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| | - Jian Liu
- Beijing National Laboratory
for Molecular Sciences, Institute of Theoretical and Computational
Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| |
Collapse
|
2
|
Hu D, Ying W, Huo P. Resonance Enhancement of Vibrational Polariton Chemistry Obtained from the Mixed Quantum-Classical Dynamics Simulations. J Phys Chem Lett 2023; 14:11208-11216. [PMID: 38055902 PMCID: PMC10726371 DOI: 10.1021/acs.jpclett.3c02985] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/25/2023] [Revised: 11/21/2023] [Accepted: 12/01/2023] [Indexed: 12/08/2023]
Abstract
We applied a variety of mixed quantum-classical (MQC) approaches to simulate the VSC-influenced reaction rate constant. All of these MQC simulations treat the key vibrational levels associated with the reaction coordinate in the quantum subsystem (as quantum states), whereas all other degrees of freedom (DOFs) are treated inside the classical subsystem. We find that, as long as we have the quantum state descriptions for the vibrational DOFs, one can correctly describe the VSC resonance condition when the cavity frequency matches the bond vibrational frequency. This correct resonance behavior can be obtained regardless of the detailed MQC methods that one uses. The results suggest that the MQC approaches can generate semiquantitative agreement with the exact results for rate constant changes when changing the cavity frequency, the light-matter coupling strength, or the cavity lifetime. The finding of this work suggests that one can use computationally economic MQC approaches to explore the collective coupling scenario when many molecules are collectively coupled to many cavity modes in the future.
Collapse
Affiliation(s)
- Deping Hu
- Center
for Advanced Materials Research, Beijing
Normal University, Zhuhai 519087, China
| | - Wenxiang Ying
- Department
of Chemistry, University of Rochester, Rochester, New York 14627, United States
| | - Pengfei Huo
- Department
of Chemistry, University of Rochester, Rochester, New York 14627, United States
- Institute
of Optics, Hajim School of Engineering, University of Rochester, Rochester, New York 14627, United States
| |
Collapse
|
3
|
Bossion D, Chowdhury SN, Huo P. Non-adiabatic ring polymer molecular dynamics in the phase space of the SU(N) Lie group. J Chem Phys 2023; 158:044123. [PMID: 36725494 DOI: 10.1063/5.0133970] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023] Open
Abstract
We derive the non-adiabatic ring polymer molecular dynamics (RPMD) approach in the phase space of the SU(N) Lie Group. This method, which we refer to as the spin mapping non-adiabatic RPMD (SM-NRPMD), is based on the spin-mapping formalism for the electronic degrees of freedom (DOFs) and ring polymer path-integral description for the nuclear DOFs. Using the Stratonovich-Weyl transform for the electronic DOFs and the Wigner transform for the nuclear DOFs, we derived an exact expression of the Kubo-transformed time-correlation function (TCF). We further derive the spin mapping non-adiabatic Matsubara dynamics using the Matsubara approximation that removes the high frequency nuclear normal modes in the TCF and derive the SM-NRPMD approach from the non-adiabatic Matsubara dynamics by discarding the imaginary part of the Liouvillian. The SM-NRPMD method has numerical advantages compared to the original NRPMD method based on the Meyer-Miller-Stock-Thoss (MMST) mapping formalism due to a more natural mapping using the SU(N) Lie Group that preserves the symmetry of the original system. We numerically compute the Kubo-transformed position auto-correlation function and electronic population correlation function for three-state model systems. The numerical results demonstrate the accuracy of the SM-NRPMD method, which outperforms the original MMST-based NRPMD. We envision that the SM-NRPMD method will be a powerful approach to simulate electronic non-adiabatic dynamics and nuclear quantum effects accurately.
Collapse
Affiliation(s)
- Duncan Bossion
- Department of Chemistry, University of Rochester, 120 Trustee Road, Rochester, New York 14627, USA
| | - Sutirtha N Chowdhury
- Department of Chemistry, University of Rochester, 120 Trustee Road, Rochester, New York 14627, USA
| | - Pengfei Huo
- Department of Chemistry, University of Rochester, 120 Trustee Road, Rochester, New York 14627, USA
| |
Collapse
|
4
|
Saller MAC, Lai Y, Geva E. An Accurate Linearized Semiclassical Approach for Calculating Cavity-Modified Charge Transfer Rate Constants. J Phys Chem Lett 2022; 13:2330-2337. [PMID: 35245071 DOI: 10.1021/acs.jpclett.2c00122] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/14/2023]
Abstract
We show that combining the linearized semiclasscial approximation with Fermi's golden rule (FGR) rate theory gives rise to a general-purpose cost-effective and scalable computational framework that can accurately capture the cavity-induced rate enhancement of charge transfer reactions that occurs when the molecular system is placed inside a microcavity. Both partial linearization with respect to the nuclear and photonic degrees of freedom and full linerization with respect to nuclear, photonic, and electronic degrees of freedom (the latter within the mapping Hamiltonian approach) are shown to be highly accurate, provided that the Wigner transforms of the product (WoP) of operators at the initial time is not replaced by the product of their Wigner transforms. We also show that the partial linearization method yields the quantum-mechanically exact cavity-modified FGR rate constant for a model system in which the donor and acceptor potential energy surfaces are harmonic and identical except for a shift in the equilibrium energy and geometry, if WoP is applied.
Collapse
Affiliation(s)
- Maximilian A C Saller
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| | - Yifan Lai
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| | - Eitan Geva
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| |
Collapse
|
5
|
Dodin A, Provazza J, Coker DF, Willard AP. Trajectory Ensemble Methods Provide Single-Molecule Statistics for Quantum Dynamical Systems. J Chem Theory Comput 2022; 18:2047-2061. [PMID: 35230105 DOI: 10.1021/acs.jctc.1c00477] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The emergence of experiments capable of probing quantum dynamics at the single-molecule level requires the development of new theoretical tools capable of simulating and analyzing these dynamics beyond an ensemble-averaged description. In this article, we present an efficient method for sampling and simulating the dynamics of the individual quantum systems that make up an ensemble and apply it to study the nonequilibrium dynamics of the ubiquitous spin-boson model. We generate an ensemble of single-system trajectories, and we analyze this trajectory ensemble using tools from classical statistical mechanics. Our results demonstrate that the dynamics of quantum coherence is highly heterogeneous at the single-system level due to variations in the initial bath configuration, which significantly affects the transient exchange of coherence between the system and its bath. We observe that single systems tend to retain coherence over time scales longer than that of the ensemble. We also compute a novel thermodynamic entanglement entropy that quantifies a thermodynamic driving force favoring system-bath entanglement.
Collapse
Affiliation(s)
- Amro Dodin
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Justin Provazza
- Department of Chemistry, Northwestern University, Evanston, Illinois 60208, United States.,Department of Chemistry, Boston University, Boston, Massachusetts 02215, United States
| | - David F Coker
- Department of Chemistry, Boston University, Boston, Massachusetts 02215, United States
| | - Adam P Willard
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| |
Collapse
|
6
|
Ananth N. Path Integrals for Nonadiabatic Dynamics: Multistate Ring Polymer Molecular Dynamics. Annu Rev Phys Chem 2022; 73:299-322. [PMID: 35081325 DOI: 10.1146/annurev-physchem-082620-021809] [Citation(s) in RCA: 13] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
This review focuses on a recent class of path-integral-based methods that simulate nonadiabatic dynamics in the condensed phase using only classical molecular dynamics trajectories in an extended phase space. Specifically, a semiclassical mapping protocol is used to derive an exact, continuous, Cartesian variable path-integral representation for the canonical partition function of a system in which multiple electronic states are coupled to nuclear degrees of freedom. Building on this exact statistical foundation, multistate ring polymer molecular dynamics methods are developed for the approximate calculation of real-time thermal correlation functions. The remarkable promise of these multistate ring polymer methods, their successful applications, and their limitations are discussed in detail.Expected final online publication date for the Annual Review of Physical Chemistry, Volume 73 is April 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.
Collapse
Affiliation(s)
- Nandini Ananth
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York, USA;
| |
Collapse
|
7
|
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
|
8
|
He X, Wu B, Gong Z, Liu J. Commutator Matrix in Phase Space Mapping Models for Nonadiabatic Quantum Dynamics. J Phys Chem A 2021; 125:6845-6863. [PMID: 34339600 DOI: 10.1021/acs.jpca.1c04429] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023]
Abstract
We show that a novel, general phase space mapping Hamiltonian for nonadiabatic systems, which is reminiscent of the renowned Meyer-Miller mapping Hamiltonian, involves a commutator variable matrix rather than the conventional zero-point-energy parameter. In the exact mapping formulation on constraint space for phase space approaches for nonadiabatic dynamics, the general mapping Hamiltonian with commutator variables can be employed to generate approximate trajectory-based dynamics. Various benchmark model tests, which range from gas phase to condensed phase systems, suggest that the overall performance of the general mapping Hamiltonian is better than that of the conventional Meyer-Miller Hamiltonian.
Collapse
Affiliation(s)
- Xin He
- Beijing National Laboratory for Molecular Sciences, Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| | - Baihua Wu
- Beijing National Laboratory for Molecular Sciences, Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| | - Zhihao Gong
- Beijing National Laboratory for Molecular Sciences, Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| | - Jian Liu
- Beijing National Laboratory for Molecular Sciences, Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| |
Collapse
|
9
|
Lang H, Vendrell O, Hauke P. Generalized discrete truncated Wigner approximation for nonadiabatic quantum-classical dynamics. J Chem Phys 2021; 155:024111. [PMID: 34266254 DOI: 10.1063/5.0054696] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Nonadiabatic molecular dynamics occur in a wide range of chemical reactions and femtochemistry experiments involving electronically excited states. These dynamics are hard to treat numerically as the system's complexity increases, and it is thus desirable to have accurate yet affordable methods for their simulation. Here, we introduce a linearized semiclassical method, the generalized discrete truncated Wigner approximation (GDTWA), which is well-established in the context of quantum spin lattice systems, into the arena of chemical nonadiabatic systems. In contrast to traditional continuous mapping approaches, e.g., the Meyer-Miller-Stock-Thoss and the spin mappings, GDTWA samples the electron degrees of freedom in a discrete phase space and thus forbids an unphysical unbounded growth of electronic state populations. The discrete sampling also accounts for an effective reduced but non-vanishing zero-point energy without an explicit parameter, which makes it possible to treat the identity operator and other operators on an equal footing. As numerical benchmarks on two linear vibronic coupling models and Tully's models show, GDTWA has a satisfactory accuracy in a wide parameter regime, independent of whether the dynamics is dominated by relaxation or by coherent interactions. Our results suggest that the method can be very adequate to treat challenging nonadiabatic dynamics problems in chemistry and related fields.
Collapse
Affiliation(s)
- Haifeng Lang
- Theoretical Chemistry, Institute of Physical Chemistry, Heidelberg University, Im Neuenheimer Feld 229, 69120 Heidelberg, Germany
| | - Oriol Vendrell
- Theoretical Chemistry, Institute of Physical Chemistry, Heidelberg University, Im Neuenheimer Feld 229, 69120 Heidelberg, Germany
| | - Philipp Hauke
- INO-CNR BEC Center and Department of Physics, University of Trento, Via Sommarive 14, I-38123 Trento, Italy
| |
Collapse
|
10
|
Hu D, Xie Y, Peng J, Lan Z. On-the-Fly Symmetrical Quasi-Classical Dynamics with Meyer-Miller Mapping Hamiltonian for the Treatment of Nonadiabatic Dynamics at Conical Intersections. J Chem Theory Comput 2021; 17:3267-3279. [PMID: 34028268 DOI: 10.1021/acs.jctc.0c01249] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The on-the-fly version of the symmetrical quasi-classical dynamics method based on the Meyer-Miller mapping Hamiltonian (SQC/MM) is implemented to study the nonadiabatic dynamics at conical intersections of polyatomic systems. The current on-the-fly implementation of the SQC/MM method is based on the adiabatic representation and the dressed momentum. To include the zero-point energy (ZPE) correction of the electronic mapping variables, we employ both the γ-adjusted and γ-fixed approaches. Nonadiabatic dynamics of the methaniminium cation (CH2NH2+) and azomethane are simulated using the on-the-fly SQC/MM method. For CH2NH2+, both ZPE correction approaches give reasonable and consistent results. However, for azomethane, the γ-adjusted version of the SQC/MM dynamics behaves much better than the γ-fixed version. Further analysis indicates that it is always recommended to use the γ-adjusted SQC/MM dynamics in the on-the-fly simulation of photoinduced dynamics of polyatomic systems, particularly when the excited state is well separated from the ground state in the Franck-Condon region. This work indicates that the on-the-fly SQC/MM method is a powerful simulation protocol to deal with the nonadiabatic dynamics of realistic polyatomic systems.
Collapse
Affiliation(s)
- Deping Hu
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety & MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou 510006, China.,School of Environment, South China Normal University, Guangzhou 510006, China
| | - Yu Xie
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety & MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou 510006, China.,School of Environment, South China Normal University, Guangzhou 510006, China
| | - Jiawei Peng
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety & MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou 510006, China.,School of Environment, South China Normal University, Guangzhou 510006, China
| | - Zhenggang Lan
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety & MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou 510006, China.,School of Environment, South China Normal University, Guangzhou 510006, China
| |
Collapse
|
11
|
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
|
12
|
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
|
13
|
Mannouch JR, Richardson JO. A partially linearized spin-mapping approach for nonadiabatic dynamics. I. Derivation of the theory. J Chem Phys 2020; 153:194109. [DOI: 10.1063/5.0031168] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
14
|
Mannouch JR, Richardson JO. A partially linearized spin-mapping approach for nonadiabatic dynamics. II. Analysis and comparison with related approaches. J Chem Phys 2020; 153:194110. [DOI: 10.1063/5.0031173] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
15
|
Mandal A, Krauss TD, Huo P. Polariton-Mediated Electron Transfer via Cavity Quantum Electrodynamics. J Phys Chem B 2020; 124:6321-6340. [PMID: 32589846 DOI: 10.1021/acs.jpcb.0c03227] [Citation(s) in RCA: 63] [Impact Index Per Article: 15.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We investigate the polariton-mediated electron transfer reaction in a model system with analytic rate constant theory and direct quantum dynamical simulations. We demonstrate that the photoinduced charge transfer reaction between a bright donor state and dark acceptor state can be significantly enhanced or suppressed by coupling the molecular system to the quantized radiation field inside an optical cavity. This is because the quantum light-matter interaction can influence the effective driving force and electronic couplings between the donor state, which is the hybrid light-matter excitation, and the molecular acceptor state. Under the resonance condition between the photonic and electronic excitations, the effective driving force can be tuned by changing the light-matter coupling strength; for an off-resonant condition, the same effect can be accomplished by changing the molecule-cavity detuning. We further demonstrate that using both the electronic coupling and light-matter coupling helps to extend the effective couplings across the entire system, even for the dark state that carries a zero transition dipole. Theoretically, we find that both the counter-rotating terms and the dipole self-energy in the quantum electrodynamics Hamiltonian are important for obtaining an accurate polariton eigenspectrum as well as the polariton-mediated charge transfer rate constant, especially in the ultrastrong coupling regime.
Collapse
Affiliation(s)
- Arkajit Mandal
- Department of Chemistry, University of Rochester, 120 Trustee Road, Rochester, New York 14627, United States
| | - Todd D Krauss
- Department of Chemistry, University of Rochester, 120 Trustee Road, Rochester, New York 14627, United States
| | - Pengfei Huo
- Department of Chemistry, University of Rochester, 120 Trustee Road, Rochester, New York 14627, United States
| |
Collapse
|
16
|
Heller ER, Richardson JO. Semiclassical instanton formulation of Marcus–Levich–Jortner theory. J Chem Phys 2020; 152:244117. [DOI: 10.1063/5.0013521] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Eric R. Heller
- Laboratory of Physical Chemistry, ETH Zürich, 8093 Zürich, Switzerland
| | | |
Collapse
|
17
|
Runeson JE, Richardson JO. Generalized spin mapping for quantum-classical dynamics. J Chem Phys 2020; 152:084110. [DOI: 10.1063/1.5143412] [Citation(s) in RCA: 27] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023] Open
Affiliation(s)
- Johan E. Runeson
- Laboratory of Physical Chemistry, ETH Zürich, 8093 Zürich, Switzerland
| | | |
Collapse
|
18
|
Abstract
Dynamics at molecule-metal interfaces are a subject of intense current interest and come in many different flavors of experiments: gas-phase scattering, chemisorption, electrochemistry, nanojunction transport, and heterogeneous catalysis, to name a few. These dynamics involve nuclear degrees of freedom entangled with many electronic degrees of freedom (in the metal), and as such there is always the possibility for nonadiabatic phenomena to appear: the nuclei do not necessarily need to move slower than the electrons to break the Born-Oppenheimer (BO) approximation. In this Feature Article, we review a set of dynamical methods developed recently to deal with such nonadiabatic phenomena at a metal surface, methods that serve as alternatives to Tully's independent electron surface hopping (IESH) model. In the weak molecule-metal coupling regime, a classical master equation (CME) can be derived and a simple surface hopping approach is proposed to propagate nuclear and electronic dynamics stochastically. In the strong molecule-metal interaction regime, a Fokker-Planck equation can be derived for the nuclear dynamics, with electronic DoFs incorporated into the overall friction and random force. Lastly, a broadened classical master equation (BCME) can interpolate between the weak and strong molecule-metal interactions. Here, we briefly review these methods and the relevant benchmarking data, showing in particular how the methods can be used to calculate nonequilibrium transport properties. We highlight several open questions and pose several avenues for future study.
Collapse
Affiliation(s)
- Wenjie Dou
- Department of Chemistry , University of California, Berkeley , Berkeley , California 94720 , United States
| | - Joseph E Subotnik
- Department of Chemistry , University of Pennsylvania , Philadelphia , Pennsylvania 19104 , United States
| |
Collapse
|
19
|
Zheng J, Xie Y, Jiang S, Long Y, Ning X, Lan Z. Initial sampling in symmetrical quasiclassical dynamics based on Li-Miller mapping Hamiltonian. Phys Chem Chem Phys 2019; 21:26502-26514. [PMID: 31777888 DOI: 10.1039/c9cp03975a] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Abstract
A symmetrical quasiclassical (SQC) dynamics approach based on the Li-Miller (LM) mapping Hamiltonian (SQC-LM) was employed to describe nonadiabatic dynamics. In principle, the different initial sampling procedures may be applied in the SQC-LM dynamics, and the results may be dependent on different initial sampling. We provided various initial sampling approaches and checked their influence. We selected two groups of models including site-exciton models for exciton dynamics and linear vibronic coupling models for conical intersections to test the performance of SQC-LM dynamics with the different initial sampling methods. The results were examined with respect to those of the accurate multiconfigurational time-dependent Hartree (MCTDH) quantum dynamics. For both the models, the SQC-LM method more-or-less gives a reasonable description of the population dynamics, while the influence of the initial sampling approaches on the final results is noticeable. It seems that the suitable initial sampling methods should be determined by the system under study. This indicates that the combination of the SQC-LM method with a suitable sampling approach may be a potential method in the description of nonadiabatic dynamics.
Collapse
Affiliation(s)
- Jie Zheng
- Industrial Research Institute of Nonwovens & Technical Textiles, College of Textiles Clothing, Qingdao University, Qingdao 266071, China.
| | | | | | | | | | | |
Collapse
|
20
|
Vaillant CL, Thapa MJ, Vaníček J, Richardson JO. Semiclassical analysis of the quantum instanton approximation. J Chem Phys 2019; 151:144111. [DOI: 10.1063/1.5123800] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Christophe L. Vaillant
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, École Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Manish J. Thapa
- Laboratory of Physical Chemistry, ETH Zürich, 8093 Zürich, Switzerland
| | - Jiří Vaníček
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, École Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | | |
Collapse
|
21
|
Mandal A, Huo P. Investigating New Reactivities Enabled by Polariton Photochemistry. J Phys Chem Lett 2019; 10:5519-5529. [PMID: 31475529 DOI: 10.1021/acs.jpclett.9b01599] [Citation(s) in RCA: 74] [Impact Index Per Article: 14.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
We perform quantum dynamics simulations to investigate new chemical reactivities enabled by cavity quantum electrodynamics. The quantum light-matter interactions between the molecule and the quantized radiation mode inside an optical cavity create a set of hybridized electronic-photonic states, so-called polaritons. The polaritonic states adapt the curvatures from both the ground and the excited electronic states, opening up new possibilities to control photochemical reactions by exploiting intrinsic quantum behaviors of light-matter interactions. With quantum dynamics simulations, we demonstrate that the selectivity of a model photoisomerization reaction can be controlled by tuning the photon frequency of the cavity mode or the light-matter coupling strength, providing new ways to manipulate chemical reactions via the light-matter interaction. We further investigate collective quantum effects enabled by coupling the quantized radiation mode to multiple molecules. Our results suggest that in the resonance case, a photon is recycled among molecules to enable multiple excited state reactions, thus effectively functioning as a catalyst. In the nonresonance case, molecules emit and absorb virtual photons to initiate excited state reactions through fundamental quantum electrodynamics processes. These results from quantum dynamics simulations reveal basic principles of polariton photochemistry as well as promising reactivities that take advantage of intrinsic quantum behaviors of photons.
Collapse
Affiliation(s)
- Arkajit Mandal
- Department of Chemistry , University of Rochester , 120 Trustee Road , Rochester , New York 14627 , United States
| | - Pengfei Huo
- Department of Chemistry , University of Rochester , 120 Trustee Road , Rochester , New York 14627 , United States
| |
Collapse
|
22
|
He X, Liu J. A new perspective for nonadiabatic dynamics with phase space mapping models. J Chem Phys 2019; 151:024105. [DOI: 10.1063/1.5108736] [Citation(s) in RCA: 30] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Xin He
- Beijing National Laboratory for Molecular Sciences, Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| | - Jian Liu
- Beijing National Laboratory for Molecular Sciences, Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| |
Collapse
|
23
|
Teh HH, Jin BY, Cheng YC. Frozen-mode small polaron quantum master equation with variational bound for excitation energy transfer in molecular aggregates. J Chem Phys 2019; 150:224110. [DOI: 10.1063/1.5096287] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Hung-Hsuan Teh
- Department of Chemistry, National Taiwan University, Taipei City 10617, Taiwan
| | - Bih-Yaw Jin
- Department of Chemistry, National Taiwan University, Taipei City 10617, Taiwan
| | - Yuan-Chung Cheng
- Department of Chemistry, National Taiwan University, Taipei City 10617, Taiwan
| |
Collapse
|
24
|
Thapa MJ, Fang W, Richardson JO. Nonadiabatic quantum transition-state theory in the golden-rule limit. I. Theory and application to model systems. J Chem Phys 2019; 150:104107. [DOI: 10.1063/1.5081108] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Manish J. Thapa
- Laboratory of Physical Chemistry, ETH Zürich, 8093 Zürich, Switzerland
| | - Wei Fang
- Laboratory of Physical Chemistry, ETH Zürich, 8093 Zürich, Switzerland
| | | |
Collapse
|
25
|
Mandal A, Sandoval C. JS, Shakib FA, Huo P. Quasi-Diabatic Propagation Scheme for Direct Simulation of Proton-Coupled Electron Transfer Reaction. J Phys Chem A 2019; 123:2470-2482. [DOI: 10.1021/acs.jpca.9b00077] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Affiliation(s)
- Arkajit Mandal
- Department of Chemistry, University of Rochester, 120 Trustee Road, Rochester, New York 14627, United States
| | - Juan S. Sandoval C.
- Department of Chemistry, University of Rochester, 120 Trustee Road, Rochester, New York 14627, United States
| | - Farnaz A. Shakib
- Department of Chemistry, University of Rochester, 120 Trustee Road, Rochester, New York 14627, United States
| | - Pengfei Huo
- Department of Chemistry, University of Rochester, 120 Trustee Road, Rochester, New York 14627, United States
| |
Collapse
|
26
|
Xie Y, Zheng J, Lan Z. Performance evaluation of the symmetrical quasi-classical dynamics method based on Meyer-Miller mapping Hamiltonian in the treatment of site-exciton models. J Chem Phys 2018; 149:174105. [DOI: 10.1063/1.5047002] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Yu Xie
- Qingdao Institute of Bioenergy and Bioprocess Technology, Chinese Academy of Sciences, Qingdao 266101, Shandong, China
- The Environmental Research Institute; MOE Key Laboratory of Theoretical Chemistry of Environment, South China Normal University, Guangzhou 510006, China
| | - Jie Zheng
- Industrial Research Institute of Nonwovens and Technical Textiles, College of Textiles and Clothing, Qingdao University, Qingdao 266071, China
| | - Zhenggang Lan
- Qingdao Institute of Bioenergy and Bioprocess Technology, Chinese Academy of Sciences, Qingdao 266101, Shandong, China
- The Environmental Research Institute; MOE Key Laboratory of Theoretical Chemistry of Environment, South China Normal University, Guangzhou 510006, China
| |
Collapse
|
27
|
Sandoval C. JS, Mandal A, Huo P. Symmetric quasi-classical dynamics with quasi-diabatic propagation scheme. J Chem Phys 2018; 149:044115. [DOI: 10.1063/1.5036787] [Citation(s) in RCA: 23] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022] Open
Affiliation(s)
| | - Arkajit Mandal
- Department of Chemistry, University of Rochester, 120 Trustee Road, Rochester, New York 14627, USA
| | - Pengfei Huo
- Department of Chemistry, University of Rochester, 120 Trustee Road, Rochester, New York 14627, USA
| |
Collapse
|
28
|
Mandal A, Shakib FA, Huo P. Investigating photoinduced proton coupled electron transfer reaction using quasi diabatic dynamics propagation. J Chem Phys 2018; 148:244102. [DOI: 10.1063/1.5030634] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022] Open
Affiliation(s)
- Arkajit Mandal
- Department of Chemistry, University of Rochester, 120 Trustee Road, Rochester, New York 14627, USA
| | - Farnaz A. Shakib
- Department of Chemistry, University of Rochester, 120 Trustee Road, Rochester, New York 14627, USA
| | - Pengfei Huo
- Department of Chemistry, University of Rochester, 120 Trustee Road, Rochester, New York 14627, USA
| |
Collapse
|
29
|
Tao X, Shushkov P, Miller TF. Path-integral isomorphic Hamiltonian for including nuclear quantum effects in non-adiabatic dynamics. J Chem Phys 2018; 148:102327. [DOI: 10.1063/1.5005544] [Citation(s) in RCA: 36] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Xuecheng Tao
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Philip Shushkov
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Thomas F. Miller
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| |
Collapse
|
30
|
Mandal A, Yamijala SSRKC, Huo P. Quasi-Diabatic Representation for Nonadiabatic Dynamics Propagation. J Chem Theory Comput 2018; 14:1828-1840. [DOI: 10.1021/acs.jctc.7b01178] [Citation(s) in RCA: 31] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Arkajit Mandal
- Department of Chemistry, University of Rochester, 120 Trustee Road, Rochester, New York 14627, United States
| | - Sharma SRKC Yamijala
- Department of Chemistry, University of Rochester, 120 Trustee Road, Rochester, New York 14627, United States
| | - Pengfei Huo
- Department of Chemistry, University of Rochester, 120 Trustee Road, Rochester, New York 14627, United States
| |
Collapse
|
31
|
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
|
32
|
Provazza J, Segatta F, Garavelli M, Coker DF. Semiclassical Path Integral Calculation of Nonlinear Optical Spectroscopy. J Chem Theory Comput 2018; 14:856-866. [DOI: 10.1021/acs.jctc.7b01063] [Citation(s) in RCA: 27] [Impact Index Per Article: 4.5] [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
- European Center for Theoretical Studies in Nuclear Physics and Related Areas (ECT*-FBK), 38123 Trento, Italy
- Dipartimento
di Chimica Industriale “Toso Montanari”, University of Bologna, Viale del Risorgimento, 4, 40136 Bologna, Italy
| | - Marco Garavelli
- 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
|
33
|
Chowdhury SN, Huo P. Coherent state mapping ring polymer molecular dynamics for non-adiabatic quantum propagations. J Chem Phys 2017; 147:214109. [DOI: 10.1063/1.4995616] [Citation(s) in RCA: 37] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Sutirtha N. Chowdhury
- Department of Chemistry, University of Rochester, 120 Trustee Road, Rochester, New York 14627, USA
| | - Pengfei Huo
- Department of Chemistry, University of Rochester, 120 Trustee Road, Rochester, New York 14627, USA
| |
Collapse
|
34
|
Castellanos MA, Huo P. Enhancing Singlet Fission Dynamics by Suppressing Destructive Interference between Charge-Transfer Pathways. J Phys Chem Lett 2017; 8:2480-2488. [PMID: 28520444 DOI: 10.1021/acs.jpclett.7b00972] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
We apply a real-time path-integral approach to investigate the charge-transfer (CT)-mediated singlet fission quantum dynamics in a model pentacene dimer. Our path-integral method gives reliable fission dynamics across various reaction regimes as well as a broad range of reorganization energies and temperatures. With this method, we investigated the destructive interference between the two CT-mediated fission pathways and discovered two mechanisms that can suppress this deleterious effect. First, increasing the energy gap between the two CT states effectively shuts down the high-lying CT pathway, leaving a better functioning low-lying CT pathway with a minimum amount of destructive interference. Second, intermolecular vibrations induce electronic coupling fluctuations, such that the destructive cancellations due to the different signs in static electronic couplings are suppressed. Our numerical results suggest that these two effects can enhance the fission rate up to three times. These findings reveal promising design principles for more efficient singlet fission materials.
Collapse
Affiliation(s)
- Maria A Castellanos
- Department of Chemistry, University of Rochester , 120 Trustee Road, Rochester, New York 14627, United States
- Department of Chemistry, Universidad Icesi , Cali, Colombia
| | - Pengfei Huo
- Department of Chemistry, University of Rochester , 120 Trustee Road, Rochester, New York 14627, United States
| |
Collapse
|
35
|
Welsch R, Song K, Shi Q, Althorpe SC, Miller TF. Non-equilibrium dynamics from RPMD and CMD. J Chem Phys 2016; 145:204118. [DOI: 10.1063/1.4967958] [Citation(s) in RCA: 42] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023] Open
Affiliation(s)
- Ralph Welsch
- Division of Chemistry and Chemical Engineering, California Institute of Technology,1200 E. California Blvd., Pasadena, California 91125, USA
| | - Kai Song
- Beijing National Laboratory for Molecular Sciences, State Key Laboratory for Structural Chemistry of Unstable and Stable Species, Institute of Chemistry, Chinese Academy of Sciences, Zhongguancun, Beijing 100190, China
| | - Qiang Shi
- Beijing National Laboratory for Molecular Sciences, State Key Laboratory for Structural Chemistry of Unstable and Stable Species, Institute of Chemistry, Chinese Academy of Sciences, Zhongguancun, Beijing 100190, China
| | - Stuart C. Althorpe
- Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Thomas F. Miller
- Division of Chemistry and Chemical Engineering, California Institute of Technology,1200 E. California Blvd., Pasadena, California 91125, USA
| |
Collapse
|
36
|
Ouyang W, Dou W, Jain A, Subotnik JE. Dynamics of Barrier Crossings for the Generalized Anderson–Holstein Model: Beyond Electronic Friction and Conventional Surface Hopping. J Chem Theory Comput 2016; 12:4178-83. [DOI: 10.1021/acs.jctc.6b00533] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Wenjun Ouyang
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| | - Wenjie Dou
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| | - Amber Jain
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| | - Joseph E. Subotnik
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| |
Collapse
|
37
|
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
|
38
|
Subotnik JE, Jain A, Landry B, Petit A, Ouyang W, Bellonzi N. Understanding the Surface Hopping View of Electronic Transitions and Decoherence. Annu Rev Phys Chem 2016; 67:387-417. [DOI: 10.1146/annurev-physchem-040215-112245] [Citation(s) in RCA: 234] [Impact Index Per Article: 29.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Affiliation(s)
- Joseph E. Subotnik
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104;
| | - Amber Jain
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104;
| | - Brian Landry
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104;
| | - Andrew Petit
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104;
| | - Wenjun Ouyang
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104;
| | - Nicole Bellonzi
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104;
| |
Collapse
|
39
|
Antipov SV, Ye Z, Ananth N. Dynamically consistent method for mixed quantum-classical simulations: A semiclassical approach. J Chem Phys 2016; 142:184102. [PMID: 25978878 DOI: 10.1063/1.4919667] [Citation(s) in RCA: 28] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We introduce a new semiclassical (SC) framework, the Mixed Quantum-Classical Initial Value Representation (MQC-IVR), that can be tuned to reproduce existing quantum-limit and classical-limit SC approximations to quantum real-time correlation functions. Applying a modified Filinov transformation to a quantum-limit SC formulation leads to the association of a Filinov parameter with each degree of freedom in the system; varying this parameter from zero to infinity controls the extent of quantization of the corresponding mode. The resulting MQC-IVR expression provides a consistent dynamic framework for mixed quantum-classical simulations and we demonstrate its numerical accuracy in the calculation of real-time correlation functions for a model 1D system and a model 2D system over the full range of quantum- to classical-limit behaviors.
Collapse
Affiliation(s)
- Sergey V Antipov
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, USA
| | - Ziyu Ye
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, USA
| | - Nandini Ananth
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, USA
| |
Collapse
|
40
|
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
|
41
|
Kretchmer JS, Miller III TF. Kinetically-constrained ring-polymer molecular dynamics for non-adiabatic chemistries involving solvent and donor–acceptor dynamical effects. Faraday Discuss 2016; 195:191-214. [DOI: 10.1039/c6fd00143b] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We investigate the performance of the recently developed kinetically-constrained ring polymer molecular dynamics (KC-RPMD) method for the description of model condensed-phase electron transfer (ET) reactions in which solvent and donor–acceptor dynamics play an important role. Comparison of KC-RPMD with results from Golden-Rule rate theories and numerically exact quantum dynamics calculations demonstrates that KC-RPMD accurately captures the combination of electronic- and nuclear-dynamical effects throughout the Marcus (intermediate solvent friction) and Zusman (large solvent friction) regimes of ET. It is also demonstrated that KC-RPMD accurately describes systems in which the magnitude of the diabatic coupling depends on the position of a dynamical donor–acceptor mode. In addition to these successes, however, we present an unsurprising failure of KC-RPMD to capture the enhancement of the ET rate in the low solvent friction regime associated with nuclear coherence effects. In this analysis, we re-visit several aspects of the original KC-RPMD formulation, including the form of the kinetic constraint and the choice of the mass of the auxiliary electronic variable. In particular, we introduce a Langevin bath for the auxiliary electronic variable to correct for its unphysically low coupling with the nuclear degrees of freedom. This work demonstrates that the KC-RPMD method is well suited for the direct simulation of non-adiabatic donor–acceptor chemistries associated with many complex systems, including those for which solvent dynamics plays an important role in the reaction mechanism.
Collapse
Affiliation(s)
- Joshua S. Kretchmer
- Division of Chemistry and Chemical Engineering
- California Institute of Technology
- Pasadena
- USA
| | - Thomas F. Miller III
- Division of Chemistry and Chemical Engineering
- California Institute of Technology
- Pasadena
- USA
| |
Collapse
|
42
|
Hele TJH, Ananth N. Deriving the exact nonadiabatic quantum propagator in the mapping variable representation. Faraday Discuss 2016; 195:269-289. [DOI: 10.1039/c6fd00106h] [Citation(s) in RCA: 26] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Abstract
We derive an exact quantum propagator for nonadiabatic dynamics in multi-state systems using the mapping variable representation, where classical-like Cartesian variables are used to represent both continuous nuclear degrees of freedom and discrete electronic states. The resulting Liouvillian is a Moyal series that, when suitably approximated, can allow for the use of classical dynamics to efficiently model large systems. We demonstrate that different truncations of the exact Liouvillian lead to existing approximate semiclassical and mixed quantum–classical methods and we derive an associated error term for each method. Furthermore, by combining the imaginary-time path-integral representation of the Boltzmann operator with the exact Liouvillian, we obtain an analytic expression for thermal quantum real-time correlation functions. These results provide a rigorous theoretical foundation for the development of accurate and efficient classical-like dynamics to compute observables such as electron transfer reaction rates in complex quantized systems.
Collapse
Affiliation(s)
| | - Nandini Ananth
- Department of Chemistry and Chemical Biology
- Cornell University
- Ithaca
- USA
| |
Collapse
|
43
|
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
|
44
|
Jain A, Subotnik JE. Does Nonadiabatic Transition State Theory Make Sense Without Decoherence? J Phys Chem Lett 2015; 6:4809-4814. [PMID: 26631360 DOI: 10.1021/acs.jpclett.5b02148] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/05/2023]
Abstract
We analyze thermal rate constants as computed with surface hopping dynamics and resolve certain inconsistencies that have permeated the literature. On one hand, according to Landry and Subotnik (J. Chem. Phys. 2012, 137, 22A513), without decoherence, direct dynamics with surface hopping overestimates the rate of relaxation for the spin-boson Hamiltonian. On the other hand, according to Jain and Subotnik (J. Chem. Phys. 2015, 143, 134107), without decoherence, a transition state theory with surface hopping underestimates spin-boson rate constants. In this Letter, we resolve this apparent contradiction. We show that, without decoherence, direct dynamics and transition state theory should not agree; agreement is guaranteed only with decoherence. We also show that, even though the effects of decoherence may be hidden for isoenergetic reactions, these decoherence failures are exposed for exothermic reactions. We believe these lessons are essential when interpreting surface hopping papers published in the literature without any decoherence corrections.
Collapse
Affiliation(s)
- Amber Jain
- Department of Chemistry, University of Pennsylvania , 231 South 34th Street, Philadelphia, Pennsylvania 19104, United States
| | - Joseph E Subotnik
- Department of Chemistry, University of Pennsylvania , 231 South 34th Street, Philadelphia, Pennsylvania 19104, United States
| |
Collapse
|
45
|
Sun X, Geva E. Equilibrium Fermi’s Golden Rule Charge Transfer Rate Constants in the Condensed Phase: The Linearized Semiclassical Method vs Classical Marcus Theory. J Phys Chem A 2015; 120:2976-90. [DOI: 10.1021/acs.jpca.5b08280] [Citation(s) in RCA: 38] [Impact Index Per Article: 4.2] [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
|
46
|
Richardson JO. Ring-polymer instanton theory of electron transfer in the nonadiabatic limit. J Chem Phys 2015; 143:134116. [DOI: 10.1063/1.4932362] [Citation(s) in RCA: 34] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Jeremy O. Richardson
- Institut für Theoretische Physik und Interdisziplinäres Zentrum für Molekulare Materialien, Friedrich-Alexander-Universität Erlangen-Nürnberg (FAU), Staudtstraße 7/B2, 91058 Erlangen, Germany
| |
Collapse
|
47
|
Jain A, Subotnik JE. Surface hopping, transition state theory, and decoherence. II. Thermal rate constants and detailed balance. J Chem Phys 2015; 143:134107. [DOI: 10.1063/1.4930549] [Citation(s) in RCA: 50] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Amber Jain
- Department of Chemistry, University of Pennsylvania, 231 South 34th Street, Philadelphia, Pennsylvania 19104, USA
| | - Joseph E. Subotnik
- Department of Chemistry, University of Pennsylvania, 231 South 34th Street, Philadelphia, Pennsylvania 19104, USA
| |
Collapse
|
48
|
Kondorskiy AD, Nanbu S. Electronically nonadiabatic wave packet propagation using frozen Gaussian scattering. J Chem Phys 2015; 143:114103. [DOI: 10.1063/1.4930923] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Alexey D. Kondorskiy
- P. N. Lebedev Physical Institute of the Russian Academy of Sciences, 53, Leninsky Prospekt, Moscow 119991, Russia
| | - Shinkoh Nanbu
- Department of Materials and Life Sciences, Faculty of Science and Technology, Sophia University, 7-1 Kioi-cho, Chiyoda-ku, Tokyo 102-8554, Japan
| |
Collapse
|
49
|
Hele TJH, Willatt MJ, Muolo A, Althorpe SC. Boltzmann-conserving classical dynamics in quantum time-correlation functions: “Matsubara dynamics”. J Chem Phys 2015; 142:134103. [DOI: 10.1063/1.4916311] [Citation(s) in RCA: 71] [Impact Index Per Article: 7.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Timothy J. H. Hele
- Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Michael J. Willatt
- Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Andrea Muolo
- Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Stuart C. Althorpe
- Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| |
Collapse
|
50
|
Kelly A, Brackbill N, Markland TE. Accurate nonadiabatic quantum dynamics on the cheap: Making the most of mean field theory with master equations. J Chem Phys 2015; 142:094110. [DOI: 10.1063/1.4913686] [Citation(s) in RCA: 44] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023] Open
Affiliation(s)
- Aaron Kelly
- Department of Chemistry, Stanford University, Stanford, California 94305, USA
| | - Nora Brackbill
- Department of Physics, Stanford University, Stanford, California 94305, USA
| | - Thomas E. Markland
- Department of Chemistry, Stanford University, Stanford, California 94305, USA
| |
Collapse
|