Finite temperature effects on Majorana bound states in chiral $p$-wave superconductors

We study Majorana fermions bound to vortex cores in a chiral $p$-wave superconductor at temperatures non-negligible compared to the superconducting gap. Thermal occupation of Caroli de Gennes-Matricon states, below the full gap, causes the free energy difference between the two fermionic parity sectors to decay algebraically with increasing temperature. The power law acquires an additional factor of $T^{-1}$ for each bound state thermally excited. The zero-temperature result is exponentially recovered well below the minigap (lowest-lying CdGM level). Our results suggest that temperatures larger than the minigap may not be disastrous for topological quantum computation. We discuss the prospect of precision measurements of pinning forces on vortices as a readout scheme for Majorana qubits.

In a p-wave superconductor two vortex bound states of finite separation R hybridize and split in energy 34 .The splitting effectively produces a finite-energy two-level system of the Majorana qubit, with the levels distinguished by the fermion parity: two Majoranas can 'fuse' into a state with either zero or one fermions (two states of different parity), like two spin-1/2 particles can combine into either a spin-0 or spin-1 state.In some cases, a barrier to exploiting the vortex bound states in topological superconductors is the mixing with the tower of excited states in the vortex cores, known as the Caroli de Gennes-Matricon (CdGM) states [35][36][37] .The CdGM states are characterized by a level spacing, the 'minigap', of size δ ε ≈ ∆ 2 0 /E F , with ∆ 0 the full superconducting gap and E F the Fermi energy.Exciting these CdGM states does not lead to loss of quantum information, but it can hide the information so that it is very difficult to manipulate or to measure 38 .
Consider a qubit made from two vortices each having one Majorana zero mode.To measure the state of this qubit, one might, in principle, measure the force be-tween the vortices as we move them close together.At zero temperature there would be a difference in forces for the two different qubit states.At finite temperatures, excitation of the CdGM states reduces the force difference between the two different qubit states, although the quantum information is not lost until temperatures high enough that a bulk quasiparticle is excited that can carry away the fermionic parity 38 .It is thus generally assumed that temperature should be minimized as far as possible in realistic Majorana setups.
In this paper we examine the impact of thermal occupation of the CdGM states on the Majorana energy splitting in a two-vortex system.In a simple analytical model we find that the free energy difference between fermionic parity sectors lies exponentially close to the zero-temperature result for temperatures below the minigap, and decreases only as T −2 with increasing temperature T just above this temperature (and well below the superconducting gap).More generally, the contrast in free energy between parity sectors as two Majoranas are fused decreases with an additional factor of T −1 for each thermally-occupied bound state.The polynomial dependence suggests that temperature need not be an immediately limiting factor for topological quantum computation.This result relies, however, on having a 'sizable' minigap.If the vortex core is swamped with a continuum of in-gap states, we find the parity contrast to be exponentially suppressed in temperature.We discuss experimental aspects, including limiting time scales from quasiparticle poisoning and thermal vortex motion, towards the end of the paper.

