Majorana Edge Modes in Topological Superconductors

Updated 12 April 2026

Majorana edge modes are emergent, self-conjugate quasiparticles localized at superconductor interfaces with a zero-energy signature and non-Abelian statistics.
Their existence is controlled by topological invariants like the Chern number and is highly sensitive to geometry, pairing symmetry, and external flux factors.
These modes yield distinct transport phenomena such as quantized charge and thermal conductance, underpinning potential applications in fault-tolerant quantum computing.

Majorana edge modes are emergent, non-Abelian, charge-neutral quasiparticle excitations localized at the boundaries or interfaces of topological superconductors. They are characterized by their equivalence to their own antiparticle (self-conjugation), protected zero-energy crossing, and robust nonlocal correlations, reflecting the underlying topological order of the parent bulk phase. The physics of Majorana edge modes is governed by both local and global features—such as pairing symmetry, geometric/topological properties of the host manifold, and external perturbations—yielding a diversity of transport, spectroscopic, and quantum-information signatures.

1. Theoretical Foundations: BdG Formalism and Topological Superconductivity

In gapful two-dimensional spinless $p_x+ip_y$ superconductors, the Bogoliubov–de Gennes (BdG) Hamiltonian captures the essential structure of Majorana edge and vortex modes: $H = \int d^2r\,\left[ \psi^\dagger\left(-\frac{\nabla^2}{2m}-\mu\right)\psi + \frac12\left(\Delta(p_x+ip_y)\psi^\dagger\psi^\dagger + \mathrm{h.c.}\right) \right]$ where $\psi(\mathbf{r})$ is the spinless fermion operator, $\mu$ the chemical potential, and $\Delta$ the $p$ -wave pairing amplitude (Quelle et al., 2016).

A more general surface with nontrivial geometry or topology is described by introducing local frames, the metric, and the spin connection. The pairing field $\phi$ transforms as a vector and satisfies the covariant constancy condition $D_\mu\phi=0$ in the ground state, with $D_\mu$ including both electromagnetic and geometric contributions. The BdG Hamiltonian on such surfaces includes these geometric effects explicitly: $\mathcal{H} = \begin{pmatrix} -\mu & \frac{1}{2\sqrt{g}}\{\sqrt{g}\,\phi,\partial_{-}\} \ -\frac{1}{2\sqrt{g}}\{\sqrt{g}\,\phi^*,\partial_{+}\} & \mu \end{pmatrix}$ Under appropriate boundary conditions, these Hamiltonians admit localized edge states whose existence and properties are dictated as much by global topology as by local structure (Quelle et al., 2016).

The appearance and number of Majorana edge modes are controlled by a bulk topological invariant, typically the Chern number: $H = \int d^2r\,\left[ \psi^\dagger\left(-\frac{\nabla^2}{2m}-\mu\right)\psi + \frac12\left(\Delta(p_x+ip_y)\psi^\dagger\psi^\dagger + \mathrm{h.c.}\right) \right]$ 0 where the Hamiltonian is recast as $H = \int d^2r\,\left[ \psi^\dagger\left(-\frac{\nabla^2}{2m}-\mu\right)\psi + \frac12\left(\Delta(p_x+ip_y)\psi^\dagger\psi^\dagger + \mathrm{h.c.}\right) \right]$ 1 with $H = \int d^2r\,\left[ \psi^\dagger\left(-\frac{\nabla^2}{2m}-\mu\right)\psi + \frac12\left(\Delta(p_x+ip_y)\psi^\dagger\psi^\dagger + \mathrm{h.c.}\right) \right]$ 2 for quasiparticle bands (Tsutsumi et al., 2010, Rachel et al., 2017). For chiral $H = \int d^2r\,\left[ \psi^\dagger\left(-\frac{\nabla^2}{2m}-\mu\right)\psi + \frac12\left(\Delta(p_x+ip_y)\psi^\dagger\psi^\dagger + \mathrm{h.c.}\right) \right]$ 3 pairing, $H = \int d^2r\,\left[ \psi^\dagger\left(-\frac{\nabla^2}{2m}-\mu\right)\psi + \frac12\left(\Delta(p_x+ip_y)\psi^\dagger\psi^\dagger + \mathrm{h.c.}\right) \right]$ 4, ensuring the presence of a single (per boundary) chiral Majorana edge mode.

2. Geometric and Topological Dependence: Annulus, Cylinder, Möbius, and Cone

The presence or absence of Majorana edge modes is sensitive to both the global geometry and topology of the host:

Annulus: With inner/outer edges and no intrinsic spin connection, the winding of the pairing field $H = \int d^2r\,\left[ \psi^\dagger\left(-\frac{\nabla^2}{2m}-\mu\right)\psi + \frac12\left(\Delta(p_x+ip_y)\psi^\dagger\psi^\dagger + \mathrm{h.c.}\right) \right]$ 5 is constrained by threading flux. In the absence of external flux, the necessary even pairing-winding condition for zero-mode existence fails ( $H = \int d^2r\,\left[ \psi^\dagger\left(-\frac{\nabla^2}{2m}-\mu\right)\psi + \frac12\left(\Delta(p_x+ip_y)\psi^\dagger\psi^\dagger + \mathrm{h.c.}\right) \right]$ 6), leading to the absence of edge Majorana zero modes in the ground state (Quelle et al., 2016).
Cylinder: The geometry allows $H = \int d^2r\,\left[ \psi^\dagger\left(-\frac{\nabla^2}{2m}-\mu\right)\psi + \frac12\left(\Delta(p_x+ip_y)\psi^\dagger\psi^\dagger + \mathrm{h.c.}\right) \right]$ 7 at zero flux, admitting an $H = \int d^2r\,\left[ \psi^\dagger\left(-\frac{\nabla^2}{2m}-\mu\right)\psi + \frac12\left(\Delta(p_x+ip_y)\psi^\dagger\psi^\dagger + \mathrm{h.c.}\right) \right]$ 8 (zero angular quasi-momentum) solution for the edge Majorana. Thus, each boundary component hosts a chiral Majorana mode at zero energy in the ground state.
Truncated Cone: Interpolates between annulus ( $H = \int d^2r\,\left[ \psi^\dagger\left(-\frac{\nabla^2}{2m}-\mu\right)\psi + \frac12\left(\Delta(p_x+ip_y)\psi^\dagger\psi^\dagger + \mathrm{h.c.}\right) \right]$ 9) and cylinder ( $\psi(\mathbf{r})$ 0). The holonomy of the spin connection $\psi(\mathbf{r})$ 1 shifts the winding, so the Majorana-supporting ground state is discontinuous across certain $\psi(\mathbf{r})$ 2. Adiabatic tuning of $\psi(\mathbf{r})$ 3 changes the current-carrying sector and illustrates the role of global holonomy.
Möbius Band: Being non-orientable, a globally consistent $\psi(\mathbf{r})$ 4 pairing is not possible. Instead, a domain wall (line defect) where the chirality of the pairing reverses is required; this domain wall itself binds a robust Majorana mode. The Möbius band, with only one physical boundary, supports exactly two Majoranas: one at the (sole) edge and one at the domain wall. The domain-wall Majorana's existence is protected by topology rather than boundary counting (Quelle et al., 2016).

These results illustrate that local flatness and gap protection do not fully determine Majorana edge mode configurations; the interplay of holonomy, orientability, and flux quantization controls the detailed pattern of edge and defect-bound Majoranas.

3. Microscopic Structure, Symmetry Classes, and Boundary Conditions

Majorana edge modes are realized as zero-energy, self-conjugate solutions to the BdG equations subject to boundary conditions ensuring self-adjointness and topological protection. The relevant symmetry class for most realizations is class D or BDI, supporting a $\psi(\mathbf{r})$ 5 or $\psi(\mathbf{r})$ 6 classification.

On nontrivial geometries, the quantization of angular (or linear) momentum is affected by the winding of the pairing field and the presence of fluxes or domain walls. For instance, on an annulus, the BdG ground state only allows edge zero modes for even winding; generically, the allowed values of angular quasi-momentum $\psi(\mathbf{r})$ 7 preclude zero energy solutions unless additional flux is threaded to modify $\psi(\mathbf{r})$ 8 to be even.

Self-adjoint boundary conditions for the spinless $\psi(\mathbf{r})$ 9-wave BdG operator take the form: $\mu$ 0 with $\mu$ 1 for parabolic-band parent superconductors, and the topological phase requiring $\mu$ 2 (Quelle et al., 2016). These conditions are essential for realizing protected, localized Majorana edge modes.

4. Transport Signatures and Experimental Detection

Majorana edge modes are associated with universal and quantized transport phenomena:

Quantized Charge Conductance: In 2D chiral topological superconductors (e.g. nano-scale “Shiba islands” of magnetic adatoms on superconductors), each chiral Majorana branch at the edge leads to a quantized zero-bias differential conductance $\mu$ 3, where $\mu$ 4 is the Chern number. This conductance doubling arises from perfect Andreev reflection and is a key spectroscopic observable (Rachel et al., 2017).
Chiral Supercurrent: Chiral Majorana modes also carry an equilibrium edge supercurrent whose circulation direction is set by $\mu$ 5. Its magnitude and sense can be mapped using scanning magnetometry, substantiating the topological nature of the edge (Rachel et al., 2017).
Thermal Conductance: Each chiral Majorana edge contributes a quantized thermal conductance $\mu$ 6, half that of a Dirac channel (Tsutsumi et al., 2010). Chirality inversion of Majorana edge modes (effected by supercurrent-induced Doppler boost) leads to a doubling of thermal conductance, a unique signature distinguishing spinful and spinless superconductors (Vela et al., 2021).
Fractional Josephson Effect: Majorana edge modes at boundaries of topological exciton condensates or at Josephson junctions with engineered barriers yield $\mu$ 7-periodic Josephson relations $\mu$ 8, identifiable via missing odd Shapiro steps and characteristic current-phase relations (Seradjeh, 2012, Subhadarshini et al., 24 Mar 2025).
Zero-Bias Peaks: Tunneling spectroscopy reveals robust, quantized zero-bias peaks at locations of Majorana edge modes (boundaries, domain walls), robust against weak disorder and system parameter variations. For systems with nonzero-energy edge modes (e.g., extended or trimerized Kitaev chains), finite-bias peaks at symmetric energies provide a means to distinguish topological from trivial edge states (Ghuneim et al., 4 Sep 2025, Zhang et al., 12 May 2025).

