## Persistent Topological Negativity in a High-Temperature Mixed-State

Yonna Kim, 1 Ali Lavasani, 2 and Sagar Vijay 1

1 Department of Physics, University of California, Santa Barbara, CA 93106, USA 2 Kavli Institute for Theoretical Physics, University of California, Santa Barbara, CA, 93106

We study the entanglement structure of the Greenberger-Horne-Zeilinger (GHZ) state as it thermalizes under a strongly-symmetric quantum channel describing the Metropolis-Hastings dynamics for the d -dimensional classical Ising model at inverse temperature β . This channel outputs the classical Gibbs state when acting on a product state in the computational basis. When applying this channel to a GHZ state in spatial dimension d &gt; 1, the resulting mixed-state changes character at the Ising phase transition temperature from being long-range entangled to short-range-entangled as temperature increases. Nevertheless, we show that the topological entanglement negativity of a large region is insensitive to this transition and takes the same value as that of the pure GHZ state at any finite temperature β &gt; 0. We establish this result by devising a local operations and classical communication (LOCC) 'decoder' that provides matching lower and upper bounds on the negativity in the thermodynamic limit which may be of independent interest. This perspective connects the negativity to an error-correction problem on the ( d -1)-dimensional bipartitioning surface and explains the persistent negativity in certain correlated noise models found in previous studies. Numerical results confirm our analysis.

Introduction: The nature and fate of quantum entanglement in mixed-states of quantum many-body systems -arising in thermal quantum matter or after a pure state experiences local decoherence - have been recently investigated in order to understand the universal quantum correlations that can persist in open quantum systems [1-18]. Patterns of long-range-entanglement (LRE) in mixed quantum matter are intricately related to quantum many-body phases which can act as quantum error-correcting codes [19-21]. Thermal states of quantum many-body systems are also now known to be efficiently classically-preparable at sufficiently high temperatures [22].

The entanglement negativity provides a well-studied diagnostic of bipartite entanglement in a mixed-state, which vanishes if the state can be prepared from a product state by local operations and classical communication (LOCC) across the bipartition [23]. While the negativity of a region grows with the area of its boundary in a generic Gibbs state [24], contributions to the negativity which are independent to coarse-grained details of the bipartition signal the presence of LRE and can probe the presence of mixed-state topological quantum order [20]. Diagnostics of non-local mixed-state entanglement, such as the topological entanglement negativity [20], are challenging to study analytically, and the subtle behavior of mixed-state LRE under a finite-depth local quantum channel (FDLC) remains to be fully understood [8, 25].

In this work, we investigate the fate of mixed-state quantum entanglement under a Metropolis-Hastings [26, 27] channel describing the equilibration of a thermal classical system, specifically that of a d -dimensional Ising model at inverse temperature β . When acting on a pure state in the computational basis, this channel prepares a classical thermal state of the d -dimensional Ising model.

When the initial state is chosen to be a Green- berger-Horne-Zeilinger (GHZ) state, we show that the negativity of a subregion of the resulting mixed-state ρ β at any finite temperature β &gt; 0 is asymptotically exactly that of the pure GHZ state in the limit that the boundary of the sub-region becomes large. In spatial dimensions d &gt; 1, the negativity is thus a constant, independent of geometric details of the bipartitioning surface and insensitive to the thermal phase transition , as summarized in Eq. (8). Nevertheless, the entanglement properties of the mixed-state do change across this transition: we show that ρ β cannot be expressed as a convex sum of shortrange-entangled (SRE) pure states in the ordered phase but is preparable using an FDLC acting on a product initial state in the disordered phase. This observed behavior is in stark contrast to the negativity of the Gibbs state of a local, quantum many-body Hamiltonian which necessarily vanishes above a finite temperature [22].

This result is established using the fact that the entanglement negativity is an entanglement monotone [23, 28], so that the negativity of a subregion A of any state ρ , denoted N A ( ρ ), cannot increase under an LOCC operation R across the bipartitioning surface,

$$\mathcal { N } _ { A } ( \mathcal { R } [ \rho ] ) \leq \mathcal { N } _ { A } ( \rho ). \text{ \quad \ \ } ( 1 )$$

We explicitly construct an LOCC operation that recovers the GHZ state at any finite temperature and for any bipartition in the asymptotic limit that the boundary of A becomes large. We show that for a contiguous bipartition, this channel is the decoder for a ( d -1)-dimensional repetition code, which the two parties in the LOCC protocol must use to determine the precise unitary circuit that will recover the GHZ state. The well-known success of the ( d -1)-dimensional repetition code up to a maximal error strength in dimensions d &gt; 1 is intimately related to the recovery of the GHZ state, and thus the constant entanglement negativity, at any finite temperature.

We comment on broader utility of this perspective for understanding mixed-state negativity by demonstrating that this can provide a simple understanding of why the topological entanglement negativity can remain a constant in correlated noise models (e.g. in [21]).

Setup: Consider N qubits on a d -dimensional hypercubic lattice, in a state described by the density matrix

$$\rho _ { \beta } \equiv \frac { 1 } { \mathcal { Z } _ { \beta } } \sum _ { \sigma } e ^ { - \beta H } \left | \psi _ { \sigma } \right \rangle \left \langle \psi _ { \sigma } \right |. \quad \quad ( 2 ) \quad \text{fig}$$

Here,

$$H \equiv - J \sum _ { \langle i, j \rangle } Z _ { i } Z _ { j } \quad \quad ( 3 ) \quad \text{$\text{te}$} \\ \text{th}$$

is the energy of a classical Ising model on the d -dimensional hypercubic lattice with ferromagnetic ( J &gt; 0) nearest-neighbor interactions, while Z β ≡ Tr ( e -βH ) is the corresponding partition function at inverse temperature β . Furthermore, the state | ψ σ ⟩ ≡ ( | σ ⟩ + U σ | ⟩ ) / √ 2, where | σ ⟩ ≡ | σ , . . . , σ 1 N ⟩ is a product state in the Pauli Z basis, and U ≡ ∏ j X j is a unitary transformation which generates the Ising symmetry of Eq. (3).

The mixed-state ρ β is naturally obtained by starting with a GHZ state | ψ + ⟩ = ( |↑ · · ·⟩ + |↓ · · ·⟩ ) / √ 2, which is an equal-amplitude superposition of the macroscopicallydistinct all-up and all-down qubit configurations, and then applying a channel Φ β which describes a MetropolisHastings algorithm [26, 27] in which the energy H is measured, and then a unitary operation X j is applied on a random qubit with probability min { 1 exp( , -β ∆ E j ) } , where ∆ E j ≡ E new -E old is the energy difference between the new and old qubit configurations.

Repeated measurements and feedback, starting from any product state in the Pauli Z basis, will eventually produce a Gibbs ensemble of the d -dimensional classical Ising model at inverse temperature β [26, 27]. We specifically define Φ β to be repeated application of the Metropolis-Hastings algorithm so as to produce a thermal steady-state when acting on the all-up state ρ 0 ≡ |↑ · · ·⟩ ⟨↑ · · ·| , i.e.

$$\Phi _ { \beta } [ \rho _ { 0 } ] = \frac { 1 } { \mathcal { Z } _ { \beta } } \sum _ { \sigma } e ^ { - \beta H } \left | \sigma \right \rangle \left \langle \sigma \right |. \quad \quad \left ( 4 \right ) \, \left ( \, \text{sam} \, \\ \text{is st} \,. \,.$$

Each operation that comprises the channel Φ β manifestly commutes with the Ising symmetry transformation U . Therefore, U is a strong symmetry of this channel [29, 30], i.e. U Φ [ β ρ ] = Φ [ β Uρ ] ∀ ρ , so that

$$\Phi _ { \beta } [ | \psi _ { + } \rangle \, \langle \psi _ { + } | ] = \frac { 1 + U } { \sqrt { 2 } } \Phi _ { \beta } [ \rho _ { 0 } ] \frac { 1 + U } { \sqrt { 2 } } = \rho _ { \beta }. \quad ( 5 )$$

Mixed-State Entanglement: The state ρ β witnesses a thermal phase transition with increasing temperature in spatial dimensions d &gt; 1 between a ferromagnetically-ordered and a disordered state. At any

Figure 1: The ABC partitioning used for computing the conditional mutual information I ( A : C B | ).

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00066v2-figure-bce5c379-73fd-4f20-9f92-cb8b73ff0a2d.jpg)

