---
title: 'Helical Shiba Chain: Topological Superconductivity'
url: https://www.emergentmind.com/topics/helical-shiba-chain
type: topic
---

# Helical Shiba Chain: Topological Superconductivity

A helical Shiba chain is a one-dimensional array of magnetic adatoms deposited on a conventional \(s\)-wave superconductor in which the impurity moments form a spiral texture along the chain. Each adatom binds a Yu–Shiba–Rusinov (YSR) state, and the overlap of these bound states produces long-range, oscillatory subgap bands that can realize one-dimensional topological superconductivity and Majorana end modes. In the canonical planar geometry, with adatoms at \(x_n = n a\), the texture may be written as
\[
S_n = S[\cos(Qna)\,\hat x + \sin(Qna)\,\hat y],
\]
so that the helix itself supplies the spin-momentum locking needed for an effective \(p\)-wave channel, without requiring Rashba spin–orbit coupling in the host [1604.02040].

## 1. Geometry, spin texture, and defining mechanisms

The basic object is a chain of classical magnetic moments on a superconducting surface. Several equivalent parametrizations are used in the literature. A planar helix is often written as
\[
\mathbf S_n = S[\cos(Qna)\,\hat x + \sin(Qna)\,\hat y],
\]
while a more general conical texture is
\[
\mathbf S_n=\{\sin\theta\cos[\phi(n)],\sin\theta\sin[\phi(n)],\cos\theta\},\qquad \phi(n+1)-\phi(n)=g,
\]
with \(Q\equiv g\) in lattice units and \(\theta=\pi/2\) recovering the in-plane helix [1604.02040], [2507.21446]. Earlier tight-binding formulations also used a planar \(x\text{–}z\) helix,
\[
\mathbf B_n = B_0\{\sin(n\theta)\hat x+\cos(n\theta)\hat z\},
\]
which is equivalent up to spin-axis conventions [1308.6108].

In the helical setting, the magnetic texture mixes spin sectors spatially. The resulting YSR band behaves as a spinless \(p\)-wave-like superconductor once projected to low energies, which is why Majorana end states can appear. This mechanism is conceptually close to semiconductor-nanowire constructions, but the spin-momentum locking is generated by the helix itself rather than by an externally imposed Rashba term plus Zeeman field [1604.02040], [1308.3969]. In microscopic discussions of chain formation, the spiral is commonly associated with RKKY exchange, sometimes supplemented by Dzyaloshinskii–Moriya interactions or substrate-induced spin–orbit coupling [2412.11784].

The helical description is not the only route to the same low-energy structure. A ferromagnetic chain on a Rashba-coupled superconducting surface can generate an analogous topological YSR band because spin-flip correlations in the host Green’s function induce odd-parity pairing even for uniform magnetization [1407.6345]. This suggests that “helical Shiba chain” denotes both a literal spiral-moment chain and, in a broader effective sense, a class of impurity-band systems whose low-energy YSR sector acquires helical spin structure.

## 2. Microscopic formulation and effective Shiba-band theory

A standard microscopic starting point is the Bogoliubov–de Gennes Hamiltonian
\[
H = \left[\frac{p^2}{2m}-\mu\right]\tau_z + \Delta \tau_x - J\sum_n \mathbf S_n\cdot \boldsymbol\sigma\,\delta(\mathbf r-\mathbf R_n),
\]
with \(\mathbf R_n=(x_n,0,0)\), Nambu matrices \(\tau_i\), spin matrices \(\sigma_i\), exchange coupling \(J\), and \(s\)-wave gap \(\Delta\) [1604.02040]. For a single classical impurity, the YSR energy is
\[
E_{\mathrm{YSR}}=\Delta\frac{1-\alpha^2}{1+\alpha^2},\qquad \alpha=\pi \nu_0 J S,
\]
so the deep-YSR regime corresponds to \(\alpha\approx 1\), where the bound state approaches midgap and hybridization between impurities becomes strong [1604.02040], [1308.3969].

In the dilute-chain limit, projection onto the impurity YSR basis yields an effective one-dimensional BdG problem with long-range hopping and pairing,
\[
H_{\mathrm{eff}}=\sum_{m,n}\Psi_m^\dagger\,[h_{mn}\tau_z+\Delta_{mn}\tau_x]\,\Psi_n.
\]
For a three-dimensional host superconductor, the couplings decay as \(e^{-r_{mn}/\xi_E}/(k_F r_{mn})\) and oscillate at \(k_F\),
\[
h_{mn}\propto -\Delta \frac{\sin(k_F r_{mn})}{k_F r_{mn}}e^{-r_{mn}/\xi_E},\qquad
\Delta_{mn}\propto +\Delta \frac{\cos(k_F r_{mn})}{k_F r_{mn}}e^{-r_{mn}/\xi_E},
\]
with \(r_{mn}=|x_m-x_n|\) and \(\xi_E=\xi_0/\sqrt{1-E^2/\Delta^2}\) [1604.02040]. In two-dimensional host models the asymptotics become \(e^{-r/\xi}/\sqrt{k_F r}\), and the same host-mediated structure produces long-range odd-parity pairing after projection [1407.6345], [1505.05862].

The helical texture enters through spinor overlaps. In the planar case,
\[
\langle\uparrow|\uparrow\rangle_{mn}
=
\cos^2\frac{\theta}{2}\,e^{+ik_H r_{mn}}
+
\sin^2\frac{\theta}{2}\,e^{-ik_H r_{mn}},
\qquad
\langle\uparrow|\downarrow\rangle_{mn}
=
i\sin\theta\,\sin(k_H r_{mn}),
\]
so \(k_F\)-oscillations and helix-induced phases jointly determine the YSR bandwidth and the effective pairing symmetry [1604.02040]. In the deep-dilute regime this reduces to the long-range tight-binding model originally emphasized for helical chains, where hopping and pairing decay only algebraically over distances shorter than the coherence length [1308.3969], [1312.5723].

Beyond the deep-dilute approximation, the subgap problem can be formulated exactly as a nonlinear eigenvalue problem. One form used for helical chains is
\[
\tilde G^{-1}(E)\Psi=0,
\]
with a \(4N\times 4N\) matrix depending nonlinearly on
\[
\lambda = \frac{\Delta+E}{\sqrt{\Delta^2-E^2}}.
\]
This approach keeps the full energy dependence of the decay length and avoids truncation to short-range couplings [1604.02040], [1410.5639]. Complementary semi-analytical methods treat the infinite chain as a line impurity and extract the effective Green’s function by T-matrix techniques, thereby removing finite-size artifacts from direct tight-binding simulations [2102.02214].

## 3. Topological structure, invariants, and Majorana end states

For a generic helical chain, particle–hole symmetry is present and the system is in class \(D\). For a planar helix, an additional chiral symmetry appears, and the classification is upgraded to class \(BDI\) [1604.02040], [1308.6108]. In the planar case the literature frequently uses either a winding number or a \(Z_2\) parity sufficient for the single-Majorana phase. A convenient construction for the exact nonlinear formulation is the “topological Hamiltonian”
\[
\tilde H\equiv \tilde G^{-1}(0),
\]
from which one evaluates a real-space Pfaffian under periodic and antiperiodic boundary conditions; in clean periodic chains, topological transitions coincide with gap closings at \(ka=0,\pi\) [1604.02040]. For planar helices, one compact phase-boundary expression is
\[
\alpha_{0,\pi}
=
\frac{1}{\sqrt{(h_k^{\uparrow\uparrow})^2+(d_k^{\uparrow\uparrow})^2}}\Big|_{ka=0,\pi},
\]
which tracks the change of the \(Z_2\) invariant in the exact treatment [1604.02040].

In the original long-range Shiba-chain formulation, the Bloch Hamiltonian takes the form
\[
\mathcal H(k)=h_k\tau_z+\Delta_k\tau_x,
\]
with \(h_k\) and \(\Delta_k\) given by slowly decaying oscillatory sums. The planar helix then realizes a long-range analogue of a spinless \(p\)-wave chain. Unlike the nearest-neighbor Kitaev model, however, both hopping and pairing are power-law over distances \(r\lesssim \xi_0\), and for nonplanar helices the hopping amplitudes become complex, so \(h_k\neq h_{-k}\) and the pairing between \(k\) and \(-k\) is suppressed [1308.3969].

Long-range coupling also alters the critical theory. At the “Bragg point” \(k_F=k_h\), the transition between single-channel and two-channel regimes is unconventional: the critical point supports exponentially localized Majorana bound states with a short localization length unrelated to the topological gap, and away from that point the exponential core develops a power-law tail. In the notation of the exact solution,
\[
\xi_{\mathrm{eff}}=\frac{a}{\ln|\beta^{-1}|},
\]
while away from criticality the asymptotic tail behaves as \(\sin y/[y\ln^2 y]\) in the long-coherence-length limit [1312.5723]. This sharply distinguishes helical Shiba chains from short-range class-\(D\) wires.

Finite open chains in the topological phase host Majorana zero modes at their ends. In current-biased or FFLO generalizations, these modes remain diagnosed by bulk quantities such as the many-body polarization \(P_x\), with \(P_x=0.5\) marking the topological phase in the self-consistent FFLO studies [2412.11784]. A plausible implication is that “Majorana end mode” remains the unifying observable across static, current-driven, and Floquet helical-chain variants even though the bulk invariants and band constructions differ.

## 4. Disorder, vacancies, and anti-Shiba physics

The disorder problem in helical Shiba chains is atypical because missing adatoms do not merely broaden the YSR band; they generate their own localized subgap states. In a topological chain, a vacancy acts analogously to a magnetic impurity in a clean \(s\)-wave superconductor and binds a low-lying “anti-Shiba” state below the band edge of the regular chain [1604.02040]. Formally, the defect is treated by a T-matrix,
\[
T(E)=[I-VG_0(E)]^{-1}V,
\]
where \(V\) is the removal of the local exchange potential and poles of \(T(E)\) give the vacancy-bound-state energy \(E_{\mathrm{vac}}\) [1604.02040].

The physical interpretation is that a missing magnetic atom creates a weak link inside the topological YSR band, binding a localized pair of hybridized Majorana modes. A single vacancy therefore becomes a direct local probe of the nontrivial bulk phase. In the helical chain, vacancy-bound states are generically present throughout the nontrivial \(Z_2\) region, and their local density of states is concentrated near the defect and decays over \(\xi_E\) with Friedel oscillations at \(k_F\) [1604.02040].

At finite vacancy density, the hybridization of these defect states produces an anti-Shiba band. Its bandwidth scales parametrically as
\[
W_{\mathrm{anti}}
\sim
\Delta\times [\cos(k_F\ell)\ \text{or}\ \sin(k_F\ell)]\,
\frac{e^{-\ell/\xi_E}}{k_F\ell},
\]
up to texture-dependent projector factors, with \(\ell\approx a/p\) for vacancy concentration \(p\) [1604.02040]. The resulting deterioration of topology is not monotonic in the clean gap. Instead, dilute chains exhibit stripe-like regions of “unusual fragility” in the \((k_Fa,\alpha)\) plane whenever a single-vacancy level is tuned close to zero energy. In those regions, even a few percent vacancy concentration can split the nontrivial domain into disconnected islands [1604.02040].

This vacancy mechanism differs from random fluctuations of the Shiba coupling \(\alpha=\pi\nu_0 J S\). In both ferromagnetic and helical models, \(\alpha\)-disorder mainly erodes topology from phase boundaries, whereas vacancy disorder nucleates gapless phases in the middle of otherwise robust topological regions [1604.02040]. This suggests that defect spectroscopy is not only a diagnostic tool but also a stringent constraint on materials design.

## 5. Broader static variants: ferromagnetic, multichannel, antiferromagnetic, and multichain systems

The helical-chain mechanism belongs to a broader family of YSR-band topological superconductors. One important variant replaces the helical texture by a ferromagnetic chain on a superconducting surface with Rashba spin–orbit coupling. In that setting, odd-parity pairing arises from SOC-induced spin-flip propagators in the host rather than from noncollinear magnetic order, and the effective chain again maps to a spinless \(p\)-wave system with long-range couplings [1407.6345]. Momentum-resolved STM on atom-by-atom Mn chains on Nb(110) later provided evidence for multi-orbital Shiba bands, including a nondegenerate \(\alpha\)-derived band with a \(p\)-wave gap \(2\Delta_{\mathrm{ind}}\approx 360\,\mu\text{eV}\), while a \(\delta\)-derived band remained effectively gapless within the experimental resolution [2104.11497].

A second generalization keeps the impurity-band language but includes several angular-momentum scattering channels. In the multichannel YSR chain, simultaneous \(l=0\) and \(l=\pm1\) channels mixed by Rashba SOC generate a multiband Shiba structure. Depending on whether the deep band is \(s\)-like or \(p\)-like, the effective theory is either a single-band or two-band topological superconductor, and the inclusion of higher channels can enlarge the topological phase space relative to single-channel models [1505.05862]. This is directly relevant to realistic transition-metal adatoms, for which multiple YSR resonances are often observed experimentally.

Related static platforms can also be organized around altered magnetic order. In an antiferromagnetic chain on a superconductor, the \(\pi\)-pitch spin pattern may be regarded as a special helical limit. There, a supercurrent together with a weak Zeeman field generates a staggered spin-current term and drives the YSR band into a topological phase with spin-filtered Majorana edge states; the edge-spin polarization depends only on the parity of the number of magnetic moments [1402.5901]. In monolayer transition-metal dichalcogenide superconductors, a ferromagnetic adatom chain can realize the same effective spinless \(p\)-wave structure through spin–valley locking and parity-mixed intervalley pairing; the topological phase then depends not only on adatom spacing and moment direction but also on the orientation of the chain relative to the crystal axes [1604.02134].

The one-dimensional chain can also be coupled laterally into ladders. For planar helices, chirality-reversing domain walls bind two protected Majorana states, and multichain ladders exhibit a transverse-mode structure richer than a simple even–odd \(Z_2\) rule. In particular, a ladder of trivial chains can become topological because the transverse couplings renormalize the effective longitudinal parameters of each channel [1308.6108]. This multichain behavior is one of the clearest indications that helical Shiba systems are better understood as long-range, mode-resolved topological bands than as simple arrays of paired end Majoranas.

## 6. Supercurrent, FFLO, and Floquet extensions

The static helical Shiba chain has been generalized in two main nonequilibrium directions: phase-biased condensates and periodic drives. A uniform supercurrent in the host superconductor can be encoded by a phase gradient \(\Delta(\mathbf r)=|\Delta|e^{i\phi(\mathbf r)}\), with \(\phi(\mathbf r)=2\mathbf q\cdot\mathbf r\). After a gauge transformation, electrons and holes acquire opposite momentum shifts, and the effective Shiba-chain couplings gain terms linear in \(\varepsilon=v_F|\mathbf q|\cos\beta\), where \(\beta\) is the angle between current and chain direction [1406.4288]. For nonplanar helices this can tune a chain from gapless to topological gapped, or between trivial and topological gapped phases; for planar helices it mainly proliferates gapless regions [1406.4288].

A more elaborate current-biased extension is the self-consistent FFLO helical Shiba chain, where the order parameter takes the Fulde–Ferrell form
\[
\Delta(x)=\Delta e^{iqx}.
\]
With an out-of-plane Zeeman field, the spiral-induced effective SOC and the field-induced band asymmetry stabilize a finite-\(q\) superconducting state that supports end Majorana zero modes and a superconducting diode effect [2412.11784]. In that formulation the bulk topology is diagnosed by the polarization
\[
P_x=\frac{1}{2\pi}\operatorname{Im}\left[\operatorname{Tr}\{\ln(U^\dagger WU)\}\right],
\]
with \(P_x=0.5\) in the topological phase [2412.11784]. A field-free variant proximitizes the helical chain by a \(d\)-wave altermagnet, where induced altermagnetic spin splitting stabilizes topological FFLO superconductivity and strong nonreciprocal supercurrents without external fields [2507.21446].

Periodic driving introduces a distinct Floquet topology. In the driven helical Shiba chain, a sinusoidal modulation of the chemical potential,
\[
V(t)=V_0\cos(\Omega t),
\]
creates quasienergy gaps at both \(\varepsilon=0\) and \(\varepsilon=\pi/T\), enabling regular \(0\)-Majorana end modes and anomalous \(\pi\)-Majorana end modes [2304.02352]. Their topology is characterized by dynamical winding numbers computed from periodized evolution operators with twisted boundary conditions [2304.02352]. Transport theory for the same Floquet Shiba chain shows that the sideband-summed differential conductance obeys a Floquet sum rule, yielding \(n_M\times 2e^2/h\) for \(n_M\) end-localized Floquet Majorana modes at a given quasienergy in the end-resolved limit [2407.01135].

These driven and current-biased constructions do not replace the static helical chain; they enlarge its phase space. A plausible synthesis is that the chain has become a flexible platform for engineering effective \(p\)-wave superconductivity under conditions where the control parameter is no longer only \((k_Fa,\alpha,Q)\), but also superfluid momentum \(q\), drive frequency \(\Omega\), and nonequilibrium symmetry breaking.

## 7. Experimental signatures, materials, and current status

Scanning tunneling microscopy and spectroscopy remain the principal probes. In a conventional static topological phase, one expects end-localized zero-bias peaks separated from the bulk Shiba bands by a minigap, with spatial decay governed by the Majorana localization length. Vacancy spectroscopy adds a complementary bulk-sensitive probe: a single vacancy produces a pronounced LDOS peak at \(E_{\mathrm{vac}}\) below the clean band edge, and in fragility stripes this resonance moves close to zero bias [1604.02040]. Two vacancies produce bonding–antibonding splitting that decreases with separation, directly imaging the long-range hybridization of anti-Shiba states [1604.02040].

The clearest momentum-resolved bulk evidence so far comes from atom-by-atom Mn chains on Nb(110), where Bogoliubov quasiparticle interference revealed multi-orbital Shiba bands and identified one nondegenerate band with a gap shape and particle–hole asymmetry consistent with topological \(p\)-wave pairing [2104.11497]. In that system, the fitted \(\alpha\)-band parameters were
\[
\Delta_s=1.5\,\text{meV},\quad A=3.1,\quad B=2.35,\quad \xi=0.77\,\text{nm},\quad k_F=0.69\pi/a,\quad k_h=0.14\pi/a,
\]
and the inferred Majorana localization length from the gapped \(\alpha\) band was \(\xi_M\approx 3.8\,\text{nm}\) [2104.11497]. The same measurements also showed that a gapless \(\delta\) band can complicate the edge phenomenology, underscoring the multiband nature of realistic Shiba chains [2104.11497].

Materials proposals span conventional three-dimensional superconductors such as Nb and Pb with adatoms such as Fe, Co, or Mn, atomically assembled chains on heavy-element surfaces with strong Rashba effects, monolayer TMD superconductors, and heterostructures with \(d\)-wave altermagnets [2104.11497], [1604.02134], [2507.21446]. Across these platforms, the experimentally relevant tuning variables are impurity spacing \(a\), effective exchange \(\alpha\), helix pitch \(Q\), chain orientation, and defect density. The collected literature suggests a consistent strategy: maximize the clean topological gap by favorable interference of long-range hopping and pairing, while avoiding vacancy-tuned stripe regions and uncontrolled gapless bands.

Helical Shiba chains therefore occupy a distinctive position within topological-superconductivity research. They are impurity-band systems rather than proximitized semiconductor bands; they are intrinsically long-range rather than nearest-neighbor; and their topology is often most transparent in momentum-resolved bulk observables or defect-bound-state spectroscopy rather than in zero-bias end peaks alone.

Source: https://www.emergentmind.com/topics/helical-shiba-chain