Nonthermal entanglement dynamics in a dipole-facilitated glassy model with disconnected subspaces

We construct a dipole-facilitated kinetic constraint to partition the Hilbert space into three disconnected subspaces, two of which are nonthermal and the other acts as an intrinsic thermal bath. The resulting glassy system freely oscillates in nonthermal subspaces, making the quantum entanglement perform like a substantial qubit. The spatially spreading entanglement, quantified by concurrence, fidelity and 2-R\'{e}nyi entropy, is found to be spontaneously recovered which is absent in other reference models. Under low-frequency random flip noise, this reversible hydrodynamics of entanglement holds high fidelity and volume law, while at high frequency thermalization unusually occurs leading to a strange phase transition. Our work offers an elaborate space structure for realizing ergodicity breaking and controllable entanglement dynamics.

In traditional thermodynamic perspective, the entropy simply increases when more states in the Hilbert space are populated.Recent progresses in random quantum circuits however indicate the measurement can induce breakdown of entropic volume law, as a unitary makes the spreading entanglement entropy increase by one in each step, implying its particle nature [1,2].It is then pretty interesting whether the entanglement can be manipulated as a substantial particle, but just like in superconducing circuits [3,4], the uniform levels of particle hinder efficient quantum control and it has to be reduced to two-level system enabled by specific anharmonicity.Subsequently, it serves as the main object of the present work to introduce well-designed kinetic constraints to disconnect the subspaces to which the entanglement may spread.
In glassy systems, kinetic constraints give rise to weak ergodic breaking in presence of the fragmentation of Hilbert space with outlying nonthermal states [5][6][7].For example, in a triangular antiferromagnetic lattice, the dynamics of a local spin is frustrated due to the constraint from neighboring spins [8,9].In the spin glass phase, therefore, long-range orders are absent and the featured timescale of thermalization becomes ultraslow.This interesting local dynamics is then described by kinetically constrained models (KCMs), such as the Fredrickson-Andersen (FA) model and the East model [10][11][12].Based on the disorder-free quantum East model, it has been proven that a large number of nonthermal states can be constructed to manifest area-law entanglement entropy [13,14].In this Letter, we introduce a new dipole-facilitated kinetic constraint, which results in unusual ergodicity breaking by two disconnected subspaces.Initially from the Bell state and the Greenberger-Horne-Zeilinger (GHZ) state, we will study the exotic dynamics of entanglement within this elaborately constructed block structure of Hilbert space.
Our basic idea is to constrain the spin flip by some specific neighboring state, which is motivated by the CNOT gate in quantum computation with the control circuit being instead composed of two spins or a dipole.We start from cutting a triangular spin lattice into a quasi-one-dimensional chain as sketched in Fig. 1(a) [15,16].Odd and even sites are separated into lower and upper sides labeled by different colors.The constraint rule is described as follows.First, two neighboring spins favor antiparallel configuration due to antiferromagnetic interaction.Second, a third spin in the triangle formed with these two spins is allowed to flip only if its neighboring spin is in down state.Of course we can parallel consider the constraint from up state but that is surely equivalent.Subsequently, this newly-proposed KCM can be named as dipole-facilitated model (DFM).The model Hamiltonian is written as where X i is the x-Pauli operator on i-th site, Q i = | ↑ ↑ | i and P i = | ↓ ↓ | i are projectors of spin-up and spin-down, respectively.Throughout this work, the periodic boundary condition (PBC) is adopted for this model.On the potential experimental realization of the model, we first notice that the Rydberg atoms, an ideal experimental platform for KCMs, have been used to study nonequilibrium and slow dynamics [17][18][19].The excited and ground state of atom can be mapped as spin up and down.Under some controllable experimental conditions, the nearestneighbour Rydberg blockade effect can be well described by the so-called PXP model [20,21], which has got very similar constraint with the present DFM.In this context, we think that the DFM should not be difficult to be realized in a Rydberg atom array.Let us first discuss the Hilbert space of DFM in block representation [22].We consider four sites as instance.Two successive sites are grouped into one "group-site", and we denote (↓↓), (↑↓), (↓↑) and (↑↑) as •, ⊲, ⊳ and •, respectively.Under the action of H df , we notice that these configurations (••), (⊳⊲), (⊳•), (•⊲) and (••) are annihilated and the rest can be categorized into three subspaces.The transformation rules are These three reaction paths can be straightforwardly extended to any longer chains and subsequently construct a large thermal subspace labeled by T and two nonthermal subspaces labeled by L (all even sites are spin-down) and R (all odd sites are spin-down), as sketched in Fig. 1(b).Most interestingly, these three disconnected subspaces totally decide the active spatial state of systems.For example, if the initial state of a system without perturbation is | ↑↓ . . .↑↓ , even sites will persistently stay in spin-down state, implying that the time evolution of the system is confined to the subspace L. The reverse holds true as well, namely states initiated from | ↓↑ . . .↓↑ stays in the R subspace.Now, if we have an initial state |0 = | ↓↓ . . .↓↓ , it will evolve into completely different subspace and stay there depending on how we add spin-up into the lattice.Adding a single spin-up to even (odd) site activates R (L) subspace, respectively, and two successive spin-up's break the confinement of above two nonthermal subspaces leading to the T subspace which may act as an intrinsic thermal bath as discussed below.
One may ask if these disconnected subspaces possess similar feature with that in the Hilbert space fragmentation.Breaking the Hilbert space into disconnected sectors is the nontrivial feature of generic KCM, but in most cases the number of Krylov subspaces is exponentially dependent of the system size [19].In our DFM the number of subspaces is fixed no matter how large the system is, which is as stated essential to produce nontrivial entanglement dynamics.
In the following, we use time-evolving block decimation (TEBD) to compute the dynamics on the chain [23,24].We first show the time evolution of expectation value of z-Pauli operator Z i with different initial states for L = 24.Fig. 2(a) and (b) show the results from the initial states | ↓↑ . . .↓↑ and | ↓ . . .↓↑↓ . . .↓ , respectively.In two cases, the spins are continuously flipped on the even sites, but all odd sites keep spin-down without any flipping.That is, the system freely oscillates in the R subspace resulting from ergodicity breaking.The parallel case in L subspace is not shown.For a comparison, we calculate the initial state with two successive spin-up sites in the middle of the chain, namely | ↓ . . .↓↑↑↓ . . .↓ as shown in Fig. 2(c).This will transform the subspace L and R to T and the system becomes ergodic, as all sites are being flipped in the time evolution.
Let us now categorize all sites into two for convenience with regard to quantum computations.The sites always being spin-down in a product state are called idle sites and the others as work sites.As long as the system is sufficiently large, the probability that all work sites simultaneously flip to spin-down is extremely low.Therefore, we can denote idle sites to be 0, work sites together as 1, to form a large and substantial logical qubit.In superconducting circuits, quantum non-demolition parity measurements can be realized by assigning some qubits as control nodes to detect errors of adjacent data qubits [25].Analogously, herein, if we introduce a spin-up at a left end work site, the right end work site will tell us the initial odevity, and errors caused by noise can be detected by the appearance of successive spin-up's.A single logical qubit can thus be in a superposition state of product states in the L and R subspace, labeled by |L or |R respectively, paving a novel way for fault tolerance with this "big" qubit.Essentially, this logical qubit can be promisingly applied to various entangled states.By a two-qubit gate, it is easy to obtain a local Bell state in many-body system (L = 8), i.e., |Bell = 1 √ 2 (| ↓↓↓↓↑↓↓↓ + | ↓↓↓↑↓↓↓↓ ), which will generally be evolving into |L + |R .Then, we focus on the entanglement between middle two sites to see the influence of the disconnected subspaces.To this end, we calculate the concurrence to quantize the bipartite entanglement between them [26,27], which is defined as where λ k 's are the eigenvalues in descending order of the Hermitian matrix r = √ ρρ √ ρ and ρ = (σ y ⊗σ y )ρ * (σ y ⊗σ y ).σ y and ρ are y-Pauli matrix and reduced density matrix of the two middle sites by partially tracing others [28].
As displayed in Fig. 3(a), we observe a surprising periodicity of the entanglement, especially its maximum, which is very like a Newton's cradle.This novel nonthermal dynamic behavior stems from that two superposed bases of |Bell are restricted in two disconnected nontrivial subspaces, and the time evolution in these two subspaces are completely symmetric allowing the entanglement to be spontaneously recovered after spreading.To be comparisons, three other typical models of spin chain are also calculated, namely the East model, the PXP model and the antiferromagnetic Heisenberg model as shown in Fig. 3(c), (d) and (e).There are no conserved quantities in the PXP model which was first introduced in the quantum many-body scars and also in the classical East model for describing spin glasses [14].The Heisenberg model merely conserves the total spin.Remarkably, concurrences in these three models do not perform any visible cradle-like periodicity but just irregular oscillations.
To see how the nonthermal dynamics is broken down, we then mimic the inevitable noise which can be described as (off-diagonal) bit flip and the (diagonal) projection, respectively.In the first case if an idle site is flipped from spindown to spin-up by off-diagonal noise, one of subspaces L or R will be transferred into T.As a result, the entangled state is destroyed irreversibly by the T thermal bath.We then add a random measurement into the Hamiltonian as a noisy source, following the form: where j is a random site to be measured and the noisy term occurs every t = nT X .In Fig. 3(b), we observe that although the concurrence still oscillates, the envelope decays gradually under the noise T X = 1 and will not recover, manifesting the irreversible breaking of periodicity.It implies the measurement has induced a transition of the dynamic feature.In the case of diagonal noise, we can also set a noisy Hamiltonian similar to Eq. ( 4) with projectors Q i or P i .Compared with the off-diagonal noise which connects subspaces, these diagonal projectors do not hybridize the subspaces unless all work sites are simultaneously set to be 0 which is almost impossible, so their influence is much easier to be eliminated and is not considered here.In order to further comprehend the measurement-induced phase transition in DFM, we move on to consider the 8-site GHZ state |GHZ = 1 √ 2 (| ↑↓↑↓↑↓↑↓ + | ↓↑↓↑↓↑↓↑ ) as the initial state which also belongs to |L + |R and has many-body entangled properties distinguished from the Bell state.We calculate the fidelity F = | GHZ|φ(t) | between evolving and GHZ state to characterize the influence of wrong actions under noise, which actually quantifies the fluctuation of initial state.As shown in Fig. 4(a), for the system without random actions, the fidelity is able to come back to 1 periodically, implying perfect many-body revivals.By adding a noisy term, the fidelity is gradually decreasing over time, i.e. the distance between evolving state and |GHZ is increasing.If setting the time axis in the unit of nanosecond as for usual case, it is estimated that even in a very noisy environment with frequency being about 125 MHz, the fidelity of |GHZ can keep above 0.8 after 400 ns, suggesting a fairly good robustness [29].
We can also calculate the entanglement entropy between even and odd sites.The reduced 2-Rényi entropy with a normalized logarithmic base for the reduced density matrix ρ of even sites is then defined as which also characterizes the localization of wave function [30].As depicted in Fig. 4(b), without noise the entanglement remains a fixed value S min = log 16 2 = 0.25, with respect to two nonzero diagonal elements in the reduced density matrix.In other cases, the Rényi entropy keeps growing before saturated.There is a significant phase transition: For T X > 1 the entropy is in the early stage linearly dependent of time suggesting the system stays in an entangling phase with volume law, and for T X < 1 the entropy rapidly saturates indicating the measurement induces an area-law disentangling phase and thermalization takes place instead of localization in common cases.More interestingly, the area-law entropy is even larger than the volume-law entropy.In the small nonthermal subspaces the entanglement is periodic and localized so the increase is extremely slow, while when the localization is broken down, a large number of states will instantly evolve into large thermal T subspace to make the entropy rapidly increase, which is very similar with the entanglement tsunami [31].Consequently, the T subspace itself performs an intrinsic thermal bath for nonthermal subspaces which is the most appealing structure of DFM.Before ending, we make a simple analogy between our model and the disorder-free localization model which is composed of spinless fermions and spin-1/2 on the bond [32].As a fermionic matter field, the fermion tunneling is related to coupled spins.We take one triangle on the lattice as instance.It is easy to derive that where σ + i = | ↑ ↓ | and σ − i = | ↓ ↑ | are spin upper and lower operators.One can that, the last two terms are nothing but the constrained hopping term in the disorder-free localization model, namely the hopping is determined by spin on bond.As a result, our model H df just decomposes the localization into two subspaces and thus enables the nonthermal entanglement.
In summary, originated from triangular frustrated kinetic constraints, the DFM manifests exotic disconnected subspace structure.The entanglement exhibits appealing nonthermal periodic dynamics acting as a substantial qubit.Under low-frequency flip noise, the entanglement ballistically increases while at high frequency it saturates in an instant, indicating a phase transition to thermalization takes place stemming from the intrinsic thermal subspace.

