Multiple Andreev reflection and critical current in topological superconducting nanowire junctions

We study transport in a voltage biased superconductor-normal-superconductor (SNS) junction made of semiconducting nanowires with strong spin-orbit coupling, as it transitions into a topological superconducting phase for increasing Zeeman field. Despite the absence of a fractional steady-state ac Josephson current in the topological phase, the dissipative multiple Andreev reflection (MAR) current I_dc at different junction transparencies is particularly revealing. It exhibits unique features related to topology, such as the gap inversion, the formation of Majorana bound states, and fermion-parity conservation. Moreover, the critical current I_c, which remarkably does not vanish at the critical point where the system becomes gapless, provides direct evidence of the topological transition.


Introduction
Semiconducting nanowires (NWs) with a strong spin-orbit (SO) coupling in the proximity of s-wave superconductors and in the presence of an external Zeeman magnetic field B are a promising platform to study Majorana physics. Theory predicts that above a critical field B c ≡ µ 2 + ∆ 2 , defined in terms of the Fermi energy µ and the induced s-wave pairing ∆, the wire undergoes a topological transition into a phase hosting zero energy Majorana bound states (MBSs) at the ends of the wire [1,2]. Recent experiments have reported measurements of differential conductance dI/dV that support the existence of such MBSs at normal-superconductor (NS) junctions in InSb [3,4] and InAs [5] NWs. The main result of these experiments is an emergent zero-bias anomaly (ZBA) in dI/dV as B increases. In this context, the ZBA results from tunnelling into the MBS [6,7,8]. Although these experiments are partially consistent with the MBS interpretation [9,10,11,12], some important features such as the expected superconducting gap inversion were not observed. Moreover, other mechanisms that give rise to ZBAs, such as disorder [12,13,14,15], Kondo physics [16], or Andreev bound states (ABSs) [17,18], cannot be completely ruled out. In particular Zeeman-resolved ABSs in nanowires with charging effects [18] can give magnetic field dependencies essentially undistinguishable from some of the claimed Majorana experiments [5]. Furthermore, it has been recently pointed out that even ZBAs similar to the 0.7 anomaly in quantum point contacts may play a role in single barrier structures [19].
Stronger evidence could be provided by the observation of non-Abelian interference (braiding) [20], or by transport in phase-sensitive superconductor-normalsuperconductor (SNS) junctions. The latter approach, which typically involves the measurement of an anomalous "fractional" 4π-periodic ac Josephson effect [21,22,23], is much less demanding than performing braiding. Realistically, however, the fractional effect, detected through, e.g., the absence of odd steps in Shapiro experiments [23,24,25,26], may be difficult to measure (dissipation is expected to destroy it in the steady state), or may even develop without relation to topology [27]. Although it has been shown that the 4π periodicity survives in the dynamics, such as noise and transients [28,29,30], simpler experimental probes of MBSs are extremely desirable.
Here we propose the multiple Andreev reflection (MAR) current in voltage-biased SNS junctions made of NWs [31,32] as an alternative, remarkably powerful, yet simple tool to study the topological transition. This is made possible by the direct effect that gap inversion, MBS formation and fermion-parity conservation have on the MAR current I dc (V ) at various junction transparencies T N . For tunnel junctions, I dc (V ) traces the closing and reopening of the superconducting gap at B c , ∆ eff ∼ |B−B c |. This gap inversion can be shown to be a true topological transition by tuning the junction to perfect transparency T N = 1. In this regime, the limiting current I dc (V → 0) shows signatures of the parity conservation effects that are responsible for the fractional Josephson current in the presence of MBSs, but which, in contrast to the latter, survive in the steady state limit. Moreover, the detailed dependence of MAR as a function of T N has the fundamental advantage over NS junctions in that it contains information about the peculiar dependence of MBS hybridization with superconductor phase difference φ, despite not requiring any external control on it. Similarly, we show that another important phase-insensitive quantity, the critical current I c , remains unexpectedly finite for all B due to a significant continuum contribution, and exhibits an anomaly at the topological transition.
This paper is organized as follows. In section 2, we review the model for Rashba nanowires in the presence of both s-wave superconducting pairing and an external Zeeman field and describe how a spinless p-wave superconductor regime can be achieved. In particular, we discuss how the problem can be understood in terms of two independent p-wave superconductors, originated from the Rashba helical bands, and weakly coupled by an interband pairing term. This two-band description is very useful in order to understand the main results of this paper. Such results are discussed in section 3 which is divided in two parts. The first part (subsection 3.1) is devoted to the ABSs which are confined in the junction. The detailed evolution of these ABSs as the system undergoes a topological transition has not been discussed in the literature, to the best of our knowledge, and becomes essential in order to gain a deep understanding of transport across the junction, which is discussed in subsection 3.2. In this part we study the ac Josephson effect in nanowire SNS junctions, with focus on how the MAR currents reflect the topological transition in the nanowires as the Zeeman field increases. In particular, we present a thorough analysis of MAR transport in topological SNS junctions for arbitrary transparency of the normal part. Finally, we discuss in section 4 how the critical current I c does not vanish at the critical point where the system becomes gapless and, importantly, how I c provides direct evidence of the topological transition.

