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 Physical Chemistry Letters

Small Matrix Quantum-Classical Path Integral

Sohang Kundu, Nancy Makri

Published Apr 13, 2022• 8 citations
Publisher

Abstract

The quantum-classical path integral (QCPI) is a rigorous formulation of nonadiabatic dynamics, where the dynamical interaction between a quantum system and its environment is captured consistently through classical trajectories driven by forces along quantum paths of the system. In this Letter, we develop a small matrix decomposition (SMatQCPI) that eliminates the tensor storage requirements of the iterative QCPI algorithm. In the case of a system coupled to a harmonic bath, SMatQCPI provides fully quantum mechanical propagation, which also reduces the computational cost to that of a single QCPI step. Further, the SMatQCPI matrices only need to account for quantum contributions to decoherence, allowing high efficiency in challenging regimes of incoherent dynamics. Overall, this new composite algorithm combines the best features of two powerful path integral formulations and offers a versatile tool for simulating condensed phase quantum dynamics.

DOI
10.1021/acs.jpclett.2c00668

Related papers

Small Matrix Path Integral for Driven Dissipative Dynamics

Nancy Makri

The small matrix decomposition of the quasi-adiabatic propagator path integral (SMatPI) for a system coupled to a harmonic bath, which accounts for multitime memory correlations in the influence functional without the use of tensors, is extended to include a time-dependent term that drives the system. In the case of a periodic field, the algorithm requires the construction of SMatPI matrices initialized over a short time interval. The SMatPI algorithm circumvents the large array storage of tensor-based iterative path integral decompositions and, in the case of a periodic field, also eliminates the demanding tensor multiplication at each time step, leading to dramatic savings which allow the fully quantum mechanical treatment of multistate systems and long-memory environments.

The Journal of Physical Chemistry A•Nov 23, 2021•29 citations

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

Small Matrix Path Integral with Extended Memory

Nancy Makri

The small matrix decomposition of the path integral (SMatPI) for a discrete system coupled to a harmonic bath expresses the reduced density matrix in terms of matrices whose size is given by the number of states comprising the system, circumventing the large storage requirements of iterative tensor-based algorithms. The present work extends the SMatPI methodology to account for residual memory that exceeds the entanglement length without an increase in computational effort.

Journal of Chemical Theory and Computation•Dec 2, 2020•55 citations

Tree-Based Implementation of the Small Matrix Path Integral for System-Bath Dynamics

Geshuo Wang, Zhenning Cai

The small matrix path integral (SMatPI) method is an efficient numerical approach to simulating the evolution of a quantum system coupled to a harmonic bath. The method relies on a sequence of kernel matrices that defines the non-Markovian dynamics of the quantum system. In the original SMatPI method, these kernels are computed indirectly through the QuAPI method. Instead, we focus on the definition of the kernel matrices and reveal a recurrence relation in these matrices. Using such a relationship, a tree based algorithm (t-SMatPI) is developed, which is shown to be much faster than straightforward computation of the kernel matrices based on their definitions. This algorithm bypasses the step to compute the SMatPI matrices by other path integral methods and provides more understanding of the SMatPI matrices themselves. Meanwhile, it keeps the memory cost and computational cost low. Numerical experiments show that the t-SMatPI algorithm gives exactly the same result as i-QuAPI and SMatPI. In spite of this, our method may indicate some new properties of open quantum systems, and has the potential to be generalized to higher-order numerical schemes.

arXiv:2207.11830•Communications in Computational Physics•Sep 5, 2024•10 citations

Small matrix path integral propagation for long-time quantum dynamics of multistate systems in one and two dimensions

Nancy Makri

The small matrix path integral (SMatPI) methodology allows fully quantum mechanical calculations on system-bath Hamiltonians under a variety of conditions with minimal array storage and is thus applicable to multistate systems. This work focuses on the SMatPI iteration algorithm, which can be computationally prohibitive in systems with many states. By transforming the SMatPI propagation matrices to the desired initial state, the cost of the iterative multiplication is reduced by a factor of n2, where n is the number of system states. This leads to dramatic savings when n≫10 and allows the simulation of large systems over very long propagation times. Applications are presented on models with up to 50 system states, including exciton transfer dynamics in one-dimensional aggregates and a square lattice.

The Journal of Chemical Physics•Feb 17, 2026•1 citation

Path integral methods for quantum dynamics in condensed phases

Tuseeta Banerjee

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.

•Dec 2, 2015•0 citations
The quantum-classical path integral (QCPI) is a rigorous formulation of nonadiabatic dynamics, where the dynamical interaction between a quantum system and its environment is captured consistently through classical trajectories driven by forces along quantum paths of the system. In this Letter, we develop a small matrix decomposition (SMatQCPI) that eliminates the tensor storage requirements of the iterative QCPI algorithm. In the case of a system coupled to a harmonic bath, SMatQCPI provides fully quantum mechanical propagation, which also reduces the computational cost to that of a single QCPI step. Further, the SMatQCPI matrices only need to account for quantum contributions to decoherence, allowing high efficiency in challenging regimes of incoherent dynamics. Overall, this new composite algorithm combines the best features of two powerful path integral formulations and offers a versatile tool for simulating condensed phase quantum dynamics.
The small matrix decomposition of the quasi-adiabatic propagator path integral (SMatPI) for a system coupled to a harmonic bath, which accounts for multitime memory correlations in the influence functional without the use of tensors, is extended to include a time-dependent term that drives the system. In the case of a periodic field, the algorithm requires the construction of SMatPI matrices initialized over a short time interval. The SMatPI algorithm circumvents the large array storage of tensor-based iterative path integral decompositions and, in the case of a periodic field, also eliminates the demanding tensor multiplication at each time step, leading to dramatic savings which allow the fully quantum mechanical treatment of multistate systems and long-memory environments.
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.
The small matrix path integral (SMatPI) method is an efficient numerical approach to simulating the evolution of a quantum system coupled to a harmonic bath. The method relies on a sequence of kernel matrices that defines the non-Markovian dynamics of the quantum system. In the original SMatPI method, these kernels are computed indirectly through the QuAPI method. Instead, we focus on the definition of the kernel matrices and reveal a recurrence relation in these matrices. Using such a relationship, a tree based algorithm (t-SMatPI) is developed, which is shown to be much faster than straightforward computation of the kernel matrices based on their definitions. This algorithm bypasses the step to compute the SMatPI matrices by other path integral methods and provides more understanding of the SMatPI matrices themselves. Meanwhile, it keeps the memory cost and computational cost low. Numerical experiments show that the t-SMatPI algorithm gives exactly the same result as i-QuAPI and SMatPI. In spite of this, our method may indicate some new properties of open quantum systems, and has the potential to be generalized to higher-order numerical schemes.
The small matrix path integral (SMatPI) methodology allows fully quantum mechanical calculations on system-bath Hamiltonians under a variety of conditions with minimal array storage and is thus applicable to multistate systems. This work focuses on the SMatPI iteration algorithm, which can be computationally prohibitive in systems with many states. By transforming the SMatPI propagation matrices to the desired initial state, the cost of the iterative multiplication is reduced by a factor of n2, where n is the number of system states. This leads to dramatic savings when n≫10 and allows the simulation of large systems over very long propagation times. Applications are presented on models with up to 50 system states, including exciton transfer dynamics in one-dimensional aggregates and a square lattice.
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.