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.
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.
The contact process is a paradigmatic example of nonequilibrium dynamics, with broad applications ranging from chemistry to sociology. Its quantum counterpart, the quantum contact process (QCP), extends the classical model to include coherent processes. Despite sustained interest, the nature of the transition in the one-dimensional (1D) QCP remains debatable. Here, combining Liouvillian spectral analysis, the tensor jump method, exact quantum jump Monte Carlo, and truncated Wigner simulations, we show that 1D QCP undergoes a continuous absorbing-state phase transition, with critical exponents distinct from the classical case. We further find Liouvillian gap closes well below the critical point, highlighting that spectral gap analysis alone cannot distinguish a phase transition from metastability in the QCP. Crucially, the 1D QCP is weakly entangled, yet even this weak entanglement is indispensable for capturing the correct critical behavior, whereas semiclassical methods artificially stabilize the active state and predict a spurious first-order transition. Our work establishes the quantum origin of the phase transition in the 1D QCP and underscores the essential role of entanglement in dissipative quantum many-body systems.
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.
Phase transitions to absorbing states are among the simplest examples of critical phenomena out of equilibrium. The characteristic feature of these models is the presence of a fluctuationless configuration that the dynamics cannot leave, which has proved a rather stringent requirement in experiments. Recently, a proposal to seek such transitions in highly tunable systems of cold atomic gases offers to probe this physics and, at the same time, to investigate the robustness of these transitions to quantum coherent effects. Here, the authors specifically focus on the interplay between classical and quantum fluctuations in a simple driven open quantum model which, in the classical limit, reproduces a contact process, which is known to undergo a continuous transition in the ``directed percolation'' universality class. The authors derive an effective long-wavelength field theory for the present class of open spin systems and show that, due to quantum fluctuations, the nature of the transition changes from second to first order, passing through a critical point which appears to belong instead to the ``tricritical directed percolation'' class.
The role of quantum fluctuations in modifying the critical behavior of nonequilibrium phase transitions is a fundamental but unsolved question. In this study, we examine the absorbing state phase transition of a 1D chain of qubits undergoing a contact process that involves both coherent and classical dynamics. We adopt a discrete-time quantum model with states that can be described in the stabilizer formalism, and therefore allows for an efficient simulation of large system sizes. The extracted critical exponents indicate that the absorbing state phase transition of this Clifford circuit model belongs to the directed percolation universality class. This suggests that the inclusion of quantum fluctuations does not necessarily alter the critical behavior of nonequilibrium phase transitions of purely classical systems. Finally, we extend our analysis to a non-Clifford circuit model, where a tentative scaling analysis in small systems reveals critical exponents that are also consistent with the directed percolation universality class. Published by the American Physical Society 2024
The absorbing-state transition in the three-dimensional contact process with and without quenched randomness is investigated by means of Monte Carlo simulations. In the clean case, a reweighting technique is combined with a careful extrapolation of the data to infinite time to determine with high accuracy the critical behavior in the three-dimensional directed percolation universality class. In the presence of quenched spatial disorder, our data demonstrate that the absorbing-state transition is governed by an unconventional infinite-randomness critical point featuring activated dynamical scaling. The critical behavior of this transition does not depend on the disorder strength, i.e., it is universal. Close to the disordered critical point, the dynamics is characterized by the nonuniversal power laws typical of a Griffiths phase. We compare our findings to the results of other numerical methods, and we relate them to a general classification of phase transitions in disordered systems based on the rare region dimensionality.
Stochastic processes with absorbing states feature examples of nonequilibrium universal phenomena. While the classical regime has been thoroughly investigated in the past, relatively little is known about the behavior of these nonequilibrium systems in the presence of quantum fluctuations. Here, we theoretically address such a scenario in an open quantum spin model which, in its classical limit, undergoes a directed percolation phase transition. By mapping the problem to a nonequilibrium field theory, we show that the introduction of quantum fluctuations stemming from coherent, rather than statistical, spin flips alters the nature of the transition such that it becomes first order. In the intermediate regime, where classical and quantum dynamics compete on equal terms, we highlight the presence of a bicritical point with universal features different from the directed percolation class in a low dimension. We finally propose how this physics could be explored within gases of interacting atoms excited to Rydberg states.