We numerically analyze a discretized version of the quantum contact process, a non-equilibrium model having a competition between a coherent infection process and a dissipative recovery process. Previously, the continuous-time version has been discussed to have a phase transition in a novel quantum universality class, while experimental investigations of a related process have pointed towards the transition being in the classical directed percolation class. We demonstrate that these differences stem from details of the discretization of the process and construct a systematic way to continuously interpolate between these results. Finally, we provide evidence that the modification of the universality class is due to destructive interference effects, showing a clear quantum-mechanical origin.
We numerically analyze a discretized version of the quantum contact process, a non-equilibrium model having a competition between a coherent infection process and a dissipative recovery process. Previously, the continuous-time version has been discussed to have a phase transition in a novel quantum universality class, while experimental investigations of a related process have pointed towards the transition being in the classical directed percolation class. We demonstrate that these differences stem from details of the discretization of the process and construct a systematic way to continuously interpolate between these results. Finally, we provide evidence that the modification of the universality class is due to destructive interference effects, showing a clear quantum-mechanical origin.
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.
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 investigate a nonequilibrium phase transition in a dissipative and coherent quantum spin system using the quantum Langevin equation and mean-field theory. Recently, the quantum contact process (QCP) was theoretically investigated using the Rydberg antiblockade effect, in particular, when the Rydberg atoms were excited in s states so that their interactions were regarded as being between the nearest neighbors. However, when the atoms are excited to d states, the dipole-dipole interactions become effective, and long-range interactions must be considered. Here we consider a quantum spin model with a long-range QCP, where the branching and coagulation processes are allowed not only for the nearest-neighbor pairs but also for long-distance pairs, coherently and incoherently. Using the semiclassical approach, we show that the mean-field phase diagram of our long-range model is similar to that of the nearest-neighbor QCP, where the continuous (discontinuous) transition is found in the weak (strong) quantum regime. However, at the tricritical point, we find a new universality class, which was neither that of the QCP at the tricritical point nor that of the classical directed percolation model with long-range interactions. Implementation of the long-range QCP using interacting cold gases is discussed.
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 contact process is a simple infection spreading model showcasing an out-of-equilibrium phase transition between a macroscopically active and an inactive phase. Such absorbing state phase transitions are often sensitive to the presence of quenched disorder. Traditionally, a phase transition in the disordered contact process is either triggered by dilution or by locally varying the infection rate. However, when both factors play an important role, a multicritical point emerges that remains poorly understood. Here, we study the multicritical contact process by large-scale Monte Carlo simulations in two and three dimensions. The multicritical behavior is found to be universal and exhibits ultra-slow, activated dynamical scaling, with exponents consistent with those predicted by the strong disorder renormalization group method. This finding indicates that the multicritical contact process belongs to the same universality class as the multicritical quantum Ising model, in both two and three dimensions.
Abstract Non-Hermitian (NH) quantum system recently have attracted a lots of attentions theoretically and experimentally. However, the results based on the single-particle picture may not apply to understand the property of NH many-body system. How the property of quantum many-body system especially the phase transition will be affected by the non-Hermiticity remains unclear. Here we study NH quantum contact process (QCP) model, whose effective Hamiltonian is derived from Lindbladian master equation. We show that there is a continuous phase transition induced by the non-Hermiticity in QCP. We also determine the critical exponents β of order parameter, γ of susceptibility and study the correlation and entanglement near phase transition point. We observe that the order parameter and susceptibility display infinitely singularity even for finite size system, since non-Hermiticity endow many-body system with different singular behavior from classical phase transition. Moreover our results show that the phase transition have no counterpart in Hermitian case and belongs to completely different universality class.