1D Majorana Goldstinos and partial supersymmetry breaking in quantum wires

Realizing Majorana modes in topological superconductors, i.e., the condensed-matter counterpart of Majorana fermions in particle physics, may lead to a major advance in the field of topologically-protected quantum computation. Here, we introduce one-dimensional, counterpropagating, and dispersive Majorana modes as bulk excitations of a periodic chain of partially-overlapping, zero-dimensional Majorana modes in proximitized nanowires via periodically-modulated fields. This system realizes centrally-extended quantum-mechanical supersymmetry with spontaneous partial supersymmetry breaking. The massless Majorana modes are the Nambu-Goldstone fermions (Goldstinos) associated with the spontaneously broken supersymmetry. Their experimental fingerprint is a dip-to-peak transition in the zero-bias conductance, which is generally not expected for Majorana modes overlapping at a finite distance. Moreover, the Majorana modes can slide along the wire by applying a rotating magnetic field, realizing a"Majorana pump". This may suggest new braiding protocols and implementations of topological qubits.

In this work, we propose the realization of centrally extended quantum-mechanical SUSY, i.e., extended SUSY with central charges, in an experimentally-accessible condensed matter system, by employing a periodic array of partially-overlapping 0D quasi-Majorana modes to realize a Majorana chain with dispersive 1D Majorana fermions, and show how this can be achieved in proximitized semiconducting nanowires [75][76][77][78][79] via periodically-modulated magnetic fields [80][81][82][83][84][85] with large variations of the field intensity. We find that, for strong variations of the magnetic field amplitude, the topological mass gap M can assume alternatively positive and negative values along the wire, corresponding to topologically trivial and nontrivial segments, with partially-overlapping 0D quasi-Majorana modes localized at their boundaries and forming a 1D periodic lattice. This is the first realistic proposal to realize the Majorana chain model [32][33][34][35]86] in proximitized nanowires with strong spin-orbit coupling. The system exhibits a pair of dispersive and counterpropagating 1D Majorana modes delocalized along the whole wire and separated from the higher-energy bulk states, with a mass gap that can be tuned by externally applied fields. To characterize these emergent 1D Majorana modes, we introduce the concept of pseudohelicity and unveil the existence of an extended SUSY algebra with central charges. In the massless case, we found indeed that the 1D Majorana modes are pseudohelical, i.e., have opposite Majorana pseudospin, and exhibit centrally extended SUSY, with a finite zero-energy density of states and zero-bias peak delocalized along the whole wire. We find that the signature of the emergent SUSY can be revealed by the transition from a dip G = 0 to a quantized peak G = 2e 2 /h in the zero-bias conductance, or by inducing an adiabatic Majorana pumping in a sliding lattice of 0D quasi-Majorana modes, with quantized transport of one quasi-Majorana mode per a half cycle. Note that zero-bias peaks are generally not expected in the presence of several 0D quasi-Majorana modes localized at a finite distance. Moreover, we identify the massless Majorana fermion with a Goldstino, i.e., the Nambu-Goldstone fermion [87] associated with spontaneously broken SUSY from N = 4 to N = 2. Such a partial SUSY breaking of the extended SUSY is known to be possible only in the presence of central charges. This is perfectly compatible with our setup, in which we explicitly identify the central charges of the extended superalgebra. While the extended SUSY was proposed in several condensed matter systems . Lattice of 0D Majorana modes in a proximitized nanowire. A semiconducting nanowire with Rashba spin-orbit coupling α covered by a superconducting shell, with a periodically-modulated magnetic field induced by a regular array of nanomagnets and an externally applied field Ba. The periodic modulation of the field induces topologically nontrivial (NT) segments alternating to trivial ones, which correspond respectively to negative and positive values of the local Majorana mass M. 0-dimensional quasi-Majorana modes γAj and γBj localize at the boundaries between trivial and nontrivial segments, with overlaps w and v. 90], to our knowledge, this is the first real-world realization of centrally-extended quantum-mechanical SUSY, which plays essential roles in non-perturbative aspects of quantum field theory in high energy physics [91][92][93], and partial breaking of the extended SUSY algebras.

