Majorana Zero Modes

A Majorana zero mode is the most economical object in the quantum-computing toolbox: half of a fermion, pinned at zero energy, storing a qubit no local instrument can read. Ordinary matter is built from Dirac fermions, whose creation operator c†c^\dagger is emphatically not their annihilation operator cc. In 1937 Ettore Majorana asked whether a relativistic fermion could instead be real, its own antiparticle, with γ=γ†\gamma=\gamma^\dagger. In a metal with fixed electron number the answer is no; but a superconductor's condensate swallows and emits charge freely, breaking particle-number conservation down to a single robust quantity — fermion parity. In that setting a Dirac electron splits cleanly into two self-conjugate Majorana halves that a topological phase drags to opposite ends of a wire. The payoff is where quantum information theory enters. One Dirac fermion carries a single occupation number, empty or filled — a qubit. Cut it into two Majorana halves separated by a length LL, and that qubit becomes nonlocal: neither end alone knows whether the fermion is present, so no local perturbation can corrupt it, and the residual error scales as e−L/ξe^{-L/\xi} — protection bought geometrically rather than by active correction. Dynamically the modes are Ising (σ\sigma) anyons: exchanging them enacts logic gates fixed by topology alone, immune to the calibration drift that plagues every conventional qubit — though braiding yields only the Clifford group, so universality demands a non-Clifford resource injected from outside. Experimentally the hunt turns on distinguishing genuine, nonlocal Majoranas from trivial Andreev states that mimic the same zero-bias signature — the question that has defined the field for a decade.

Key takeaways
  • A Majorana zero mode (MZM) is a localized, zero-energy quasiparticle in a topological superconductor whose operator is its own Hermitian conjugate, γ=γ†\gamma=\gamma^\dagger — a real fermion with γ2=1\gamma^2=1.
  • MZMs come in separated pairs; two distant modes fuse into one fermion c=12(γ1+iγ2)c=\tfrac12(\gamma_1+i\gamma_2) whose occupation is a nonlocal qubit, split only by δE∼e−L/ξ\delta E\sim e^{-L/\xi} in the canonical Kitaev pp-wave chain.
  • MZMs are Ising anyons: fusion σ×σ=1+ψ\sigma\times\sigma=1+\psi gives a 2n−12^{n-1}-dimensional protected ground space for 2n2n modes, and exchanging them realizes the braid group.
  • Braiding generates only the Clifford group {H,S,CNOT}\{H,S,\mathrm{CNOT}\} — not universal; a non-Clifford TT gate via magic-state distillation is required. Signatures (a 2e2/h≈7.748×10−5 S2e^2/h\approx7.748\times10^{-5}\ \mathrm{S} zero-bias peak, a 4π4\pi-periodic Josephson effect) are confounded by disorder-induced Andreev states.
  • Key idea: one delocalized fermion splits into two spatially separated Majorana halves; the qubit lives in the joint parity of the pair and is invisible to any local probe, so braiding rotates it fault-tolerantly through Clifford gates.
The core insight

One delocalized fermion is split into two spatially separated Majorana halves; the qubit lives in the joint parity of the pair and is invisible to any local probe, so braiding — a topological operation — rotates it fault-tolerantly through the Clifford gates.

Where this lives today

Majorana zero modes are the highest-risk, highest-reward bet in quantum hardware. The mainstream platform is the proximitized semiconductor nanowire — InAs or InSb with epitaxial aluminium, driven across VZ=Δ2+μ2V_Z=\sqrt{\Delta^2+\mu^2} by field and gate voltage — but the road has been rough: an early quantized-conductance claim (Zhang et al., Nature 2018) was retracted in 2021 after trivial Andreev states were shown to reproduce the signal, sharpening the insistence on nonlocal tests. Microsoft has since pursued a "topological gap protocol" as a statistical gatekeeper and, in early 2025, announced its "Majorana 1" processor on measurement-based tetron devices — a claim still under active scrutiny. Parallel platforms include magnetic-atom chains and iron-based superconductors (FeTex_xSe1−x_{1-x}) with near-quantized vortex-core peaks. The decisive experiments now target parity: measurement-based braiding, the 4π4\pi-periodic Josephson effect via missing odd Shapiro steps, and nonlocal parity readout. Whether any has produced a genuinely protected, braidable qubit remains open.

The physical picture

Picture an ordinary electron mode as a complex number cc; algebraically you can always write it as two real parts, c=12(γ1+iγ2)c=\tfrac12(\gamma_1+i\gamma_2). Normally those halves sit on top of one another as a single local degree of freedom. A topological superconductor tears them apart and pins one at each end of a wire, so the fermion's occupation — the stored bit — becomes a global property: the joint parity of two objects a micron apart, which no local jolt can read or flip. The qubit is not stored at either end; it is whether the two distant halves, brought together, would fuse to a filled or an empty fermion — a question with no local answer. The first figure makes this exact on Kitaev's lattice, where at the sweet spot the bulk pairs up every Majorana except two, stranded at the ends with nothing to couple to.

