Pseudomode approaches allow for an exact and unapproximated description of a quantum system interacting arbitrarily strongly with a bath. In general, a system coupled to pseudomodes will not thermalize to the system's Gibbs state: This is to be expected when the system-bath coupling is non-perturbative, but conflicts with common thermodynamic intuition when the system-bath coupling is asymptotically weak. We explore under which circumstances pseudomode models satisfy detailed balance (and subsequently thermalize to the system Gibbs state) and how specific choices of parameters can force "weak" detailed balance that is restricted to a limited frequency range. A combination of Hermitian and non-Hermitian pseudomodes that yields a flat effective-temperature profile is also considered. The results and criteria established here are relevant for the construction of pseudomode models in contexts where thermodynamic consistency is required.
Pseudomode approaches allow for an exact and unapproximated description of a quantum system interacting arbitrarily strongly with a bath. In general, a system coupled to pseudomodes will not thermalize to the system's Gibbs state: This is to be expected when the system-bath coupling is non-perturbative, but conflicts with common thermodynamic intuition when the system-bath coupling is asymptotically weak. We explore under which circumstances pseudomode models satisfy detailed balance (and subsequently thermalize to the system Gibbs state) and how specific choices of parameters can force "weak" detailed balance that is restricted to a limited frequency range. A combination of Hermitian and non-Hermitian pseudomodes that yields a flat effective-temperature profile is also considered. The results and criteria established here are relevant for the construction of pseudomode models in contexts where thermodynamic consistency is required.
The pseudomode framework provides an exact description of the dynamics of an open quantum system coupled to a non-Markovian environment. Using this framework, the influence of the environment on the system is studied in an equivalent model, where the open system is coupled to a finite number of unphysical pseudomodes that follow a time-local master equation. Building on the insight that this master equation does not need to conserve the hermiticity of the pseudomode state, we here ask for the most general conditions on the master equation that guarantee the correct reproduction of the system's original dynamics. We demonstrate that our generalized approach decreases the number of pseudomodes that are required to model, for example, underdamped environments at finite temperature. We also provide an unraveling of the master equation into quantum jump trajectories of non-Hermitian states, which further facilitates the utilization of the pseudomode technique for numerical calculations by enabling the use of easily parallelizable Monte Carlo simulations. Finally, we show that pseudomodes, despite their unphysical nature, provide a natural picture in which physical processes, such as the creation of system-bath correlations or the exchange of heat, can be studied. Hence, our results pave the way for future investigations of the system-environment interaction leading to a better understanding of open quantum systems far from the Markovian weak-coupling limit. Published by the American Physical Society 2024
Abstract The treatment of quantum thermodynamic systems beyond weak coupling is of increasing relevance, yet extremely challenging. The evaluation of thermodynamic quantities in strong-coupling regimes requires a nonperturbative knowledge of the bath dynamics, which in turn relies on heavy numerical simulations. To tame these difficulties, considering thermal bosonic baths linearly coupled to the open system, we derive expressions for heat, work, and average system-bath interaction energy that only involve the autocorrelation function of the bath and two-time expectation values of system operators. We then exploit the pseudomode approach, which replaces the physical continuous bosonic bath with a small finite number of damped, possibly interacting, modes, to numerically evaluate these relevant thermodynamic quantities. We show in particular that this method allows for an efficient numerical evaluation of thermodynamic quantities in terms of one-time expectation values of the open system and the pseudomodes. We apply this framework to the investigation of two paradigmatic situations. In the first instance, we study the entropy production for a two-level system (TLS) coupled to an ohmic bath, simulated via interacting pseudomodes, allowing for the presence of time-dependent driving. Secondly, we consider a quantum thermal machine composed of a TLS interacting with two thermal baths at different temperatures, showing that an appropriate sinusoidal modulation of the coupling with the cold bath only is enough to obtain work extraction.
We investigate steady states of macroscopic quantum systems under dissipation not obeying the detailed balance condition. We argue that the Gibbs state at an effective temperature gives a good description of the steady state provided that the system Hamiltonian obeys the eigenstate thermalization hypothesis (ETH) and the perturbation theory in the weak system-environment coupling is valid in the thermodynamic limit. We derive a criterion to guarantee the validity of the perturbation theory, which is satisfied in the thermodynamic limit for sufficiently weak dissipation when the Liouvillian is gapped for bulk-dissipated systems, while the perturbation theory breaks down in boundary-dissipated chaotic systems due to the presence of diffusive transports. We numerically confirm these theoretical predictions. This paper suggests a connection between steady states of macroscopic open quantum systems and the ETH.
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.
Starting from a microscopic description of weak system-bath interactions, we derive from first principles a quantum master equation that does not rely on the well-known rotating wave approximation. This includes generic many-body systems, with Hamiltonians with vanishingly small energy spacings that forbid that approximation. The equation satisfies a general form of detailed balance, called KMS (Kubo-Martin-Schwinger), which ensures exact convergence to the many-body Gibbs state. Unlike the more common notion of GNS (Gelfand-Naimark-Segal) detailed balance, this notion is compatible with the absence of the rotating wave approximation. We show that the resulting Lindbladian dynamics not only reproduces the thermal equilibrium point up to a small renormalization of the system Hamiltonian, but it also approximates the true system evolution with an error that grows at most linearly in time, giving an exponential improvement upon previous estimates. This master equation has quasilocal jump operators, can be efficiently simulated on a quantum computer, and reduces to the usual Davies dynamics in the limit of a coarse-graining time much larger than the inverse of the smallest frequency difference. With it, we provide a rigorous model of many-body thermalization relevant to both open quantum systems and quantum algorithms.
Learning quantum interactions from finite-temperature many-body systems is a central task in emerging quantum platforms. Recently, the problem of learning from lattice quantum Gibbs states has found rigorous, efficient protocols. Nevertheless, exact Gibbs states, as the input premise, are in fact computationally intractable to prepare and may not faithfully represent generic finite-temperature quantum systems. In contrast, a system coupled to a heat bath can be stuck at an approximate stationary state (metastable state) long before it truly equilibrates. Here, we formulate a physically and algorithmically consistent alternative: learning from such metastable states of detailed-balanced master equations (Lindbladians) arising from system-bath interactions. We distill the algorithmic mechanism and structural condition underlying Gibbs-state learning and extend it in full to metastable states, attaining nearly optimal sample and computational complexity (in the system size and the precision). More broadly, we sharpen notions of metastability and develop a unified framework for finite-temperature learning.