of 7
1
Supplemental Material for “Symmetry breaking and error correction in open
quantum systems”
Simon Lieu,
1
,
2
Ron Belyansky,
1
,
2
Jeremy T. Young,
1
Rex Lundgren,
1
,
2
Victor V. Albert,
3
,
4
Alexey V. Gorshkov
1
,
2
1
Joint Quantum Institute, NIST/University of Maryland, College Park, MD 20742, USA
2
Joint Center for Quantum Information and Computer Science,
NIST/University of Maryland, College Park, Maryland 20742 USA
3
Institute for Quantum Information and Matter and Walter Burke Institute for
Theoretical Physics, California Institute of Technology, Pasadena, CA 91125, USA
4
National Institute of Standards and Technology, Gaithersburg, MD 20899, USA
In Sec. 1, we analytically show that the dissipative gap closes at the critical point by utilizing an exact solution
for the Lindblad spectrum [Fig. 1(b) in the main text]. Sec. 2 exhibits numerical evidence for a noiseless subsystem
steady state in the strong-broken phase (away from
L
0
). Sec. 3 tracks the evolution of the state throughout the error
protocol in the main text. We show numerical evidence for the state structure defined in Eq. (5) of the main text
for errors which keep the model in the strong-broken phase. Sec. 4 uses the asymptotic projection method to confirm
perfect fidelity recovery in the thermodynamic limit, in agreement with the direct numerical evolution discussed in
the main text. In Sec. 5 we show that a
classical bit
encoded into the steady state is protected against perturbations
which keep the Lindbladian in the
weak-broken
phase.
1. CLOSING OF THE DISSIPATIVE GAP AT THE CRITICAL POINT
We show that an extensive number of spectral eigenvalues touch zero at the critical boundary [Fig. 1(a) in the main
text] when approaching from the unbroken phase in the thermodynamic limit. We utilize Prosen’s “third quantization”
technique which allows us to fully diagonalize a quadratic Lindbladian [S1, S2]. For the Hamiltonian (3) in the presence
of one-photon loss only (i.e. the weak transition), the Lindbladian can be expressed as
L
=

+
β
+
β
+
+

β
β
,
where
β
are bosonic superoperators satisfying generalized commutation relations [
β
i
j
] =
δ
ij
. These excite a quantum
of “complex energy”

±
=
κ
1
±
4
λ
2
ω
2
, where the (unique) steady state is annihilated by all quasiparticles
β
±
ρ
ss
= 0, and the many-body spectrum is built from these single-particle excitations
L
[(
β
+
)
n
(
β
)
m
ρ
ss
] = (
n
+
+
m
)[(
β
+
)
n
(
β
)
m
ρ
ss
]. The single-particle spectrum touches zero at
κ
1
=
4(
λ/ω
)
2
1, which coincides with the
emergence of a non-zero order parameter (see main text). This implies that an infinite number of eigenvalues of
L
are
zero at the critical point of the weak transition from 1 steady state to 2 steady states. We plot both the single-particle
spectrum and match it with many-body numerics in Fig. S1. [Fig. S1(a) and Fig. 1(b) are equivalent; here we plot
the real and imaginary parts side by side.] The numerical spectrum deviates from analytical predictions only near
the critical boundary due to truncation of the Hilbert space dimension. Note that the analytical and numerical plots
are only valid in the unbroken phase. The steady state has an infinite number of photons in the broken phase, hence
any finite-size Hilbert space will not produce a converged spectrum. Finite-size scaling [Fig. 1(d)] suggests that two
eigenvalues are exponentially close to zero in the weak-broken phase with a dissipative gap to the rest of the modes.
0
1
λ/ω
0
4
Re[
Λ
]
(a)
0
1
λ/ω
2
0
2
Im[Λ
]
(b)
FIG. S1. Analytical single-particle spectrum (red lines) and numerical many-body spectrum (black dots) with
κ
1
= 2
2
=
κ
d
= 0. The many-body spectrum comes in integer multiples of the single-particle excitations. As the system approaches
the critical point from the unbroken phase, the single-particle spectrum touches zero at the phase boundary
λ/ω
1
.
1. The
numerical spectrum starts to deviate from the analytical predictions near the transition due to truncation of the Hilbert space
dimension
d
Hilbert
= 70. We plot up to 25 eigenvalues closest to zero for clarity.
2
2. NOISELESS SUBSYSTEM IN THE STRONG-BROKEN PHASE
We demonstrate that the model described in the main text possesses a qubit steady-state structure in the thermo-
dynamic limit of the strong-broken phase. In particular, we will show that the four right eigenoperators with zero
eigenvalue can be written in the form
r
μν
=
|
μ
〉〈
ν
|⊗
z
with (
μ,ν
)
(+
,
). This is called a noiseless subsystem (NS)
if
z
is mixed, and a decoherence-free subspace (DFS) if
z
is pure [S3–S5].
The four steady-state right eigenoperators belonging to the different parity sectors are
r
F
++
=
(
s
++
0
0
0
)
, r
F
−−
=
(
0
0
0
s
−−
)
, r
F
+
=
(
0
s
+
0
0
)
, r
F
+
=
(
0
0
s
+
0
)
(S1)
in the Fock basis [
|
0
,
|
2
,
|
4
,...,
|
1
,
|
3
,
|
5
,...
]
T
. They each satisfy
L
(
r
) = 0 (in the thermodynamic limit).
Since
s
++
,s
−−
are guaranteed to be Hermitian matrices, we can diagonalize them via a unitary transformation
U
= Diag[
U
+
,U
] which relates the Fock basis to the diagonal basis
r
d
i
=
U
r
F
i
U
. In this new basis, the eigenoperators
are
r
d
++
=
(
z
++
0
0
0
)
, r
d
−−
=
(
0
0
0
z
−−
)
, r
d
+
=
(
0
z
+
0
0
)
, r
d
+
=
(
0
0
z
+
0
)
,
(S2)
where
z
++
,z
−−
are diagonal by construction, and
z
+
,z
+
are diagonal in the thermodynamic limit. We will show
that
z
++
=
z
−−
=
z
+
=
z
+
in this limit, which implies that the system hosts a NS or a DFS.
In the special limit
ω
=
κ
d
=
κ
1
= 0
6
= 0
2
6
= 0, any pure superposition of even and odd cat states remains
steady, as discussed in the main text. Thus
z
++
=
z
−−
=
z
+
=
z
+
= Diag[1
,
0
,
0
,
0
,...
], which implies a DFS.
We now consider a parameter regime away from this limit but within the strong-broken phase. We start by adding
dephasing:
ω
=
κ
1
= 0
d
6
= 0
6
= 0
2
6
= 0. We will numerically show that the
z
matrices are equal and not
pure. For the matrix distance, we choose the trace distance
D
t
(
A,B
) = Tr[
(
A
B
)
2
]
/
2. In Fig. S2(a,b), we plot
D
t
(
z
++
,z
−−
) and
D
t
(
z
++
,z
+
) as the system approaches the thermodynamic limit
λ/κ
2
=
N
→ ∞
. Indeed, we
find that the matrices
z
++
,z
−−
,z
+
all converge to a single matrix as
N
is increased. (
z
+
and
z
+
are related
by Hermiticity.) Additionally, in Fig. S2(c), we show that
z
++
is a non-pure matrix with elements that fall off as
(
z
++
)
ii
exp [
i
]. The purity of
z
++
degrades
with
N
(not shown). We conclude that the system tends to a noiseless
subsystem in the thermodynamic limit, since the
z
±±
all converge to a single non-pure matrix. [For completeness, in
Fig. S2(d), we show that the smallest eigenvalue in the off-diagonal sector indeed tends to zero exponentially quickly
with
N
. The steady-state degeneracy is split by an exponentially small factor, characteristic of symmetry-breaking
transitions.]
We repeat this analysis in the limit of no dephasing but non-zero
ω
:
κ
d
=
κ
1
= 0
6
= 0
6
= 0
2
6
= 0. Fig. S3
shows that the
z
±±
converge to a single non-pure matrix in the thermodynamic limit, similar to the case of dephasing.
We therefore conclude that a generic model in the strong-broken phase possesses a noiseless subsystem, whilst a
decoherence-free subspace exists at a special point
L
0
in the phase diagram.
10
15
20
N
10
9
8
7
6
(a)
log
10
[
|
D
t
(
z
++
,z
−−
)
|
]
10
15
20
N
7
6
5
4
(b)
log
10
[
|
D
t
(
z
++
,z
+
)
|
]
10
15
20
N
10
9
8
7
(d)
log
10
[
Re(Λ
+
)]
0
3
6
9
i
8
6
4
2
0
(c)
log
10
[(
z
++
)
ii
]
FIG. S2.
Parameters:
λ/κ
2
=
N,κ
d
= 0
.
03
=
κ
1
= 0, i.e. non-zero dephasing. (a,b) The trace norm
D
t
(
A,B
) =
Tr[
(
A
B
)
2
]
/
2 between the different right eigenoperators with zero eigenvalue goes to zero in the thermodynamic limit
N
→∞
. (c) Diagonal matrix elements of
z
++
for
N
= 20. The matrix is not pure, with elements scaling as (
z
++
)
ii
exp [
ci
]
for some
c >
0. (d) The off-diagonal symmetry sector of the Lindbladian acquires an eigenvalue of zero as
N
→∞
. Here Λ
+
is the smallest eigenvalue in the off-diagonal sector.
3
10
15
20
N
10
8
6
(a)
log
10
[
|
D
t
(
z
++
,z
−−
)
|
]
10
15
20
N
4
.
1
4
.
0
3
.
9
3
.
8
3
.
7
(b)
log
10
[
|
D
t
(
z
++
,z
+
)
|
]
10
15
20
N
10
8
6
(d)
log
10
[
Re(Λ
+
)]
0
3
6
9
i
15
10
5
0
(c)
log
10
[(
z
++
)
ii
]
FIG. S3.
Parameters:
λ/κ
2
=
N,ω/λ
= 0
.
5
d
=
κ
1
= 0. (a,b) The trace norm between the different right eigenoperators
with zero eigenvalue goes to zero in the thermodynamic limit
N
→∞
. (c) Diagonal matrix elements of
z
++
for
N
= 20. The
matrix is not pure, with elements scaling as (
z
++
)
ii
exp [
ci
] for some
c >
0. (d) The off-diagonal symmetry sector of the
Lindbladian acquires an eigenvalue of zero as
N
→∞
. Here Λ
+
is the smallest eigenvalue in the off-diagonal sector.
3. EVOLUTION FROM DECOHERENCE-FREE SUBSPACE TO NOISELESS SUBSYSTEM
We now track the state throughout the error protocol described in the main text for both dephasing errors and
Hamiltonian-frequency errors. Our analysis will confirm that the state can be written as a qubit tensored with a mixed
state thoughout the entire quench protocol, i.e. the structure described in Eq. (5) in the main text. An analytical
understanding of this mechanism requires an exact solution for the steady states in the entire strong-broken phase—an
important direction for future work.
We prepare the system in a pure steady state of
L
0
:
ρ
i
=
(
|
c
e
|
2
c
e
c
o
c
e
c
o
|
c
o
|
2
)
(S3)
in the basis of even and odd cat states
|
α
e
,
|
α
o
,
where
|
c
e
|
2
+
|
c
o
|
2
= 1 and
L
0
(
ρ
i
) = 0. We evolve this initial state
with an error to a “middle” state
ρ
m
(
τ
q
) =
e
(
L
0
+
L
)
τ
q
ρ
i
.
(S4)
We wish to show that this middle state can be written in the form
ρ
m
(
τ
q
) =
(
|
c
e
|
2
c
e
c
o
c
e
c
o
|
c
o
|
2
)
M
(S5)
for some
M
which is not necessarily pure.
We numerically solve for
ρ
m
(
τ
q
) via Eq. (S4) for arbitrary quench times and
L
0
+
L
in the strong-broken phase.
We then split the matrix up into symmetry sectors in the Fock basis
ρ
m
=
|
c
e
|
2
ρ
F
++
+
|
c
o
|
2
ρ
F
−−
+ (
c
e
c
o
ρ
F
+
+
h.c.
).
The four operators belonging to the different parity sectors are
ρ
F
++
=
(
x
++
0
0
0
)
, ρ
F
−−
=
(
0
0
0
x
−−
)
, ρ
F
+
=
(
0
x
+
0
0
)
, ρ
F
+
=
(
0
0
x
+
0
)
(S6)
in the Fock basis [
|
0
,
|
2
,
|
4
,...,
|
1
,
|
3
,
|
5
,...
]
T
. Since
x
++
,x
−−
are guaranteed to be Hermitian matrices, we can
diagonalize them via a unitary transformation
V
= Diag[
V
+
,V
] which relates the Fock basis to the diagonal basis
ρ
d
i
=
V
ρ
F
i
V
. In this new basis, the eigenoperators are
ρ
d
++
=
(
M
++
0
0
0
)
, ρ
d
−−
=
(
0
0
0
M
−−
)
, ρ
d
+
=
(
0
M
+
0
0
)
, ρ
d
+
=
(
0
0
M
+
0
)
,
(S7)
where all the
M
s are diagonal by construction. We now show that all
M
s converge to a single matrix in the
thermodynamic limit, confirming the form of Eq. (S5).
We plot the trace distance between the different
M
s for both short and long quench times
τ
q
λ
[10
2
,
10
2
]. In
Fig. S4, we consider a quench in the dephasing strength. Indeed, the trace distance between the different
M
s goes to
zero exponentially fast as a function of
N
, which suggests that the ansatz in Eq. (S5) is correct in the limit
N
→∞
.
4
5
10
15
N
4
8
12
(a)
log
10
[
|
D
t
(
M
++
,M
−−
)
|
]
log
10
τ
q
λ
=2
log
10
τ
q
λ
=1
log
10
τ
q
λ
=0
log
10
τ
q
λ
=-1
log
10
τ
q
λ
=-2
5
10
15
N
4
8
12
(b)
log
10
[
|
D
t
(
M
++
,M
+
)
|
]
log
10
τ
q
λ
=2
log
10
τ
q
λ
=1
log
10
τ
q
λ
=0
log
10
τ
q
λ
=-1
log
10
τ
q
λ
=-2
4
2
0
2
log
10
[
τ
q
λ
]
0
.
85
1
.
00
(c)
Tr[
M
2
++
]
FIG. S4. Parameters:
λ/κ
2
=
N,κ
d
= 0
.
03
=
κ
1
= 0
,c
e
= 1
/
2
,c
o
=
i/
2. (a) The trace distance between
M
++
and
M
−−
goes to zero exponentially fast in
N
. (b) Analogous behavior is observed for
M
++
and
M
+
. (c)
N
= 15
,
the red line
is the time scale set by the inverse dissipative gap
τ
g
= ∆
1
g
of
L
0
+
L
. The state is approximately pure for short quenches
compared to this time scale, while it settles to its (mixed) steady-state value for quenches longer than this timescale.
5
10
15
N
4
8
12
(a)
log
10
[
|
D
t
(
M
++
,M
−−
)
|
]
log
10
τ
q
λ
=2
log
10
τ
q
λ
=1
log
10
τ
q
λ
=0
log
10
τ
q
λ
=-1
log
10
τ
q
λ
=-2
5
10
15
N
4
8
12
(b)
log
10
[
|
D
t
(
M
++
,M
+
)
|
]
log
10
τ
q
λ
=2
log
10
τ
q
λ
=1
log
10
τ
q
λ
=0
log
10
τ
q
λ
=-1
log
10
τ
q
λ
=-2
4
2
0
2
log
10
[
τ
q
λ
]
0
.
99
1
.
00
(c)
Tr[
M
2
++
]
FIG. S5.
Parameters:
λ/κ
2
=
N,ω/λ
= 0
.
5
d
=
κ
1
= 0
,c
e
= 1
/
2
,c
o
=
i/
2. (a) The trace distance between
M
++
and
M
−−
goes to zero exponentially fast in
N
. (b) Analogous behavior is observed for
M
++
and
M
+
. (c)
N
= 15
,
the red line
is the time scale set by the inverse dissipative gap
τ
g
= ∆
1
g
of
L
0
+
L
. The state is approximately pure for short quenches
compared to this time scale, while it settles to its (mixed) steady-state value for quenches longer than this timescale.
We also track the purity of this matrix: At quench times that are short compared to the timescale set by the dissipative
gap (red line), the middle state remains approximately pure, whilst longer quenches imply that the system settles into
its new steady state, which is mixed (see previous section). Analogous behavior is observed for a quench in frequency
(Fig. S5).
4. ASYMPTOTIC PROJECTION
We verify the perfect recovery of the fidelity observed in Fig. 2 of the main text via the asymptotic projection
method [S6]. Fig. 2 shows that qubit cat states will self correct via the environment if
L
0
+
L
remains in the strong
symmetry-broken phase. This behavior can be understood via perturbation theory for short quenches (compared to
the time scale set by the dissipative gap) [S7]. Here, we consider long quench times where the system evolves into
the steady state of
L
0
+
L
. Remarkably, such a drastic error can still be passively corrected via the environment
L
0
.
We provide simple expressions relating the initial, intermediate, and final states by projecting onto the corresponding
steady state manifolds.
Defining our initial state as
ρ
i
, we evolve it with an error (
L
0
+
L
) to a “middle” state
ρ
m
(
τ
q
) =
e
(
L
0
+
L
)
τ
q
ρ
i
. We
then evolve the state with
L
0
for an infinite time to reach the final state
ρ
f
(
τ
q
) = lim
τ
→∞
e
L
0
τ
ρ
m
(
τ
q
)
.
We will discuss
how
ρ
i,m,f
relate to one another in this protocol when
τ
q
is much longer than the inverse dissipative gap of
L
0
+
L
.
We first prepare the system in a pure steady state of
L
0
,
ρ
i
=
|
a
|
2
r
0
++
+
|
b
|
2
r
0
−−
+
a
br
0
+
+
ab
r
0
+
,
(S8)
where
r
0
++
=
|
α
e
α
|
e
,
r
0
−−
=
|
α
o
α
|
o
,
r
0
+
=
|
α
e
α
|
o
,
r
0
+
=
|
α
o
α
|
e
;
|
α
e/o
is the even/odd cat state, and
L
0
(
r
0
±±
) = 0. To find
ρ
m
, it is useful to define the right and left eigenoperators of the error:
(
L
0
+
L
)( ̃
r
j
) =
̃
Λ
j
( ̃
r
j
)
,
(
L
0
+
L
′†
)(
̃
l
j
) =
̃
Λ
j
(
̃
l
j
)
,
(S9)
5
10
15
20
N
8
10
log
10
[
|
1
γ
f
|
]
(a)
10
15
20
N
10
8
6
log
10
[
|
1
γ
f
|
]
(b)
FIG. S6.
Scaling of
|
1
γ
f
|
as a function of
N
for (a) a dephasing error
λ/κ
2
=
N,κ
d
= 0
.
03 and (b) a frequency error
λ/κ
2
=
N,ω/λ
= 0
.
03.
γ
f
approaches one exponentially fast in
N
for both cases.
where the spectrum
{
̃
Λ
}
and eigenoperators determine the dynamics under
L
0
+
L
. Assuming that the error keeps
the system in the strong-broken phase, we know that two eigenvalues will be exactly zero
̃
Λ
0
++
=
̃
Λ
0
−−
= 0 and two
eigenvalues will be exponentially close to zero
̃
Λ
0
+
= (
̃
Λ
0
+
)
e
N
. We label the eigenvalue of the first “excited”
state (above these four) as
̃
Λ
g
, which sets the dissipative gap in the thermodynamic limit. The exact expression for
ρ
m
(
τ
q
) reads
ρ
m
(
τ
q
) =
j
exp[
̃
Λ
j
τ
q
]Tr[
̃
l
j
ρ
i
] ̃
r
j
,
(S10)
where we have used the orthogonality relation Tr[
̃
l
j
̃
r
k
] =
δ
jk
.
Re[
̃
Λ
1
j
] sets the lifetime of each eigenoperator.
Consider a quench time that obeys
Re[
̃
Λ
1
g
]

