Nonequilibrium steady states in multibath quantum collision models

Collision models provide a simple and versatile setting to capture the dynamics of open quantum systems. The standard approach to thermalization in this setting involves an environment of independent and identically prepared thermal qubits, interacting sequentially for a finite duration <a:math xmlns:a="http://www.w3.org/1998/Math/MathML"> <a:mrow> <a:mi mathvariant="normal">Δ</a:mi> <a:mi>t</a:mi> </a:mrow> </a:math> with the system. We compare this to a two-bath scenario in which collisional qubits are prepared in either their ground or excited states and the environment temperature is encoded in system-environment couplings. The system reaches the same thermal steady state for both settings, although even in this limit they describe fundamentally different physical processes, with the two-bath setup yielding a nonequilibrium state with finite heat currents. Non-Markovian dynamics arises when intraenvironment interactions in either setting are introduced. Here the system in the single-bath setup again reaches a steady state at the canonical temperature of the bath, but the nonequilibrium steady state of the two-bath setup tends to a different temperature due to the generation of strong system-environment and intraenvironment correlations. The two-bath setting is particularly suited to studying quantum trajectories, which are well defined also for the non-Markovian case. We showcase this with a trajectory analysis of the heat currents within a two-point measurement scheme. Finally, we consider how our results are impacted when the system-environment interaction leads to strict homogenization. Our results provide insights into the dynamics and thermodynamics of thermalization towards nonequilibrium steady states and the role of non-Markovian interactions.