-
Semiclassical Dynamics in Wigner Phase Space II : nonadiabatic Hybrid Wigner Dynamics
Authors:
Shreyas Malpathak,
Nandini Ananth
Abstract:
We present an approximate semiclassical (SC) framework for mixed quantized dynamics in Wigner phase space in a two-part series. In the first article, we introduced the Adiabatic Hybrid Wigner Dynamics (AHWD) method that allows for a few important `system' degrees of freedom to be quantized using high-level Double Herman-Kluk SC theory while describing the rest (the `bath') using classical-limit Li…
▽ More
We present an approximate semiclassical (SC) framework for mixed quantized dynamics in Wigner phase space in a two-part series. In the first article, we introduced the Adiabatic Hybrid Wigner Dynamics (AHWD) method that allows for a few important `system' degrees of freedom to be quantized using high-level Double Herman-Kluk SC theory while describing the rest (the `bath') using classical-limit Linearized SC theory. In this second article, we extend our hybrid Wigner dynamics to nonadiabatic processes. The resulting Nonadiabatic Hybrid Wigner Dynamics (NHWD) has two variants that differ in the choice of degrees of freedom to be quantized. Specifically, we introduce NHWD(E) where only the electronic state variables are quantized and the NHWD(V) where both electronic state variables and a handful of strongly coupled nuclear modes are quantized. We show that while NHWD(E) proves accurate for a wide range of scattering models and spin-boson models, systems where a few nuclear modes are strongly coupled to electronic states require NHWD(V) to accurately capture the long-time dynamics. Taken together, we show that AHWD and NHWD represent a new framework for SC simulations of high-dimensional systems with significant quantum effects.
△ Less
Submitted 12 June, 2024; v1 submitted 30 May, 2024;
originally announced May 2024.
-
Semiclassical Dynamics in Wigner Phase Space I : Adiabatic Hybrid Wigner Dynamics
Authors:
Shreyas Malpathak,
Nandini Ananth
Abstract:
The Wigner phase space formulation of quantum mechanics is a complete framework for quantum dynamic calculations that elegantly highlights connections with classical dynamics. In this series of two articles, building upon previous efforts, we derive the full hierarchy of approximate semiclassical (SC) dynamic methods for adiabatic and non-adiabatic problems in Wigner phase space. In paper I, focus…
▽ More
The Wigner phase space formulation of quantum mechanics is a complete framework for quantum dynamic calculations that elegantly highlights connections with classical dynamics. In this series of two articles, building upon previous efforts, we derive the full hierarchy of approximate semiclassical (SC) dynamic methods for adiabatic and non-adiabatic problems in Wigner phase space. In paper I, focusing on adiabatic single surface processes, we derive the well-known Double Herman-Kluk (DHK) approximation for real-time correlation functions in Wigner phase space, and connect it to the Linearized SC (LSC) approximation through a stationary phase approximation. We exploit this relationship to introduce a new hybrid SC method, termed Adiabatic Hybrid Wigner Dynamics (AHWD) that allows for a few important `system' degrees of freedom (dofs) to be treated at the DHK level, while treating the rest of the dofs (the `bath') at the LSC level. AHWD is shown to accurately capture quantum interference effects in models of coupled oscillators, and the decoherence of vibrational probability density of a model I$_2$ Morse oscillator coupled to an Ohmic thermal bath. We show that AHWD significantly mitigates the sign-problem and employs a reduced dimensional prefactor bringing calculations of complex system-bath problems within the reach of SC methods. Paper II focuses on extending this hybrid SC dynamics to nonadiabatic processes.
△ Less
Submitted 12 June, 2024; v1 submitted 30 May, 2024;
originally announced May 2024.
-
Hybrid Quantum Algorithm for Simulating Real-Time Thermal Correlation Functions
Authors:
Elliot C. Eklund,
Nandini Ananth
Abstract:
We present a hybrid Path Integral Monte Carlo (hPIMC) algorithm to calculate real-time quantum thermal correlation functions and demonstrate its application to open quantum systems. The hPIMC algorithm leverages the successes of classical PIMC as a computational tool for high-dimensional system studies by exactly simulating dissipation using the Feynman-Vernon influence functional on a classical c…
▽ More
We present a hybrid Path Integral Monte Carlo (hPIMC) algorithm to calculate real-time quantum thermal correlation functions and demonstrate its application to open quantum systems. The hPIMC algorithm leverages the successes of classical PIMC as a computational tool for high-dimensional system studies by exactly simulating dissipation using the Feynman-Vernon influence functional on a classical computer. We achieve a quantum speed-up over the classical algorithm by computing short-time matrix elements of the quantum propagator on a quantum computer. We show that the component of imaginary-time evolution can be performed accurately using the recently developed Probabilistic Imaginary-Time Evolution (PITE) algorithm, and we introduce a novel low-depth circuit for approximate real-time evolution under the kinetic energy operator using a Discrete Variable Representation (DVR). We test the accuracy of the approximation by computing the position-position thermal correlation function of a proton transfer reaction.
△ Less
Submitted 29 May, 2024;
originally announced May 2024.
-
Mean-Field Ring Polymer Rates Using a Population Dividing Surface
Authors:
Nathan London,
Siyu Bu,
Britta Ann Johnson,
Nandini Ananth
Abstract:
Mean-field Ring Polymer Molecular Dynamics (MF-RPMD) offers a computationally efficient method for the simulation of reaction rates in multi-level systems. Previous work has established that, to model a nonadiabatic state-to-state reaction accurately, the dividing surface must be chosen to explicitly sample kinked ring polymer configurations where at least one bead is in a different electronic sta…
▽ More
Mean-field Ring Polymer Molecular Dynamics (MF-RPMD) offers a computationally efficient method for the simulation of reaction rates in multi-level systems. Previous work has established that, to model a nonadiabatic state-to-state reaction accurately, the dividing surface must be chosen to explicitly sample kinked ring polymer configurations where at least one bead is in a different electronic state than the others. Building on this, we introduce a population difference coordinate and a kink-constrained dividing surface, and we test the accuracy of the resulting mean-field rate theory on a series of linear vibronic coupling model systems as well as spin-boson models. We demonstrate that this new MF-RPMD rate approach is efficient to implement and quantitatively accurate for models over a wide range of driving forces, coupling strengths, and temperatures.
△ Less
Submitted 7 May, 2024;
originally announced May 2024.
-
Nonadiabatic simulations of photoisomerization and dissociation in ethylene using ab initio classical trajectories
Authors:
Ken Miyazaki,
Nandini Ananth
Abstract:
We simulate the nonadiabatic dynamics of photo-induced isomerization and dissociation in ethylene using ab initio classical trajectories in an extended phase space of nuclear and electronic variables. This is achieved by employing the Linearized Semiclassical Initial Value Representation (LSC-IVR) method for nonadiabatic dynamics where discrete electronic states are mapped to continuous classical…
▽ More
We simulate the nonadiabatic dynamics of photo-induced isomerization and dissociation in ethylene using ab initio classical trajectories in an extended phase space of nuclear and electronic variables. This is achieved by employing the Linearized Semiclassical Initial Value Representation (LSC-IVR) method for nonadiabatic dynamics where discrete electronic states are mapped to continuous classical variables using either the Meyer-Miller Stock-Thoss representation or a more recently introduced spin mapping approach. Trajectory initial conditions are sampled by constraining electronic state variables to a single initial excited state, and by drawing nuclear phase space configurations from a Wigner distribution at finite temperature. An ensemble of classical ab initio trajectories are then generated to compute thermal population correlation functions and to analyze the mechanisms of isomerization and dissociation. Our results serve as a demonstration that this parameter-free semiclassical approach is computationally efficient and accurate, identifying mechanistic pathways in agreement with previous theoretical studies, and also uncovering dissociation pathways observed experimentally.
△ Less
Submitted 2 August, 2023; v1 submitted 16 June, 2023;
originally announced June 2023.
-
Non-Linear Correlation Functions and Zero-Point Energy Flow in Mixed Quantum-Classical Semiclassical Dynamics
Authors:
Shreyas Malpathak,
Nandini Ananth
Abstract:
Mixed Quantum Classical (MQC)-IVR is a recently introduced semiclassical framework that allows for selective quantization of the modes of a complex system. In the quantum limit, MQC reproduces the semiclassical Double Herman-Kluk IVR results, accurately capturing nuclear quantum coherences and conserving zero-point energy. However, in the classical limit, while MQC mimics the Husimi-IVR for real-t…
▽ More
Mixed Quantum Classical (MQC)-IVR is a recently introduced semiclassical framework that allows for selective quantization of the modes of a complex system. In the quantum limit, MQC reproduces the semiclassical Double Herman-Kluk IVR results, accurately capturing nuclear quantum coherences and conserving zero-point energy. However, in the classical limit, while MQC mimics the Husimi-IVR for real-time correlation functions with linear operators, it is significantly less accurate for non-linear correlation functions with errors even at time zero. Here, we identify the origin of this discrepancy in the MQC formulation and propose a modification. We analytically show that the modified MQC approach is exact for all correlation functions at time zero, and in a study of zero-point energy (ZPE) flow, we numerically demonstrate that it correctly obtains the quantum and classical limits as a function of time. Interestingly, while classical-limit MQC simulations show the expected, unphysical ZPE leakage, we find it is possible to predict and even modify the direction of ZPE flow through selective quantization of the system, with the quantum-limit modes accepting energy additions but preserving the minimum quantum mechanically required energy.
△ Less
Submitted 4 November, 2022;
originally announced November 2022.
-
A Semiclassical Framework for Mixed Quantum Classical Dynamics
Authors:
Shreyas Malpathak,
Matthew S. Church,
Nandini Ananth
Abstract:
Semiclassical approximations for quantum dynamic simulations in complex chemical systems range from rigorously accurate methods that are computationally expensive to methods that exhibit near-classical scaling with system size but are limited in their ability to describe quantum effects. In practical studies of high-dimensional reactions, neither extreme is the best choice: frequently a high-level…
▽ More
Semiclassical approximations for quantum dynamic simulations in complex chemical systems range from rigorously accurate methods that are computationally expensive to methods that exhibit near-classical scaling with system size but are limited in their ability to describe quantum effects. In practical studies of high-dimensional reactions, neither extreme is the best choice: frequently a high-level quantum mechanical description is only required for a handful of modes, while the majority of environment modes that do not play a key role in the reactive event of interest are well served with a lower level of theory. In this feature we introduce Modified Filinov filtration as a powerful tool for mixed quantum-classical simulations in a uniform semiclassical framework.
△ Less
Submitted 21 July, 2022;
originally announced July 2022.
-
A Skew Dividing Surface for Accurate Nonadiabatic Mean-Field Ring Polymer Rates
Authors:
Britta Ann Johnson,
Nandini Ananth
Abstract:
Mean-Field Ring Polymer Molecular Dynamics (MF-RPMD) is a powerful, efficient, and accurate method for approximate quantum dynamic simulations of multi-level system dynamics. Initial efforts to compute nonadiabatic reaction rates using MF-RPMD were not successful; recent work showed that this can be remedied by including a simple, if adhoc, correction term that accounts for the formation of `kinke…
▽ More
Mean-Field Ring Polymer Molecular Dynamics (MF-RPMD) is a powerful, efficient, and accurate method for approximate quantum dynamic simulations of multi-level system dynamics. Initial efforts to compute nonadiabatic reaction rates using MF-RPMD were not successful; recent work showed that this can be remedied by including a simple, if adhoc, correction term that accounts for the formation of `kinked' or mixed electronic state ring polymer configurations. Here, we build on this idea, introducing a electronic state population based reaction coordinate and novel skew dividing surface that constrains nuclear positions to configurations where the reactant and product state potentials are near-degenerate and that samples kinked electronic state configurations. We then demonstrate the numerical accuracy of this method in computing rates for a series of nonadiabatic model systems.
△ Less
Submitted 28 April, 2021; v1 submitted 26 February, 2021;
originally announced March 2021.
-
Multistate ring polymer instantons and nonadiabatic reaction rates
Authors:
Srinath Ranya,
Nandini Ananth
Abstract:
We present two multistate ring polymer instanton (RPI) formulations, both obtained from an exact path integral representation of the quantum canonical partition function for multistate systems. The two RPIs differ in their treatment of the electronic degrees of freedom; whereas the Mean-Field (MF)-RPI averages over the electronic state contributions, the Mapping Variable (MV)-RPI employs explicit…
▽ More
We present two multistate ring polymer instanton (RPI) formulations, both obtained from an exact path integral representation of the quantum canonical partition function for multistate systems. The two RPIs differ in their treatment of the electronic degrees of freedom; whereas the Mean-Field (MF)-RPI averages over the electronic state contributions, the Mapping Variable (MV)-RPI employs explicit continuous Cartesian variables to represent the electronic states.
We compute both RPIs for a series of model two-state systems coupled to a single nuclear mode with electronic coupling values chosen to describe dynamics in both adiabatic and nonadiabatic regimes. We show that the MF-RPI for symmetric systems are in good agreement with previous literature, and we show that our numerical techniques are robust for systems with non-zero driving force. The nuclear MF-RPI and the nuclear MV-RPI are similar, but the MV-RPI uniquely reports on the changes in the electronic state populations along the instanton path. In both cases, we analytically demonstrate the existence of a zero-mode and we numerically and that these solutions are true instantons with a single unstable mode as expected for a first order saddle point. Finally, we use the MF-RPI to accurately calculate rate constants for adiabatic and nonadiabatic model systems with the coupling strength varying over three orders of magnitude.
△ Less
Submitted 26 February, 2020; v1 submitted 10 August, 2019;
originally announced August 2019.
-
Semiclassical dynamics in the mixed quantum-classical limit
Authors:
Matthew S. Church,
Nandini Ananth
Abstract:
The semiclassical Double Herman-Kluk Initial Value Representation is an accurate approach to computing quantum real time correlation functions, but its applications are limited by the need to evaluate an oscillatory integral. In previous work, we have shown that this `sign problem' can be mitigated using the modified Filinov filtration technique to control the extent to which individual modes of t…
▽ More
The semiclassical Double Herman-Kluk Initial Value Representation is an accurate approach to computing quantum real time correlation functions, but its applications are limited by the need to evaluate an oscillatory integral. In previous work, we have shown that this `sign problem' can be mitigated using the modified Filinov filtration technique to control the extent to which individual modes of the system contribute to the overall phase of the integrand. Here we follow this idea to a logical conclusion: we analytically derive a general expression for the mixed quantum-classical limit of the semiclassical correlation function - AMQC-IVR, where the phase contributions from the `classical' modes of the system are filtered while the `quantum' modes are treated in the full semiclassical limit. We numerically demonstrate the accuracy and efficiency of the AMQC-IVR formulation in calculations of quantum correlation functions and reaction rates using three model systems with varied coupling strengths between the classical and quantum subsystems. We also introduce a separable prefactor approximation that further reduces the computational cost, but is only accurate in the limit of weak coupling between the quantum and classical subsystems.
△ Less
Submitted 29 June, 2019;
originally announced July 2019.
-
Anticipating acene-based chromophore spectra with molecular orbital arguments
Authors:
Timothy J. H. Hele,
Eric G. Fuemmeler,
Samuel N. Sanders,
Elango Kumarasamy,
Matthew Y. Sfeir,
Luis M. Campos,
Nandini Ananth
Abstract:
Recent synthetic studies on the organic molecules tetracene and pentacene have found certain dimers and oligomers to exhibit an intense absorption in the visible region of the spectrum which is not present in the monomer or many previously-studied dimers. In this article we combine experimental synthesis with electronic structure theory and spectral computation to show that this absorption arises…
▽ More
Recent synthetic studies on the organic molecules tetracene and pentacene have found certain dimers and oligomers to exhibit an intense absorption in the visible region of the spectrum which is not present in the monomer or many previously-studied dimers. In this article we combine experimental synthesis with electronic structure theory and spectral computation to show that this absorption arises from an otherwise dark charge-transfer excitation 'borrowing intensity' from an intense UV excitation. Further, by characterizing the role of relevant monomer molecular orbitals, we arrive at a design principle that allows us to predict the presence or absence of an additional absorption based on the bonding geometry of the dimer. We find this rule correctly explains the spectra of a wide range of acene derivatives and solves an unexplained structure-spectrum phenomenon first observed seventy years ago. These results pave the way for the design of highly absorbent chromophores with applications ranging from photovoltaics to liquid crystals.
△ Less
Submitted 18 December, 2018;
originally announced December 2018.
-
Nonadiabatic semiclassical dynamics in the mixed quantum-classical initial value representation
Authors:
Matthew S. Church,
Timothy J. H. Hele,
Gregory S. Ezra,
Nandini Ananth
Abstract:
We extend the Mixed Quantum-Classical Initial Value Representation (MQC-IVR), a semiclassical method for computing real-time correlation functions, to electronically nonadiabatic systems using the Meyer-Miller-Stock-Thoss (MMST) Hamiltonian to treat electronic and nuclear degrees of freedom (dofs) within a consistent dynamic framework. We introduce an efficient symplectic integration scheme, the M…
▽ More
We extend the Mixed Quantum-Classical Initial Value Representation (MQC-IVR), a semiclassical method for computing real-time correlation functions, to electronically nonadiabatic systems using the Meyer-Miller-Stock-Thoss (MMST) Hamiltonian to treat electronic and nuclear degrees of freedom (dofs) within a consistent dynamic framework. We introduce an efficient symplectic integration scheme, the MInt algorithm, for numerical time-evolution of the nuclear and electronic phase space variables as well as the Monodromy matrix, under the non-separable MMST Hamiltonian. We then calculate the probability of transmission through a curve-crossing in model two-level systems and show that in the quantum limit MQC-IVR is in good agreement with the exact quantum results, whereas in the classical limit the method yields results in keeping with mean-field approaches like the Linearized Semiclassical IVR. Finally, exploiting the ability of MQC-IVR to quantize different dofs to different extents, we present a detailed study of the extents to which quantizing the nuclear and electronic dofs improves numerical convergence properties without significant loss of accuracy.
△ Less
Submitted 21 September, 2017;
originally announced September 2017.
-
Validating and Implementing Modified Filinov Phase Filtration in Semiclassical Dynamics
Authors:
Matthew S. Church,
Sergey V. Antipov,
Nandini Ananth
Abstract:
The Mixed Quantum-Classical Initial Value Representation (MQC-IVR) is a recently introduced approximate semiclassical (SC) method for the calculation of real-time quantum correlation functions. MQC-IVR employs a modified Filinov filtration (MFF) scheme to control the overall phase of the SC integrand, extending the applicability of SC methods to complex systems while retaining their ability to acc…
▽ More
The Mixed Quantum-Classical Initial Value Representation (MQC-IVR) is a recently introduced approximate semiclassical (SC) method for the calculation of real-time quantum correlation functions. MQC-IVR employs a modified Filinov filtration (MFF) scheme to control the overall phase of the SC integrand, extending the applicability of SC methods to complex systems while retaining their ability to accurately describe quantum coherence effects. Here, we address questions regarding the effectiveness of the MFF scheme in combination with SC dynamics. Previous work showed that this filtering scheme is of limited utility in the context of semiclassical wavepacket propagation, but we find the MFF is extraordinarily powerful in the context of correlation functions. By examining trajectory phase and amplitude contributions to the real-time SC correlation function in a model system, we clearly demonstrate that the MFF serves to reduce noise by damping amplitude only in regions of highly oscillatory phase leading to a reduction in computational effort while retaining accuracy. Further, we introduce a novel and efficient MQC-IVR formulation that allows for linear scaling in computational cost with the total simulation length, a significant improvement over the more-than quadratic scaling exhibited by the original method.
△ Less
Submitted 18 May, 2017;
originally announced May 2017.
-
An Almost Analytical Approach to Simulating 2D Electronic Spectra
Authors:
Pallavi Bhattacharyya,
Nandini Ananth
Abstract:
We introduce an almost analytical method to simulate 2D electronic spectra as a double Fourier transform of the the non-linear response function (NRF) corresponding to a particular optical pulse sequence. We employ a unitary transformation to represent the total system Hamiltonian in a stationary basis that allows us to separate contributions from decoherence and phonon-mediated population relaxat…
▽ More
We introduce an almost analytical method to simulate 2D electronic spectra as a double Fourier transform of the the non-linear response function (NRF) corresponding to a particular optical pulse sequence. We employ a unitary transformation to represent the total system Hamiltonian in a stationary basis that allows us to separate contributions from decoherence and phonon-mediated population relaxation to the NRF. Previously, one of us demonstrated the use of an analytic, cumulant expansion approach to calculate the decoherence term. Here, we extend this idea to obtain an accurate expression for the population relaxation term, a significant improvement over standard quantum master equation-based approximations. We numerically demonstrate the accuracy of our method by computing the photon echo spectrum of a two-level system coupled to a thermal bath, and we highlight the mechanistic insights obtained from our simulation.
△ Less
Submitted 16 May, 2017; v1 submitted 1 May, 2017;
originally announced May 2017.
-
State Space Path Integrals for Electronically Nonadiabatic Reaction Rates
Authors:
Jessica Ryan Duke,
Nandini Ananth
Abstract:
We present a state-space-based path integral method to calculate the rate of electron transfer (ET) in multi-state, multi-electron condensed-phase processes. We employ an exact path integral in discrete electronic states and continuous Cartesian nuclear variables to obtain a transition state theory (TST) estimate to the rate. A dynamic recrossing correction to the TST rate is then obtained from re…
▽ More
We present a state-space-based path integral method to calculate the rate of electron transfer (ET) in multi-state, multi-electron condensed-phase processes. We employ an exact path integral in discrete electronic states and continuous Cartesian nuclear variables to obtain a transition state theory (TST) estimate to the rate. A dynamic recrossing correction to the TST rate is then obtained from real-time dynamics simulations using mean field ring polymer molecular dynamics. We employ two different reaction coordinates in our simulations and show that, despite the use of mean field dynamics, the use of an accurate dividing surface to compute TST rates allows us to achieve remarkable agreement with Fermi's golden rule rates for nonadiabatic ET in the normal regime of Marcus theory. Further, we show that using a reaction coordinate based on electronic state populations allows us to capture the turnover in rates for ET in the Marcus inverted regime.
△ Less
Submitted 4 June, 2016;
originally announced June 2016.
-
Deriving the exact nonadiabatic quantum propagator in the mapping variable representation
Authors:
Timothy J. H. Hele,
Nandini Ananth
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 expression is a Moyal series that, when suitably approximated, can allow for the use of classical dynamics to efficien…
▽ More
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 expression 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 propagator 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 propagator, 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.
△ Less
Submitted 2 May, 2016;
originally announced May 2016.
-
Mapping Variable Ring Polymer Molecular Dynamics: A Path-Integral Based Method for Nonadiabatic Processes
Authors:
Nandini Ananth
Abstract:
We introduce mapping-variable ring polymer molecular dynamics (MV-RPMD), a model dynamics for the direct simulation of multi-electron processes. An extension of the RPMD idea, this method is based on an exact, imaginary time path-integral representation of the quantum Boltzmann operator using continuous Cartesian variables for both electronic states and nuclear degrees of freedom. We demonstrate t…
▽ More
We introduce mapping-variable ring polymer molecular dynamics (MV-RPMD), a model dynamics for the direct simulation of multi-electron processes. An extension of the RPMD idea, this method is based on an exact, imaginary time path-integral representation of the quantum Boltzmann operator using continuous Cartesian variables for both electronic states and nuclear degrees of freedom. We demonstrate the accuracy of the MV-RPMD approach in calculations of real-time, thermal correlation functions for a range of two-state single-mode model systems with different coupling strengths and asymmetries. Further, we show that the ensemble of classical trajectories employed in these simulations preserves the Boltzmann distribution and provides a direct probe into real-time coupling between electronic state transitions and nuclear dynamics.
△ Less
Submitted 17 September, 2013;
originally announced September 2013.
-
Flux-correlation approach to characterizing reaction pathways in quantum systems: A study of condensed-phase proton-coupled electron transfer
Authors:
Nandini Ananth,
Thomas F. Miller III
Abstract:
We introduce a simple method for characterizing reactive pathways in quantum systems. Flux auto- correlation and cross-correlation functions are employed to develop a quantitative measure of dynamical coupling in quantum transition events, such as reactive tunneling and resonant energy transfer. We utilize the method to study condensed-phase proton-coupled electron transfer (PCET) reactions and to…
▽ More
We introduce a simple method for characterizing reactive pathways in quantum systems. Flux auto- correlation and cross-correlation functions are employed to develop a quantitative measure of dynamical coupling in quantum transition events, such as reactive tunneling and resonant energy transfer. We utilize the method to study condensed-phase proton-coupled electron transfer (PCET) reactions and to determine the relative importance of competing concerted and sequential reaction pathways. Results presented here include numerically exact quantum dynamics simulations for model condensed-phase PCET reactions. This work demonstrates the applicability of the new method for the analysis of both approximate and exact quantum dynamics simulations.
△ Less
Submitted 15 April, 2012;
originally announced April 2012.
-
Direct simulation of electron transfer using ring polymer molecular dynamics: Comparison with semiclassical instanton theory and exact quantum methods
Authors:
Artur R. Menzeleev,
Nandini Ananth,
Thomas F. Miller III
Abstract:
The use of ring polymer molecular dynamics (RPMD) for the direct simulation of electron transfer (ET) reaction dynamics is analyzed in the context of Marcus theory, semiclassical instanton theory, and exact quantum dynamics approaches. For both fully atomistic and system-bath representations of condensed-phase ET, we demonstrate that RPMD accurately predicts both ET reaction rates and mechanisms t…
▽ More
The use of ring polymer molecular dynamics (RPMD) for the direct simulation of electron transfer (ET) reaction dynamics is analyzed in the context of Marcus theory, semiclassical instanton theory, and exact quantum dynamics approaches. For both fully atomistic and system-bath representations of condensed-phase ET, we demonstrate that RPMD accurately predicts both ET reaction rates and mechanisms throughout the normal and activationless regimes of the thermodynamic driving force. Analysis of the ensemble of reactive RPMD trajectories reveals the solvent reorganization mechanism for ET that is anticipated in the Marcus rate theory, and the accuracy of the RPMD rate calculation is understood in terms of its exact description of statistical fluctuations and its formal connection to semiclassical instanton theory for deep-tunneling processes. In the inverted regime of the thermodynamic driving force, neither RPMD nor a related formulation of semiclassical instanton theory capture the characteristic turnover in the reaction rate; comparison with exact quantum dynamics simulations reveals that these methods provide inadequate quantization of the real-time electronic-state dynamics in the inverted regime.
△ Less
Submitted 25 July, 2011;
originally announced July 2011.