Skip to content

Topological Superconductors

A topological superconductor is a superconducting phase whose gapped bulk Bogoliubov quasiparticles cannot be deformed to a topologically trivial superconducting reference while the relevant symmetries and bulk gap are preserved. Its bulk obstruction is expressed through a Bogoliubov–de Gennes (BdG) invariant, and a boundary or topological defect can consequently support Majorana modes.

This definition has three important qualifications.

  1. The topology belongs to a bulk phase, not to an isolated zero-bias feature.
  2. The standard classification discussed here is a classification of gapped quadratic BdG Hamiltonians. Nodes, strong interactions, crystalline symmetries, and disorder can require refined tools.
  3. A Majorana zero mode is a self-conjugate emergent quasiparticle operator. It is not evidence that a fundamental relativistic Majorana particle has been discovered.

The mathematical framework is well established. Experimental identification is not comparably settled. As of July 2026, many systems display ingredients or signatures compatible with topological superconductivity, but no solid-state platform has achieved a universally accepted demonstration of topologically protected Majorana fusion and braiding. Individual claims should therefore be reported together with their assumptions, alternative explanations, and independent checks.

This page is the canonical home for BdG topology, particle–hole redundancy, the class-D and class-DIII invariants, and Majorana boundary and vortex modes. Bogoliubov Quasiparticles owns quasiparticle branch counting and coherence factors; BCS Mean-Field Theory owns the pairing saddle and gap equation; the Kitaev Chain is the compact canonical model entry; Anyons and Braiding owns the braid-group consequences; and Topological Qubits owns complete Majorana modules, parity readout, control, metrics, and architecture evidence.

Use Superfluidity and Superconductivity to establish the generic superconducting object and evidence ledger before adding topology; this page owns the topological obstruction and its consequences.

Required background. BCS Mean-Field Theory supplies the pairing saddle, Bogoliubov Quasiparticles supplies Nambu branch counting, and Topology in Quantum Matter supplies the gapped-deformation and evidence framework.

Helpful background. Symmetry Classification Preview supplies tenfold-way language, the Kitaev Chain supplies the minimal class-D model, and the superconductivity gateway supplies the material and proximity-effect preparation.

A conventional band insulator is described by occupied electron bands separated from empty electron bands. A superconductor is different: its condensate mixes electron annihilation and creation operators, so a quasiparticle generally carries neither a definite electron number nor a definite electric charge. The natural single-particle description is therefore written in Nambu space, where particles and holes appear together.

The mean-field Hamiltonian does not conserve the electron number operator N^\hat N. It does preserve fermion parity,

P^f=(−1)N^,\hat P_f = (-1)^{\hat N},

because pairing creates or annihilates electrons in pairs. In an isolated exact system, fermion parity is a genuine superselection label. In an open device, quasiparticle poisoning can exchange a single electron with the environment and spoil the effective parity lifetime.

Topology requires a specified gap. For a fully gapped superconductor, the positive BdG spectrum obeys

En(k)≥Δqp>0E_n(\mathbf k) \geq \Delta_{\mathrm{qp}} > 0

throughout the Brillouin zone. A nodal superconductor can possess topological charges attached to points or lines of nodes, as well as protected surface flat bands or arcs, but it is not classified by simply inserting its Hamiltonian into the table for fully gapped phases.

The condensate itself need not have an exotic microscopic origin. Conventional ss-wave pairing can produce an effective topological superconductor when combined with spin–orbit coupling, Zeeman splitting, a topological-insulator surface, or an appropriately structured junction. Conversely, unconventional pairing does not automatically imply nontrivial topology. Pairing symmetry, normal-state band structure, dimensionality, and protecting symmetries must be considered together.

Let h(k)h(\mathbf k) be an N×NN\times N normal-state Bloch Hamiltonian, including orbital and spin labels, and let Δ(k)\Delta(\mathbf k) be its pairing matrix. For the Nambu spinor

Ψk=(ck1⋮ckNc−k1†⋮c−kN†),\Psi_{\mathbf k} = \begin{pmatrix} c_{\mathbf k1} \\ \vdots \\ c_{\mathbf kN} \\ c^\dagger_{-\mathbf k1} \\ \vdots \\ c^\dagger_{-\mathbf kN} \end{pmatrix},

the quadratic mean-field Hamiltonian can be written

H^MF=12∑kΨk†HBdG(k)Ψk+E0,\hat H_{\mathrm{MF}} = \frac{1}{2} \sum_{\mathbf k} \Psi_{\mathbf k}^\dagger \mathcal H_{\mathrm{BdG}}(\mathbf k) \Psi_{\mathbf k} + E_0,

with

HBdG(k)=(h(k)−μΔ(k)Δ†(k)−hT(−k)+μ).\mathcal H_{\mathrm{BdG}}(\mathbf k) = \begin{pmatrix} h(\mathbf k)-\mu & \Delta(\mathbf k) \\ \Delta^\dagger(\mathbf k) & -h^T(-\mathbf k)+\mu \end{pmatrix}.

The factor 1/21/2 compensates for Nambu doubling. One may instead sum over an independent half of momentum space with a corresponding boundary prescription. Counting every eigenvalue of the doubled matrix as an independent excitation would count the same physical degree of freedom twice.

Fermi statistics impose

Δ(k)=−ΔT(−k).\Delta(\mathbf k) = -\Delta^T(-\mathbf k).

This is the compact form of the antisymmetry of a Cooper pair under exchange. A spin-singlet pair can be even in momentum because its spin tensor is antisymmetric; a single-component or spin-polarized pair must be odd in momentum.

Under a global electron rephasing c↦eiχ/2cc\mapsto e^{i\chi/2}c, the Nambu spinor transforms as

Ψ↦eiχτz/2Ψ,\Psi \mapsto e^{i\chi\tau_z/2}\Psi,