Figure
Kitaev chain at the sweet spot μ = 0, t = Δ: each site carries two Majoranas, every bulk pair bonds γ2ⱼ to γ2ⱼ₊₁ in H = it Σ γ2ⱼ γ2ⱼ₊₁, and the end modes γ₁, γ2N drop out — one nonlocal qubit P = iγ₁γ2N = ±1 at zero energy, split only by δE ~ Δ e^(−L/ξ)
One fermion = two Majoranas: the sweet-spot chain strands γ₁ and γ2N at zero energy H = it Σ γ2ⱼ γ2ⱼ₊₁ sweet spot μ = 0, t = Δ site 1 site 2 site 3 site N γ₁ γ2N bulk pair nonlocal qubit P = iγ₁γ2N = ±1

Braiding is choreography. Slide the endpoint Majoranas slowly around one another and the operation writes a unitary onto the ground space that depends only on the topology of the world lines — not on speed or exact path, the way a knot stays a knot however you wiggle the rope. That is why a topological gate is calibration-free, and why the single-qubit exchanges realize SS and HH but nothing finer. The second figure is the map an experimentalist actually navigates: not the abstract lattice, but the phase diagram of a real proximitized nanowire, where field and gate voltage must be tuned past a hard boundary before a single Majorana appears.

Figure
Nanowire topological phase diagram with induced gap Δ = 0.2 meV: the boundary V_Z = √(Δ² + μ²) separates a trivial region from a topological one hosting an end MZM whose ideal signature is a zero-bias conductance peak quantized at G = 2e²/h ≈ 7.748×10⁻⁵ S
Cross the boundary V_Z = √(Δ² + μ²), Δ = 0.2 meV, to make an end Majorana topological — one MZM per end trivial — no end modes V_Z = √(Δ² + μ²) G = 2e²/h ≈ 7.748×10⁻⁵ S raise field 0.2 0.4 0.6 −0.6 0 0.6 chemical potential μ (meV) Zeeman energy V_Z (meV)

Underneath both pictures is a Z2\mathbb{Z}_2 invariant. As the chemical potential tunes through the bulk gap closing, the band touches zero and the system flips between trivial (every Majorana paired within a unit cell) and topological (one unpaired half per end). Self-conjugacy γ=γ†\gamma=\gamma^\dagger is exactly why this needs a superconductor: a Majorana carries no definite charge, being an equal superposition of electron and hole, so it can exist only where the condensate makes charge ambiguous and leaves fermion parity as the sole good quantum number.

The core equations

The starting point is the algebra defining a Majorana operator: Hermitian, self-conjugate, and squaring to one, with a normalized anticommutator that replaces the usual {ci,cj†}=δij\{c_i,c_j^\dagger\}=\delta_{ij}.

γi=γi†,{γi,γj}=2δij(1)\gamma_i=\gamma_i^\dagger,\qquad \{\gamma_i,\gamma_j\}=2\delta_{ij} \qquad (1)

where: γi\gamma_i — the ii-th Hermitian Majorana operator; {⋅,⋅}\{\cdot,\cdot\} — anticommutator; δij\delta_{ij} — Kronecker delta. Setting i=ji=j forces γi2=1\gamma_i^2=1, so a Majorana has eigenvalues ±1\pm1 and no number operator of its own — it is genuinely half a degree of freedom. Because one Majorana carries no occupation number, they are assembled in pairs — two combining into one complex fermion, always possible algebraically but useful only when the halves can be pulled far apart.

c=12(γ1+iγ2),c†=12(γ1−iγ2)(2)c=\tfrac12(\gamma_1+i\gamma_2),\qquad c^\dagger=\tfrac12(\gamma_1-i\gamma_2) \qquad (2)

where: c,c†c,c^\dagger — annihilation and creation operators of one ordinary (Dirac) fermion obeying {c,c†}=1\{c,c^\dagger\}=1; γ1,γ2\gamma_1,\gamma_2 — the two Majorana halves. Inverting gives γ1=c+c†\gamma_1=c+c^\dagger, γ2=−i(c−c†)\gamma_2=-i(c-c^\dagger). The occupation of that fermion is the encoded qubit; its parity is the logical ZZ — the operator a local probe cannot access once the halves are separated.

n^=c†c=12(1+iγ1γ2),P=iγ1γ2=2n^−1=±1(3)\hat n=c^\dagger c=\tfrac12\bigl(1+i\gamma_1\gamma_2\bigr),\qquad P=i\gamma_1\gamma_2=2\hat n-1=\pm1 \qquad (3)

