pith. sign in

arxiv: 2411.05785 · v2 · pith:HXXXEWMEnew · submitted 2024-11-08 · 🪐 quant-ph · cond-mat.dis-nn· cond-mat.stat-mech

Phases of decodability in the surface code with unitary errors

Pith reviewed 2026-05-23 17:13 UTC · model grok-4.3

classification 🪐 quant-ph cond-mat.dis-nncond-mat.stat-mech
keywords surface codeunitary errorsmaximum likelihood decodingstatistical mechanics modelentanglement scalingtransfer matrixdecodability phasestopological codes
0
0 comments X

The pith

Volume-law entanglement coexists with ferromagnetic order in the surface code under unitary errors, yielding a phase where encoded information survives but maximum-likelihood decoding fails.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper models maximum-likelihood decoding for the surface code with generic unitary errors via a statistical mechanics model whose complex weights are contracted using a (1+1)D transfer matrix. Information loss appears as a ferromagnetic-to-paramagnetic transition in this contraction, while entanglement scaling undergoes a separate area-to-volume-law transition. These transitions do not coincide, so volume-law entanglement can appear inside the ferromagnetic phase. Simulations for single- and two-qubit X-rotations locate the three phases, and tilting the rotation axis couples the X and Z sectors to produce the ferromagnetic volume-law regime. In that regime the logical information remains in principle but becomes effectively undecodable.

Core claim

The transfer-matrix contraction of the complex-weight statistical mechanics model reveals a ferromagnetic volume-law phase in which the encoded information is retained yet is effectively undecodable because entanglement scales with volume rather than area.

What carries the argument

The (1+1)D transfer-matrix contraction of the complex-weight statistical mechanics model for maximum-likelihood decoding success probability.

If this is right

  • A ferromagnetic area-law phase exists in which decoding succeeds with high probability.
  • A paramagnetic volume-law phase exists in which logical information is lost.
  • Tilting the rotation axis away from the X direction couples the X and Z statistical-mechanics models and can stabilize the ferromagnetic volume-law phase.
  • In the ferromagnetic volume-law phase, Z errors remain correctable in principle while recovery of the encoded classical information becomes hard.
  • An isometric tensor-network algorithm enables efficient syndrome sampling for these unitary-error simulations.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Decodability thresholds based solely on ferromagnetic order may overestimate practical performance when unitary errors produce volume-law entanglement.
  • The separation between order and entanglement transitions suggests that decoder design must track both quantities rather than error rate alone.
  • The same transfer-matrix framework could be applied to other topological codes whose error models map to complex-weight statistical mechanics.

Load-bearing premise

The (1+1)D transfer-matrix contraction of the complex-weight statistical mechanics model exactly captures the maximum-likelihood decoding success probability for generic unitary errors on the surface code.

What would settle it

Direct numerical comparison of the actual maximum-likelihood decoding success rate on small surface-code patches against the order parameter and entanglement scaling extracted from the same transfer-matrix contraction under identical unitary error rates.

Figures

Figures reproduced from arXiv: 2411.05785 by Sajant Anand, Yimu Bao.

Figure 1
Figure 1. Figure 1: (a) Two-dimensional surface code on the square lat [PITH_FULL_IMAGE:figures/full_fig_p001_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Decoding transition in the surface code under single- and two-qubit Pauli-X rotations. (a) Defect free energy as a [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Decoding in the surface code under single-qubit [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
Figure 1
Figure 1. Figure 1: Tensor network representation of the partition function of RBIM associated with single-qubit unitary X-rotation. [PITH_FULL_IMAGE:figures/full_fig_p008_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Tensor network representation of the partition function of RBIM associated with single-qubit unitary X-rotation. [PITH_FULL_IMAGE:figures/full_fig_p009_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Tensor network representation of the complex weight partition function associated with general single-qubit rotation [PITH_FULL_IMAGE:figures/full_fig_p011_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: (a) Isometric tensor network representation for [PITH_FULL_IMAGE:figures/full_fig_p013_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Sampling two-qubit observables from a one-dimensional quantum state, represented by a sequential unitary state. [PITH_FULL_IMAGE:figures/full_fig_p015_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Defect free energy as a function of θ in the surface code with Pauli-X rotation. Curves with increasing opacity represents system sizes L = 8, 16, 24. We choose a fixed aspect ratio Ly = 4Lx. The single qubit errors can be applied on every tensor without affecting the bond dimension or the isometric conditions on each tensor. The two-qubit error on primal rows is contained entirely within each row and thus… view at source ↗
Figure 7
Figure 7. Figure 7: Defect free energy as a function of θ for fixed ϕ = 0.02π, 0.04π, 0.06π, 0.08π, 0.12π, 0.16π, 0.2π. Data collapse are obtained using the ansatz ∆FX = f((θ − θc)L 1/ν) [PITH_FULL_IMAGE:figures/full_fig_p018_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Entanglement entropy as a function of contraction step [PITH_FULL_IMAGE:figures/full_fig_p019_8.png] view at source ↗
read the original abstract

The maximum likelihood (ML) decoder in the two-dimensional surface code with generic unitary errors is governed by a statistical mechanics model with complex weights, which can be simulated via (1+1)D transfer matrix contraction. Information loss with an increasing error rate manifests as a ferromagnetic-to-paramagnetic transition in the contraction dynamics. In this work, we establish entanglement as a separate obstruction to decoding; it can undergo a transition from area- to volume-law scaling in the transfer matrix contraction with increasing unitary error rate. In particular, the volume-law entanglement can coexist with ferromagnetic order, giving rise to a phase in which the encoded information is retained yet is effectively undecodable. We numerically simulate the ML decoding in the surface code subject to both single- and two-qubit Pauli-X rotations and obtain a phase diagram that contains a ferromagnetic area-law, a paramagnetic volume-law, and a potential ferromagnetic volume-law phase. We further show that, starting from the paramagnetic volume-law phase, tilting the single-qubit rotation away from the X-axis couples the stat-mech models for X and Z errors and can lead to a ferromagnetic volume-law phase in which, although Z errors remain in principle correctable, the encoded classical information is hard to recover. To perform numerical simulations, we develop an algorithm for syndrome sampling based on an isometric tensor network representation of the surface code.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

1 major / 2 minor

Summary. The paper maps maximum-likelihood decoding of the surface code under generic unitary errors to a (1+1)D transfer-matrix contraction of a complex-weight statistical mechanics model. Ferromagnetic order in the contraction signals retention of encoded information while entanglement scaling provides an independent diagnostic; the central claim is the existence of a ferromagnetic volume-law phase in which information is retained yet effectively undecodable. Numerical contractions for single- and two-qubit X-rotations produce a phase diagram containing ferromagnetic area-law, paramagnetic volume-law, and a potential ferromagnetic volume-law regime. Tilting the rotation axis couples the X and Z models and can induce the ferromagnetic volume-law phase. An isometric tensor-network construction is introduced for efficient syndrome sampling.

Significance. If the central mapping and numerical observations hold, the work supplies a concrete, simulable diagnostic that separates information retention from decodability via entanglement scaling. The explicit isometric tensor-network representation of syndrome sampling is a technical strength that enables direct extraction of both order parameters from the same contraction. The reported coexistence phase offers a falsifiable prediction for coherent-error models and may inform decoder design beyond standard Pauli ML thresholds.

major comments (1)
  1. [Numerical results section] The numerical phase diagram (section describing results for X-rotations) asserts the existence of a ferromagnetic volume-law regime, yet no system sizes, bond dimensions for the transfer-matrix contraction, or finite-size scaling analysis are reported. Without these details the distinction between a true phase and a finite-size crossover cannot be assessed, which is load-bearing for the central claim of a distinct undecodable-yet-retaining phase.
minor comments (2)
  1. [Abstract] The abstract states that the volume-law phase is 'potential'; the main text should clarify whether this regime is observed in the contractions or remains a conjecture based on the order-parameter definitions.
  2. [Model definition] Notation for the complex weights in the statistical mechanics model and the precise definition of the entanglement entropy extracted from the transfer matrix should be introduced with an equation reference early in the manuscript.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their careful reading of the manuscript and for highlighting the importance of finite-size analysis in supporting the central claim. We address the major comment below.

read point-by-point responses
  1. Referee: [Numerical results section] The numerical phase diagram (section describing results for X-rotations) asserts the existence of a ferromagnetic volume-law regime, yet no system sizes, bond dimensions for the transfer-matrix contraction, or finite-size scaling analysis are reported. Without these details the distinction between a true phase and a finite-size crossover cannot be assessed, which is load-bearing for the central claim of a distinct undecodable-yet-retaining phase.

    Authors: We agree that the distinction between a true phase and a finite-size crossover requires explicit reporting of system sizes, bond dimensions, and scaling analysis. The current manuscript presents numerical contractions but does not include these details in the main text or supplementary material. In the revised version we will add the lattice sizes employed (up to L=16), the bond dimensions used in the transfer-matrix contractions (χ up to 128), and a finite-size scaling collapse of the order parameters and entanglement entropy to locate the phase boundaries and confirm the volume-law ferromagnetic regime. revision: yes

Circularity Check

0 steps flagged

No significant circularity identified

full rationale

The paper constructs an explicit isometric tensor-network representation of the surface code, derives the (1+1)D transfer-matrix contraction of the complex-weight statistical mechanics model directly from that representation, and extracts both ferromagnetic order and entanglement scaling via numerical contraction of the same object. No load-bearing step reduces by definition, by fitting a parameter then relabeling the output as a prediction, or by a self-citation chain; the phase diagram is obtained from direct simulation under the paper's own mapping.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The mapping from ML decoding to a complex-weight statistical mechanics model and the use of transfer-matrix contraction are taken as given; no explicit free parameters or invented entities are stated in the abstract.

axioms (1)
  • domain assumption The maximum-likelihood decoder for unitary errors is exactly equivalent to the partition function of a (1+1)D statistical mechanics model with complex weights.
    Invoked in abstract paragraph 1 to justify the transfer-matrix simulation.

pith-pipeline@v0.9.0 · 5773 in / 1274 out tokens · 41162 ms · 2026-05-23T17:13:13.116900+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Forward citations

Cited by 2 Pith papers

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Decoding coherent errors in toric codes on honeycomb and square lattices: duality to Majorana monitored dynamics and symmetry classes

    cond-mat.stat-mech 2026-04 unverdicted novelty 8.0

    Toric code decodability under coherent X/Z errors is dual to Majorana monitored dynamics whose symmetry class (D or DIII) dictates whether the generic transition is a measurement-induced entanglement transition or a t...

  2. Coherent error induced phase transition

    quant-ph 2025-05 unverdicted novelty 5.0

    Coherent unitary errors on stabilizer codes trigger a phase transition at critical rate pc, below which the syndrome state keeps the original logical information and above which it shifts to a different logical state.

Reference graph

Works this paper leans on

55 extracted references · 55 canonical work pages · cited by 2 Pith papers · 2 internal anchors

  1. [1]

    A. Y. Kitaev, Fault-tolerant quantum computation by anyons, Annals of physics303, 2 (2003)

  2. [2]

    Dennis, A

    E. Dennis, A. Kitaev, A. Landahl, and J. Preskill, Topological quantum memory, Journal of Mathematical Physics 43, 4452 (2002)

  3. [3]

    Semeghini, H

    G. Semeghini, H. Levine, A. Keesling, S. Ebadi, T. T. Wang, D. Bluvstein, R. Verresen, H. Pichler, M. Kali- nowski, R. Samajdar,et al., Probing topological spin liq- uidsonaprogrammablequantumsimulator,Science 374, 1242 (2021)

  4. [4]

    Bluvstein, H

    D. Bluvstein, H. Levine, G. Semeghini, T. T. Wang, S. Ebadi, M. Kalinowski, A. Keesling, N. Maskara, H. Pichler, M. Greiner,et al., A quantum processor based on coherent transport of entangled atom arrays, Nature 604, 451 (2022)

  5. [5]

    Satzinger, Y.-J

    K. Satzinger, Y.-J. Liu, A. Smith, C. Knapp, M. New- man, C. Jones, Z. Chen, C. Quintana, X. Mi, A. Dunsworth, et al. , Realizing topologically ordered states on a quantum processor, Science374, 1237 (2021)

  6. [6]

    Suppressingquantum errorsby scaling asurface codelog- ical qubit, Nature614, 676 (2023)

  7. [7]

    Andersen, Y

    T. Andersen, Y. Lensky, K. Kechedzhi, I. Drozdov, A. Bengtsson, S. Hong, A. Morvan, X. Mi, A. Opremcak, E.-A. Kim, et al., Observation of non-abelian exchange statistics on a superconducting processor, Bulletin of the American Physical Society68 (2023)

  8. [8]

    Bluvstein, S

    D. Bluvstein, S. J. Evered, A. A. Geim, S. H. Li, H. Zhou, T. Manovitz, S. Ebadi, M. Cain, M. Kalinowski, D. Hangleiter,et al., Logical quantum processor based on reconfigurable atom arrays, Nature , 1 (2023)

  9. [9]

    B. M. Terhal, Quantum error correction for quantum memories, Reviews of Modern Physics87, 307 (2015)

  10. [10]

    Bravyi, M

    S. Bravyi, M. Englbrecht, R. König, and N. Peard, Cor- recting coherent errors with surface codes, npj Quantum Information 4, 55 (2018)

  11. [11]

    A. J. Ferris and D. Poulin, Tensor networks and quan- tum error correction, Physical review letters113, 030501 (2014)

  12. [12]

    Bravyi, M

    S. Bravyi, M. Suchara, and A. Vargo, Efficient algorithms for maximum likelihood decoding in the surface code, Physical Review A90, 032326 (2014)

  13. [13]

    D. K. Tuckett, S. D. Bartlett, and S. T. Flammia, Ultra- high error threshold for surface codes with biased noise, Physical review letters120, 050505 (2018)

  14. [14]

    D. K. Tuckett, A. S. Darmawan, C. T. Chubb, S. Bravyi, S.D.Bartlett,andS.T.Flammia,Tailoringsurfacecodes for highly biased noise, Physical Review X 9, 041031 (2019). 6

  15. [15]

    C. T. Chubb and S. T. Flammia, Statistical mechanical models for quantum codes with correlated noise, Annales de l’Institut Henri Poincaré D8, 269 (2021)

  16. [16]

    C. T. Chubb, General tensor network decoding of 2d pauli codes, arXiv preprint arXiv:2101.04125 (2021)

  17. [17]

    Behrends, F

    J. Behrends, F. Venn, and B. Béri, Surface codes, quan- tum circuits, and entanglement phases, Physical Review Research 6, 013137 (2024)

  18. [18]

    Y. Li, X. Chen, and M. P. Fisher, Quantum zeno ef- fect and the many-body entanglement transition, Physi- cal Review B98, 205136 (2018)

  19. [19]

    Skinner, J

    B. Skinner, J. Ruhman, and A. Nahum, Measurement- induced phase transitions in the dynamics of entangle- ment, Physical Review X9, 031009 (2019)

  20. [20]

    Y. Li, X. Chen, and M. P. Fisher, Measurement-driven entanglement transition in hybrid quantum circuits, Physical Review B100, 134306 (2019)

  21. [21]

    S. Choi, Y. Bao, X.-L. Qi, and E. Altman, Quantum er- ror correction in scrambling dynamics and measurement- induced phase transition, Physical Review Letters125, 030505 (2020)

  22. [22]

    M. J. Gullans and D. A. Huse, Dynamical purifica- tion phase transition induced by quantum measurements, Physical Review X10, 041020 (2020)

  23. [23]

    A. C. Potter and R. Vasseur, Entanglement dynamics in hybrid quantum circuits, in Entanglement in Spin Chains: From Theory to Quantum Technology Applica- tions (Springer, 2022) pp. 211–249

  24. [24]

    M. P. Fisher, V. Khemani, A. Nahum, and S. Vijay, Random quantum circuits, Annual Review of Condensed Matter Physics 14, 335 (2023)

  25. [25]

    Sang and T

    S. Sang and T. H. Hsieh, Measurement-protected quan- tum phases, Physical Review Research3, 023200 (2021)

  26. [26]

    Alberton, M

    O. Alberton, M. Buchhold, and S. Diehl, Entanglement transition in a monitored free-fermion chain: From ex- tended criticality to area law, Physical Review Letters 126, 170602 (2021)

  27. [27]

    M. P. Zaletel and F. Pollmann, Isometric tensor network states in two dimensions, Physical review letters 124, 037201 (2020)

  28. [28]

    C. Wang, J. Harrington, and J. Preskill, Confinement- higgs transition in a disordered gauge theory and the ac- curacy threshold for quantum memory, Annals of Physics 303, 31 (2003)

  29. [29]

    H. G. Katzgraber, H. Bombín, and M. A. Martin- Delgado, Error threshold for color codes and random three-body ising models, Physical review letters 103, 090501 (2009)

  30. [30]

    Bombin, R

    H. Bombin, R. S. Andrist, M. Ohzeki, H. G. Katzgraber, and M. A. Martin-Delgado, Strong resilience of topologi- cal codes to depolarization, Physical Review X2, 021004 (2012)

  31. [31]

    Kubica, M

    A. Kubica, M. E. Beverland, F. Brandão, J. Preskill, and K. M. Svore, Three-dimensional color code thresh- olds via statistical-mechanical mapping, Physical review letters 120, 180501 (2018)

  32. [32]

    H. Song, J. Schönmeier-Kromer, K. Liu, O. Viyuela, L. Pollet, and M. A. Martin-Delgado, Optimal thresh- olds for fracton codes and random spin models with sub- system symmetry, Physical Review Letters129, 230502 (2022)

  33. [33]

    F. Venn, J. Behrends, and B. Béri, Coherent error thresh- old for surface codes from majorana delocalization, arXiv preprint arXiv:2211.00655 (2022)

  34. [34]

    Verstraete and J

    F. Verstraete and J. I. Cirac, Matrix product states represent ground states faithfully, Physical Review B—Condensed Matter and Materials Physics73, 094423 (2006)

  35. [35]

    Schuch, M

    N. Schuch, M. M. Wolf, F. Verstraete, and J. I. Cirac, En- tropy scaling and simulability by matrix product states, Physical review letters100, 030504 (2008)

  36. [36]

    See supplementary online material for details

  37. [37]

    Vidal, Efficient simulation of one-dimensional quan- tum many-body systems, Physical review letters 93, 040502 (2004)

    G. Vidal, Efficient simulation of one-dimensional quan- tum many-body systems, Physical review letters 93, 040502 (2004)

  38. [38]

    Y. Bao, S. Choi, and E. Altman, Symmetry enriched phases of quantum circuits, Annals of Physics 435, 168618 (2021)

  39. [39]

    C.-M. Jian, B. Bauer, A. Keselman, and A. W. Lud- wig, Criticality and entanglement in nonunitary quantum circuits and tensor networks of noninteracting fermions, Physical Review B106, 134206 (2022)

  40. [40]

    M. Fava, L. Piroli, T. Swann, D. Bernard, and A. Nahum, Nonlinear sigma models for monitored dynamics of free fermions, Physical Review X13, 041045 (2023)

  41. [41]

    C.-M. Jian, H. Shapourian, B. Bauer, and A. W. Ludwig, Measurement-induced entanglement transi- tions in quantum circuits of non-interacting fermions: Born-rule versus forced measurements, arXiv preprint arXiv:2302.09094 (2023)

  42. [42]

    Li and M

    Y. Li and M. P. Fisher, Decodable hybrid dynamics of open quantum systems with z 2 symmetry, Physical Re- view B 108, 214302 (2023)

  43. [43]

    M. A. Nielsen and I. L. Chuang,Quantum computation and quantum information (Cambridge university press, 2010)

  44. [44]

    Y. Bao, R. Fan, A. Vishwanath, and E. Altman, Mixed- state topological order and the errorfield double formu- lation of decoherence-induced transitions, arXiv preprint arXiv:2301.05687 (2023)

  45. [45]

    R. Fan, Y. Bao, E. Altman, and A. Vishwanath, Diag- nostics of mixed-state topological order and breakdown of quantum memory, arXiv preprint arXiv:2301.05689 (2023)

  46. [46]

    J. Y. Lee, C.-M. Jian, and C. Xu, Quantum criticality under decoherence or weak measurement, PRX quantum 4, 030317 (2023)

  47. [47]

    Pollmann, S

    F. Pollmann, S. Mukerjee, A. M. Turner, and J. E. Moore, Theory of finite-entanglement scaling at one- dimensional quantum critical points, Physical review let- ters 102, 255701 (2009)

  48. [48]

    Behrends and B

    J. Behrends and B. Béri, Statistical mechanical mapping and maximum-likelihood thresholds for the surface code under generic single-qubit coherent errors, arXiv preprint arXiv:2410.22436 (2024)

  49. [49]

    Petravick, B

    D. Petravick, B. Kramer, D. Olson, M. Livny, A. Roy, P. Avery, K. Blackburn, T. Wenaus, F. Würthwein, I. Foster, et al. , The open science grid, in Journal of Physics: Conference Series , Vol. 78 (IOP Publishing,

  50. [50]

    Sfiligoi, D

    I. Sfiligoi, D. C. Bradley, B. Holzman, P. Mhashilkar, S. Padhi, and F. Wurthwein, The pilot way to grid re- sources using glideinwms, in 2009 WRI World congress on computer science and information engineering , Vol. 2 (IEEE, 2009) pp. 428–432

  51. [51]

    OSG, Ospool, OSG 10.21231/906P-4D78 (2006)

  52. [52]

    Phases of decodability in the surface code with unitary errors

    OSG, Open science data federation, OSG 10.21231/0KVZ-VE57 (2015). Supplementary Material for “Phases of decodability in the surface code with unitary errors” Yimu Bao1 and Sajant Anand2 1Kavli Institute for Theoretical Physics, University of California, Santa Barbara, CA 93106, USA 2Department of Physics, University of California, Berkeley, CA 94720, USA ...

  53. [53]

    Due to the isometric structure of the network, this RDM can be found efficiently as rows above the ones we sample trace to the identity

    Sample the stabilizer contained on either two (boundary) or three (bulk) rows by finding the RDM for the involved rows. Due to the isometric structure of the network, this RDM can be found efficiently as rows above the ones we sample trace to the identity

  54. [54]

    We apply the projectors within the rows being sampled and represent the result as an MPS

    Project|Ψ⟩into the eigenspace of the stabilizer with the desired measurement outcome. We apply the projectors within the rows being sampled and represent the result as an MPS

  55. [55]

    Fermionic isometric tensor network states in two dimensions,

    Measure out qubits in the computational basis after all stabilizers in which they are involved are measured. This ensures that the number of qubits in the MPS does now grow as one continues to sample As we proceed row-by-row upwards sampling syndromes, the bond dimension of the MPS representing the projected rows will grow and require truncation. Thus, th...