and the pairing matrix changes by Δ↦eiχΔ\Delta\mapsto e^{i\chi}\Delta. A uniform condensate phase can therefore be moved between the off-diagonal matrix elements and the Nambu basis. Relative phases across a junction and phase winding around a vortex remain physical.

The Pauli matrices τi\tau_i below act in particle–hole space. Spin Pauli matrices are written sis_i. In a multiorbital effective model, the meaning of either set must be stated; a symbol appearing like physical spin can instead label a pseudospin doublet.

Every BdG Hamiltonian has an antiunitary particle–hole relation

CHBdG(k)C−1=−HBdG(−k).\mathcal C \mathcal H_{\mathrm{BdG}}(\mathbf k) \mathcal C^{-1} = -\mathcal H_{\mathrm{BdG}}(-\mathbf k).

In the canonical Nambu basis one may choose

C=τxK,C2=+1,\mathcal C = \tau_xK, \qquad \mathcal C^2=+1,

where KK denotes complex conjugation. Thus, if

HBdG(k)Φnk=EnkΦnk,\mathcal H_{\mathrm{BdG}}(\mathbf k) \Phi_{n\mathbf k} = E_{n\mathbf k} \Phi_{n\mathbf k},

then

τxΦnk∗\tau_x\Phi_{n\mathbf k}^*

is an eigenvector at −k-\mathbf k with energy −Enk-E_{n\mathbf k}. In second-quantized language, the quasiparticle associated with the negative-energy partner is the creation operator for the positive-energy mode:

α−E,−k=αE,k†.\alpha_{-E,-\mathbf k} = \alpha^\dagger_{E,\mathbf k}.

This relation is often called particle–hole symmetry, but it is first of all a redundancy introduced by the Nambu description. It is not ordinary charge conservation, not a transformation that maps the laboratory to a distinct state of equal positive energy, and not relativistic charge conjugation. Additional physical symmetries, such as time reversal, inversion, mirror symmetry, or spin rotation, constrain the BdG Hamiltonian independently.

At exactly zero energy, the particle–hole partner lies in the same eigenspace. For an isolated nondegenerate zero mode, its phase may be chosen so that

CΦ0=Φ0.\mathcal C\Phi_0 = \Phi_0.

The corresponding quasiparticle operator obeys

γ†=γ.\gamma^\dagger = \gamma.

This self-adjoint operator is a Majorana zero mode. Particle–hole redundancy permits such an operator, but it does not force every superconductor to have one. Existence and protection require topology, symmetry, or fine tuning.

The stable free-fermion classification asks which gapped BdG Hamiltonians can be continuously deformed into one another after adding trivial bands. Three antiunitary or unitary constraints organize the problem:

TH(k)T−1=H(−k),CH(k)C−1=−H(−k),SH(k)S−1=−H(k).\begin{aligned} \mathcal T \mathcal H(\mathbf k) \mathcal T^{-1} &= \mathcal H(-\mathbf k), \\ \mathcal C \mathcal H(\mathbf k) \mathcal C^{-1} &= -\mathcal H(-\mathbf k), \\ \mathcal S \mathcal H(\mathbf k) \mathcal S^{-1} &= -\mathcal H(\mathbf k). \end{aligned}

Here T\mathcal T is physical time reversal, C\mathcal C is BdG particle–hole conjugation, and chiral symmetry S\mathcal S can arise from their product. The most common symmetry classes for topological superconductors are:

ClassTime reversalParticle–holeChiral1D2D3D
DabsentC2=+1\mathcal C^2=+1absentZ2\mathbb Z_2Z\mathbb Z00
DIIIT2=−1\mathcal T^2=-1C2=+1\mathcal C^2=+1presentZ2\mathbb Z_2Z2\mathbb Z_2Z\mathbb Z
BDIT2=+1\mathcal T^2=+1C2=+1\mathcal C^2=+1presentZ\mathbb Z0000

Class D describes a generic superconductor without time-reversal symmetry. Class DIII describes a time-reversal-invariant spinful superconductor. A Kitaev chain with all parameters real has an additional effective time-reversal symmetry and lies in class BDI, allowing an integer winding number. Generic complex perturbations that preserve only BdG particle–hole conjugation reduce the robust distinction to the class-D parity.

The entries in the table are stable, noninteracting, internal-symmetry classifications. Interactions can identify phases that are distinct in the quadratic table; for example, the one-dimensional BDI integer reduces to a Z8\mathbb Z_8 distinction and the three-dimensional DIII integer reduces to a Z16\mathbb Z_{16} distinction under suitable symmetry-preserving interactions. Crystalline symmetry can add weak, mirror, rotation, hinge, or higher-order indices. Disorder can preserve a real-space or scattering invariant even after momentum ceases to be a good quantum number.

Whatever invariant is used, it cannot change under a continuous deformation unless at least one defining condition fails. Typically the quasiparticle gap closes:

Δqp⟶0.\Delta_{\mathrm{qp}} \longrightarrow 0.

A gap closing is therefore necessary at a generic quadratic topological phase transition, but a measured gap closing by itself does not identify which invariant changed. Topological Phase Transitions derives the ideal nanowire criterion and separates quasiparticle closings from first-order, disordered, and interacting routes.

For a translation-invariant one-dimensional class-D superconductor with C=τxK\mathcal C=\tau_xK, define at the particle–hole-invariant momenta k⋆=0,πk_\star=0,\pi

B(k⋆)=HBdG(k⋆)τx.B(k_\star) = \mathcal H_{\mathrm{BdG}}(k_\star)\tau_x.

Particle–hole conjugation makes B(k⋆)B(k_\star) antisymmetric in this basis. When its Pfaffian is real and nonzero, the Majorana number can be written

M=sgn⁡[Pf⁡B(0)Pf⁡B(π)].\mathcal M = \operatorname{sgn} \left[ \operatorname{Pf}B(0) \operatorname{Pf}B(\pi) \right].

