The Hamiltonian of mean force is a widely used concept to describe the modification of the usual canonical Gibbs state for a quantum system whose coupling strength with the thermal bath is non-negligible. Here, we perturbatively derive general approximate expressions for the Hamiltonians of mean force in the weak-coupling approximation and in the high-temperature one. We numerically analyze the accuracy of the corresponding expressions and show that the precision of the Bloch–Redfield equantum master equation can be improved if we replace the original system Hamiltonian by the Hamiltonian of mean force.
The Hamiltonian of mean force is a widely used concept to describe the modification of the usual canonical Gibbs state for a quantum system whose coupling strength with the thermal bath is non-negligible. Here, we perturbatively derive general approximate expressions for the Hamiltonians of mean force in the weak-coupling approximation and in the high-temperature one. We numerically analyze the accuracy of the corresponding expressions and show that the precision of the Bloch–Redfield equantum master equation can be improved if we replace the original system Hamiltonian by the Hamiltonian of mean force.
When the coupling of a quantum system to its environment is nonnegligible, its steady state is known to deviate from the textbook Gibbs state. The Bloch-Redfield quantum master equation, one of the most widely adopted equations for solving the open quantum dynamics, cannot predict all the deviations of the steady state of a quantum system from the Gibbs state. In this paper, for a generic spin-boson model, we use a higher-order quantum master equation (in system-environment coupling strength) to analytically calculate all the deviations of the steady state of the quantum system up to second order in coupling strength. We also show that this steady state is exactly identical to the corresponding generalized Gibbs state, the so-called quantum mean force Gibbs state, at arbitrary temperature. All these calculations are highly general, making them immediately applicable to a wide class of systems well modeled by the spin-Boson model, spanning a diverse range of topics, from nanomaterials to various condensed-phase processes, and quantum computing (e.g., environment-induced corrections to the steady state of a superconducting qubit). As an example, we use our results to study the dynamics and the steady state of a solid-state double-quantum-dot system under physically relevant choices of parameters.
The Lindblad form guarantees complete positivity of a Markovian quantum master equation (QME). However, its microscopic derivation for a quantum system weakly interacting with a thermal bath requires several approximations, which may result in inaccuracies in the QME. Recently, various Lindbladian QMEs were derived without resorting to the secular approximation from the Redfield equation which does not guarantee the complete positivity. Here we explicitly calculate, in a perturbative manner, the equilibrium steady states of these Lindbladian QMEs. We compare the results with the steady state of the Redfield equation obtained from an analytic continuation method, which coincides with the so-called mean force Gibbs (MFG) state. The MFG state is obtained by integrating out the bath degrees of freedom for the Gibbs state of the total Hamiltonian. We explicitly show that the steady states of the Lindbladian QMEs are different from the MFG state. Our results indicate that manipulations of the Redfield equation needed to enforce complete positivity of a QME drives its steady state away from the MFG state. We also find that, in the high-temperature regime, both the steady states of the Lindbladian QMEs and MFG state reduce to the same Gibbs state of a system Hamiltonian under certain conditions.
When the strength of interaction between a quantum system and bath is non-negligible, the equilibrium state can deviate from the Gibbs state. But the expression of such a mean force Gibbs state in an arbitrary parameter regime is unknown and is numerically challenging to determine. In this work, we first review the local harmonic approximation to this problem [Maier et al., Phys. Rev. E 81, 021107 (2010)], which can accurately determine the mean force Gibbs state when either the system-bath coupling or the temperature is large, or when the third and higher derivatives of the potential are small compared to certain system-bath specific parameters. In the appropriate limit, we show that the local harmonic approximation reduces to the ultra-strong coupling and high temperature results recently derived in the literature. After deriving an estimate for the error induced by this method, we apply it to study some systems, like a quartic oscillator and a particle in a quartic double-well potential. We also apply this method to analyze the proton tunneling problem in a DNA recently studied in literature [Slocombe et al., Comm. Phys., vol. 5, no. 1, p. 109, 2022], where our results suggest the equilibrium value of the probability of mutation to be orders of magnitude lower than the steady state value obtained there (10−8 vs 10−4).
We put forth a new class of quantum master equations that correctly reproduce the asymptotic state of an open quantum system beyond the infinitesimally weak system-bath coupling limit. Our method is based on incorporating the knowledge of the reduced steady state into its dynamics. The correction not only steers the reduced system toward a correct steady state but also improves the accuracy of the dynamics, thereby refining the archetypal Born-Markov weak-coupling second-order master equations. In case of equilibrium, we use the exact mean-force Gibbs state to correct the Redfield quantum master equation. By benchmarking it with the exact solution of the damped harmonic oscillator, we show that our method also helps correct the long-standing issue of positivity violation, albeit without complete positivity. Our method of a canonically consistent quantum master equation opens a new perspective in the theory of open quantum systems leading to a reduced density matrix accurate beyond the commonly used Redfield and Lindblad equations, while retaining the same conceptual and numerical complexity.
The Gibbs state is widely taken to be the equilibrium state of a system in contact with an environment at temperature T. However, non-negligible interactions between system and environment can give rise to an altered state. Here we derive general expressions for this mean force Gibbs state, valid for any system that interacts with a bosonic reservoir. First, we derive the state in the weak coupling limit and find that, in general, it maintains coherences with respect to the bare system Hamiltonian. Second, we develop a new expansion method suited to investigate the ultrastrong coupling regime. This allows us to derive the explicit form for the mean force Gibbs state, and we find that it becomes diagonal in the basis set by the system-reservoir interaction instead of the system Hamiltonian. Several examples are discussed including a single qubit, a three-level V-system and two coupled qubits all interacting with bosonic reservoirs. The results shed light on the presence of coherences in the strong coupling regime, and provide key tools for nanoscale thermodynamics investigations.
The Gibbs state is widely taken to be the equilibrium state of a system in contact with an environment at temperature T. However, non-negligible interactions between system and environment can give rise to an altered state. Here, we derive general expressions for this mean force Gibbs state, valid for any system that interacts with a bosonic reservoir. First, we derive the state in the weak coupling limit and find that, in general, it maintains coherences with respect to the bare system Hamiltonian. Second, we develop a new expansion method suited to investigate the ultrastrong coupling regime. This allows us to derive the explicit form for the mean force Gibbs state, and we find that it becomes diagonal in the basis set by the system-reservoir interaction instead of the system Hamiltonian. Several examples are discussed including a single qubit, a three-level V-system, and two coupled qubits all interacting with bosonic reservoirs. The results shed light on the presence of coherences in the strong coupling regime, and provide key tools for nanoscale thermodynamics investigations.