Abstract We study the dynamics of a qubit system interacting with thermalized bath-ancilla spins via a repeated interaction scheme. Considering generic initial conditions for the system and employing a Heisenberg-type interaction between the system and the ancillas, we analytically prove the following: (i) the populations and coherences of the qubit system evolve independently toward a nonequilibrium steady-state solution, which is diagonal in the qubit’s energy eigenbasis. The population relaxes to this state geometrically, whereas the coherences decay through a more compound behavior. (ii) In the long time limit, the system approaches a steady state that generally differs from the thermal state of the ancilla. We derive this steady-state solution and show its dependence on the interaction parameters and collision frequency. (iii) We bound the number of interaction steps required to achieve the steady state within a specified error tolerance, and we evaluate the energetic cost associated with the process. Our key finding is that deterministic system-ancilla interactions do not typically result in the system thermalizing to the thermal state of the ancilla. Instead, they generate a distinct nonequilibrium steady state, which we explicitly derive. However, we also identify an operational regime that leads to thermalization with a few long and possibly randomized collisions.
Abstract We study the dynamics of a qubit system interacting with thermalized bath-ancilla spins via a repeated interaction scheme. Considering generic initial conditions for the system and employing a Heisenberg-type interaction between the system and the ancillas, we analytically prove the following: (i) the populations and coherences of the qubit system evolve independently toward a nonequilibrium steady-state solution, which is diagonal in the qubit’s energy eigenbasis. The population relaxes to this state geometrically, whereas the coherences decay through a more compound behavior. (ii) In the long time limit, the system approaches a steady state that generally differs from the thermal state of the ancilla. We derive this steady-state solution and show its dependence on the interaction parameters and collision frequency. (iii) We bound the number of interaction steps required to achieve the steady state within a specified error tolerance, and we evaluate the energetic cost associated with the process. Our key finding is that deterministic system-ancilla interactions do not typically result in the system thermalizing to the thermal state of the ancilla. Instead, they generate a distinct nonequilibrium steady state, which we explicitly derive. However, we also identify an operational regime that leads to thermalization with a few long and possibly randomized collisions.
We investigate the relaxation dynamics and steady state of a qubit interacting with a coherent environment through an energy-conserving flip-flop interaction within the Markovian repeated-interaction framework. We derive an exact analytical solution for the evolving and stationary state, nonperturbative in the environmental coherence. We show that environmental coherence provides a resource for dissipative state preparation by expanding the set of accessible steady states. In particular, coherent ancillas can increase the system ground- state population beyond that of the bath, corresponding to an effective temperature below the bath temperature, and enhance its steady-state purity, with these effects controlled by the interplay between the system energy scale and the system-ancilla coupling. Nevertheless, we prove that the purity of the system can never exceed the purity of the ancilla. We further analyze the thermodynamics of the protocol and show that environmental coherence can modify and even reverse the natural direction of heat transfer, giving rise to anomalous heat flow. Remarkably, for fixed interaction parameters, the direction of the coherence-induced heat transfer is governed by the relative phase between the system and ancilla coherence, whereas their magnitudes determine the amount of heat exchanged. Our results establish environmental coherence as a resource for controlling steady-state preparation and energy transfer in repeated-interaction dynamics.
We study how memory affects the thermalization of a qubit in a repeated-interaction model. In our model, a memory qubit interacts repeatedly with the system, while both qubits undergo collisions with fresh thermal ancillas. We derive exact closed evolution equations for their populations and mutual coherence, and demonstrate memory effects using a causal-break operation. Although the system converges to the target thermal state, the thermalization behavior shows non-monotonic approach. We also find a Mpemba-like reversal, in which the number of collisions required for thermalization depends non-monotonically on the temperature of the target state. Within the initial-state family considered, both phenomena require initial coherence between the system and memory or an imbalance in their excited-state populations.
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 Δ t 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.
This paper compares the cooling performance of Markovian and non-Markovian quantum collisional models for nanoscale thermodynamic devices. A minimal setup is considered: a single qubit system repeatedly collides with qubit ancillas acting as a heat bath, while a local cooling “demon” operates on the system after each collision. Two regimes are analyzed: (i) a Markovian protocol where the system interacts with a fresh thermal ancilla at every step and the ancilla is subsequently rethermalized, and (ii) a non-Markovian protocol where the system repeatedly interacts with the same initially thermal environment qubit with no reset, allowing correlations (memory) to build up. A central contribution is an efficiency definition that explicitly includes all external energetic resources, counting both the energy removed by the demon and the otherwise hidden work cost of rethermalizing the ancillas in the Markovian implementation. The demon is modeled as an amplitude damping channel; via an explicit Stinespring dilation, the extracted energy is interpreted as heat dumped into an auxiliary cold reservoir, and the paper discusses how additional demon memory/erasure (Landauer) costs would further reduce operational efficiency. Numerical simulations show that the non-Markovian protocol can cool the system to much lower final energies and achieve substantially higher efficiency than the Markovian protocol across a wide range of collision (coupling) and damping strengths. The enhancement is traced to the finite environment’s initial energy and its capacity to store correlations as usable thermodynamic resources. The results emphasize that bounds and benchmarks derived under Markovian assumptions may significantly underestimate performance in finite, memoryful environments, and that reported efficiencies should be interpreted as upper bounds when demon implementation costs are not fully modeled.
In this work we investigate a model of thermalization wherein a single ancillary qubit randomly interacts with the system to be thermalized. This not only sheds light on the emergence of Gibbs states in nature, but also provides a routine for preparing arbitrary thermal states on a digital quantum computer. For desired β and random interaction G the routine boils down to time independent Hamiltonian simulation and is represented by the channel Φ:ρ↦EGTrEnv[e−i(H+αG)t(ρ⊗Ze−βHE)ei(H+αG)t]. We rigorously prove that these dynamics reduce to a Markov chain process in the weak-coupling regime with the thermal state as the approximate fixed point. We upper bound the total simulation time required in terms of the Markov chain spectral gap λ⋆, which we compute exactly in the ground state limit. These results are independent of any eigenvalue knowledge of the system, but we are further able to show that with knowledge of eigenvalue differences λS(i)−λS(j), then the total simulation time is dramatically reduced. The ratio of the complete ignorance simulation cost to the perfect knowledge simulation cost scales as O(δmin7ε3.5λ⋆(β)3.5∥HS∥7), where δmin is related to the eigenvalue differences of the system. Additionally, we provide more specific results for single qubit and harmonic oscillator systems as well as numeric experiments with hydrogen chains. In addition to the algorithmic merits, these results can be viewed as broad extensions of the Repeated Interactions model to generic Hamiltonians with unknown interactions, giving a complete picture of the thermalization process for quantum systems.
Providing the microscopic behavior of a thermalization process has always been an intriguing issue. There are several models of thermalization, which often require interaction of the system under consideration with the microscopic constituents of the macroscopic heat bath. With the aim to simulate such a thermalization process, here we look at the thermalization of a two-level quantum system under the action of a Markovian master equation corresponding to memory-less action of a heat bath, kept at a certain temperature, using a single-qubit ancilla. A two-qubit interaction Hamiltonian ( H th , say) is then designed —with a single-qubit thermal state as the initial state of the ancilla— which gives rise to thermalization of the system qubit in the infinite time limit. Further, we study the general form of Hamiltonian, of which ours is a special case, and look for the conditions for thermalization to occur. We also derive a Lindblad-type non-Markovian master equation for the system dynamics under the general form of system-ancilla Hamiltonian.