temperature, the state ρ β has a strong Ising symmetry Uρ β = ρ β . However, the ferromagnetically-ordered phase is also characterized by long-range-order in local operators which are charged under this symmetry | i - |→∞ j

̸

Tr( ρ β Z Z i j ) = 0. These conditions are sufficient to establish [9] that this mixed-state cannot be described as an ensemble of SRE pure-states. In Ref.[5] it has been argued that ρ β in the disordered phase can be written as an ensemble of SRE pure states. We prove a potentially stronger result in the Supplemental Material [31], by showing that in the disordered phase, ρ β can be prepared by a FDLC acting on a product state and is thus SRE based on the even more restrictive definition of SRE mixed-states in Refs.[32-34].

Conditional Mutual Information: We now investigate entanglement properties of ρ β and their relation to the Ising phase transition. The conditional mutual information (CMI) is a measure of tripartite correlations in a state and has been used to study various aspects of entanglement in quantum systems [35-40], including probing mixed state quantum phases of matter [25]. Given subsystems A B , , and C , the CMI is defined as

$$\text{assi-} \quad I ( A \colon C | B ) = S ( A B ) + S ( B C ) - S ( B ) - S ( A B C ). \ ( 6 )$$

Consider the partitioning shown in Fig. 1 with r &gt; 0, where C is taken to be the complement of AB . Note that for any proper subset A of qubits, ρ β,A = tr ¯ A ( ρ β ) = ρ Cl β,A , where ρ Cl β = e -βH / Z β is the classical Gibbs state [41]. This shows that the first three terms in Eq. (6) take the same value when computed for ρ β as for ρ Cl β . Moreover, it is straightforward to use Eq.(2) to show S ρ ( β ) = S ρ ( Cl β ) -log(2). Therefore I ( A : C B | ) = log(2) + I Cl ( A : C B | ), where I Cl ( A : C B | ) is the CMI computed for ρ Cl β . On the other hand, since B is assumed to be non-empty, there is no interaction term in H with support on both A and C and hence I Cl ( A : C B | ) = 0 [42], and

$$I ( A \, \colon C | B ) = \log ( 2 ), \text{ \quad \ \ } ( 7 )$$

independent of r , temperature β , and the system size. It is worth mentioning that in contrast, for the thermal Gibbs state of local quantum Hamiltonians, CMI decays exponentially with r [43-45].

&lt;

l

a

t

e

x

i

s

h

1

\_

b

6

4

=

"

R

g

n

q

v

Z

U

+

f

X

H

T

0

J

c

A

&gt;

B

9

V

D

L

S

M

F

3

j

W

Y

C

u

K

N

G

p

E

O

r

Q

/

o

7

8

w

P

d

2

5

m

z

y

k

I

&lt;

l

a

t

e

x

i

s

h

1

\_

b

6

4

=

"

/

r

y

G

S

D

j

P

7

5

U

R

X

o

B

8

9

3

Q

&gt;

A

n

c

V

L

T

g

J

E

O

z

F

+

I

p

v

q

Y

Z

H

w

d

M

2

C

W

k

u

f

N

K

0

m

&lt;

l

a

t

e

x

i

s

h

1

\_

b

6

4

=

"

E

9

r

z

H

M

N

P

5

S

f

V

Z

q

D

W

j

k

Q

p

O

c

&gt;

A

B

X

L

g

F

3

Y

C

u

K

G

R

J

U

/

o

7

8

w

n

+

d

2

0

T

m

v

y

I

&lt;

l

a

t

e

x

i

s

h

1

\_

b

6

4

=

"

S

Y

C

X

q

H

D

o

L

k

M

A

9

2

8

Q

5

+

U

&gt;

B

n

c

V

T

g

J

E

O

z

F

I

p

r

v

Z

w

d

0

/

P

G

W

3

u

7

m

N

f

j

y

K

R

&lt;

l

a

t

e

x

i

s

h

1

\_

b

6

4

=

"

Y

Z

p

d

y

B

v

U

w

3

z

J

L

5

A

u

c

P

N

&gt;

n

V

D

S

g

E

O

G

M

r

H

C

W

j

Q

9

I

o

T

2

m

k

8

+

f

/

F

R

0

7

q

X

K

&lt;

l

a

t

e

x

i

s

h

1

\_

b

6

4

=

"

P

m

Z

w

8

0

C

T

d

5

p

k

S

r

O

J

7

o

&gt;

A

B

9

X

c

V

D

L

g

I

F

3

/

G

u

n

j

Q

N

y

U

+

2

f

R

E

z

q

W

K

M

Y

H

v

&lt;

l

a

t

e

x

i

s

h

1

\_

b

6

4

=

"

R

g

n

q

v

Z

U

+

f

X

H

T

0

J

c

A

&gt;

B

9

V

D

L

S

M

F

3

j

W

Y

C

u

K

N

G

p

E

O

r

Q

/

o

7

8

w

P

d

2

5

m

z

y

k

I

&lt;

l

a

t

e

x

i

s

h

1

\_

b

6

4

=

"

R

g

n

q

v

Z

U

+

f

X

H

T

0

J

c

A

&gt;

B

9

V

D

L

S

M

F

3

j

W

Y

C

u

K

N

G

p

E

O

r

Q

/

o

7

8

w

P

d

2

5

m

z

y

k

I

&lt;

l

a

t

e

x

i

s

h

1

\_

b

6

4

=

"

/

r

y

G

S

D

j

P

7

5

U

R

X

o

B

8

9

3

Q

&gt;

A

n

c

V

L

T

g

J

E

O

z

F

+

I

p

v

q

Y

Z

H

w

d

M

2

C

W

k

u

f

N

K

0

m

&lt;

l

a

t

e

x

i

s

h

1

\_

b

6

4

=

"

/

r

y

G

S

D

j

P

7

5

U

R

X

o

B

8

9

3

Q

&gt;

A

n

c

V

L

T

g

J

E

O

z

F

+

I

p

v

q

Y

Z

H

w

d

M

2

C

W

k

u

f

N

K

0

m

&lt;

l

a

t

e

x

i

s

h

1

\_

b

6

4

=

"

Q

p

U

S

8

v

F

q

f

2

+

r

u

g

P

O

j

E

&gt;

A

C

D

X

c

V

N

M

w

H

3

n

5

Y

9

K

d

o

R

W

k

J

m

G

/

B

7

y

L

Z

I

T

0

&lt;

l

a

t

e

x

i

s

h

1

\_

b

6

4

=

"

p

D

g

/

9

z

E

u

v

H

f

I

S

k

K

F

Z

J

X

w

&gt;

A

C

c

V

N

8

M

3

n

5

Y

d

O

U

0

L

r

Q

q

T

B

P

j

R

m

W

2

o

y

G

7

+

&lt;

l

a

t

e

x

i

s

h

1

\_

b

6

4

=

"

E

9

r

z

H

M

N

P

5

S

f

V

Z

q

D

W

j

k

Q

p

O

c

&gt;

A

B

X

L

g

F

3

Y

C

u

K

G

R

J

U

/

o

7

8

w

n

+

d

2

0

T

m

v

y

I

&lt;

l

a

t

e

x

i

s

h

1

\_

b

6

4

=

"

S

Y

C

X

q

H

D

o

L

k

M

A

9

2

8

Q

5

+

U

&gt;

B

n

c

V

T

g

J

E

O

z

F

I

p

r

v

Z

w

d

0

/

P

G

W

3

u

7

m

N

f

j

y

K

R

&lt;

l

a

t

e

x

i

s

h

1

\_

b

6

4

=

"

E

9

r

z

H

M

N

P

5

S

f

V

Z

q

D

W

j

k

Q

p

O

c

&gt;

A

B

X

L

g

F

3

Y

C

u

K

G

R

J

U

/

o

7

8

w

n

+

d

2

0

T

m

v

y

I

&lt;

l

a

t

e

x

i

s

h

1

\_

b

6

4

=

"

S

Y

C

X

q

H

D

o

L

k

M

A

9

2

8

Q

5

+

U

&gt;

B

n

c

V

T

g

J

E

O

z

F

I

p

r

v

Z

w

d

0

/

P

G

W

3

u

7

m

N

f

j

y

K

R

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00066v2-figure-263c3db1-96f8-4ff0-a952-e262d76a9a71.jpg)

-

Figure 2: Measurements of domain walls that lie exclusively within A or A lead to configurations such as the one shown in the left panel, with domain walls indicated by red lines on the dual lattice. There are two possibilities for how these domain wall connect across this interface, corresponding to the choices σ σ 1 1 = ± 1.

Entanglement Negativity: We now investigate the entanglement negativity of a subregion A of the mixedstate (2), N A ( ρ β ) ≡ || ( ρ T A β || 1 -1) / 2. We will show that for a region A ,

$$\stackrel { \cdot } { \mathcal { N } } _ { A } ( \rho _ { \beta } ) = \begin{cases} \begin{array} { c c c } 1 \\ \frac { 1 } { 2 } - O ( \exp ( - | \partial A | ) ) & \beta > 0 \\ 0 & \beta = 0 \end{array} \end{cases}$$

where | ∂A | is the number of qubits on the boundary of A . The second result is easily established. At infinite temperature, ρ β =0 = (1 + U / ) 2 N , so that ρ T A β = ρ β and N A ( ρ β =0 ) = 0.

The first result is determined as follows. The negativity is a convex function [23, 28] so that

$$\mathcal { N } _ { A } ( \rho _ { \beta } ) & \leq \frac { 1 } { \mathcal { Z } _ { \beta } } \sum _ { \sigma } e ^ { \beta J \sum _ { \langle i, j \rangle } \sigma _ { i } \sigma _ { j } } \mathcal { N } _ { A } ( | \psi _ { \sigma } \rangle \langle \psi _ { \sigma } | ) = \frac { 1 } { 2 } \quad ( 9 ) \quad \text{$\text{perat}$} \\ \hat { \ }. \quad \hat { \ }. \quad \hat { \ }. \quad \hat { \ }. \quad \hat { \ }. \quad \hat { \ }. \quad \sim \hat { \ }. \quad \hat { \ }. \quad \hat { \ }. \quad \hat { \ }. \quad \hat { \ }. \quad \hat { \ }. \quad \hat { \ }. \quad \hat { \ }. \quad \hat { \ }.$$

for any region A , since | ψ σ ⟩ is related to the GHZ state | ψ + ⟩ by single qubit gates, so that N A ( | ψ σ ⟩⟨ ψ σ | ) = N A ( | ψ + ⟩⟨ ψ + | ). The latter is exactly 1 / 2 for any region A . The negativity is also an entanglement monotone [23, 28] (see Eq. (1)). Here, we show that given the state ρ β for any β &gt; 0, there is an LOCC operation R with respect to any region A and its complement ¯ that recovers the A GHZ state with a fidelity that approaches unity exponentially quickly in the size of the boundary of A . Combining this with (9), allows us to establish Eq. (8).

We construct the LOCC channel R for a contiguous bipartition before generalizing this construction for a non-contiguous subregion A . The state ρ β describes the ensemble of pure-states | ψ σ ⟩ drawn with probability e βJ ∑ ⟨ i,j ⟩ σ σ i j / Z β . Given a state from this ensemble, we measure the operators Z Z i j along all nearest-neighbor bonds lying exclusively within A and A . The resulting pattern of domain walls (see e.g. Fig. 2) - is consistent with having drawn the state (i) | ψ σ ⟩ or the state (ii) ∏ j ∈ A X j | ψ σ ⟩ and leaves one classical bit's worth of ambiguity in determining the precise state.

The preparation of the GHZ state thus maps onto a decoding problem in a classical ( d -1)-dimensional repetition code living on the interface between A and A . The classical bit of information is the relative alignment of any two qubits across the bipartition in the Pauli Z basis (e.g. τ 1 ≡ σ σ 1 1 in Fig. 2). The measured 'syndromes' of the repetition code (e.g. τ 1 τ 2 = ( σ σ 1 1 )( σ σ 2 2 )) and the Boltzmann weights for the d -dimensional Ising model are used to perform maximum-likelihood decoding to guess the correct state. Let w 1 and w 2 denote the Boltzmann weights corresponding to states ( ) and ( i ii ) respectively. We guess that the state we have is ( ) with probability i p ≡ w / w 1 ( 1 + w 2 ) and ( ii ) with probability 1 -p . We subsequently apply the appropriate single-site unitary gates to turn this into the GHZ state. The LOCC channel is thus

- 1. Measure Z Z i j along all nearest-neighbor bonds lying exclusively within A or A .
- 2. Let | σ ⟩ be a product state in the Pauli Z basis which is consistent with the measurement outcomes. Apply the unitary U 1 = ∏ j X (1 -σ j ) / 2 j with probability p , or the unitary U 2 = ∏ j ∈ A X j ∏ j X (1 -σ j ) / 2 j with probability 1 -p .

At any finite temperature, the relative probability p/ (1 -p ) = w /w 1 2 scales exponentially in the size of the boundary of A , due to the fact that the energy of an Ising domain wall is extensive in its length. Thus, the decoder for the repetition code should be successful at any finite temperature and fail exactly at infinite temperature when w 1 = w 2 = 1 2. / This is consistent with the known behavior of the repetition code with singlequbit bit-flip errors with probability q , for which perfect decoding is possible for any q ≤ 1 / 2.

We may generalize this channel for the case where A is a non-contiguous region. Step 1 (above) is replaced by measurements of Z Z i j across all pairs of sites entirely within A and entirely within ¯ . A As before, only one classical bit's worth of ambiguity remains after these measurements. The quantum circuit that takes this state onto the GHZ state is then determined by a decoder for a repetition code with | ∂A | sites. Thus, decoding the correct quantum circuit should again be successful at any β &gt; 0 as | ∂A | becomes large, regardless of the geometry of the bipartition.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00066v2-figure-20a659b6-92d7-4036-ad74-dbbe347176e0.jpg)

(a)

Figure 3: a) Entanglement negativity N A versus temperature in the state ρ β in d = 2, where subregion A is taken to be half the system. The inset shows 1 / 2 -N A at T = 5 versus L in a semi-log plot, which shows that N A saturates exponentially quickly in the size of the boundary. b) Derivative of entanglement negativity of a single site with respect to temperature, showing a singularity at the Ising phase transition temperature.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00066v2-figure-a9f44437-33f2-4a4c-8f4b-4f8bb6e56db9.jpg)

For any bipartition, we show that the fidelity of the recovered state R ρ ( β ) with the GHZ state is [31]

$$\mathcal { F } \left ( \mathcal { R } ( \rho _ { \beta } ), \left | \psi _ { + } \right \rangle \left \langle \psi _ { + } \right | \right ) = \frac { 1 } { 2 } \left [ 1 - \left \langle \tanh ( \beta H _ { \partial } ) \right \rangle \right ] \quad ( 1 0 ) \quad \text{ aries,}$$

where H ∂ denotes terms in the Ising Hamiltonian acting exclusively along bonds crossing the bipartition between A and A and the expectation value is taken in the thermal state for the classical Ising model. At any finite temperature, the energy H ∂ is negative, and extensively large in the boundary of A , so that 〈 tanh( βH ∂ ) 〉 → 1 as A becomes large.

Numerical Results: To numerically confirm this result, we obtain an alternate expression for the entanglement negativity of any bi-partition

$$\mathcal { N } _ { A } [ \rho _ { \beta } ] = \frac { 1 } { 2 } \left \langle \left | \tanh ( \beta H _ { \partial } [ \sigma ] ) \right | \right \rangle _ { q } \quad \ \ ( 1 1 ) \quad \ \text{only} \quad \ \text{only}$$

Figure 4: Dephasing of the toric code by the local operators Z j + δ X j at each lattice site j as studied in Ref. [21]. Correlations in the syndromes permit perfect preparation of a mixed-state within the ground-space of the toric code using LOCC across A and A .

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00066v2-figure-31bfe7a9-171c-4a2c-97c8-0d9a1c941ed1.jpg)

where the final expectation value is taken with respect to a distribution q σ [ ], which is derived from the Boltzmann distribution for the Ising model, and is given in the Supplemental Material [31]. We note that a similar result has been found before for a closely related set of mixed states [46].

We present the numerical results of classical Monte Carlo simulations for a 2D system of qubits on an L × L square lattice with periodic boundary conditions. Fig. 3a shows the negativity N A ( ρ β ) when A is taken to be the L × L 2 cylinder. As is clear from Fig. 3a, in the thermodynamic limit, the negativity is equal to 1 / 2 even above the Ising critical temperature (marked by the dashed line in Fig.3a), as expected from the expression in (43). The fact that the negativity of subregions with extensive boundary is 1 / 2 irrespective of temperature, means that if one defines the topological negativity as the constant contribution to the negativity as in Ref. [2], one would find that the resulting mixed state has N topo = -1 / 2 even above the critical temperature where ρ β is an SRE mixed-state. [5, 6, 32, 47, 48] (cf. Ref. [20]).

Although the thermal phase transition at T c has no effect on the negativity of subregions with extensive boundaries, the negativity of finite subregions (while continuous) becomes singular at T = T c . Fig. 3b shows the derivative of the negativity of a single site with respect to temperature which shows a singularity at T = T c , which is consistent with the singularity of the heat capacity at the d = 2 classical Ising phase transition.

Discussion and Outlook: Using LOCC to recover states with a known pattern of LRE may provide another perspective that can be used to investigate mixed-state LRE under decoherence. For correlated decoherence, this is apparent in specific examples. The two-dimensional toric code dephased by local operators Z j + δ X j at each site j , with the vector δ given in Fig. 4, was studied in Ref. [21]. Quantum state trajectories under this channel only contain patterns of Z 2 charge and flux which are separated by δ . Therefore, for any dephasing strength, it

is possible to use LOCC across the bipartition in Fig. 4 to deterministically recover the toric code ground-state. The Z 2 charge and flux lying on the bipartitioning surface, which cannot be directly measured in the LOCC protocol, can be perfectly inferred from stabilizer measurements in A and A due to the noise correlations. The negativity for this bipartition and for any dephasing strength [21] is thus exactly that of the toric-code ground-state [49]. Another example for which the LOCC perspective presented here could be insightful is the study of entanglement structure in the toric code with singlesite Y decoherence [50, 51]. In this case, while the noise model itself is not correlated, the LOCC protocol can benefit from the correlations that are present in the syndromes due to the extensively many strong symmetries of this channel [50, 52]. The detailed study of this model and other correlated noise models is left to a future work.

We identify other directions implied by our results. First, it is interesting to consider the fate of mixedstate entanglement under thermalization by more generic quantum channels which are strongly-symmetric. More concretely, we may consider thermalizing the GHZ state via a strongly symmetric implementation of the quantum Metropolis-Hastings dynamics [53, 54] for the d -dimensional transverse field quantum Ising model. The resulting mixed state would be the projection of the quantum Gibbs state onto the symmetric sector. A similar LOCC protocol to recover the GHZ state can be used to relate the entanglement negativity to the performance of a ( d -1)-dimensional repetition code under a coherent correlated noise model. Energetic considerations can be used to argue about the success of the LOCC protocol in turning the resulting mixed state into the GHZ state. This lower bounds the negativity of the resulting mixedstate and suggests that the negativity is non-vanishing at high temperatures, though this remains to be confirmed.

Second, while the persistence of topological negativity at high temperatures in our work can be understood as a consequence of the repetition code possessing maximum error threshold, it is interesting to search for settings where a similar line of argument would lead to a different code with a finite threshold. Such a system would exhibit a finite temperature phase transition in negativity, which may not necessarily coincide with the standard order-disorder thermal phase transition.

Lastly, given that the topological negativity is unable to detect the phase transition in the complexity of the mixed-state ρ β (from LRE to SRE), it is interesting to look for other quantities that would be able to capture the complexity of the quantum phase transition in our system as well as in more general settings. A possible candidate would be the convex-roof extension of the conditional mutual information as introduced in Ref. [55] using the same ABC partitioning used above.

Acknowledgments: We thank Yimu Bao, Matthew Fisher, Tarun Grover, Tibor Rakovszky, Subhayan Sahu,

Shengqi Sang and Yaodong Li for useful discussions. This research was supported in part by grant NSF PHY2309135 to the Kavli Institute for Theoretical Physics (KITP). Use was made of computational facilities purchased with funds from the National Science Foundation (CNS-1725797) and administered by the Center for Scientific Computing (CSC). The CSC is supported by the California NanoSystems Institute and the Materials Research Science and Engineering Center (MRSEC; NSF DMR 2308708) at UC Santa Barbara. SV is partly supported as an Alfred P. Sloan Research Fellow.

- [1] K.-H. Wu, T.-C. Lu, C.-M. Chung, Y.-J. Kao, and T. Grover, Entanglement renyi negativity across a finite temperature transition: a monte carlo study, Physical Review Letters 125 , 140603 (2020).
- [2] T.-C. Lu and T. Grover, Structure of quantum entanglement at a finite temperature critical point, Physical Review Research 2 , 043345 (2020).
- [3] T.-C. Lu and S. Vijay, Characterizing long-range entanglement in a mixed state through an emergent order on the entangling surface, Physical Review Research 5 , 033031 (2023).
- [4] R. Fan, Y. Bao, E. Altman, and A. Vishwanath, Diagnostics of mixed-state topological order and breakdown of quantum memory, PRX Quantum 5 , 020343 (2024).
- [5] Y.-H. Chen and T. Grover, Symmetry-enforced many-body separability transitions, arXiv preprint arXiv:2310.07286 (2023).
- [6] Y.-H. Chen and T. Grover, Separability transitions in topological states induced by local decoherence, Physical Review Letters 132 , 170602 (2024).
- [7] Y.-H. Chen and T. Grover, Unconventional topological mixed-state transition and critical phase induced by selfdual coherent errors, arXiv preprint arXiv:2403.06553 (2024).
- [8] S. Sang, Y. Zou, and T. H. Hsieh, Mixed-state quantum phases: Renormalization and quantum error correction, arXiv preprint arXiv:2310.08639 (2023).
- [9] T.-C. Lu, Z. Zhang, S. Vijay, and T. H. Hsieh, Mixedstate long-range order and criticality from measurement and feedback, PRX Quantum 4 , 030318 (2023).
- [10] L. A. Lessa, M. Cheng, and C. Wang, Mixed-state quantum anomaly and multipartite entanglement, arXiv preprint arXiv:2401.17357 (2024).
- [11] Z. Zhang, U. Agrawal, and S. Vijay, Quantum communication and mixed-state order in decohered symmetry-protected topological states, arXiv preprint arXiv:2405.05965 (2024).
- [12] J. Y. Lee, C.-M. Jian, and C. Xu, Quantum criticality under decoherence or weak measurement, PRX quantum 4 , 030317 (2023).
- [13] C. de Groot, A. Turzillo, and N. Schuch, Symmetry protected topological order in open quantum systems, Quantum 6 , 856 (2022).
- [14] G.-Y. Zhu, N. Tantivasadakarn, A. Vishwanath, S. Trebst, and R. Verresen, Nishimori's cat: stable longrange entanglement from finite-depth unitaries and weak measurements, Physical Review Letters 131 , 200201

(2023).

