Non-equilibrium absorbing state phase transitions in discrete-time quantum cellular automaton dynamics on spin lattices

We introduce a discrete-time quantum dynamics on a two-dimensional lattice that describes the evolution of a 1 + 1-dimensional spin system. The underlying quantum map is constructed such that the reduced state at each time step is separable. We show that for long times this state becomes stationary and displays a continuous phase transition in the density of excited spins. This phenomenon can be understood through a connection to the so-called Domany–Kinzel automaton, which implements a classical non-equilibrium process that features a transition to an absorbing state. Near the transition density–density correlations become long-ranged, and interestingly the same is the case for quantum correlations despite the separability of the stationary state. We quantify quantum correlations through the local quantum uncertainty and show that in some cases they may be determined experimentally solely by measuring expectation values of classical observables. This work is inspired by recent experimental progress in the realization of Rydberg lattice quantum simulators, which—in a rather natural way—permit the realization of conditional quantum gates underlying the discrete-time dynamics discussed here.


Introduction
Recent years have witnessed breakthroughs in the realization of quantum simulator platforms based on cold atomic systems [1][2][3][4][5]. One of the most recent generations of these quantum simulators is based on Rydberg atoms and offers freely programmable and addressable spin arrays [2][3][4]6]. When excited to (high-lying) Rydberg states atoms interact strongly, thereby offering a versatile platform for the study of quantum matter in an out of equilibrium.
Strong interactions between Rydberg atoms are moreover at the heart of implementations of quantum information processing protocols [7] where they allow the realization of conditional gates [8,9] that generate entangling operations. Digital quantum simulators [10] employ such gates-similar to the circuit-based approach to quantum computing-and represent a route towards emulating quantum dynamics with exotic interactions. The possibility of digitally simulating open and closed many-body systems with Rydberg lattice systems was theoretically explored in [11] and the capability of this platform for preparing exotic many-body systems and states was highlighted. While their experimental realization has not yet been achieved, first proofof-principle demonstrations of the feasibility of this idea were demonstrated within a trapped ion quantum simulator [12] and superconducting circuits [13,14].
Here we introduce a class of spin models with discrete-time quantum dynamics that lends itself rather naturally to the implementation on a Rydberg quantum simulator. The dynamics takes place within a 1+1dimensional lattice in which the directions can be thought of representing time and space, respectively. Propagation between time slices proceeds via the successive application of three-body gates that perform conditional unitary rotations. Despite the fact that the dynamics of the whole system is unitary, entanglement between time slices leads to the convergence of the reduced state on the final time slice. This stationary state may display a non-equilibrium phase transition and features non-classical correlations that become long-ranged in the vicinity of the transition point. We illustrate our idea using an example that is efficiently solvable in the sense that it permits the mapping onto the non-equilibrium process of a classical cellular automaton for site percolation.
Our work highlights the emergence of stationary behavior in closed quantum systems and introduces a new aspect in extending the concept of a cellular automaton into the quantum domain-for a few examples see [15][16][17][18][19][20][21]. Moreover, it connects to quantum generalizations of perceptrons in neural networks [22,23]. Our proposed setting provides a natural test bed for assessing the capabilities of current Rydberg lattice quantum simulators: it possesses non-trivial features, such as a phase transition and long-ranged quantum correlations, but yet can be efficiently solved. It can thus be used for the certification of a Rydberg simulator in a regime (two dimensions, strong interactions, long times) which is usually numerically intractable.

The setting
The two-dimensional lattice system we are considering is depicted in figure 1. Each row consists of N sites, with a spin-1/2 degree of freedom per site. The horizontal and vertical directions we consider as space and time, respectively. The dynamics starts from a state where all spins are in the state ñ | except for the first time slice (first row) which is prepared in the desired initial configuration. The evolution then proceeds by applying a sequence of elementary gates linking the time slice at time t to the time slice at time t+1.
For the case we are mainly interested in this work these elementary gates are unitary operators that act on three spins simultaneously-two consecutive ones on time slice t (control spins) and one on time slice t+1 (target spin), as shown in figure 1. These gates perform a rotation of the state of the target spin, conditioned on the presence or absence of excited spins (in state ñ | ) among the two-control spins. We consider this type of gate here because it can be rather naturally implemented in Rydberg lattice quantum simulators as is discussed further below. Formally, we can write the gate as . Note, that we dropped the time slice index t in the explicit forms of both the projectors and the unitary in order not to make the notation too contrived.
The rule (1) can be considered as imposing a kinetic constraint in the dynamics, reminiscent of facilitated models of glasses [24]. Constrained dynamics-which has been experimentally shown to take place e.g. in interacting Rydberg gases [25]-can give rise to complex evolution both in classical [26] and closed and open quantum systems [21,[27][28][29][30][31][32]. In particular, a rule akin to (1), of at least one nearest neighbour in the excited | . An initial state is prepared on the first time slice and propagated towards future times, i.e. lower rows, by a sequence of gates that connect subsequent time slices. In the example here we use three-body gates ( ) G m t which can be implemented for example in Rydberg lattice quantum simulators, where the spin degree of freedoms are encoded in two electronic levels. state required to allow for local evolution, is known in classical facilitated models to lead to an effective dynamics of the reaction-diffusion kind [33,34].
In our model the propagation from time slice t−1 to t is achieved via the concatenation of gates, , where we assume periodic boundary conditions when applying ( ) G N t . Note, that due to the specific choice made in equation (2) the actual order of the gates is not important since the projectors commute. The successive application of the gate  ( ) t to subsequent time slices propagates the initial state and creates a pure state (provided that the initial state has been pure) on the entire lattice.
The reduced state ρ t on time slice t is linked to the reduced state of the previous time slice by a recurrence relation: . The state ρ t is separable as it is formed by a convex superposition of product states of the form r r The weight of each state is given by the expectation value of the product of projection operators taken in the state of the previous time slice, rt 1 . In our protocol local quantum operations, such as ( ) U m t , are conditioned by a measurement result i m , that can be communicated 'classically'. Such scheme cannot produce an entangled state on time slice t. Nevertheless, ρ t can exhibit non-classical correlations as we show later.

Mean field approximation
For a first understanding of the discrete-time dynamics we conduct a mean field study. To this end we consider the evolution of the local density on site m under the gate (1), which yields . We take the expectation value over the t−1-time slice, make use of the form (2) of the projector + P m m , 1 and perform the mean field approximation (decoupling of pair correlators and assumption of homogeneity). This yields a recurrence relation, connecting the mean field densities ν at time slices t and t−1: To make progress we turn the recurrence relation into a differential equation ]. Choosing the initial condition n(0)=1, we obtain the solution which has an interesting limiting behavior at long times: for x x crit the excitation density assumes the non-zero stationary value Thus, x crit defines a critical rotation angle a p = 4 crit which in the limit  ¥ t separates two qualitatively different states. At x=x crit we find n = + ( ) ( ) t t 2 2 and thus the density displays an algebraic approach to stationarity. This result is reminiscent of mean field calculations of classical reaction-diffusion problems that feature absorbing state phase transitions [35].

Mapping to a classical non-equilibrium process
Further insight into this phase transition behavior is obtained by exploiting a link to a classical stochastic process: due to the separability of the reduced density matrices ρ t and the structure of the projectors (2), the probabilities , which appear in the reduced state (3), can be generated via a classical discrete time dynamics. This process takes place on a two-dimensional lattice, as depicted in figure 1, that contains classical spins (either up or down), initially prepared in the state   ⟫. The discrete time evolution proceeds via the classical maps Note, that instead of taking the trace, expectation values in this classical description are calculated by applying the desired operator to the probability vector and multiplying from the left with a (flat) reference state: for N spins this is The structural resemblance between the reduced state (3) and the probability vector (7) is evident. The local quantum states ρ ( k) and classical states  ⟫ s k are constructed such that they yield the same expectation values for classical observables, e.g.
x Tr m m m ,1 . Thus, also the states (3) and (7) yield identical expectation values of classical observables, and in this sense the discrete time quantum dynamics is mapped onto a classical process.
The connection to the DK cellular automaton provides an explanation for the phase transition behavior observed in the mean field calculation: it is known that the cellular automaton dynamics (6) leads to a nonequilibrium stationary state which displays a continuous (absorbing state) phase transition between a so-called inactive phase-in which the expectation value of the average density á ñ = å á ñ n n N m N m 1 is zero-and an active phase in which á ñ ¹ n 1. This transition occurs at x≈0.7 and is in the directed percolation universality class. The corresponding numerical data is shown in figures 2(a)-(c).