II. BACKGROUND: MAJORANA BOUND STATES IN THE p + ip MODEL
We consider an effective spinless p x + ip y superconductor in two dimensions, described in the Bogoliubov-de Gennes (BdG) formalism in terms of a coupled eigensys-tem in the particle-hole basis 39 , where u n (v n ) is the particle (hole) component of the eigenstate, and (2) Here ∆(r) is the pairing function at position r = (x, y), k F and E F are the Fermi wavevector and Fermi energy, m is the effective electron mass, and ε n /2 is the energy of level n of the BdG spectrum.The BdG Hamiltonian in Eq. ( 2) features a particle-hole-symmetry which ensures that the resulting spectrum is symmetric around zero, such that ε n > 0 denote the energy differences between the particle and the hole states.Formally, the model belongs to class D in the classification of non-interacting topological superconductors and insulators 40,41 .
The quasiparticle annihilation (creation) operator γ( †) n associated with level ε n of the BdG Hamiltonian is a superposition of the spinless electron annihilation (creation) operators ψ( †) (r), The operators associated with localized Majorana (zero) modes obey the defining criterion γn = γ † n and the anticommutation relations {γ n , γm } = 2δ nm .
Vortices in the superconductor are regions around which the complex phase of the order parameter ∆ winds through 2π, with the gap magnitude going to zero in the vortex center.At position r, defining r = |r| and ϕ = arg (r), this is described by the order parameter ∆(r) = ∆ 0 e i ϕ f (r).Here ∆ 0 is the full (asymptotic) gap, the vorticity which is fixed to = 1 in the following, and f (r) is the radial profile, where f (r) ∼ r | | close to the vortex core.Assuming a profile f (r) = tanh (r/ξ), where ξ = v F /∆ 0 is the coherence length and v F is the Fermi velocity, the Majorana zero mode solution of Eq. (1) for a vortex in an infinite system with 2mE manifestly satisfying the Majorana condition u 0 = v * 0 .Here, N is a normalization constant, and J 1 is a Bessel function of the first kind.With N vortices centred at positions R j , the order parameter may be approximated by when assuming no spatial or phase fluctuations of the vortices, limiting the validity to the (zero temperature) mean field result for the gap profile 34 .Within the ground state manifold defined by a collection of Majorana zero modes, the pairwise exchanges of excitations constitute a higher-dimensional representation of the (non-Abelian) braid group generators, which is encoded by the unitary operators Ûn,n+1 = exp − π 4 γn γn+1 10,42 .The braiding operators do not span all unitary gates needed for universal quantum computation, but the remaining set of gates can be implemented in a non-topological way with arbitrarily small errors 11,43 .
With two vortices of finite separation R = |R 1 − R 2 | in a topological superconductor, each containing a Majorana mode, the energy splitting between the Majorana modes is given by 34 : This expression holds in the regime R ξ k −1 F , where k F is the Fermi wavevector.The hybridization energy above, ε 0 = ε + − ε − , is the energy difference between the two fermionic parity sectors, which, when the vortices are well-separated, are associated with the two-vortex wavefunctions , where Ψ n = u n (r), v n (r) (with denoting the transpose) is the wavefunction associated with vortex n.

III. EFFECT OF HEATING ON CDGM STATES
The lowest-energy levels, the CdGM levels, in the vortex core occur symmetrically about the Fermi energy, and their level spacing is given by the minigap, δ ε .The presence of such in-gap states, which derive from the finite density of states at the Fermi level in the normal state of the core, can lead to reduced distinguishability of the Majorana parity sectors 34,35 .In this section we consider the quantitative effect of occupying these states at finite temperature.More details on the CdGM levels can be found in Appendix A.

A. The partition function
In the superconducting condensate, where the number of particles is not conserved, we employ the grand canonical ensemble in evaluations of thermal expectation values.Since the (total) fermionic parity, i.e. the number of fermions modulo 2, is conserved the proper partition function is projected onto the even (+) and odd (−) par-ity sectors [44][45][46][47] : where Z 0 = n 2 cosh (βε n /2) is the partition function without parity restrictions, β = (k B T ) −1 , and all the products run over the positive energy levels in the BdG spectrum.We ignore states with energies far above k B T as they are suppressed by Boltzmann factors at low temperatures.The free energy in the ± parity sectors can be found in the usual way: B. Low-temperature model We consider two vortices at finite separation such that the Majorana modes in the vortex cores split according to Eq. ( 6).In the same fashion the excited CdGM states of the cores will also split (see Sec. III D).As a general and effective description at low temperatures, where only the lowest levels in the tower of CdGM states are thermally activated, we assume that the system has six levels: where δ ε is the minigap.Based on the splitting of the zero modes (Eq.( 6)) and numerics the deviation of the excited states from the minigap, w 1 and w 2 as defined above, are expected to decay with the vortex-vortex separation as w i ∼ exp(−R/ξ) for i = 1, 2.
The two parity sectors P = ±1 are treated separately, with the associated particle configurations being denoted by n ± .In the BdG-spectrum these configurations are depicted in Fig. 1.Notice that the number of particle excitations (disks above zero energy in the figure) is therefore even or odd in the respective parity channels.We consider the hierarchy of energy scales ε 0 |w 1 |, |w 2 | < δ ε , which we expect to apply when the vortices are further apart than about 3ξ (see Sec. III D).The partition function is given by The observable consequence of switching fermionic parity sector is that the free energy involved in bringing two vortices together changes.We term this difference in free energy between the two parity sectors the parity disparity, FIG. 1. Occupancy configurations, labeled by n± in columns, of the three-level system in the (a) even (P = +1) and (b) odd (P = −1) parity sector.The single-particle energy levels in the particle-hole symmetric BdG spectrum are denoted by ±εn/2.This depicts a low energy model for a system with finite distance between the two vortices as investigated numerically in Sec.III D.
which in the current case evaluates to: .
(11) As expected, the parity disparity decreases monotonically with increasing temperature.It reflects the feasibility of directly probing state of a Majorana qubit.In the limit of low temperature, k B T δ ε , the leading correction to the zero-temperature result is suppressed as ∼ exp(−βδ ε ): (12) At temperatures above the minigap we find F ± = − 1 β log 4 + O(βδ ε ) to leading order in the hightemperature expansion, simply reflecting thermal occupation of the four available configurations in each parity sector shown in Fig. 1.The linear corrections in βδ ε are independent of parity, such that the parity disparity acquires a leading order correction at second order in β, ) This inverse square law decay in the parity disparity with increasing temperature is weak enough that it will likely not preclude direct measurement at temperatures above the minigap.However, this toy model can only describe the regime k B T 2δ ε : at higher temperatures further CdGM states will be thermally excited.
Differences between the two parity channels are expected to be washed out at high temperatures.At low temperatures, the fact that the configurations 2 + and 3 + have excitation energies on the order of ε 0 greater than 2 − and 3 − (Fig. 1), respectively, distinguishes the two sectors.

