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:2507.07648 · Physical review. B./Physical review. B

Summing real time Feynman paths of lattice polarons with matrix product states

Qi Gao, Yuan Wan

Published Oct 1, 2025• 1 citation
arXivHTMLPDF

Abstract

We study numerically the real time dynamics of lattice polarons by combining the Feynman path integral and the matrix product state (MPS) approach. By constructing and solving a flow equation, we show that the integrand, viewed as a multivariable function of polaron world line parameters, can be compressed as a low bond dimension MPS, thereby allowing for efficient evaluation of various dynamical observables. We establish the effectiveness of our method by benchmarking the calculated polaron spectral function in one dimension against available results. We further demonstrate its potential by presenting the polaron spectral function in two dimensions and simulating polaron diffusion in both one and two dimensions.

DOI
10.1103/5sdf-fnl5

Related papers

Finite-temperature spectral functions of lattice polarons with process tensors

Qi Gao, Yuan Wan

We develop a process tensor approach for calculating the finite-temperature spectral functions of lattice polarons. We show that the finite-temperature phonon influence functional can be expressed as the elementwise exponential of a low-rank matrix product state, enabling an efficient summation of electron worldlines. An embedding construction embeds worldlines of different propagation times into a single matrix product state. This construction yields the propagator over the full simulated time interval in a single calculation. We benchmark this approach against available finite-temperature density-matrix renormalization group results for the one-dimensional Holstein polaron. We further obtain numerically converged spectra for the two-dimensional Holstein polaron. The spectra reveal a progressive transfer of spectral weight from coherent spectral features into an incoherent background with increasing temperature.

arXiv:2610.27582•Oct 20, 2026•0 citations

Real-time impurity solver using Grassmann time-evolving matrix product operators

Ruofan Chen, Xiansong Xu, Chu Guo

An emergent and promising tensor-network-based impurity solver is to represent the Feynman-Vernon influence functional as a matrix product state, where the bath is integrated out analytically. Here we present an approach to calculate the equilibrium impurity spectral function based on the recently proposed Grassmann time-evolving matrix product operators method. The central idea is to perform a quench from a separable impurity-bath initial state as in the nonequilibrium scenario. The retarded Green's function $G(t+{t}_{0},{t}^{\ensuremath{'}}+{t}_{0})$ is then calculated after an equilibration time ${t}_{0}$ such that the impurity and bath are approximately in thermal equilibrium. There are two major advantages of this method. First, since we focus on real-time dynamics, we do not need to perform the numerically ill-posed analytic continuation as in imaginary- time evolution-based methods. Second, the required bond dimension of the matrix product state in real-time calculations is observed to be much smaller than that in imaginary-time calculations, leading to a significant improvement in numerical efficiency. The accuracy of this method is demonstrated using the single-orbital Anderson impurity model and benchmarked against the continuous-time quantum Monte Carlo method.

arXiv:2401.04880•Physical review. B./Physical review. B•Apr 2, 2024•10 citations

Solving the Peierls-Boltzmann transport equation with matrix product states

Anonymous, Hirad Alipanah, Juan José Mendoza-Arenas

The Peierls-Boltzmann transport equation (PBE), which governs non-equilibrium phonon transport, suffers from the curse of dimensionality due to its high-dimensional phase space including both real and modal spaces. We explore the use of matrix product states (MPS) for numerical simulation of the PBE. We show that an MPS configuration based on scattering events combined with a dimensionless form of the solution can drastically increase the locality of correlations between tensors in the MPS representation, enhancing its effectiveness in dimension reduction. We further examine the effects of index ordering in an MPS and find that the highest locality is achieved when tensor chains associated with both real and modal spaces are connected from the coarsest grid to each other in the center of the MPS. Using this optimal configuration and a solver inspired by the density matrix renormalization group, we solve the PBE discretized by a finite volume method (FVM). The solution is obtained for crystalline silicon across ballistic, quasi-ballistic, and diffusive transport regimes. An MPS truncated to the compression ratio of $10^{-3}$ suffices to reproduce reference solutions with high fidelity. The computational cost scales sublinearly with the number of grid points in both real and modal spaces, achieving roughly an order of magnitude reduction in computational time compared to the FVM with sparse matrix operation.

arXiv:2604.06153•Physical review. B./Physical review. B•Aug 18, 2026•0 citations

Momentum-resolved time evolution with matrix product states

Maarten Van Damme, Laurens Vanderstraeten