where: n^\hat n — occupation number in {0,1}\{0,1\}; PP — fermion parity of the pair with eigenvalues ±1\pm1. A logical error is precisely a flip of PP, which requires simultaneously touching both distant modes. The canonical microscopic model is Kitaev's chain of spinless fermions with pp-wave pairing — the minimal Hamiltonian supporting end Majoranas.

H=−∑j(t cj†cj+1+Δ cjcj+1+h.c.)−μ∑jcj†cj(4)H=-\sum_j\Bigl(t\,c_j^\dagger c_{j+1}+\Delta\,c_j c_{j+1}+\text{h.c.}\Bigr)-\mu\sum_j c_j^\dagger c_j \qquad (4)

where: tt — nearest-neighbour hopping; Δ\Delta — pp-wave pairing amplitude (taken real); μ\mu — chemical potential; cjc_j — spinless fermion on site jj. Rewriting each cjc_j in Majoranas exposes the physics. At the fine-tuned sweet spot the Hamiltonian collapses to a chain of inter-site Majorana bonds, with two operators conspicuously absent.

H=it∑j=1N−1γ2j γ2j+1(μ=0, t=Δ)(5)H=it\sum_{j=1}^{N-1}\gamma_{2j}\,\gamma_{2j+1}\qquad(\mu=0,\ t=\Delta) \qquad (5)

where: γ2j−1,γ2j\gamma_{2j-1},\gamma_{2j} — the two Majoranas on site jj. The end operators γ1\gamma_1 and γ2N\gamma_{2N} appear in no term, so they commute with HH and cost zero energy — the unpaired end modes. Away from the sweet spot the topological phase is delimited by where the bulk gap closes, which the momentum-space Bogoliubov spectrum locates.

E±(k)=±(2tcos⁡k+μ)2+4Δ2sin⁡2k(6)E_\pm(k)=\pm\sqrt{(2t\cos k+\mu)^2+4\Delta^2\sin^2 k} \qquad (6)

where: kk — crystal momentum; E±E_\pm — quasiparticle energies. The gap vanishes only at k=0k=0 (μ=−2t\mu=-2t) or k=πk=\pi (μ=2t\mu=2t); between them, ∣μ∣<2∣t∣|\mu|<2|t| with Δ≠0\Delta\neq0, the chain is topological. The modes are Ising anyons, whose fusion algebra makes the ground space degenerate and non-Abelian.

σ×σ=1+ψ,σ×ψ=σ,ψ×ψ=1(7)\sigma\times\sigma=1+\psi,\qquad \sigma\times\psi=\sigma,\qquad \psi\times\psi=1 \qquad (7)

where: σ\sigma — the Ising anyon (an MZM); ψ\psi — the fermion (Bogoliubov quasiparticle); 11 — the vacuum channel. The two outcomes of σ×σ\sigma\times\sigma are exactly the two parity states of a pair. Exchanging two modes acts on the ground space as a definite unitary — the topological square root of parity.

Bij=exp⁡ ⁣(π4γiγj)=12(1+γiγj)(8)B_{ij}=\exp\!\Bigl(\tfrac{\pi}{4}\gamma_i\gamma_j\Bigr)=\tfrac{1}{\sqrt2}\bigl(1+\gamma_i\gamma_j\bigr) \qquad (8)

where: BijB_{ij} — the unitary exchanging modes ii and jj. Conjugation gives BijγiBij†=γjB_{ij}\gamma_iB_{ij}^\dagger=\gamma_j, BijγjBij†=−γiB_{ij}\gamma_jB_{ij}^\dagger=-\gamma_i, a π/2\pi/2 rotation in the (γi,γj)(\gamma_i,\gamma_j) plane; on the qubit it is the phase gate SS. For 2n2n modes the degeneracy grows exponentially, split into two superselected parity sectors.

dim⁡H2n=2n,dim⁡H2n(fixed P)=2n−1(9)\dim\mathcal H_{2n}=2^{n},\qquad \dim\mathcal H_{2n}^{(\text{fixed }P)}=2^{n-1} \qquad (9)

where: 2n2n — number of MZMs; PP — total (superselected) fermion parity. Four Majoranas therefore encode exactly one protected logical qubit. Finally, the protection is quantified by the residual splitting as the modes' tails overlap.

δE  ∼  Egap e−L/ξcos⁡(kFL)(10)\delta E\;\sim\;E_{\mathrm{gap}}\,e^{-L/\xi}\cos(k_F L) \qquad (10)

where: δE\delta E — residual splitting of the protected pair; LL — separation of the two MZMs; ξ\xi — superconducting coherence length; kFk_F — Fermi wavevector, source of the oscillations. Doubling L/ξL/\xi squares the smallness of the error — protection bought by geometry.

Derivation AFrom the Kitaev chain to unpaired end Majoranas

The goal is to show that at the sweet spot μ=0, t=Δ>0\mu=0,\ t=\Delta>0 the open Kitaev chain leaves exactly two zero-energy Majorana modes at its ends, encoding one protected fermion, and how that exact zero mode softens to δE∼e−L/ξ\delta E\sim e^{-L/\xi} off the fine-tuned point. We rewrite the Hamiltonian in Majorana variables, read off which operators drop out, and assemble the survivors into a nonlocal qubit.

Step 1. Decompose each lattice fermion into two Majoranas, cj=12(γ2j−1+iγ2j)c_j=\tfrac12(\gamma_{2j-1}+i\gamma_{2j}), with {γa,γb}=2δab\{\gamma_a,\gamma_b\}=2\delta_{ab}. Why: the Majorana basis is the natural language of a superconductor, where pairing terms cjcj+1c_jc_{j+1} are as fundamental as hopping and both become simple bilinears γaγb\gamma_a\gamma_b.

Step 2. Substitute at the sweet spot. The chemical-potential term −μ cj†cj=−μ2(1+iγ2j−1γ2j)-\mu\,c_j^\dagger c_j=-\tfrac{\mu}{2}(1+i\gamma_{2j-1}\gamma_{2j}) vanishes at μ=0\mu=0; the combined hopping and pairing at t=Δt=\Delta reduce to −t(cj†cj+1+cjcj+1+h.c.)=it γ2jγ2j+1-t(c_j^\dagger c_{j+1}+c_jc_{j+1}+\text{h.c.})=it\,\gamma_{2j}\gamma_{2j+1}. Why: at t=Δt=\Delta the intra-site bond γ2j−1γ2j\gamma_{2j-1}\gamma_{2j} cancels and only inter-site bonds survive, dimerizing the chain across, not within, its sites.

Step 3. Collect the survivors into H=it∑j=1N−1γ2jγ2j+1H=it\sum_{j=1}^{N-1}\gamma_{2j}\gamma_{2j+1}. Every Majorana is paired except γ1\gamma_1 (left end) and γ2N\gamma_{2N} (right end), which appear in no term. Why: an operator absent from the Hamiltonian commutes with it and costs zero energy — that is the exact zero mode, localized on the terminal site with ξ→0\xi\to0.

Step 4. Assemble the two stranded ends into one delocalized fermion f=12(γ1+iγ2N)f=\tfrac12(\gamma_1+i\gamma_{2N}), with parity P=iγ1γ2N=2f†f−1=±1P=i\gamma_1\gamma_{2N}=2f^\dagger f-1=\pm1 labelling a two-fold degenerate ground space — one qubit whose logical Z=PZ=P joins opposite ends. Why: the logical operator has no single-end support, so no local perturbation can equal it or, at first order, cause a logical error.

Step 5. Leave the sweet spot into the topological bulk ∣μ∣<2∣t∣|\mu|<2|t|. The modes acquire tails ∝e−x/ξ\propto e^{-x/\xi} whose overlap splits the degeneracy by δE∼Egap e−L/ξcos⁡(kFL)\delta E\sim E_{\mathrm{gap}}\,e^{-L/\xi}\cos(k_FL), recovering an exact zero only as L→∞L\to\infty. Why: real wires are never at the sweet spot, so protection is exponential, not exact — a long enough wire is arbitrarily good. *As a consistency check: the μ\mu term maps to iγ2j−1γ2ji\gamma_{2j-1}\gamma_{2j} and vanishes at μ=0\mu=0; at t=Δt=\Delta only γ2jγ2j+1\gamma_{2j}\gamma_{2j+1} survives, γ1\gamma_1 and γ2N\gamma_{2N} drop out, and the bulk spectrum (6) closes at ∣μ∣=2∣t∣|\mu|=2|t|.*

Go deeper · Derivation B · 1 minBraiding two Majoranas gives the Clifford phase gate
Derivation BBraiding two Majoranas gives the Clifford phase gate

The goal is to derive the unitary generated by adiabatically exchanging two Majoranas γ1,γ2\gamma_1,\gamma_2 and to show it acts on the encoded qubit as the Clifford phase gate S=diag(1,i)S=\mathrm{diag}(1,i) up to a global phase — the concrete sense in which topology alone performs a quantum gate. We fix the exchange rule from the world-line geometry, find the generator that implements it by conjugation, exponentiate, and read the action in the parity basis.

Step 1. Fix the exchange rule. A counterclockwise exchange carries one mode into the other's position through a branch cut, imparting a relative sign: γ1→γ2\gamma_1\to\gamma_2, γ2→−γ1\gamma_2\to-\gamma_1. Why: doing it twice must equal a full encircling, γi→−γi\gamma_i\to-\gamma_i — the 2π2\pi monodromy that distinguishes an anyon from a boson or fermion, so the sign is the topological content.

Step 2. Seek a unitary B=eθγ1γ2B=e^{\theta\gamma_1\gamma_2} reproducing this by conjugation. Using (γ1γ2)2=−1(\gamma_1\gamma_2)^2=-1, eθγ1γ2γ1e−θγ1γ2=γ1cos⁡2θ+γ2sin⁡2θe^{\theta\gamma_1\gamma_2}\gamma_1 e^{-\theta\gamma_1\gamma_2}=\gamma_1\cos2\theta+\gamma_2\sin2\theta, so 2θ=π/22\theta=\pi/2 fixes θ=π/4\theta=\pi/4. Why: the exchange is a rigid π/2\pi/2 rotation in the (γ1,γ2)(\gamma_1,\gamma_2) plane, and only the generator γ1γ2\gamma_1\gamma_2 rotates the two operators into one another.

Step 3. Exponentiate. Because (γ1γ2)2=−1(\gamma_1\gamma_2)^2=-1 the series collapses: B12=eπγ1γ2/4=cos⁡π4+γ1γ2sin⁡π4=12(1+γ1γ2)B_{12}=e^{\pi\gamma_1\gamma_2/4}=\cos\tfrac\pi4+\gamma_1\gamma_2\sin\tfrac\pi4=\tfrac{1}{\sqrt2}(1+\gamma_1\gamma_2). Why: the closed form is a discrete, once-and-for-all matrix with no continuous knob to miscalibrate — the whole point of a topological gate.

Step 4. Read the action on the qubit. Since P=iγ1γ2P=i\gamma_1\gamma_2, B12=e−iπP/4=diag(e−iπ/4,e+iπ/4)B_{12}=e^{-i\pi P/4}=\mathrm{diag}(e^{-i\pi/4},e^{+i\pi/4}) in the parity basis; factoring the global phase, B12=e−iπ/4 diag(1,i)=e−iπ/4SB_{12}=e^{-i\pi/4}\,\mathrm{diag}(1,i)=e^{-i\pi/4}S. Why: this is exactly the Clifford phase gate S=ZS=\sqrt Z — a bare exchange is one specific Clifford generator, not an arbitrary rotation.

Step 5. Note what it does not give. Exchanging a third mode produces a Hadamard-like gate, and the two generate the single-qubit Clifford group — but never a phase finer than π/2\pi/2. Why: the braid matrices are eighth roots of unity times Paulis, so no sequence approaches the π/4\pi/4 phase of the TT gate — the gap Derivation 3 makes precise. *As a consistency check: B122=−iZB_{12}^2=-iZ (two exchanges give ZZ) and B128=1B_{12}^8=\mathbb 1 (eight exchanges close the loop), confirming the eighth-root-of-unity structure of Ising anyons.*

Go deeper · Derivation C · 1 minWhy braiding stops at the Clifford group
Derivation CWhy braiding stops at the Clifford group

The goal is to argue, via the Gottesman–Knill theorem, that Ising braiding generates exactly the Clifford group and is therefore efficiently classically simulable — hence not universal — and to identify what resource restores universality. We catalogue the braid gates, recognize them as the Clifford generators, invoke the simulability theorem, and exhibit the one gate that lies outside.

Step 1. Catalogue the braid gates: B12∼SB_{12}\sim S, B23∼HB_{23}\sim H, and inter-qubit exchanges give an entangler equivalent to CNOT\mathrm{CNOT}. Each maps encoded Paulis to encoded Paulis, g P g†=Pg\,\mathcal P\,g^\dagger=\mathcal P. Why: every braid is a Majorana bilinear rotation, which only permutes and signs the Pauli bilinears iγaγbi\gamma_a\gamma_b — it can never take a Pauli outside the Pauli group.

Step 2. Recognize the group. Unitaries normalizing the Paulis, Cn={U: UPnU†=Pn}\mathcal C_n=\{U:\,U\mathcal P_nU^\dagger=\mathcal P_n\}, are by definition the Clifford group, generated precisely by {H,S,CNOT}\{H,S,\mathrm{CNOT}\} — exactly the braid output. Why: this fixes the reachable set as Cn\mathcal C_n and no larger, closing the door on any gate that fails to conjugate Paulis to Paulis.

Step 3. Invoke Gottesman–Knill. A Clifford circuit on computational-basis inputs with Pauli measurements is classically simulable in polynomial time, tracking how Paulis transform rather than the full 2n2^n amplitude vector. Why: if braiding alone were universal it would contradict this theorem, so it cannot deliver a quantum speedup — a structural limitation, not a practical one.