C. Arbitrary numbers of core states
Now consider an arbitrary number of CdGM states thermally active in the two vortex cores.If the vortices are well-separated it is reasonable to assume that ε 0 < ε 1 ≈ ε 2 < ε 3 ≈ ε 4 < . . ., where ε 0 is the Majorana energy, and ε i>0 are the excited CdGM levels.The approximate equalities become exact as R/ξ → ∞.If the temperature is much greater than at least one of the levels, ε n k B T ∆ 0 for n 0, the parity disparity can be approximated by since tanh (βε m /2) ≈ βε m /2 1 ∀m n.Thus, the temperature decay above level ε n is algebraic: ∆F ∼ T −n .The above decay suggests that temperatures on the order of the minigap can be tolerated but that exciting additional bound states causes the parity disparity to decrease rapidly.The derivation holds in the limit of a well-defined and sizable minigap.If instead the vortex core contains a continuum of in-gap levels, the parity disparity of Eq. ( 14) is recast 16 and approximated 48 in the high-temperature limit as Here ρ 0 ≡ 1/δ ε is the density of states in the vortex core.The first term in the middle row of Eq. ( 15) appears when replacing the sum of equally spaced levels with an integral using the trapezoidal rule.The exponential suppression in k B T implies that vortex cores with densely packed ingap states, as on topological insulator surfaces 20 or in superconducting Pb monolayers 49 , would likely remain unworkable in the context of readout schemes and topological quantum computation.

D. Numerical results
Consider two vortices with separation R in a spinless p + ip superconductor.We first assume that T = 0 and that the two vortices have the same (positive) phase winding, as in Eq. ( 5).It should be noted that this is formally different from the vortex-antivortex and the antivortex-antivortex configuration in a p + ip superconductor.
We approach the problem numerically by solving the BdG equations with a finite element method for a large range of inter-vortex distances on a slab no smaller than [−10ξ, 10ξ] × [−9ξ, 9ξ] with Dirichlet boundary conditions.In Appendix B the calculation is repeated with one vortex replaced by a hole, with a flux quantum penetrating the hole.We fix the (dimensionless) parameters of the p + ip model to E F = 3, ∆ 0 = 1, and mk The CdGM wavefunctions are exponentially localized (with decay length ξ) around the two vortex cores, with spatial oscillations set by 2πk −1 F .Close to the vortex cores the expected small argument behaviour of the appropriate Bessel functions is recovered (e.g.Eq. ( 4)) 36 .In the well-separated limit the two vortex cores each host a single Majorana zero mode, and the CdGM levels become doubly degenerate as reflected in Fig. 2. Using the Majorana basis 39 for intermediate separation reveals that the wavefunctions have started to separate into two localized states at separations as small as about R/ξ ≈ 3.At larger separations the states Ψ ± , corresponding to opposite fermionic parity, are reasonable approximations to the true ground state.
The near alignment of the crossing of the energy levels seen in Fig. 2 is linked to the radial profiles of the vortex states which are given by Bessel functions of the first kind 36 , with an argument that increases slightly as a function of the excitation number.The Bessel-type wavefunctions enforce the energy splittings, which derive from the overlaps between the respective states 39 , to acquire similar oscillations.
We note that on a finite slab with Dirichlet boundary conditions, the net angular momentum induced by the two vortices enforces boundary states peaked around the edges of the sample, with an energy spacing set by the boundary length, v F /L.When the vortices are located far from the edge compared to ξ, these states are not affected by the presence of the vortices or their separation (see Appendix C).Considering instead periodic boundary conditions in both directions (solving the system on a torus) 1 and replacing one vortex by an antivortex 50 , removes the boundary states completely.
FIG. 3. The free energy F± (Eq.( 8)) normalized by the full gap ∆0 in the two parity sectors P = ±1 as a function of the inter-vortex separation R/ξ for a range of temperatures, with labels referring to the values of kBT /∆0.Here F± is the free energy averaged over R/ξ.
In Figs. 3 and 4 we display the influence of temperature in the presence of in-gap states for the two-vortex system, when Eq. ( 8) is applied to the numerically found energy levels of Fig. 2. As the temperature surpasses the minigap, the oscillations begin to smear.Fig. 4 displays the numerically exact temperature dependence of the parity disparity at R/ξ = 3.0 in the main figure.The black straight lines represent the simple power law derived below Eq. ( 14), when replacing tanh(βε n /2) by βε n /2 for k B T > ε n and 1 for k B T < ε n .Measuring the fusion channel of the corresponding qubit as defined by the two Majoranas, and hence the effect of braiding, thus becomes correspondingly difficult at temperatures well above the minigap.k B T /" 1 < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > 0.05 0.10 0.50 1 F/" 0 < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > 0.5 1.0 1.5 2.0 2.5 3.0 3.5 k B T/" 1 < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > 0.02 0.04 0.06 0.08 F/ 0 < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > 1.8 < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t >

