Structure of quantum mean force Gibbs states for coupled harmonic systems

An open quantum system interacting with a heat bath at a given temperature is expected to reach the mean force Gibbs (MFG) state as a steady state. The MFG state is given by tracing out the bath degrees of freedom from the equilibrium Gibbs state of the total system plus bath. When the interaction between the system and the bath is not negligible, it is different from the usual system Gibbs state obtained from the system Hamiltonian only. Using the path integral method, we present the exact MFG state for a coupled system of quantum harmonic oscillators in contact with multiple thermal baths at the same temperature. We develop a nonperturbative method to calculate the covariances with respect to the MFG state. By comparing them with those obtained from the system Gibbs state, we find that the effect of coupling to the bath decays exponentially as a function of the distance from the system-bath boundary. This is similar to the skin effect found recently for a quantum spin chain interacting with an environment. Using the exact results, we also investigate the ultrastrong coupling limit where the coupling between the system and the bath gets arbitrarily large and makes a connection with the recent result found for a general quantum system.