Fermionic Toric Code

• The fermionic toric code is an exactly solvable 2+1D lattice model that generalizes Kitaev’s toric code using Z2-graded fusion and fermion parity constraints.
• It features a unique Hamiltonian structure with commuting projectors that enforce closed-string conditions and proper fermionic occupancy.
• Its topological order is characterized by modified anyon braiding statistics, a four-fold degenerate ground state, and tensor network formulations.

The fermionic toric code is an exactly solvable lattice model that generalizes Kitaev’s original toric code to the fermionic setting. It realizes (2+1)D intrinsically fermionic topological order by incorporating $Z_2$-graded fusion rules, nontrivial associativity relations sensitive to fermion parity, and a Hamiltonian structure that depends crucially on the anticommuting nature of fermion operators. This model cannot be fully captured by bosonic string-net or gauge theory constructions, as reflected in its fundamentally distinct locality, braiding statistics, and field-theoretic description.

1. Construction and Hamiltonian Structure

The canonical construction places qubits on lattice edges (typically a honeycomb or trivalent lattice) and spinless fermions on the vertices. The Hilbert space is

$$L_\Gamma^{\mathrm{fTC}} = \bigoplus_{I \subset V} \left( \prod_{v \in I} c_v^\dagger |0_V\rangle \otimes \bigotimes_{e \in E(\Gamma)} \mathbb{C}^2 \right)$$

where $c_v^\dagger$ is the creation operator for a fermion at vertex $v$ and the rest is the standard qubit space of the bosonic toric code (Gu et al., 2013).

The Hamiltonian is a sum of commuting projectors:

$H_{f\mathrm{TC}} = - \sum_{v} Q_v - \sum_{p} Q_p$

where the vertex term $Q_v$ enforces both the closed-string constraint and proper local fermion occupancy:

$$Q_v = \frac{1}{2}\left(1 + \prod_{i \in v} \sigma_i^z\right) \left\{ 1 - \left[ c_v^\dagger c_v - N^f(\{ \sigma_{a \in v}^z \}) \right]^2 \right\}$$

with $N^f$ a function determined by the $Z_2$-graded fusion rule:

$N^f(1,1,0) = 1,\quad N^f = 0\ \text{otherwise}$

The plaquette term flips the spins on the boundary and updates adjacent fermion decorations based on the fusion rules and local associativity relations:

$$L_\Gamma^{\mathrm{fTC}} = \bigoplus_{I \subset V} \left( \prod_{v \in I} c_v^\dagger |0_V\rangle \otimes \bigotimes_{e \in E(\Gamma)} \mathbb{C}^2 \right)$$0

where $$L_\Gamma^{\mathrm{fTC}} = \bigoplus_{I \subset V} \left( \prod_{v \in I} c_v^\dagger |0_V\rangle \otimes \bigotimes_{e \in E(\Gamma)} \mathbb{C}^2 \right)$$1 involves products of fermion creation/annihilation operators and local phase factors determined by the graded pentagon equations.

2. $$L_\Gamma^{\mathrm{fTC}} = \bigoplus_{I \subset V} \left( \prod_{v \in I} c_v^\dagger |0_V\rangle \otimes \bigotimes_{e \in E(\Gamma)} \mathbb{C}^2 \right)$$2-Graded Fusion Rules and Associativity

The key difference from the bosonic toric code is the introduction of $$L_\Gamma^{\mathrm{fTC}} = \bigoplus_{I \subset V} \left( \prod_{v \in I} c_v^\dagger |0_V\rangle \otimes \bigotimes_{e \in E(\Gamma)} \mathbb{C}^2 \right)$$3-graded fusion coefficients. In the bosonic case, $$L_\Gamma^{\mathrm{fTC}} = \bigoplus_{I \subset V} \left( \prod_{v \in I} c_v^\dagger |0_V\rangle \otimes \bigotimes_{e \in E(\Gamma)} \mathbb{C}^2 \right)$$4 (two strings fuse to vacuum). In the fermionic variant,

• $$L_\Gamma^{\mathrm{fTC}} = \bigoplus_{I \subset V} \left( \prod_{v \in I} c_v^\dagger |0_V\rangle \otimes \bigotimes_{e \in E(\Gamma)} \mathbb{C}^2 \right)$$5: two strings fusing at a vertex can optionally decorate that vertex with a fermion (odd parity).
• $$L_\Gamma^{\mathrm{fTC}} = \bigoplus_{I \subset V} \left( \prod_{v \in I} c_v^\dagger |0_V\rangle \otimes \bigotimes_{e \in E(\Gamma)} \mathbb{C}^2 \right)$$6: as in the bosonic model. All other channels vanish.

