Recognition: unknown
Bootstrapping Open Quantum Many-body Systems with Absorbing Phase Transitions
Pith reviewed 2026-05-10 02:43 UTC · model grok-4.3
The pith
Combining density matrix positivity with steady-state conditions produces a bootstrap hierarchy for bounding observables in open quantum many-body systems on infinite lattices.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The positivity of the density matrix together with the condition that its time derivative vanishes under the Lindblad generator yields a hierarchy of inequalities that can be systematically truncated to bound physical quantities in open quantum systems on infinite lattices that undergo absorbing phase transitions.
What carries the argument
The positivity-plus-steady-state hierarchy: the set of linear and semidefinite constraints obtained by requiring that the density matrix be positive semidefinite and that its commutator with the Liouvillian superoperator be zero.
If this is right
- Bootstrap bounds are obtained on steady-state expectation values of the quantum contact process.
- An estimate follows for the critical coupling that separates the absorbing and active phases.
- Certain ratios of expectation values become bounded in the nontrivial steady state above criticality.
- An upper bound is placed on the Liouvillian spectral gap in the subcritical phase.
Where Pith is reading between the lines
- The same hierarchy could be applied to other Lindblad models that possess absorbing states to obtain comparable bounds.
- Convergence rate of the bounds with truncation depth may indicate how rapidly the method captures infinite-lattice physics.
- Direct comparison of the bootstrap intervals with tensor-network or Monte-Carlo data would test consistency on larger but finite lattices.
Load-bearing premise
A finite truncation of the positivity and steady-state hierarchy produces bounds that remain meaningful and tighten toward the exact infinite-lattice values.
What would settle it
If successively larger truncations of the hierarchy produce bounds that fail to enclose independently computed values or fail to tighten for the quantum contact process.
read the original abstract
We demonstrate that combining the positivity of density matrices with steady-state conditions yields a systematic bootstrap method for studying open quantum many-body systems governed by Lindblad master equations on infinite lattices, which exhibit absorbing phase transitions. As a concrete example, we apply this method to the quantum contact process with an absorbing state. We obtain bootstrap bounds on steady-state expectation values, the critical coupling, certain ratios of expectation values in the nontrivial steady state in the supercritical phase, and the Liouvillian spectral gap in the subcritical phase.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript develops a bootstrap method for open quantum many-body systems governed by Lindblad master equations on infinite lattices that exhibit absorbing phase transitions. It combines positivity constraints on (reduced) density matrices with the steady-state condition obtained by setting the expectation value of the Lindblad superoperator applied to local operators to zero. The resulting semidefinite program is truncated at finite order and applied to the quantum contact process, producing bounds on steady-state expectation values, the location of the absorbing transition, certain ratios of observables in the active phase, and the Liouvillian gap in the absorbing phase.
Significance. If the truncation errors can be controlled, the approach supplies a new, parameter-light route to rigorous bounds on non-equilibrium steady states that avoids finite-size scaling artifacts. The explicit construction from positivity plus the Lindblad steady-state condition, together with the numerical demonstration on the quantum contact process, would constitute a useful addition to the toolbox for open quantum systems, especially near criticality where tensor-network or Monte-Carlo methods face severe challenges.
major comments (2)
- [Method section (hierarchy truncation)] The truncation of the positivity-plus-steady-state hierarchy at finite order N (described in the method section) is introduced without an a-priori error estimate, monotonicity argument, or numerical convergence study showing that the SDP feasible set shrinks to the true infinite-lattice steady state as N increases. This is load-bearing for the central claim that the reported bounds on the critical coupling and the Liouvillian gap are systematic and meaningful for the infinite system, particularly when the correlation length diverges at the transition.
- [Results on critical coupling and gap] The numerical bounds on the critical coupling and on the gap (presented in the results tables/figures) are obtained at a single truncation level; without evidence that these bounds are monotonic or that the discarded higher-order correlators cannot shift the location of the transition, it is unclear whether the quoted intervals are rigorous or contaminated by uncontrolled truncation bias.
minor comments (2)
- [Abstract] The abstract states that bounds are obtained but does not indicate the truncation order used for the quoted numerical values; adding this information would improve reproducibility.
- [Method section] Notation for the moment matrices or correlator hierarchy should be defined more explicitly (e.g., which operators are retained at each level) to allow readers to reproduce the SDP.
Simulated Author's Rebuttal
We thank the referee for the careful and constructive report. The concerns about the truncation scheme and its implications for the rigor of the reported bounds are well taken. We address each major comment below and will revise the manuscript to incorporate additional numerical evidence of convergence.
read point-by-point responses
-
Referee: The truncation of the positivity-plus-steady-state hierarchy at finite order N (described in the method section) is introduced without an a-priori error estimate, monotonicity argument, or numerical convergence study showing that the SDP feasible set shrinks to the true infinite-lattice steady state as N increases. This is load-bearing for the central claim that the reported bounds on the critical coupling and the Liouvillian gap are systematic and meaningful for the infinite system, particularly when the correlation length diverges at the transition.
Authors: We agree that the manuscript would benefit from a more explicit discussion of the truncation. While a general a-priori error estimate for the N to infinity limit is not derived here, each finite truncation yields rigorous bounds because the included positivity and steady-state constraints are necessary conditions satisfied by any physical steady state. In the revision we will add a dedicated subsection with numerical results at successive orders (N=3 to N=6) that demonstrate stabilization of the bounds on the critical coupling and gap, together with observed monotonic tightening of the feasible intervals. A complete convergence proof remains an open theoretical question. revision: partial
-
Referee: The numerical bounds on the critical coupling and on the gap (presented in the results tables/figures) are obtained at a single truncation level; without evidence that these bounds are monotonic or that the discarded higher-order correlators cannot shift the location of the transition, it is unclear whether the quoted intervals are rigorous or contaminated by uncontrolled truncation bias.
Authors: The original submission emphasized the highest computationally feasible truncation. The revision will include a new table and accompanying figure that compare bounds across multiple truncation orders. These data will show that the location of the critical coupling remains stable and that the intervals tighten monotonically, indicating that higher-order correlators do not shift the transition outside the reported precision. The intervals at each finite N are rigorously valid outer bounds; the multi-order comparison will quantify the truncation uncertainty. revision: yes
- A rigorous a-priori error estimate or mathematical proof of convergence of the finite-N hierarchy to the exact infinite-lattice steady state is not available and would require substantial further theoretical work.
Circularity Check
No circularity: derivation rests on standard positivity and Lindblad steady-state axioms
full rationale
The paper constructs its bootstrap by imposing positivity of finite-subsystem density matrices together with the steady-state condition obtained from the expectation value of the Lindblad superoperator. These are independent physical axioms, not defined in terms of the output bounds. The hierarchy is truncated at finite order to obtain an SDP; the resulting bounds are presented as consequences of the truncated constraints rather than as inputs that are renamed or fitted to themselves. No self-citation is invoked to justify uniqueness or to smuggle an ansatz, and no known empirical pattern is merely relabeled. The truncation is an approximation whose convergence is not proven in the excerpt, but that is a question of rigor, not circularity. The central claim therefore remains self-contained against external benchmarks.
Axiom & Free-Parameter Ledger
free parameters (1)
- truncation level of the operator hierarchy
axioms (2)
- standard math Density matrices are positive semi-definite
- domain assumption The system possesses a steady state obeying the Lindblad master equation
Reference graph
Works this paper leans on
-
[1]
W. Casteels, R. Fazio and C. Ciuti,Critical dynamical properties of a first-order dissipative phase transition,Phys. Rev. A95(2017) 012128, [1608.00717]
-
[2]
F. Minganti, A. Biella, N. Bartolo and C. Ciuti,Spectral theory of Liouvillians for dissipative phase transitions,Phys. Rev. A98(2018) 042118, [1804.11293]
-
[3]
B. Skinner, J. Ruhman and A. Nahum,Measurement-Induced Phase Transitions in the Dynamics of Entanglement,Phys. Rev. X9(2019) 031009, [1808.05953]. 14
- [4]
- [5]
-
[6]
Hinrichsen,Critical phenomena in nonequilibrium systems,Adv
H. Hinrichsen,Critical phenomena in nonequilibrium systems,Adv. Phys.49(2000) 815–958, [cond-mat/0001070]
-
[7]
M. Marcuzzi, M. Buchhold, S. Diehl and I. Lesanovsky,Absorbing State Phase Transition with Competing Quantum and Classical Fluctuations,”Phys. Rev. Lett.”116(June, 2016) 245701, [1601.07305]
-
[8]
F. Carollo, E. Gillman, H. Weimer and I. Lesanovsky,Quantum contact process,Phys. Rev. Lett.123(2019) 100604, [1902.04515]
-
[9]
E. Gillman, F. Carollo and I. Lesanovsky,Numerical Simulation of Critical Dissipative Non-Equilibrium Quantum Systems with an Absorbing State,1907.02433
- [10]
-
[11]
T. Barthel and R. H¨ ubener,Solving condensed-matter ground-state problems by semidefinite relaxations,Phys. Rev. Lett.108(2012) 200404, [1106.4966]
- [12]
- [13]
- [14]
- [15]
- [16]
- [17]
-
[18]
M. Cho, B. Gabai, H. W. Lin, J. Yeh and Z. Zheng,Bootstrapping Euclidean Two-point Correlators,2511.08560
work page internal anchor Pith review Pith/arXiv arXiv
-
[19]
G. Robichon and A. Tilloy,Bootstrapping the stationary state of bosonic open quantum systems, 2410.07384
-
[20]
L. Mortimer, D. Farina, G. Di Bello, D. Jansen, A. Leitherer, P. Mujal et al.,Certifying steady-state properties of open quantum systems,Phys. Rev. Res.7(2025) 033237, [2410.13646]
-
[21]
Cho,Nonequilibrium Phase Transitions in LargeNMatrix Quantum Mechanics, 2508.04764
M. Cho,Nonequilibrium phase transitions in Large N matrix quantum mechanics,JHEP03 (2026) 092, [2508.04764]. 15
- [22]
-
[23]
Anderson and M
P. Anderson and M. Kruczenski,Loop equation in Lattice gauge theories and bootstrap methods, EPJ Web Conf.175(2018) 11011
2018
-
[24]
V. Kazakov and Z. Zheng,Analytic and numerical bootstrap for one-matrix model and “unsolvable” two-matrix model,JHEP06(2022) 030, [2108.04830]
-
[25]
V. Kazakov and Z. Zheng,Bootstrap for lattice Yang-Mills theory,Phys. Rev. D107(2023) L051501, [2203.11360]
- [26]
- [27]
-
[28]
Cho,Bootstrapping Nonequilibrium Stochastic Processes,SciPost Phys.19(2025) 124, [2505.13609]
M. Cho,Bootstrapping Nonequilibrium Stochastic Processes,SciPost Phys.19(2025) 124, [2505.13609]
-
[29]
D. A. Mazziotti,Contracted schr¨ odinger equation: Determining quantum energies and two-particle density matrices without wave functions,Phys. Rev. A57(Jun, 1998) 4219–4234
1998
-
[30]
Nakata, H
M. Nakata, H. Nakatsuji, M. Ehara, M. Fukuda, K. Nakata and K. Fujisawa,Variational calculations of fermion second-order reduced density matrices by semidefinite programming algorithm,The Journal of Chemical Physics114(05, 2001) 8282–8292
2001
-
[31]
Bounding scalar operator dimensions in 4D CFT
R. Rattazzi, V. S. Rychkov, E. Tonni and A. Vichi,Bounding scalar operator dimensions in 4D CFT,JHEP12(2008) 031, [0807.0004]
work page Pith review arXiv 2008
-
[32]
R. Rattazzi, S. Rychkov and A. Vichi,Central Charge Bounds in 4D Conformal Field Theory, Phys.Rev.D83(2011) 046011, [1009.2725]
-
[33]
R. Rattazzi, S. Rychkov and A. Vichi,Bounds in 4D Conformal Field Theories with Global Symmetry,J.Phys.A44(2011) 035402, [1009.5985]
-
[34]
O. Hern´ andez-Lerma and J. Lasserre,Markov Chains and Invariant Probabilities. Progress in Mathematics. Birkh¨ auser Basel, 2012, https://doi.org/10.1007/978-3-0348-8024-4
-
[35]
G. Fantuzzi, D. Goluskin, D. Huang and S. I. Chernyshenko,Bounds for deterministic and stochastic dynamical systems using sum-of-squares optimization,SIAM Journal on Applied Dynamical Systems15(2016) 1962–1988, [https://doi.org/10.1137/15M1053347]
-
[36]
I. Tobasco, D. Goluskin and C. R. Doering,Optimal bounds and extremal trajectories for time averages in nonlinear dynamical systems,Physics Letters A382(Feb., 2018) 382–386, [1705.07096]
-
[37]
Korda, D
M. Korda, D. Henrion and I. Mezi´ c,Convex computation of extremal invariant measures of nonlinear dynamical systems and markov processes,Journal of Nonlinear Science31(2021) 14
2021
-
[38]
J.-B. Lasserre, D. Henrion, C. Prieur and E. Tr´ elat,Nonlinear optimal control via occupation measures and LMI-relaxations,arXiv Mathematics e-prints(Mar., 2007) math/0703377, [math/0703377]
-
[39]
Stochastic optimal control via local occupation measures,
F. Holtorf, A. Edelman and C. Rackauckas,Stochastic Optimal Control via Local Occupation Measures,arXiv e-prints(Nov., 2022) arXiv:2211.15652, [2211.15652]. 16
-
[40]
Lindblad,On the Generators of Quantum Dynamical Semigroups,Commun
G. Lindblad,On the Generators of Quantum Dynamical Semigroups,Commun. Math. Phys.48 (1976) 119
1976
-
[41]
Gorini, A
V. Gorini, A. Kossakowski and E. C. G. Sudarshan,Completely Positive Dynamical Semigroups of N Level Systems,J. Math. Phys.17(1976) 821
1976
-
[42]
Manzano,A short introduction to the Lindblad master equation,AIP Adv.10(2020) 025106, [1906.04478]
D. Manzano,A short introduction to the Lindblad master equation,AIP Adv.10(2020) 025106, [1906.04478]
-
[43]
T. E. Harris,Contact Interactions on a Lattice,The Annals of Probability2(1974) 969 – 988
1974
-
[44]
Liggett,Interacting Particle Systems
T. Liggett,Interacting Particle Systems. Classics in Mathematics. Springer Berlin, Heidelberg, 2004, https://doi.org/10.1007/978-1-4613-8542-4
- [45]
- [46]
-
[47]
M. Reehorst, S. Rychkov, D. Simmons-Duffin, B. Sirois, N. Su and B. van Rees,Navigator Function for the Conformal Bootstrap,SciPost Phys.11(2021) 072, [2104.09518]. 17
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.