- [15] J. Y. Lee, Y.-Z. You, and C. Xu, Symmetry protected topological phases under decoherence, arXiv preprint arXiv:2210.16323 (2022).
- [16] Z. Wang, X.-D. Dai, H.-R. Wang, and Z. Wang, Topologically ordered steady states in open quantum systems, arXiv preprint arXiv:2306.12482 (2023).
- [17] Z. Wang and L. Li, Anomaly in open quantum systems and its implications on mixed-state quantum phases, arXiv preprint arXiv:2403.14533 (2024).
- [18] J. Y. Lee, W. Ji, Z. Bi, and M. Fisher, Decoding measurement-prepared quantum phases and transitions: from ising model to gauge theory, and beyond, arXiv preprint arXiv:2208.11699 (2022).
- [19] Y. Bao, R. Fan, A. Vishwanath, and E. Altman, Mixedstate topological order and the errorfield double formulation of decoherence-induced transitions, arXiv preprint arXiv:2301.05687 (2023).
- [20] T.-C. Lu, T. H. Hsieh, and T. Grover, Detecting topological order at finite temperature using entanglement negativity, Physical Review Letters 125 , 116801 (2020).
- [21] Z. Wang, Z. Wu, and Z. Wang, Intrinsic mixedstate topological order without quantum memory, arXiv preprint arXiv:2307.13758 (2023).
- [22] A. Bakshi, A. Liu, A. Moitra, and E. Tang, Hightemperature gibbs states are unentangled and efficiently preparable, arXiv preprint arXiv:2403.16850 (2024).
- [23] G. Vidal and R. F. Werner, Computable measure of entanglement, Physical Review A 65 , 032314 (2002).
- [24] N. E. Sherman, T. Devakul, M. B. Hastings, and R. R. Singh, Nonzero-temperature entanglement negativity of quantum spin models: Area law, linked cluster expansions, and sudden death, Physical Review E 93 , 022128 (2016).
- [25] S. Sang and T. H. Hsieh, Stability of mixed-state quantum phases via finite markov length, arXiv preprint arXiv:2404.07251 (2024).
- [26] N. Metropolis, A. W. Rosenbluth, M. N. Rosenbluth, A. H. Teller, and E. Teller, Equation of state calculations by fast computing machines, The journal of chemical physics 21 , 1087 (1953).
- [27] W. K. Hastings, Monte carlo sampling methods using markov chains and their applications, (1970).
- [28] M. B. Plenio, Logarithmic negativity: a full entanglement monotone that is not convex, Physical review letters 95 , 090503 (2005).
- [29] V. V. Albert and L. Jiang, Symmetries and conserved quantities in lindblad master equations, Physical Review A 89 , 022118 (2014).
- [30] S. Lieu, R. Belyansky, J. T. Young, R. Lundgren, V. V. Albert, and A. V. Gorshkov, Symmetry breaking and error correction in open quantum systems, Physical Review Letters 125 , 240405 (2020).
- [31] See supplemental material.
- [32] R. Ma and C. Wang, Average symmetry-protected topological phases, Physical Review X 13 , 031016 (2023).
- [33] R. Ma, J.-H. Zhang, Z. Bi, M. Cheng, and C. Wang, Topological phases with average symmetries: the decohered, the disordered, and the intrinsic, arXiv preprint arXiv:2305.16399 (2023).
- [34] A. Anshu and C. Nirkhe, Circuit lower bounds for low-energy states of quantum code hamiltonians, arXiv preprint arXiv:2011.02044 (2020).
- [35] M. Levin and X.-G. Wen, Detecting topological order in
- a ground state wave function, Physical review letters 96 , 110405 (2006).
- [36] K. Kato, F. Furrer, and M. Murao, Informationtheoretical analysis of topological entanglement entropy and multipartite correlations, Physical Review A 93 , 022317 (2016).
- [37] D. Ding, P. Hayden, and M. Walter, Conditional mutual information of bipartite unitaries and scrambling, Journal of High Energy Physics 2016 , 1 (2016).
- [38] O. Fawzi and R. Renner, Quantum conditional mutual information and approximate markov chains, Communications in Mathematical Physics 340 , 575 (2015).
- [39] S.-u. Lee, C. Oh, Y. Wong, S. Chen, and L. Jiang, Universal spreading of conditional mutual information in noisy random circuits, arXiv preprint arXiv:2402.18548 (2024).
- [40] Y. Zhang and S. Gopalakrishnan, Nonlocal growth of quantum conditional mutual information under decoherence, arXiv preprint arXiv:2402.03439 (2024).
- [41] This in particular shows that for mutual information we have I ( A : C ) = I Cl ( A : C ) provided that r &gt; 0, where I Cl is computed with respect to ρ Cl β .
- [42] P. Clifford and J. Hammersley, Markov fields on finite graphs and lattices, (1971).
- [43] K. Kato and F. G. Brandao, Quantum approximate markov chains are thermal, Communications in Mathematical Physics 370 , 117 (2019).
- [44] T. Kuwahara, K. Kato, and F. G. Brand˜o, Clustering of a conditional mutual information for quantum gibbs states above a threshold temperature, Physical review letters 124 , 220601 (2020).
- [45] T. Kuwahara, Clustering of conditional mutual information and quantum markov structure at arbitrary temperatures, arXiv preprint arXiv:2407.05835 (2024).
- [46] T. Grover and T.-C. Lu, Quantum correlations at finite temperature transitions (2018), conference talk, UC Riverside.
- [47] R. F. Werner, Quantum states with einstein-podolskyrosen correlations admitting a hidden-variable model, Physical Review A 40 , 4277 (1989).
- [48] M. B. Hastings, Topological order at nonzero temperature, Physical review letters 107 , 210501 (2011).
- [49] C. Castelnovo, Negativity and topological order in the toric code, Physical Review A-Atomic, Molecular, and Optical Physics 88 , 042319 (2013).
- [50] D. K. Tuckett, A. S. Darmawan, C. T. Chubb, S. Bravyi, S. D. Bartlett, and S. T. Flammia, Tailoring surface codes for highly biased noise, Physical Review X 9 , 041031 (2019).
- [51] T. Ellison and M. Cheng, Towards a classification of mixed-state topological orders in two dimensions, arXiv preprint arXiv:2405.02390 (2024).
- [52] J. Hauser, Y. Bao, S. Sang, A. Lavasani, U. Agrawal, and M. Fisher, Information dynamics in decohered quantum memory with repeated syndrome measurements: a dual approach, arXiv preprint arXiv:2407.07882 (2024).
- [53] K. Temme, T. J. Osborne, K. G. Vollbrecht, D. Poulin, and F. Verstraete, Quantum metropolis sampling, Nature 471 , 87 (2011).
- [54] M.-H. Yung and A. Aspuru-Guzik, A quantum-quantum metropolis algorithm, Proceedings of the National Academy of Sciences 109 , 754 (2012).
- [55] T.-T. Wang, M. Song, Z. Y. Meng, and T. Grover, An analog of topological entanglement entropy for mixed

states, arXiv preprint arXiv:2407.20500 (2024).

- [56] Another way to derive this pre-factor is to observe that for a given configuration of µ j satisfying the constraint
- ∏ j µ j = 1, there are 2 N b / 2 N -1 different configurations of σ ij such that µ j = A σ j [ ].
- [57] T. Cacoullos, Exercises in probability (Springer Science &amp; Business Media, 2012).

## Supplemental Material

## Entanglement of ρ β

Here, we show that the state (Eq. 2) is short-range-entangled (SRE) throughout the disordered phase at sufficiently high temperatures, and in any spatial dimension. First, observe that this state may be written as

$$\rho _ { \beta } = \sum _ { \mu } p _ { \mu } \left | \phi _ { \mu } \right \rangle \left \langle \phi _ { \mu } \right |$$

where | ϕ µ ⟩ ≡ ρ 1 / 2 β | µ / ⟩ √ p µ and p µ ≡ ⟨ µ ρ | β | µ ⟩ , while | µ ⟩ is a product state in the Pauli X basis. Because of the strong symmetry of this state, Tr( ρ β U ) = 1, where U ≡ ∏ j X j . As a result, each pure state appearing in (12) must satisfy U ϕ | µ ⟩ = | ϕ µ ⟩ .

Furthermore, observe that

$$\rho _ { \beta } ^ { 1 / 2 } & = \frac { 1 } { \sqrt { Z _ { \beta } } } \sum _ { \sigma } e ^ { - \beta H / 2 } \left | \psi _ { \sigma } \right \rangle \left \langle \psi _ { \sigma } \right |.$$

It is evident from this expression that [ Z Z , ρ i j 1 / 2 β ] = 0 for any i , j . As a result, any state | ϕ µ ⟩ appearing in (12) can be obtained by applying single-site unitary gates on the reference state | ϕ + ⟩ = ρ 1 / 2 β | + · · · + ⟩ / √ p + , where | + · · · + is ⟩ a product state in which X j = +1 for each qubit:

$$| \phi _ { \mu } \rangle = \prod _ { j } Z _ { j } ^ { \frac { 1 - \mu _ { j } } { 2 } } \, | \phi _ { + } \rangle \,.$$

We note that p µ = ⟨ µ ρ | β | µ ⟩ = + ⟨ · · · + | ∏ j Z 1 -µ j 2 j ρ β ∏ j Z 1 -µ j 2 j | + · · · + = + ⟩ ⟨ · · · + | ρ β | + · · · + = ⟩ p + . As a result, p µ = 1 2 / N -1 where N is the total number of qubits. To summarize, we may write (12) as

$$\rho _ { \beta } = \frac { 1 } { 2 ^ { N - 1 } } \sum _ { \mu \text{ s.t. } \prod _ { j } \mu _ { j } = + 1 } \prod _ { j } Z _ { j } ^ { \frac { 1 - \mu _ { j } } { 2 } } \left | \phi _ { + } \right \rangle \left \langle \phi _ { + } \right | \prod _ { j } Z _ { j } ^ { \frac { 1 - \mu _ { j } } { 2 } }.$$

We now re-cast the summation by introducing Ising degrees of freedom σ ij along bonds of the hypercubic lattice, and letting µ j ≡ A σ j [ ] ≡ ∏ i ∈ NN j ( ) σ ij , where the product is taken over sites that are nearest-neighbor to j . With periodic boundary conditions, ∏ j A σ j [ ] = +1, and as a result, the constraint on the summation ∏ j µ j = +1 is naturally satisfied by this replacement. With this replacement, we may write

$$\rho _ { \beta } = \frac { 1 } { 2 ^ { N _ { b } } } \sum _ { \sigma } \prod _ { j } Z _ { j } ^ { \frac { 1 - A _ { j } [ \sigma ] } { 2 } } \left | \phi _ { + } \right \rangle \left \langle \phi _ { + } \right | \prod _ { j } Z _ { j } ^ { \frac { 1 - A _ { j } [ \sigma ] } { 2 } }.$$

Here, N b = d · N is the total number of bonds on the d -dimensional hypercubic lattice. The pre-factor ensures the normalization of the density matrix [56]. Using the fact that Z 1 -Aj σ [ ] 2 j = ∏ i ∈ NN( ) j Z 1 -σ ij 2 j , we may write ∏ j Z 1 -Aj σ [ ] 2 j = ∏ ⟨ i,j ⟩ ( Z Z i j ) 1 -σ ij 2 , where the product in the second expression is taken over distinct, nearest-neighbor bonds. As a result,

$$\rho _ { \beta } = \frac { 1 } { 2 ^ { N _ { b } } } \sum _ { \sigma } \prod _ { \langle i, j \rangle } ( Z _ { i } Z _ { j } ) ^ { \frac { 1 - \sigma _ { i j } } { 2 } } \left | \phi _ { + } \rangle \left \langle \phi _ { + } \right | \prod _ { \langle i, j \rangle } ( Z _ { i } Z _ { j } ) ^ { \frac { 1 - \sigma _ { i j } } { 2 } }.$$

We conclude that ρ β may be obtained by dephasing Z Z i j on each bond of the lattice, with maximal strength, a manifestly finite-depth, local quantum channel.

Finally, we may argue that the state | ϕ + ⟩ is SRE when the system is in a disordered phase. We do this by first writing

$$| \phi _ { + } \rangle = \sqrt { 2 ^ { N - 1 } } \rho _ { \beta } ^ { 1 / 2 } \, | + \cdots + \rangle = \sqrt { \frac { 2 ^ { N - 1 } } { Z _ { \beta } } } \sum _ { \sigma } e ^ { - \beta H / 2 } \, | \psi _ { \sigma } \rangle \, \langle \psi _ { \sigma } | + \cdots + \rangle = \frac { 1 } { \sqrt { Z _ { \beta } } } \sum _ { \sigma } e ^ { - \beta H / 2 } \, | \sigma \rangle \,.$$

We now observe that (18) is the ground-state of the Hamiltonian

$$H = - \sum _ { j } Q _ { j }$$

where

$$Q _ { j } \equiv - X _ { j } + \prod _ { i \in \mathrm N N ( j ) } e ^ { - \beta J Z _ { i } Z _ { j } }.$$

## LOCC Operations and Recovering the GHZ State

Here we will show the existence of an LOCC protocol R with respect to bi-partition A and ¯ , A which transforms state ρ β back to the GHZ state in the thermodynamic limit. The LOCC protocol consists of two steps:

- 1. Measure Z Z i j on each link that is entirely in A or in ¯ subregions, finding the outcome A µ ij = ± 1. This does not include the links across the boundary, ∂A .
- 2. For any given set of measurement outcomes µ = { µ i,j } , there exist two possible states which have domain wall configuration consistent with µ . We pick one of the two states based on a coin toss weighted with their respective Boltzmann weight in ρ β (Eq.(2)) and apply a product of Pauli X i operators to eliminate the domain walls in that assumed configuration.

In the following, we will derive Eq.(10) for the fidelity between the GHZ state, | ψ + ⟩ , and the output state R ( ρ β ).

## Fidelity

Let { M µ } be the projective measurements that measure whether there is a domain wall or not on the links which are entirely in region A or ¯ , A

$$M _ { \mu } \equiv \prod _ { < i, j > \in A, \bar { A } } \frac { 1 + \mu _ { i, j } Z _ { i } Z _ { j } } { 2 },$$

where µ ij is the measurement outcome of measuring Z Z i j and the product only includes nearest-neighbor sites which are both in A or both in ¯ . A Note that,

$$M _ { \mu } \left | \psi _ { \sigma } \right \rangle = \prod _ { < i, j > \in A, \bar { A } } \delta _ { \mu _ { i, j }, \sigma _ { i } \sigma _ { j } } \left | \psi _ { \sigma } \right \rangle.$$

A set of measurement outcomes µ specifies a certain domain wall configuration in A and ¯ , A with some domain walls ending on the boundary ∂A . Let 1 ∈ A and ¯ 1 ∈ ¯ be two adjacent sites across the boundary and define A τ = ± 1 to be a marker for whether or not there is a domain wall across the edge connecting 1 and ¯ , i.e. 1 τ ≡ σ σ 1 ¯ 1 (see Fig. 2). Note that µ and τ together uniquely specify the global domain wall configuration, because there should be even number of domain walls across each plaquette of the square lattice. In particular if we label the sites along the boundary as shown in Fig.5, we have,

$$\sigma _ { i } \sigma _ { i } = \tau \prod _ { j = 1 } ^ { i - 1 } \mu _ { j, j + 1 } \mu _ { j, j + 1 }.$$

Let X µ,τ be the appropriate product of single qubit PauliX s that corresponds to the domain wall configuration specified by µ and τ . In particular, if M µ | ψ σ ⟩ = | ψ σ ⟩ , we have

$$X _ { \mu, \tau } \left | \psi _ { \sigma } \right \rangle \equiv \begin{cases} \left | \psi _ { + } \right \rangle & \text{if $\sigma_{1}\sigma_{\bar{1}}=\tau$} \\ \prod _ { j \in A } X _ { j } \left | \psi _ { + } \right \rangle & \text{if $\sigma_{1}\sigma_{\bar{1}}\neq\tau$}. \end{cases}$$

̸

Figure 5: Labeling the boundary sites in A and ¯ A

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00066v2-figure-bbed4101-bec4-4fde-a69d-a957a7553447.jpg)

Given the measurement outcomes µ , the decoder R applies X µ,τ with τ = ± 1 chosen randomly with probability P µ,τ . The probability P µ,τ is chosen to be proportional to the Boltzmann weight of the spin configuration specified by µ and τ and according to Eq.(23) can be written as,

$$P _ { \mu, \tau } = \frac { e ^ { \beta J \tau } \sum _ { i \geq 1 } \prod _ { j = 1 } ^ { i - 1 } \mu _ { j, j + 1 } \mu _ { j, j + 1 } } { \sum _ { \tau ^ { \prime } = \pm 1 } e ^ { \beta J \tau ^ { \prime } } \sum _ { i \geq 1 } \prod _ { j = 1 } ^ { i - 1 } \mu _ { j, j + 1 } \mu _ { j, j + 1 } }.$$