The phase is topological for M=−1\mathcal M=-1 and trivial for M=+1\mathcal M=+1. A change of sign requires a Pfaffian to vanish, which is precisely a bulk gap closing at one of the invariant momenta. Gauge-covariant formulations or scattering-matrix invariants should be used when the convenient real Pfaffian convention is unavailable.

For the real, nearest-neighbor Kitaev-chain convention

H(k)=(−μ−2tcos⁡k)τz+2Δsin⁡k τy,\mathcal H(k) = \left( -\mu-2t\cos k \right)\tau_z + 2\Delta\sin k\,\tau_y,

the pairing term vanishes at k=0,πk=0,\pi. The invariant reduces to

M=sgn⁡[(μ+2t)(μ−2t)].\mathcal M = \operatorname{sgn} \left[ (\mu+2t)(\mu-2t) \right].

Hence

∣μ∣<2∣t∣⟹M=−1.\lvert\mu\rvert < 2\lvert t\rvert \quad\Longrightarrow\quad \mathcal M=-1.

An open chain in this regime has one exponentially localized Majorana mode at each end. At ∣μ∣=2∣t∣\lvert\mu\rvert=2\lvert t\rvert, the bulk gap closes and the localization length diverges. The Kitaev Chain entry develops the lattice Hamiltonian and solvable limits; the point here is the general correspondence between a bulk Pfaffian parity and an odd number of protected boundary Majoranas.

A normalized Majorana operator localized near a boundary has the form

γa=∑j(uajcj+uaj∗cj†),\gamma_a = \sum_j \left( u_{aj}c_j + u_{aj}^*c_j^\dagger \right),

with algebra

γa†=γa,{γa,γb}=2δab.\gamma_a^\dagger=\gamma_a, \qquad \{\gamma_a,\gamma_b\} = 2\delta_{ab}.

Two separated Majoranas form one ordinary complex fermion:

f=γL+iγR2,iγLγR=2f†f−1.f = \frac{\gamma_L+i\gamma_R}{2}, \qquad i\gamma_L\gamma_R = 2f^\dagger f-1.

The occupation nf=0,1n_f=0,1 changes the fermion parity associated with the pair. A finite overlap produces the effective Hamiltonian

H^split=iε2γLγR=ε(f†f−12).\hat H_{\mathrm{split}} = \frac{i\varepsilon}{2} \gamma_L\gamma_R = \varepsilon \left( f^\dagger f-\frac{1}{2} \right).

For a clean long wire, the splitting commonly behaves as

ε(L)∼Ae−L/ξcos⁡(kFL+φ),\varepsilon(L) \sim A e^{-L/\xi} \cos(k_FL+\varphi),

where ξ\xi is the Majorana localization length. The oscillation is not itself a topological invariant; it reflects microscopic phase accumulation and can be strongly modified by disorder, multiple bands, and smooth confinement.

Local perturbations near one end cannot directly implement the nonlocal bilinear iγLγRi\gamma_L\gamma_R when the two wavefunctions do not overlap. This is the origin of exponential protection in an idealized wire. It is not absolute protection:

  • finite overlap splits the pair;
  • a local perturbation can close the protecting gap;
  • accidental low-energy Andreev states can couple to an end mode;
  • quasiparticle poisoning changes total fermion parity;
  • measurement leads can broaden or decohere the state;
  • a single pair does not provide two usable states within a fixed total-parity sector.

Encoding a parity-conserving qubit requires at least four Majoranas or an equivalent multi-island architecture, plus controlled measurement and initialization. Braiding separated Majorana defects implements noncommuting operations only when the evolution remains adiabatic relative to the bulk gap, fast relative to poisoning, and insensitive to residual splittings. Those operations are developed in Anyons and Braiding.

BdG energy pairing, Majorana end modes, and vortex-core level comparison

Three diagnostics that must not be conflated. Particle–hole redundancy pairs every finite-energy BdG state with a state at −E-E; a topological open wire can leave self-conjugate modes γL\gamma_L and γR\gamma_R near zero with exponentially small overlap; and an ordinary vortex has a half-shifted Caroli–de Gennes–Matricon ladder, whereas an appropriate topological vortex can carry an unpaired zero mode.

A minimal continuum model for a spin-polarized chiral px+ipyp_x+ip_y superconductor is

H(k)=ξkτz+ΔkF(kxτx−kyτy),\mathcal H(\mathbf k) = \xi_{\mathbf k}\tau_z + \frac{\Delta}{k_F} \left( k_x\tau_x-k_y\tau_y \right),

where ξk=k2/(2m)−μ\xi_{\mathbf k}=k^2/(2m)-\mu. After a physically necessary ultraviolet regularization, the negative-energy BdG band has Chern number

C=12π∑En<0∫BZΩn,xy(k) d2k.C = \frac{1}{2\pi} \sum_{E_n<0} \int_{\mathrm{BZ}} \Omega_{n,xy}(\mathbf k)\,d^2k.

The weak-pairing regime μ>0\mu>0 is topological, with ∣C∣=1\lvert C\rvert=1 for one chiral species; the strong-pairing regime μ<0\mu<0 is trivial. The transition at μ=0\mu=0 closes the bulk gap. This is the Read–Green distinction: two states with the same broken symmetries can differ by their BdG topology.

A boundary of a class-D phase with Chern number CC carries net chiral Majorana central charge

c−=C2.c_- = \frac{C}{2}.

In the ideal low-temperature limit its thermal Hall response is

κxyT=Cπ2kB26h.\frac{\kappa_{xy}}{T} = C \frac{\pi^2k_B^2}{6h}.

The factor of one half relative to a chiral complex fermion is characteristic of a Majorana channel. BdG quasiparticles do not carry a conserved U(1)U(1) charge, and the electromagnetic response also includes the condensate. Thus CC fixes the ideal chiral-Majorana thermal response, but it does not fix an analogous universal quantized electric Hall conductance.