Effective model: 1D lattice of 0D Majorana modes
We consider a bipartite 1D lattice of 2N 0D Majorana modes where γ Aj , γ Bj are the Majorana operators corresponding to a single Dirac operator c j = (γ Bj + iγ Aj )/2 per unit cell, and with w, v ∈ R. This model is a special case of the Kitaev chain model [11] if µ = 2w and t = ∆ = −v. In momentum space we get up to a constant term, with H eff (k) = (0, v sin k, v cos k − w), and τ the vector of Pauli matrices. The energy dispersion is with a topological mass gap M eff = |w| − |v|. In the continuum limit (and assuming vw > 0), the Hamiltonian coincides with a 1D Dirac equation H = vk τ y + mv 2 − v 2 k 2 τ z with mass gap mv 2 = v − w and a quadratic correction in the momentum. The covariant form is obtained by multiplying the Hamiltonian by τ z . The Dirac equation is topologically trivial or nontrivial, respectively, for mv < 0 and mv > 0, i.e., for M eff > 0 and M eff < 0 [94].
In the massless case |v| = |w| the zero-energy eigenstates are doubly degenerate at gapless points and described by the fermionic operator d M = ( γ A + i γ B )/2 and its hermitian conjugate d † M , where the nonlocal Majorana operators are γ A = (1/ √ N ) j γ Aj , and γ B = (1/ √ N ) j γ Bj . The gapless state v = ±w separates two topologically inequivalent phases described by the topological invariant sgn M eff = ±1 where M eff = |w| − |v|. The Hamiltonian also exhibits 0D Majorana end modes in the nontrivial phase |v| > |w| in the case of open boundary conditions. The two end states are γ L ∝ j (w/v) j γ Aj and γ R ∝ j (w/v) N +1−j γ Bj localized at the opposite ends of the chain with localization length ξ eff = 1/| log |w/v||. The massless Majorana fields (|v| = |w|) describe a 1D free Majorana fermion in a 1+1D conformal field theory [7], which coincides with a pair of counterpropagating 1D Majorana modes. The solutions of the Dirac equation for a relativistic massless particle m = 0 form a helical pair of counterpropagating modes with opposite spin. However, our effective model is spinless: The role played by the spin degree of freedom is now played by the particle-hole degree of freedom. Hence, to characterize the properties of the 1D Majorana modes in our model, we introduce the Majorana pseudospin operator as τ /2 in analogy to the spin operator σ/2 (we use natural units). We obtain that the expectation values of the Majorana pseudospin have opposite directions for the two modes near the gapless point, being τ = sgn (v sin k)ŷ at k → 0, π for v = ±w, respectively. Analogously to the notion of helical modes [14], i.e., a pair of counterpropagating modes having opposite spin, we introduce the notion of pseudohelical modes, as a pair of counterpropagating modes having opposite Majorana pseudospin. Hence, the two modes form a pseudohelical pair near the gapless point. In this case, elastic backscattering is suppressed since the two modes crossing at zero energy are orthogonal, analogously to the case of helical modes [14]. The Majorana pseudospin introduced here generalizes the Majorana polarization [74,95]: The expectation values of the two components τ x,y /2 coincide with the Majorana polarization (up to prefactors).

Centrally extended superalgebra and partial supersymmetry breaking
In the massless case, we find that Eq. (1) exhibits extended N = 4 quantum mechanical SUSY given by the combined algebra defined by the supercharges which satisfy the superalgebra Here, H SUSY = H eff + 2h|v| (with h > 1) is the many-body Hamiltonian having nonnegative energy levels, P = N j=1 iγ Aj γ Bj the fermion parity, and T the translation defined by T γ Aj T † = γ Bj , T γ Bj T † = γ Aj+1(modN ) , which satisfies {T, P } = 0 and [T, H SUSY ] = 0 for |v| = |w|. (The supercharges Q 2 and Q 1 were introduced separately in previous works [32,33], but the existence of the combined N = 4 superalgebra has not been previously demonstrated in Majorana chain models, up to our knowledge.) All many-body eigenstates, including the groundstate, have superpartners with opposite parity. Thus, the Witten index is zero, and SUSY is spontaneously broken [27]. The two degenerate groundstates are the vacuum |0 and the state |1 = d † M |0 , respectively with even and odd fermion parity. The supercharge Q 1 |0 annihilates both groundstates However, the two groundstates are superpartners with respect to Q 2 : Since [T, H SUSY ] = 0 and {T, P } = 0, the eigenstates |0 and T |0 have same energy but opposite parities, which mandates T |0 = |1 and T |1 = |0 (up to a complex phase). This yields Hence, the supersymmetry Q 1 is unbroken whereas Q 2 is spontaneously broken: The N = 4 superalgebra (Q 1 , Q 2 ) is spontaneously broken down into the N = 2 superalgebra Q 1 . This mandates the presence of a Goldstino [87], which we identify with the massless Majorana fermion. The zero mass gap is protected by SUSY, i.e., the gap closes if and only if the Hamiltonian exhibits SUSY. We note that the supercharge Q 1 can be defined even in the presence of disorder and broken translational symmetry [33,96]: In this case, the corresponding Goldstino has a mass gap which closes if and only if the Hamiltonian exhibits SUSY. We note that the no-go theorem by Witten [27,97] forbids partial supersymmetry breaking in extended superalgebras with zero central charges: Either all supersymmetries Q i are broken, or they are all unbroken. However, the no-go theorem can be evaded in the presence of nonzero central charges [98], which allow the partial breaking of the extended SUSY.
Proximitized semiconducting nanowire with spin-orbit coupling and periodic magnetic field The Hamiltonian in Eq. (1) may describe the low-energy effective theory of a 1D topological superconductor with spatiallymodulated fields. Specifically, we consider a semiconducting nanowire with Rashba spin-orbit coupling and coated with a conventional superconductor, as in Fig. 1. A periodically-modulated magnetic field in the zx plane B nm (x) is induced by an array of nanomagnets [80][81][82][83]99] with magnetic moments parallel to the z-axis. In addition, we consider a uniform applied field B a in the z-direction, which can be used to control the lengths of the nontrivial and trivial segments. The wire is described by the Hamiltonian density where σ and τ are the vectors of Pauli matrices in spin and particle-hole space, m the effective mass, α the spin-orbit coupling, µ the chemical potential, b(x) = (gµ B /2)B(x) the Zeeman field in the zx plane (B(x) = B nm (x) + B aẑ is the total magnetic field), and ∆(x) the proximization-induced superconducting pairing. We require the wavelength λ of the periodically-modulated field to be comparable to the Majorana localization length, which is ξ M ≈ (b/E SO )α/∆ and α/∆ respectively for E SO = mα 2 /2 2 ∆ and ∆ (weak and strong spin-orbit coupling regimes) [9,100,101]. It is essential to our proposal to consider magnetic fields with large variations of the field intensity, contrarily to the well-known regime of periodically rotating fields with constant amplitude (or negligible amplitude variations), considered before [80,81,102,103]. Notice that periodically-modulated magnetic fields can also be induced by employing magnetic textures [104,105] or domain walls [106,107], while periodicallymodulated chemical potentials can be obtained by a periodic modulation of the width of the superconducting coating [96,108] in epitaxial 1D semiconductor-superconductor heterostructures [109][110][111][112].
If all fields are uniform, the Hamiltonian above reduces to the Oreg-Lutchyn minimal model [12,13]. The sign of the Majorana mass gap M = µ 2 + ∆ 2 − |b| characterize the trivial (M > 0) and nontrivial phases (M < 0) with M = 0 at the closing of the particle-hole gap. If the chemical potential, Zeeman field, or superconducting pairing is not uniform, one can define a local Majorana mass 1), where the coupling parameters v, w coincide with the Hamiltonian matrix elements between contiguous 0D quasi-Majorana modes separated by trivial and nontrivial segments, respectively. For |v| = |w|, the overlaps between contiguous 0D quasi-Majorana modes across trivial and nontrivial segments become equal, and the 1D Majorana modes become massless (gapless). We calculate the magnetic field B nm of the nanomagnets array via the finite-element method, see Supplementary Figure 1. We use the resulting Zeeman field in all numerical calculations. For reference, we find that a reasonable approximation of the Zeeman field is where b nm coincides with the magnitude of the average field along the wire, which gives |b nm (x)| ≈ |b nm | (2 − 2 cos (2πx/λ)).
We then discretize the Hamiltonian and take realistic values for the model parameters [44,71]. Moreover, we assume where ∆ 0 is the superconducting gap at zero field and B c the critical field, in order to account for the suppression of the superconducting pairing induced by the magnetic field. By numerically diagonalizing the discretized Hamiltonian, we obtain the energy E n and Nambu spinor Ψ n (x) of each eigenstate (see Methods section). Figures 2(a) and 2(b) show the local density of states (LDOS) at zero energy (x) = 1 π Im n |Ψ n (x)| 2 /(E n − iΓ) with finite broadening Γ = 0.1∆ 0 (to simulate the experimental conditions), and energy spectra as a function of the uniform applied magnetic field, in the case of periodic boundary conditions. The LDOS shows the presence of a periodic lattice of 0D quasi-Majorana modes with finite overlap, localized at the nodes of the Majorana mass gap M(x) = 0 [see also Fig. 2(e)]. This lattice corresponds to dispersive 1D Majorana modes with energy below the particle-hole gap and separated from the higherenergy bulk states, highlighted in Fig. 2(b), and with a Majorana mass gap equal to M eff = |w| − |v|. When the overlaps v, w between 0D quasi-Majorana modes across the nontrivial and trivial segments become equal when b a = b SUSY (which, in first approximation, occurs when L AB ≈ L BA ), the periodic lattice becomes invariant up to translations T [see Eq. (5)]. Hence, the wire exhibits SUSY and the dispersion becomes gapless (massless) with maximum LDOS at zero energy.
When the applied field increases above the threshold b NT , such that |b(x)| ≥ µ(x) 2 + ∆(x) 2 ∀x, the Majorana mass becomes negative on the whole wire and the trivial segments disappear. Conversely, when the applied field decreases below the threshold b T , such that |b(x)| ≤ µ(x) 2 + ∆(x) 2 ∀x, the Majorana mass becomes positive on the whole wire and the nontrivial segments disappear. In these two cases, the 0D quasi-Majorana modes at the ends of the trivial (or nontrivial) segments fuse into finite-energy Andreev-like fermionic modes. The continuous crossover between Majorana and Andreev-like modes is realized by increasing the overlaps between contiguous 0D quasi-Majorana modes at the ends of either the trivial or nontrivial segments, such that |v| |w| or |w| |v|, without closing the particle-hole gap. This crossover also occurs when the wavelength λ becomes smaller than the Majorana localization length λ ξ M : This results in larger overlaps v, w between contiguous 0D quasi-Majorana modes fusing into fermionic Andreev-like modes. In the opposite regime λ ξ M one has v, w → 0, which corresponds to decoupled quasi-Majorana modes with flat dispersion E k ≈ 0.
Figures 2(c) and 2(d) show the LDOS and the energy spectra in the case of open boundary conditions. For b a < b T and b a > b NT the local Majorana mass gap M(x) has the same sign along the wire, realizing a gapped trivial or nontrivial phase. In the latter case, 0D Majorana end modes γ L , γ R localize at the opposite ends. However, for b T < b a < b NT , the local Majorana mass gap M(x) changes its sign along the wire, with 0D quasi-Majorana modes described by the effective Hamiltonian in Eq. (1). This effective Hamiltonian can be trivial or nontrivial with Majorana mass gap M eff = |w| − |v|, as analyzed before. Hence, the nontrivial phase with 0D Majorana end modes γ L , γ R is also realized for b SUSY < b a < b NT (i.e., |v| > |w|). The LDOS integrated over the whole energy dispersion of the 1D Majorana mode is shown in Supplementary Figure 3. The LDOS as a function of energy below the bulk gap is shown in the Supplementary Movie.
If b nm (x) is approximately harmonic as in Eq. (9) one has |b(x) . Thus, assuming µ 2 + ∆ 2 ≈ b 2 nm + b 2 a the local Majorana mass gap becomes M(x) = (b nm b a / b 2 nm + b 2 a ) sin (2πx/λ + θ) which has equally-spaced nodes at x n /λ = θ/2π + n/2. Slowly varying the applied field direction θ induces the adiabatic sliding of the 1D lattice of 0D quasi-Majorana modes, corresponding to the pumping of one 0D quasi-Majorana mode every half-turn θ → θ+π and one full fermionic state every full turn θ → θ+2π of the applied field direction. In the case of periodic and closed boundary conditions, a half-turn of the applied field direction corresponds to the translation T of 0D quasi-Majorana modes entering the definition of the supercharges in Eq. (5). Figure 3(d) shows the intensity of the total magnetic field in Eq. (10) as a function of θ (see also Supplementary Figure 2). Figure 3(a) shows the evolution of the LDOS when the applied field direction θ turns around in the zx-plane. As the field rotates, 0D quasi-Majorana modes slide along the wire, resulting in an adiabatic pumping of 0D quasi-Majorana modes, as shown in Fig. 3(a). 0D quasi-Majorana modes translate by T at each half-turn θ → θ + π. Figure 3(b) shows the energy of the dispersive 1D Majorana modes below the gap, which corresponds to the sliding of the lattice of 0D quasi-Majorana modes. For reference, Fig. 3(c) shows the local Majorana mass M and its nodes M = 0 as a function of θ. This result may be generalized to the dynamical Floquet regime, with the possible realization of a finite "Majorana current" through the wire.