The local “F-symbols” (associativity relations) are $$L_\Gamma^{\mathrm{fTC}} = \bigoplus_{I \subset V} \left( \prod_{v \in I} c_v^\dagger |0_V\rangle \otimes \bigotimes_{e \in E(\Gamma)} \mathbb{C}^2 \right)$$7-graded and potentially pick up a phase $$L_\Gamma^{\mathrm{fTC}} = \bigoplus_{I \subset V} \left( \prod_{v \in I} c_v^\dagger |0_V\rangle \otimes \bigotimes_{e \in E(\Gamma)} \mathbb{C}^2 \right)$$8 when a fermion is transposed in a fusion tree. Specifically, consistency of the associativity relations under moves that generate a fermion loop constrain $$L_\Gamma^{\mathrm{fTC}} = \bigoplus_{I \subset V} \left( \prod_{v \in I} c_v^\dagger |0_V\rangle \otimes \bigotimes_{e \in E(\Gamma)} \mathbb{C}^2 \right)$$9 (Gu et al., 2013). These Grassmann-valued data lead to nontrivial phase factors ($c_v^\dagger$0) in processes involving fermion exchange, and more generally to “graded pentagon equations.”

3. Ground State Manifold and Topological Order

The ground state is a superposition of all closed string configurations, each decorated according to the $c_v^\dagger$1-graded fusion rule:

$c_v^\dagger$2

where the sum runs over closed string coverings and $c_v^\dagger$3 denotes the spin configuration (Gu et al., 2013). On a torus, this leads to a four-fold degenerate ground state, as in the bosonic toric code.

The $c_v^\dagger$4-matrix description is modified from the bosonic case:

$c_v^\dagger$5

The presence of an odd diagonal entry in $c_v^\dagger$6 signals fermionic topological order; no local bosonic system realizes this $c_v^\dagger$7-matrix (Gu et al., 2013).

4. Anyon Types, Fusion, and Braiding

Both bosonic and fermionic toric codes exhibit four anyon species with abelian braiding, but their statistics differ. In $c_v^\dagger$8TC:

• The fermion ($c_v^\dagger$9) is a superselection sector with odd parity.
• Two bosonic anyons ($v$0) correspond to parity-even electric and magnetic sectors.
• Fusion rules:

$v$1

• Mutual statistics:

$v$2

Thus, the $v$3 and $v$4-like excitations exhibit mutual semionic statistics, while the $v$5 sector is fermionic.

The modular $v$6 and $v$7 matrices explicitly encode the distinction via nontrivial $v$8 phase factors:

$v$9

with $H_{f\mathrm{TC}} = - \sum_{v} Q_v - \sum_{p} Q_p$0 from the associativity relations (Gu et al., 2013).

5. Mapping to Conventional Fermion and Spin Systems

Via the Jordan–Wigner transformation, the fermionic toric code can be mapped to spin or Majorana systems. For example, in the toric honeycomb model,

• Spins are mapped to pairs of Majorana operators ($H_{f\mathrm{TC}} = - \sum_{v} Q_v - \sum_{p} Q_p$1) per site,
• Complex fermions are constructed as combinations:

$H_{f\mathrm{TC}} = - \sum_{v} Q_v - \sum_{p} Q_p$2

• Plaquette operators translate to products of local fermion parity, e.g.,

$H_{f\mathrm{TC}} = - \sum_{v} Q_v - \sum_{p} Q_p$3

• The Hamiltonian becomes:

$H_{f\mathrm{TC}} = - \sum_{v} Q_v - \sum_{p} Q_p$4

Extra boundary couplings are introduced to recover the toric code’s fourfold degeneracy on a torus in the fermionic chain mapping (Liang et al., 2010).

6. Tensor Network and Category Theoretic Formulations

The ground state and excitation spectrum can be formulated via fermionic tensor networks, in particular, fPEPS and fermionic MPO-injectivity (Wille et al., 2016, Williamson et al., 2016). Here,

• Local tensors $H_{f\mathrm{TC}} = - \sum_{v} Q_v - \sum_{p} Q_p$5 carry both bosonic and Grassmann-valued fermionic indices.
• The virtual layer supports an fMPO symmetry constructed from tensors $H_{f\mathrm{TC}} = - \sum_{v} Q_v - \sum_{p} Q_p$6 with a branch-dependent structure and explicit Grassmann parity constraints.
• fMPO-injectivity is enforced via the existence of a (pseudo-)inverse tensor $H_{f\mathrm{TC}} = - \sum_{v} Q_v - \sum_{p} Q_p$7 such that $H_{f\mathrm{TC}} = - \sum_{v} Q_v - \sum_{p} Q_p$8, $H_{f\mathrm{TC}} = - \sum_{v} Q_v - \sum_{p} Q_p$9 projecting onto the fMPO-symmetric subspace.
• The topological data, including fusion and associativity ($Q_v$0-symbols), is encoded via group cohomology elements ($Q_v$1, $Q_v$2, supercohomology $Q_v$3).

