Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2604.19862 · SciPost Physics

Bootstrapping Open Quantum Many-body Systems with Absorbing Phase Transitions

Minjae Cho, Colin Oscar Nancarrow, Petar Tadić, Yuan Xin

quant-phcond-mat.stat-mechcond-mat.str-elhep-thSubmitted Apr 21, 2026• Announced Apr 23, 2026• 0 citations
arXivHTMLPDF

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.

Comments
17 pages

Related papers

Criticality and Phase Classification for Quadratic Open Quantum Many-Body Systems

Yikang Zhang, Thomas Barthel

We study the steady states of translation-invariant open quantum many-body systems governed by Lindblad master equations, where the Hamiltonian is quadratic in the ladder operators, and the Lindblad operators are either linear or quadratic and Hermitian. These systems are called quasifree and quadratic, respectively. We find that steady states of one-dimensional systems with finite-range interactions necessarily have exponentially decaying Green's functions. For the quasifree case without quadratic Lindblad operators, we show that fermionic systems with finite-range interactions are noncritical for any number of spatial dimensions and provide bounds on the correlation lengths. Quasifree bosonic systems can be critical in D>1 dimensions. Last, we address the question of phase transitions in quadratic systems and find that, without symmetry constraints beyond invariance under single-particle basis and particle-hole transformations, all gapped Liouvillians belong to the same phase.

arXiv:2204.05346•Physical Review Letters•Sep 12, 2022•29 citations

Critical Behavior of the Quantum Contact Process in One Dimension

Federico Carollo, Edward Gillman, Hendrik Weimer +1 more

The contact process is a paradigmatic classical stochastic system displaying critical behavior even in one dimension. It features a nonequilibrium phase transition into an absorbing state that has been widely investigated and shown to belong to the directed percolation universality class. When the same process is considered in a quantum setting, much less is known. So far, mainly semiclassical studies have been conducted and the nature of the transition in low dimensions is still a matter of debate. Also, from a numerical point of view, from which the system may look fairly simple-especially in one dimension-results are lacking. In particular, the presence of the absorbing state poses a substantial challenge, which appears to affect the reliability of algorithms targeting directly the steady state. Here we perform real-time numerical simulations of the open dynamics of the quantum contact process and shed light on the existence and on the nature of an absorbing state phase transition in one dimension. We find evidence for the transition being continuous and provide first estimates for the critical exponents. Beyond the conceptual interest, the simplicity of the quantum contact process makes it an ideal benchmark problem for scrutinizing numerical methods for open quantum nonequilibrium systems.

arXiv:1902.04515•Physical Review Letters•Sep 6, 2019•57 citations

Local Certification of Many-Body Steady States

Miguel Frías Pérez, Antonio Acín

We present a relaxation-based method to bound expectation values on the steady state of dissipative many-body quantum systems described by master equations of the Lindblad form. Instead of targeting to represent the entire state, we promote the reduced density matrices to our variables and enforce the constraints that are imposed on them by consistency with a global steady state. The resulting constraints have the form of a semidefinite program, which allows us to efficiently bound the values a given expectation value can take in the steady state. Our results show fast convergence of the bounds with the size of the reduced density matrices, giving very competitive predictions for the steady state of several one- and two-dimensional models for an arbitrary number of particles.

quant-pharXiv:2602.06487•Feb 9, 2026•1 citation

Power-law approach to steady state in open lattices of noninteracting electrons

Mariya V. Medvedyeva, Stefan K. Kehrein

We address the question of how a nonequilibrium steady state (NESS) is reached in the Lindbladian dynamics of an open quantum system. We develop an expansion of the density matrix in terms of the NESS excitations, each of which has its own (exponential) decay rate. However, when the decay rates tend to zero for many NESS excitations (the spectral gap of the Liouvillian is closed in the thermodynamic limit), the long-time dynamics of the system can exhibit a power-law behavior. This relaxation to NESS expectation values is determined by the density of states close to zero spectral gap and the value of the operator in these states. We illustrate this main idea on the example of the lattice of noninteracting fermions coupled to Markovian leads at infinite bias voltage. The current comes towards its NESS value starting from a typical initial state as ${\ensuremath{\tau}}^{\ensuremath{-}3/2}$. This behavior is universal and independent of the space dimension.

arXiv:1406.1408•Physical Review B•Nov 11, 2014•26 citations

Universal Lindblad equation for open quantum systems

Frederik Nathan, Mark Spencer Rudner

We develop a Markovian master equation in the Lindblad form that enables the efficient study of a wide range of open quantum many-body systems that would be inaccessible with existing methods. The validity of the master equation is based entirely on properties of the bath and the system-bath coupling, without any requirements on the level structure within the system itself. The master equation is derived using a Markov approximation that is distinct from that used in earlier approaches. We provide a rigorous bound for the error induced by this Markov approximation; the error is controlled by a dimensionless combination of intrinsic correlation and relaxation timescales of the bath. Our master equation is accurate on the same level of approximation as the Bloch-Redfield equation. In contrast to the Bloch-Redfield approach, our approach ensures preservation of the positivity of the density matrix. As a result, our method is robust, and can be solved efficiently using stochastic evolution of pure states (rather than density matrices). We discuss how our method can be applied to static or driven quantum many-body systems, and illustrate its power through numerical simulation of a spin chain that would be challenging to treat by existing methods.

arXiv:2004.01469•Physical review. B./Physical review. B•Sep 3, 2020•212 citations

Numerical simulation of critical dissipative non-equilibrium quantum systems with an absorbing state

Edward Gillman, Federico Carollo, Igor Lesanovsky

Abstract The simulation of out-of-equilibrium dissipative quantum many body systems is a problem of fundamental interest to a number of fields in physics, ranging from condensed matter to cosmology. For unitary systems, tensor network methods have proved successful and extending these to open systems is a natural avenue for study. In particular, an important question concerns the possibility of approximating the critical dynamics of non-equilibrium systems with tensor networks. Here, we investigate this by performing numerical simulations of a paradigmatic quantum non-equilibrium system with an absorbing state: the quantum contact process. We consider the application of matrix product states and the time-evolving block decimation algorithm to simulate the time-evolution of the quantum contact process at criticality. In the Lindblad formalism, we find that the Heisenberg picture can be used to improve the accuracy of simulations over the Schrödinger approach, which can be understood by considering the evolution of operator-space entanglement. Furthermore, we also consider a quantum trajectories approach, which we find can reproduce the expected universal behaviour of key observables for a significantly longer time than direct simulation of the average state. These improved results provide further evidence that the universality class of the quantum contact process is not directed percolation, which is the class of the classical contact process.

arXiv:1907.02433•New Journal of Physics•Sep 1, 2019•35 citations
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.
We study the steady states of translation-invariant open quantum many-body systems governed by Lindblad master equations, where the Hamiltonian is quadratic in the ladder operators, and the Lindblad operators are either linear or quadratic and Hermitian. These systems are called quasifree and quadratic, respectively. We find that steady states of one-dimensional systems with finite-range interactions necessarily have exponentially decaying Green's functions. For the quasifree case without quadratic Lindblad operators, we show that fermionic systems with finite-range interactions are noncritical for any number of spatial dimensions and provide bounds on the correlation lengths. Quasifree bosonic systems can be critical in D>1 dimensions. Last, we address the question of phase transitions in quadratic systems and find that, without symmetry constraints beyond invariance under single-particle basis and particle-hole transformations, all gapped Liouvillians belong to the same phase.
The contact process is a paradigmatic classical stochastic system displaying critical behavior even in one dimension. It features a nonequilibrium phase transition into an absorbing state that has been widely investigated and shown to belong to the directed percolation universality class. When the same process is considered in a quantum setting, much less is known. So far, mainly semiclassical studies have been conducted and the nature of the transition in low dimensions is still a matter of debate. Also, from a numerical point of view, from which the system may look fairly simple-especially in one dimension-results are lacking. In particular, the presence of the absorbing state poses a substantial challenge, which appears to affect the reliability of algorithms targeting directly the steady state. Here we perform real-time numerical simulations of the open dynamics of the quantum contact process and shed light on the existence and on the nature of an absorbing state phase transition in one dimension. We find evidence for the transition being continuous and provide first estimates for the critical exponents. Beyond the conceptual interest, the simplicity of the quantum contact process makes it an ideal benchmark problem for scrutinizing numerical methods for open quantum nonequilibrium systems.
We present a relaxation-based method to bound expectation values on the steady state of dissipative many-body quantum systems described by master equations of the Lindblad form. Instead of targeting to represent the entire state, we promote the reduced density matrices to our variables and enforce the constraints that are imposed on them by consistency with a global steady state. The resulting constraints have the form of a semidefinite program, which allows us to efficiently bound the values a given expectation value can take in the steady state. Our results show fast convergence of the bounds with the size of the reduced density matrices, giving very competitive predictions for the steady state of several one- and two-dimensional models for an arbitrary number of particles.
We address the question of how a nonequilibrium steady state (NESS) is reached in the Lindbladian dynamics of an open quantum system. We develop an expansion of the density matrix in terms of the NESS excitations, each of which has its own (exponential) decay rate. However, when the decay rates tend to zero for many NESS excitations (the spectral gap of the Liouvillian is closed in the thermodynamic limit), the long-time dynamics of the system can exhibit a power-law behavior. This relaxation to NESS expectation values is determined by the density of states close to zero spectral gap and the value of the operator in these states. We illustrate this main idea on the example of the lattice of noninteracting fermions coupled to Markovian leads at infinite bias voltage. The current comes towards its NESS value starting from a typical initial state as ${\ensuremath{\tau}}^{\ensuremath{-}3/2}$. This behavior is universal and independent of the space dimension.
We develop a Markovian master equation in the Lindblad form that enables the efficient study of a wide range of open quantum many-body systems that would be inaccessible with existing methods. The validity of the master equation is based entirely on properties of the bath and the system-bath coupling, without any requirements on the level structure within the system itself. The master equation is derived using a Markov approximation that is distinct from that used in earlier approaches. We provide a rigorous bound for the error induced by this Markov approximation; the error is controlled by a dimensionless combination of intrinsic correlation and relaxation timescales of the bath. Our master equation is accurate on the same level of approximation as the Bloch-Redfield equation. In contrast to the Bloch-Redfield approach, our approach ensures preservation of the positivity of the density matrix. As a result, our method is robust, and can be solved efficiently using stochastic evolution of pure states (rather than density matrices). We discuss how our method can be applied to static or driven quantum many-body systems, and illustrate its power through numerical simulation of a spin chain that would be challenging to treat by existing methods.
Abstract The simulation of out-of-equilibrium dissipative quantum many body systems is a problem of fundamental interest to a number of fields in physics, ranging from condensed matter to cosmology. For unitary systems, tensor network methods have proved successful and extending these to open systems is a natural avenue for study. In particular, an important question concerns the possibility of approximating the critical dynamics of non-equilibrium systems with tensor networks. Here, we investigate this by performing numerical simulations of a paradigmatic quantum non-equilibrium system with an absorbing state: the quantum contact process. We consider the application of matrix product states and the time-evolving block decimation algorithm to simulate the time-evolution of the quantum contact process at criticality. In the Lindblad formalism, we find that the Heisenberg picture can be used to improve the accuracy of simulations over the Schrödinger approach, which can be understood by considering the evolution of operator-space entanglement. Furthermore, we also consider a quantum trajectories approach, which we find can reproduce the expected universal behaviour of key observables for a significantly longer time than direct simulation of the average state. These improved results provide further evidence that the universality class of the quantum contact process is not directed percolation, which is the class of the classical contact process.