We study the thermalization properties of one-dimensional open quantum systems coupled to baths at their boundary. The baths are driven to their thermal states via Lindblad operators, while the system undergoes Hamiltonian dynamics. We specifically consider multisite baths and investigate the extent to which the late-time steady state resembles a Gibbs state at some controllable temperature set by the baths. We study three models: a noninteracting fermion model accessible via free-fermion technology, and two interacting models, the XZ model and the chiral clock model, which are accessible via tensor network methods. We show that, by tuning towards the weak coupling and slow relaxation limits, one can engineer low temperatures in the bulk of the system provided the bath size is big enough. We use this capability to study energy transport in the XZ model at lower temperatures than previously reported. Our work paves the way for future studies of interacting open quantum systems at low temperatures.
We study the thermalization properties of one-dimensional open quantum systems coupled to baths at their boundary. The baths are driven to their thermal states via Lindblad operators, while the system undergoes Hamiltonian dynamics. We specifically consider multisite baths and investigate the extent to which the late-time steady state resembles a Gibbs state at some controllable temperature set by the baths. We study three models: a noninteracting fermion model accessible via free-fermion technology, and two interacting models, the XZ model and the chiral clock model, which are accessible via tensor network methods. We show that, by tuning towards the weak coupling and slow relaxation limits, one can engineer low temperatures in the bulk of the system provided the bath size is big enough. We use this capability to study energy transport in the XZ model at lower temperatures than previously reported. Our work paves the way for future studies of interacting open quantum systems at low temperatures.
Many physical phenomena, including thermalization in open quantum systems and quantum Gibbs sampling, are modeled by Lindbladians approximating a system weakly coupled to a bath. Understanding the convergence speed of these Lindbladians to their steady states is crucial for bounding algorithmic runtimes and thermalization timescales. We study two such families of processes: one characterizing a repeated-interaction Gibbs sampling algorithm, and another modeling open many-body quantum thermalization. We prove that both converge in polynomial time for several non-commuting systems, including high-temperature local lattices, weakly interacting fermions, and 1D spin chains. These results demonstrate that simple dissipative quantum algorithms can prepare complex Gibbs states and that Lindblad dynamics accurately capture thermal relaxation. Our proofs rely on a novel technical result that extrapolates spectral gap lower bounds from quasi-local Lindbladians to the non-local generators governing these dynamics.
We numerically study two non-interacting fermion models, a quantum wire model and a Chern insulator model, governed by open system Lindblad dynamics. The physical setup consists of a unitarily evolving "bulk" coupled via its boundaries to two dissipative "leads". The open system dynamics is chosen to drive the leads to thermal equilibrium, and by choosing different temperatures and chemical potentials for the two leads we may drive the bulk into a non-equilibrium current carrying steady state. We report two main results in this context. First, we show that for an appropriate choice of dynamics of the leads, the bulk state is also driven to thermal equilibrium even though the open system dynamics does not act directly on it. Second, we show that the steady state which emerges at late time, even in the presence of currents, is lightly entangled in the sense of having small mutual information and conditional mutual information for appropriate regions. We also report some results for the rate of approach to the steady state. These results have bearing on recent attempts to formulate a numerically tractable method to compute currents in strongly interacting models; specifically, they are relevant for the problem of designing simple leads that can drive a target system into thermal equilibrium at low temperature.
We study the temperature dependence of energy diffusion in two chaotic gapped quantum spin chains, a tilted-field Ising model and an XZ model, using an open-system approach. We introduce an energy imbalance by coupling the chain to thermal baths at its boundary and study the nonequilibrium steady states of the resulting Lindblad dynamics using a matrix product operator ansatz for the density matrix. We define an effective local temperature profile by comparing local reduced density matrices in the steady state with those of a uniform thermal state. We then measure the energy current for a variety of driving temperatures and extract the temperature dependence of the energy diffusion constant. For the Ising model, we are able to study temperatures well below the energy gap and find a regime of dilute excitations where rare three-body collisions control energy diffusion. A kinetic model correctly predicts the observed exponential increase of the energy diffusion constant at low temperatures. For the XZ model, we are only able to access intermediate to high temperatures relative to the energy gap, and we show that the data are well described by an expansion around the infinite temperature limit. We also discuss the limitations of the particular driving scheme and suggest that lower temperatures can be accessed using larger baths.
We present a methodology to simulate the quantum thermodynamics of thermal machines which are built from an interacting working medium in contact with fermionic reservoirs at a fixed temperature and chemical potential.Our method works at a finite temperature, beyond linear response and weak systemreservoir coupling, and allows for nonquadratic interactions in the working medium.The method uses mesoscopic reservoirs, continuously damped toward thermal equilibrium, in order to represent continuum baths and a novel tensor-network algorithm to simulate the steady-state thermodynamics.Using the example of a quantum-dot heat engine, we demonstrate that our technique replicates the well-known Landauer-Büttiker theory for efficiency and power.We then go beyond the quadratic limit to demonstrate the capability of our method by simulating a three-site machine with nonquadratic interactions.Remarkably, we find that such interactions lead to power enhancement, without being detrimental to the efficiency.Furthermore, we demonstrate the capability of our method to tackle complex many-body systems by extracting the superdiffusive exponent for high-temperature transport in the isotropic Heisenberg model.Finally, we discuss transport in the gapless phase of the anisotropic Heisenberg model at a finite temperature and its connection to charge conjugation parity, going beyond the predictions of single-site boundary driving configurations.
We study thermalisation and the possible occurrence of the Mpemba effect in a closed quantum setting that mimics the interaction of a system with thermal reservoirs coupled only at its boundaries. Specifically, we consider a tripartite geometry in which a finite chain, initially prepared at a finite temperature, is suddenly connected on both sides to two semi-infinite chains of the same nature held at a different temperature. These outer chains act as thermal baths, while the full system evolves unitarily under a homogeneous Hamiltonian. This setup provides a simple quantum realisation of a temperature quench and closely resembles the original scenario in which the classical Mpemba effect was first observed. We focus on two paradigmatic free-fermion models, the XX chain and the transverse-field Ising chain, which respectively preserve and break the global U(1) particle-number symmetry. As a probe of relaxation, we consider the Frobenius distance between the time-evolved reduced density matrix of the central subsystem and its stationary state, which is the thermal state at the bath temperature. Exploiting the free-fermionic structure of both models, the dynamics remains Gaussian and the Frobenius distance can be expressed exactly in terms of two-point correlation functions. Combining this representation with generalised hydrodynamics, we derive analytical predictions for the Frobenius distance in the hydrodynamic limit, providing a complete characterisation of the thermalisation process. Using these results, we investigate the possible occurrence of the Mpemba effect. We find that, despite the genuine non-equilibrium dynamics displayed by the system, no Mpemba effect arises in this setting. Our analysis identifies a broad class of boundary-driven thermalisation protocols in which relaxation is fully characterised analytically and exhibits no anomalous acceleration of equilibration.
Understanding the dynamics of dissipative quantum systems, particularly beyond the weak coupling approximation, is central to various quantum applications. While numerically exact methods provide accurate solutions, they often lack the analytical insight provided by theoretical approaches. In this study, we employ the recently developed method dubbed the effective Hamiltonian theory to understand the dynamics of system-bath configurations without resorting to a perturbative description of the system-bath coupling energy. Through a combination of mapping steps and truncation, the effective Hamiltonian theory offers both analytical insights into signatures of strong couplings in open quantum systems and a straightforward path for numerical simulations. To validate the accuracy of the method, we apply it to two canonical models: a single spin immersed in a bosonic bath and two noninteracting spins in a common bath. In both cases, we study the transient regime and the steady state limit at nonzero temperature and spanning system-bath interactions from the weak to the strong regime. By comparing the results of the effective Hamiltonian theory with numerically exact simulations, we show that although the former overlooks non-Markovian features in the transient equilibration dynamics, it correctly captures non-perturbative bath-generated couplings between otherwise non-interacting spins, as observed in their synchronization dynamics and correlations. Altogether, the effective Hamiltonian theory offers a powerful approach for understanding strong coupling dynamics and thermodynamics, capturing the signatures of such interactions in both relaxation dynamics and in the steady state limit.