Class DIII preserves spinful time reversal with T2=−1\mathcal T^2=-1. Its protected boundaries appear in dimension-dependent forms:

  • a one-dimensional phase has a Kramers pair of Majorana end modes;
  • a two-dimensional phase has a helical Majorana edge;
  • a three-dimensional phase can have a gapless Majorana surface cone.

After flattening a three-dimensional chiral-symmetric BdG Hamiltonian into an off-diagonal unitary block q(k)q(\mathbf k), one convention for the integer winding number is

ν3=124π2∫BZd3k ϵijkTr⁡[q−1∂iq q−1∂jq q−1∂kq].\nu_3 = \frac{1}{24\pi^2} \int_{\mathrm{BZ}} d^3k\, \epsilon^{ijk} \operatorname{Tr} \left[ q^{-1}\partial_iq\, q^{-1}\partial_jq\, q^{-1}\partial_kq \right].

Sign conventions depend on the orientation and block definition, but ∣ν3∣\lvert\nu_3\rvert counts the stable surface-cone multiplicity in the minimal free description. Superfluid 3He^3\mathrm{He}-B is a benchmark class-DIII topological neutral superfluid. It demonstrates the phase structure without being an electronic superconducting material.

For a centrosymmetric, fully gapped, time-reversal-invariant superconductor, odd-parity pairing and the pattern of Fermi surfaces around time-reversal-invariant momenta can provide a useful sufficient topological criterion. It is not a substitute for determining the pairing symmetry: uncertainty about the order parameter propagates directly into uncertainty about the topological label.

A vortex carries condensate phase winding

∮∇ϕ⋅dℓ=2πn.\oint\nabla\phi\cdot d\boldsymbol\ell = 2\pi n.

The order parameter is suppressed in its core, so subgap states are expected even in a conventional superconductor. In the semiclassical regime EF≫ΔE_F\gg\Delta, an ordinary singly quantized vortex has Caroli–de Gennes–Matricon levels of order

Em≃(m+12)ω0,ω0∼Δ2EF,E_m \simeq \left( m+\frac{1}{2} \right)\omega_0, \qquad \omega_0 \sim \frac{\Delta^2}{E_F},

with integer mm in a simple spinless notation. These levels can lie very close to zero when Δ/EF\Delta/E_F is small. A near-zero vortex peak is therefore not automatically a Majorana mode.

In a two-dimensional chiral topological superconductor with odd Chern parity, an odd-vorticity defect can bind an unpaired Majorana zero mode. Spectral flow effectively removes the half-level offset for one branch, leaving a state pinned to zero by particle–hole conjugation unless it couples to another Majorana or the protecting assumptions fail. A proximitized topological-insulator surface provides another route: spin-momentum locking plus ordinary ss-wave pairing yields an effective spinless pairing structure, and a vortex can bind a Fu–Kane Majorana mode.

For a finite collection of vortices, zero modes hybridize pairwise:

H^M=i2∑a<btabγaγb.\hat H_{\mathrm{M}} = \frac{i}{2} \sum_{a<b} t_{ab}\gamma_a\gamma_b.

The amplitudes tabt_{ab} depend exponentially and oscillatory on separation in a clean system. Disorder, anisotropy, nearby surfaces, other vortex-core levels, and three-dimensional vortex-line dispersion complicate the spectrum. A credible vortex claim should therefore establish localization, particle–hole structure, robustness over a parameter region, separation from ordinary core levels, and consistency with the bulk topology.

An ideal single-band model for a Rashba nanowire with Zeeman energy VZV_Z and induced ss-wave pairing is

HNW(k)=(ℏ2k22m∗−μ)τz+αk syτz+VZsx+Δindτx.\begin{aligned} \mathcal H_{\mathrm{NW}}(k) ={}& \left( \frac{\hbar^2k^2}{2m^*}-\mu \right)\tau_z \\ &+ \alpha k\,s_y\tau_z + V_Zs_x + \Delta_{\mathrm{ind}}\tau_x. \end{aligned}

At k=0k=0, the ideal bulk gap closes when

VZ2=μ2+Δind2.V_Z^2 = \mu^2+\Delta_{\mathrm{ind}}^2.

The single-band topological regime lies on the side

VZ2>μ2+Δind2,V_Z^2 > \mu^2+\Delta_{\mathrm{ind}}^2,

provided a reopened gap survives. Spin–orbit coupling converts the spin-singlet proximity effect into an effective odd-momentum pairing channel in the helical band.

This criterion is a model result, not a device certification. Real devices have several transverse bands, electrostatic inhomogeneity, orbital magnetic effects, field-dependent parent pairing, interface disorder, smooth quantum dots, finite length, dissipative leads, and interactions. A large Zeeman field can satisfy the algebraic inequality while simultaneously destroying the induced gap.

RouteIntended mechanismWhat has been demonstratedWhat remains decisive
Semiconductor–superconductor nanowires and planar junctionsSpin–orbit coupling, Zeeman splitting, and induced pairingHard gaps, subgap states, end spectroscopy, nonlocal transport, parity-sensitive measurementsA robust topological gap tied to nonlocal end modes, controlled fusion, and braiding
Topological-insulator surface or nanowire hybridsSpin-momentum locking plus proximity pairingJosephson and Andreev phenomena, induced gaps, vortex and bound-state candidatesUnambiguous Majorana identification and scalable control
Magnetic atom chains on superconductorsHelical or ferromagnetic order, spin–orbit coupling, and proximity pairingEnd-localized zero-energy spectral weight and increasingly controlled atom-by-atom chainsBulk-gap and end-mode consistency, exclusion of trivial Shiba states, and manipulation
Iron-based superconductorsTopological surface bands combined with intrinsic pairingVortex-core zero-bias candidates and surface-state spectroscopyReproducible topological bulk/defect ledger and non-Abelian operations
Odd-parity and chiral candidate materialsIntrinsic unconventional pairingEvidence for unconventional order in several compoundsSettled pairing symmetry, a full protecting gap, and boundary response matching one invariant
Quantum-dot Kitaev chainsTunable crossed Andreev and elastic tunnelingTwo- and three-site spectra resembling short Kitaev chainsScaling to a long gapped chain and demonstrating topological, rather than sweet-spot, protection
Quantum simulatorsEncoded Kitaev or Majorana HamiltoniansControlled spectra and dynamics of target modelsThese validate a simulator, not an intrinsic electronic topological superconductor