Step 4. Exhibit the missing resource. The T=diag(1,eiπ/4)T=\mathrm{diag}(1,e^{i\pi/4}) gate sends X→12(X+Y)∉P1X\to\tfrac{1}{\sqrt2}(X+Y)\notin\mathcal P_1, so T∉C1T\notin\mathcal C_1 and no braid produces it. Supplementing braids with a noisy magic state ∣T⟩∝∣0⟩+eiπ/4∣1⟩|T\rangle\propto|0\rangle+e^{i\pi/4}|1\rangle, distilled and injected by teleportation, promotes the set to universal {H,S,CNOT,T}\{H,S,\mathrm{CNOT},T\}. Why: this fixes the Majorana architecture — cheap protected Clifford braids plus expensive distilled magic — placing the true resource cost in the distillation factories. *As a consistency check: T X T†=12(X+Y)T\,X\,T^\dagger=\tfrac{1}{\sqrt2}(X+Y) by direct matrix multiplication, confirming that TT maps XX outside the Pauli group and therefore lies outside the Clifford group.*

Limiting cases

The regimes below decide whether Majoranas are real, protected, and useful.

CaseResultMeaning
Long-wire / topologicalL≫ξL\gg\xiδE∼e−L/ξ→0\delta E\sim e^{-L/\xi}\to0; the parity qubit is well-defined and protected — the operating regime.
Short-wire / hybridizedL≲ξL\lesssim\xiOverlapping tails split the modes to finite energy; the degeneracy lifts and the qubit dephases — no protection.
Kitaev sweet spotμ=0, t=Δ\mu=0,\ t=\DeltaExactly decoupled end modes, zero splitting at any length; transparent but fine-tuned, ξ→0\xi\to0.
Adiabatic braidingℏ/τ≪Egap\hbar/\tau\ll E_{\mathrm{gap}}Slow gates avoid diabatic excitation out of the ground space — a speed limit on braiding.
Finite temperaturekBT≪Egapk_BT\ll E_{\mathrm{gap}}Thermal quasiparticle poisoning flips parity unless TT is well below the gap (millikelvin).

Worked examples

1 · Protection versus poisoning in a real nanowire (fully worked)

An InAs/Al wire has induced gap Δ≈0.2\Delta\approx0.2 meV and coherence length ξ≈130\xi\approx130 nm. (a) Estimate the residual splitting δE\delta E of the parity qubit at length L=1.3 μL=1.3\ \mum. (b) Compare it to thermal quasiparticle poisoning at fridge base T=20T=20 mK. Which limits the qubit?

Step 1. Set the exponent. With L=1300L=1300 nm and ξ=130\xi=130 nm, L/ξ=10L/\xi=10, so e−10=4.54×10−5e^{-10}=4.54\times10^{-5}. Taking Egap≈ΔE_{\mathrm{gap}}\approx\Delta and dropping the order-one cos⁡(kFL)\cos(k_FL), the splitting is δE∼Δ e−L/ξ=0.2 meV×4.54×10−5≈9.1 neV\delta E\sim\Delta\,e^{-L/\xi}=0.2\ \mathrm{meV}\times4.54\times10^{-5}\approx\mathbf{9.1\ neV}.

Step 2. Turn the splitting into a timescale. A nonzero δE\delta E makes the two parity states precess, dephasing on a scale t∼ℏ/δEt\sim\hbar/\delta E. With ℏ=6.58×10−16\hbar=6.58\times10^{-16} eV·s and δE=9.1×10−9\delta E=9.1\times10^{-9} eV, t∼72 μst\sim\mathbf{72\ \mu s} — the coherence budget this length buys.

Step 3. Estimate thermal poisoning. A stray quasiparticle flipping the parity must cross the gap, with Boltzmann weight e−Δ/kBTe^{-\Delta/k_BT}. Using kB=0.0862k_B=0.0862 meV/K and T=0.020T=0.020 K, Δ/kBT=0.2/0.001724≈116\Delta/k_BT=0.2/0.001724\approx116, so the weight is e−116∼10−50e^{-116}\sim10^{-50} — utterly negligible in equilibrium.

Step 4. Read the verdict. Equilibrium thermal poisoning is astronomically suppressed (e−116e^{-116}), so the practical limit is the residual hybridization splitting δE≈9.1\delta E\approx9.1 neV — the finite wire length. Doubling LL to 2.6 μ2.6\ \mum sends L/ξ→20L/\xi\to20, cutting δE\delta E to ∼0.4\sim0.4 peV. (Real devices are limited well before this by non-equilibrium quasiparticles, which do not obey the Boltzmann weight.) Answer: δE ≈ 9.1 neV (t ~ 72 μs); thermal poisoning is negligible, so wire length, not temperature, sets the protection.

2 · What a braid actually does to the qubit (fully worked)

Four Majoranas encode one qubit with Z=iγ1γ2Z=i\gamma_1\gamma_2, X=iγ2γ3X=i\gamma_2\gamma_3. (a) Which gate is the exchange B12B_{12}? (b) What is B122B_{12}^2? (c) How many exchanges return the identity?

