2
of quantum many-body models in general, whose key
properties of quantum entanglement and superposition
are particularly sensitive to interactions with the envi-
ronment. Prominent examples of the detrimental effects
of decoherence on quantum many-body systems include
1/f noise in superconducting quantum interference de-
vices (SQUIDs) that undermines superconducting qubits
[58–63]. Given that superconducting qubits, as well as
other platforms, have been of great recent interest in the
quantum simulation of gauge theories [26,27], suppress-
ing 1/f noise sources using efficient and experimentally
feasible schemes becomes of central importance.
In this work, using exact diagonalization calculations
and time-dependent perturbation theory, we demonstrate
how the principle of linear gauge protection, initially de-
vised to control coherent gauge-breaking errors [51], can
be employed to suppress the growth of the gauge vio-
lations due to incoherent errors with spectral form 1/fβ
(β > 0) as 1/V β, where Vis the protection strength. The
rest of this paper is organized as follows: We briefly re-
view the concept of linear gauge protection in Sec. II, and
1/fβnoise and the corresponding Bloch–Redfield formal-
ism in Sec. III. We present our main numerical results in
Sec. IV. We finally conclude and provide an outlook in
Sec. V. We include Appendix Afor a derivation of the
Bloch–Redfield equation employed for our analysis, Ap-
pendix Bfor our derivations in time-dependent perturba-
tion theory, in addition to Appendix Cwhere we provide
supplemental numerical results.
II. LINEAR GAUGE PROTECTION
Let us consider an Abelian gauge theory described by
the Hamiltonian ˆ
H0, and whose gauge symmetry is gen-
erated by the operator ˆ
Gj, where jdenotes a site on a
lattice of length L. The gauge invariance of ˆ
H0is encoded
in the commutation relations ˆ
H0,ˆ
Gj= 0,∀j. The set
of gauge-invariant states {|ψi} is defined as the simulta-
neous eigenstates of the generators: ˆ
Gj|ψi=gj|ψi,∀j.
A set of these eigenvalues g= (g1, g2, . . . , gL) over the
volume of the system defines a unique gauge superselec-
tion sector, the projector onto which is ˆ
Pg. One can
further define a target or physical gauge superselection
sector gtar = (gtar
1, gtar
2, . . . , gtar
L) in which one wishes to
restrict the dynamics in an experiment, for example.
In experimental implementations of gauge theories, ˆ
H0
is mapped onto the microscopic degrees of freedom of
a quantum simulator. In general, unavoidable gauge
symmetry-breaking errors λˆ
H1at strength λwill arise
in this process either due to higher orders in the per-
turbation theory used to perform the mapping, or in ex-
perimental imperfections in equipment. Even when per-
turbative, these errors can generate gauge violations that
grow as λ2t2over evolution time t, which in turn lead to a
complete departure from faithful gauge-theory dynamics
beyond timescales t∝1/λ [31].
In order to suppress these errors in a controlled way,
the concept of linear gauge protection was introduced in
Ref. [51]. It entails adding the protection term
Vˆ
HG=VX
j
cjˆ
Gj,(1)
where Vis the protection strength. The sequence cj
can be chosen to be rational and satisfying the condition
Pjcjg−gtar
j= 0 ⇐⇒ gj=gtar
j,∀j. In this case,
the sequence is said to be compliant, and, for a volume-
independent and sufficiently large V, the gauge violation
is controlled up to times exponential in V[51,64]. Al-
though Vis volume-independent, the sequence cjwould
have to grow (not faster than) exponentially with system
size in order to satisfy the compliance condition. This
renders the compliant sequence somewhat inconvenient
for large-scale gauge-theory quantum simulators such as
those realized in recent cold-atom setups [23,24].
However, reality turns out to be more forgiving, and
even simple noncompliant sequences such as cj= (−1)j
can give excellent protection in the target sector against
gauge errors up to all accessible evolution times in both fi-
nite systems [51] and the thermodynamic limit [65]. This
can be explained through the coherent quantum Zeno ef-
fect [66–69], which guarantees that upon adding the pro-
tection term (1) an effective Zeno Hamiltonian ˆ
HZ=
ˆ
H0+λˆ
Pgtar ˆ
H1ˆ
Pgtar emerges that faithfully reproduces
the dynamics of the faulty gauge theory ˆ
H0+λˆ
H1+Vˆ
HG
up to timescales linear in Vin a worst-case scenario [51].
For certain gauge theories, the full local generator ˆ
Gj
may be too challenging to realize in an experiment [19],
in which case the linear gauge protection as given in
Eq. (1) becomes impractical. Nevertheless, a powerful
workaround exists based on local pseudogenerators ˆ
Wj,
which are identical to the full local generators ˆ
Gjin the
target sector, but not necessarily outside of it [52]. For-
mally, they satisfy the relation
ˆ
Wj|φi=gtar
j|φi ⇐⇒ ˆ
Gj|φi=gtar
j|φi.(2)
One can then extend the principle of linear gauge protec-
tion to one in terms of the local pseudogenerator, with
protection term
Vˆ
HW=VX
j
cjˆ
Wj,(3)
where the same rules apply for the sequence cjas in
the case of Eq. (1). Note that even though ˆ
H0com-
mutes with ˆ
Gj, it generally does not commute with
ˆ
Wj, with the latter associated with a local symmetry
richer than that generated by ˆ
Gj[70]. The result-
ing Zeno Hamiltonian when protecting with Eq. (3) is
ˆ
HZ=ˆ
Pgtar ˆ
H0+λˆ
H1ˆ
Pgtar , under which the dynam-
ics of the faulty gauge theory ˆ
H0+λˆ
H1+Vˆ
HWcan be
faithfully reproduced up to times at least linear in V[52].
In terms of purely unitary errors, extensive numeri-
cal simulations in exact diagonalization (ED) and infi-