We introduce a method based on matrix product states (MPS) for computing spectral functions of (quasi-) one-dimensional spin chains, working directly in momentum space in the thermodynamic limit. We simulate the time evolution after applying a momentum operator to an MPS ground state by working with the momentum superposition of a window MPS. We show explicitly for the spin-1 Heisenberg chain that the growth of entanglement is smaller in momentum space, even inside a two-particle continuum, such that we can attain very accurate spectral functions with relatively small bond dimension. We apply our method to compute spectral lineshapes of the gapless XXZ chain and the square-lattice ${J}_{1}\ensuremath{-}{J}_{2}$ Heisenberg model on a six-leg cylinder.

arXiv:2201.07314•Physical review. B./Physical review. B•May 23, 2022•7 citations

Grassmann time-evolving matrix product operators: An efficient numerical approach for fermionic path integral simulations

Xiansong Xu, Chu Guo, Ruofan Chen

Developing numerical exact solvers for open quantum systems is a challenging task due to the non-perturbative and non-Markovian nature when coupling to structured environments. The Feynman-Vernon influence functional approach is a powerful analytical tool to study the dynamics of open quantum systems. Numerical treatments of the influence functional including the quasi-adiabatic propagator technique and the tensor-network-based time-evolving matrix product operator method have proven to be efficient in studying open quantum systems with bosonic environments. However, the numerical implementation of the fermionic path integral suffers from the Grassmann algebra involved. In this work, we present a detailed introduction to the Grassmann time-evolving matrix product operator method for fermionic open quantum systems. In particular, we introduce the concepts of Grassmann tensor, signed matrix product operator, and Grassmann matrix product state to handle the Grassmann path integral. Using the single-orbital Anderson impurity model as an example, we review the numerical benchmarks for structured fermionic environments for real-time nonequilibrium dynamics, real-time and imaginary-time equilibration dynamics, and its application as an impurity solver. These benchmarks show that our method is a robust and promising numerical approach to study strong coupling physics and non-Markovian dynamics. It can also serve as an alternative impurity solver to study strongly correlated quantum matter with dynamical mean-field theory.

arXiv:2410.11541•The Journal of Chemical Physics•Oct 15, 2024•6 citations

A first introduction to Matrix Product State algorithms for the integration of Lindblad equation

Christophe Chatelain

In this introductory review, we present and compare four algorithms for the numerical integration of the Lindblad equation for one-dimensional quantum lattice systems. All four methods are based on Matrix Product State representations and can be viewed as extensions of the Time-Evolving Block Decimation (TEBD) algorithm to open quantum systems. Two approaches directly integrate the vectorized Lindblad equation, one of them explicitly enforcing the positivity of the density matrix. The other two rely on stochastic unravelings of the Lindblad equation, namely the quantum trajectory and quantum state diffusion approaches. We discuss the principles, numerical implementation, accuracy, and computational efficiency of the different methods, and benchmark them against an exactly solvable free fermion model.

