Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

The Journal of Chemical Physics

Quantum corrections to classical evaluation of nonadiabatic transition rates

Hyojoon Kim, Peter J. Rossky

Published Aug 8, 2006• 7 citations
Publisher

Abstract

A recently developed quantum correction approach is applied to evaluating the nonadiabatic quantum-mechanical transition rate between Born-Oppenheimer states of a subsystem embedded in a thermal bath of harmonic oscillators. In the first-order perturbation theory, the nonadiabatic rate can be expressed in terms of a quantum-mechanical correlation function, which can be estimated directly from classical data. Application to a popular spin-boson model shows that our results are in excellent agreement with the exact quantum-mechanical results.

DOI
10.1063/1.2217443

Related papers

Numerical sampling of nonadiabatic dynamics of quantum-classical systems.

Daniel A. Uken

The simulation of the dynamics of quantum systems is very difficult, due to the fact that, in general, it cannot be calculated exactly for interacting many-body systems. Brute force simulations of quantum dynamics are simply not feasible, and approximations need to be made. In many instances a quantum system can be approximated as a quantum-classical system, where only a subsystem of interest is treated quantum mechanically, and the rest is considered as a classical bath. When energy is free to be exchanged between the subsystem and its environment, the dynamics that occur is said to be nonadiabatic. This type of dynamics is challenging to calculate on a computer, as it can lead to large statistical errors at long times. Hence, there is a need for improved algorithms for nonadiabatic dynamics. In this thesis, a recently introduced nonadiabatic sampling scheme [A. Sergi and F. Petruccione, Phys. Rev. E 81, 032101 (2010)] is used to calculate the long-time dynamics of a model system comprising a quantum spin coupled to a bath of harmonic oscillators. Also, various technical aspects of the algorithm are investigated.

•Jan 1, 2010•0 citations

Surface-hopping dynamics of a spin-boson system

Dónal MacKernan, Giovanni Ciccotti, Raymond Kapral

The spin-boson model is solved within the framework of quantum-classical dynamics using our recently-developed surface-hopping scheme. The quantum-classical equation of motion is expressed in an adiabatic basis and its solution is constructed from an ensemble of trajectories which undergo nonadiabatic transitions and evolve coherently on the adiabatic surfaces. Details of the algorithm for the simulation of the dynamics are presented and the method of simple Monte Carlo sampling used to evaluate the expectation values of observables is discussed. The simulation method is applied to a spin-boson system with a harmonic bath composed of ten oscillators with an Ohmic spectral density. For the spin-boson model the present implementation of quantum-classical dynamics is exact and the results of our surface-hopping simulations are in accord with previous numerically exact results for this model.

The Journal of Chemical Physics•Feb 8, 2002•101 citations

Trotter-Based Simulation of Quantum-Classical Dynamics

Dónal MacKernan, Giovanni Ciccotti, Raymond Kapral

Quantum rate processes in condensed phase systems are often computed by combining quantum and classical descriptions of the dynamics. An algorithm for simulating the quantum-classical Liouville equation, which describes the dynamics of a quantum subsystem coupled to a classical bath, is presented in this paper. The algorithm is based on a Trotter decomposition of the quantum-classical propagator, in conjunction with Monte Carlo sampling of quantum transitions, to yield a surface-hopping representation of the dynamics. An expression for the nonadiabatic propagator that is responsible for quantum transitions and associated bath momentum changes is derived in a form that is convenient for Monte Carlo sampling and exactly conserves the total energy of the system in individual trajectories. The expectation values of operators or quantum correlation functions can be evaluated by initial sampling of quantum states and use of quantum-classical Liouville dynamics for the time evolution. The algorithm is tested by calculations on the spin-boson model, for which exact quantum results are available, and is shown to reproduce the exact results for stronger nonadiabatic coupling and much longer times using fewer trajectories than other schemes for simulating quantum-classical Liouville dynamics.

The Journal of Physical Chemistry B•Dec 22, 2007•105 citations

A Study of Quantum-classical Dynamics in the Mapping Basis

