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

arXiv:2308.10540 · The Journal of Chemical Physics

Quantum correlation functions through tensor network path integral

Amartya Bose

Published Dec 5, 2023• 7 citations
arXivHTMLPDF

Abstract

Tensor networks have historically proven to be of great utility in providing compressed representations of wave functions that can be used for the calculation of eigenstates. Recently, it has been shown that a variety of these networks can be leveraged to make real time non-equilibrium simulations of dynamics involving the Feynman-Vernon influence functional more efficient. In this work, a tensor network is developed for non-perturbatively calculating the equilibrium correlation function for open quantum systems using the path integral methodology. These correlation functions are of fundamental importance in calculations of rates of reactions, simulations of response functions and susceptibilities, spectra of systems, etc. The influence of the solvent on the quantum system is incorporated through an influence functional, whose unconventional structure motivates the design of a new optimal matrix product-like operator that can be applied to the so-called path amplitude matrix product state. This complex-time tensor network path integral approach provides an exceptionally efficient representation of the path integral, enabling simulations for larger systems strongly interacting with baths and at lower temperatures out to longer time. The derivation, design, and implementation of this method are discussed along with a wide range of illustrations ranging from rate theory and symmetrized spin correlation functions to simulation of response of the Fenna-Matthews-Olson complex to light.

DOI
10.1063/5.0174338

Related papers

A tensor network representation of path integrals: Implementation and analysis

Amartya Bose, Peter Walters

Tensors with finite correlation afford very compact tensor network representations. A novel tensor network-based decomposition of real-time path integral simulations involving Feynman-Vernon influence functional is introduced. In this tensor network path integral (TNPI) technique, the finite temporarily non-local interactions introduced by the influence functional can be captured very efficiently using matrix product state representation for the path amplitude (PA) tensor. We illustrate this particular TNPI method through various realistic examples, including a charge transfer reaction and an exciton transfer in a dimer. We also show how it is readily applied to systems with greater than two states by simulating a 7-site model of FMO and a molecular wire model. The augmented propagator (AP) TNPI utilizes the symmetries of the problem, leading to accelerated convergence and dramatic reductions of computational effort. We also introduce an approximate method that speeds up propagation beyond the non-local memory length. Furthermore, the structure imposed by the tensor network representation of the PA tensor naturally suggests other factorizations that make simulations for extended systems more efficient. These factorizations would be the subject of future explorations. The flexibility of the AP-TNPI framework makes it a promising new addition to the family of path integral methods for non-equilibrium quantum dynamics.

arXiv:2106.12523•Jun 23, 2021•4 citations

A multisite decomposition of the tensor network path integrals

Amartya Bose, Peter Walters

Tensor network decompositions of path integrals for simulating open quantum systems have recently been proven to be useful. However, these methods scale exponentially with the system size. This makes it challenging to simulate the non-equilibrium dynamics of extended quantum systems coupled with local dissipative environments. In this work, we extend the tensor network path integral (TNPI) framework to efficiently simulate such extended systems. The Feynman-Vernon influence functional is a popular approach used to account for the effect of environments on the dynamics of the system. In order to facilitate the incorporation of the influence functional into a multisite framework (MS-TNPI), we combine a matrix product state (MPS) decomposition of the reduced density tensor of the system along the sites with a corresponding tensor network representation of the time axis to construct an efficient 2D tensor network. The 2D MS-TNPI network, when contracted, yields the time-dependent reduced density tensor of the extended system as an MPS. The algorithm presented is independent of the system Hamiltonian. We outline an iteration scheme to take the simulation beyond the non-Markovian memory introduced by solvents. Applications to spin chains coupled to local harmonic baths are presented; we consider the Ising, XXZ, and Heisenberg models, demonstrating that the presence of local environments can often dissipate the entanglement between the sites. We discuss three factors causing the system to transition from a coherent oscillatory dynamics to a fully incoherent dynamics. The MS-TNPI method is useful for studying a variety of extended quantum systems coupled with solvents.

arXiv:2109.09723•The Journal of Chemical Physics•Jan 10, 2022•51 citations

Exploiting the Causal Tensor Network Structure of Quantum Processes to Efficiently Simulate Non-Markovian Path Integrals