The categories should not be collapsed. An engineered proximity device, an intrinsic bulk material, and a quantum simulator answer different scientific questions. A short quantum-dot chain can test the Kitaev Hamiltonian beautifully while remaining too small to establish asymptotic topological protection. A candidate intrinsic material can exhibit unconventional pairing while its topological invariant remains unknown.

No single local conductance curve measures a bulk invariant. The strongest case combines a bulk-gap diagnosis, a boundary or defect response, nonlocal consistency, parameter evolution through a transition, and controls that exclude known trivial mechanisms.

ObservationIdeal topological expectationImportant alternatives or caveats
Zero-bias end peakResonant Andreev reflection from an isolated Majorana can approach 2e2/h2e^2/h at zero temperatureSmooth-confinement Andreev states, quantum dots, disorder, heating, dissipation, and peak merging
Gap closing and reopeningRequired when a one-parameter path changes a bulk invariantA local probe can miss the bulk gap; ordinary level crossings and electrostatic rearrangements can mimic reopening
Correlated signals at both endsSeparated end Majoranas arise from one nonlocal topological segmentExtended trivial states and inhomogeneous multi-dot spectra can correlate the ends
Nonlocal conductanceCan reveal an extended bulk gap and distinguish local from nonlocal processesElastic cotunneling, crossed Andreev reflection, disorder, contacts, and charge redistribution require a full conductance matrix
Fractional Josephson responseFixed-parity branches can be 4π4\pi periodicParity relaxation restores 2π2\pi periodicity; Landau–Zener transitions and nonequilibrium occupation can mimic missing Shapiro steps
Vortex-core zero modeAn appropriate topological vortex binds an unpaired MajoranaOrdinary Caroli–de Gennes–Matricon states, impurity states, finite resolution, and vortex-line dispersion
Thermal boundary transportA chiral Majorana edge contributes a half-integer channelPhonons, additional edges, domains, contact thermalization, and non-topological neutral modes
Fusion or braid protocolOutcomes depend on fermion parity and braid order as predicted for Ising anyonsRequires initialized separated modes, a maintained gap, parity control, and controls against ordinary coherent dynamics

Fractional Josephson effect and 4π periodicity

Section titled “Fractional Josephson effect and 4π periodicity”

The fractional Josephson effect is also described as 4π Josephson periodicity, or as 4pi periodicity in plain-text literature searches. The notation refers to a parity-resolved spectral branch, not automatically to the equilibrium free energy.

For a topological Josephson junction with two coupled Majoranas, the fixed-parity energies can take the form

E±(ϕ)=±EMcos⁡ϕ2.E_\pm(\phi) = \pm E_M\cos\frac{\phi}{2}.

Following one parity branch through ϕ↦ϕ+2π\phi\mapsto\phi+2\pi changes its sign, so that branch returns only after 4π4\pi. In thermal equilibrium, however, the system can switch parity at the crossing and recover a 2π2\pi-periodic free energy. The phrase “4π4\pi Josephson effect” must therefore specify the drive frequency, parity lifetime, occupation dynamics, and alternative nonequilibrium mechanisms.

Current status and the topological gap protocol critique

Section titled “Current status and the topological gap protocol critique”

The experimental literature has advanced well beyond the first zero-bias peaks. Hybrid devices have been tested with a multistep topological-gap protocol, interferometric parity readout, three-terminal conductance, and short artificial Kitaev chains. Those are significant control achievements. They do not yet amount to a consensus demonstration of protected Majorana braiding.

The status is also actively contested. A 2023 study reported InAs–Al devices passing a specified topological gap protocol, and a 2025 experiment reported single-shot parity-sensitive interferometry in protocol-selected devices. A 2026 published critique argued that the underlying operating regions were substantially disordered and apparently gapless, challenging the topological interpretation. This disagreement is precisely why raw gap, locality, and calibration evidence must accompany the phase label.

The defensible hierarchy is:

  1. Ingredient demonstrated: proximity pairing, spin–orbit coupling, Zeeman control, a hard gap, or tunable Andreev processes.
  2. Candidate Majorana-compatible state: a zero mode with some expected stability or localization.
  3. Candidate topological phase: bulk and boundary evidence pass a declared protocol.
  4. Non-Abelian operation demonstrated: initialized modes exhibit fusion or braid-order statistics with trivial explanations excluded.
  5. Protected logical operation demonstrated: error suppression scales with separation, gap, or code distance in the intended architecture.

Claims should state the highest completed rung without borrowing the language of the next one.

Treating every BdG eigenvalue as an independent quasiparticle. The EE and −E-E branches are related by Nambu redundancy. Independent excitations may be counted using positive energies with zero modes handled separately.

Calling particle–hole conjugation ordinary charge symmetry. A BdG quasiparticle is generally a coherent electron–hole mixture, and electric charge is not the conserved quantity represented by C\mathcal C.

Equating a zero-energy state with a Majorana zero mode. Any tuned Andreev level can cross zero. A topological Majorana requires self-conjugacy plus a protecting bulk or defect invariant and suitable localization.

Applying the gapped tenfold-way table to a nodal spectrum. Nodal superconductors require invariants on loops or surfaces surrounding the nodes and may have different boundary phenomena.

