---
title: 3D Fermionic Toric Code
url: https://www.emergentmind.com/topics/3d-fermionic-toric-code
type: topic
---

# 3D Fermionic Toric Code

The **3D fermionic toric code** is a \(3{+}1\)-dimensional \(\mathbb Z_2\) topological order in which the deconfined point excitation is an emergent fermion and the flux sector is loop-like rather than anyonic. In its canonical form it is a variant of the usual 3D toric code, defined on a cubic lattice with a single qubit on each edge and commuting vertex and plaquette terms, but with operator algebra chosen so that the point-like charge has fermionic self-statistics rather than bosonic self-statistics [2503.02928]. The phase has also been described as the ordinary fermionic \(\mathbb Z_2\) gauge theory with fermionic charges and bosonic \(\mathbb Z_2\) flux loops, distinguishing it from an anomalous variant with fermionic loops [2110.14654]. Recent work places the 3D fermionic toric code at the intersection of higher-form symmetry, finite-temperature topological order, symmetry enrichment, and Floquet quantum error correction [2503.02928][1511.02563][2602.12685].

## 1. Canonical formulation as a 3D \(\mathbb Z_2\) gauge theory

A standard presentation uses a cubic lattice with periodic boundary conditions, i.e. a 3-torus, and a single qubit on each edge. The Hamiltonian has the same commuting-projector structure as the conventional 3D toric code,
\[
H = - \sum_v A_v - \sum_p B_p,
\]
with mutually commuting vertex terms \(A_v\) and plaquette terms \(B_p\) [2503.02928]. The stabilizers satisfy the same global relations as in the bosonic 3D toric code,
\[
\prod_v A_v = 1, \qquad \prod_{p \in c} B_p = 1,
\]
where the second product is over the six plaquettes of an elementary cube \(c\) [2503.02928].

The shared stabilizer skeleton obscures the essential distinction. The local form of \(A_v\), \(B_p\), and of the short string operators is described as a deformation of the usual 3D toric-code operators, arranged so that the point-like violations of \(A_v\) are fermions rather than bosons [2503.02928]. A short string operator \(W_e\) on an edge \(e\) anticommutes with the two adjacent vertex terms and commutes with all other stabilizers, so products of \(W_e\) along a path create and move point-like excitations at the endpoints. Flux excitations arise from violations of \(B_p\), and membrane operators
\[
M_\Sigma = \prod_{e \in \Sigma} X_e
\]
create loop excitations on the boundary of a surface \(\Sigma\) [2503.02928].

This formulation makes the 3D fermionic toric code a qubit Hamiltonian with intrinsically fermionic emergent content. The microscopic Hilbert space is bosonic, but the deconfined charge sector is fermionic, and this distinction controls both the braiding structure and the higher-form symmetry of the phase [2503.02928].

## 2. Excitations, braiding, and loop self-statistics

The excitation content consists of point-like charges and loop-like fluxes. Vertex violations \(A_v=-1\) are point-like and are created in pairs at the ends of string operators built from the short string pieces \(W_e\). Plaquette violations \(B_p=-1\) form closed loops on the dual lattice because \(\prod_{p\in c} B_p = 1\) forbids isolated endpoints of flux lines [2503.02928]. As in ordinary \(\mathbb Z_2\) gauge theory, taking a charge around a flux loop yields a minus sign, so the charge–loop mutual statistics is nontrivial [2503.02928].

A crucial clarification is that the ordinary 3D fermionic toric code does **not** assign fermionic self-statistics to the flux loop. The phase analyzed as the ordinary 3D fermionic toric code is the **FcBl** phase: fermionic charge, bosonic loop. By contrast, a distinct **FcFl** phase has fermionic charges and fermionic loops, and that phase is anomalous rather than a stand-alone \(3{+}1\)-dimensional bosonic lattice topological order [2110.14654].

| Phase | Point excitation | Flux loop |
|---|---|---|
| **BcBl** | bosonic charge | bosonic loop |
| **FcBl** | fermionic charge | bosonic loop |
| **FcFl** | fermionic charge | fermionic loop |

The distinction between **FcBl** and **FcFl** is diagnosed by a loop self-statistics invariant \(\mu=\pm1\), defined from a 36-step membrane process on a tetrahedral geometry. The ordinary 3D fermionic toric code has \(\mu=+1\), while the anomalous fermionic-loop phase has \(\mu=-1\) [2110.14654]. The same work argues that **FcFl** can only exist at the boundary of a non-trivial \(4{+}1\)-dimensional invertible bosonic phase with action
\[
S=\frac{1}{2}\int w_2 w_3,
\]
so the phrase “3D fermionic toric code” conventionally refers to **FcBl**, not to the anomalous fermionic-loop theory [2110.14654].

This corrects a common compression of terminology. In \(3{+}1\) dimensions, “fermionic” may refer either to the charge sector or to loop self-statistics; the ordinary fermionic toric code is fermionic in the first sense and bosonic in the second [2110.14654].

## 3. Anomalous 2-form symmetry and finite-temperature quantum topological order

The 3D fermionic toric code has a distinguished higher-form symmetry generated by closed Wilson loops of the emergent fermion. For a contractible closed path \(C\), the associated string operator \(W(C)\) commutes with the Hamiltonian and detects the parity of linked flux loops; in the loop-less sector it acts trivially, so these operators define a \(2\)-form symmetry [2503.02928]. The key property is that this \(2\)-form symmetry is **anomalous** when viewed in a purely bosonic \(3\)-dimensional setting. The anomaly disappears if physical fermions are available, because physical fermions can be bound to the emergent fermions and condensed [2503.02928].

This anomalous \(2\)-form symmetry is the mechanism behind the phase’s finite-temperature behavior. A 2025 analysis identified the 3D fermionic toric code as the first explicit example of a three-dimensional local Hamiltonian system whose equilibrium thermal states exhibit quantum topological order at sufficiently small but nonzero temperature [2503.02928]. The argument constructs a quasi-local channel that removes all loop excitations shorter than \(\ln(L)^2\), yielding a cleaned state close to a loop-less state \(\rho_\varnothing\). The fidelity bound is
\[
F(C(\rho_{\rm fTC}), \rho_\varnothing) \ge 1 - L^3 e^{-(2\beta - \ln 5)\ln(L)^2},
\]
so for
\[
\beta > \frac{\ln 5}{2}, \qquad T_0 = \frac{2}{\ln 5} \approx 1.24,
\]
the cleaned thermal state approaches \(\rho_\varnothing\) faster than any inverse polynomial in \(L\) [2503.02928].

The loop-less state realizes the anomalous \(2\)-form symmetry as a strong symmetry, and results on anomalous higher-form symmetries then imply that no short-range entangled mixed state can approximate it. Pulling this conclusion back through the quasi-local channel gives the main theorem: below \(T_0\), the fidelity between the thermal state \(\rho_{\rm fTC}\) and any short-range entangled state decays as \(O(L^{-\infty})\), equivalently the trace distance approaches \(1\) up to \(O(L^{-\infty})\) corrections [2503.02928].

The contrast with the bosonic 3D toric code is sharp. There the corresponding \(2\)-form symmetry is anomaly-free, and Hastings’ finite-temperature short-range-entanglement result applies. In the fermionic toric code, low-temperature thermal states remain long-range entangled precisely because the fermionic Wilson-loop symmetry is anomalous [2503.02928].

## 4. Symmetry enrichment, loop fractionalization, and fermionic SPT connections

When a global symmetry \(G\) is imposed, the 3D fermionic toric code becomes a symmetry-enriched \(\mathbb Z_2\) topological order with fermionic charges. The point-like charge \(e\) carries symmetry fractionalization classified by
\[
[w_e] \in H^2(G,\mathbb Z_2),
\]
exactly as in two-dimensional anyon fractionalization. Flux loops \(m\) introduce two additional layers of structure: loop–membrane intersection fractionalization
\[
[w_m] \in H^2(G,\mathbb Z_2),
\]
and, when \([w_e]=[w_m]=1\), an intrinsic loop fractionalization class
\[
[\nu] \in H^3(G,\mathbb Z_2),
\]
interpreted as the edge data of a \(2\)-dimensional SPT phase living on the loop [1511.02563].

This framework is directly tied to fermionic symmetry-protected phases. Gauging fermion parity in a \(3\)-dimensional fermionic SPT produces a \(\mathbb Z_2\) gauge theory with fermionic gauge charges, i.e. a 3D fermionic toric code. In that interpretation, the Gu–Wen supercohomology datum
\[
n_3 \in H^3(G,\mathbb Z_2)
\]
is proposed to be the loop fractionalization class \([\nu]\) of the fermion-parity flux loop in the gauged theory [1511.02563].

The classification is only partial because \([w_m]\) and \([\nu]\) are not fully independent when \([w_m]\neq 1\); the paper emphasizes equivalences and ambiguities in that regime [1511.02563]. Even so, several structural consequences follow. For \(G=\mathbb Z_2\), candidate nontrivial loop fractionalization patterns in the fermionic toric code are argued to be anomalous, leading to the conclusion that there is no nontrivial \(3\)-dimensional \(\mathbb Z_2\) fermionic SPT of that type. For \(G=\mathbb Z_2^T\), by contrast, a \(4\)-dimensional layer construction provides evidence for a nontrivial interacting \(3\)-dimensional fermionic SPT whose gauged form is a symmetry-enriched fermionic toric code with nontrivial \([\nu]\) on the flux loops [1511.02563].

Within this viewpoint, the 3D fermionic toric code is not merely “\(\mathbb Z_2\) gauge theory with a fermionic charge.” It is also the natural gauged endpoint of a class of fermionic SPT constructions, with loop fractionalization providing the bridge between higher-form topological order and fermionic symmetry protection [1511.02563].

## 5. Floquet realizations, logical structure, and monitored dynamics

The 3D fermionic toric code has a particularly important operational realization in Floquet quantum error correction. A 2023 construction introduced a 3D Floquet fermionic toric code by extending the 2D honeycomb Floquet-code framework to a trivalent lattice in three dimensions, with instantaneous stabilizer codes having the same topological order as the 3D fermionic toric code [2307.13668]. That work framed the construction in terms of condensation of topological excitations and “rewinding” of measurement schedules, using periodic sequences of non-commuting local measurements to engineer desired instantaneous stabilizer groups [2307.13668].

A subsequent 2026 construction made this program explicit for the full logical code space. It identifies a 3D Kekulé–Kitaev lattice, a tricoordinated and 3-edge-colored geometry in which deleting any one edge color yields a two-color subgraph that decomposes into short, closed loops rather than homologically nontrivial chains [2602.12685]. This loop property prevents sequential color measurements from collapsing logical information. On a 3-torus, the corresponding \(3\)D fermionic toric-code phase encodes three logical qubits, with round-independent inner logical line operators and round-dependent membrane representatives for the complementary logicals [2602.12685].

The basic backbone is the three-color cycle
\[
z \rightarrow x \rightarrow y \rightarrow z,
\]
implemented using only two-body Pauli measurements. However, on the 3D Kekulé–Kitaev lattice a simple 3-round color cycle does not expose the full plaquette-syndrome set, because some plaquettes involve all three colors and some syndrome information remains in the \(p_2,p_4,p_6\) sector [2602.12685]. The full schedule is therefore extended to
\[
(z,\ x,\ y,\ z_{\rm intra}^{\rm even},\ x,\ z_{p_2},\ y_{p_2},\ x,\ z_{\rm inter}^{\rm even},\ y),
\]
which reconstructs all plaquette stabilizer eigenvalues without disturbing the logical subspace [2602.12685]. The same work reports that the relevant logical-operator group remains invariant over the entire 10-round sequence, so the protocol implements a trivial logical automorphism while preserving all three logical qubits [2602.12685].

The same lattice geometry also supports a family of monitored Kitaev models with random measurements of the non-commuting bond parities. In the measurement-probability simplex \((p_x,p_y,p_z)\), the corner regions near \((1,0,0)\), \((0,1,0)\), and \((0,0,1)\) exhibit area-law entanglement, while a central critical region shows
\[
S(L/2) \propto L^{d-1}\log L
\]
in \(d=3\) [2602.12685]. The absence of edge critical points on the Kekulé–Kitaev lattice is again tied to the finite-loop property. This suggests that Floquet protocols preserving logical information are closely linked to trajectories in measurement space that remain within area-law topological regimes rather than crossing measurement-induced criticality [2602.12685].

## 6. Formal frameworks, antecedents, and research frontiers

The 3D fermionic toric code inherits much of its conceptual vocabulary from \(2{+}1\)-dimensional fermionic topological order. The original fermionic toric-code lattice model was introduced in 2013 as an exactly soluble \(2{+}1\)-dimensional fermionic version of the toric code, built from \(\mathbb Z_2\)-graded fusion rules and fermionic associativity data with \(\alpha=\pm i\), and described at low energies by spin Chern–Simons theories with
\[
K^{fTC}=\begin{pmatrix}0&2\\2&1\end{pmatrix}, \qquad
K^{\overline{fTC}}=\begin{pmatrix}0&2\\2&3\end{pmatrix}
\]
[1309.7032]. In \(3{+}1\) dimensions the excitation content changes from anyons to point charges and flux loops, but the same theme persists: fermionic locality is stricter than bosonic locality, and intrinsically fermionic topological data cannot be reduced to a bosonic toric-code deformation [1309.7032].

Two formal developments sharpen this point. First, fermionic MPO-injective tensor networks provide an exact PEPS/fPEPS description of the \(2{+}1\)-dimensional fermionic toric code and of fermionic twisted quantum doubles, but their explicit extension to \(3\)D remains open; the natural \(3\)D analogue would require a higher-dimensional version of fermionic MPO symmetry, i.e. membrane-level graded virtual symmetries [1609.02574]. Second, Majorana–Pauli stabilizer codes furnish an exact stabilizer realization of the \(2{+}1\)-dimensional intrinsically fermionic toric code using \(\mathbb Z_8\) Pauli operators coupled to Majorana modes, and organize it within a duality web generated by anyon condensation and gauging of bosonic or fermion-parity symmetries [2606.25048]. A plausible implication is that comparable hybrid stabilizer descriptions in \(3{+}1\) dimensions would have to encode loop condensation and higher-form braiding directly at the stabilizer-algebra level rather than by a straightforward lift of ordinary Pauli stabilizer codes.

Current frontiers therefore separate into three partially connected directions. One is the **field-theoretic and symmetry-based** direction, centered on anomalous \(2\)-form symmetry and finite-temperature long-range entanglement [2503.02928]. A second is the **symmetry-enriched and higher-categorical** direction, in which loop fractionalization and anomaly constraints organize possible fermionic \(\mathbb Z_2\) gauge theories [1511.02563]. The third is the **operational and coding-theoretic** direction, where Floquet protocols and monitored dynamics produce instantaneous 3D fermionic toric-code order using only two-body measurements while preserving a nontrivial logical subspace [2307.13668][2602.12685].

Taken together, these developments establish the 3D fermionic toric code as a central example of a \(3{+}1\)-dimensional spin topological order: simple enough to admit explicit lattice, Floquet, and symmetry-based descriptions, yet rich enough to exhibit anomalous higher-form symmetry, boundary-only variants, and low-temperature quantum topological order unavailable in the bosonic 3D toric code.

Source: https://www.emergentmind.com/topics/3d-fermionic-toric-code