Rashba nanowire model and effective p-wave pairing
A single one-dimensional NW in the normal state is described by the Hamiltonian where m * is the effective mass, α so the SO coupling, µ the Fermi energy and σ i the spin Pauli matrices. An external magnetic field B along the wire produces a Zeeman splitting B = 1 2 gµ B B, where µ B is the Bohr magneton and g is the wire g-factor. The

Nambu Hamiltonian
models the NW in the presence of an induced s-wave superconducting pairing ∆ (here assumed real without loss of generality). The essential ingredient for a topological superconductor is an effective p-wave pairing acting on a single ("spinless") fermionic species [21]. SO coupling splits NW states into two subbands of opposite helicity at B = 0. At finite B, these two subbands, which we label + and − [black and orange lines in Fig. 1(a)], have spins canted away from the SO axis. The s-wave pairing ∆, expressed in the ± basis, takes the form of an intraband p-wave ∆ ++/−− p (p) = ±ip∆α so / B 2 + (α so p) 2 , plus an interband s-wave pairing ∆ +− s (p) = ∆B/ B 2 + (α so p) 2 [33]. Without the latter, the problem decouples into two independent p-wave superconductors, while ∆ +− s acts as a weak coupling between them. Each quasi-independent ± sector has a different (B-dependent) gap, which we call ∆ − (at small p) and ∆ + (large p), see Fig. 1(a,b). While ∆ + remains roughly constant with B (for strong SO coupling [9,11]), ∆ − vanishes linearly as B approaches the critical field, ∆ − ≈ |B − B c | ‡. This closing and reopening (gap inversion) signals a topological ‡ Note that, in general, ∆ − is at a small but finite momentum. However, as B approaches B c , ∆ − becomes centered at p = 0 and is approximately equal to |E 0 |, where E 0 is the zero momentum transition, induced by the effective removal of the − sector away from the low-energy problem. Below B c the NW is composed of two spinless p-wave superconductors, and is therefore topologically trivial. Above B c , ∆ − is no longer a p-wave gap, but rather a normal (Zeeman) spectral gap already present in the normal state, transforming the wire into a single-species p-wave superconductor with non-trivial topology. This phase contains MBSs, protected by the effective gap ∆ eff = Min(∆ + , ∆ − ), at the wire ends. Above a certain field B (2) c , the gap ∆ eff saturates at ∆ + and the physics of superconducting helical edge states in spin-Hall insulators is recovered [22,28].

Nanowire SNS junctions
In the previous section we described how a semiconducting nanowire with a strong SO coupling in the proximity of an s-wave superconductor and in the presence of an external Zeeman magnetic field B behaves as a topological superconductor above a critical field B c . Here we are concerned with the effects of this topological transition on the MAR current I dc (V ) across junctions formed with such nanowires. In particular, we consider SNS junctions of different normal transparencies T N . Experimentally, such geometry can be fabricated by partially covering a single NW with two superconducting leads and leaving an uncovered normal region in the middle. The coupling of the normal part of the NW to the superconducting leads can be tuned by local control of the electron density with a gated constriction. This can be realized by using, e.g., bottom-gates forming a quantum point contact, see Fig. 2. Such geometry has been successfully implemented experimentally in Ref. [19] for NS junctions, where control of the coupling between the superconducting and normal sections from near pinchoff (tunneling limit) to the multichannel regime is demonstrated. For simplicity, here energy of the lowest subband, E 0 = B − B c , and is related to the topological charge of the lowest superconducting band. [34] Out[429]=

