Overcoming the Lamb Shift in System-Bath Interaction Models via KMS Detailed Balance: High-Accuracy Thermalization with Time-Bounded Interactions
Pith reviewed 2026-05-19 17:04 UTC · model grok-4.3
The pith
Engineering the system-bath interaction so its transition rates satisfy KMS detailed balance makes the steady state arbitrarily close to the Gibbs state in the weak-coupling limit, no matter how the Lamb shift is structured.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
If the system-bath interaction is engineered so that the transition part of the approximate Lindbladian generator satisfies the KMS detailed balance condition, then the unique fixed point of the dynamics can be made arbitrarily close to the Gibbs state in the weak-coupling limit, regardless of the structure of the Lamb shift term. This remains true even when the approximate Lindbladian differs substantially from the ideal Davies generator and the Lamb shift term does not commute with the thermal state.
What carries the argument
The KMS detailed balance condition imposed only on the transition rates of the approximate Lindbladian, which fixes the steady state independently of the Lamb shift contribution.
If this is right
- High-accuracy thermalization becomes possible without forcing the Lamb shift to commute with the target thermal state.
- An end-to-end complexity of O(ε^{-1}) follows for preparing the Gibbs state to accuracy ε.
- The guarantees extend to any Hamiltonian whose KMS-detailed-balance Lindbladian mixes rapidly.
- The same fixed-point control holds for time-bounded interactions in the system-bath model.
Where Pith is reading between the lines
- The KMS condition may serve as a design rule for engineering interactions in experiments where realizing an exact Davies generator is impractical.
- Similar separation of transition rates from shift terms could be tested in other approximate dissipative models outside thermalization.
- Numerical checks on small qubit systems with explicit non-commuting Lamb shifts would provide direct evidence for the claimed convergence.
Load-bearing premise
The Hamiltonian must belong to the class for which the corresponding KMS-detailed-balance Lindbladian is already known to be fast mixing; the paper supplies no new mixing-time proof.
What would settle it
Take a small finite-dimensional system whose ideal Davies generator mixes rapidly but whose Lamb shift does not commute with the Gibbs state; numerically integrate the engineered weak-coupling dynamics and verify whether the distance from the steady state to the Gibbs state tends to zero as the coupling parameter approaches zero.
read the original abstract
We investigate quantum thermal state preparation algorithms based on system-bath interactions and uncover a surprising phenomenon in the weak-coupling regime. We rigorously prove that, if the system-bath interaction is engineered so that the transition part of the approximate Lindbladian generator satisfies the Kubo--Martin--Schwinger (KMS) detailed balance condition, then the unique fixed point of the dynamics can be made arbitrarily close to the Gibbs state in the weak-coupling limit, regardless of the structure of the Lamb shift term. Importantly, this remains true even when the approximate Lindbladian differs substantially from the ideal Davies generator and the Lamb shift term does not commute with the thermal state. Our result shows that the role of the KMS detailed balance condition extends well beyond standard Lindbladian dynamics, serving as a general principle for a broader class of dissipative systems. Furthermore, by combining this with a general perturbation framework, we bound the mixing time of the dynamics and establish an end-to-end complexity of $\mathcal{O}(\varepsilon^{-1})$ for Gibbs state preparation. These guarantees apply to any Hamiltonian whose associated KMS-detailed-balance Lindbladian is known to be fast mixing.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript claims that engineering the system-bath interaction so the transition part of the approximate Lindbladian satisfies the KMS detailed balance condition makes the unique fixed point of the dynamics arbitrarily close to the Gibbs state in the weak-coupling limit, even when the Lamb shift does not commute with the thermal state and the generator differs from the ideal Davies form. It further combines this fixed-point result with a general perturbation framework to bound the mixing time, yielding an end-to-end complexity of O(ε^{-1}) for Gibbs state preparation, but only for Hamiltonians whose associated ideal KMS-detailed-balance Lindbladian is already known to be fast mixing.
Significance. If the central fixed-point result holds, the work is significant for extending the applicability of the KMS condition as a design principle beyond standard Lindbladian generators to a broader class of approximate dissipative dynamics. The ability to achieve high-accuracy thermalization despite non-commuting Lamb shifts and the derivation of a concrete complexity bound (conditional on fast mixing of the ideal case) represent a useful contribution to quantum thermal state preparation algorithms with time-bounded interactions.
major comments (1)
- The O(ε^{-1}) complexity claim in the abstract and main results section is load-bearing on the external fast-mixing assumption for the ideal KMS Lindbladian; the manuscript supplies no new mixing-time analysis or proof for this ideal case and treats the property as an input, so the end-to-end guarantee does not apply when that assumption fails even if the fixed-point closeness result remains valid.
Simulated Author's Rebuttal
We thank the referee for their careful reading of the manuscript and for the positive overall assessment. We address the major comment below.
read point-by-point responses
-
Referee: The O(ε^{-1}) complexity claim in the abstract and main results section is load-bearing on the external fast-mixing assumption for the ideal KMS Lindbladian; the manuscript supplies no new mixing-time analysis or proof for this ideal case and treats the property as an input, so the end-to-end guarantee does not apply when that assumption fails even if the fixed-point closeness result remains valid.
Authors: We agree that the O(ε^{-1}) end-to-end complexity bound is conditional on the fast-mixing property of the ideal KMS-detailed-balance Lindbladian, which we treat as an external input drawn from prior literature rather than re-proving in this work. This conditional character is already stated explicitly in the abstract and main text. The primary technical contribution remains the fixed-point result, which holds independently of mixing times. To address the concern, we will add a clarifying sentence in the abstract and introduction emphasizing that the complexity guarantee applies only when the ideal generator is known to be fast mixing. revision: partial
Circularity Check
No significant circularity; result is conditional on external fast-mixing input
full rationale
The paper explicitly qualifies its end-to-end complexity claim by requiring that the ideal KMS-detailed-balance Lindbladian already be known to be fast mixing, treating this as an external input rather than deriving it internally. The fixed-point closeness result (unique fixed point arbitrarily close to Gibbs state when transition part satisfies KMS condition) is proved independently of mixing times and does not reduce to a fitted parameter, self-definition, or self-citation chain. No equations or steps in the abstract or described derivation exhibit the patterns of self-definitional reduction, fitted inputs renamed as predictions, or ansatz smuggling. The perturbation framework is applied only conditionally, preserving independence from the target result.
Axiom & Free-Parameter Ledger
axioms (2)
- domain assumption The dynamics possesses a unique fixed point in the weak-coupling limit.
- domain assumption The KMS-detailed-balance Lindbladian is fast mixing for the Hamiltonians under consideration.
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
if the system-bath interaction is engineered so that the transition part of the approximate Lindbladian generator satisfies the KMS detailed balance condition, then the unique fixed point ... can be made arbitrarily close to the Gibbs state ... regardless of the structure of the Lamb shift term
-
IndisputableMonolith/Foundation/AlphaCoordinateFixation.leancostAlphaLog_high_calibrated_iff unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
L(ρ) = −i[H_Lamb, ρ] + L_KMS(ρ) ... L_KMS satisfies the KMS detailed balance condition
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Forward citations
Cited by 1 Pith paper
-
Rigorous error bounds for dissipative thermal state preparation from weak system-bath coupling
The unitary contribution from weak system-bath coupling in collision-model thermal state preparation tightens the fixed-point error bound, scaling rigorously as J² where J is the coupling strength.
Reference graph
Works this paper leans on
-
[1]
End-to-end Complexity Analysis of End-to-End Efficient Quantum Thermal State Preparation [12] Following the discussion in Section II D, we give additional assumptions to the Setup II.6: Assumption III.3.•We set the Gaussian widthσ= Ω(β)and the coupling strengthα=o(1). •We set ˜β= 2β 2−β 2/(4σ2) , g(ω) = βp 2π(2−β 2/(4σ2)) exp − (βω+ 1) 2 2 (2−β 2/(4σ2)) ,...
-
[2]
Applications to the algorithms in [15, 21] In this section, we focus on the algorithms proposed in [15, 21], reviewed in Section II D. In both setups, it is straightforward to verify that the transition part ofDVAS ,f,T (ρ)exactly satisfiesthe KMS detailed balance condition, and therefore admits an expansion of the same form as Eq. (3). Consequently, Theo...
-
[3]
Rapid thermalization of spin chain commuting Hamiltonians.Phys
Ivan Bardet, ´Angela Capel, Li Gao, Angelo Lucia, David P´ erez-Garc´ ıa, and Cambyse Rouz´ e. Rapid thermalization of spin chain commuting Hamiltonians.Phys. Rev. Lett., 130(6):060401, 2023. 18
work page 2023
-
[4]
Thiago Bergamaschi and Chi-Fang Chen. Fast mixing of quantum spin chains at all temperatures.arXiv/2510.08533, 2026
- [5]
-
[6]
Chi-Fang Chen, Michael Kastoryano, Fernando G. S. L. Brand˜ ao, and Andr´ as Gily´ en. Efficient quantum thermal simulation. Nature, 646(8085):561–566, 10 2025
work page 2025
-
[7]
Quantum Thermal State Preparation
Chi-Fang Chen, Michael J Kastoryano, Fernando GSL Brand˜ ao, and Andr´ as Gily´ en. Quantum thermal state preparation. arXiv:2303.18224, 2023
work page internal anchor Pith review Pith/arXiv arXiv 2023
-
[8]
An efficient and exact noncommutative quantum Gibbs sampler
Chi-Fang Chen, Michael J Kastoryano, and Andr´ as Gily´ en. An efficient and exact noncommutative quantum Gibbs sampler.arXiv:2311.09207, 2023
-
[9]
Efficient quantum algorithms for simulating Lindblad evolution
Richard Cleve and Chunhao Wang. Efficient quantum algorithms for simulating Lindblad evolution. InICALP 2017, volume 80, pages 17:1–17:14, 2017
work page 2017
-
[10]
Single-ancilla ground state preparation via Lindbladians.Phys
Zhiyan Ding, Chi-Fang Chen, and Lin Lin. Single-ancilla ground state preparation via Lindbladians.Phys. Rev. Research, 6:033147, 2024
work page 2024
-
[11]
Efficient quantum Gibbs samplers with Kubo–Martin–Schwinger detailed balance condition.Commun
Zhiyan Ding, Bowen Li, and Lin Lin. Efficient quantum Gibbs samplers with Kubo–Martin–Schwinger detailed balance condition.Commun. Math. Phys., 406(3):67, 2025
work page 2025
-
[12]
Polynomial-time preparation of low-temperature Gibbs states for 2D Toric Code.arXiv:2410.01206, 2024
Zhiyan Ding, Bowen Li, Lin Lin, and Ruizhe Zhang. Polynomial-time preparation of low-temperature Gibbs states for 2D Toric Code.arXiv:2410.01206, 2024
-
[13]
Simulating open quantum systems using Hamiltonian simulations.PRX Quantum, 5:020332, 2024
Zhiyan Ding, Xiantao Li, and Lin Lin. Simulating open quantum systems using Hamiltonian simulations.PRX Quantum, 5:020332, 2024
work page 2024
-
[14]
End-to-end efficient quantum thermal and ground state preparation made simple.arXiv:2508.05703, 2025
Zhiyan Ding, Yongtao Zhan, John Preskill, and Lin Lin. End-to-end efficient quantum thermal and ground state preparation made simple.arXiv:2508.05703, 2025
-
[15]
Mixing time of open quantum systems via hypocoercivity.Phys
Di Fang, Jianfeng Lu, and Yu Tong. Mixing time of open quantum systems via hypocoercivity.Phys. Rev. Lett., 134:140405, Apr 2025
work page 2025
-
[16]
Matthew Hagan and Nathan Wiebe. The thermodynamic cost of ignorance: Thermal state preparation with one ancilla qubit.arXiv:2502.03410, 2025
- [17]
-
[18]
Quantum logarithmic Sobolev inequalities and rapid mixing.J
Michael J Kastoryano and Kristan Temme. Quantum logarithmic Sobolev inequalities and rapid mixing.J. Math. Phys., 54(5):1–34, 2013
work page 2013
-
[19]
Rapid thermalization of dissipative many-body dynamics of commuting Hamiltonians.Commun
Jan Kochanowski, Alvaro M Alhambra, Angela Capel, and Cambyse Rouz´ e. Rapid thermalization of dissipative many-body dynamics of commuting Hamiltonians.Commun. Math. Phys., 2024
work page 2024
- [20]
-
[21]
Speeding up quantum markov processes through lifting.arXiv:2505.12187, 2025
Bowen Li and Jianfeng Lu. Speeding up quantum markov processes through lifting.arXiv:2505.12187, 2025
-
[22]
Simulating Markovian open quantum systems using higher-order series expansion
Xiantao Li and Chunhao Wang. Simulating Markovian open quantum systems using higher-order series expansion. In ICALP 2023, volume 261, pages 87:1–87:20, 2023
work page 2023
- [23]
-
[24]
Quantum simulation of lindbladian dynamics via repeated interactions
Matthew Pocrnic, Dvira Segal, and Nathan Wiebe. Quantum simulation of lindbladian dynamics via repeated interactions. Journal of Physics A: Mathematical and Theoretical, 58(30):305302, jul 2025
work page 2025
-
[25]
Thermal state preparation via rounding promises.Quantum, 7:1132, 2023
Patrick Rall, Chunhao Wang, and Pawel Wocjan. Thermal state preparation via rounding promises.Quantum, 7:1132, 2023
work page 2023
-
[26]
Carlos Ramon-Escandell, Alessandro Prositto, and Dvira Segal. Thermal state preparation by repeated interactions at and beyond the Lindblad limit.arXiv:2506.12166, 2025
-
[27]
Optimal quantum algorithm for Gibbs state preparation
Cambyse Rouz´ e, Daniel Stilck Fran¸ ca, and´Alvaro M Alhambra. Optimal quantum algorithm for Gibbs state preparation. arXiv:2411.04885, 2024
- [28]
- [29]
-
[30]
Trace inequalities for matrices.Bulletin of the Australian Mathematical Society, 87(1):139–148, 2013
Khalid Shebrawi and Hussien Albadawi. Trace inequalities for matrices.Bulletin of the Australian Mathematical Society, 87(1):139–148, 2013
work page 2013
-
[31]
Oles Shtanko and Ramis Movassagh. Preparing thermal states on noiseless and noisy programmable quantum processors. arXiv:2112.14688, 2021
-
[32]
Alhambra, Daniel Stilck Fran¸ ca, and Cambyse Rouz´ e
Samuel Slezak, Matteo Scandi, ´Alvaro M. Alhambra, Daniel Stilck Fran¸ ca, and Cambyse Rouz´ e. Polynomial-time ther- malization and gibbs sampling from system-bath couplings.arXiv/2601.16154, 2026
-
[33]
Theχ 2- divergence and mixing times of quantum Markov processes.J
Kristan Temme, Michael James Kastoryano, Mary Beth Ruskai, Michael Marc Wolf, and Frank Verstraete. Theχ 2- divergence and mixing times of quantum Markov processes.J. Math. Phys., 51(12), 2010
work page 2010
-
[34]
Yu Tong and Yongtao Zhan. Fast mixing of weakly interacting fermionic systems at any temperature.PRX Quantum, 6:030301, Jul 2025
work page 2025
-
[35]
Ke Wang and Zhiyan Ding. Beyond lindblad dynamics: Rigorous guarantees for thermal and ground state preservation under system bath interactions.arXiv/2512.03457, 2025
-
[36]
Rapid quantum ground state preparation via dissipative dynamics
Yongtao Zhan, Zhiyan Ding, Jakob Huhn, Johnnie Gray, John Preskill, Garnet Kin-Lic Chan, and Lin Lin. Rapid quantum 19 ground state preparation via dissipative dynamics.arXiv/2503.15827, 2025
-
[37]
ˇStˇ ep´ anˇSm´ ıd, Richard Meister, Mario Berta, and Roberto Bondesan. Polynomial time quantum Gibbs sampling for Fermi-Hubbard model at any temperature.arXiv:2501.01412, 2025
-
[38]
Rapid Mixing of Quantum Gibbs Samplers for Weakly-Interacting Quantum Systems
ˇStˇ ep´ anˇSm´ ıd, Richard Meister, Mario Berta, and Roberto Bondesan. Rapid mixing of quantum gibbs samplers for weakly- interacting quantum systems.arXiv/2510.04954, 2025. APPENDIX. The appendix collects the technical details deferred from the main text and is organized as follows: •In Section A, we give a self-contained proof of the general integer mi...
work page internal anchor Pith review Pith/arXiv arXiv 2025
-
[39]
Auxiliary operator construction and approximation to the Gibbs state We first explain how the auxiliary operatorρ ∗ is constructed. The asymptotic expansion from the informal discus- sion Section IV suggests looking for a correction of the form ρ∗ =ρ β +α 2E. We will first define the correction termE, then use it to control∥ρ ∗ −ρ β∥1, and finally show th...
-
[40]
Auxiliary operator approximating the fixed point We next give a rigorous verification that the same construction also makesρ ∗ an approximate fixed point of the full channel. The informal asymptotic expansion already indicates that the averaged order-α 2 term should cancel, so the task here is to turn that heuristic cancellation into a quantitative estima...
-
[41]
Fixed point approximating the Gibbs state It remains to show that the fixed point of the quantum channel Φ α is close to the Gibbs stateρ β. According to Corollary C.7, we have constructed an auxiliary operatorρ ∗ that is close toρ β, and we have shown that one application of Φ α movesρ ∗ only by a higher-order term: ∥Φα(ρ∗)−ρ ∗∥1 ≤ O σβlog(σ)α 4 ,(C8) Th...
-
[42]
To do so, we begin by defining the following weighted L2-distance between quantum states
Mixing time of the ideal channel We first prove the first part (D1) of Proposition D.1. To do so, we begin by defining the following weighted L2-distance between quantum states. Definition D.2.For two quantum statesρandσ, define theρ β-weighted distance as dβ(ρ, σ) := ρ−1/4 β (ρ−σ)ρ −1/4 β 2 . 29 Note the spectral gap ofL KMS implies the contraction ofL K...
-
[43]
From the ideal channel to the implemented channel Next, we transfer the contraction of Φα stated in Lemma D.3 to the mixing time of the implemented channel Φ α, which will establish the second part of Proposition D.1. We will use the following stability argument for the mixing time of quantum channels, which is part of the [12, Theorem 8]. For completenes...
-
[44]
Proof of Theorem III.4.LetL ε := log 8∥ρ−1/2 β ∥2 ε =O(β∥H∥+ log(1/ε))
Proof of Theorem III.4 and Corollary III.5 We now combine Propositions C.1 and D.1 to prove Theorem III.4 and Corollary III.5. Proof of Theorem III.4.LetL ε := log 8∥ρ−1/2 β ∥2 ε =O(β∥H∥+ log(1/ε)). By Proposition D.1, one can choose σ= ckβ2 λgap , α 2 =c α ελgap σβlog(σ) log−1 4∥ρ−1/2 β ∥2 ε ! , T 0 ≥2σ p log((α2βlog(σ)) −1), so that tmix,Φα(2ε)≤ O 1 λga...
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.