Step 1. Identify the single braid. From equation (8), B12=e−iπZ/4=diag(e−iπ/4,e+iπ/4)B_{12}=e^{-i\pi Z/4}=\mathrm{diag}(e^{-i\pi/4},e^{+i\pi/4}) in the ZZ-eigenbasis; pulling out the global phase, B12=e−iπ/4 diag(1,i)=e−iπ/4SB_{12}=e^{-i\pi/4}\,\mathrm{diag}(1,i)=e^{-i\pi/4}S — the phase gate S=ZS=\sqrt Z.

Step 2. Square it. B122=e−iπZ/2=diag(−i,i)=−i diag(1,−1)=−i ZB_{12}^2=e^{-i\pi Z/2}=\mathrm{diag}(-i,i)=-i\,\mathrm{diag}(1,-1)=\mathbf{-i\,Z}. So two exchanges of the same pair implement logical ZZ (up to a phase) — a full braid, not just an exchange.

Step 3. Find the order. B12k=e−iπkZ/4B_{12}^k=e^{-i\pi k Z/4}, so B128=e−i2πZ=1B_{12}^8=e^{-i2\pi Z}=\mathbb 1 exactly: eight exchanges close the loop, the eighth-root-of-unity signature of Ising anyons.

Step 4. Check non-commutativity, the fingerprint of non-Abelian statistics. B23=e−iπX/4∼HB_{23}=e^{-i\pi X/4}\sim H-type, and since {X,Z}=0\{X,Z\}=0 the two braids fail to commute: B12B23≠B23B12B_{12}B_{23}\neq B_{23}B_{12}. Their products sweep out the single-qubit Clifford group — and stop there. (Quick corollary: a Kitaev chain with t=Δt=\Delta is topological iff ∣μ∣<2∣t∣|\mu|<2|t|, so μ=1.5\mu=1.5 hosts end modes while μ=2.5\mu=2.5 does not.) Answer: B₁₂ = S (up to phase); B₁₂² = −iZ; eight exchanges give the identity.

Physics of the portrait

The Majorana zero mode carries the name of a physicist who imagined its parent particle as fundamental — and vanished before it was ever found, emergent, inside a superconductor.

Origins · 1 minEttore Majorana and how this was discovered
EM
Ettore Majorana
1906–1938 (disappeared) · self-conjugate fermions and their revival in topological superconductors

Ettore Majorana was, in Enrico Fermi's estimation, a physicist of the rank of Galileo and Newton — a prodigy of Rome's Via Panisperna group who published rarely and reluctantly. In 1937, in his final paper, he found a real-valued alternative to Dirac's complex equation: a relativistic fermion that is its own antiparticle, with γ=γ†\gamma=\gamma^\dagger. He proposed it might describe the neutrino, a question still open in the search for neutrinoless double-beta decay. In March 1938 Majorana withdrew his savings, boarded a ferry between Palermo and Naples, and was never seen again — accident, suicide, or deliberate disappearance has never been settled. His self-conjugate fermion lay dormant until Alexei Kitaev's 2001 toy model showed that a one-dimensional pp-wave superconductor pins exactly such an object — now at zero energy, a Majorana zero mode — at each end. The idea detonated after the 2008 non-Abelian-anyon review and the 2010 nanowire proposals. The particle Majorana imagined as fundamental turned out to matter most as an emergent, engineered quasiparticle — a fitting afterlife for a man who left behind an equation and then left the world.

History
1937
— Ettore Majorana publishes his real-valued relativistic fermion equal to its own antiparticle
2001
— Alexei Kitaev models a 1D p-wave chain with unpaired zero-energy Majoranas at its ends
2008
— Nayak, Simon, Stern, Freedman and Das Sarma review non-Abelian anyons and topological quantum computation
2010
— Lutchyn et al. and Oreg et al. propose the semiconductor-superconductor nanowire recipe for engineered Majoranas
2012
— Mourik et al. report the first zero-bias conductance peaks in InSb nanowires

Going further

Texts and reviews. Kitaev, "Unpaired Majorana fermions in quantum wires", Phys.-Usp. 44, 131 (2001) · Nayak, Simon, Stern, Freedman, Das Sarma, "Non-Abelian anyons and topological quantum computation", Rev. Mod. Phys. 80, 1083 (2008) · Alicea, "New directions in the pursuit of Majorana fermions", Rep. Prog. Phys. 75, 076501 (2012) · Preskill, Lecture Notes on Quantum Computation, Ch. 9 · Nielsen & Chuang, Quantum Computation and Quantum Information, Ch. 10. Key relations to keep in mind: γ2=1\gamma^2=1 and {c,c†}=1\{c,c^\dagger\}=1; the sweet-spot Hamiltonian is H=it∑γ2jγ2j+1H=it\sum\gamma_{2j}\gamma_{2j+1} with γ1,γ2N\gamma_1,\gamma_{2N} absent; the topological phase boundary follows from (6); B122=−iZB_{12}^2=-iZ and B128=1B_{12}^8=\mathbb 1; and TXT†=12(X+Y)∉P1T X T^\dagger=\tfrac{1}{\sqrt2}(X+Y)\notin\mathcal P_1. The quantized conductance 2e2/h≈7.748×10−52e^2/h\approx7.748\times10^{-5} S and 4π4\pi Josephson periodicity are standard theoretical predictions; the worked-example parameters (Δ=0.2\Delta=0.2 meV, ξ=130\xi=130 nm) are representative InAs/Al values, not asserted device data.

