Abstract We generalize time-evolving matrix product operators method to nonequilibrium quantum transport problems. The nonequilibrium current is obtained via numerical differentiation of the generating functional which is represented as a tensor network. The approach is numerically exact and the non-Markovian effects are fully taken into account. In the transport process, a part of the heat that flows out from a bath flows into the system and other baths, and the rest is stored in the system-bath coupling part. We take the spin-boson model as a demonstration to show the details of this heat flowing and the establishment of a steady current between two baths.
Abstract We generalize time-evolving matrix product operators method to nonequilibrium quantum transport problems. The nonequilibrium current is obtained via numerical differentiation of the generating functional which is represented as a tensor network. The approach is numerically exact and the non-Markovian effects are fully taken into account. In the transport process, a part of the heat that flows out from a bath flows into the system and other baths, and the rest is stored in the system-bath coupling part. We take the spin-boson model as a demonstration to show the details of this heat flowing and the establishment of a steady current between two baths.
Simulating the non-perturbative and non-Markovian dynamics of open quantum systems is a very challenging many body problem, due to the need to evolve both the system and its environments on an equal footing. Tensor network and matrix product states (MPS) have emerged as powerful tools for open system models, but the numerical resources required to treat finite-temperature environments grow extremely rapidly and limit their applications. In this study we use time-dependent variational evolution of MPS to explore the striking theory of Tamascelli et al. (Phys. Rev. Lett. 2019, 123, 090402.) that shows how finite-temperature open dynamics can be obtained from zero temperature, i.e., pure wave function, simulations. Using this approach, we produce a benchmark dataset for the dynamics of the Ohmic spin-boson model across a wide range of coupling strengths and temperatures, and also present a detailed analysis of the numerical costs of simulating non-equilibrium steady states, such as those emerging from the non-perturbative coupling of a qubit to baths at different temperatures. Despite ever-growing resource requirements, we find that converged non-perturbative results can be obtained, and we discuss a number of recent ideas and numerical techniques that should allow wide application of MPS to complex open quantum systems.
We present a numerically exact method to compute the full counting statistics of heat transfer in non-Markovian open quantum systems, which is based on the time-evolving matrix product operator algorithm.This approach is applied to the paradigmatic spin-boson model in order to calculate the mean and fluctuations of the heat transferred to the environment during thermal equilibration.We show that system-reservoir correlations make a significant contribution to the heat statistics at low temperature and present a variational theory that quantitatively explains our numerical results.We also demonstrate a fluctuation-dissipation relation connecting the mean and variance of the heat distribution at high temperature.Our results reveal that system-bath interactions make a significant contribution to heat transfer even when the dynamics of the open system is effectively Markovian.The method presented here provides a flexible and general tool to predict the fluctuations of heat transfer in open quantum systems in nonperturbative regimes.
We present a numerically exact method to compute the full counting statistics of heat transfer in non-Markovian open quantum systems, which is based on the time-evolving matrix product operator (TEMPO) algorithm. This approach is applied to the paradigmatic spin-boson model in order to calculate the mean and fluctuations of the heat transferred to the environment during thermal equilibration. We show that system-reservoir correlations make a significant contribution to the heat statistics at low temperature and present a variational theory that quantitatively explains our numerical results. We also demonstrate a fluctuation-dissipation relation connecting the mean and variance of the heat distribution at high temperature. Our results reveal that system-bath interactions make a significant contribution to heat dissipation even when the dynamics of the open system is effectively Markovian. The method presented here provides a flexible and general tool to predict the fluctuations of heat transfer in open quantum systems in non-perturbative regimes.
Open quantum systems play an important role in developing quantum sciences, and therefore the study of corresponding numerical method is of great significance. For the open quantum systems, the quasi-adiabatic propagator path integral invented in 1990s is one of the few numerically exact methods. However, its computational complexity scales exponentially with system size and correlation length, and therefore its application is limited in practical calculation. In recent years, the study and application of tensor network have made rapid progress. Representing the path integral by tensor network makes the computational complexity increase polynomially, thus greatly improving the computational efficiency. Such a new method is called time-evolving matrix product operator. At the very beginning, the reduced density matrix is represented as a matrix product state. Then the time evolution of the system can be achieved by iteratively applying matrix product operators to the matrix product state. The iterative process is amenable to the standard matrix product states compression algorithm, which keeps the computational cost on a polynomial scale. The time-evolving matrix product operator is an efficient, numerically exact and fully non-Markovian method, which has a broad application prospect in the study of quantum open systems. For instance, it is already used in the study of the thermalization, heat statistic, heat transfer and optimal control of the quantum open systems, and conversely it can be also used to investigate the effect of the system on the environment. In addition, the TEMPO method is naturally related to the process tensor, and can be used to calculate the correlation function of the system efficiently. In this article we review this method and its applications. We give a brief introduction of the path integral formalism of Caldeira-Leggett model. According to the path integral formalism, we demonstrate the usage of quasi-adiabatic propagator path integral method. we give the basic idea of matrix product states, and we show how to recast quasi-adiabatic propagator path integral method into time-evolving matrix product operators method by employing the concept of matrix product states and matrix product operators, and give a review of its applications. In addition, we use the calculation results of physical quantities, correlation functions and heat currents in the spin-boson model to illustrate the applications of the time-evolving matrix product operator method.
In order to model realistic quantum devices it is necessary to simulate quantum systems strongly coupled to their environment. To date, most understanding of open quantum systems is restricted either to weak system-bath couplings or to special cases where specific numerical techniques become effective. Here we present a general and yet exact numerical approach that efficiently describes the time evolution of a quantum system coupled to a non-Markovian harmonic environment. Our method relies on expressing the system state and its propagator as a matrix product state and operator, respectively, and using a singular value decomposition to compress the description of the state as time evolves. We demonstrate the power and flexibility of our approach by numerically identifying the localisation transition of the Ohmic spin-boson model, and considering a model with widely separated environmental timescales arising for a pair of spins embedded in a common environment.
Dynamics of open quantum systems with structured reservoirs can often be simulated efficiently with tensor network influence functionals. The standard variants of the time-evolving matrix product operator (TEMPO) method are applicable when the systems is coupled to Gaussian bosonic baths via hermitian coupling operators that mutually commute. In this work we introduce a generalization to cases where the system is coupled to a single reservoir through multiple non-commuting operators, representing the most general form of linear system-bath coupling. We construct a Gaussian influence functional that properly handles Trotter errors arising from a finite evolution time step, thus ensuring convergence for long evolution times. Based on this result, the uniform TEMPO scheme can be employed to obtain a matrix product operator form of the influence functional, enabling efficient simulations of the real-time dynamics of the open system. As a demonstration, we simulate the time evolution of driven two-level emitters coupled to a bosonic lattice at different lattice sites.