HaL
HbL HcL HdL HeL H f L we focus on short SNS junctions § with single channel nanowires. For computation purposes, we consider a discretisation of the continuum model Eq. (1) for the Rashba nanowire into a tight-binding lattice with a small lattice spacing a. This transforms terms containing the momentum operator p into nearest-neighbour hopping matrices v.
are matrices in spin space, and t = 2 /2m * a 2 . The pairing is incorporated like in Eq. (1). A short SNS junction is modelled by suppressing the hopping matrix v 0 = νv between two sites in the middle of the wire, which represent the junction. The dimensionless factor ν ∈ [0, 1] controls the junction's normal transparency at B = 0, which we denote T N (ν). A phase difference φ across the junction is implemented by multiplying ∆ to the left and right of the junction by e ∓iφ/2 , respectively. Despite the § Results for the long junction limit will be published elsewhere [35].
simplicity of this description, it contains the relevant physics of a short SNS junction. As it has been shown for standard junctions [36], such physics essentially depend on the contact normal transmission as well as the voltage drop across it . Thus, we expect that a more detailed modeling, including e.g. a spatially-dependent voltage drop, would only modify the effective transmission T N (ν) which defines the different regimes we shall explain in the following.

Andreev bound states
In such short SNS junction, an ABS should form for each of the two p-wave sectors described in section 2 for B < B c , while only one, associated to ∆ + , should remain for B > B c . To support this picture, we present calculations of the local density of states (LDOS) at the junction in the transparent limit (T N = 1). This LDOS is peaked at the energy ± of the ABS, which is a function of the phase difference φ across the junction. For B = 0 ( Fig. 3a) the two ABSs are degenerate and confined within the gap ∆ ¶. As the Zeeman field increases, Fig. 3b, the two ABS split and the system develops the two distinct gaps ∆ + and ∆ − described in Section 2. Note that both ABSs are truly bound at energies below the lowest gap ∆ − , but only quasibound in the energy window ∆ − < < ∆ + . This is readily apparent in the plot as a broadening of the ABS resonances (see, e.g. Fig. 3c). As B approaches the critical field B c , ∆ − gets reduced, and becomes exactly zero at B = B c . Note that at this point the upper ABS has reached zero energy at φ = π and is quasibound for all energies, Fig. 3d. Upon entering the topological phase (B ≥ B c ), ∆ − reopens but one of the ABSs of the problem has disappeared (Fig. 3e). The surviving ABS associated to ∆ + arises due to the hybridization of the two emerging MBSs across the junction. Global fermionparity conservation protects the φ = π level crossing. Due to the residual ∆ +− s coupling between the two sectors, the ∆ + Andreev state is once more quasibound in the energy window ∆ − < < ∆ + . At high enough magnetic fields, ∆ + is the smallest gap of the problem and hence the Majorana ABS is truly bound, Fig. 3f. In long junctions, more ABSs can be confined in the junction [9,35,38]. These extra ABSs coexist with the ones described here and may, for example, anti-cross with the Majorana-like ∆ + Note that, for the sake of simplicity, we do not include the possibility of junctions containing resonant levels or quantum dots. A study of such junctions, including Coulomb blockade effects, is beyond the scope of this paper but might be useful in order to analyse the possibility of Majorana physics arising in experiments with short SNS junctions containing quantum dot nanowires, such as the ones reported in Ref. [4] ¶ Note that even this non-topological case is anomalous as the ABS energies do not reach zero at φ = π, unlike predicted by the standard theory for a transparent channel T N = 1 within the Andreev approximation µ ∆ [37]. We have checked that the energy minimum δ π does indeed vanish as µ/∆ grows, see Fig. C1 in Appendix C Andreev level, affecting its character near zero energy [9].