Ali M. Nassimi

Solving quantum dynamics is an exponentially difficult problem. Thus, an exact numerical solution is inaccessible for any condensed matter system.\nA promising approach is to divide the system into a quantum subsystem containing degrees of freedom which are of greater interest or those which have more profound quantum character (e.g., have smaller mass) and a classical bath containing the rest of the system.\nImposing such a partition and treating the bath classically results in quantum-classical dynamics. The quantum-classical Liouville equation is a general equation in the Hilbert space of quantum degrees of freedom while it resides in the phase space of the classical degrees of freedom.\nAny numerical solution to this equation requires representation of the quantum subsystem in some basis. Solutions to this equation have been already proposed in the subsystem, adiabatic and force bases, each with its own cons and pros.\n\nIn this work, the quantum-classical equations of motion are cast in the subsystem basis and subsequently mapped to a number of fictitious harmonic oscillators.\nThe result is quantum-classical dynamics in the mapping basis which treats both quantum and classical degrees of freedom on the same footing, i.e., in phase space. Neglecting a portion of the back reaction of the quantum-subsystem to classical bath results in an expression for the time evolution of an operator (density matrix) equal to its Poisson bracket with the Hamiltonian.\nThis equation can be solved in terms of characteristics to provide a computationally tractable method for calculating quantum-classical dynamical properties. The expressions for expectation values and correlation functions in this formalism are derived.\nCalculations on spin-boson system, barrier crossing models---the so called Tully models---and the Fenna-Mathews-Olson pigments show very good agreement between the results of this method and numerical solutions to the Schrödinger equation.

•Aug 31, 2011•0 citations

Variance of the energy of a quantum system in a time-dependent perturbation: Determination by nonadiabatic transition probabilities

Anirban Mandal, Katharine L. C. Hunt

For a quantum system in a time-dependent perturbation, we prove that the variance in the energy depends entirely on the nonadiabatic transition probability amplitudes bk(t). Landau and Lifshitz introduced the nonadiabatic coefficients for the excited states of a perturbed quantum system by integrating by parts in Dirac’s expressions for the coefficients ck(1)(t) of the excited states to first order in the perturbation. This separates ck(1)(t) for each state into an adiabatic term ak(1)(t) and a nonadiabatic term bk(1)(t). The adiabatic term follows the adiabatic theorem of Born and Fock; it reflects the adjustment of the initial state to the perturbation without transitions. If the response to a time-dependent perturbation is entirely adiabatic, the variance in the energy is zero. The nonadiabatic term bk(1)(t) represents actual excitations away from the initial state. As a key result of the current work, we derive the variance in the energy of the quantum system and all of the higher moments of the energy distribution using the values of |bk(t)|2 for each of the excited states along with the energy differences between the excited states and the ground state. We prove that the same variance (through second order) is obtained in terms of Dirac’s excited-state coefficients ck(t). We show that the results from a standard statistical analysis of the variance are consistent with the quantum results if the probability of excitation Pk is set equal to |bk(t)|2, but not if the probability of excitation is set equal to |ck(t)|2. We illustrate the differences between the variances calculated with the two different forms of Pk for vibration–rotation transitions of HCl in the gas phase.

The Journal of Chemical Physics•Mar 12, 2020•2 citations

Nonadiabatic transition probabilities for quantum systems in electromagnetic fields: Dephasing and population relaxation due to contact with a bath

Sara D. Jovanovski, Anirban Mandal, Katharine L. C. Hunt