5. Generalizations and Topological Quantum Information

Majorana edge modes arise and are stabilized under a variety of generalizations:

Temporal Periodicity and Floquet Engineering: Periodically driven (Floquet) p-wave superconductors admit not only zero-energy (MZM) but also quasienergy $\mu$ 9 (MPM) Majorana edge modes, characterized by independent winding numbers, greatly enriching the phase diagram—allowing for multiple protected edge modes and “Floquet-enriched” nontrivial phases (Thakurathi et al., 2013, Wang et al., 2019, Bomantara, 2023, Zhou, 2019).
Number-Conserving and Interacting Models: Even in fully particle-number conserving systems with long-range interactions—without mean-field superconductivity—Majorana edge modes emerge as robust, nonlocal parity-carrying excitations, identified via ground-state degeneracy and spatially resolved edge correlation functions (Thomas-Markarian et al., 29 Aug 2025, Lisandrini et al., 2022, Zhang et al., 2016).
Quantum Information Applications: The nonlocality and braiding properties of Majorana edge modes support quantum memory encoding and the implementation of topologically fault-tolerant quantum gates. In hybrid architectures, combinations of Majorana edge modes and vortex-core Majoranas allow for high-fidelity Clifford operations and non-Abelian statistics, with the edge modes functioning as quantum memory while vortices carry out gates (Bedow et al., 13 May 2025).
Geometric Engineering and Domain Walls: The existence and robustness of edge Majoranas can be enforced or enhanced by introducing sharp geometric features (edges, corners) or topological defects (domain walls, Möbius band flips), yielding localized zero modes distinguished by their spatial profiles and response to flux insertion (Quelle et al., 2016, Ghuneim et al., 4 Sep 2025).

These properties render Majorana edge modes a central component in both the fundamental classification of topological phases and in technological applications, including quantum computation.

6. Robustness, General Classification, and Open Directions

Majorana edge modes are stabilized by a confluence of bulk gap protection, global geometric/topological constraints, and symmetry-enforced boundary conditions:

Protection by Topological Invariants: The bulk-boundary correspondence ensures that the total number of robust Majorana edge modes equals the change in topological invariant (Chern number or winding number) across the boundary. On non-orientable or holonomy-twisted manifolds, one must track the effective winding of the pairing around nontrivial cycles; even winding guarantees an unpaired zero-mode per boundary or domain wall (Quelle et al., 2016).
Stability under Perturbations and Disorder: Majorana edge modes remain spectrally isolated and pinned at zero energy under generic symmetry-preserving and even many types of symmetry-breaking perturbations, including onsite disorder, as long as the bulk gap remains open and particle–hole (or at least chiral) symmetry is preserved (Ghuneim et al., 4 Sep 2025, Zhang et al., 12 May 2025). In particle-number-conserving and interacting systems, the edge parity remains well-defined up to finite-size corrections (Thomas-Markarian et al., 29 Aug 2025, Lisandrini et al., 2022).
Distinguishing Topological from Trivial or Weak-Protection Regimes: Not all edge-localized states are robust Majorana modes; in weak topological or trivial phases, multiple unprotected near-zero modes may arise but can be readily split or hybridized by perturbations or dissipation (Zhang et al., 12 May 2025, Ghuneim et al., 4 Sep 2025). Mode classification and localization can be sharply determined by analytical solutions of the characteristic equations governing edge profiles, and by diagnostics such as the behavior of the energy splitting, spectral isolation, and spatial decay.
Impact of Geometry and Flux: The presence or absence of Majorana edge modes is nonuniversal and, in general, cannot be inferred from local properties alone; it is sensitive to the holonomy, flux quantization, and even nonorientability of the manifold, as demonstrated in the smooth but topologically distinct interpolation between cylinder, annulus, cone, and Möbius band (Quelle et al., 2016).

Majorana edge modes thus represent a central paradigm of boundary and defect-induced fractionalization in topological phases, unifying aspects of geometry, global topology, symmetry, and quantum information.

Key References: