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

Journal of Chemical Theory and Computation

Small Matrix Path Integral for System-Bath Dynamics

Nancy Makri

Published Jun 15, 2020• 95 citations
Publisher

Abstract

A small matrix decomposition of the path integral expression (SMatPI) that yields the reduced density matrix of a system interacting with a dissipative harmonic bath is obtained by recursively spreading the entangled influence functional terms over longer time intervals while simultaneously decreasing their magnitude until these terms become negligible. This allows summation over the path integral variables one by one through multiplication of small matrices with dimension equal to that of the bare system. The theoretical framework of the decomposition is described using a diagrammatic approach. Analytical and numerical calculations show that the necessary time length for the temporal entanglement to become negligible is practically the same as the bath-induced memory. The properties and structure of the propagator matrices are discussed, and applications to multistate systems are presented.

DOI
10.1021/acs.jctc.0c00039

Related papers

Small matrix disentanglement of the path integral: Overcoming the exponential tensor scaling with memory length

Nancy Makri

The discretized path integral expression for the reduced density matrix (RDM) of a system interacting with a dissipative harmonic bath is fully entangled because of influence functional terms that couple the variables at different time points. The iterative decomposition of the path integral, which exploits the finite length of influence functional memory, involves a tensor propagator whose size grows exponentially with the memory length. The present Communication disentangles the path integral by recursively spreading the temporal entanglement over longer path segments, while decreasing its contribution. Eventually, the entangled term becomes sufficiently small and may be neglected, leading to iterative propagation of the RDM through simple multiplication of matrices whose size is equal to that of the bare system. It is found that the temporal entanglement length is practically equal to the bath-induced memory length. The small matrix decomposition of the path integral (SMatPI) is stable and very efficient, extending the applicability of numerically exact real-time path integral methods to multi-state systems.

The Journal of Chemical Physics•Jan 27, 2020•101 citations

Kink Sum for Long-Memory Small Matrix Path Integral Dynamics

Nancy Makri

The small matrix decomposition of the real-time path integral (SMatPI) allows for numerically exact and efficient propagation of the reduced density matrix (RDM) for system-bath Hamiltonians. Its high efficiency lies in the small size of the SMatPI matrices employed in the iterative algorithm, whose size is equal to that of the full RDM. By avoiding the storage and multiplication of large tensors, the SMatPI algorithm is applicable in multistate systems under a variety of conditions. The main computational effort is the evaluation of path sums within the entangled memory length to construct the SMatPI matrices. A number of methods are available for this task, each with its own favorable parameter regime, but calculations with strong system-bath coupling and long memory at low temperatures remain out of reach. The present paper evaluates the path sums by binning the paths (in forward time only) based on their amplitudes, which depend on the number and type of kinks they contain. The algorithm is very efficient, leading to a dramatic acceleration of path sums and significantly extending the accessible memory length in the most challenging regimes.

The Journal of Physical Chemistry B•Mar 4, 2024•12 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

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

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 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
A small matrix decomposition of the path integral expression (SMatPI) that yields the reduced density matrix of a system interacting with a dissipative harmonic bath is obtained by recursively spreading the entangled influence functional terms over longer time intervals while simultaneously decreasing their magnitude until these terms become negligible. This allows summation over the path integral variables one by one through multiplication of small matrices with dimension equal to that of the bare system. The theoretical framework of the decomposition is described using a diagrammatic approach. Analytical and numerical calculations show that the necessary time length for the temporal entanglement to become negligible is practically the same as the bath-induced memory. The properties and structure of the propagator matrices are discussed, and applications to multistate systems are presented.
The discretized path integral expression for the reduced density matrix (RDM) of a system interacting with a dissipative harmonic bath is fully entangled because of influence functional terms that couple the variables at different time points. The iterative decomposition of the path integral, which exploits the finite length of influence functional memory, involves a tensor propagator whose size grows exponentially with the memory length. The present Communication disentangles the path integral by recursively spreading the temporal entanglement over longer path segments, while decreasing its contribution. Eventually, the entangled term becomes sufficiently small and may be neglected, leading to iterative propagation of the RDM through simple multiplication of matrices whose size is equal to that of the bare system. It is found that the temporal entanglement length is practically equal to the bath-induced memory length. The small matrix decomposition of the path integral (SMatPI) is stable and very efficient, extending the applicability of numerically exact real-time path integral methods to multi-state systems.
The small matrix decomposition of the real-time path integral (SMatPI) allows for numerically exact and efficient propagation of the reduced density matrix (RDM) for system-bath Hamiltonians. Its high efficiency lies in the small size of the SMatPI matrices employed in the iterative algorithm, whose size is equal to that of the full RDM. By avoiding the storage and multiplication of large tensors, the SMatPI algorithm is applicable in multistate systems under a variety of conditions. The main computational effort is the evaluation of path sums within the entangled memory length to construct the SMatPI matrices. A number of methods are available for this task, each with its own favorable parameter regime, but calculations with strong system-bath coupling and long memory at low temperatures remain out of reach. The present paper evaluates the path sums by binning the paths (in forward time only) based on their amplitudes, which depend on the number and type of kinks they contain. The algorithm is very efficient, leading to a dramatic acceleration of path sums and significantly extending the accessible memory length in the most challenging regimes.
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 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 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.