These features enable a systematic classification of fermionic topological phases and direct implementation in tensor-network-based simulations (Williamson et al., 2016).

7. Physical Realizations and Applications

The fermionic toric code underpins the classification of 2+1D fermionic topological orders, provides intuition for the physics of fractional quantum Hall states, and serves as a target phase for quantum simulation platforms:

• Majorana island arrays, as in two-dimensional superconducting networks (Terhal et al., 2012), where fourth-order perturbation theory yields an effective toric code Hamiltonian.
• Platform-agnostic approaches, such as superconducting qubit arrays or custom fermion codes on quantum hardware, can implement the toric code Hamiltonian or its fermion-decorated variant (Chien et al., 2020).
• Quantum error correction using Majorana qubit codes that specifically detect and correct for odd- and even-weight errors, leveraging the fermionic code’s unique stabilizer structure (Kundu et al., 2023).

Experimental protocols for state preparation exploit the Clifford circuit structure inherent in toric code-like stabilizer Hamiltonians, with measurement-based "gluing" for arbitrary lattice geometries (Chen et al., 2022).

---

Summary Table: Fermionic Toric Code vs. Bosonic Toric Code

Aspect	Bosonic Toric Code	Fermionic Toric Code
Hilbert space	$Q_v$4	Direct sum over fermionic Fock sectors on vertices and qubits on edges
Fusion rules	$Q_v$5	$Q_v$6, $Q_v$7, graded by fermion parity
Ground state degeneracy	Four (on torus)	Four (on torus), but structure distinguished by fermionic parity
Effective field theory	$Q_v$8	$Q_v$9 (or dual with 3)
Anyon data	Abelian, mutual semions with bosonic $$Q_v = \frac{1}{2}\left(1 + \prod_{i \in v} \sigma_i^z\right) \left\{ 1 - \left[ c_v^\dagger c_v - N^f(\{ \sigma_{a \in v}^z \}) \right]^2 \right\}$$0, $$Q_v = \frac{1}{2}\left(1 + \prod_{i \in v} \sigma_i^z\right) \left\{ 1 - \left[ c_v^\dagger c_v - N^f(\{ \sigma_{a \in v}^z \}) \right]^2 \right\}$$1 matrices	Abelian, phases $$Q_v = \frac{1}{2}\left(1 + \prod_{i \in v} \sigma_i^z\right) \left\{ 1 - \left[ c_v^\dagger c_v - N^f(\{ \sigma_{a \in v}^z \}) \right]^2 \right\}$$2 in braiding; $$Q_v = \frac{1}{2}\left(1 + \prod_{i \in v} \sigma_i^z\right) \left\{ 1 - \left[ c_v^\dagger c_v - N^f(\{ \sigma_{a \in v}^z \}) \right]^2 \right\}$$3, $$Q_v = \frac{1}{2}\left(1 + \prod_{i \in v} \sigma_i^z\right) \left\{ 1 - \left[ c_v^\dagger c_v - N^f(\{ \sigma_{a \in v}^z \}) \right]^2 \right\}$$4 matrices encode fermion parity
Locality	Strict tensor product structure	Fock space, non-trivial anti-commutation, non-local in bosonic mapping
Associativity (F-symbols)	Trivial pentagon relation	Graded pentagon, $$Q_v = \frac{1}{2}\left(1 + \prod_{i \in v} \sigma_i^z\right) \left\{ 1 - \left[ c_v^\dagger c_v - N^f(\{ \sigma_{a \in v}^z \}) \right]^2 \right\}$$5 symbols with $$Q_v = \frac{1}{2}\left(1 + \prod_{i \in v} \sigma_i^z\right) \left\{ 1 - \left[ c_v^\dagger c_v - N^f(\{ \sigma_{a \in v}^z \}) \right]^2 \right\}$$6 required for consistency

---

The fermionic toric code thus exemplifies a class of exactly solvable 2+1D models with intrinsic fermionic topological order, not realizable in purely local bosonic systems. Its defining features—graded fusion, nontrivial associativity, modified $$Q_v = \frac{1}{2}\left(1 + \prod_{i \in v} \sigma_i^z\right) \left\{ 1 - \left[ c_v^\dagger c_v - N^f(\{ \sigma_{a \in v}^z \}) \right]^2 \right\}$$7-matrix, and unique braiding statistics—establish it as the paradigmatic model for investigating fermionic topologically ordered phases, their classification, and potential application in quantum simulation and error correction.