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

Path integral methods for quantum dynamics in condensed phases

Tuseeta Banerjee

Published Dec 2, 2015• 0 citations
OpenAlex

Abstract

Quantum mechanical effects play an important role in dynamics of condensed phases. It is also well known that the difficulty of solving the full Schrödinger equation grows exponentially with time. One approach to cope up with this difficulty is to restrict quantum treatment to a few particles (system) and treat the dynamics of the rest of particles (bath) using classical trajectories. Path integral methods due to their trajectory like nature provide excellent tool to develop quantum-classical methods, but to capture long time dynamics, the number of classical trajectories grow exponentially. A part of this work focuses on improving Quantum-Classical Path Integral (QCPI) treatments that allow for larger path integral time step by building “smarter” propagators. These improved solvent-driven reference propagators is developed by incorporating physically motivated approximations to the solvent. When used in QCPI expression they allow convergence with larger time steps leading to an exponential reduction of the number of trajectories required. Further advantages of these propagators include improvements in path filtering techniques, smaller number of path integral steps for achieving the memory decoherence time, and smoothing of the integrand which leads to convergence with fewer Monte Carlo sample points. These ideas have been validated on the spin boson model - a prototypical model to study condensed phase dynamics, which consists of a two-level system coupled to a harmonic bath. The new approach of building the reference propagators combined with their iterative evaluation and filtering is validated using parameters that mimic the first electron transfer in wild-type photosynthetic reaction centres. A real-world dynamical simulation contains various anharmonic effects, Ferrocene (donor) - Ferrocenium (acceptor) system in liquid Benzene was studied where the anharmonic effects of the bath of liquid Benzene on the dynamics of the system were treated using Ensemble Averaged Classical Path. The process of electron or proton transfer occurs between two species having different charges and fluctuations in the solvent plays a key role in determining the rates of such reaction dynamics. To correctly account for the effects of initial distribution of the solvent on the dynamics of the quantum system, the concept of equilibrating the bath to the donor state of the system is used. Two different approaches have been used within the influence functional framework.

Related papers

Quantum‐classical path integral: A rigorous approach to condensed phase dynamics

Nancy Makri

A rigorous methodology that combines a fully quantum mechanical treatment of a small system with a classical trajectory description of a large number of degrees of freedom is discussed in connection with the mechanism of decoherence. The need for mean field approximations in traditional quantum‐classical calculations is removed by abandoning the delocalized wavefunction description of the quantum particle in favor of Feynman's formulation, which is based on local paths. The effects of the environment enter the quantum‐classical path integral (QCPI) expression in terms of phase factors along classical trajectories, whose number grows exponentially with the number of time steps. It is argued that memory loss allows termination of trajectory branching and that the main contribution to decoherence arises from a single classical trajectory (from each sampled initial condition). Exploiting these ideas allows a dramatic reduction in the required number of trajectories, making QCPI calculations feasible. © 2015 Wiley Periodicals, Inc.

International Journal of Quantum Chemistry•Jul 29, 2015•71 citations

Numerical path integral techniques for long time dynamics of quantum dissipative systems

Nancy Makri

Recent progress in numerical methods for evaluating the real-time path integral in dissipative harmonic environments is reviewed. Quasi-adiabatic propagators constructed numerically allow convergence of the path integral with large time increments. Integration of the harmonic bath leads to path integral expressions that incorporate the exact dynamics of the quantum particle along the adiabatic path, with an influence functional that describes nonadiabatic corrections. The resulting quasi-adiabatic propagator path integral is evaluated by efficient system-specific quadratures in most regimes of parameter space, although some cases are handled by grid Monte Carlo sampling. Exploiting the finite span of nonlocal influence functional interactions characteristic of broad condensed phase spectra leads to an iterative scheme for calculating the path integral over arbitrary time lengths. No uncontrolled approximations are introduced, and the resulting methodology converges to the exact quantum result with modest amounts of computational power. Applications to tunneling dynamics in the condensed phase are described.

Journal of Mathematical Physics•May 1, 1995•377 citations

Quantum dynamics for a system coupled to slow baths: On-the-fly filtered propagator method

Eunji Sim

An on-the-fly filtered propagator functional path integral scheme is introduced as an efficient way of calculating the iterative dynamics of complex condensed systems. Time evolution of the reduced density matrix of a dissipative quantum system is evaluated iteratively by filtering the negligible propagator elements at each propagation step. This on-the-fly filtering along with the finiteness of the bath memory dramatically reduces the configuration space to be integrated without losing numerical accuracy of the results. The required computational storage space for the configuration and the weight of the survived path segments increases linearly with the bath memory length. Up to the bath memory time, it is found that a strikingly small fraction of the configurations survives the on-the-fly filtering process and the number of surviving configurations increases algebraically with the propagation time. At times longer than the bath memory time, the number of surviving configurations required for numerically accurate results is essentially saturated and is less than 0.1% of the total number of the configurations. This new scheme is extremely useful for the problems of low-frequency solvents for which the bath memory spans many time steps of the propagation.

The Journal of Chemical Physics•Sep 8, 2001•92 citations

Quantum phenomena in condensed phase systems

Emily R. Dunkel

This thesis develops mixed quantum-classical methods to model dynamical nonequilibrium quantum phenomena in condensed phase systems, and also explores properties of 2-d systems near a quantum critical point using field theoretic techniques and purely classical simulation. Mixed quantum-classical methods involve dividing the system into a quantum sub-system coupled to a classical bath, which is useful for examining the dynamics of systems in the condensed phase since the full quantum treatment of a large number of degrees of freedom is often computationally infeasible. Although the bath is treated classically, its coupling to the quantum system precludes a Newtonian description. We implement an approach that uses the quantum-classical Liouville equation to propagate the quantum subsystem density matrix, which results in an ensemble of surface-hopping trajectories. We test this method on model systems and compare with a linearized path-integral approach. We also develop a novel approach to propagating the reduced density matrix that incorporates ideas from both methods while

•Jan 1, 2008•0 citations

Path integral renormalization for quantum dissipative dynamics with multiple timescales

Nancy Makri

The reduced density matrix of a system coupled to a harmonic dissipative bath can be obtained accurately using an iterative path integral methodology that involves a propagator tensor, whose rank is roughly equal to the number of time steps that span the bath-induced memory. The present work shows that even in situations where the bath is characterized by high-frequency modes (which necessitate a small time step) and slow oscillators (which give rise to very long memory), efficient propagation that takes into account the memory on all time scales is feasible by renormalizing the propagator matrix to include successively larger-length influence functional correlations. At its crudest level, the renormalized propagator is an ordinary matrix, thus the computational demands of the method are very modest, similar to those for propagating the bare system. Yet, test applications on dissipative two-level systems show that the method produces accurate results in demanding situations. By including memory on successively increasing time scales, the method is capable of capturing nonexponential kinetics induced by sluggish solvents.

Molecular Physics•May 10, 2012•30 citations

Path integral formulation of dissipative quantum dynamics

Alexey Novikov

In this thesis the path integral formalism is applied to the calculation of the dynamics of dissipative quantum systems. The time evolution of a system of bilinearly coupled bosonic modes is treated using the real-time path integral technique in coherent-state representation. This method is applied to a damped harmonic oscillator within the Caldeira-Leggett model. In order to get the stationary trajectories the corresponding Lagrangian function is diagonalized and then the path integrals are evaluated by means of the stationary-phase method. The time evolution of the reduced density matrix in the basis of coherent states is given in simple analytic form for weak system-bath coupling, i.e. the so-called rotating-wave terms can be evaluated exactly but the non-rotating-wave terms only in a perturbative manner. The validity range of the rotating-wave approximation is discussed from the viewpoint of spectral equations. In addition, it is shown that systems without initial system-bath correlations can exhibit initial jumps in the population dynamics even for rather weak dissipation. Only with initial correlations the classical trajectories for the system coordinate can be recovered. The path integral formalism in a combined phase-space and coherent-state representation is applied to the problem of curve-crossing dynamics. The system of interest is described by two coupled one-dimensional harmonic potential energy surfaces interacting with a heat bath. The mapping approach is used to rewrite the Lagrangian function of the electronic part of the system. Using the Feynman-Vernon influence-functional method the bath is eliminated whereas the non-Gaussian part of the path integral is treated using the perturbation theory in the small coordinate shift between potential energy surfaces. The vibrational and the population dynamics is considered in a lowest order of the perturbation. The dynamics of a Gaussian wave packet is analyzed along a one-dimensional reaction coordinate. Also the damping rate of coherence in the electronic part of the relevant system is evaluated within the ordinary and variational perturbation theory. The analytic expressions for the rate functions are obtained in the low and high temperature regimes.

•Jan 20, 2005•0 citations
Quantum mechanical effects play an important role in dynamics of condensed phases. It is also well known that the difficulty of solving the full Schrödinger equation grows exponentially with time. One approach to cope up with this difficulty is to restrict quantum treatment to a few particles (system) and treat the dynamics of the rest of particles (bath) using classical trajectories. Path integral methods due to their trajectory like nature provide excellent tool to develop quantum-classical methods, but to capture long time dynamics, the number of classical trajectories grow exponentially. A part of this work focuses on improving Quantum-Classical Path Integral (QCPI) treatments that allow for larger path integral time step by building “smarter” propagators. These improved solvent-driven reference propagators is developed by incorporating physically motivated approximations to the solvent. When used in QCPI expression they allow convergence with larger time steps leading to an exponential reduction of the number of trajectories required. Further advantages of these propagators include improvements in path filtering techniques, smaller number of path integral steps for achieving the memory decoherence time, and smoothing of the integrand which leads to convergence with fewer Monte Carlo sample points. These ideas have been validated on the spin boson model - a prototypical model to study condensed phase dynamics, which consists of a two-level system coupled to a harmonic bath. The new approach of building the reference propagators combined with their iterative evaluation and filtering is validated using parameters that mimic the first electron transfer in wild-type photosynthetic reaction centres. A real-world dynamical simulation contains various anharmonic effects, Ferrocene (donor) - Ferrocenium (acceptor) system in liquid Benzene was studied where the anharmonic effects of the bath of liquid Benzene on the dynamics of the system were treated using Ensemble Averaged Classical Path. The process of electron or proton transfer occurs between two species having different charges and fluctuations in the solvent plays a key role in determining the rates of such reaction dynamics. To correctly account for the effects of initial distribution of the solvent on the dynamics of the quantum system, the concept of equilibrating the bath to the donor state of the system is used. Two different approaches have been used within the influence functional framework.
A rigorous methodology that combines a fully quantum mechanical treatment of a small system with a classical trajectory description of a large number of degrees of freedom is discussed in connection with the mechanism of decoherence. The need for mean field approximations in traditional quantum‐classical calculations is removed by abandoning the delocalized wavefunction description of the quantum particle in favor of Feynman's formulation, which is based on local paths. The effects of the environment enter the quantum‐classical path integral (QCPI) expression in terms of phase factors along classical trajectories, whose number grows exponentially with the number of time steps. It is argued that memory loss allows termination of trajectory branching and that the main contribution to decoherence arises from a single classical trajectory (from each sampled initial condition). Exploiting these ideas allows a dramatic reduction in the required number of trajectories, making QCPI calculations feasible. © 2015 Wiley Periodicals, Inc.
Recent progress in numerical methods for evaluating the real-time path integral in dissipative harmonic environments is reviewed. Quasi-adiabatic propagators constructed numerically allow convergence of the path integral with large time increments. Integration of the harmonic bath leads to path integral expressions that incorporate the exact dynamics of the quantum particle along the adiabatic path, with an influence functional that describes nonadiabatic corrections. The resulting quasi-adiabatic propagator path integral is evaluated by efficient system-specific quadratures in most regimes of parameter space, although some cases are handled by grid Monte Carlo sampling. Exploiting the finite span of nonlocal influence functional interactions characteristic of broad condensed phase spectra leads to an iterative scheme for calculating the path integral over arbitrary time lengths. No uncontrolled approximations are introduced, and the resulting methodology converges to the exact quantum result with modest amounts of computational power. Applications to tunneling dynamics in the condensed phase are described.
An on-the-fly filtered propagator functional path integral scheme is introduced as an efficient way of calculating the iterative dynamics of complex condensed systems. Time evolution of the reduced density matrix of a dissipative quantum system is evaluated iteratively by filtering the negligible propagator elements at each propagation step. This on-the-fly filtering along with the finiteness of the bath memory dramatically reduces the configuration space to be integrated without losing numerical accuracy of the results. The required computational storage space for the configuration and the weight of the survived path segments increases linearly with the bath memory length. Up to the bath memory time, it is found that a strikingly small fraction of the configurations survives the on-the-fly filtering process and the number of surviving configurations increases algebraically with the propagation time. At times longer than the bath memory time, the number of surviving configurations required for numerically accurate results is essentially saturated and is less than 0.1% of the total number of the configurations. This new scheme is extremely useful for the problems of low-frequency solvents for which the bath memory spans many time steps of the propagation.
This thesis develops mixed quantum-classical methods to model dynamical nonequilibrium quantum phenomena in condensed phase systems, and also explores properties of 2-d systems near a quantum critical point using field theoretic techniques and purely classical simulation. Mixed quantum-classical methods involve dividing the system into a quantum sub-system coupled to a classical bath, which is useful for examining the dynamics of systems in the condensed phase since the full quantum treatment of a large number of degrees of freedom is often computationally infeasible. Although the bath is treated classically, its coupling to the quantum system precludes a Newtonian description. We implement an approach that uses the quantum-classical Liouville equation to propagate the quantum subsystem density matrix, which results in an ensemble of surface-hopping trajectories. We test this method on model systems and compare with a linearized path-integral approach. We also develop a novel approach to propagating the reduced density matrix that incorporates ideas from both methods while
The reduced density matrix of a system coupled to a harmonic dissipative bath can be obtained accurately using an iterative path integral methodology that involves a propagator tensor, whose rank is roughly equal to the number of time steps that span the bath-induced memory. The present work shows that even in situations where the bath is characterized by high-frequency modes (which necessitate a small time step) and slow oscillators (which give rise to very long memory), efficient propagation that takes into account the memory on all time scales is feasible by renormalizing the propagator matrix to include successively larger-length influence functional correlations. At its crudest level, the renormalized propagator is an ordinary matrix, thus the computational demands of the method are very modest, similar to those for propagating the bare system. Yet, test applications on dissipative two-level systems show that the method produces accurate results in demanding situations. By including memory on successively increasing time scales, the method is capable of capturing nonexponential kinetics induced by sluggish solvents.
In this thesis the path integral formalism is applied to the calculation of the dynamics of dissipative quantum systems. The time evolution of a system of bilinearly coupled bosonic modes is treated using the real-time path integral technique in coherent-state representation. This method is applied to a damped harmonic oscillator within the Caldeira-Leggett model. In order to get the stationary trajectories the corresponding Lagrangian function is diagonalized and then the path integrals are evaluated by means of the stationary-phase method. The time evolution of the reduced density matrix in the basis of coherent states is given in simple analytic form for weak system-bath coupling, i.e. the so-called rotating-wave terms can be evaluated exactly but the non-rotating-wave terms only in a perturbative manner. The validity range of the rotating-wave approximation is discussed from the viewpoint of spectral equations. In addition, it is shown that systems without initial system-bath correlations can exhibit initial jumps in the population dynamics even for rather weak dissipation. Only with initial correlations the classical trajectories for the system coordinate can be recovered. The path integral formalism in a combined phase-space and coherent-state representation is applied to the problem of curve-crossing dynamics. The system of interest is described by two coupled one-dimensional harmonic potential energy surfaces interacting with a heat bath. The mapping approach is used to rewrite the Lagrangian function of the electronic part of the system. Using the Feynman-Vernon influence-functional method the bath is eliminated whereas the non-Gaussian part of the path integral is treated using the perturbation theory in the small coordinate shift between potential energy surfaces. The vibrational and the population dynamics is considered in a lowest order of the perturbation. The dynamics of a Gaussian wave packet is analyzed along a one-dimensional reaction coordinate. Also the damping rate of coherence in the electronic part of the relevant system is evaluated within the ordinary and variational perturbation theory. The analytic expressions for the rate functions are obtained in the low and high temperature regimes.