3.0
< l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t >

4.2
< l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " ( n u l l ) " > ( n u l l ) < / l a t e x i t > Appendix B contains the same finite-temperature calculation with one vortex replaced by a hole.Qualitatively, the results are very similar to the vortex-vortex case.However, when one vortex is replaced by a hole, and reflection symmetry about (R 1 + R 2 )/2 in the order parameter is lost, the previous double degeneracy of the levels in the R/ξ → ∞ limit is lifted.The level spacing of the hole states is reduced with the hole circumference (cf.Ref. 51), generally making the issue of a small minigap worse when the hole radius is larger than ξ.We note also that the frequency of the CdGM level oscillations with the vortex-hole separation roughly doubles when comparing to the vortex-vortex case.

IV. MEASUREMENT CONSIDERATIONS
In practice, the topological splitting of the Majoranas constitutes a small contribution compared to the vortexvortex repulsion (see Appendix C for a listing of the expressions).The topological contribution will in principle remain clear under the Fourier transform of the free energy.The (Friedel-like) oscillations of frequency k F from the topological splitting manifest as a bump, with distinguishable features in the two parity channels, revealed for instance in Re{ F− − F+ }, with being the spatial Fourier transform of F ± .We show this in Fig. 5.By subtracting the free energy before and after a parity flip one is left with a finite result only if the vortices contain non-Abelian anyons.Although the k F bump in the energy difference would constitute strong evidence of successful braiding, the subtraction of the two curves will require a sensitivity set by ε 0 in the measured free energy.The measurement series of F − − F + can in principle be derived from the corresponding pinning force when one vortex is brought close to another, as measured before (F ± ) and after (F ∓ ) braiding with a third vortex.
In the following we consider the putative topological superconductor (Li 1−x Fe x )OHFeSe to evaluate the scope of applicability for force measurements to readout of Majorana qubits 28 .The material has a small Fermi energy E F 50 − 60 meV and a superconducting gap on the order of 2∆ 0 20 meV.Tunneling measurements on superconducting vortices identify a zero-bias peak, attributed to a Majorana zero mode, clearly separated from excited CdGM levels, with the unusually large minigap δ ε 1 meV ≈ 11 K.
Numerically, we find the parity disparity at temperatures below the minigap to be roughly ε 0 ≈ 0.04∆ 0 at R/ξ = 3.To measure the corresponding pinning force difference between the parity sectors would require a force sensitivity of δF ε 0 /ξ ≈ 0.05 pN, with ξ 1.4 nm 28 .This precision is at least an order of magnitude below previously reported force measurement thresholds [53][54][55] , although optimizing these measurements was not the primary goal of those studies.