Inferring topology from unconventional pairing alone. Odd parity or broken time reversal can be helpful, but the Fermi surfaces, full gap, dimensionality, and symmetry representation still matter.

Claiming quantized electric Hall conductance from a BdG Chern number. BdG quasiparticles do not carry a conserved U(1)U(1) charge, and the condensate contributes to electromagnetic response. A BdG Chern number fixes ideal chiral-Majorana thermal response, not a universal quantized electric Hall coefficient.

Calling an artificial few-site chain an intrinsic material realization. It may be a precise Hamiltonian emulator and an important device milestone without establishing a long-chain topological phase.

Calling a Majorana pair a protected qubit. Fixed total parity removes the two-state freedom of a single pair, and poisoning, finite overlap, control errors, and above-gap excitations remain.

For

H(k)=ξkτz+Δk τx,\mathcal H(k) = \xi_k\tau_z + \Delta k\,\tau_x,

where ξ−k=ξk\xi_{-k}=\xi_k and Δ\Delta is real, verify particle–hole conjugation with C=τxK\mathcal C=\tau_xK. If Φk\Phi_k has energy EE, construct its partner.

Solution

The Hamiltonian is real, so

CH(k)C−1=τxH(k)τx.\mathcal C\mathcal H(k)\mathcal C^{-1} = \tau_x\mathcal H(k)\tau_x.

Using τxτzτx=−τz\tau_x\tau_z\tau_x=-\tau_z and τx3=τx\tau_x^3=\tau_x gives

τxH(k)τx=−ξkτz+Δk τx.\tau_x\mathcal H(k)\tau_x = -\xi_k\tau_z + \Delta k\,\tau_x.

Meanwhile,

−H(−k)=−ξkτz+Δk τx,-\mathcal H(-k) = -\xi_k\tau_z + \Delta k\,\tau_x,

so the relation holds. If

H(k)Φk=EΦk,\mathcal H(k)\Phi_k = E\Phi_k,

then

τxΦk∗\tau_x\Phi_k^*

is an eigenvector of H(−k)\mathcal H(-k) with energy −E-E. It represents the creation operator of the same positive-energy quasiparticle, not a second independent excitation.

Use

M=sgn⁡[(μ+2t)(μ−2t)]\mathcal M = \operatorname{sgn} \left[ (\mu+2t)(\mu-2t) \right]

to classify μ=0\mu=0, μ=t\mu=t, μ=3t\mu=3t, and μ=−3t\mu=-3t for t>0t>0. Where must the bulk gap close?

Solution

For μ=0\mu=0,

M=sgn⁡(−4t2)=−1.\mathcal M = \operatorname{sgn}(-4t^2) = -1.

For μ=t\mu=t,

M=sgn⁡(−3t2)=−1.\mathcal M = \operatorname{sgn}(-3t^2) = -1.

Both are topological. For μ=±3t\mu=\pm3t, the product is 5t2>05t^2>0, so M=+1\mathcal M=+1 and the phase is trivial. The sign changes at

μ=±2t.\mu=\pm2t.

At those values the normal term vanishes at k=πk=\pi or k=0k=0, respectively, while the odd pairing term also vanishes. The BdG gap must therefore close.

Given {γa,γb}=2δab\{\gamma_a,\gamma_b\}=2\delta_{ab}, prove that

f=γ1+iγ22f = \frac{\gamma_1+i\gamma_2}{2}

obeys the canonical fermion algebra and derive iγ1γ2=2f†f−1i\gamma_1\gamma_2=2f^\dagger f-1.

Solution

Self-adjointness gives

f†=γ1−iγ22.f^\dagger = \frac{\gamma_1-i\gamma_2}{2}.

Using γ12=γ22=1\gamma_1^2=\gamma_2^2=1 and γ1γ2=−γ2γ1\gamma_1\gamma_2=-\gamma_2\gamma_1,

{f,f†}=1,f2=(f†)2=0.\{f,f^\dagger\} = 1, \qquad f^2=(f^\dagger)^2=0.

The number operator is

f†f=14(γ1−iγ2)(γ1+iγ2)=12(1+iγ1γ2).\begin{aligned} f^\dagger f &= \frac{1}{4} \left( \gamma_1-i\gamma_2 \right) \left( \gamma_1+i\gamma_2 \right) \\ &= \frac{1}{2} \left( 1+i\gamma_1\gamma_2 \right). \end{aligned}

Rearranging yields

iγ1γ2=2f†f−1.i\gamma_1\gamma_2 = 2f^\dagger f-1.

The bilinear measures the parity of the complex fermion, up to the stated convention.

Suppose

ε(L)=0.2 meV e−L/(0.25 μm)\varepsilon(L) = 0.2\,\mathrm{meV}\, e^{-L/(0.25\,\mu\mathrm m)}

after tuning to an oscillation maximum. How long must the wire be for ε<1 μeV\varepsilon<1\,\mu\mathrm{eV}?

Solution

Since 0.2 meV=200 μeV0.2\,\mathrm{meV}=200\,\mu\mathrm{eV}, the condition is

200e−L/(0.25 μm)<1.200e^{-L/(0.25\,\mu\mathrm m)} < 1.

Therefore

L>0.25 μm ln⁡200≃1.32 μm.L > 0.25\,\mu\mathrm m\, \ln 200 \simeq 1.32\,\mu\mathrm m.

This estimates only overlap splitting. It says nothing about poisoning, disorder-induced states, lead broadening, or whether a hard topological gap exists.

A two-dimensional class-D superconductor has C=3C=3. Find its chiral central charge and ideal low-temperature κxy/T\kappa_{xy}/T. Why is the coefficient not 3π2kB2/(3h)3\pi^2k_B^2/(3h)?

Solution

Each chiral Majorana mode has central charge 1/21/2, so

c−=C2=32.c_- = \frac{C}{2} = \frac{3}{2}.

Thus