ac Josephson effect and MAR currents
Under a constant voltage bias V , the pairings ∆ to the left and right of the junction acquire an opposite and time-dependent phase difference, φ(t) = 2eV t/ . This induces Landau-Zener transitions between the ABSs and into the continuum, thereby developing a time dependent Josephson current with both I dc and I ac components. Such is the point of view in e.g. Refs. [39,29]. Alternatively, φ(t) can be gauged away into the hopping across the junction, v 0 (t) = νe −i eV tτz σσ c + rσ v σσ c lσ +h.c., where τ z is the z-Pauli matrix in Nambu space. By employing Keldysh-Floquet theory [36,40], we obtain the stationary-state time-dependent ac Josephson current I(t) = n e in eV t I n (note that only even harmonics survive, see Appendix A for full details). Here, we concentrate on the dc-current I dc = I 0 . The results for I dc (V ) at small, intermediate and full transparency are summarised in Fig. 4(a-c) for increasing values of B spanning the topological transition.
3.2.1. Tunneling regime. For non-topological tunnel junctions, dc-transport vanishes below an abrupt threshold voltage V t = 2∆ eff /e = 2∆ − /e (Fig. 4(a), blue curves). This well known result follows from the fact that there are no quasiparticle excitations in the decoupled wires for energy ∈ (−∆ eff , ∆ eff ) if B < B c . Indeed, to second order in perturbation theory in ν, the MAR current takes the form of a convolution between A 0 (ω) and A 0 (ω ± eV / ), where A 0 is the decoupled (ν = 0) spectral density at each side of the junction (Appendix B). [The trace of A 0 ( ), proportional to the local density of states (LDOS), is shown in Fig. 5(a)]. Hence, as B increases, the tunnelling current threshold follows the closing of the gap in the LDOS, until V t vanishes and I dc becomes linear in small V at B c (black curve). As B > B c , the gap reopens, but the threshold is now halved to V t = ∆ eff /e (red curves) + . The change, easily detectable as a halving of the slope of the threshold dV t /dB across B c , is due to the emergence of an intra-gap zero-energy MBS in the topological phase [see zero energy peak in Fig. 5(a)], which opens a tunnelling transport channel from or into the new zero energy state. Moreover, when B = B (2) c , ∆ − surpasses ∆ + , and ∆ eff saturates at ∆ + . This is directly visible in V t (B) as a kink at B 3.2.2. Intermediate transparency regime. As transparency increases, subharmonic MAR steps develop at voltages V t /n = 2∆ eff /en (n = 2, 3, 4, . . .), see Fig. 4(b). The specific profile of each step with V still contains information on the LDOS of the junction at energies around ∆ eff . At B = 0, the power-law LDOS for | | > ∆ results in a staircase-like curve I dc (V ) [blue line in Fig. 4(d)]. This shape is roughly preserved up to B = B c . For B > B c the MAR profile changes qualitatively, however. The subharmonic threshold voltages V t /n are halved (since V t is halved), and the MAR current profile becomes oscillatory instead of step-like. A blowup of the oscillations is presented in Fig. 4(d) (red curve), together with guidelines for the corresponding V t /n in gray. The emergence of oscillatory MAR steps, which here is connected to the formation of zero energy peaks in the LDOS owing to the localized MBSs, is well known in Josephson junctions containing a resonant level [42,43,44]. Note, however, that the oscillations in a topologically trivial system, such as for instance a quantum dot between two superconductors, arise at odd fractions of 2∆ eff , i.e. at voltages 2∆ eff /(2n − 1)e, instead of the ∆ eff /en of the Majorana case. Interestingly, this difference is ultimately due to the fact that a resonant level spatially localised within the junction cannot carry current directly into the reservoirs, while a zero energy MBS (essentially half a non local fermion) can. This same situation arises in d-wave Josephson junctions, which also exhibit oscillatory ∆ eff /en MAR subharmonics owing to the presence of mid gap states [45].
3.2.3. Transparent limit. In the limit T N → 1, ABS energies ± (φ) [Fig. 3(c)] touch the continuum at φ = 0. This has an important consequence. From the Landau-Zener point of view of the ac Josephson effect [39], the time dependence of φ(t) = 2eV t/ for an arbitrarily small V will induce the escape of any quasiparticle occupying an ABS into the continuum after a single φ(t) cycle. A given ABS becomes occupied with high probability in each cycle around φ = π if the rate dφ(t)/dt = 2eV exceeds its energy minimum (π) ≡ δ π . (Recall this energy is finite, since the Andreev approximation does not apply, see Appendix C.) One quasiparticle is then injected into the continuum per cycle, and a finite I dc (V δ π /e) arises. Below such voltage, however, the ABS remains empty, so that if δ π is finite, as is the case of a realistic non-topological junction [see Fig. 3(a-c)], one obtains I dc (V → 0) = 0 (valid for any transparency at B < B c ). This is in contrast to the conventional B = 0, T N = 1 result I(V → 0) = 4∆G 0 /e, predicted within the Andreev approximation (G 0 = e 2 /h).
After the topological transition, this picture changes dramatically. The two MBSs at each side of the junction hybridise for a given φ into a single ABS. This seemingly innocent change has a notable consequence. Since fermion parity in the superconducting wires is globally preserved, an anticrossing at φ = π, which would represent a mixing of a state with one and zero fermions in the lone ABS, is forbidden. Parity conservation therefore imposes δ π = 0 in the presence of MBSs, irrespective of T N or µ/∆ . This is a true topologically protected property of the junction, and gives rise to a finite I dc (V → 0) = 2∆ eff G 0 /e, i.e. half the value expected for the nontopological junction in the Andreev approximation. This abrupt change is shown in Fig. 4(c,e). The I dc (V → 0) MAR current in transparent junctions, therefore, directly probes the emergence of parity protection.

Critical current
In the transparent limit, a supercurrent peak [31,32,46] may hinder the experimental identification of the I dc (V → 0) limit, but itself holds valuable information about the transition. The critical current I c may be computed in general by maximizing the V = 0 (time-independent) current I(φ) respect to φ (including the contribution from the continuum). For a short transparent junction at B = 0, I c is maximum, and equal to I 0 c ≡ e∆/ in the Andreev approximation. Fig. 5(b) shows I c for increasing values of B. Naively, one may expect that a junction without a superconducting gap should not carry a finite supercurrent, but this is not the case here. At B = B c , I c is finite † †, while ∆ eff = 0 [the junction LDOS at criticality is also gapless, see Fig.  5(a)]. This gapless supercurrent comes from the + (φ) quasi-bound Andreev state in the continuum, which contributes almost as if it were a subgap ABS. It is thus a reasonable approximation to write I c as the sum of the critical current from each ABS. For B < B c , Note that residual splitting may survive in the topological phase for finite length nanowires, for which a finite (albeit exponentially small) coupling between four MBSs exist. † † Note also that in junctions with trivial superconductors, I c → 0 as the nanowire becomes helical at B = µ < B c [47]. Our result is therefore another nontrivial consequence of topology in the junction. I c ≈ 1 2 I 0 c (∆ + + ∆ − )/∆. The ∆ − contribution, however, should not be included for B > B c , leading to a discontinuity in ∂I c /∂B. This simple model gives a qualitative fit [dotted line in Fig. 5(b)] to the exact numerics (solid line), with deviations coming from corrections to the Andreev approximation, and contributions above ∆ + . Additional deviations in experiments, coming e.g. from the finite impedance of the electromagnetic environment, are not expected to alter the discontinuity in ∂I c /∂B, which remains a signature of the topological transition.

Conclusions
In conclusion, we have shown that the dc-current in voltage biased Josephson junctions is a flexible experimental probe into the various aspects of the topological superconducting transition in semiconducting nanowires.
Tuning the junction transparency one may obtain evidence of MBS formation as conclusive as a fractional Josephson effect, without requiring control of the junction phase. Moreover, we have found that the critical current in the wire does not vanish at the transition due to above-gap contributions, although its derivative with B exhibits a discontinuity as a result of the disappearance of one ABS. This behavior of I c provides a direct evidence of the topological transition. MAR spectroscopy and critical current measurements in nanowires similar to the ones studied here have already been reported [31,32].
Although we have focused here on the simplest case (single-band, short junction limit) we expect the main features of the topological transition to remain robust under more general conditions. Preliminary results in the quasi-one dimensional multiband case show that I c is a non-monotonic function for increasing magnetic fields. For weak interband SO mixing [48], the behavior discussed in Fig. 5(b) can be generalized and I c presents a series of minima at different fields corresponding to the topological transition of each subband.
Importantly, the alternative physical scenarios, such as, e. g., disorder [12,13,14,15] or Andreev bound states [18], that produce ZBAs in NS junctions (and thus mimic Majorana physics), cannot give the distinct features associated to global parity that were discussed here for SNS junctions. We therefore believe that experiments along the lines discussed in this paper could provide the first unambiguous report of a topological transition in nanowires, and the emergence of Majorana bound states.