DISCUSSION
We theoretically proposed the realization of extended quantum-mechanical SUSY with central charges and dispersive 1D Majorana fermions in condensed matter, specifically, in a proximitized semiconducting nanowire via spatially-modulated magnetic fields. As shown in previous studies, a chain of Majorana modes exhibits both quantum mechanical SUSY [32,33] and space-time SUSY [34,35,86]. In our work, we unveiled the presence of an additional highly-nontrivial structure, i.e., an extended quantum-mechanical SUSY with central charges. This structure emerges as the combination of two coexisting N = 2 superalgebras: To our knowledge, the properties of these two coexisting superalgebras have not been previously demonstrated in the Majorana chain context. Extended SUSY with central charges has been one of the most important notions in quantum field theory and string theory over the decades since the second string revolution in the 90s and the revolution of quantum field theory by Seiberg and Witten [91][92][93]. In spite of its great importance in formal aspects of quantum field theory and string theory, all high energy theorists regard it as a useful tool which is not directly related to reality. This is because the extended N = 2 SUSY does not allow chiral fermions relevant for elementary particles such as quarks and leptons. In high-energy phenomenology, only N = 1 SUSY and its breaking are considered. There are many proposals to realize SUSY in condensed matter, e.g., in Bose-Fermi mixtures of ultracold atoms [113][114][115][116][117][118][119][120][121][122], Majorana Cooper-pair boxes [123], and at the boundaries of topological superconductors or insulators [30,124,125]. However, our model describes the first accessible example of extended SUSY with central charges and partial spontaneous SUSY breaking realized in nature. The action of the Goldstino is usually accompanied by higher derivative corrections determined thoroughly by a symmetrybreaking pattern, as found by Volkov and Akulov [87] (see also recent works on extended SUSY [126,127]). Our theory should be considered as the leading order of the derivative expansion. Thus, higher derivative correction terms, if one could obtain them, should be summed up to the Volkov-Akulov type action for 1D Majorana fermions.
While the N = 4 SUSY conformal field theories in 1D are known to be characterized by a central charge c = 2 [128], the central charge in our case is c = 1, due to the presence of two helical pair of counterpropagating Majorana fermions, implying only N = 2 SUSY instead of N = 4 SUSY. This is compatible with the fact that N = 4 SUSY is spontaneously broken and only N = 2 SUSY remains, giving strong evidence of the existence of the unbroken extended N = 2 SUSY in our model.
We notice that electronic interactions with the external environment may contribute to the pinning of the Majorana modes to zero energy in a typical nanowire setup. Without interactions, the SUSY regime corresponds to a single point of the parameter space, which coincides with the closing of the particle-hole gap, and occurs when the externally applied magnetic field is exactly b a = b SUSY . In the presence of electronic interactions, however, the lowest energy level may become pinned to zero energy, and therefore the particle-hole gap may remain closed in an extended window of the parameter space. This may occur due to the self-interaction of the charge distribution of the Majorana modes mediated by the external environment [129] in the case of a pair of Majorana modes localized at the edges of the nanowire. Perfectly spatially-separated Majorana modes are charge neutral. However, in our case, contiguous 0D quasi-Majorana modes have a finite overlap, which gives rise to a finite charge density [130,131]. This finite charge induces a screening charge distribution in the dielectric environment, which acts back onto the Majorana modes. This results in a self-interaction term that pushes the energy of the Majorana modes back to zero, as long as the nanowire has a larger dielectric constant of the external environment (e.g., the nanowire substrate) [129]. This mechanism can stabilize the SUSY by pinning the Majorana mass M eff ≡ 0 over an extended parameter space.
To experimentally realize our proposal, the wire must be much longer than the field periodicity, which must be comparable with the Majorana localization length, i.e., L λ ξ M ≈ α/∆. Moreover, variations of the gate and spin-orbit coupling fields must be negligible at length scales larger than λ, to guarantee an unbroken translational invariance at the mesoscopic level. Conversely, the physics described here is not affected by perturbations having a length scale shorter than λ (e.g., disorder). A different approach to realizing localized quasi-Majorana modes is by employing noncollinear magnetic textures or domain walls in complex magnet-superconductor heterostructures [106,107]. However, our proposal does not require the presence of a magnetic substrate and has the advantage of using proximitized semiconducting nanowires, which by far are the most extensively studied platform for topological superconductivity [75][76][77][78][79].
In a finite wire with 2N 0D Majorana modes with open boundary conditions and at zero temperature, the differential conductance exhibits 2N − 1 zeros and 2N quantized peaks G = 2e 2 /h, and the zero-bias conductance is zero [132]. In an infinite wire with identical couplings (v = w) instead, the conductance shows a zero-bias peak G = 2e 2 /h [132], which corresponds to the tunneling into the delocalized Majorana mode at zero energy. For v = w, the energy of the Majorana modes is lifted by the broken SUSY, and therefore there is no Majorana mode available at zero energy: In this case, the zero-bias conductance is zero. Hence, for a sufficiently long Majorana chain or equivalently in the case of closed boundary conditions (i.e., in a loop geometry), the transition between a zero-bias dip G = 0 to a peak G = 2e 2 /h signals the onset of SUSY at v = w. This transition should be observable in sufficiently long nanowires or in a setup with loop geometry (i.e., with closed boundary conditions) by varying the applied magnetic field close to the SUSY point b a = b SUSY . However, these signatures may be difficult to distinguish from conductance peaks induced by disorder, impurities, or finite-size effects. The characterization of the signatures of SUSY in the conductance will be the subject of future work. Stronger experimental signatures are the presence of spatially-periodic peaks in the LDOS along the whole wire, and their adiabatic evolution in the Majorana pump regime, obtained by varying the applied magnetic field direction. These signatures can be obtained by locally probing the differential conductivity in multiterminal setups [56,[133][134][135], or by scanning tunneling microscopy (STM) in epitaxial 1D heterostructures [109][110][111][112]. Other nonlocal fingerprints of SUSY can be provided by the signatures of the closing of the bulk particle-hole gap with a finite density of states on the whole wire. In particular, these signatures can be revealed via a tunneling probe placed at the bulk of the nanowire [56], by the quantized peak of the thermal conductance and electrical shot noise at the transition, and the doubling of the magnetoconductance oscillations in an Aharonov-Bohm ring geometry [136], or by the peak of the 4π component of the Josephson current in a superconducting ring geometry [137].
We also mention that the proposed experimental protocol to realize an adiabatic Majorana pumping in a sliding lattice of 0D quasi-Majorana modes may suggest new methods of braiding Majorana modes in 1D nanowire networks, which is the next milestone in the route to topological quantum computation [25, 26, 37-39, 42, 43].
Concluding, we proposed an experimentally accessible realization of a Majorana chain with emergent SUSY and with dispersive 1D Majorana fermions in proximitixed nanowires via spatially-modulated magnetic fields. In this system, we demonstrated the presence of an extended N = 4 SUSY with central charges and we identified the massless 1D Majorana fermions as the Nambu-Goldstone fermions (Goldstinos) associated with the spontaneous partial breaking of SUSY. Their experimental signatures are the finite LDOS at zero energy (zero-bias peak) delocalized on the whole length of the wire, and the dip-to-peak transition in the zero-bias conductance. This has to be contrasted with zero-bias peaks of 0D Majorana end modes, localized only at the ends of the wire, and to the general case of 0D quasi-Majorana modes, whose energy is lifted by their finite overlap. We finally showed how to realize an adiabatic Majorana pump by varying the applied magnetic field direction, which induces a sliding lattice of 0D quasi-Majorana modes with quantized transport of one quasi-Majorana mode per a half cycle. The manipulation of Majorana modes via spatially-modulated fields may lead to the realization of alternative non-abelian braiding protocols.

METHODS
The numerical results were obtained by discretizing Eq. (8) into a tight-binding model and diagonalizing the resulting Hamiltonian. The LDOS in the main article and in the Supplementary Information was calculated from the spectra as a function of energy. The momentum dispersion was obtained by Fourier-transform the tight-binding Hamiltonian and diagonalizing the resulting Hamiltonian in the momentum basis. In agreement with previous works [44,71], we consider an InSb nanowire proximitized by a conventional superconductor, and take m = 0.015 m e , α = 1 eV Å, b/B = 1.5 meV T −1 (g ≈ 50), B c = 3 T, and ∆ 0 = 1 meV at zero magnetic field. In this regime, the Majorana localization length is estimated to be ξ M ≈ 200-300 nm depending on the applied field b a . The magnetic field induced by the array of nanomagnets was calculated by numerical integration over a finite mesh. We considered a periodic array of 250 nm×250 nm×250 nm nanomagnets placed at a mutual distance 750 nm and at d = 250 nm from the wire, with remnant magnetic field B R ≈ 1 T and parallel magnetic moments in the z direction. See Supplementary Note 1 and Supplementary Figures 1 and 2 for more details.

DATA AVAILABILITY
The code used for the numerical simulations within this paper and the resulting data are available on Zenodo [138].