κxyT=c−π2kB23h=π2kB22h.\frac{\kappa_{xy}}{T} = c_- \frac{\pi^2k_B^2}{3h} = \frac{\pi^2k_B^2}{2h}.

The larger proposed coefficient would count three chiral complex fermions. A complex channel contains two real Majorana channels, so it carries twice the thermal central charge of one Majorana mode.

At k=0k=0, the Rashba term vanishes. Show that the positive excitation energies are

E±(0)=∣VZ±μ2+Δind2∣E_\pm(0) = \left| V_Z \pm \sqrt{\mu^2+\Delta_{\mathrm{ind}}^2} \right|

and find the gap-closing condition.

Solution

At k=0k=0,

H(0)=−μτz+Δindτx+VZsx.\mathcal H(0) = -\mu\tau_z + \Delta_{\mathrm{ind}}\tau_x + V_Zs_x.

The spin operator sxs_x commutes with the Nambu matrices. In an sxs_x eigenstate with eigenvalue σ=±1\sigma=\pm1, the Nambu block has eigenvalues

E=σVZ±μ2+Δind2.E = \sigma V_Z \pm \sqrt{\mu^2+\Delta_{\mathrm{ind}}^2}.

Taking positive magnitudes gives the two displayed branches. The smaller branch reaches zero when

∣VZ∣=μ2+Δind2,\lvert V_Z\rvert = \sqrt{\mu^2+\Delta_{\mathrm{ind}}^2},

or

VZ2=μ2+Δind2.V_Z^2 = \mu^2+\Delta_{\mathrm{ind}}^2.

The derivation locates an ideal bulk transition; it does not prove that a finite disordered device has entered a robust topological regime.

A conventional superconductor has Δ=1.5 meV\Delta=1.5\,\mathrm{meV} and EF=150 meVE_F=150\,\mathrm{meV}. Estimate the Caroli–de Gennes–Matricon spacing. Could a spectrometer with 30 μeV30\,\mu\mathrm{eV} resolution distinguish the lowest ordinary level from zero?

Solution

The characteristic spacing is

ω0∼Δ2EF=(1.5 meV)2150 meV=0.015 meV=15 μeV.\omega_0 \sim \frac{\Delta^2}{E_F} = \frac{(1.5\,\mathrm{meV})^2} {150\,\mathrm{meV}} = 0.015\,\mathrm{meV} = 15\,\mu\mathrm{eV}.

The lowest half-shifted level is of order

ω02∼7.5 μeV.\frac{\omega_0}{2} \sim 7.5\,\mu\mathrm{eV}.

That is well below the stated resolution, so it can appear as a zero-bias feature. Spatial structure, temperature dependence, field evolution, higher levels, and a bulk topological diagnosis are needed before assigning a Majorana interpretation.

A finite nanowire shows a stable zero-bias peak at one end over a range of magnetic field. No second end contact is available, the induced gap softens continuously with field, and no gap reopening is resolved. The peak is labeled “proof of braidable Majoranas.” Evaluate the claim and propose four stronger checks.

Solution

The observation is compatible with a Majorana end mode, but it is far below proof of a topological phase and says nothing direct about braiding. A partially separated Andreev bound state, a smooth quantum dot, disorder, or a dissipative near-zero level can produce similar local conductance. The absence of a resolved reopened gap and of a second-end measurement removes two central consistency checks.

Stronger tests include:

  • map the bulk or nonlocal gap through a closing and reopening;
  • add a second end contact and test correlated end behavior while ruling out extended trivial states;
  • quantify peak height and temperature scaling against resonant-Andreev theory;
  • vary length or segment boundaries to test exponential end-mode overlap;
  • measure parity lifetimes and controlled fusion outcomes;
  • implement a braid-order protocol with initialized modes and ordinary coherent controls;
  • repeat across devices with a preregistered analysis and disorder diagnostics.