Appendix A. Floquet-Keldysh formalism
Consider a mesoscopic system composed of two semi-infinite leads (labeled L and R), each in thermal equilibrium at the same temperature T and with the same chemical potential µ = 0. Each lead has a finite s-wave superconducting pairing ∆ α , where α = L, R. A central system (α = S), which may or may not be superconducting, is coupled to both leads through operator v. In its Nambu form, the Hamiltonian of the system readsĤ where the Nambu Hamiltonian matrix takes the general form Here h α is the normal Hamiltonian for each section of the system. The blocks delimited by lines denote the Nambu particle, hole and pairing sectors. If we apply a left-right voltage bias V through the junction, the Bardeen-Cooper-Schrieffer (BCS) pairing of the leads will become time dependent, ∆ L/R → e ±iV t ∆ L/R , while h L/R → h L/R ± V /2 (we take e = = 1). Both these changes can be gauged away from the leads and into the system by properly redefining c + i → c + i (t) = e ±iV t/2 c + i . This transformation is done also inside the system S, thereby effectively dividing it into two, the portion with an e iV t/2 phase (denoted S L ), and the portion with the opposite phase (denoted S R ). This restores H to its unbiased form, save for a new time dependence in h S → h S (V t), which is constrained to the coupling between the S L and S R , It is important to note that H(t) is periodic, with angular frequency ω 0 = V . In the steady state limit (at long times t after switching on the potential V ) all response functions and observables will exhibit the same time periodicity (all transient effects are assumed to be completely damped away). In particular, the steady state current I(t) = I(t + 2π/ω 0 ), so that for some harmonic amplitudes I n , in general complex, that satisfy I n = I * −n since I(t) is real. This current can be computed using the Keldysh Green's function formalism. [49] The standard expression for I(t) is computed starting from the definition of I(t) = ∂ t N L , where N L is the total number of fermions in the left lead. By using Heisenberg equation and the Keldysh-Dyson equation, one arrives at where The z-Pauli matrix τ z above acts on the Nambu particle-hole sector, The self energy from the left lead is defined as Σ a,< L (t , t) = v g a,< L (t , t)v + , where g L (t , t) = g L (t − t) stands for the left lead's propagator, when decoupled from the system (this propagator depends only on the time difference since the decoupled lead is time independent in this gauge). We define the Fourier transform of g as The retarded propagator in Fourier space is while the advanced g a L (ω) = [g r L (ω)] + . One can compute g < L (ω) = if (ω)A L (ω), where f (ω) = 1/(e ω/k B T + 1) is the Fermi distribution in the leads, and A L (ω) = i(g r L (ω) − g a L (ω)) is the Nambu spectral function. The g r L/R [and in particular A L/R (ω)] is assumed known, or at least easily obtainable from h L/R and v. Finally, the Green functions G r (t , t) and G < (t , t) correspond to the propagator for the full system, including the coupling to the leads. (Note that, in practice, since G is inside a trace in J(t), only matrix elements of G inside the S portion of the full system are needed). The retarded G r satisfies the equation of motion while G < (when projected onto the finite-dimensional system S) satisfies the Keldysh relation (t 2 , t). Since h S in H is time dependent, G propagators depend on two times; unlike Σ L/R or g L/R they are not Fourier diagonal. Instead, we can exploit the steady-state condition, which reads to expand the system's G as a Fourier transform in t − t and a Fourier series in t. We define The natural question is how the equation of motion is expressed in terms of the harmonics G n ( ). It takes the most convenient form if we redefine G n ( ) (where is unbounded) in terms of the quasienergy˜ ∈ [0, ω 0 ], i.e. =˜ + mω 0 This has the advantage that the equation of motion translates to a matrix equation analogous to that of a static system in Fourier space This is known as the Floquet description of the steady state dynamics in terms of sidebands, which appear formally as a new quantum number n. Time dependent portions of H(t) act as a coupling between different sidebands. The effective Hamiltonian for the n-th sideband is the static portion of H(t), shifted by −nω 0 . One therefore sometimes defines the Floquet "Hamiltonian" of the system as where, as before, h S n n = dte i(n −n)t h S (t). Likewise, one may define the Floquet self-energies as Σ Lnm (˜ ) = δ nm Σ L/R (˜ + nω 0 ), (since the leads are static, Σ is sideband-diagonal).
The Floquet equation of motion for G r nm (˜ ) can be solved like in the case of a static system. Within the S portion of the system, we have Boldface denotes the sideband structure implicit in all the above matrices. Similarly, the Keldysh relation takes the simple form (˜ ). Finally, the time averaged current I dc ≡ I 0 takes the form where the trace includes the sideband index. In a practical computation, the number of sidebands that must be employed is finite, and depends on the applied voltage bias V (the typical number scales as n max ∼ v 0 /V ). We employ an adaptive scheme that increases the number of sidebands recursively until convergence for each value of V .

