We apply finite-basis duality to a statistical-mechanical model introduced by Chubb and Flammia for the surface code under spatially correlated bit-flip noise, for which their mapping gives a random-bond Ising model with both two-body edge interactions and four-body face interactions, and obtain several estimations of the error thresholds via the duality analysis in conjunction with the replica method.
We apply finite-basis duality to a statistical-mechanical model introduced by Chubb and Flammia for the surface code under spatially correlated bit-flip noise, for which their mapping gives a random-bond Ising model with both two-body edge interactions and four-body face interactions, and obtain several estimations of the error thresholds via the duality analysis in conjunction with the replica method.
We apply the single-equation duality criterion to the square-octagonal random-bond Ising model recently obtained from an exact error-edge map for a surface code with nearest-neighbor correlated errors. The calculation is performed for the minimal cell after the error-edge reduction. For the symmetric case p1=p2=p3=p, this gives pc=0.0288427147, in close agreement with the reported numerical threshold of about 3%.
Nishimori Threshold Estimation for Bayesian Inference and Zq Surface Code Decoding
In quantum error correction, the error threshold provides essential quantitative guidance for the ability to bring about fault-tolerance through decoding the effects of incoherent noise, weak measurement or inference. However, the numerical value of an error threshold is typically only accessible through large-scale numerical simulations of the underlying noise model. Here we introduce an analytical estimate of error thresholds falling into the Nishimori universality class via a Fourier--Walsh projection scheme that maps the critical point of the underlying disorder-free statistical-mechanics model to the Born-disordered Nishimori critical point. Using a minimal replica theory approach, this closed-form estimate is obtained from a projection of the exact replicated single-bond weight which we find to reproduce (within a percentage point) the known numerical thresholds of random-bond and random-plaquette Ising models / Z2 stabilizer codes in spatial dimensions d=2−5, and extends to Potts variables with q≤4. The main application of our projection scheme is to Zq surface codes, whose decoding problem maps to the disordered q-state clock model. For q≥5 the clean clock model has \textit{two} Berezinskii--Kosterlitz--Thouless transitions, which the projection maps to two Nishimori temperatures that bound an intermediate information-critical phase. The resulting threshold values not only accurately agree with recent decohered-Zq-toric-code numerics, but are found to satisfy the Gilbert--Varshamov self-dual entropy relation lnq≃Hq(T1∗)+Hq(T2∗), although no duality condition is imposed in the construction. Our approach thereby points to a deeper connection between the clean and Born-disordered models, while allowing for instant analytical estimates of error thresholds for a variety of stabilizer codes.
The surface code represents a promising candidate for fault-tolerant quantum computation due to its high error threshold and experimental accessibility with nearest-neighbor interactions. However, current exact surface code threshold analyses are based on the assumption of independent and identically distributed (i.i.d.) errors. Though there are numerical studies for the threshold with correlated error, they are only the lower bound rather than the exact value, which offers potential for higher error thresholds. Here, we establish an error-edge map, which allows for the mapping of quantum error correction to a square-octagonal random bond Ising model. We then present the maximum-likelihood threshold under a realistic noise model that combines independent single-qubit errors with correlated errors between nearest-neighbor data qubits. Our method is applicable for any ratio of nearest-neighbor correlated errors to i.i.d. errors. We present analytical constraints for the error correction threshold of surface codes, and by applying these constraints, we numerically obtain the error threshold higher than the existing values. Our error threshold is both an upper bound and achievable.
We investigate how the lattice geometry influences the error thresholds of two-dimensional toric codes. In the presence of bit- or phase-flip errors, the toric code maps onto the two-dimensional random-bond Ising model (RBIM). We determine the critical behaviour of the RBIM using replica-exchange Monte Carlo simulations. While previous studies have explored thresholds for general lattice geometries under suboptimal decoders such as minimum-weight perfect matching, the rapid development of near-optimal decoders makes resolving the ultimate, maximum-likelihood code capacity relevant. We compute these optimal error thresholds on the Nishimori line for the square, honeycomb, triangular, dice, and kagome lattices. Owing to lattice duality, the phase-flip threshold on a given lattice is equivalent to the bit-flip threshold on its dual. We find that the optimal thresholds of dual-lattice pairs display a characteristic duality-driven splitting around the self-dual square-lattice, mirroring the qualitative behaviour observed for suboptimal decoders. While the average coordination number has the strongest impact on the error threshold, our results demonstrate that the detailed arrangement of vertices and plaquettes also plays a significant role in determining its precise value.
The code-capacity threshold of a scalable quantum error correcting stabilizer code can be expressed as a thermodynamic phase transition of a corresponding random-bond Ising model. Here we study the XY and XZZX surface codes under phase-biased noise, px=py=pz/(2η) , with η≥1/2 , and total error rate p=px+py+pz . By appropriately formulating the boundary conditions, in the rotated code geometry, we find exact solutions at a special disordered point, p=1+η−12+η−1≳0.5 , for arbitrary odd code distance d , where the codes reduce to one-dimensional Ising models. The total logical failure rate is given by Pf=34−14e−2dZartanh(1/2η) , where dZ=d2 and d for the two codes respectively, is the effective code distance for pure phase-flip noise. As a consequence, for code distances d≪η , and error rates near the threshold, the XZZX code is effectively equivalent to the phase-flip correcting repetition code over d qubits. The large finite size corrections for dZ<η also make threshold extractions, from numerical calculations at moderate code distances, unreliable. We show that calculating thresholds based not only on the total logical failure rate, but also independently on the phase- and bit-flip logical failure rates, gives a more confident estimate. Using this method for the XZZX code with a tensor-network based decoder and code distances up to d≈100 , we find that the thresholds converge to a single value at moderate bias ( η=30,100 ), at an error rate above the hashing bound.
It is known that noisy topological quantum codes are related to random-bond Ising models where the order-disorder phase transition in the classical model is mapped to the error threshold of the corresponding topological code. On the other hand, there is a dual mapping between classical spin models and quantum Calderbank-Shor-Steane (CSS) states where the partition function of a classical model defined on a hypergraph H is written as an inner product of a product state and a CSS state on dual hypergraph $\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{H}$. It is then interesting to see what the interpretation is of the classical phase transition in the random-bond Ising model within the framework of the above duality, and whether such an interpretation has any connection to the error threshold of the corresponding topological CSS code. In this paper, we consider the above duality relation specifically for a two-dimensional random-bond Ising model. We show that the order parameter of this classical system is mapped to a coherence order parameter in a noisy Toric code model. In particular, a quantum phase transition from a coherent phase to a noncoherent phase occurs when the initial coherent state is affected by two sequences of bit-flip quantum channels where a quenched disorder is induced by measurement of the errors after the first channel. On the other hand, the above transition is directly related to the error threshold of the Toric code model. Accordingly, and since the noisy process can be applied to other topological CSS states, we conclude that the dual correspondence can also provide a useful tool for the study of error thresholds in different topological CSS codes.