The justified label is “a robust local zero-bias feature compatible with several mechanisms,” not proof of a braidable excitation.

  • Unconventional Superconductivity owns pairing representation, nodes, relative signs, and time-reversal-breaking evidence; this page begins only when a BdG invariant, protecting gap, boundary or defect mode, and topological evidence are the claim.
  • BCS Theory develops Cooper pairing, coherence factors, thermodynamics, and conventional superconducting observables before topology is imposed.
  • Vortex Matter, Pinning, and Flux Flow owns entry, barriers, pinning, creep, collective order, and driven vortex motion; none of those material responses supplies the BdG invariant or Majorana evidence required here.
  • Pair-Density Waves and Exotic Orders distinguishes finite pair momentum and sign modulation from a protected Bogoliubov topological invariant.
  • Superconducting Proximity Effect owns the nontopological bulk induced correlations, inverse suppression, and spatial gap profile that must be established before a heterostructure is assigned a BdG invariant.
  • Proximity and Andreev Physics owns the ordinary NS-interface, BTK, few-mode Andreev-level, and hybrid-device baselines that candidate Majorana devices must exceed.
  • Symmetry Classification Preview owns the general tenfold-way vocabulary and explains the difference between physical time reversal and BdG particle–hole conjugation.
  • Topology in Quantum Matter gives the general bulk-gap, deformation, and evidence framework used here.
  • Topological Insulators is the charge-conserving comparison; proximity to its spin-momentum-locked boundary provides one route to effective topological pairing.
  • Edge and Surface States compares charged and neutral chiral boundaries, helical modes, Dirac cones, finite-size gaps, and the probes that can or cannot identify them.
  • Bulk–Boundary Correspondence explains stable boundary indices, removable mode pairs, regulator dependence, and the interacting limits of a free BdG count.
  • Symmetry-Protected Topological Phases explains stable interacting equivalence and the reductions Z→Z8\mathbb Z\to\mathbb Z_8 for class BDI and Z→Z16\mathbb Z\to\mathbb Z_{16} for class DIII.
  • Operator Identities collects the fermionic parity and Majorana-bilinear identities used in low-energy descriptions.
  • Superselection Sectors Preview explains why coherent control is restricted by total fermion parity.
  • Entanglement Spectrum develops a virtual-boundary diagnostic for free and interacting topological phases.
  • Topological Quantum Computation Bridge turns Majorana modes and parity measurements into an explicit encoding, gate, protection, and experimental-milestone ledger.
  • Topological Qubits expands that encoding into a complete tetron or anyon module with readout, reset, controller, outer code, universal-gate resources, and dated device evidence.
  • M. Sato and Y. Ando, “Topological Superconductors: A Review,” Reports on Progress in Physics 80, 076501 (2017), doi:10.1088/1361-6633/aa6ac7.
  • J. Alicea, “New Directions in the Pursuit of Majorana Fermions in Solid State Systems,” Reports on Progress in Physics 75, 076501 (2012), doi:10.1088/0034-4885/75/7/076501.
  • C. W. J. Beenakker, “Search for Majorana Fermions in Superconductors,” Annual Review of Condensed Matter Physics 4, 113–136 (2013), doi:10.1146/annurev-conmatphys-030212-184337.
  • B. A. Bernevig and T. L. Hughes, Topological Insulators and Topological Superconductors, Princeton University Press, 2013.
  1. N. Read and D. Green, “Paired States of Fermions in Two Dimensions with Breaking of Parity and Time-Reversal Symmetries and the Fractional Quantum Hall Effect,” Physical Review B 61, 10267–10297 (2000), doi:10.1103/PhysRevB.61.10267.
  2. A. Y. Kitaev, “Unpaired Majorana Fermions in Quantum Wires,” Physics-Uspekhi 44, 131–136 (2001), doi:10.1070/1063-7869/44/10S/S29.
  3. A. P. Schnyder, S. Ryu, A. Furusaki, and A. W. W. Ludwig, “Classification of Topological Insulators and Superconductors in Three Spatial Dimensions,” Physical Review B 78, 195125 (2008), doi:10.1103/PhysRevB.78.195125.
  4. A. Kitaev, “Periodic Table for Topological Insulators and Superconductors,” AIP Conference Proceedings 1134, 22–30 (2009), doi:10.1063/1.3149495.
  5. X.-L. Qi, T. L. Hughes, S. Raghu, and S.-C. Zhang, “Time-Reversal-Invariant Topological Superconductors and Superfluids in Two and Three Dimensions,” Physical Review Letters 102, 187001 (2009), doi:10.1103/PhysRevLett.102.187001.
  6. L. Fu and C. L. Kane, “Superconducting Proximity Effect and Majorana Fermions at the Surface of a Topological Insulator,” Physical Review Letters 100, 096407 (2008), doi:10.1103/PhysRevLett.100.096407.
  7. L. Fu and E. Berg, “Odd-Parity Topological Superconductors: Theory and Application to CuxBi2Se3\mathrm{Cu_xBi_2Se_3},” Physical Review Letters 105, 097001 (2010), doi:10.1103/PhysRevLett.105.097001.
  8. R. M. Lutchyn, J. D. Sau, and S. Das Sarma, “Majorana Fermions and a Topological Phase Transition in Semiconductor–Superconductor Heterostructures,” Physical Review Letters 105, 077001 (2010), doi:10.1103/PhysRevLett.105.077001.
  9. Y. Oreg, G. Refael, and F. von Oppen, “Helical Liquids and Majorana Bound States in Quantum Wires,” Physical Review Letters 105, 177002 (2010), doi:10.1103/PhysRevLett.105.177002.
  10. C. Caroli, P. G. de Gennes, and J. Matricon, “Bound Fermion States on a Vortex Line in a Type II Superconductor,” Physics Letters 9, 307–309 (1964), doi:10.1016/0031-9163(64)90375-0.
  11. K. T. Law, P. A. Lee, and T. K. Ng, “Majorana Fermion Induced Resonant Andreev Reflection,” Physical Review Letters 103, 237001 (2009), doi:10.1103/PhysRevLett.103.237001.
  12. V. Mourik et al., “Signatures of Majorana Fermions in Hybrid Superconductor–Semiconductor Nanowire Devices,” Science 336, 1003–1007 (2012), doi:10.1126/science.1222360.
  13. S. Nadj-Perge et al., “Observation of Majorana Fermions in Ferromagnetic Atomic Chains on a Superconductor,” Science 346, 602–607 (2014), doi:10.1126/science.1259327.
  14. D. Wang et al., “Evidence for Majorana Bound States in an Iron-Based Superconductor,” Science 362, 333–335 (2018), doi:10.1126/science.aao1797.
  15. M. Aghaee et al., “InAs–Al Hybrid Devices Passing the Topological Gap Protocol,” Physical Review B 107, 245423 (2023), corrected 2024, doi:10.1103/PhysRevB.107.245423.
  16. Microsoft Azure Quantum et al., “Interferometric Single-Shot Parity Measurement in InAs–Al Hybrid Devices,” Nature 638, 651–655 (2025), doi:10.1038/s41586-024-08445-2.
  17. H. F. Legg, “On the Robustness of Topological Gap Detection via Transport,” Nature 654, E22–E26 (2026), doi:10.1038/s41586-026-10567-8.
  18. S. L. D. ten Haaf et al., “Observation of Edge and Bulk States in a Three-Site Kitaev Chain,” Nature 641, 890–895 (2025), doi:10.1038/s41586-025-08892-5.
  19. J. Feng et al., “Long-Range Crossed Andreev Reflection in a Topological Insulator Nanowire Proximitized by a Superconductor,” Nature Physics 21, 708–715 (2025), doi:10.1038/s41567-025-02806-y.
  20. H. Jang et al., “Robustness of Majorana Modes to a Disorder Potential in Atomic Chains on a Superconducting Rashba Alloy,” Nature Physics (2026), doi:10.1038/s41567-026-03322-3.