Prerequisites. Anyons and Fractional Statistics · Fibonacci and Ising Anyons · Non-Abelian Braiding and Fusion Rules · Topological Quantum Computation. Leads to. Magic State Distillation · The Cost of Magic · Toric Code and Topological Order · Measurement-Based Quantum Computation.

Ask the tutor: Want to see exactly how four Majoranas encode one logical qubit and why braiding them stays stuck inside the Clifford group — or how magic-state distillation breaks out of it?

Is a Majorana zero mode a particle or a mode?

It is a mode — a single-particle Bogoliubov–de Gennes solution at exactly zero energy — not a particle in the everyday sense. The confusion is the name: Majorana's 1937 idea was a fundamental particle equal to its own antiparticle, whereas a superconductor hosts an emergent excitation obeying the same self-conjugacy γ=γ†\gamma=\gamma^\dagger. It is localized, immobile, and carries no definite charge or spin. Calling it a "Majorana fermion" is loose; "Majorana zero mode" is precise — half of an ordinary fermionic mode cut in two in space.

What exactly is fermion parity, and why does it matter?

Fermion parity is the operator P=iγ1γ2=(−1)n^P=i\gamma_1\gamma_2=(-1)^{\hat n} asking whether the number of fermions in a pair is even or odd. In a superconductor total particle number is not conserved — Cooper pairs are freely created and destroyed — so parity, not number, is the good quantum number. That is why Majoranas can exist: a self-conjugate mode has no definite charge and needs an environment where only parity survives. The encoded qubit is the parity of a mode pair, and a logical error is an unwanted parity flip — "quasiparticle poisoning," a stray electron changing the parity, is the dominant real-world decoherence channel.

Why does braiding give only Clifford gates and not universal computation?

Every braid is a rotation generated by a Majorana bilinear, and such rotations map Paulis to Paulis. Unitaries that do this form the Clifford group, generated by {H,S,CNOT}\{H,S,\mathrm{CNOT}\} — exactly what the exchanges produce. Gottesman–Knill then says any Clifford circuit on stabilizer inputs is classically simulable in polynomial time, so braiding alone cannot beat a classical computer. The missing ingredient is any non-Clifford gate, standardly T=diag(1,eiπ/4)T=\mathrm{diag}(1,e^{i\pi/4}); because the Ising braid matrices are only eighth roots of unity times Paulis, no braid reaches the π/4\pi/4 phase TT needs.

How is universality restored if braiding is not enough?

By injecting a non-Clifford resource from outside the topological protection. The standard route is magic-state distillation: prepare many noisy copies of ∣T⟩∝∣0⟩+eiπ/4∣1⟩|T\rangle\propto|0\rangle+e^{i\pi/4}|1\rangle with non-topological operations, then run a Clifford-only circuit that consumes them to output one nearly perfect ∣T⟩|T\rangle. Teleporting it in with Clifford braids realizes the TT gate, completing the universal set {H,S,CNOT,T}\{H,S,\mathrm{CNOT},T\}. The trade is deliberate: the fragile, expensive part is confined to the distillation factories, while the bulk of the logic runs on cheap, protected braids.

Why is a zero-bias conductance peak not proof of a Majorana?

Because trivial physics fakes it. A Majorana should give a zero-bias peak quantized to 2e2/h2e^2/h, but a smooth end potential or disorder can bind an ordinary Andreev state near zero energy that yields the same peak — even one close to 2e2/h2e^2/h over a parameter range. This ambiguity is exactly why a 2018 claim of quantized Majorana conductance was retracted in 2021. Distinguishing them requires nonlocality: correlated peaks at both ends, a 4π4\pi-periodic Josephson current (missing odd Shapiro steps), or a nonlocal parity measurement. A single local probe can never settle it.

How do measurement-based schemes braid without moving anything?

Physically dragging Majoranas around one another is hard and risks diabatic errors, so modern architectures (tetron and hexon designs) never move them. Instead they run a scripted sequence of two- and four-Majorana parity measurements via charge sensors. By the logic of measurement-based quantum computation, a well-chosen sequence of projective parity measurements teleports the encoded information through the array and enacts the same Clifford unitary a physical braid would — a "braid by measurement". This trades the control problem of transport for the readout problem of fast, high-fidelity parity measurement, the paradigm most current Majorana processors pursue.

Before this

Where this leads