arxiv: Statistical mechanical models for quantum codes subject to correlated noise

Transcription

1 Statistical mechanical models for quantum codes subject to correlated noise Christopher Chubb (EQUS USyd), Steve Flammia (EQUS USyd/YQI Yale) arxiv: Perimeter Institute

2 Quantum codes Two important questions about quantum codes: How do I decode? What is the threshold? The statistical mechanical mapping 1 allows us to address both questions codes subject to Pauli noise. 1 Dennis, Kitaev, Landahl, Preskill, JMP 2002, doi:10/cs2mtf, arxiv:quant-ph/

3 Quantum codes Two important questions about quantum codes: How do I decode optimally? What is the threshold? The statistical mechanical mapping 1 allows us to address both questions codes subject to Pauli noise. 1 Dennis, Kitaev, Landahl, Preskill, JMP 2002, doi:10/cs2mtf, arxiv:quant-ph/

4 Quantum codes Two important questions about quantum codes: How do I decode optimally? What is the fault-tolerant threshold? The statistical mechanical mapping 1 allows us to address both questions codes subject to Pauli noise. 1 Dennis, Kitaev, Landahl, Preskill, JMP 2002, doi:10/cs2mtf, arxiv:quant-ph/

5 Quantum codes Two important questions about quantum codes: How do I decode optimally? What is the fault-tolerant threshold? The statistical mechanical mapping 1 allows us to address both questions codes subject to Pauli noise. 1 Dennis, Kitaev, Landahl, Preskill, JMP 2002, doi:10/cs2mtf, arxiv:quant-ph/

6 Statistical mechanical mapping The idea here it to construct a family of statistical mechanical models, whose thermodynamic properties reflect the error correction properties of the code. X X X X Z Z Z Z Toric code Eight-vertex model This will allow us to use the analytic and numerical tools developed to study stat mech systems to study quantum codes.

7 Statistical mechanical mapping Stabiliser code & Pauli noise Disordered statistical mechanical model Threshold Phase transition Decoding Calculating partition functions Allows us to reappropriate techniques for studying stat. mech. systems to study quantum codes: Threshold approximation Optimal decoding Monte Carlo simulation Partition function calculation

8 Statistical mechanical mapping Stabiliser code & Pauli noise Disordered statistical mechanical model Threshold Phase transition Decoding Calculating partition functions Allows us to reappropriate techniques for studying stat. mech. systems to study quantum codes: Threshold approximation Optimal decoding Monte Carlo simulation Partition function calculation

9 Statistical mechanical mapping Stabiliser code & Pauli noise Disordered statistical mechanical model Threshold Phase transition Decoding Calculating partition functions Allows us to reappropriate techniques for studying stat. mech. systems to study quantum codes: Threshold approximation Optimal decoding Monte Carlo simulation Partition function calculation

10 Our results Generalise mapping to correlated noise for arbitrary subsystem codes Numerically show that mild correlations can lower the threshold of the toric code considerably Show how to apply our mapping to circuit noise via the history code, allowing us to approximate fault tolerant thresholds Show that the stat. mech. mapping gives tensor network maximum likelihood decoders which generalises the MPS decoder of Bravyi, Suchara and Vargo

11 Our results Generalise mapping to correlated noise for arbitrary subsystem codes Numerically show that mild correlations can lower the threshold of the toric code considerably Show how to apply our mapping to circuit noise via the history code, allowing us to approximate fault tolerant thresholds Show that the stat. mech. mapping gives tensor network maximum likelihood decoders which generalises the MPS decoder of Bravyi, Suchara and Vargo

12 Our results Generalise mapping to correlated noise for arbitrary subsystem codes Numerically show that mild correlations can lower the threshold of the toric code considerably Show how to apply our mapping to circuit noise via the history code, allowing us to approximate fault tolerant thresholds Show that the stat. mech. mapping gives tensor network maximum likelihood decoders which generalises the MPS decoder of Bravyi, Suchara and Vargo

13 Our results Generalise mapping to correlated noise for arbitrary subsystem codes Numerically show that mild correlations can lower the threshold of the toric code considerably Show how to apply our mapping to circuit noise via the history code, allowing us to approximate fault tolerant thresholds Show that the stat. mech. mapping gives tensor network maximum likelihood decoders which generalises the MPS decoder of Bravyi, Suchara and Vargo