τ
q
 −
Re[(
̃
Λ
0
+
)
1
]
e
N
. This quench is long enough for the
system to relax into the new steady state but not so long that coherences are lost. In this regime,
ρ
m
will tend to the
following matrix
t
m
lim
N
→∞
ρ
m
(
τ
q
) =
t
m
, t
m
=
|
a
|
2
̃
r
0
++
+
|
b
|
2
̃
r
0
−−
+
[
a
m
̃
r
0
+
+
h.c.
]
, γ
m
= Tr
[
(
̃
l
0
+
)
r
0
+
]
.
(S11)
If
τ
q
is longer than
Re[
̃
Λ
1
g
], then all excitations will vanish and we will be left with the projection onto
the steady-state manifold of the error. We have confirmed this numerically by doing the full time evolution
ρ
m
= exp [(
L
+
L
)
τ
q
]
ρ
i
and comparing the resulting matrix with
t
m
. Indeed, the trace distance
D
t
(
ρ
m
,t
m
) =
Tr[
(
ρ
m
t
m
)
2
]
/
2) goes to zero exponentially quickly in
N
. We have thus found a simple expression for
ρ
m
(
τ
q
) for
this range of
τ
q
.
Having understood the structure of this intermediate state,
ρ
m
t
m
, we now project this state back onto the
steady-state manifold of
L
0
. Without any additional approximations, the resulting state is
lim
N
→∞
ρ
f
=
|
a
|
2
r
0
++
+
|
b
|
2
r
0
−−
+
γ
f
a
br
0
+
+
γ
f
ab
r
0
+
, γ
f
= Tr[(
̃
l
0
+
)
r
0
+
]Tr[(
l
0
+
)
̃
r
0
+
]
.
(S12)
We see that the final state is very simply related to the initial state via the
γ
f
parameter in Eq. (S12). Moreover,
numerically we observe that
γ
f
approaches 1 exponentially fast in the thermodynamic limit, depicted in Fig. S6 for
both the case of (a)
κ
d
6
= 0 and (b)
ω
6
= 0. (We have also checked that
γ
m
approaches 1 in the same limit.) This
implies that the final state
ρ
f
is indeed expected to return to its initial (pure) state
ρ
i
in the thermodynamic limit.
Structure of the left eigenoperators
̃
l
In Sec. 3 and earlier in this Section, we saw that the initial state settles into the noiseless subsystem of the
intermediate Lindbladian
L
0
+
L
without losing any coherences as
N
→∞
. We would like to find a simple explanation
for this behavior. This evolution would be accounted for (in the limit
N
→ ∞
) if the left eigenoperators of
L
0
+
L
with zero eigenvalue are equal to the identity in each symmetry sector, since, in this case,
γ
m
= Tr[(
l
0
+
)
̃
r
0
+
] =
Tr[ ̃
s
0
+
] = Tr[ ̃
z
0
+
] = 1 where in the last step we have used Tr[ ̃
z
0
+
] = Tr[ ̃
z
0
++
] = 1. (See Sec. 2 for definitions of
r,s,z
.) We will show that this is indeed true. Splitting up the left eigenoperators into symmetry sectors, we have
̃
l
F
++
=
(
y
++
0
0
0
)
,
̃
l
F
−−
=
(
0
0
0
y
−−
)
,
̃
l
F
+
=
(
0
y
+
0
0
)
,
̃
l
F
+
=
(
0
0
y
+
0
)
.
(S13)
6
0
9
0
9
(a)
q
+
,N
= 10
1
.
0
0
.
5
0
.
0
0
.
5
1
.
0
0
9
0
9
(b)
q
+
,N
= 20
1
.
0
0
.
5
0
.
0
0
.
5
1
.
0
0
9
0
9
(c)
q
+
,N
= 30
1
.
0
0
.
5
0
.
0
0
.
5
1
.
0
FIG. S7. Plot of a 10
×
10 block of
q
+
; all elements are real. Parameters:
κ
2
= 1
/N,κ
d
= 0
.
03
=
κ
1
= 0. As the
system approaches the thermodynamic limit, the matrix tends to the identity.
0
9
0
9
(a)
q
+
,N
= 10
1
.
0
0
.
5
0
.
0
0
.
5
1
.
0
0
9
0
9
(b)
q
+
,N
= 20
1
.
0
0
.
5
0
.
0
0
.
5
1
.
0
0
9
0
9
(c)
q
+
,N
= 30
1
.
0
0
.
5
0
.
0
0
.
5
1
.
0
FIG. S8. Plot of a 10
×
10 block of
q
+
; all elements are real. Parameters:
κ
2
= 1
/N,ω/λ
= 0
.
5
d
=
κ
1
= 0. As the system
approaches the thermodynamic limit, the matrix tends to the identity.
As before, we are in the Fock basis [
|
0
,
|
2
,
|
4
,...,
|
1
,
|
3
,
|
5
,...
]
T
. Then
y
++
=
y
−−
=
I
since any arbitrary
initial state must have unit overlap with the steady-state solutions with non-zero trace. Now we switch from the Fock
basis to the diagonal basis of
r
,
r
d
i
=
U
r
F
i
U
,
̃
l
d
i
=
U
̃
l
F
i
U
, and obtain
̃
l
d
++
=
(
q
++
0
0
0
)
,
̃
l
d
−−
=
(
0
0
0
q
−−
)
,
̃
l
d
+
=
(
0
q
+
0
0
)
,
̃
l
d
+
=
(
0
0
q
+
0
)
.
(S14)
Again,
q
++
=
q
−−
=
I
; we shall now probe the structure of the off-diagonal matrix
q
+
.
In this basis, the four
right
eigenoperators
r
of
L
0
+
L
with zero eigenvalue are just a single diagonal matrix
z
in each
of the four symmetry quadrants in the thermodynamic limit (see Sec. 2). This matrix
z
is not pure, and in principle
has infinite rank although its eigenvalues fall off exponentially quickly as a function of the index, i.e.
z
jj
e
cj
for
some
c >
0. In the case of a noiseless subsystem with full rank
z
, Ref. [S8] proved that the corresponding conserved
quantity must be the identity in each symmetry sector for a finite-dimensional Hilbert space. Since our bosonic model
has an infinite-dimensional Hilbert space, these results do not immediately apply. Nevertheless, we numerically show
that the conserved quantities approach the identity in the thermodynamic limit.
In Fig. S7, we plot the elements of a 10
×
10 block of the matrix
q
+
for the case of non-zero dephasing. Indeed,
we find that the matrix tends to the identity as we approach the thermodynamic limit. The matrix
q
+
acquires
off-diagonal terms at entries where the corresponding matrix elements
z
jj
are small, i.e. we are limited by numerical
precision. Analogous behavior is observed for the case of non-zero
ω
, depicted in Fig. S8. So indeed we expect
lim
N
→∞
q
+
=
I
for the full rank noiseless subsystem. This explains why
ρ
i
does not lose coherences when relaxing
into the steady state of
L
0
+
L
.
5. PASSIVE PROTECTION OF A CLASSICAL BIT IN THE WEAK-BROKEN PHASE
Up to now, we have focused on describing the quantum error correcting properties for a qubit encoded in the steady
state of
L
0
subject to errors that keep it in the
strong-broken
phase. Here, we show that a classical bit encoded
into the steady state of
L
0
will be passively protected against any error which keeps the model in the
weak-broken
phase. This region of parameter space includes single-photon loss, which is non-negligible in experimental setups and
represents the dominant decoherence mechanism for photonic cat qubits.
We consider the same protocol as in the main text, but with a restricted initial state: Initialize the system in the
state
ρ
i
=
c
|
α
〉〈
α
|
+ (1
c
)
|−
α
〉〈−
α
|
where
c
= 0 or 1, which represents the classical bit and satisfies
L
0
(
ρ
i
) = 0.
(Any choice of
c
[0
,
1] will also exhibit protection.) Then quench the state with an “error” for an arbitrary time
τ
q
to obtain
ρ
m
= exp [(
L
0
+
L
)
τ
q
](
ρ
i
). Finally, turn off the error and evolve the system with
L
0
for a long time such
7
5
10
N
0
.
0
0
.
5
1
.
0
F
(
ρ
i
f
)
(a)
κ
1
= 0.50
κ
1
= 2.50
2
6
10
N
1
2
3
log
10
[1
F
(
ρ
i
f
)]
(b)
FIG. S9. Fidelity of the initial and final
classical
-bit state with
c
= 1
,λ/κ
2
=
N,κ
d
=
ω
= 0
q
λ
= 10,
F
(
ρ
i
f
) =
Tr[
ρ
i
ρ
f
ρ
i
]
2
. Quenches to the weak-broken phase (black dots) have a fidelity that tends to one in the thermodynamic
limit, while quenches to the weak-unbroken phase (red dots) do not. (b) Same parameters as in (a) with
κ
1
= 0
.
5; the fidelity
tends to one exponentially fast in
N
.
that it reaches its steady state:
ρ
f
= lim
t
→∞
exp [
L
0
t
](
ρ
m
). We show that any error which keeps the model in the
weak-broken
phase is correctable passively.
In Fig. S9, we plot the fidelity between the initial and final states after a long quench of single-photon-loss error
L
1
=
κ
1
a
. Only if the error keeps the model in the weak-broken phase (i.e.
κ
1
/λ <
2) does the classical bit recover
its initial state in the thermodynamic limit.
We can understand this behavior by considering the steady-state structure for a generic system in the weak-broken
phase. In a parity basis, it assumes the form
ρ
ss
=
(
1
/
2
c
1
/
2
c
1
/
2
1
/
2
)
z,
(S15)
where
c
[0
,
1] is a real variable parameterizing the classical-bit manifold, and
z
is some (generically mixed) state.
This structure is in agreement with Table 2 in the main text. This suggests that any time after the introduction of
the error, the state has the form
ρ
m
(
τ
q
) =
(
1
/
2
c
1
/
2
c
1
/
2
1
/
2
)
M
(
τ
q
)
,
(S16)
where the classical bit
c
remain perfectly stored in the off-diagonal sector and only
M
changes. If the error is large
enough to move the system to the weak-unbroken phase (e.g.
κ
1
/λ >
2), then the state will evolve toward a unique
steady state and the classical information will be lost. This agrees with our numerical results.
Going back to cat qubit superpositions encoded in the steady state of
L
0
,
i.e.
,
|
ψ
=
c
1
|
α
+
c
0
|−
α
, our analysis
implies that “small” single-photon loss (
κ
1
/λ <
2) can induce phase-flip errors (in the basis above) that are not
passively correctable (while keeping bit-flip errors passively correctable), while large single-photon loss (
κ
1
/λ >
2)
can induce both phase-flip and bit-flip errors that are not passively correctable. Errors that cannot be corrected
passively must be actively corrected, e.g. via redundantly encoding into ancilla qubits.
[S1] T. Prosen, New Journal of Physics
10
, 043026 (2008).
[S2] T. Prosen and T. H. Seligman, Journal of Physics A: Mathematical and Theoretical
43
, 392004 (2010).
[S3] V. V. Albert and L. Jiang, Phys. Rev. A
89
, 022118 (2014).
[S4] E. Knill, R. Laflamme, and L. Viola, Phys. Rev. Lett.
84
, 2525 (2000).
[S5] D. A. Lidar, I. L. Chuang, and K. B. Whaley, Phys. Rev. Lett.
81
, 2594 (1998).
[S6] V. V. Albert, B. Bradlyn, M. Fraas, and L. Jiang, Phys. Rev. X
6
, 041031 (2016).
[S7] M. Mirrahimi, Z. Leghtas, V. V. Albert, S. Touzard, R. J. Schoelkopf, L. Jiang, and M. H. Devoret, New Journal of
Physics
16
, 045014 (2014).
[S8] R. Blume-Kohout, H. K. Ng, D. Poulin, and L. Viola, Phys. Rev. A
82
, 062306 (2010).