cond-mat.stat-mechCross-listedarXiv:2609.18841•Sep 17, 2026•0 citations
We study numerically the real time dynamics of lattice polarons by combining the Feynman path integral and the matrix product state (MPS) approach. By constructing and solving a flow equation, we show that the integrand, viewed as a multivariable function of polaron world line parameters, can be compressed as a low bond dimension MPS, thereby allowing for efficient evaluation of various dynamical observables. We establish the effectiveness of our method by benchmarking the calculated polaron spectral function in one dimension against available results. We further demonstrate its potential by presenting the polaron spectral function in two dimensions and simulating polaron diffusion in both one and two dimensions.
We develop a process tensor approach for calculating the finite-temperature spectral functions of lattice polarons. We show that the finite-temperature phonon influence functional can be expressed as the elementwise exponential of a low-rank matrix product state, enabling an efficient summation of electron worldlines. An embedding construction embeds worldlines of different propagation times into a single matrix product state. This construction yields the propagator over the full simulated time interval in a single calculation. We benchmark this approach against available finite-temperature density-matrix renormalization group results for the one-dimensional Holstein polaron. We further obtain numerically converged spectra for the two-dimensional Holstein polaron. The spectra reveal a progressive transfer of spectral weight from coherent spectral features into an incoherent background with increasing temperature.
An emergent and promising tensor-network-based impurity solver is to represent the Feynman-Vernon influence functional as a matrix product state, where the bath is integrated out analytically. Here we present an approach to calculate the equilibrium impurity spectral function based on the recently proposed Grassmann time-evolving matrix product operators method. The central idea is to perform a quench from a separable impurity-bath initial state as in the nonequilibrium scenario. The retarded Green's function $G(t+{t}_{0},{t}^{\ensuremath{'}}+{t}_{0})$ is then calculated after an equilibration time t0{t}_{0}t0​ such that the impurity and bath are approximately in thermal equilibrium. There are two major advantages of this method. First, since we focus on real-time dynamics, we do not need to perform the numerically ill-posed analytic continuation as in imaginary- time evolution-based methods. Second, the required bond dimension of the matrix product state in real-time calculations is observed to be much smaller than that in imaginary-time calculations, leading to a significant improvement in numerical efficiency. The accuracy of this method is demonstrated using the single-orbital Anderson impurity model and benchmarked against the continuous-time quantum Monte Carlo method.
The Peierls-Boltzmann transport equation (PBE), which governs non-equilibrium phonon transport, suffers from the curse of dimensionality due to its high-dimensional phase space including both real and modal spaces. We explore the use of matrix product states (MPS) for numerical simulation of the PBE. We show that an MPS configuration based on scattering events combined with a dimensionless form of the solution can drastically increase the locality of correlations between tensors in the MPS representation, enhancing its effectiveness in dimension reduction. We further examine the effects of index ordering in an MPS and find that the highest locality is achieved when tensor chains associated with both real and modal spaces are connected from the coarsest grid to each other in the center of the MPS. Using this optimal configuration and a solver inspired by the density matrix renormalization group, we solve the PBE discretized by a finite volume method (FVM). The solution is obtained for crystalline silicon across ballistic, quasi-ballistic, and diffusive transport regimes. An MPS truncated to the compression ratio of 10−310^{-3}10−3 suffices to reproduce reference solutions with high fidelity. The computational cost scales sublinearly with the number of grid points in both real and modal spaces, achieving roughly an order of magnitude reduction in computational time compared to the FVM with sparse matrix operation.
We introduce a method based on matrix product states (MPS) for computing spectral functions of (quasi-) one-dimensional spin chains, working directly in momentum space in the thermodynamic limit. We simulate the time evolution after applying a momentum operator to an MPS ground state by working with the momentum superposition of a window MPS. We show explicitly for the spin-1 Heisenberg chain that the growth of entanglement is smaller in momentum space, even inside a two-particle continuum, such that we can attain very accurate spectral functions with relatively small bond dimension. We apply our method to compute spectral lineshapes of the gapless XXZ chain and the square-lattice ${J}_{1}\ensuremath{-}{J}_{2}$ Heisenberg model on a six-leg cylinder.
Developing numerical exact solvers for open quantum systems is a challenging task due to the non-perturbative and non-Markovian nature when coupling to structured environments. The Feynman-Vernon influence functional approach is a powerful analytical tool to study the dynamics of open quantum systems. Numerical treatments of the influence functional including the quasi-adiabatic propagator technique and the tensor-network-based time-evolving matrix product operator method have proven to be efficient in studying open quantum systems with bosonic environments. However, the numerical implementation of the fermionic path integral suffers from the Grassmann algebra involved. In this work, we present a detailed introduction to the Grassmann time-evolving matrix product operator method for fermionic open quantum systems. In particular, we introduce the concepts of Grassmann tensor, signed matrix product operator, and Grassmann matrix product state to handle the Grassmann path integral. Using the single-orbital Anderson impurity model as an example, we review the numerical benchmarks for structured fermionic environments for real-time nonequilibrium dynamics, real-time and imaginary-time equilibration dynamics, and its application as an impurity solver. These benchmarks show that our method is a robust and promising numerical approach to study strong coupling physics and non-Markovian dynamics. It can also serve as an alternative impurity solver to study strongly correlated quantum matter with dynamical mean-field theory.
In this introductory review, we present and compare four algorithms for the numerical integration of the Lindblad equation for one-dimensional quantum lattice systems. All four methods are based on Matrix Product State representations and can be viewed as extensions of the Time-Evolving Block Decimation (TEBD) algorithm to open quantum systems. Two approaches directly integrate the vectorized Lindblad equation, one of them explicitly enforcing the positivity of the density matrix. The other two rely on stochastic unravelings of the Lindblad equation, namely the quantum trajectory and quantum state diffusion approaches. We discuss the principles, numerical implementation, accuracy, and computational efficiency of the different methods, and benchmark them against an exactly solvable free fermion model.