Topological BdG Superconductors and Boundary Modes

A topological BdG superconductor is a gapped mean-field Hamiltonian whose symmetry class and bulk invariant obstruct deformation to an atomic limit while preserving the gap and symmetries. For a specified boundary or defect, the bulk invariant predicts protected particle–hole-related modes. This is a theorem about a model under stated hypotheses; realizing the Hamiltonian, maintaining the gap, and observing non-Abelian operations are separate experimental claims.

Required background. Inhomogeneous BdG theory supplies particle–hole redundancy, boundaries, and normalization.

Helpful background. The free-fermion periodic table and protected interfaces supply the classification framework.

The Kitaev-chain invariant

For spinless lattice fermions with real nearest-neighbor pairing,

$$H=\sum_j\left[-t(c_j^\dagger c_{j+1}+\mathrm{h.c.}) -\mu\left(c_j^\dagger c_j-\frac12\right) +\Delta(c_jc_{j+1}+\mathrm{h.c.})\right].$$

In the Nambu basis $(c_k,c_{-k}^\dagger)^T$,

$$H(k)=(-\mu-2t\cos k)\tau_z+2\Delta\sin k\,\tau_y.$$

The explicitly real model also has spinless time reversal $\mathcal T=K$ and the resulting chiral symmetry, so it lies in class BDI and has an integer winding number. If generic complex perturbations preserve particle–hole symmetry but break those additional symmetries, the classification reduces to class D. The parity of the BDI winding is the class-D $\mathbb Z_2$ invariant

$$\nu=\operatorname{sgn}[(-\mu-2t)(-\mu+2t)].$$

$\nu=-1$, equivalently $\lvert\mu\rvert<2\lvert t\rvert$ for $\Delta\ne0$, is the odd-winding topological phase. The sign formula is basis-dependent in appearance, but its change at a bulk gap closing is invariant. In the single-chain phase written here the BDI winding is $\pm1$, with its sign set by orientation conventions, while $\nu$ retains only its odd parity.

At $t=\Delta>0$ and $\mu=0$, introduce $c_j=(\gamma_{A,j}+i\gamma_{B,j})/2$ with $\gamma=\gamma^\dagger$ and $\{\gamma_m,\gamma_n\}=2\delta_{mn}$. The Hamiltonian pairs $\gamma_{B,j}$ with $\gamma_{A,j+1}$, leaving $\gamma_{A,1}$ and $\gamma_{B,L}$ absent. These are exact boundary zero modes. Away from the solvable point but within the phase, their wavefunctions decay over a localization length and their finite-chain splitting behaves schematically as

$$\delta E\sim e^{-L/\xi_M}\cos(k_FL+\phi).$$

Kitaev 2001, §§2–3 gives the construction and parity structure.

Bulk–boundary hypotheses

Protection requires all of the following:

• the relevant bulk or mobility gap remains open;
• perturbations preserve the symmetry class used by the invariant;
• the boundary separates phases with different invariants;
• distant Majoranas remain sufficiently separated and parity remains meaningful; and
• interactions do not invalidate the free-fermion classification or introduce a different low-energy topological order.

A zero-energy BdG vector obeying particle–hole self-conjugacy can be normalized as a Majorana operator, but numerical zero within tolerance is not proof of protection. Track spectral flow under every symmetry-allowed perturbation, system length, and boundary condition, and compute the bulk invariant independently.

In two-dimensional chiral $p_x+ip_y$ pairing, a Chern number controls chiral edge modes, while odd-vorticity defects can bind zero modes. Read and Green 2000, §§II–IV develops the weak- and strong-pairing distinction. Crystal symmetries can protect additional phases, but surface disorder that breaks the protecting symmetry changes the boundary claim.

From model to engineered platform

Spin–orbit-coupled semiconductor wires with Zeeman energy $V_Z$, induced pairing $\Delta_{\mathrm{ind}}$, and chemical potential $\mu$ provide one effective route. The clean single-band criterion

$$V_Z^2>\mu^2+\Delta_{\mathrm{ind}}^2$$

marks a model bulk transition. Multiband occupancy, orbital magnetic effects, parent-gap suppression, electrostatics, interactions, and disorder modify the inference. Lutchyn, Sau, and Das Sarma 2010 and Oreg, Refael, and von Oppen 2010 derive the ideal platform. Whether a device satisfies it belongs to the evidence page.

The validity map keeps model existence and platform identification on different sides of an evidence gate.

Bulk topology licenses protected modes for a specified model and interface. Device calibration, nonlocal spectroscopy, parity control, and adversarial trivial models are additional obligations. Original schematic, not to scale.

The paired-matter claim test matrix separates these levels.

Exercise

Find the phase boundaries. Determine where the Kitaev-chain spectrum closes.

Solution

$E(k)^2=(-\mu-2t\cos k)^2+4\Delta^2\sin^2k$. For $\Delta\ne0$, both terms vanish only at $k=0$ with $\mu=-2t$ or $k=\pi$ with $\mu=2t$. The gapped interval between them has opposite signs of the normal-state mass at the two particle–hole-invariant momenta and hence $\nu=-1$.

References

• Kitaev, A. Y. (2001). “Unpaired Majorana fermions in quantum wires.” Physics-Uspekhi 44(10S), 131–136. doi:10.1070/1063-7869/44/10S/S29.
• Lutchyn, R. M., Sau, J. D., and Das Sarma, S. (2010). “Majorana fermions and a topological phase transition in semiconductor–superconductor heterostructures.” Physical Review Letters 105, 077001. doi:10.1103/PhysRevLett.105.077001.
• Oreg, Y., Refael, G., and von Oppen, F. (2010). “Helical liquids and Majorana bound states in quantum wires.” Physical Review Letters 105, 177002. doi:10.1103/PhysRevLett.105.177002.
• Read, N., and Green, D. (2000). “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. doi:10.1103/PhysRevB.61.10267.