Putting everything together, we can compute the average 'decoding' fidelity,

$$F = \langle \psi _ { + } | \mathcal { R } ( \rho _ { \beta } ) | \psi _ { + } \rangle = \sum _ { \sigma, \mu, \tau } \frac { e ^ { - \beta H [ \sigma ] } } { \mathcal { Z } _ { \beta } } P _ { \mu, \tau } | \langle \psi _ { + } | \, X _ { \mu, \tau } M _ { \mu } \, | \psi _ { \sigma } \rangle | ^ { 2 }.$$

Using Eqs.(22) and (24) this simplifies to

$$F = \sum _ { \sigma, \mu, \tau } \frac { e ^ { - \beta H [ \sigma ] } } { \mathcal { Z } _ { \beta } } P _ { \mu, \tau } \prod _ { < i, j > \in A, \bar { A } } \delta _ { \mu _ { i, j }, \sigma _ { i } \sigma _ { j } } \delta _ { \sigma _ { 1 } \sigma _ { \bar { i } }, \tau }.$$

Plugging the expression for P µ,τ into Eq. (27), we get

$$F = & \sum _ { \sigma, \tau } \frac { e ^ { - \beta H [ \sigma ] } } { \mathcal { Z } _ { \beta } } \frac { e ^ { \beta J \tau \sum _ { i \geq 1 } \sigma _ { 1 } \sigma _ { 1 } \sigma _ { 1 } \sigma _ { 1 } \sigma _ { i } } } { \sum _ { \tau ^ { \prime } = \pm 1 } e ^ { \beta J \tau ^ { \prime } \sum _ { i \geq 1 } \sigma _ { 1 } \sigma _ { 1 } \sigma _ { 1 } \sigma _ { 1 } \sigma _ { i } } } \delta _ { \sigma _ { 1 } \sigma _ { 1 }, \tau } \\ & _ { - \, \, \,$$

$$= & \sum _ { \sigma } ^ { \cdot \cdot } \frac { e ^ { - \beta H [ \sigma ] } } { \mathcal { Z } _ { \beta } } \frac { e ^ { \beta J \sum _ { i \geq 1 } \sigma _ { i } \sigma _ { i } } } { \sum _ { \tau ^ { \prime } = \pm 1 } e ^ { \beta \tau ^ { \prime } \sigma _ { 1 } \sigma _ { 1 } J \sum _ { i \geq 1 } \sigma _ { i } \sigma _ { i } } } } \\ & \dots \cdot \cdot \cdot / \quad / \quad \ \ \ \ \ \$$

$$= & \frac { 1 } { 2 } \sum _ { \sigma } \frac { e ^ { - \beta H [ \sigma ] } } { \mathcal { Z } _ { \beta } } \left ( 1 - \tanh \left ( \beta J \sum _ { i \geq 1 } \sigma _ { i } \sigma _ { \bar { i } } \right ) \right )$$

$$= & \frac { 1 } { 2 } \left [ 1 - \left \langle \tanh ( \beta H _ { \partial } ) \right \rangle \right ]$$

where H ∂ is the boundary Hamiltonian between A and ¯ . A

## Mapping to the d -1 repetition code

The problem of choosing which domain wall configuration to assume can be understood as a decoding problem for a ( d -1)-dimensional classical repetition code. To see this, note that the decoding task in a classical error correcting

code is essentially figuring out what error has happened based on a given set of syndromes. In particular the decoding task in the classical repetition code on N bits can be described as follows. Consider an error vector e ∈ { 1 , - } 1 N such that e i = -1 if a bit flip has happened on the i 'th bit and e i = 1 otherwise. The goal of the decoder is to figure out e , knowing only the syndrome values, i.e. e e i j for adjacent i and j on the appropriate graph, e.g. the d -dimensional square lattice in the standard d -dimensional repetition code. Moreover, the decoder assumes an underlying noise model, i.e. a random distribution from which the vector e is drawn. For example, the independent and identically distributed (iid) noise model assumes each component e i is -1 with probability p and 1 with probability 1 -p , independent of other components.

In our setting finding the domain wall configuration is clearly the same task, with e i defined to be e i ≡ σ σ i ¯ i . Note that we cannot know the value of e i , but we know the value of e e i j = µ i,j µ ¯ ¯ i,j for adjacent i and j . At low temperatures, the domain walls are rare so most e i s would be 1, and as one increases the temperature more e i s would be -1. The underlying noise model is given by P µ,τ as described in the previous section. This noise model is not iid, but rather a correlated noise model. However, as we argue in the next section, the repetition code is able to guess e correctly with probability 1 in the thermodynamic limit unless β = 0, at which point each e i = ± 1 completely at random.

## Threshold for the Repetition Code

Here we argue that the repetition code decoding under the noise distribution P µ,τ is successful with probability 1 for any β &gt; 0 in the thermodynamic limit. First we review the argument for maximum error correction threshold in the presence of the iid noise and then we argue that the same should remain true for more generic noise models. We use quantum notation for convenience although the problem we are considering is inherently classical.

Consider an initial state | Ψ undergoing some errors, leaving it in the mixed state ⟩

$$\rho = \sum _ { e } p _ { e } \ e \left | \Psi \right \rangle \left \langle \Psi \right | e$$

where e is a product of X bit flips and p e is the probability of having that error. We rewrite the sum to distinguish between the errors that leave ρ with fewer than N/ 2 bit flips and errors with higher weights:

$$\rho = \sum _ { e \colon | e | < N / 2 } p _ { e } \ e \left | \Psi \right \rangle \left \langle \Psi \right | e + \sum _ { e \colon | e | \geq N / 2 } p _ { e } \ e \left | \Psi \right \rangle \left \langle \Psi \right | e.$$

The error correction channel, R , measures the syndromes and between two possible error configuration consistent with a syndrome, applies the one with the smaller weight. Therefore it always successfully corrects the errors when there are less than N/ 2 bit flips. Hence,

$$\mathcal { R } ( \rho ) & = \sum _ { e \colon | e | < N / 2 } p _ { e } \ \mathcal { R } ( e \, | \Psi \rangle \, \langle \Psi | \, e ) + \sum _ { e \colon | e | \geq N / 2 } p _ { e } \ \mathcal { R } ( e \, | \Psi \rangle \, \langle \Psi | \, e ) \\ & = \sum _ { e \colon | e | < N / 2 } p _ { e } \, | \Psi \rangle \, \langle \Psi | + \sum _ { e \colon | e | \geq N / 2 } p _ { e } \ \mathcal { R } ( e \, | \Psi \rangle \, \langle \Psi | \, e ).$$

The fidelity of this post-error correction state with the original state is lower bounded as,

$$F & = \langle \Psi | \, \mathcal { R } ( \rho ) \, | \Psi \rangle \\ & = \sum _ { e \colon | e | < N / 2 } p _ { e } + \sum _ { e \colon | e | \geq N / 2 } p _ { e } \, \langle \Psi | \, \mathcal { R } ( e | \Psi \rangle \, \langle \Psi | \, e ) \, | \Psi \rangle \geq \sum _ { e \colon | e | < N / 2 } p _ { e }$$

where the lower bound is because the second term in the expression above is non-negative.

Now we will show that this lower bound is equal to one in the thermodynamic limit. First, we rearrange the sum to be over | e | = k , the number of bit flips error, and because p k is a binomial distribution, p k = ( N k ) p k (1 -p ) N -p . In the limit of large N , we can approximate the binomial distribution as a Gaussian distribution with mean Np and standard deviation √ Np (1 -p ) via the central limit theorem. Since N/ 2 is N/ 2 -Np √ Np (1 -p ) = O ( √ N ) standard deviations away from the peak of the Gaussian, we see that lim N →∞ ∑ N/ 2 -1 k =0 p k = 1 when p &lt; 1 / 2. Since fidelity is upper-bounded by one, the fidelity is indeed equal to one in the thermodynamic limit. Thus, the repetition code will work successfully when p &lt; 1 / 2.