14 Stabiliser codes and Pauli noise For qubits, the Paulis P := {I, X, Y, Z} are defined ( ) ( ) ( ) ( ) i 1 0 I :=, X :=, Y :=, Z := i We will be considering stabiliser codes, which are specified by an Abelian subgroup of the Paulis S, and whose code space C is the joint +1 eigenspace, { } C = ψ S ψ = ψ, S S. Any two errors which differ by a stabiliser are logically equivalent, so the logical classes of errors are E := {ES S S}

15 Independent case: Hamiltonian Let A, B be the scalar commutator of two Paulis, such that AB =: A, B BA. For a stabiliser code generated by {S k } k, and an error Pauli E, the (disordered) Hamiltonian H E is defined H E ( s) := i σ P i Coupling {}}{ J i (σ) DoF Disorder {}}{{ }}{ σ, E k: σ,s k = 1 s k for s k = ±1, and coupling strengths J i (σ) R. Take-aways: Ising-type, with interactions corresponding to single-site Paulis σ Disorder E flips some interactions (Ferro Anti-ferro) Local code = local stat. mech. model

16 Independent case: Hamiltonian Let A, B be the scalar commutator of two Paulis, such that AB =: A, B BA. For a stabiliser code generated by {S k } k, and an error Pauli E, the (disordered) Hamiltonian H E is defined H E ( s) := i σ P i Coupling {}}{ J i (σ) DoF Disorder {}}{{ }}{ σ, E k: σ,s k = 1 s k for s k = ±1, and coupling strengths J i (σ) R. Take-aways: Ising-type, with interactions corresponding to single-site Paulis σ Disorder E flips some interactions (Ferro Anti-ferro) Local code = local stat. mech. model

17 Independent case: Hamiltonian Let A, B be the scalar commutator of two Paulis, such that AB =: A, B BA. For a stabiliser code generated by {S k } k, and an error Pauli E, the (disordered) Hamiltonian H E is defined H E ( s) := i σ P i Coupling {}}{ J i (σ) DoF Disorder {}}{{ }}{ σ, E k: σ,s k = 1 s k for s k = ±1, and coupling strengths J i (σ) R. Take-aways: Ising-type, with interactions corresponding to single-site Paulis σ Disorder E flips some interactions (Ferro Anti-ferro) Local code = local stat. mech. model

18 Independent case: Gauge symmetry H E ( s) = i J i (σ) σ, E σ P i k: σ,s k = 1 s k Using A, B A, C = A, BC, we see this system has a gauge symmetry s k s k and E ES k. This gauge symmetry will capture the logical equivalence of errors, Z E = Z ESk.

19 Independent case: Gauge symmetry H E ( s) = i J i (σ) σ, E σ P i k: σ,s k = 1 s k Using A, B A, C = A, BC, we see this system has a gauge symmetry s k s k and E ES k. This gauge symmetry will capture the logical equivalence of errors, Z E = Z ESk.

20 Independent case: Nishimori conditon Suppose we have an independent error model Pr(E) = i p i (E i ), we now want Z E = Pr(E). Using the gauge symmetry we have that the partition function can be written as a sum stabiliser-equivalent errors Z E = s e βh E ( s) = S e βh ES ( 1) = F E e βh F ( 1). If we select the coupling strength such that e βh E ( 1) = Pr(E), then Z E = Pr(E) will follow.

21 Independent case: Nishimori conditon Suppose we have an independent error model Pr(E) = i p i (E i ), we now want Z E = Pr(E). Using the gauge symmetry we have that the partition function can be written as a sum stabiliser-equivalent errors Z E = s e βh E ( s) = S e βh ES ( 1) = F E e βh F ( 1). If we select the coupling strength such that e βh E ( 1) = Pr(E), then Z E = Pr(E) will follow.

22 Independent case: Nishimori conditon Suppose we have an independent error model Pr(E) = i p i (E i ), we now want Z E = Pr(E). Using the gauge symmetry we have that the partition function can be written as a sum stabiliser-equivalent errors Z E = s e βh E ( s) = S e βh ES ( 1) = F E e βh F ( 1). If we select the coupling strength such that e βh E ( 1) = Pr(E), then Z E = Pr(E) will follow.

23 Independent case: Nishimori condition We now want to pick our couplings such that e βh E ( 1) = Pr(E). Expanding this out, we get log p i (E) = βj i (σ) σ, E. i i σ Using the Fourier-like orthogonality relation 1 4 σ σ, τ = δ τ,i, this becomes Nishimori condition: βj i (σ) = 1 log p i (τ) σ, τ, 4 τ P which implies e βh E ( 1) = Pr(E), and therefore Z E = Pr(E). This intrinsically links the error correcting behaviour of the code to the thermodynamic behaviour of the model (along the Nishimori line).

24 Independent case: Nishimori condition We now want to pick our couplings such that e βh E ( 1) = Pr(E). Expanding this out, we get log p i (E) = βj i (σ) σ, E. i i σ Using the Fourier-like orthogonality relation 1 4 σ σ, τ = δ τ,i, this becomes Nishimori condition: βj i (σ) = 1 log p i (τ) σ, τ, 4 τ P which implies e βh E ( 1) = Pr(E), and therefore Z E = Pr(E). This intrinsically links the error correcting behaviour of the code to the thermodynamic behaviour of the model (along the Nishimori line).

25 Independent case: Nishimori condition We now want to pick our couplings such that e βh E ( 1) = Pr(E). Expanding this out, we get log p i (E) = βj i (σ) σ, E. i i σ Using the Fourier-like orthogonality relation 1 4 σ σ, τ = δ τ,i, this becomes Nishimori condition: βj i (σ) = 1 log p i (τ) σ, τ, 4 τ P which implies e βh E ( 1) = Pr(E), and therefore Z E = Pr(E). This intrinsically links the error correcting behaviour of the code to the thermodynamic behaviour of the model (along the Nishimori line).

26 Toric code and the random-bond Ising model Step 0: Code and noise model Toric code with iid bit-flips Step 1: Degrees of freedom s v = ±1 on each vertex v Step 2: Interactions H I = v v J s v s v Step 3: Disorder H E = v v Je vv s v s v where e vv = { +1 Evv = I, 1 E vv = X.

27 Toric code and the random-bond Ising model Step 0: Code and noise model Toric code with iid bit-flips Step 1: Degrees of freedom s v = ±1 on each vertex v Step 2: Interactions H I = v v J s v s v Step 3: Disorder H E = v v Je vv s v s v where e vv = { +1 Evv = I, 1 E vv = X.

28 Toric code and the random-bond Ising model Step 0: Code and noise model Toric code with iid bit-flips Step 1: Degrees of freedom s v = ±1 on each vertex v Step 2: Interactions H I = v v J s v s v Step 3: Disorder X X X X H E = v v Je vv s v s v where e vv = { +1 Evv = I, 1 E vv = X.

29 Toric code and the random-bond Ising model Step 0: Code and noise model Toric code with iid bit-flips Pr(X e ) = p, Pr(I e ) = 1 p. Step 1: Degrees of freedom s v = ±1 on each vertex v Step 2: Interactions H I = J s v s v v v Step 3: Disorder X X X X H E = v v Je vv s v s v where e vv = { +1 Evv = I, 1 E vv = X.

30 Toric code and the random-bond Ising model Step 0: Code and noise model Toric code with iid bit-flips Pr(X e ) = p, Pr(I e ) = 1 p. Step 1: Degrees of freedom s v = ±1 on each vertex v Step 2: Interactions H I = v v J s v s v Step 3: Disorder H E = v v Je vv s v s v where e vv = { +1 Evv = I, 1 E vv = X.

31 Toric code and the random-bond Ising model Step 0: Code and noise model Toric code with iid bit-flips Pr(X e ) = p, Pr(I e ) = 1 p. Step 1: Degrees of freedom s v = ±1 on each vertex v Step 2: Interactions H I = v v J s v s v Step 3: Disorder H E = v v Je vv s v s v where e vv = { +1 Evv = I, 1 E vv = X.

32 Toric code and the random-bond Ising model Step 0: Code and noise model Toric code with iid bit-flips Pr(X e ) = p, Pr(I e ) = 1 p. Step 1: Degrees of freedom s v = ±1 on each vertex v Step 2: Interactions H I = J s v s v v v I Y X Step 3: Disorder Z H E = v v Je vv s v s v where e vv = { +1 Evv = I, 1 E vv = X. σ e, S v = { 1 if σ = X, Y, v e, +1 else.

33 Toric code and the random-bond Ising model Step 0: Code and noise model Toric code with iid bit-flips Pr(X e ) = p, Pr(I e ) = 1 p. Step 1: Degrees of freedom s v = ±1 on each vertex v Step 2: Interactions H I = J s v s v v v I Y X Step 3: Disorder Z H E = v v Je vv s v s v where e vv = { +1 Evv = I, 1 E vv = X. σ e, S v = { 1 if σ = X, Y, v e, +1 else.

34 Toric code and the random-bond Ising model Step 0: Code and noise model Toric code with iid bit-flips Pr(X e ) = p, Pr(I e ) = 1 p. Step 1: Degrees of freedom s v = ±1 on each vertex v Step 2: Interactions H I = v v J s v s v Step 3: Disorder H E = v v Je vv s v s v where e vv = { +1 Evv = I, 1 E vv = X.

35 Toric code and the random-bond Ising model Step 0: Code and noise model Toric code with iid bit-flips Pr(X e ) = p, Pr(I e ) = 1 p. Step 1: Degrees of freedom s v = ±1 on each vertex v Step 2: Interactions H I = v v J s v s v Step 3: Disorder H E = v v Je vv s v s v where e vv = { +1 Evv = I, 1 E vv = X. σ, S k σ, ES k

36 Toric code and the random-bond Ising model Step 0: Code and noise model Toric code with iid bit-flips Pr(X e ) = p, Pr(I e ) = 1 p. Step 1: Degrees of freedom s v = ±1 on each vertex v Step 2: Interactions H I = v v J s v s v Step 3: Disorder H E = v v Je vv s v s v where e vv = { +1 Evv = I, 1 E vv = X. σ, S k σ, ES k

37 Toric code and the random-bond Ising model Step 0: Code and noise model Toric code with iid bit-flips Pr(+J) = p, Pr( J) = 1 p. Step 1: Degrees of freedom s v = ±1 on each vertex v Step 2: Interactions H I = v v J s v s v Step 3: Disorder H E = v v Je vv s v s v where e vv = { +1 Evv = I, 1 E vv = X. ±J Random-bond Ising Model

38 Other independent examples Toric code Bit-flip Random-bond Ising 1 Indep. X &Z 2 Random-bond Ising Depolarising Random 8-vertex model 2 Colour code Bit-flip Random 3-spin Ising Indep. X &Z 2 Random 3-spin Ising Depolarising Random interacting 8-vertex 2 X X Z Y 1 Dennis et.al., JMP 2002, doi:10/cs2mtf, arxiv:quant-ph/ Bombin et.al., PRX 2012, doi:10/crz5, arxiv:

39 Other independent examples Toric code Bit-flip Random-bond Ising 1 Indep. X &Z 2 Random-bond Ising Depolarising Random 8-vertex model 2 Colour code Bit-flip Random 3-spin Ising Indep. X &Z 2 Random 3-spin Ising Depolarising Random interacting 8-vertex 2 X X Z Z Y 1 Dennis et.al., JMP 2002, doi:10/cs2mtf, arxiv:quant-ph/ Bombin et.al., PRX 2012, doi:10/crz5, arxiv:

40 Other independent examples Toric code Bit-flip Random-bond Ising 1 Indep. X &Z 2 Random-bond Ising Depolarising Random 8-vertex model 2 Colour code Bit-flip Random 3-spin Ising Indep. X &Z 2 Random 3-spin Ising Depolarising Random interacting 8-vertex 2 X Y X Z Z Y 1 Dennis et.al., JMP 2002, doi:10/cs2mtf, arxiv:quant-ph/ Bombin et.al., PRX 2012, doi:10/crz5, arxiv:

41 Other independent examples Toric code Bit-flip Random-bond Ising 1 Indep. X &Z 2 Random-bond Ising Depolarising Random 8-vertex model 2 Colour code Bit-flip Random 3-spin Ising Indep. X &Z 2 Random 3-spin Ising Depolarising Random interacting 8-vertex 2 X Y X Z Z Y 1 Dennis et.al., JMP 2002, doi:10/cs2mtf, arxiv:quant-ph/ Bombin et.al., PRX 2012, doi:10/crz5, arxiv:

42 Error correction threshold as a quenched phase transition Consider the free energy cost of a logical error L, E (L) = 1 β log Z EL + 1 β log Z E. Along the Nishimori line which implies E (L) = 1 β Below threshold : log Pr(E) Pr(EL), E (L) (in mean) Above threshold : E (L) 0 (in prob.)

43 Error correction threshold as a quenched phase transition Consider the free energy cost of a logical error L, E (L) = 1 β log Z EL + 1 β log Z E. Along the Nishimori line which implies E (L) = 1 β Below threshold : log Pr(E) Pr(EL), E (L) (in mean) Above threshold : E (L) 0 (in prob.)

44 Phase diagram sketch Temperature Ordered EC possible E (L) Nishimori condition Disordered EC impossible E (L) 0 ML threshold ME threshold Noise parameter

45 Correlated case The key point independence gave us was the ability to factor our noise model Pr(E) = i p i (E i ). We can generalise this to correlated models: Factored distribution An error model factors over regions {R j } j if there exist φ j : P Rj Pr(E) = ( ) φ j ERj j R such that This model includes many probabilistic graphical models, such as Bayesian Networks and Markov/Gibbs Random Fields.

46 Correlated case The key point independence gave us was the ability to factor our noise model Pr(E) = i p i (E i ). We can generalise this to correlated models: Factored distribution An error model factors over regions {R j } j if there exist φ j : P Rj Pr(E) = ( ) φ j ERj j R such that This model includes many probabilistic graphical models, such as Bayesian Networks and Markov/Gibbs Random Fields.

47 Correlated case By construction, we can extend to the correlated case by changing σ P i to σ P Rj : H E ( s) := j J j (σ) σ, E σ P Rj k: σ,s k = 1 s k Nishimori condition: βj j (σ) = 1 P Rj τ P Rj log φ j (τ) σ, τ, As before we get that Z E = Pr(E), and so the threshold manifests as a phase transition.

48 Correlated case By construction, we can extend to the correlated case by changing σ P i to σ P Rj : H E ( s) := j J j (σ) σ, E σ P Rj k: σ,s k = 1 s k Nishimori condition: βj j (σ) = 1 P Rj τ P Rj log φ j (τ) σ, τ, As before we get that Z E = Pr(E), and so the threshold manifests as a phase transition.

49 Correlated example Toric code with correlated bit-flips Correlations induce longer-range interactions

50 Correlated example Toric code with correlated bit-flips Correlations induce longer-range interactions

51 Correlated example Toric code with correlated bit-flips Correlations induce longer-range interactions

52 Correlated example Toric code with correlated bit-flips Correlations induce longer-range interactions

53 Across plaquette correlated bit-flips This error model is entirely specified by the conditional error probabilities Pr(I e I e ) Pr(I e X e ) Pr(X e I e ) Pr(X e X e ) e e for all neighbouring edges e and e. For our purposes, it will convenient to parameterise things by p := Pr(X e ), η := Pr(X e X e ) Pr(X e I e ). Here p is the marginal error rate, and η is a measure of the correlations.

54 Across plaquette correlated bit-flips This error model is entirely specified by the conditional error probabilities Pr(I e I e ) Pr(I e X e ) Pr(X e I e ) Pr(X e X e ) e e for all neighbouring edges e and e. For our purposes, it will convenient to parameterise things by p := Pr(X e ), η := Pr(X e X e ) Pr(X e I e ). Here p is the marginal error rate, and η is a measure of the correlations.

55 Monte Carlo simulations 2.5 η = 1 (indep.) η = 2 (corr.) 2 Temperature T Marginal error probability p(%) Thresholds Indep.: p th = (2)% 1,2 Corr.: p th = 10.04(6)% 1 Dennis et.al., JMP 2002, doi:10/cs2mtf, arxiv:quant-ph/ Toldin et.al., JSP 2009, doi:10/c3r2kc, arxiv:

56 Circuit noise Noise followed by ideal measurements is unrealistic. In reality, circuits will be faulty. Noisy circuits give faulty outcomes, and will spread around and correlate noise: X X X X X

57 Circuit noise We will consider measurement circuits of the form M U where U is a Clifford and M is a Pauli. For convenience we only consider independent noise on each circuit. We also will push noise through until after the unitary: I I I X Y X Y Z I X = X Y

58 Circuit noise We will consider measurement circuits of the form M U where U is a Clifford and M is a Pauli. For convenience we only consider independent noise on each circuit. We also will push noise through until after the unitary: I I I X Y X Y Z I X = X Y

59 Circuit noise Ideal measurements Time

60 Circuit noise Ideal measurements Time

61 Circuit noise Ideal measurements Time

62 Circuit noise Ideal measurements Time

63 Circuit noise Ideal measurements Pr(Error history) = j,t Time Pr ( Error after meas. (x, t) Error before meas. (x, t) )

64 Circuit noise Ideal measurements Pr(Error history) = j,t Time Pr ( Error after meas. (x, t) Error before meas. (x, t) )

65 Circuit noise Ideal measurements Pr(Error history) = j,t Not time Pr ( Error after meas. (x, t) Error before meas. (x, t) )

66 History code Ideal measurements History code: Qubits placed at points in space-time (including ancillae) Stabiliers correspond to measurements and stabilisers at final time Logicals correspond to logicals at final time Circuit noise on original code = Spatially correlated noise on history code FT decoding of original code = Decoding the history code FT threshold of original code = Threshold of history code

67 History code Ideal measurements History code: Qubits placed at points in space-time (including ancillae) Stabiliers correspond to measurements and stabilisers at final time Logicals correspond to logicals at final time Circuit noise on original code = Spatially correlated noise on history code FT decoding of original code = Decoding the history code FT threshold of original code = Threshold of history code

68 Decoding Can the stat. mech. model give us a decoder? If an error E occurs, a decoder needs to select one of the degenerate logical error classes E EL 1 EL 2 EL 3... The optimal (maximum likelihood) decoder selects the most likely class D ML = EL l where l = arg max Pr ( ) EL l. l

69 Decoding Can the stat. mech. model give us a decoder? If an error E occurs, a decoder needs to select one of the degenerate logical error classes E EL 1 EL 2 EL 3... The optimal (maximum likelihood) decoder selects the most likely class D ML = EL l where l = arg max Pr ( ) EL l. l

70 Decoding Can the stat. mech. model give us a decoder? If an error E occurs, a decoder needs to select one of the degenerate logical error classes E EL 1 EL 2 EL 3... The optimal (maximum likelihood) decoder selects the most likely class D ML = EL l where l = arg max Pr ( ) EL l. l

71 Decoding from partition functions Along the Nishimori line, the maximum likelihood condition corresponds to maximising the partition function l = arg max Z ELl. l Approximating Z ELl therefore allows us to approximate the ML decoder. Step 1: Measure the syndrome s Step 2: Construct an arbitrary error C s which has syndrome s Step 3: Approximate Z Cs L l = Pr(C s L l ) for each logical l Step 4: Find the l such that Z Cs L l is maximised Step 5: Apply (C s L l ) 1

72 Decoding from partition functions Along the Nishimori line, the maximum likelihood condition corresponds to maximising the partition function l = arg max Z ELl. l Approximating Z ELl therefore allows us to approximate the ML decoder. Step 1: Measure the syndrome s Step 2: Construct an arbitrary error C s which has syndrome s Step 3: Approximate Z Cs L l = Pr(C s L l ) for each logical l Step 4: Find the l such that Z Cs L l is maximised Step 5: Apply (C s L l ) 1

73 Decoding from (approximate) tensor network contraction Partition functions can be expressed as tensor networks 2,3, allowing us to use approximate tensor network contraction schemes. For 2D codes and locally correlated noise, this tensor network is also 2D. Here we can use the MPS-MPO approximation contraction scheme considered by Bravyi, Suchara and Vargo 4 : 4 Verstraete et. al., PRL 2006, doi:10/dfgcz8, arxiv:quant-ph/ Bridgeman and Chubb, JPA 2017, doi:10/cv7m, arxiv: Bravyi, Scuhara, Vargo, PRA 2014, doi:10/cv7n, arxiv:

74 Decoding from (approximate) tensor network contraction Partition functions can be expressed as tensor networks 2,3, allowing us to use approximate tensor network contraction schemes. For 2D codes and locally correlated noise, this tensor network is also 2D. Here we can use the MPS-MPO approximation contraction scheme considered by Bravyi, Suchara and Vargo 4 : 4 Verstraete et. al., PRL 2006, doi:10/dfgcz8, arxiv:quant-ph/ Bridgeman and Chubb, JPA 2017, doi:10/cv7m, arxiv: Bravyi, Scuhara, Vargo, PRA 2014, doi:10/cv7n, arxiv:

75 Decoding from (approximate) tensor network contraction This gives an algorithm for (approximate) maximum likelihood decoding for any 2D code, subject to any locally correlated noise, generalising BSV. Indeed, applying this to iid noise in the surface code reproduces BSV: i b i r j j i r j r = i b j b = p(e e X i b +j b Z ir +jr ) i k = i k = δ i,j,k,l j b j r l l

76 Conclusions and further work Extended the stat. mech. mapping to correlated models Can apply stat. mech. mapping to circuit noise via the history code Numerically evaluated the threshold of correlated bit-flips in the toric code Stat. mech. mapping gives tensor network maximum likelihood decoders Can we apply this to experimentally relevant correlated models? Construction extends to all Abelian quantum doubles. Can we extend non-abelian models? Non-stochastic noise? Non-Pauli noise? Coherent noise? Can we use the decoders to understand to better understand the connection between correlation and the threshold (ongoing work with David Tuckett and Benjamin Brown). Thank you! ArXiv: