REVIEW 3 major objections 4 minor 2 cited by
Equilibrium and nonequilibrium steady states with the repeated interaction protocol: Relaxation dynamics and energetic cost
T0 review · 3 major / 4 minor · reviewed 2026-08-10 · deepseek-v4-flash
Pith's one-line read For a qubit repeatedly colliding with thermal ancilla spins, the system settles into a diagonal nonequilibrium steady state that generally differs from the ancilla's thermal state, and the paper derives its population, relaxation rate…
desk verdict The population results are exact and worth having, but the abstract overclaims the coherence decay proof; the paper itself concedes no general proof exists. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing object is the repeated-interaction (collision) map: each step tensors the system qubit with a fresh thermal ancilla, applies the unitary $\hat U(\tau)=e^{-i\hat H_{\rm tot}\tau}$ generated by the free qubit and ancilla Hamiltonians plus the Heisenberg interaction, and traces out the ancilla. Because this Hamiltonian has a block structure separating single-excitation sectors from the zero/two-excitation sector, the population and coherence components of the qubit density matrix evolve independently. The analytic engine is the linear recursive ansatz used for the population, $p_S^{(n+1)} - p_S^{(\infty)} = \eta\,(p_S^{(n)} - p_S^{(\infty)})$, together with the analogous one-step recursion for the coherence amplitude with coefficient $\psi$ (or $\tilde\psi$ when $J_{zz}\neq0$); identifying $\eta$ from the exact one-step map and solving for the fixed point $p_S^{(\infty)}$ gives all later results, including the collision count and work formulas.
What would settle it
A concrete search: evaluate the exact coherence factor $\psi$ from Eq. (20), or $\tilde\psi$ from Eq. (29), over a fine grid or random sample of $J_{xx}, J_{yy}, J_{zz}, \omega_S, \omega_A, \tau, \beta$ and look for any point with $|\psi|\ge 1$; alternatively, iterate the exact collision map for such a candidate and check whether a nonzero coherence survives in the long-time limit. A single example with persistent coherences would falsify the diagonal-steady-state claim, and a long-time population that depends on the initial state would falsify the global-fixed-point claim.
Extended reading notes
Core claim
At the core is an exact one-step recursion for the qubit's ground-state population under the $J_{xx}$-$J_{yy}$ interaction: $p_S^{(n+1)} - p_S^{(\infty)} = \eta\,(p_S^{(n)} - p_S^{(\infty)})$, with $\eta = 1 - \frac{4(J_{xx}+J_{yy})^2}{\theta^2}\sin^2(\theta\tau/2) - \frac{4(J_{xx}-J_{yy})^2}{\phi^2}\sin^2(\phi\tau/2)$, where $\theta = \sqrt{4(J_{xx}+J_{yy})^2+(\omega_A-\omega_S)^2}$ and $\phi = \sqrt{4(J_{xx}-J_{yy})^2+(\omega_A+\omega_S)^2}$. Because $\eta$ obeys $-1<\eta<1$ except on a countable set of resonant collision times, the population converges geometrically at a rate independent of the ancilla temperature. The fixed point is $p_S^{(\infty)} = \frac{\theta^2(1-\cos\phi\tau)(J_{xx}-J_{yy})^2(1-p_A)+\phi^2(1-\cos\theta\tau)(J_{xx}+J_{yy})^2 p_A}{\phi^2(1-\cos\theta\tau)(J_{xx}+J_{yy})^2+\theta^2(1-\cos\phi\tau)(J_{xx}-J_{yy})^2}$, which equals the ancilla population $p_A$ only in the energy-conserving limit $J_{xx}=J_{yy}$. Coherences obey $c_S^{(n+1)} = \psi\,|c_S^{(n)}|$ with the complex factor $\psi$ of Eq. (20), and vanish when $|\psi|<1$, a condition proved for short collision times and verified by extensive numerical sampling; adding a $J_{zz}\hat\sigma^S_z\otimes\hat\sigma^A_z$ term leaves the steady-state population unchanged and modifies only the coherence factor. The paper also derives the trace-distance lower bound $n^* \ge \ln(\epsilon/|p_S^{(0)}-p_S^{(\infty)}|)/\ln|\eta|$ for diagonal initial states and explicit per-collision work and heat expressions.
Load-bearing premise
The load-bearing premise is that the coherence decay factor satisfies $|\psi|<1$ (and $|\tilde\psi|<1$ when $J_{zz}\neq0$) for all parameter values, so that coherences vanish and the steady state is diagonal; the paper proves this analytically only in the short-collision-time limit and otherwise relies on numerical sampling.
Editorial extensions
If this is right
- For generic parameters the steady state is initial-condition independent and nonthermal, so a repeated-interaction collider does not by itself act as a thermal bath unless the couplings are symmetric or the long-weak-collision protocol is used.
- The convergence count $n^*$ from Eq. (37) grows logarithmically with the required accuracy and, for short collisions, scales like $(J\tau)^{-2}$, so the runtime and gate count of a circuit implementation can be predicted from the interaction parameters.
- In steady state the system's energy change vanishes and the work spent switching the interaction on and off equals the heat delivered to the ancillas, meaning a nonequilibrium steady state requires continuous housekeeping work whenever $J_{xx}\neq J_{yy}$.
- Erasing coherences in these models costs no energy: the per-collision work and heat expressions depend only on the population difference, so coherence decay is thermodynamically silent.
- A protocol with equal qubit and ancilla frequencies, weak (random) couplings, and long collisions brings the qubit to the ancilla's thermal state in about five steps, achieving thermalization rather than a nonthermal steady state.
Reading between the lines
- One testable extension is to use the temperature-independence of $\eta$ as a diagnostic: measuring the population relaxation rate at several ancilla temperatures should give the same value, and any observed temperature dependence would indicate that one of the model assumptions (refreshed uncorrelated ancillas, fixed collision unitary) is being violated.
- The separation between population and coherence dynamics suggests a resource-allocation principle for collision-model quantum circuits: population engineering can be optimized using only $\eta$ and $p_S^{(\infty)}$, while coherence engineering can be optimized separately through $\psi$ (or $\tilde\psi$); the paper does not develop this two-track optimization.
- The few-long-randomized-collisions thermalization route hints that randomized coupling strengths could serve as a built-in error-averaging mechanism on noisy hardware, but that application is not explored here and would need a noise analysis.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper analyzes a repeated-interaction (collision) model in which a qubit system sequentially interacts with thermal spin ancillas under XX+YY and XX+YY+ZZ Heisenberg-type couplings. It derives an exact affine map for the ground-state population, from which it obtains the relaxation rate eta (Eq. 13) and the long-time population fixed point p_S^(infinity) (Eq. 14), showing that this fixed point generally differs from the ancilla thermal population. It further writes the coherence evolution in the form c^{(n+1)} = psi |c^{(n)}| (Eqs. 19-20 and 28-29) and claims that coherences decay to zero, yielding a diagonal nonequilibrium steady state. The paper also derives a lower bound on the number of collisions needed to reach a target accuracy for diagonal initial states (Eq. 37), presents explicit work/heat expressions (Eqs. 46-48), and identifies a long-collision weak-coupling protocol that thermalizes the qubit in a few steps.
Significance. The population-sector results are a genuine contribution: the fixed point Eq. (14) and rate Eq. (13) are derived analytically with no fitted parameters, and they are corroborated by direct simulations. The explicit work/heat decomposition and the scaling bound n* ~ (J tau)^-2 are also valuable and clearly presented. The claimed few-collision thermalization protocol is interesting. The main weakness is that the diagonal-steady-state claim rests on the unproven contraction of the coherence map |psi|<1 for all reachable phases; the paper's own summary concedes this gap. If this gap is filled, or the claims are appropriately qualified, the paper would be a solid contribution to the collision-model literature.
major comments (3)
- [II.D / III.A / Appendix B] The central claim that coherences decay to zero and the steady state is diagonal (abstract item (i); Sec. III.B) requires the coherence map c^{(n+1)} = psi(chi_n)|c^{(n)}| to be a contraction for every phase chi_n reachable under the dynamics, i.e. sup_chi |psi(chi)| < 1 for Eq. (20) and sup_chi |tilde_psi(chi)| < 1 for Eq. (29). The manuscript proves this only in special cases: short collision times with J_yy ≠ 0 (Appendix B, Eq. B3) and the symmetric case J_xx = J_yy (Eq. 25). Since chi evolves nontrivially at each step, the general case is not covered, and Sec. VI explicitly states that the authors were 'unable to provide a general proof of their decay.' This gap is load-bearing: without it, the diagonal steady state and the global-fixed-point statement in Sec. III.B are not established. The 10^8-point numerical sampling in Appendix B is useful evidence but is not an analytic proof.
- [Appendix B, Eq. (B1)] The short-time expansion in Eq. (B1) appears to be incorrect for generic parameters. Taylor-expanding Eq. (20) at fixed couplings gives psi = e^{i chi}[1 + (i tau/2)((omega_A+omega_S)/phi - (omega_A-omega_S)/theta)] + O(tau^2), not e^{i chi}(1 - i omega_S tau) as written. The coefficient stated in Eq. (B1) requires an additional unstated assumption, for example small J_xx and J_yy, that is not part of the claimed 'limit of small collision time, tau -> 0+.' Consequently the bound |psi|^2 <= 1 - 4 J_yy^2 tau^2 in Eq. (B3) is not established by the derivation as presented. This is a concrete flaw in the only analytic support for coherence decay in the general J_xx != J_yy case.
- [IV, Eqs. (30)-(31)] The few-collision thermalization protocol relies not only on the population rate eta but also on the coherence factor psi being small in the regime of weak couplings and long collision times. The argument after Eq. (30) states heuristically that 'the first term in Eq. (20) is small' and that |psi|^2 << 1, but no bound is derived for this regime. Since this protocol is presented as one of the paper's main contributions, it should be supported either by an explicit estimate for |psi| in the parameter regime specified by conditions (i)-(iii), or by the same sup-norm bound needed for the general coherence-decay claim.
minor comments (4)
- [II.C after Eq. (14)] There is a typographical error: 'Eq. (14' should read 'Eq. (14)', and shortly after, 'our ansatz yields p^{(n+1)}_S -> p^{(n)}_S' should be p^{(n+1)}_S -> p^{(n)}_S as tau -> 0^+.
- [II.D, Fig. 4 caption] The caption says 'Eqs. (19)-20)' with an extra closing parenthesis; it should read 'Eqs. (19)-(20)'.
- [V.A.1, Eq. (37)] The derivation of n* assumes logarithms of positive arguments and 0 < |eta| < 1; this is stated, but it would be helpful to note how the special cases eta = 0 and |p_S^(0) - p_S^(infinity)| <= epsilon are handled in Figs. 8-12.
- [References] Reference [30] appears to duplicate reference [28]; please consolidate or distinguish them.
Circularity Check
No circularity: the population and coherence recursions are algebraic consequences of the collision unitary, and the steady state is the fixed point of that map, not a fitted input.
full rationale
The derivation chain is self-contained. The population update in Eq. (10) is linear in p(n) for a fixed ancilla state, so the ansatz Eq. (11) is not an imported assumption but a rearrangement: η is read off as the coefficient of p(n) and pS(∞) is determined by the fixed-point condition C = pS(∞)(1 − η). Equations (13) and (14) follow algebraically, with no free parameters fit to data. The coherence evolution in Eqs. (19)–(20) is likewise obtained directly from the unitary map; the unproved bound |ψ| < 1 is a mathematical gap in the paper's claim that coherences always vanish, but it is not circularity because the steady state is not defined in terms of this bound and the numerical sampling is supporting evidence rather than a fitted parameter. The work and heat expressions in Eqs. (46)–(48) are traces of commutator-type formulas with the same unitary and satisfy the first law by construction, and the proof that ΔES(∞) = 0 follows from the explicit fixed point, not from an assumed answer. There are no load-bearing self-citations: prior work is cited for context and numerical observations, not to justify the derived formulas. Therefore the central claims reduce to exact algebra from the stated model, and no circular step is present.
Assumptions & free parameters
assumptions (5)
- domain assumption Ancillas are prepared fresh and uncorrelated before each collision, and are discarded after interacting once.
- domain assumption The ancilla is initialized in a Gibbs thermal state at inverse temperature beta.
- domain assumption The system-ancilla interaction has the Heisenberg form with fixed coupling constants Jxx, Jyy, and optionally Jzz.
- domain assumption The total Hamiltonian is time-independent during each collision, and collisions are back-to-back with no idle time.
- domain assumption The system and ancilla frequencies are positive (omega_A + omega_S > 0), so the bounds on eta hold.
Cite this review
Pith. "Pith review of Equilibrium and nonequilibrium steady states with the repeated interaction protocol: Relaxation dynamics and energetic cost." pith.science (2026). https://pith.science/paper/GKX7N65M
@misc{pith2026250105392,
author = {Pith},
title = {Pith review of: Equilibrium and nonequilibrium steady states with the repeated interaction protocol: Relaxation dynamics and energetic cost},
year = {2026},
howpublished = {\url{https://pith.science/paper/GKX7N65M}},
note = {Machine review of arXiv:2501.05392}
}
read the original 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 population and coherences of the system qubit 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.
Figures
Figures from the paper (9 more)
Forward citations
Cited by 2 Pith papers
-
Thermodynamically consistent collisional master equation in a low-density gas with internal structure
A 3D scattering-derived master equation thermalizes an internal quantum system in a thermal gas, while a two-temperature gas acts as a single non-equilibrium reservoir that can build up ergotropy.
-
The Thermodynamic Cost of Ignorance: Thermal State Preparation with One Ancilla Qubit
A single-ancilla random-interaction channel is claimed to prepare thermal states with provable simulation-time bounds, but the proof of the central remainder bound is flawed.
Reference graph
Works this paper leans on
-
[14]
Strasberg P, Schaller G, Brandes T and Esposito M 2017 Quantum and Information Thermodynamics: A Unifying Frame- work Based on Repeated Interactions Phys. Rev. X 7 021003
work page 2017
-
[1]
Alternatively, the RI model can be interpreted as involving a single ancilla, which is refreshed (initialized) after each collision with the system, effectively resetting the ancilla to its initial state before every subsequent collision. The basic RI model [5–9] has been extended in various directions to address a range of fundamental questions in quantu...
arXiv 2025
-
[2]
(32) For the models under investigation, the quantity of interest is ||ρ(n) S − ρ(∞) S ||Tr
Trace Distance The Trace Distance between two matrices, ρ and σ, within the same vector space, is defined as [66] ||ρ − σ||Tr = 1 2 Tr q (ρ − σ)(ρ − σ)† . (32) For the models under investigation, the quantity of interest is ||ρ(n) S − ρ(∞) S ||Tr. Observing that (ρ(n) S − ρ(∞) S )(ρ(n) S − ρ(∞) S )† = p(n) S − p(∞) S c(n) S c(n)∗ S p(∞) S − p(n) S !2 (33)...
-
[3]
Fidelity The Fidelity F (ρ, σ) between two matrices, ρ and σ, within the same vector space, is defined according to [66] F (ρ, σ) = Tr q√ρσ√ρ 2 . (38) If the state of the system after n steps is diagonal, namely, ρ(n) S = diag(p(n) S , 1 − p(n) S ), and the steady state of the system, ρ(∞) S is diagonal, then the definition of fidelity, given in Eq. (38) ...
-
[4]
Jxx − Jyy model with diagonalρ(0) S In Fig. 8, we assess the cost of reaching steady state assuming that the initial condition of the system is diagonal: the effective initial inverse temperature of the system is βS = 1.866, corresponding to the ground state population p(0) S = 0.866. Varying the temperature of the ancilla, either a cooling or heating pro...
-
[5]
Quantum Software Consortium: Exploring Distributed Quantum Solutions for Canada
Jxx − Jyy model with nondiagonalρ(0) S We present here the energetic analysis for a system initialized in a state with coherences and evolving under the RI protocol up to its nonequilibrium steady state. In particular, we prepare the initial state with c(0) S = 0.459−0.152i and p(0) S = 0.627. The relaxation process, along with its associated energetic co...
-
[6]
Breuer H P and Petruccione F 2002 The Theory of Open Quantum Systems (Oxford: Oxford University Press)
work page 2002
-
[7]
Ciccarello F, Lorenzo S, Giovannetti V and Palma G M 2022 Quantum collision models: Open system dynamics from repeated interactions Physics Reports 954 1–70
work page 2022
Show all 71 references
-
[8]
Cattaneo M, Giorgi G L, Zambrini R and Maniscalco S 2022 A Brief Journey through Collision Models for Multipartite Open Quantum Dynamics Open Syst. Inf. Dyn. 29 2250015
2022
-
[9]
Campbell S and Vacchini B 2021 Collision models in open system dynamics: A versatile tool for deeper insights? EPL 133 60001
2021
-
[10]
Ziman M, ˇStelmachoviˇ c P, Buˇ zek V, Hillery M, Scarani V and Gisin N 2002 Diluting quantum information: An analysis of information transfer in system-reservoir interactions Phys. Rev. A 65 042105
2002
-
[11]
Scarani V, Ziman M, ˇStelmachoviˇ c P, Gisin N and Buˇ zek V 2002 Thermalizing Quantum Machines: Dissipation and Entanglement Phys. Rev. Lett. 88 097905
2002
-
[12]
Ziman M, ˇStelmachoviˇ c P and Buˇ zek V 2005 Description of Quantum Dynamics of Open Systems Based on Collision-Like Models Open Syst. Inf. Dyn. 12 81–91
2005
-
[13]
Bruneau L, Joye A and Merkli M 2014 Repeated interactions in open quantum systems Journal of Mathematical Physics 55 075204
2014
-
[15]
Strasberg P 2019 Repeated Interactions and Quantum Stochastic Thermodynamics at Strong Coupling Phys. Rev. Lett. 123 180604
2019
-
[16]
Seah S, Nimmrichter S and Scarani V 2019 Nonequilibrium dynamics with finite-time repeated interactions Phys. Rev. E 99 042103 24 FIG. 12: Cost-accuracy tradeoff in the RI relaxation process of the Jxx − Jyy − Jzz model using a non diagonal random initial state the same used f...
2019
-
[17]
O’Connor E, Vacchini B and Campbell S 2021 Stochastic Collisional Quantum Thermometry Entropy 23 1634
2021
-
[18]
Filippov S N, Semin G N and Pechen A N 2020 Quantum master equations for a system interacting with a quantum gas in the low-density limit and for the semiclassical collision model Phys. Rev. A 101 012114
2020
-
[19]
Guarnieri G, Morrone D, C ¸ akmak B, Plastina F and Campbell S 2020 Non-equilibrium steady-states of memoryless quantum collision models Physics Letters A 384 126576
2020
-
[20]
Rodrigues F L S, De Chiara G, Paternostro M and Landi G T 2019 Thermodynamics of Weakly Coherent Collisional Models Phys. Rev. Lett. 123 140601
2019
-
[21]
Comar N E and Landi G T 2021 Correlations breaking homogenization Phys. Rev. A 104 032217
2021
-
[22]
Filippov S 2022 Multipartite Correlations in Quantum Collision Models Entropy 24 508
2022
-
[23]
Grimmer D, Layden D, Mann R B and Mart ´ ın-Mart ´ ınez E 2016 Open dynamics under rapid repeated interaction Phys. Rev. A 94 032126
2016
-
[24]
Lorenzo S, Ciccarello F and Palma G M 2017 Composite quantum collision models Phys. Rev. A 96 032107
2017
-
[25]
Guarnieri G, Kol´ aˇ r M and Filip R 2018 Steady-State Coherences by Composite System-Bath Interactions Phys. Rev. Lett. 121 070401
2018
-
[26]
Cattaneo M, De Chiara G, Maniscalco S, Zambrini R and Giorgi G L 2021 Collision Models Can Efficiently Simulate Any Multipartite Markovian Quantum Dynamics Phys. Rev. Lett. 126 130403
2021
-
[27]
Layden D, Mart ´ ın-Mart ´ ınez E and Kempf A 2016 Universal scheme for indirect quantum control Phys. Rev. A 93 040301
2016
-
[28]
24(12) 1808
Shafir U and Kosloff R 2022 Wave Function Realization of a Thermal Collision Model Entropy. 24(12) 1808
2022
-
[29]
Fiusa G and Landi G T 2024 Queued quantum collision models arXiv:2403.19408
2024 arXiv
-
[30]
Myers N M, Sable H and Scarola V W 2024 Unifying Collisional Models and the Monte Carlo Metropolis Method: Algorithms for Dynamics of Open Quantum Systems arXiv:2403.00197
2024 arXiv
-
[31]
Lacroix T, Cilluffo D, Huelga S F and Plenio M B 2024 Making Quantum Collision Models Exact arXiv:2411.13166
2024 arXiv
-
[32]
Barra F 2015 The thermodynamic cost of driving quantum systems by their boundaries Sci Rep 5 14873 25
2015
-
[33]
Wang Y, Man Z-X, Zhang Y-J and Xia Y-J 2022 Work costs and operating regimes for different manners of system-reservoir interactions via collision model New J. Phys. 24 053030
2022
-
[34]
Tabanera J, Luque I, Jacob S L, Esposito M, Barra F and Parrondo J M R 2022 Quantum collisional thermostats New J. Phys. 24 023018
2022
-
[35]
Wang Y, Man Z-X, Zhang Y-J and Xia Y-J 2022 Work costs and operating regimes for different manners of system-reservoir interactions via collision model New Journal of Physics, 24 053030
2022
-
[36]
Haack G and Joye A 2024 Entropy Production of Quantum Reset Models arXiv:2401.10022
2024
-
[37]
Jacob S L, Esposito M, Parrondo J M R and Barra F 2022 Quantum scattering as a work source Quantum 6, 750
2022
-
[38]
Tabanera-Bravo J, Parrondo J M R, Esposito M and Barra F 2023 Thermalization and Dephasing in Collisional Reservoirs Phys. Rev. Lett. 130 200402
2023
-
[39]
Jacob S L, Esposito M, Parrondo J M R and Barra F 2021 Thermalization Induced by Quantum Scattering PRX Quantum 2 020312
2021
-
[40]
Arısoy O, Campbell S and M¨ ustecaplıo˘ glu¨O E 2019 Thermalization of Finite Many-Body Systems by a Collision Model Entropy 21 1182
2019
-
[41]
Karevski D and Platini T 2009 Quantum Nonequilibrium Steady States Induced by Repeated Interactions Phys. Rev. Lett. 102 207207
2009
-
[42]
Rom´ an-Ancheyta R, Kol´ aˇ r M, Guarnieri G and Filip R 2021 Enhanced steady-state coherence via repeated system-bath interactions Phys. Rev. A 104 062209
2021
-
[43]
Cusumano S and De Chiara G 2024 Structured quantum collision models: generating coherence with thermal resources New J. Phys. 26 023001
2024
-
[44]
Seah S, Nimmrichter S, Grimmer D, Santos J P, Scarani V and Landi G T 2019 Collisional Quantum Thermometry Phys. Rev. Lett. 123 180602
2019
-
[45]
Shu A, Seah S and Scarani V 2020 Surpassing the thermal Cram´ er-Rao bound with collisional thermometry Phys. Rev. A 102 042417
2020
-
[46]
Seah S, Perarnau-Llobet M, Haack G, Brunner N and Nimmrichter S 2021 Quantum Speed-Up in Collisional Battery Charging Phys. Rev. Lett. 127 100601
2021
-
[47]
Barra B 2022 Efficiency Fluctuations in a Quantum Battery Charged by a Repeated Interaction Process Entropy 24(6) 820
2022
-
[48]
De Chiara G and Antezza M 2020 Quantum machines powered by correlated baths Phys. Rev. Research 2 033315
2020
-
[49]
Hewgill A, Gonz´ alez J O, Palao J P, Alonso D, Ferraro A and De Chiara G 2020 Three-qubit refrigerator with two-body interactions Phys. Rev. E 101 012109
2020
-
[50]
Purkayastha A, Guarnieri G, Campbell S, Prior J and Goold J 2022 Periodically refreshed quantum thermal machines Quantum 6 801
2022
-
[51]
Bettmann L P, Kewming M J and Goold J 2023 Thermodynamics of a continuously monitored double-quantum-dot heat engine in the repeated interactions framework Phys. Rev. E 107 044102
2023
-
[52]
Erbanni R, Xu X, Demarie T F and Poletti D 2023 Simulating quantum transport via collisional models on a digital quantum computer Phys. Rev. A 108 032619
2023
-
[53]
Di Bartolomeo G, Vischi M, Feri T, Bassi A and Donadi S 2024 Efficient quantum algorithm to simulate open systems through a single environmental qubit Phys. Rev. Res. 6 043321
2024
-
[54]
Yosifov A, Iyer A Ebler D and Vedral V 2024 Quantum homogenization as a quantum steady-state protocol on noisy intermediate-scale quantum hardware Phys. Rev. A 109, 032624
2024
-
[55]
Pocrnic M, Segal D and Wiebe N 2023 Quantum Simulation of Lindbladian Dynamics via Repeated Interactions
2023
-
[56]
Wootters W K and Zurek W H 1982 A single quantum cannot be cloned Nature 299 802–3
1982
-
[57]
Wootters W K and Zurek W H 2009 The no-cloning theorem Physics Today 62 76–7
2009
-
[58]
Purkayastha A, Guarnieri G, Campbell S, Prior J and Goold J 2021 Periodically refreshed baths to simulate open quantum many-body dynamics Phys. Rev. B 104 045417
2021
-
[59]
W´ ojtowicz G, Purkayastha A, Zwolak M and Rams M M 2023 Accumulative reservoir construction: Bridging continuously relaxed and periodically refreshed extended reservoirs Phys. Rev. B 107 035150
2023
-
[60]
Beever A, Violaris M, Marletto C and Vedral V 2024 Comparing coherent and incoherent models for quantum homoge- nization Phys. Rev. A 110 012464
2024
-
[61]
Giovannetti V and Palma G M 2012 Master Equations for Correlated Quantum Channels Phys. Rev. Lett. 108 040401
2012
-
[62]
Ciccarello F, Palma G M and Giovannetti V 2013 Collision-model-based approach to non-Markovian quantum dynamics Phys. Rev. A 87 040103
2013
-
[63]
Lorenzo S, Ciccarello F and Palma G M 2016 Class of exact memory-kernel master equations Phys. Rev. A 93 052111
2016
-
[64]
Kretschmer S, Luoma K and Strunz W T 2016 Collision model for non-Markovian quantum dynamics Phys. Rev. A 94 012106
2016
-
[65]
Saha T, Das A and Ghosh S 2024 Quantum homogenization in non-Markovian collisional model, New J. Phys. 26 023011
2024
-
[66]
Fan H, Roychowdhury V and Szkopek T 2005 Optimal two-qubit quantum circuits using exchange interactions Phys. Rev. A 72 052323
2005
-
[67]
Auff` eves A 2022 Quantum Technologies Need a Quantum Energy Initiative PRX Quantum 3 020101
2022
-
[68]
thermalization
We refer to the process of the system reaching a diagonal state with populations equal to that of the ancilla’s, pS = pA, as “thermalization”. In this steady state solution, the system achieves the same diagonal state as the bath’s, and, as we show in Eq. (46), there is no hou...
-
[69]
We perform the partial trace using the function Partial Trace (PT) implemented in Mathematica by Roger Colbeck[65]. 26
-
[70]
https://github.com/rogercolbeck/QI
-
[71]
Nielsen M A and Chuang I L 2010 Quantum Computation and Quantum Information: 10th Anniversary Edition (Cam- bridge University Press)
2010
Reviewed August 10, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.