Mathias R. Jørgensen, Felix A. Pollock

In the path integral formulation of the evolution of an open quantum system coupled to a Gaussian, noninteracting environment, the dynamical contribution of the latter is encoded in an object called the influence functional. Here, we relate the influence functional to the process tensor-a more general representation of a quantum stochastic process-describing the evolution. Then, we use this connection to motivate a tensor network algorithm for the simulation of multitime correlations in open systems, building on recent work where the influence functional is represented in terms of time evolving matrix product operators. By exploiting the symmetries of the influence functional, we are able to use our algorithm to achieve orders-of-magnitude improvement in the efficiency of the resulting numerical simulation. Our improved algorithm is then applied to compute exact phonon emission spectra for the spin-boson model with strong coupling, demonstrating a significant divergence from spectra derived under commonly used assumptions of memorylessness.

arXiv:1902.00315•Physical Review Letters•Dec 9, 2019•161 citations

Tensor network simulations of open quantum systems

Dainius Kilda

The nonequilibrium effects of dissipation and drive play a key role in an immense variety of nanoscale and mesoscale quantum systems. To understand the behaviour of open quantum systems, we need accurate methods that capture the influence of the environment on the system, while managing the exponentially large Hilbert space required to describe the system. Tensor network algorithms offer an efficient way to approach this challenge. In this thesis, we develop and apply tensor network techniques to study the dynamics and steady states of various open quantum systems. The first part of the thesis focuses on the driven dissipative many body physics in coupled cavity arrays described by Born-Markov master equations. We extend transfer matrix product operator methods to Liouvillian dynamics, and utilize them to compute dynamical correlation functions and fluorescence spectrum of an infinite coupled cavity array in 1D. We also investigate thermalization, and observe the emergence of a quasi-thermal steady state with a negative effective temperature. In another study, we use infinite projected entangled pair state (iPEPS) methods to compute steady states of coupled cavity lattices in 2D. We find that a straightforward adaptation of iPEPS to Liouvillian dynamics is unstable, contradicting a recent publication in the field. The second part investigates more general systems involving strong couplings and structured environments that induce non-Markovian dynamics. We develop a powerful time-evolving matrix product operator (TEMPO) algorithm that builds on Feynman-Vernon influence functional formalism, and uses matrix product states (MPS) to represent the temporal non-Markovian correlations efficiently. We apply TEMPO to study the localization phase transition of the spin-boson model and the dynamics of two spatially separated two-level systems coupled to a common environment. Finally, we propose the Toblerone TEMPO algorithm, which extends TEMPO to many-body systems interacting with general bosonic environments.

•Mar 3, 2020•0 citations

Tensor network renormalization: Application to dynamic correlation functions and non-Hermitian systems

Ying-Jie Wei, Zheng‐Cheng Gu

In recent years, tensor network renormalization (TNR) has emerged as an efficient and accurate method for studying (1+1)D quantum systems or 2D classical systems using real-space renormalization group (RG) techniques. One notable application of TNR is its ability to extract central charge and conformal scaling dimensions for critical systems. In this paper, we present the implementation of the Loop-TNR algorithm, which allows for the computation of dynamical correlation functions. Our algorithm goes beyond traditional approaches by not only calculating correlations in the spatial direction, where the separation is an integer, but also in the temporal direction, where the time difference can contain decimal values. Our algorithm is designed to handle both imaginary-time and real-time correlations, utilizing a tensor network representation constructed from a path-integral formalism. Additionally, we highlight that the Loop-TNR algorithm can also be applied to investigate critical properties of non-Hermitian systems, an area that was previously inaccessible using density matrix renormalization group(DMRG) and matrix product state(MPS) based algorithms.

arXiv:2311.18785•Physical review. B./Physical review. B•Mar 26, 2026•1 citation

Pairwise connected tensor network representation of path integrals

Amartya Bose

