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

Modular path integral methodology for real-time quantum dynamics

Nancy Makri

Published Dec 6, 2018• 67 citations
Publisher

Abstract

, 101101 (2018)], it was shown that the locality of interactions in many systems of interest allows a decomposition of the path integral and its evaluation via sequential linking of the paths of relatively small "modules" (e.g., chemical groups or monomers). The present paper describes the modular path integral methodology for simulating dynamical properties by propagating the density matrix in real time. The procedure is first presented for the simple topology of a single-file arrangement of units interacting via nearest neighbor couplings and subsequently extended to the calculation of two-particle correlations in arrays that may also contain some long-range interactions, to the treatment of systems with side chains or cyclic structures, to the simulation of internal dynamics in long organic molecules, and to the modifications required for coupling of one or several units of a system to dissipative environments. Illustrative applications to the dynamics of interacting two-level-systems are presented.

DOI
10.1063/1.5058223

Related papers

Efficient matrix factorisation of the modular path integral for extended systems

Sohang Kundu, Nancy Makri

The modular path integral (MPI) formulation offers a numerically exact, versatile, and highly efficient approach to the quantum dynamics of extended systems characterised by a single-file arrangement of units with short-range interactions, such as hydrocarbons, molecular aggregates with exciton couplings, or spin chains. Rather than propagating the many-particle wavefunction or density matrix in time in the usual stepwise fashion, the MPI scheme proceeds by sequentially integrating over each unit after linking its quantum paths to those of its neighbour and leads to linear scaling with system length. This paper shows that the linking matrix can be further decomposed into a diagonal matrix and a product of low-dimensional matrices, which can be applied sequentially. This factorisation changes the cost scaling from L2 to Llog⁡L, where L is the number of quantum paths, leading to dramatic savings that are analogous to those attained by the fast Fourier transform algorithm and allowing application of the MPI procedure to systems with larger units and to longer propagation times. Applications to spin chains and to electron-vibration dynamics in a large molecular aggregate illustrate the efficiency of the algorithm.

Molecular Physics•Aug 7, 2020•21 citations

Small matrix modular path integral: iterative quantum dynamics in space and time

Nancy Makri

The modular path integral (MPI) formulation for one-dimensional extended systems, such as spin arrays or molecular aggregates, allows evaluation of spin- or exciton-vibration dynamics with effort that scales linearly with the number of units. This work presents a small matrix decomposition of the modular path integral (SMatMPI), which eliminates tensor storage and enables iterative long-time propagation.

Physical Chemistry Chemical Physics•Jan 1, 2021•20 citations

Path integral density matrix dynamics: A method for calculating time-dependent properties in thermal adiabatic and non-adiabatic systems

Scott Habershon

We introduce a new approach for calculating quantum time-correlation functions and time-dependent expectation values in many-body thermal systems; both electronically adiabatic and non-adiabatic cases can be treated. Our approach uses a path integral simulation to sample an initial thermal density matrix; subsequent evolution of this density matrix is equivalent to solution of the time-dependent Schrödinger equation, which we perform using a linear expansion of Gaussian wavepacket basis functions which evolve according to simple classical-like trajectories. Overall, this methodology represents a formally exact approach for calculating time-dependent quantum properties; by introducing approximations into both the imaginary-time and real-time propagations, this approach can be adapted for complex many-particle systems interacting through arbitrary potentials. We demonstrate this method for the spin Boson model, where we find good agreement with numerically exact calculations. We also discuss future directions of improvement for our approach with a view to improving accuracy and efficiency.

The Journal of Chemical Physics•Sep 12, 2013•8 citations

PathSum : A C++ and Fortran suite of fully quantum mechanical real-time path integral methods for (multi-)system + bath dynamics

Sohang Kundu, Nancy Makri

This paper reports the release of PathSum, a new software suite of state-of-the-art path integral methods for studying the dynamics of single or extended systems coupled to harmonic environments. The package includes two modules, suitable for system-bath problems and extended systems comprising many coupled system-bath units, and is offered in C++ and Fortran implementations. The system-bath module offers the recently developed small matrix path integral (SMatPI) and the well-established iterative quasi-adiabatic propagator path integral (i-QuAPI) method for iteration of the reduced density matrix of the system. In the SMatPI module, the dynamics within the entanglement interval can be computed using QuAPI, the blip sum, time evolving matrix product operators, or the quantum-classical path integral method. These methods have distinct convergence characteristics and their combination allows a user to access a variety of regimes. The extended system module provides the user with two algorithms of the modular path integral method, applicable to quantum spin chains or excitonic molecular aggregates. An overview of the methods and code structure is provided, along with guidance on method selection and representative examples.

The Journal of Chemical Physics•Jun 9, 2023•47 citations

Modular path integral for finite-temperature dynamics of extended systems with intramolecular vibrations

Sohang Kundu, Nancy Makri