We contrast Dirac's theory of transition probabilities and the theory of nonadiabatic transition probabilities, applied to a perturbed system that is coupled to a bath. In Dirac's analysis, the presence of an excited state |k0⟩ in the time-dependent wave function constitutes a transition. In the nonadiabatic theory, a transition occurs when the wave function develops a term that is not adiabatically connected to the initial state. Landau and Lifshitz separated Dirac's excited-state coefficients into a term that follows the adiabatic theorem of Born and Fock and a nonadiabatic term that represents excitation across an energy gap. If the system remains coherent, the two approaches are equivalent. However, differences between the two approaches arise when coupling to a bath causes dephasing, a situation that was not treated by Dirac. For two-level model systems in static electric fields, we add relaxation terms to the Liouville equation for the time derivative of the density matrix. We contrast the results obtained from the two theories. In the analysis based on Dirac's transition probabilities, the steady state of the system is not an equilibrium state; also, the steady-state population ρkk,s increases with increasing strength of the perturbation and its value depends on the dephasing time T2. In the nonadiabatic theory, the system evolves to the thermal equilibrium with the bath. The difference is not simply due to the choice of basis because the difference remains when the results are transformed to a common basis.

The Journal of Chemical Physics•Apr 26, 2023•2 citations
A recently developed quantum correction approach is applied to evaluating the nonadiabatic quantum-mechanical transition rate between Born-Oppenheimer states of a subsystem embedded in a thermal bath of harmonic oscillators. In the first-order perturbation theory, the nonadiabatic rate can be expressed in terms of a quantum-mechanical correlation function, which can be estimated directly from classical data. Application to a popular spin-boson model shows that our results are in excellent agreement with the exact quantum-mechanical results.
The simulation of the dynamics of quantum systems is very difficult, due to the fact that, in general, it cannot be calculated exactly for interacting many-body systems. Brute force simulations of quantum dynamics are simply not feasible, and approximations need to be made. In many instances a quantum system can be approximated as a quantum-classical system, where only a subsystem of interest is treated quantum mechanically, and the rest is considered as a classical bath. When energy is free to be exchanged between the subsystem and its environment, the dynamics that occur is said to be nonadiabatic. This type of dynamics is challenging to calculate on a computer, as it can lead to large statistical errors at long times. Hence, there is a need for improved algorithms for nonadiabatic dynamics. In this thesis, a recently introduced nonadiabatic sampling scheme [A. Sergi and F. Petruccione, Phys. Rev. E 81, 032101 (2010)] is used to calculate the long-time dynamics of a model system comprising a quantum spin coupled to a bath of harmonic oscillators. Also, various technical aspects of the algorithm are investigated.
The spin-boson model is solved within the framework of quantum-classical dynamics using our recently-developed surface-hopping scheme. The quantum-classical equation of motion is expressed in an adiabatic basis and its solution is constructed from an ensemble of trajectories which undergo nonadiabatic transitions and evolve coherently on the adiabatic surfaces. Details of the algorithm for the simulation of the dynamics are presented and the method of simple Monte Carlo sampling used to evaluate the expectation values of observables is discussed. The simulation method is applied to a spin-boson system with a harmonic bath composed of ten oscillators with an Ohmic spectral density. For the spin-boson model the present implementation of quantum-classical dynamics is exact and the results of our surface-hopping simulations are in accord with previous numerically exact results for this model.
Quantum rate processes in condensed phase systems are often computed by combining quantum and classical descriptions of the dynamics. An algorithm for simulating the quantum-classical Liouville equation, which describes the dynamics of a quantum subsystem coupled to a classical bath, is presented in this paper. The algorithm is based on a Trotter decomposition of the quantum-classical propagator, in conjunction with Monte Carlo sampling of quantum transitions, to yield a surface-hopping representation of the dynamics. An expression for the nonadiabatic propagator that is responsible for quantum transitions and associated bath momentum changes is derived in a form that is convenient for Monte Carlo sampling and exactly conserves the total energy of the system in individual trajectories. The expectation values of operators or quantum correlation functions can be evaluated by initial sampling of quantum states and use of quantum-classical Liouville dynamics for the time evolution. The algorithm is tested by calculations on the spin-boson model, for which exact quantum results are available, and is shown to reproduce the exact results for stronger nonadiabatic coupling and much longer times using fewer trajectories than other schemes for simulating quantum-classical Liouville dynamics.
Solving quantum dynamics is an exponentially difficult problem. Thus, an exact numerical solution is inaccessible for any condensed matter system.\nA promising approach is to divide the system into a quantum subsystem containing degrees of freedom which are of greater interest or those which have more profound quantum character (e.g., have smaller mass) and a classical bath containing the rest of the system.\nImposing such a partition and treating the bath classically results in quantum-classical dynamics. The quantum-classical Liouville equation is a general equation in the Hilbert space of quantum degrees of freedom while it resides in the phase space of the classical degrees of freedom.\nAny numerical solution to this equation requires representation of the quantum subsystem in some basis. Solutions to this equation have been already proposed in the subsystem, adiabatic and force bases, each with its own cons and pros.\n\nIn this work, the quantum-classical equations of motion are cast in the subsystem basis and subsequently mapped to a number of fictitious harmonic oscillators.\nThe result is quantum-classical dynamics in the mapping basis which treats both quantum and classical degrees of freedom on the same footing, i.e., in phase space. Neglecting a portion of the back reaction of the quantum-subsystem to classical bath results in an expression for the time evolution of an operator (density matrix) equal to its Poisson bracket with the Hamiltonian.\nThis equation can be solved in terms of characteristics to provide a computationally tractable method for calculating quantum-classical dynamical properties. The expressions for expectation values and correlation functions in this formalism are derived.\nCalculations on spin-boson system, barrier crossing models---the so called Tully models---and the Fenna-Mathews-Olson pigments show very good agreement between the results of this method and numerical solutions to the Schrödinger equation.
For a quantum system in a time-dependent perturbation, we prove that the variance in the energy depends entirely on the nonadiabatic transition probability amplitudes bk(t). Landau and Lifshitz introduced the nonadiabatic coefficients for the excited states of a perturbed quantum system by integrating by parts in Dirac’s expressions for the coefficients ck(1)(t) of the excited states to first order in the perturbation. This separates ck(1)(t) for each state into an adiabatic term ak(1)(t) and a nonadiabatic term bk(1)(t). The adiabatic term follows the adiabatic theorem of Born and Fock; it reflects the adjustment of the initial state to the perturbation without transitions. If the response to a time-dependent perturbation is entirely adiabatic, the variance in the energy is zero. The nonadiabatic term bk(1)(t) represents actual excitations away from the initial state. As a key result of the current work, we derive the variance in the energy of the quantum system and all of the higher moments of the energy distribution using the values of |bk(t)|2 for each of the excited states along with the energy differences between the excited states and the ground state. We prove that the same variance (through second order) is obtained in terms of Dirac’s excited-state coefficients ck(t). We show that the results from a standard statistical analysis of the variance are consistent with the quantum results if the probability of excitation Pk is set equal to |bk(t)|2, but not if the probability of excitation is set equal to |ck(t)|2. We illustrate the differences between the variances calculated with the two different forms of Pk for vibration–rotation transitions of HCl in the gas phase.
We contrast Dirac's theory of transition probabilities and the theory of nonadiabatic transition probabilities, applied to a perturbed system that is coupled to a bath. In Dirac's analysis, the presence of an excited state |k0⟩ in the time-dependent wave function constitutes a transition. In the nonadiabatic theory, a transition occurs when the wave function develops a term that is not adiabatically connected to the initial state. Landau and Lifshitz separated Dirac's excited-state coefficients into a term that follows the adiabatic theorem of Born and Fock and a nonadiabatic term that represents excitation across an energy gap. If the system remains coherent, the two approaches are equivalent. However, differences between the two approaches arise when coupling to a bath causes dephasing, a situation that was not treated by Dirac. For two-level model systems in static electric fields, we add relaxation terms to the Liouville equation for the time derivative of the density matrix. We contrast the results obtained from the two theories. In the analysis based on Dirac's transition probabilities, the steady state of the system is not an equilibrium state; also, the steady-state population ρkk,s increases with increasing strength of the perturbation and its value depends on the dephasing time T2. In the nonadiabatic theory, the system evolves to the thermal equilibrium with the bath. The difference is not simply due to the choice of basis because the difference remains when the results are transformed to a common basis.