Appendix B. Tunneling limit
It is possible to solve the I dc current explicitly in the tunnelling limit. To leading (second) order in the left-right coupling v 0 , Eq. (A.1) reduces, after some algebra, to where A α 11 is the particle-particle Nambu 2 × 2 matrix block of the spectral function of the α = L, R decoupled wire, ] * , and the trace is taken over spin space. The trace of A α (ω) is proportional to the local density of states. Fig. B1 shows results for the tunnel current using Eq. B.1 for increasing Zeeman fields. Overall, the agreement with the full numerics in Fig. 4(a) is very good and, importantly, all the relevant features such as, e. g., the closing of the gap, are captured by this tunneling approximation. regard the normal system as featureless, with constant Fermi velocity and density of states. In this case, a number of simplifications can be carried out in the computation of equilibrium transport properties. One important consequence of the approximation in the context of our work is that, in a short SNS junction with phase difference φ, and symmetric under time-reversal symmetry (B = 0 in our case), two degenerate Andreev states will appear of energy (φ) = |∆| 1 − T 2 N sin 2 (φ/2) [37]. This immediately implies that at perfect normal transparency T N = 1, the ABSs will reach zero energy at φ = π. In other words, in the Andreev approximation the ABS energy minimum in the non-topological phase will be (π) = δ π = 0. While in the topological phase, a zero δ π is a robust property, protected by parity conservation, a δ π = 0 in the non-topological phase is accidental, a direct consequence of the Andreev approximation, and is not protected by any symmetry. In fact, any deviation from the Andreev approximation will induce a finite splitting δ π . In semiconducting nanowires such as the ones considered in this work, this correction is very relevant. Indeed, for the nanowire to undergo a superconducting topological transition at a reasonable Zeeman field B, µ/∆ must remain relatively small (the  Figure C1. Minimum Andreev state energy δ π in a short SNS junction, at B = 0 and T N = 1. As µ/∆ grows, δ π decreases to zero, in agreement with predictions within the Andreev approximation. wire must be close to depletion, and far from the Andreev approximation), otherwise the critical field B c = µ 2 + ∆ 2 would be physically unreachable. The splitting δ π , therefore, remains a relevant quantity in the formation and detection of Majorana bound states.
The value of δ π in our system may be computed numerically. The most efficient way is to consider a short SNS junction with finite length superconductors, T N = 1, B = 0 and a phase difference φ = π. Since this system is closed, an exact diagonalization of the tight-binding Nambu Hamiltonian yields a minimum eigenvalue that is exactly δ π if the S leads are long enough (longer than the coherence length). We find that this quantity is finite in the case of wires close to depletion, µ ∆, and that it vanishes as one approaches the Andreev approximation regime µ ∆, see Fig.  C1. More specifically, we have found that δ π scales as δ π = c 1 ∆ 2 /(µ + ∆c 2 ) for some c 1,2 > 0, within very good precision.