The modular decomposition of the path integral is a linear-scaling, numerically exact algorithm for calculating dynamical properties of extended systems composed of multilevel units with local couplings. In a recent article, we generalized the method to wavefunction propagation in aggregates characterized by non-diagonal couplings between adjacent units. Here, we extend the method to the calculation of reduced density matrices in aggregates where each unit includes an arbitrary number of coupled harmonic bath modes, which may describe intramolecular normal mode vibrations, at finite temperature. The effects of harmonic modes are included through influence functional factors, which involve analytical expressions that we derive. Representative applications to spin arrays described by the Heisenberg Hamiltonian with dissipative interactions and to J-aggregates of perylene bisimide, where all coupled normal modes are treated explicitly, are presented.

The Journal of Chemical Physics•Jul 28, 2020•23 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
, 101101 (2018)], it was shown that the locality of interactions in many systems of interest allows a decomposition of the path integral and its evaluation via sequential linking of the paths of relatively small "modules" (e.g., chemical groups or monomers). The present paper describes the modular path integral methodology for simulating dynamical properties by propagating the density matrix in real time. The procedure is first presented for the simple topology of a single-file arrangement of units interacting via nearest neighbor couplings and subsequently extended to the calculation of two-particle correlations in arrays that may also contain some long-range interactions, to the treatment of systems with side chains or cyclic structures, to the simulation of internal dynamics in long organic molecules, and to the modifications required for coupling of one or several units of a system to dissipative environments. Illustrative applications to the dynamics of interacting two-level-systems are presented.
The modular path integral (MPI) formulation offers a numerically exact, versatile, and highly efficient approach to the quantum dynamics of extended systems characterised by a single-file arrangement of units with short-range interactions, such as hydrocarbons, molecular aggregates with exciton couplings, or spin chains. Rather than propagating the many-particle wavefunction or density matrix in time in the usual stepwise fashion, the MPI scheme proceeds by sequentially integrating over each unit after linking its quantum paths to those of its neighbour and leads to linear scaling with system length. This paper shows that the linking matrix can be further decomposed into a diagonal matrix and a product of low-dimensional matrices, which can be applied sequentially. This factorisation changes the cost scaling from L2 to Llog⁡L, where L is the number of quantum paths, leading to dramatic savings that are analogous to those attained by the fast Fourier transform algorithm and allowing application of the MPI procedure to systems with larger units and to longer propagation times. Applications to spin chains and to electron-vibration dynamics in a large molecular aggregate illustrate the efficiency of the algorithm.
We introduce a new approach for calculating quantum time-correlation functions and time-dependent expectation values in many-body thermal systems; both electronically adiabatic and non-adiabatic cases can be treated. Our approach uses a path integral simulation to sample an initial thermal density matrix; subsequent evolution of this density matrix is equivalent to solution of the time-dependent Schrödinger equation, which we perform using a linear expansion of Gaussian wavepacket basis functions which evolve according to simple classical-like trajectories. Overall, this methodology represents a formally exact approach for calculating time-dependent quantum properties; by introducing approximations into both the imaginary-time and real-time propagations, this approach can be adapted for complex many-particle systems interacting through arbitrary potentials. We demonstrate this method for the spin Boson model, where we find good agreement with numerically exact calculations. We also discuss future directions of improvement for our approach with a view to improving accuracy and efficiency.
This paper reports the release of PathSum, a new software suite of state-of-the-art path integral methods for studying the dynamics of single or extended systems coupled to harmonic environments. The package includes two modules, suitable for system-bath problems and extended systems comprising many coupled system-bath units, and is offered in C++ and Fortran implementations. The system-bath module offers the recently developed small matrix path integral (SMatPI) and the well-established iterative quasi-adiabatic propagator path integral (i-QuAPI) method for iteration of the reduced density matrix of the system. In the SMatPI module, the dynamics within the entanglement interval can be computed using QuAPI, the blip sum, time evolving matrix product operators, or the quantum-classical path integral method. These methods have distinct convergence characteristics and their combination allows a user to access a variety of regimes. The extended system module provides the user with two algorithms of the modular path integral method, applicable to quantum spin chains or excitonic molecular aggregates. An overview of the methods and code structure is provided, along with guidance on method selection and representative examples.
The modular decomposition of the path integral is a linear-scaling, numerically exact algorithm for calculating dynamical properties of extended systems composed of multilevel units with local couplings. In a recent article, we generalized the method to wavefunction propagation in aggregates characterized by non-diagonal couplings between adjacent units. Here, we extend the method to the calculation of reduced density matrices in aggregates where each unit includes an arbitrary number of coupled harmonic bath modes, which may describe intramolecular normal mode vibrations, at finite temperature. The effects of harmonic modes are included through influence functional factors, which involve analytical expressions that we derive. Representative applications to spin arrays described by the Heisenberg Hamiltonian with dissipative interactions and to J-aggregates of perylene bisimide, where all coupled normal modes are treated explicitly, are presented.
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.