Transfer matrix approach to quantum systems subject to certain Lindblad evolution

The Lindblad equation is commonly used for studying quantum dynamics in open systems that cannot be completely isolated from an environment, relevant to a broad variety of research fields, such as atomic physics, materials science, quantum biology and quantum information and computing. For electrons in condensed matter systems, the Lindblad dynamics is intractable even if their mutual Coulomb repulsion could somehow be switched off. This is because they would still be able to...

Path integrals have, over the years, proven to be an extremely versatile tool for simulating the dynamics of open quantum systems. The initial limitations of applicability of these methods in terms of the size of the system has steadily been overcome through various developments, making numerical explorations of large systems a more-or-less regular feature. However, these simulations necessitate a detailed description of the system-environment interaction through accurate spe...

In this paper, we generalize the Quantum Kinetic Monte Carlo (QKMC) method of the Schrodinger equation, which was first proposed by [Z. Cai and J. Lu. SIAM J. Sci. Comput., 40(3):B706-B722, 2018] to the Lindblad equation. This algorithm makes full use of the tensor product structure of the matrices in the Lindblad equation, thus significantly reducing the storage cost, and can calculate a more extensive system than the existing methods. We demonstrate the method in the framew...

We investigate the $t$-$W$ scheme for the anti-ferromagnetic XXX spin chain under both periodic and open boundary conditions. We propose a new parametrization of the eigenvalues of transfer matrix. Based on it, we obtain the exact solution of the system. By analyzing the distribution of zero roots at the ground state, we obtain the explicit expressions of the eigenfunctions of the transfer matrix and the associated $\mathbb{W}$ operator (see (2.8) and (3.20)) in the thermodyn...

The operator growth hypothesis (OGH) is a technical conjecture about the behaviour of operators -- specifically, the asymptotic growth of their Lanczos coefficients -- under repeated action by a Liouvillian. It is expected to hold for a sufficiently generic closed many-body system. When it holds, it yields bounds on the high frequency behavior of local correlation functions and measures of chaos (like OTOCs). It also gives a route to numerically estimating response functions....

A linearized tensor renormalization group (LTRG) algorithm is proposed to calculate the thermodynamic properties of one-dimensional quantum lattice models, that is incorporated with the infinite time-evolving block decimation technique, and allows for treating directly the two-dimensional transfer-matrix tensor network. To illustrate its feasibility, the thermodynamic quantities of the quantum XY spin chain are calculated accurately by the LTRG, and the precision is shown to ...

The time evolution of Markovian open quantum systems is governed by Lindblad master equations, whose solution can be formally written as the Lindbladian exponential acting on the initial density matrix. By expanding this Lindbladian exponential into the Taylor series, we propose a generic method for integrating Lindblad master equations. In this method, the series is truncated, retaining a finite number of terms, and the iterative actions of Lindbladian on the density matrix ...

We propose a path-integral variant of the DMRG method to calculate real-time correlation functions at arbitrary finite temperatures. To illustrate the method we study the longitudinal autocorrelation function of the $XXZ$-chain. By comparison with exact results at the free fermion point we show that our method yields accurate results up to a limiting time which is determined by the spectrum of the reduced density matrix.

The evolution of open systems, subject to both Hamiltonian and dissipative forces, is studied by writing the $nm$ element of the time ($t$) dependent density matrix in the form \ber \rho_{nm}(t)&=& \frac {1}{A} \sum_{\alpha=1}^A \gamma ^{\alpha}_n (t)\gamma^{\alpha *}_m (t) \enr The so called "square root factors", the $\gamma(t)$'s, are non-square matrices and are averaged over $A$ systems ($\alpha$) of the ensemble. This square-root description is exact. Evolution equations...

We present a new variational method, based on the matrix product operator (MPO) ansatz, for finding the steady state of dissipative quantum chains governed by master equations of the Lindblad form. Instead of requiring an accurate representation of the system evolution until the stationary state is attained, the algorithm directly targets the final state, thus allowing for a faster convergence when the steady state is a MPO with small bond dimension. Our numerical simulations...