For the more general case where e is drawn based on the P µ,τ distribution, the bound in Eq.(37) is still valid. Therefore, to show perfect recovery we only need to show that with probability one, less than half of the edges contain a domain wall. Intuitively, this is true because for β &gt; 0, not having a domain wall on an edge is energetically favorable. To make it rigorous, one can define indicator variables I i = (1 -e i ) / 2 in terms of which the number of domain walls can be written as k = ∑ i I i . When β &gt; 0 it is more probable for I i to be 0 rather than 1 and hence the expectation value of k/ ∂A | | is less than 1 / 2. Noting that lim | i - |→∞ j Cov( I , I i j ) = 0, the weak law of large numbers for weakly dependent variables [57] shows that k/ ∂A | | concentrates around its mean, which is less than 1 / 2 when β &gt; 0 in the thermodynamic limit.

## Negativity Spectrum

Here we calculate the entanglement negativity between a sub-region A and its complement ¯ , A in the state ρ β . Let σ denote the spin configuration that is related to σ by flipping every spin. Accordingly, the state | ψ σ ⟩ can be written as | ψ σ ⟩ = | σ ⟩ + | σ ⟩ √ √ 2 and Eq. (2) can be rewritten as,

$$\rho _ { \beta } = \frac { 1 } { 2 \mathcal { Z } _ { \beta } } \sum _ { \sigma } e ^ { - \beta H [ \sigma ] } ( | \sigma \rangle \left \langle \sigma | + | \bar { \sigma } \right \rangle \left \langle \bar { \sigma } | + | \sigma \right \rangle \left \langle \bar { \sigma } | + | \bar { \sigma } \right \rangle \left \langle \sigma | ),$$

where H σ [ ] = -J ∑ ⟨ i,j ⟩ σ σ i j is the classical Ising energy of the spin configuration σ . Taking the partial transpose of ρ β with respect to sub-region A yields,

$$\rho _ { \beta } ^ { T _ { A } } & = \frac { 1 } { 2 \mathcal { Z } _ { \beta } } \sum _ { \sigma } \left [ e ^ { - \beta H [ \sigma ] } \left ( | \sigma \right ) \langle \sigma | + | \bar { \sigma } \rangle \langle \bar { \sigma } | \right ) + e ^ { - \beta H [ \sigma _ { A }, \sigma _ { A } ] } \left ( | \bar { \sigma } _ { A }, \sigma _ { A } \rangle \langle \sigma _ { A }, \bar { \sigma } _ { A } | + | \sigma _ { A }, \bar { \sigma } _ { A } \rangle \langle \bar { \sigma } _ { A }, \sigma _ { A } | \right ) \right ] \\ & = \frac { 1 } { 2 \mathcal { Z } _ { \beta } } \sum _ { \sigma } \left [ e ^ { - \beta H [ \sigma ] } \left ( | \sigma \right ) \langle \sigma | + | \bar { \sigma } \rangle \langle \bar { \sigma } | \right ) + e ^ { - \beta H [ \sigma _ { A }, \bar { \sigma } _ { A } ] } \left ( | \bar { \sigma } \rangle \langle \sigma | + | \sigma \rangle \langle \bar { \sigma } | \right ) \right ] \\, \quad, \quad, \quad, \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad.$$

where σ A and σ ¯ A denote the spin configuration σ restericted to regions A and ¯ respectively. A Thus, the eigenvalues and eigenstates of ρ T A A are

$$\lambda _ { \sigma ; \pm } = \frac { e ^ { - \beta H [ \sigma ] } \pm e ^ { - \beta H [ \sigma _ { A }, \bar { \sigma } _ { \lambda } ] } } { Z _ { \beta } }, \quad \left | \psi _ { \sigma } ^ { ( \pm ) } \right \rangle \equiv \frac { | \sigma \rangle \pm | \bar { \sigma } \rangle } { \sqrt { 2 } }. \quad \left | \psi _ { \sigma } ^ { ( \pm ) } \right \rangle \equiv \frac { | \sigma \rangle \pm | \bar { \sigma } \rangle } { \sqrt { 2 } }. \quad \left | \psi _ { \sigma } ^ { ( 4 0 ) } \right \rangle$$

We can now calculate the entanglement negativity from the sum of all the negative eigenvalues. Firstly, it is only the λ σ ; -eigenvalues which can be negative. We can also see that λ σ ; -= λ ¯; σ -from the Z 2 symmetry of the Hamiltonian and they correspond to the same eigenvector, up to some phase. Moreover, note that λ σ ; -= -λ σ A ,σ ¯ ¯ A ; -so for a given spin configuration, either λ σ ; -or λ σ A ,σ ¯ ¯ A ; -would be negative. Thus, we account for the two incidences of double counting with an additional 1/4 and get

$$\mathcal { N } _ { A } [ \rho _ { \beta } ] = \frac { 1 } { 4 \mathcal { Z } _ { \beta } } \sum _ { \sigma } \left | e ^ { - \beta H [ \sigma ] } - e ^ { - \beta H [ \sigma _ { A }, \bar { \sigma } _ { A } ] } \right |.$$

Factoring out e -βH σ [ ] we get

$$\mathcal { N } _ { A } [ \rho _ { \beta } ] & = \frac { 1 } { 4 \mathcal { Z } _ { \beta } } \sum _ { \sigma } e ^ { - \beta H [ \sigma ] } \left | 1 - e ^ { - \beta ( H [ \sigma _ { A }, \vartheta _ { A } ] - H [ \sigma ] ) } \right | \\ & = \frac { 1 } { 4 } \left \langle \left | 1 - e ^ { 2 \beta H _ { \beta } ( \sigma ) } \right | \right \rangle.$$

One may now use classical Monte Carlo simulations to compute the thermal average in Eq.(42). However, care must be taken with sampling the spin configurations, since the expression inside the absolute value could be exponentially large for configurations with an exponentially small Boltzmann weight, in such a way that the overall contribution to the expectation value is O (1).

The aforementioned issue can be avoided by using a technique known as importance sampling. Let q σ [ ] denote the probability distribution defined as q σ [ A , σ ¯ A ] = 1 2 p σ [ A , σ ¯ A ] + 1 2 p σ [ A , σ ¯ A ], where p σ [ ] = e βJ ∑ ⟨ i,j ⟩ σ σ i j / Z β . To sample

from q σ [ ], one can simply sample from p σ [ ] and then with 1 / 2 probability, flip the spins in ¯ . A We may thus write Eq.(42) as,

$$\mathcal { N } _ { A } [ \rho _ { \beta } ] & = \frac { 1 } { 4 } \sum _ { \sigma } q [ \sigma ] \frac { p [ \sigma ] } { q [ \sigma ] } | 1 - e ^ { 2 \beta H _ { \theta } [ \sigma ] } | \\ & = \frac { 1 } { 2 } \left \langle | \tanh ( \beta H _ { \theta } [ \sigma ] ) | \right \rangle _ { q }. \\. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. - \.. \.. \..$$

where the final expectation value is taken with respect to the distribution q . The advantage of Eq.(43) compared to Eq.(42) is that the expression inside the expectation value is now bounded by unity everywhere, so that contributions from rare configurations are no longer a concern.
Yonna Kim Ali Lavasani Sagar Vijay
Persistent Topological Negativity in a High-Temperature Mixed-State
We study the entanglement structure of the Greenberger-Horne-Zeilinger (GHZ)
state as it thermalizes under a strongly-symmetric quantum channel describing
the Metropolis-Hastings dynamics for the $d$-dimensional classical Ising model
at inverse temperature $\beta$. This channel outputs the classical Gibbs state
when acting on a product state in the computational basis. When applying this
channel to a GHZ state in spatial dimension $d>1$, the resulting mixed state
changes character at the Ising phase transition temperature from being
long-range entangled to short-range-entangled as temperature increases.
Nevertheless, we show that the topological entanglement negativity of a large
region is insensitive to this transition and takes the same value as that of
the pure GHZ state at any finite temperature $\beta>0$. We establish this
result by devising a local operations and classical communication (LOCC)
``decoder" that provides matching lower and upper bounds on the negativity in
the thermodynamic limit which may be of independent interest. This perspective
connects the negativity to an error-correction problem on the
$(d-1)$-dimensional bipartitioning surface and explains the persistent
negativity in certain correlated noise models found in previous studies.
Numerical results confirm our analysis.