We introduce a wide class of quantum maps that arise in collisional reservoirs and are able to thermalize a system if they operate in conjunction with an additional dephasing mechanism. These maps describe the effect of collisions and induce transitions between populations that obey detailed balance, but also create coherences that prevent the system from thermalizing. We combine these maps with a unitary evolution acting during random Poissonian times between collisions and causing dephasing. We find that, at a low collision rate, the nontrivial combination of these two effects causes thermalization in the system. This scenario is suitable for modeling collisional reservoirs at equilibrium. We justify this claim by identifying the conditions for such maps to arise within a scattering theory approach and provide a thorough characterization of the resulting thermalization process.
We introduce a wide class of quantum maps that arise in collisional reservoirs and are able to thermalize a system if they operate in conjunction with an additional dephasing mechanism. These maps describe the effect of collisions and induce transitions between populations that obey detailed balance, but also create coherences that prevent the system from thermalizing. We combine these maps with a unitary evolution acting during random Poissonian times between collisions and causing dephasing. We find that, at a low collision rate, the nontrivial combination of these two effects causes thermalization in the system. This scenario is suitable for modeling collisional reservoirs at equilibrium. We justify this claim by identifying the conditions for such maps to arise within a scattering theory approach and provide a thorough characterization of the resulting thermalization process.
Abstract Collisional reservoirs are becoming a major tool for modelling open quantum systems. In their simplest implementation, an external agent switches on, for a given time, the interaction between the system and a specimen from the reservoir. Generically, in this operation the external agent performs work onto the system, preventing thermalization when the reservoir is at equilibrium. One can recover thermalization by considering an autonomous global setup where the reservoir particles colliding with the system possess a kinetic degree of freedom. The drawback is that the corresponding scattering problem is rather involved. Here, we present a formal solution of the problem in one dimension and for flat interaction potentials. The solution is based on the transfer matrix formalism and allows one to explore the symmetries of the resulting scattering map. One of these symmetries is micro-reversibility, which is a condition for thermalization. We then introduce two approximations of the scattering map that preserve these symmetries and, consequently, thermalize the system. These relatively simple approximate solutions constitute models of quantum thermostats and are useful tools to study quantum systems in contact with thermal baths. We illustrate their accuracy in a specific example, showing that both are good approximations of the exact scattering problem even in situations far from equilibrium. Moreover, one of the models consists of the removal of certain coherences plus a very specific randomization of the interaction time. These two features allow one to identify as heat the energy transfer due to switching on and off the interaction. Our results prompt the fundamental question of how to distinguish between heat and work from the statistical properties of the exchange of energy between a system and its surroundings.
We introduce energy-space quantum walks as a minimal framework to investigate equilibration, thermalization, and irreversibility from an effective-dynamics perspective. By mapping the configuration space of a walk onto a ladder of energy eigenlevels, we reinterpret thermalization as transport in energy space, independently of microscopic system--bath details. At the classical level, the resulting birth--death--lazy dynamics leads to equilibration of the energy distribution and, under suitable conditions, to a Gibbs stationary state. We then embed this dynamics into a unitary, collision-assisted model in which coherence is controlled by a single parameter. A central result is a structural decoupling between population dynamics and coherence generation: while the populations evolve according to the classical process and relax to the Gibbs distribution, the full quantum state exhibits a persistent coherence-induced deviation from the thermal manifold. This establishes a minimal scenario of thermalization without state convergence, where equilibration occurs at the level of populations but not at the level of the full density operator. We quantify this effect using the thermal distance to the Gibbs state and derive perturbative bounds that relate the long-time deviation to classical transport properties. Our results show that coherence acts as a controllable and quantitatively bounded source of nonthermal behavior, providing a clear separation between classical equilibration and genuinely quantum corrections.
Open quantum systems that comply with the master equation and detailed balance decay in a non-oscillatory manner to thermal equilibrium. Beyond the weak coupling limit, systems that break microreversibility (e.g., in the presence of magnetic fields) violate detailed balance but still thermalize. We study the thermalization of these systems and show that a temperature increase produces novel exceptional points that indicate a sharp transition in the thermalization dynamics. A further temperature increase fuels oscillations of the energy level populations, even without quantum coherences. Moreover, the violation of detailed balance introduces an energy scale that characterizes the oscillatory regime at high temperatures.
We investigate the dynamics of a qubit chain locally coupled to a thermal reservoir, modeled through repeated collisions with particles drawn from a heat bath. Under suitable conditions, the resulting Lindblad equation is thermodynamically consistent -- it drives the system toward thermal equilibrium -- while remaining local at short times. This framework reveals a crossover between the global Lindblad equation derived from the secular approximation in weak-coupling theory and the local dissipative models often employed in the literature, which generally fail to ensure thermodynamic consistency.
We study the applicability of collisional models for non-Markovian dynamics of open quantum systems. By allowing interactions between the separate environmental degrees of freedom in between collisions we are able to construct a collision model that allows us to study quantum memory effects in open system dynamics. We also discuss the possibility to embed non-Markovian collision model dynamics into Markovian collision model dynamics in an extended state space. As a concrete example we show how, using the proposed class of collision models, we can discretely model non-Markovian amplitude damping of a qubit. In the time-continuous limit, we obtain the well-known results for spontaneous decay of a two-level system into a structured zero-temperature reservoir.
Abstract Understanding thermalization and irreversibility in open quantum systems requires char-acterizing the competition between coherent unitary dynamics and dissipative environmen-tal interactions. In this work, we present a systematic numerical study of thermalization dynamics in a two-qubit open quantum system governed by a completely positive and trace-preserving Lindblad master equation. The coherent evolution rate is controlled by a real scaling parameter α multiplying the Hamiltonian commutator, while decoherence is indepen-dently controlled by a dephasing rate γ . Although such a scaling of coherent dynamics can be absorbed into a rescaling of time in isolated systems, it becomes physically meaningful in open systems where decoherence introduces an absolute timescale. We investigate how the relative strength of coherent evolution and decoherence influences observable dynamics, purity decay, entropy production, and thermalization times. Using robust local and two-site observables, we show that varying α at fixed γ produces distinct transient trajectories, even though all evolutions converge to the same steady state. We further define an operational thermalization time τ ( γ, α ) based on a purity threshold and map the resulting thermalization-time surface across a broad parameter range. Increasing α generally accelerates oscillatory dynamics and modifies the approach to equilibrium, leading to measurable shifts in thermalization times without altering asymptotic states. These findings provide a clear and physically consistent characterization of how Hamil-tonian strength and decoherence rate jointly determine thermalization behavior in open quantum systems. The framework offers a practical parameterization for numerical mod-eling and experimental studies of coherent–dissipative competition in quantum information processing, quantum control, and open-system dynamics.