A. Timescales and vortex motion
One of the leading threats to topological quantum computation comes from free electrons mixing with the Majorana mode, causing an uncontrolled qubit parity flip 11,38 .The timescale associated with this process is termed the 'poisoning time' 56 .The error rate Γ is set by the energy scale 57 Γ k B T e −δε/(k B T ) .
Using δ ε = 0.1 meV as a conservative estimate for the minigap of (Li 1−x Fe x )OHFeSe, and T = 100 mK, the expected poisoning time is t p /Γ ≈ 9 µs.Traversing a loop of radius 10ξ, say, would require a vortex speed of v 1 cm s −1 , which is in principle within reach of existing optical techniques 55 .Another source of qubit decoherence at finite temperature, potentially relevant for this setup, is phonon interactions 58 .
Thermal fluctuations change the vortex-vortex distance dynamically, causing smearing over an oscillating function (Eq.( 6)).This can greatly reduce the contrast between the parity sectors 34 unless the vortex pinning potential changes significantly on the order of the oscillation lengthscale 2πk −1 F 21 nm.This approach would therefore require the vortices to be artificially pinned in tight potentials 59 on this scale, if the temperature is on the order of the minigap.
Requiring adiabatic motion, to avoid the excitation of quasiparticles, introduces a lower time limit on the braiding operations.However, this timescale is small for the compound under consideration, t a /δ ε ≈ 0.7 ps.The braiding operations are therefore restricted to timescales t a t t p .This could be achievable in the near future if improvements continue to be made to individual vortex manipulation [53][54][55]59 for topological superconductors with sizable minigaps 27,28 .

V. CONCLUSIONS
In this paper we have demonstrated the effect, on the Majorana energy splitting of two bound states, of thermally exciting the CdGM states in vortex cores.The parity disparity (the difference in free energies between the two Majorana parity channels when the Majoranas are brought close together), which reflects the state of a Majorana qubit, lies exponentially close to the zerotemperature result for temperatures well below the minigap ∼ ∆ 2 0 /E F .Below the full gap, thermal excitation of higher CdGM states causes the amplitude to decay algebraically in temperature, with an additional factor of T −1 for each CdGM state thermally occupied on average.If the in-gap states are densely packed the suppression in temperature becomes exponential.
The relatively weak decay of the parity disparity with increasing temperature means that temperature does not necessarily need to be minimized in experimental scenarios.In fact, the Majorana modes can in principle coexist with excited states without loss of quantum information in any readout scheme based on the total fermion parity 38 .Local heating can potentially be used as another degree of freedom.This could prove useful, for example, if the magnitude of the external magnetic field is limited by other factors (working at a higher temperature can decrease the critical field H c1 ).Another possible application is to intentionally excite the CdGM modes in order to decouple the vortex motion from the state of the Majorana qubit, thus helping to increase decoherence times of the qubit.
As a promising candidate topological superconductor we draw specific attention to the intrinsic type-II superconductor (Li 1−x Fe x )OHFeSe, recently suggested as a possible platform for topological quantum computation 28 .With the results presented in the above sections, this material should have an experimentally accessible temperature range in which one can aim at probing the parity disparity.In the near-term, we hope the knowledge that finite temperature effects need not be disastrous for the measurement and control of Majoranas may open other avenues of enquiry.

2 F = 54 .
FIG. 2. The first energy levels εm/2 normalized by the full gap ∆0 of a two-vortex system.The CdGM levels as predicted with Eq. (A2) are shown in gray dotted lines.Colour maps of the core-localized wavefunctions |Ψ(r)| = (|un(r)| 2 + |vn(r)| 2 ) 1/2 for the three lowest-lying (positive energy) states are also shown.The wavefunctions are displayed for a vortex-vortex separation of R/ξ = 3.0 on a slab of size 8ξ × 6ξ.

FIG. 4 .
FIG. 4. The parity disparity ∆F = F− − F+ normalized by the ground state energy level ε0 a function of temperature normalized by the excited energy level ε1.In the main figure (log-log axes) we show the parity disparity for a vortex-vortex separation of R/ξ = 3.0.The CdGM levels are indicated with gray dotted lines, and the algebraic temperature law T −n , predicted in Eq. (14), is shown in black.The inset shows the parity disparity for three values of the inter-vortex separation, with values of R/ξ given by the labels.

FIG. 5 .
FIG. 5.The real part of the Fourier transform F± (expressions listed in Appendix C) of the parity disparity F− − F+ at zero temperature as a function of momentum k.Here, we set the coherence length to the value for Niobium 52 , ξ = 38 nm, and the Fermi momentum to kF = 0.1, 0.2 nm −1 in cyan and magenta, respectively.The labels refer to the prefactor of the topological contribution in Eq. (6), i.e. 4∆0π −3/2 .

FIG. C. 1 .
FIG. C.1.The Majorana energy oscillations as a function of the inter-vortex separation on top of the ground state energy shift (the vortex-vortex repulsion) from Eq. (C5) when including the lowest 150 states on a finite slab with Dirichlet boundary conditions.