Quantum correlations
Despite being separable and related to a classical dynamics, the state (3) possesses non-classical correlations, as we show now. Furthermore, by exploiting the mapping to the DK cellular automaton dynamics we find that it is possible to extract quantum correlations from the measurement of classical observables, which are straightforwardly accessible on Rydberg quantum simulators [3].
As a measure for quantum correlations we employ the LQU put forward in [37] which is a variant of bipartite quantum discord [38][39][40]. It quantifies how much of the fluctuations of a local measurement is due to the noncommutativity between the state and the measured local observable, which is caused only by the states coherence, not its mixedness. By minimising over the choice of the local observable, only non-local coherences -necessarily corresponding to quantum correlations-are captured. In our model, the non-local coherence in the state (3) can be thought of as being a result of an effectively classical communication of a measurement outcome on time slice t−1 and conditioned local coherent preparation, via ( ) U m t or ( ) I m t , of atoms on time slice t, see (1) and (3).
For the reduced state ρ ij of two spins the LQU is defined as . The reduced density matrix ρ ij can be obtained entirely from measuring the local density and density-density correlations between sites i and j. To see this we exploit the special structure of the reduced state (3): each term of the sum contains a product of pure states which allows to relate expectation values of off-diagonal operators to those of diagonal observables, e.g.
. Using this property, and assuming translation invariance (á ñ = á ñ = á ñ n n n i j ), one obtains Here, the second term in equation (8) vanishes and ρ ij becomes a product state without quantum correlations. This is the case away from a phase transition where correlations between two sites are decaying rapidly as a function of their distance. Near a phase transition, however, correlations are long-ranged, as is shown in figure 2(c), where we display the connected densitydensity correlation function = -á ñ C c n ij ij 2 2 . Here also finite and long-ranged quantum correlations, characterized through the LQU, emerge, as can be seen in figure 2(d).
Note, that the entries of the density matrix (8) can be determined experimentally via state tomography, i.e. the measurement of the expectation values of products of the three Pauli matrices and the identity matrix for the atoms on sites i and j. The LQU can thus be inferred directly from experimental measurements, and its scaling, see figure 2(d), requires a number of measurements which scales only linearly in the system size, as it is exactly the case when measuring classical correlations.

Implementation with Rydberg atoms
The open cellular automaton model discussed here can be implemented on Rydberg quantum simulators [3,4]. The three-body gates underlying the gate (1) are implemented by employing the blockade interaction [8] which yields conditional unitaries [41] discriminating between the cases in which at least one or none of the source atom is excited, in direct analogy with gate (1). For the experimental investigation of the non-equilibrium dynamics it is moreover not necessary to have a two-dimensional lattice. Two parallel one-dimensional arrays (or concentric rings if periodic boundary conditions are required) are sufficient for the following protocol (see also figure 3(a)): (i) the initial state is prepared on the first chain (ring) and all sites of the second chain (ring) are prepared in the state ñ | . (ii) The discrete-time propagation is performed from the first to the second chain (ring). (iii) The first chain (ring) is reset, so that all sites are in state ñ | . (iv) The process is repeated but the role of the chains (rings) is interchanged.
An interesting practical aspect of this protocol might be that it allows to investigate collective phenomenasuch as absorbing state phase transitions-that are commonly explored in the context of open quantum systems, see e.g. [42]. The discrete-time propagation presented is fully coherent, it does not involve processes such as radiative decay photon scattering. This removes typically detrimental sources of heating and thus may facilitate the observation of dynamical processes of long times. ( ) x x crit with β being the static critical exponent. Solid lines demarcate regions in which an active and inactive phase coexist.

Generalizations
Generalizations of the dynamics presented here can be achieved by extending the fundamental gate (1) to more source/target atoms and/or by introducing more conditional spin rotations. One possible extension of the gate to K source atoms and one target atom is given by Here the operators Π m (k, K ) project on the subspace containing k excitations among the K source atoms whose state conditions the state change of the mth target atom. The latter is rotated by the unitary . We anticipate two interesting cases here: (i) K=2 source sites and rotation angles are given by α 2 =π, α 1 =α and α 0 =0: the corresponding nonequilibrium process has the two absorbing states   ñ | ... and   ñ | ... . At α=π/2 the stationary state switches between these two possibilities and displays a phase transition that is in the directed compact percolation universality class [43].
(ii) α 0 =0, which ensures the presence of the absorbing state   ñ | ... : here the mean field density follows the recurrence relation k . This process features a host of absorbing state phase transitions, coexistence regions and critical lines. Moreover, a suitable choice of the rotation angles a k allows to set all terms of order smaller than K to zero which tunes the system to a multi-critical point (similar to tri-critical directed percolation [44]): . Here the mean field density displays a power-law behavior on approach to stationarity: n~-( ) ( ) t t K 1 1 . In figure 3(b) we illustrate the case K=4.

Summary and outlook
We studied examples of quantum non-equilibrium processes that can be mapped onto classical cellular automata and therefore efficiently solved. Beyond being of conceptual interesting these findings can be applied for scrutinizing current quantum simulation platforms under challenging, yet numerically tractable, conditions. An interesting subject for future investigations is the realization of non-equilibrium processes with absorbing (dark) states [45] that feature entanglement and/or phase coherence between different sites. Those can be achieved by employing projectors in the fundamental gate (9) that project for example on two-site entangled states, in conjunction with unitary operations acting on two and more target sites. Future studies may also investigate the role of imperfections in the fundamental gates and how they influence the statics and dynamics of the absorbing state phase transition discussed here.