It has been recently shown how the tensorial nature of real-time path integrals (PIs) involving the Feynman-Vernon influence functional can be utilized with matrix product states, taking advantage of the finite length of the bath-induced memory. Tensor networks (TNs) promise to provide a unified language to express the structure of a PI. Here, a generalized TN specifically incorporating the pairwise interaction structure of the influence functional and its invariance with respect to the average forward-backward position or the sojourn value in the form of the blip representation is derived and implemented. This pairwise connected TNPI (PC-TNPI) is illustrated through applications to typical spin-boson problems and explorations of the differences caused by the exact form of the spectral density. The storage and performance scalings are reported, showing the compactness of the representation and the efficiency of the contraction process. Finally, taking advantage of the compressed representation, the viability of using PC-TNPI for simulating multistate problems is demonstrated. The PC-TNPI structure can be shown to yield other TN algorithms currently in use. Consequently, it should be possible to use it as a starting point for deriving other optimized procedures.

arXiv:2106.14934•Physical review. B./Physical review. B•Jan 27, 2022•20 citations
Tensor networks have historically proven to be of great utility in providing compressed representations of wave functions that can be used for the calculation of eigenstates. Recently, it has been shown that a variety of these networks can be leveraged to make real time non-equilibrium simulations of dynamics involving the Feynman-Vernon influence functional more efficient. In this work, a tensor network is developed for non-perturbatively calculating the equilibrium correlation function for open quantum systems using the path integral methodology. These correlation functions are of fundamental importance in calculations of rates of reactions, simulations of response functions and susceptibilities, spectra of systems, etc. The influence of the solvent on the quantum system is incorporated through an influence functional, whose unconventional structure motivates the design of a new optimal matrix product-like operator that can be applied to the so-called path amplitude matrix product state. This complex-time tensor network path integral approach provides an exceptionally efficient representation of the path integral, enabling simulations for larger systems strongly interacting with baths and at lower temperatures out to longer time. The derivation, design, and implementation of this method are discussed along with a wide range of illustrations ranging from rate theory and symmetrized spin correlation functions to simulation of response of the Fenna-Matthews-Olson complex to light.
Tensors with finite correlation afford very compact tensor network representations. A novel tensor network-based decomposition of real-time path integral simulations involving Feynman-Vernon influence functional is introduced. In this tensor network path integral (TNPI) technique, the finite temporarily non-local interactions introduced by the influence functional can be captured very efficiently using matrix product state representation for the path amplitude (PA) tensor. We illustrate this particular TNPI method through various realistic examples, including a charge transfer reaction and an exciton transfer in a dimer. We also show how it is readily applied to systems with greater than two states by simulating a 7-site model of FMO and a molecular wire model. The augmented propagator (AP) TNPI utilizes the symmetries of the problem, leading to accelerated convergence and dramatic reductions of computational effort. We also introduce an approximate method that speeds up propagation beyond the non-local memory length. Furthermore, the structure imposed by the tensor network representation of the PA tensor naturally suggests other factorizations that make simulations for extended systems more efficient. These factorizations would be the subject of future explorations. The flexibility of the AP-TNPI framework makes it a promising new addition to the family of path integral methods for non-equilibrium quantum dynamics.
Tensor network decompositions of path integrals for simulating open quantum systems have recently been proven to be useful. However, these methods scale exponentially with the system size. This makes it challenging to simulate the non-equilibrium dynamics of extended quantum systems coupled with local dissipative environments. In this work, we extend the tensor network path integral (TNPI) framework to efficiently simulate such extended systems. The Feynman-Vernon influence functional is a popular approach used to account for the effect of environments on the dynamics of the system. In order to facilitate the incorporation of the influence functional into a multisite framework (MS-TNPI), we combine a matrix product state (MPS) decomposition of the reduced density tensor of the system along the sites with a corresponding tensor network representation of the time axis to construct an efficient 2D tensor network. The 2D MS-TNPI network, when contracted, yields the time-dependent reduced density tensor of the extended system as an MPS. The algorithm presented is independent of the system Hamiltonian. We outline an iteration scheme to take the simulation beyond the non-Markovian memory introduced by solvents. Applications to spin chains coupled to local harmonic baths are presented; we consider the Ising, XXZ, and Heisenberg models, demonstrating that the presence of local environments can often dissipate the entanglement between the sites. We discuss three factors causing the system to transition from a coherent oscillatory dynamics to a fully incoherent dynamics. The MS-TNPI method is useful for studying a variety of extended quantum systems coupled with solvents.
In the path integral formulation of the evolution of an open quantum system coupled to a Gaussian, noninteracting environment, the dynamical contribution of the latter is encoded in an object called the influence functional. Here, we relate the influence functional to the process tensor-a more general representation of a quantum stochastic process-describing the evolution. Then, we use this connection to motivate a tensor network algorithm for the simulation of multitime correlations in open systems, building on recent work where the influence functional is represented in terms of time evolving matrix product operators. By exploiting the symmetries of the influence functional, we are able to use our algorithm to achieve orders-of-magnitude improvement in the efficiency of the resulting numerical simulation. Our improved algorithm is then applied to compute exact phonon emission spectra for the spin-boson model with strong coupling, demonstrating a significant divergence from spectra derived under commonly used assumptions of memorylessness.
The nonequilibrium effects of dissipation and drive play a key role in an immense variety of nanoscale and mesoscale quantum systems. To understand the behaviour of open quantum systems, we need accurate methods that capture the influence of the environment on the system, while managing the exponentially large Hilbert space required to describe the system. Tensor network algorithms offer an efficient way to approach this challenge. In this thesis, we develop and apply tensor network techniques to study the dynamics and steady states of various open quantum systems. The first part of the thesis focuses on the driven dissipative many body physics in coupled cavity arrays described by Born-Markov master equations. We extend transfer matrix product operator methods to Liouvillian dynamics, and utilize them to compute dynamical correlation functions and fluorescence spectrum of an infinite coupled cavity array in 1D. We also investigate thermalization, and observe the emergence of a quasi-thermal steady state with a negative effective temperature. In another study, we use infinite projected entangled pair state (iPEPS) methods to compute steady states of coupled cavity lattices in 2D. We find that a straightforward adaptation of iPEPS to Liouvillian dynamics is unstable, contradicting a recent publication in the field. The second part investigates more general systems involving strong couplings and structured environments that induce non-Markovian dynamics. We develop a powerful time-evolving matrix product operator (TEMPO) algorithm that builds on Feynman-Vernon influence functional formalism, and uses matrix product states (MPS) to represent the temporal non-Markovian correlations efficiently. We apply TEMPO to study the localization phase transition of the spin-boson model and the dynamics of two spatially separated two-level systems coupled to a common environment. Finally, we propose the Toblerone TEMPO algorithm, which extends TEMPO to many-body systems interacting with general bosonic environments.
In recent years, tensor network renormalization (TNR) has emerged as an efficient and accurate method for studying (1+1)D quantum systems or 2D classical systems using real-space renormalization group (RG) techniques. One notable application of TNR is its ability to extract central charge and conformal scaling dimensions for critical systems. In this paper, we present the implementation of the Loop-TNR algorithm, which allows for the computation of dynamical correlation functions. Our algorithm goes beyond traditional approaches by not only calculating correlations in the spatial direction, where the separation is an integer, but also in the temporal direction, where the time difference can contain decimal values. Our algorithm is designed to handle both imaginary-time and real-time correlations, utilizing a tensor network representation constructed from a path-integral formalism. Additionally, we highlight that the Loop-TNR algorithm can also be applied to investigate critical properties of non-Hermitian systems, an area that was previously inaccessible using density matrix renormalization group(DMRG) and matrix product state(MPS) based algorithms.
It has been recently shown how the tensorial nature of real-time path integrals (PIs) involving the Feynman-Vernon influence functional can be utilized with matrix product states, taking advantage of the finite length of the bath-induced memory. Tensor networks (TNs) promise to provide a unified language to express the structure of a PI. Here, a generalized TN specifically incorporating the pairwise interaction structure of the influence functional and its invariance with respect to the average forward-backward position or the sojourn value in the form of the blip representation is derived and implemented. This pairwise connected TNPI (PC-TNPI) is illustrated through applications to typical spin-boson problems and explorations of the differences caused by the exact form of the spectral density. The storage and performance scalings are reported, showing the compactness of the representation and the efficiency of the contraction process. Finally, taking advantage of the compressed representation, the viability of using PC-TNPI for simulating multistate problems is demonstrated. The PC-TNPI structure can be shown to yield other TN algorithms currently in use. Consequently, it should be possible to use it as a starting point for deriving other optimized procedures.