FIG. 1 :
FIG. 1: (a) Schematic of an 8-site GHZ state with the circle representing that site is in spin-down state and the ball being spin-up.Red and blue zone denote the state of relevant sites are changeable.(b) A thermal subspace (T) and two nonthermal subspaces (L and R) are figured out in disconnected block representation.The upper and lower halves in each square (solid: ↑↑, white: ↓↓, gradient: ↑↓ or ↓↑) represent even and odd sites, respectively.

FIG. 2 :
FIG. 2: Time evolution of Zi on a 24-site lattice.(a) The initial state is | ↓↑ . . .↓↑ .The odd sites keep in spin-down and even sites behave periodic oscillation.(b) The initial state is | ↓ . . .↓↑↓ . . .↓ , i.e. a single spin-up is located at i = 12, which spreads to the ends of the chain but odd sites keep in spin-down persistently.(c) The initial state is | ↓ . . .↓↑↑↓ . . .↓ , i.e. the i = 12 and 13 sites are in spin-up.Differently, spin-up states will cover all sites after sufficiently long time.

FIG. 3 :
FIG. 3: Concurrence between middle two sites (i = 4 and 5) in 8-site systems.All initial states are set as the Bell state.(a) Time evolution with DFM shows recoverable periodic oscillations.(b) Adding a random flip noise with TX = 1, the bipartite entanglement gradually decreases during oscillation.This calculation result is averaged over 1000 samples for random X operators.Results of time evolution without noise under (c) the East model at Rokhsar-Kivelson point, (d) the PXP model from Rydberg blockade and (e) the antiferromagnetic Heisenberg model do not show periodic behavior.The East model is in open boundary condition and others are in PBC.

FIG. 4 :
FIG. 4: Fidelity and reduced 2-Rényi entropy evolving with HX.(a) The fidelity F = | GHZ|φ(t) | with TX = 0.5 (red), TX = 1 (blue), TX = 2 (orange), TX = 4 (green), TX = 8 (cyan) and TX = ∞ (black) is obtained after quench from 8-site |GHZ .The system behaves high-frequency oscillations and the envelope lines of higher bound are bolded.The TX = ∞ case is equal to evolving under H df , which shows periodicity without decay.(b) The reduced 2-Rényi entropy S(ρ) of even sites for different TX.At TX = ∞, the entropy is constant over time.Dashed lines indicate the linear dependency of entropy and time, namely the volume law.Above calculation results are all averaged over 500 runs for X operators on random sites.