Numerical evaluation and robustness of the quantum mean-force Gibbs state

We introduce a numerical method to determine the Hamiltonian of mean force (HMF) Gibbs state for a quantum system strongly coupled to a reservoir. The method adapts the time evolving matrix product operator (TEMPO) algorithm to imaginary-time propagation. By comparing the real-time and imaginary-time propagation for a generalized spin-boson model, we confirm that the HMF Gibbs state correctly predicts the steady state. We show that the numerical dynamics match the polaron master equation at strong coupling. We illustrate the potential of the imaginary-time TEMPO approach by exploring reservoir-induced entanglement between qubits.