
2
connected momentum sector is a two-level system. For
nonintegrable models, e.g., the long-range interacting
TFIC or two-dimensional quantum Ising models, the di-
rect connection between DQPTs and OP zeros occurs
for large quenches, where the transverse-field strength is
much larger than the coupling constant, and which can be
considered perturbatively close to classical precession in
two-level systems during the short timescales where this
behavior is prominent [10,12,15]. Indeed, for interme-
diate quenches where this perturbative treatment is no
longer valid, this connection breaks down [17–19], and
for small quenches within the ordered phase, anomalous
DQPTs can appear even without any OP zeros occurring
over all investigated timescales [10,11,20]. As such, it
has been argued that the periodic DQPT behavior seen
for large quenches may be a mere manifestation of effec-
tive two-level system dynamics [21].
Here, we address this argument by investigating
DQPTs in models exhibiting “state transfer”—a dynam-
ical process where the wave function evolves cyclically
between the states of a given manifold—due to quantum
many-body scarring [22–24]. Such models possess a small
number of anomalous eigenstates throughout their spec-
trum, leading to oscillatory dynamics from a few specific
states in an otherwise thermalizing system. We will fo-
cus on a formulation of the lattice Schwinger model [25]
known as the spin-SU(1) quantum link model (QLM)
[26,27]. This model has been shown to exhibit quan-
tum many-body scarring for massless quenches starting
in the maximal-flux (extreme) vacua for a wide range of
spin values [28,29]. We show in these models that peri-
odic DQPTs arise within complex many-body dynamics
that is beyond two-level systems, and where a DQPT sig-
nals state transfer from one vacuum to another within a
(2S+1)-fold degenerate vacuum manifold; see Fig. 1. We
find no direct connection between DQPTs and OP ze-
ros for integer S. For half-integer S, a DQPT is directly
connected to an OP zero only when the DQPT signals
a transfer between intermediate minimal-flux vacua of
opposite flux sign. We further show that models where
DQPT behavior resembles two-level system dynamics are
a special case of our general picture.
II. MODEL
We consider the spin-SU(1) QLM, given by the Hamil-
tonian [26,27,30]
ˆ
H=
L
X
j=1 J
2pS(S+ 1)ˆσ−
jˆs+
j,j+1 ˆσ−
j+1 + H.c.
+µˆσz
j+κ2
2ˆsz
j,j+12,(1)
where we have adopted particle-hole and Jordan–Wigner
transformations [31,32]. The Pauli operator ˆσz
jdescribes
the matter occupation on site jwith mass µ, and the
spin-Soperators ˆs+
j,j+1/pS(S+ 1) and ˆsz
j,j+1 represent
the gauge and electric fields, respectively, residing on the
link between sites jand j+ 1. The tunneling constant is
J, which we shall set to unity as the overall energy scale,
κis the gauge-coupling strength, and Lis the number of
sites.
The generator of the U(1) gauge symmetry of Hamil-
tonian (1) is
ˆ
Gj= (−1)jˆsz
j−1,j + ˆsz
j,j+1 +ˆσz
j+ˆ
1
2,(2)
and gauge-invariant states |ϕ⟩satisfy ˆ
Gj|ϕ⟩=gj|ϕ⟩,∀j,
where gj/(−1)j∈{−2S, . . . , 2S+1}. We will work in the
physical sector gj=0,∀j.
A building block of the U(1) QLM has been exper-
imentally realized for S→∞ in a cold-atom setup [33].
Large-scale implementations of the spin-1/2 U(1) QLM
on a Bose–Hubbard superlattice have been employed to
observe gauge invariance [34] and thermalization dynam-
ics [35].
III. QUENCH DYNAMICS
We now present time-evolution results obtained
through the infinite matrix product state (iMPS) tech-
nique based on the time-dependent variational principle
[36–39]. This technique works directly in the thermody-
namic limit, and also allows us to directly detect DQPTs
as level crossings between the logarithms of the eigen-
values of the MPS transfer matrix [10,40,41], without
any need for finite-size scaling with L. DQPTs for gauge
theories have already been studied in the context of the
spin-1/2 U(1) QLM [13,42,43], and also for S≥1/2 [44]
as well as in the Schwinger model [14,45], but not in the
context of quantum many-body scarring. For the most
stringent calculations that we have performed in iMPS
for this work, we find convergence for a maximal bond
dimension of 550 and a time-step of 0.0005/J.
We are interested in the dynamics of the experimen-
tally relevant return rate (RR)
λ(t) = min
mzλmz(t),(3a)
λmz(t) = −lim
L→∞
1
Lln
⟨ψmz
0|ψ(t)⟩
2,(3b)
which has recently been used to identify DQPTs in a
trapped-ion experiment [46]. Here, |ψmz
0⟩is the set
of vacua with mz∈−S, . . . , S, which are the (2S+1)-
fold degenerate ground states of Hamiltonian (1) for
κ/J=0 and µ/J→∞. We call a vacuum extreme when
|mz|=S, and intermediate when |mz|<S. We can rep-
resent a vacuum state on the two-site two-link unit cell
as |ψmz
0⟩=|−1, mz,−1,−mz⟩, indicating the eigenvalue
−1 of ˆσz
jat each (empty) matter site and the eigenvalue
mzof ˆsz
j,j+1 on each link, but we will often refer to this