Unveiling Signatures of Topological Phases in Open Kitaev Chains and Ladders

In this work, the general problem of the characterization of the topological phase of an open quantum system is addressed. In particular, we study the topological properties of Kitaev chains and ladders under the perturbing effect of a current flux injected into the system using an external normal lead and derived from it via a superconducting electrode. After discussing the topological phase diagram of the isolated systems, using a scattering technique within the Bogoliubov–de Gennes formulation, we analyze the differential conductance properties of these topological devices as a function of all relevant model parameters. The relevant problem of implementing local spectroscopic measurements to characterize topological systems is also addressed by studying the system electrical response as a function of the position and the distance of the normal electrode (tip). The results show how the signatures of topological order affect the electrical response of the analyzed systems, a subset of which being robust also against the effects of a moderate amount of disorder. The analysis of the internal modes of the nanodevices demonstrates that topological protection can be lost when quantum states of an initially isolated topological system are hybridized with those of the external reservoirs. The conclusions of this work could be useful in understanding the topological phases of nanowire-based mesoscopic devices.


Introduction
In the last decade, the properties of Majorana zero-energy modes (MZMs) hosted in topological superconductors have gathered considerable interest. Indeed, their non-abelian statistics has been proposed as the working principle of a fault-tolerant topological quantum computation [1]. However, a crucial challenge towards the topological quantum computer is to implement quantum operations of nearly degenerate quantum states by a dynamical process involving Majorana fermions; the main one is the braiding dynamics in superconducting nanowires. From the theoretical side, the simplest model for realizing MZMs is the one-dimensional spinless p-wave chain proposed by Kitaev [2]. A realistic implementation of the Kitaev model has been proposed lately in the pioneering work by Fu and Kane [3] who have predicted the presence of MZMs as a result of the proximity effect between an s-wave superconductor and the surface states of a strong topological insulator.
Many experimental works have been performed to realize the theoretical predictions of the Kitaev model and, in particular, Mourik and coworkers [4] have shown evidence of MZMs in the tunnel conductance of an InAs nanowire proximized by an s-wave superconductor. Two years later,

Majorana Fermions in the Kitaev Chain
In this section, we briefly discuss the main topological properties of the Kitaev chain. This discussion provides the appropriate starting point before treating the open quantum systems and serves to fix the notation. Following the Kitaev seminal work [2], we consider a system of spinless fermions confined to a one-dimensional lattice and subject to a p-wave superconducting coupling described by the Hamiltonian: where a j (a † j ) is the fermionic annihilation (creation) operator for a site j (j = 1, ..., L), t > 0 is the amplitude of nearest-neighbor hopping, ∆ > 0 is the amplitude of the superconducting pairing, while µ represents the chemical potential. For the Kitaev special case (t = ∆, µ = 0), using the Majorana basis c 2j−1 = a † j + a j and c 2j = −i(a j − a † j ), the Hamiltonian in Equation (1) becomes: where the hermitian Majorana operators fulfill the anticommutation relation {c k , c l } = 2δ kl . In this case, the Hamiltonian has two zero-energy Majorana modes, namely c 1 and c 2L , which are located at the end of the wire and combine to form a zero-energy non-local fermion f = (c 1 + ic 2L )/2. The other N − 1 fermionic modes (Bogoliubov modes)ã j = (c 2j+1 + ic 2j )/2, j ∈ [1, N − 1], are degenerate with energy 2t and connect the neighbouring sites in the bulk as depicted in Figure 1. For arbitrary values of ∆, t, µ, the Kitaev Hamiltonian H K can be diagonalized by the Bogoliubov transformation a j = ∑ m (u j,m α m + v * j,m α † m ), where the quasiparticle fermionic annihilation and creation operators α m and α † m obey canonical anticommutation relations. The Bogoliubov amplitudes u j,m and v j,m satisfy the Bogoliubov-de Gennes (BdG) equations: −µu j,m − t(u j+1,m + u j−1,m ) + ∆(v j−1,m − v j+1,m ) = E m u j,m , µv j,m + t(v j+1,m + v j−1,m ) + ∆(u j+1,m − u j−1,m ) = E m v j,m , with excitation energy E m ≥ 0. The diagonal form of H K reads: where E 0 is the energy of ground state |0 of Bogoliubov modes so that α m |0 = 0 for all α m . This model admits a topological phase for |µ| < 2t (see Figure 2) with a robust zero-energy edge mode, which is here labeled by setting m = M. The energy E M of this mode is not exactly zero, as expected in thermodynamic limit, and presents an exponential decay with the system size L (i.e., E M exp(−L/ξ), where ξ represents a characteristic decay length). The remaining non-topological modes of the spectrum are gapped and form a band.
The zero-energy Majorana modes are localised at the left/right edge of the wire and decay inside the bulk; the annihilation operators of such modes are given by: where the real valued eigenfunctions f L,j = (u j,M + v * j,M ) and f R,j = (u j,M − v * j,M ) are obtained by the explicit solution of the BdG Equation (2). Following Ref. [12], it is possible to show that f L,j = 2|A|ρ j sin(jθ) and f R, In the thermodynamic limit (L → ∞), the presence of such edge modes results in an exact degeneration of the ground state that corresponds to the presence or absence of the non-local fermion α M = (γ L + iγ R )/2. In fact, in this limit, the two ground states (zero-energy states) |+ ,|− have different parity and thus α M |− = 0 and α † M |− = |+ .

Ladder of Two Kitaev Chains
Various generalizations of the single Kitaev chain have recently appeared [13][14][15][16][17][18], the simplest one being obtained by coupling two single Kitaev chains with transversal hopping and pairing terms to form a Kitaev ladder (KL) [19]. The Hamiltonian of the system is given by: where H K 1 and H K 2 are the Hamiltonians of the isolated Kitaev chains given in Equation (1), while H K 1 ,K 2 describes the interchain coupling whose expression is given by: The labels 1, 2 in Equation (5) denote the two chains, j is the site index, t 1 represents the transversal hopping amplitude and ∆ 1 is the transversal pairing term. All the model parameters, namely t, t 1 , ∆ and ∆ 1 , are taken as real numbers. In the momentum representation, the Hamiltonian (4) can be written as where represents the Nambu spinor. The Hamiltonian is given by with k = −2t cos k − µ and ∆ k = 2∆ sin k. By construction, the Hamiltonian (8) satisfies the particle-hole symmetry (Ξ) implemented by the operator which is defined in terms of the complex conjugation operator K and of the Pauli matrix σ x , while I represents the identity operator. The Hamiltonian satisfies also the time-reversal symmetry (τ = K) and the chiral symmetry According to the mentioned symmetry properties enjoyed by H(k), the Hamiltonian model belongs to the BDI symmetry class with Z topological index [20]. The spectrum of the Hamiltonian can be easily obtained and is given by The topological phase diagram can be obtained by the calculation of the winding number [18,21]. Using the chiral basis, the Hamiltonian (8) assumes the following off-diagonal form [21] where U is the basis change matrix and A k is the 2 × 2 matrix Following Zhou et al. [18], we calculate the winding number as and compute the topological phase diagram, which is systematically analyzed in the following discussion. To provide evidence of the richness of the topological phases of the model, in the following, we present our results using an extended parameters range, which in principle could be completely explored in cold atoms experiments rather than in condensed matter systems. In Figure 3, we present the topological phase diagram in the (t 1 , µ) plane by fixing ∆ = 0.8 and setting different values of ∆ 1 (∆ 1 = 0, 0.09, 0.5, 0.8, respectively), while taking t as an energy unit. The integers in the figures correspond to the winding number and are related to the number of MZMs at every edge of the ladder. By direct inspection of Figure 3, we notice that an increase of the pairing term ∆ 1 produces an enlargement of the trivial region, labeled by 0. The expansion of the trivial region changes the phase boundaries and induces a modification of the linear phase boundaries, which are peculiar to the ∆ 1 = 0 case. A further consequence of the enlargement of trivial region is the possibility, changing the chemical potential and setting appropriate t 1 values, to have a direct transition from the two MZMs region to the trivial phase without crossing the one MZM phase. Interestingly, an additional consequence of the interchain coupling is that the critical value of the chemical potential, µ c , defining 2-0 phase boundary is progressively lowered as ∆ 1 increases. As a consequence, µ c becomes smaller than the one established by Kitaev for the 1-0 phase boundary of a single chain, i.e., µ c < 2t.
In Figure 4, we present the topological phase diagrams in the (∆ 1 , µ) plane. The different panels are realized by setting two values of t 1 , namely t 1 = 0.6 (left panel) and t 1 = 2.1 (right panel), while retaining the remaining parameters as in Figure 3. Direct observation of Figure 4 shows that, for a weak interchain hopping parameter (t 1 = 0.6), a phase with two MZMs is favored compared to the one MZM phase. On the other hand, when the transverse hopping is increased (Figure 4, right panel), the phase with one Majorana mode per edge is favored and the system behaves like an effective two-orbital single Kitaev chain. In order to verify the bulk-edge correspondence, in Figure 5, we present the low-energy part of the energy spectra of a ladder of N = 250 lattice sites as a function of the chemical potential µ. The analysis of the finite-size system confirms the presence and the multiplicity of zero-energy modes as prescribed by the winding number analysis.
Up to now, we have explored an extended topological phase diagram in order to emphasize the richness of the topological phase plane of a ladder model. Hereafter, in Figure 6, we specialize to the case of a condensed matter system and present phase diagrams showing a parameters region containing the device working point considered in the transport properties simulations (see Section 5). In particular, in Figure 6 (left panel), the asymmetric pairing case (i.e., ∆ = 0.02 and ∆ 1 = 0.09) is analyzed, while the equal pairing case (∆ = ∆ 1 = 0.09) is presented in the right panel. In the asymmetric situation (∆ = ∆ 1 ), a trivial region is found for t 1 ≤ 0.1 and arbitrary µ values, the latter region being reduced in the symmetric case (right panel). Let us finally note that the situation ∆ 1 > ∆, considered here, is realizable in multiorbital materials where the magnitude of different order parameters may well be different.

Quantum Transport through a Normal/Kitaev Chain/Superconductor Device
Using the scattering matrix approach, we study the transport properties of a Kitaev chain coupled to a normal and a superconducting p-wave lead ( Figure 7). The Kitaev chain and the superconducting lead define a T-stub configuration (i.e., a waveguide with a closed sidearm [22]) with the junction point located in the middle of the Kitaev chain, while the normal electrode position can be changed to simulate an STM measurement configuration. The t N parameter defines the hopping between the normal and the KC, while t S represents the hopping between the KC and the superconducting lead, while ∆ S defines the pairing amplitude. In the following, we study the topological signatures on the transport properties of the system considering a finite Kitaev chain of N = 121 sites. The scattering state of the normal electrode is given by (see Appendix A.1 for details) which represents an electron coming from the normal bulk having an ordinary reflection probability |r ee | 2 and an Andreev reflection probability given by |r eh | 2 . The wave vectors k e/h are obtained by the solution of the BdG equations: where E is the energy of the scattering process, µ N the chemical potential and t the hopping amplitude of the normal lead. Inside the superconducting electrode, the propagating state is given by (see Appendix A.2 for details): where t ee and t eh are transmission coefficients of an e-like and h-like quasiparticle. The BCS coherence factors u e,h , v e,h for a p-wave superconductor are given by: where ∆ q e/h = 2∆ sin q e/h is the momentum-dependent superconducting gap.
Resorting to a numerical solution of the scattering problem, we obtain the transmission and reflection coefficients and the wave functions of the resonant modes along the Kitaev chain (we call these modes Kitaev modes by analogy with the Majorana zero-energy modes introduced by Kitaev [2]). The wave functions of modes belonging to the Kitaev chain, being a byproduct of the solution of the scattering problem, are not normalized with respect to the position since scattering states are normalized with respect to the particle flux.
The conservation of probability of all scattering processes can be written in terms of the scattering coefficients implying the following relation: written under Andreev approximation k e/h ≈ q e/h ≈ k which has been adopted in numerical evaluations. Once the scattering coefficients are known, the conduction properties of the system can be computed adopting a BTK-like approach. Accordingly, the zero-temperature differential conductance can be written as: where A = |r eh | 2 and B = |r ee | 2 represent the Andreev reflection and the normal reflection probability, respectively [23]. Using the quantum transport theory presented above, we have studied the differential conductance and the internal resonant modes of the N-KC-SC topological device. In particular, in Figure 8, we study the differential conductance of the N-KC-SC device as a function of the N-KC link position and setting the model parameters inside the topological non-trivial region, where Majorana modes are expected to be relevant. Panels (a), (d), (g) show the conductance curves at selected N-KC link positions (specified within the figures inset) corresponding to the left end (n = 1), the middle (n = 61) and the right end (n = 121) of the nanowire. In Figure 8a, two peaks are clearly visible: a quantized zero-bias peak (E = 5 · 10 −5 ) and its satellite at higher energy (E = 6 · 10 −3 ). The analysis of the internal modes makes it evident that the quantized zero-bias peak, which is expected when the system undergoes a topological phase transition (the so-called zero-bias anomaly (ZBA)), corresponds to a MZM peaked at the left end of the wire (b). On the other hand, the peak at higher energy (E = 6 · 10 −3 ) corresponds to hybridized Majorana modes which originate a U-shaped wave function peaked at the two ends of the wire (panel (c)) and extended along the whole length of the Kitaev chain. The width of the resonant conductance peaks depends on the amplitude of the internal mode wave function at the N-KC linking point, the latter quantity being relevant in defining the overlap between the scattering states of the normal electrode and the internal modes of the Kitaev chain. The overlap increases with the wave function amplitude and depends on the normal electrode position and on the conductance peak considered. As a consequence, a narrow (broaden) peak is observed when the normal lead is weakly (strongly) coupled with the Kitaev chain internal modes, the latter situation corresponding to poor (strong) overlap between the normal electrode and the nanowire quantum states. These overlap effects are clearly visible in Figure 8a-c in which the two broadened conductance peaks originate from the non-vanishing values of the wave function amplitudes at the n = 1 position, i.e., the N-KC linking point. It is worth mentioning here that changing the linking point does not produce a qualitative change of the internal modes amplitude along the system. As a consequence, in (d) and (g), we observe the progressive and monotonic width reduction of the zero-bias peak, which follows the lowering of the wave function amplitude as the normal electrode is moved towards the right end of the Kitaev chain (panels (e) and (h)). The broadening of the satellite of the zero-bias peak is not monotonic as a function of the N-KC link position and follows the U-shaped behavior characteristic of the wave function amplitude (see (c), (f), (i)) of the internal Kitaev modes. From the experimental side, these observations suggest the possibility to reconstruct the wave function amplitude profile of the internal modes of the Kitaev chain by measuring the conductance peak broadening as the normal tip of an STM is moved along the system, the latter method being an interesting tool in characterizing the survival of topological properties under non-equilibrium conditions. It is worth mentioning here that pure Majorana modes are neutral excitations and thus they are not able to sustain a charge current. For this reason, a certain degree of hybridization of these topological modes is a necessary condition to have a finite charge current flowing through the system. The degree of hybridization of the Majorana modes inside the Kitaev chain can be quantified by studying the site-dependent charge density ρ n = | f n | 2 − |g n | 2 associated with the internal modes profile presented in Figures 8b, c, e, f, h, i. The complete analysis, which is reported in Appendix C, shows that the average charge densitȳ ρ = (∑ L n=1 ρ n )/L corresponding to the internal modes presented in (b), (e), (h) is much lower than the corresponding quantity computed for panels (c), (f), (i). This circumstance suggests that zero-bias conductance peak originates from a genuine Majorana mode presenting a weak hybridization with a hole-like non-topological mode, the latter conclusion being supported by the negative sign ofρ. On the other hand, satellite conductance peaks at higher energy are associated with internal states, presented in (c), (f) and (i), with a strong electron-like character (i.e.,ρ > 0).
As a final comment, we notice that, despite the wave guide spatial configurations corresponding to the (a) and (g) panels are one the mirror reflected of the other with respect to a reflection line coincident with the superconducting electrode, the Hamiltonian of the entire system does not satisfy this symmetry and thus the conductance curves pertaining to these cases are not coincident. This statement can be easily proven. Reflection transformation of the isolated Kitaev chain Hamiltonian (Equation (1)) can be implemented by the change of indices j → L + 1 − j. The reflection does not modify the structure of the Hamiltonian but transforms the initial pairing term ∆ into −∆. The transformed Kitaev Hamiltonian is equivalent to the initial one via a U(1) gauge transformation of the fermionic operators (a j → e i π 2 b j ). For the above reasons, physical properties of the isolated Kitaev chain are invariant under reflection. When the KC is connected to external reservoirs, the above procedure is still applicable with a different outcome. Indeed, the electrodes' Hamiltonians are invariant under reflection transformation, which is only implemented on the Hamiltonian of the Kitaev chain. On the other hand, the N-KC junction position is altered by the transformation (i.e., n = 1 → n = 121), while the position of the KC-SC link (n = 61) remains unaffected since it is a fixed point. Thus, the described reflection transformation not only provides a reflected spatial configuration of the initial system but also determines a phase difference between the superconducting pairing terms of the Kitaev chain and the superconducting electrode. Implementing a gauge transformation as the one used for the isolated Kitaev chain does not eliminate the phase gradient (which is now transferred to the N-KC and KC-SC hopping amplitudes), the latter producing observable effects. We have verified that conductance curve of the system described by the transformed Hamiltonian is coincident with the curve shown in Figure 8a, obtained for the initial Hamiltonian. Moreover, the transformed Hamiltonian does not correspond to the Hamiltonian used in Figure 8g, where the mentioned phase gradient is absent. These arguments prove our conclusions and can be equally applicable to the ladder case discussed in Section 5 (see Figure 11a-i).
Up to now, we have shown that the overlap between quantum states of the normal electrode and the internal modes of the Kitaev chain can be changed by varying the normal electrode position along the nanowire. From the experimental viewpoint, an alternative way to modify this coupling consists in changing the N-KC distance, e.g., by modifying the distance of an STM normal tip (normal electrode) from the topological nanowire to be measured. In this way, the transmission at the N-KC interface becomes exponentially suppressed at increasing distance due to the vacuum tunneling phenomenon and the same effect is induced on the overlap between quantum states belonging to the normal electrode and the Kitaev chain. Within the considered transport model, the above-mentioned effects can be considered by acting on the hopping amplitude t N , which is related to the transparency of the N-KC interface. In order to perform this analysis, in Figure 9, we present the differential conductance curves of the topological device for selected values of t N going from tunneling (t N = 0.1) to the metallic (t N = 1) regime. Direct observation of Figure 9 shows that the resonance broadening of the conductance peaks increases when t N increases, the latter condition being related to a reduced N-KC distance. On the other hand, a smaller hopping amplitude induces a resonance shrinking which reflects a reduced overlap between the tip and system modes, this condition being appropriate to describe an increased N-KC spatial gap. The analysis of Figure 9 also suggests that the transparent limit of the BTK theory, characterized by a constant subgap conductance plateau sustained by the Andreev reflection, is not accessible within the considered T-stub geometry in which additional scattering events are originated at the waveguide junctions. The BTK transparent limit can be approached using the geometry described in Appendix B, which, in general, does not correspond to the experimental conditions of an STM measurement and thus presents a limited relevance in our discussion.
To complete our analysis, let us note that the emergence of topological phases in closed systems is commonly discussed in terms of properties of the Hamiltonian spectrum. For the open system considered here, the mentioned classification can be done by studying the scattering matrix properties. Indeed, following the work by Beenakker et al. [24], we consider the reflection submatrix r of the scattering matrix S: r = r ee r eh r he r hh , (18) and introduce the topological number Q = sign[Det(r)]. We have verified that a numerical estimate of Q gives the same topological phase boundary in the t N − µ plane as the one of an isolated Kitaev chain, confirming that one can characterize the non-trivial phase of the system under study via the scattering matrix properties.

Quantum Transport through a Normal/Kitaev Ladder/Superconductor Device
In this section, we consider the transport properties of the N-KL-SC device depicted in Figure 10. Proceeding as was done for the N-KC-SC device, we have determined the scattering matrix elements by imposing appropriate matching conditions on the incident and transmitted scattering states and we have evaluated the differential conductance in terms of Andreev and normal reflection probability. The results of this analysis are shown in Figure 11 where the conductance curves and the internal modes of the Kitaev ladder are shown. The model parameters used in the computations are those that determine the emergence of a topological phase with two Majorana fermions per edge for the closed system [19] discussed in Section 3 ( Figure 6). The conductance curves ((a), (e), and (i)) show a multiple peaks structure containing a zero-bias conductance peak accompanied by two satellite peaks at higher energy (E = 4.5 · 10 −3 and E = 6.5 · 10 −3 ). The zero-bias conductance peak is related to a Majorana mode peaked at the left end of the ladder, while satellite peaks are related with U-shaped internal modes peaked at the ladder edges (see (c), (d), (g), (h), (m), and (n)). Internal modes pertaining to the satellite peaks can be associated with quantum states coming from hybridization processes of genuine Majorana modes. The hybridization of such states, which is favored by the opening of the system, produces a degradation of the initial topological properties and gives rise to quasi-Majorana states [25]. As discussed before, the degree of hybridization of the ladder internal states can be deduced by measuring the resonance width of the conductance peaks. In close analogy with the discussion reported for the N-KC-SC device, we do observe a clear correspondence between the internal modes wave function at the N-KL interface and the resonance broadening of the conductance peaks. In the next section, we describe disorder effects on quasi-Majorana states. (m): n = 121, E = 4.5 · 10 −3 . (n): n = 121, E = 6.5 · 10 −3 . The remaining model parameters have been fixed as: ∆ = 0.02, t N = t S = 0.2, t = 1, µ = 0.5, µ S = µ N = 0, ∆ 1 = 0.09, t 1 = 0.6.

Disorder Effects in a Normal/Kitaev Ladder/Superconductor Device
In order to test the robustness of the topological phase in the N-KL-SC device, we include a random potential in the ladder model which emulates, for instance, the intrinsic effect of impurities and inhomogeneities in semiconducting nanowires. On the other hand, when Majorana quasiparticles are realized using ordered assemblies of magnetic atoms on the surface of conventional superconductors [26,27], extrinsic disorder effects can be induced by the random strength of the atom-surface coupling. We model these effects assigning random values of the on-site potential U n at any system site; the random realizations have been obtained generating random numbers U n with uniform distribution in the interval (−δ, δ), with δ fixed as δ = 0.02, 0.04, 0.08 ( Figure 12). Accordingly, the statistical average of the on-site potential is zero, while its variance is given by U 2 n = δ 2 /3. The remaining model parameters have been fixed as in panel (a) of Figure 11. When the conductance of a single disordered realization is analyzed (panels (a), (e), (i)), we observe the progressive degradation of the three-peaks structure observed in Figure 11a, which has been obtained in the absence of disorder. Disorder has a different impact depending on the peak considered. In particular, the zero-bias peak, which corresponds to a robust Majorana mode, remains quantized and disorder only produces a broadening effect on the resonant peak. The above result reflects the fact that disorder introduces a finite mean free path (related to the variance of the random on-site potential) which renormalizes the quasi-particles lifetime and induces a resonance broadening. The two satellite peaks at quasi-zero energy (see Figure 11a) are fragile against disorder effects. Indeed, in the presence of disorder, one peak is completely suppressed, while the one at higher energy appears broadened and no more quantized (at high disorder). The presented phenomenology suggests that quasi-zero energy modes are related to the hybridization of genuine Majorana states. Such states, being distributed along the whole system, are quite fragile to disorder effects. The suppression of the intermediate resonant peak in the presence of disorder is easily explained looking at panel (c) of Figure 11 and (c), (g) and (m) of Figure 12. Indeed, the internal mode in Figure 11c (in the absence of disorder) is peaked at the system edges. Accordingly, the latter state is well coupled with the scattering states of the normal electrode. When the disorder is included, the overlap of this state with the normal electrode states goes to zero, producing the disappearance of the conductance peak. This effect originates from the disorder-assisted localization of the ladder resonant mode at a system edge which is not connected with the normal electrode (see (c), (g), and (m)). The latter effect is not observed for quasi-zero energy mode at higher energy, which starts to be severely degraded only at high disorder (δ = 0.08). This is confirmed by the internal modes analysis which has been shown in panels (d), (h), and (n). In this case, the analysis provides evidence that internal modes are peaked at both the system edges and thus are well coupled with the normal electrode states, the latter being a requisite to observe a finite conductance.
The above findings suggest that topological protection of quantum states can be lost as an effect of the opening of the isolated quantum system. Indeed, the connection of a topological system with external reservoirs produces hybridization of the internal modes, which is the prerequisite to manifest fragility against the detrimental effects of disorder.

Conclusions
We have addressed the problem of the topological phase characterization of open quantum systems by studying the paradigmatic cases of a Kitaev chain and a Kitaev ladder. We have characterized the topological phase diagram for the closed systems and we have found that, differently from the case of a single Kitaev chain, a Kitaev ladder displays topological phases with two Majorana states per edge. Concerning the topological phases with multiple Majorana states, the question of their stability against the opening of the system arises. In order to prove the stability of the topological phases of open systems subject to measuring processes (e.g., tunneling spectroscopy measurements), we have studied the quantum transport through a Kitaev chain and a Kitaev ladder by coupling them to a normal and a superconducting electrode. Using a lattice version of the usual Bogoliubov-de Gennes scattering theory, we have studied the differential conductance of these devices and we have found a correspondence between the conductance peaks and the resonant states confined inside the topological scattering region (i.e., the Kitaev chain or ladder). It has been proven, both for the Kitaev chain and ladder, that zero-energy peaks of the differential conductance correspond to robust Majorana states, while quasi-Majorana quantum states are signaled by almost zero-energy resonant peaks. Quasi-Majorana states originate from the hybridization of Majorana states, the latter being favored by the coupling with the scattering states of the electrodes. By studying the differential conductance of a topological device based on a Kitaev ladder, we have demonstrated that quasi-Majorana states and the associated conductance peaks are fragile against disorder effects. These findings are relevant in studying topological systems with multiple Majorana states, which are supposed to be important in multimode quantum wires or in multi-orbital topological systems. from the Fermi level. Requiring translational invariance, we can write the Nambu spinor in the following form: Solving the Equation (A1) using the ansatz in Equation (A2), we obtain the electron-and hole-like independent solutions, namely ψ e/h (j), with the corresponding energy E: Using the dispersion relations for electron and hole solutions of Equation (A3), we obtain the wave vectors of main text k e/h = arccos (µ ± E)/(−2t) (Equation (13)). An electron-like (hole-like) solution describes a particle with group velocity v g = 1 h ∂E k ∂k parallel (anti-parallel) to the particle wavevector k e (k h ).
Using ξ q = ± E 2 − ∆ 2 q , Equation (A8) can be presented in the form: Introducing the BCS coherence factors recalled in the main text, we obtain four different solutions for the BdG equations corresponding to e-like and h-like excitations, namely (u e , ±iv e ) t and (v h , ±iu h ) t , characterized by different propagation direction (±). The above arguments justify the form of the scattering states used in Equation (14) of the main text.

Appendix B. Conductance Lowering Effects in Branched Quantum Waveguides
The conduction properties of the device described in the main text (see Figure 7) are affected by the t N value, describing the interface opacity, and by the network geometry. Indeed, the so-called T-stub configuration in quantum wave guides enhances scattering events and may induce a lowering of the device conductance. According to these considerations, one easily argues that the transparency of the scattering region (i.e., the Kitaev chain) not only depends on the t N value, which plays the same role of the dimensionless barrier strength Z of the BTK theory [23], but it is also affected by the wave guide geometry. As a consequence, it is plausible that the transparent limit, in which the conductance is maximized, cannot be obtained within the relevant experimental configuration described in the main text.
For the sake of completeness, in order to study the transparent limit of the proposed model, we define a new system geometry presenting an in-line configuration which avoids branching points along the current path (see Figure A1).
We numerically solve the tight binding equations for this model and compute the differential conductance. The results of the analysis are reported in Figure A2, where analogies with the conductance plots of the BTK theory [23] are present. Differently from the usual BTK conductance curves derived for the s-wave superconducting pairing, a zero-bias peak, reminiscent of the p-wave coupling, is always present in our simulations. A similar analysis can be performed for the N-KL-SC device ( Figure A3). Compared to the differential conductance of the single Kitaev chain (see Figure A2), a peculiar two-peaks structure is present in the ladder differential conductance curve. The two-peaks structure is more pronounced for opaque N-KL interfaces (i.e., for low t N values), while, for higher t N values, the two-peaks structure is substituted by a peculiar conductance minimum originated by the incomplete coalescence of two broadened resonant peaks. Moreover, the zero-bias peak and its satellite at higher energy, constituting the two-peaks structure mentioned before, are related to the presence of a quantum state with Majorana character accompanied by a quasi-Majorana quantum state, as discussed in the main text.  Figure A2. N-KC-SC device. Zero-temperature differential conductance of the model depicted in Figure A1. The model parameters have been fixed as: ∆ = 0.02, t S = t = 1, µ = 0.5, µ S = µ N = 0, while t N = 0.1, 0.4, 0.8, 1 has been used in obtaining the different plots.

Appendix C. Charge Neutrality of Quasi-Majorana Modes
Majorana modes in closed systems are neutral excitations. Electrical neutrality implies that they cannot transport charge current unless a certain degree of hybridization with other quantum states occurs. These hybridization phenomena are supposed to be relevant when the topological system is connected with external reservoirs and a current is forced to flow through it. Under this circumstance, Majorana modes are expected to be subject to hybridization with the reservoir quantum states in order to warrant the current conservation. The degree of hybridization is expected to be variable with the system geometry, which defines the path followed by the charge current along the device. On the light of these general arguments, in Figure A4, we have investigated the site-dependent charge density ρ n = | f n | 2 − |g n | 2 (with f n and g n the electron-like and hole-like component of the Nambu spinor) corresponding to the internal modes presented in Figure 8b,c,e,f,h,i. Different ρ n versus n curves in Figure A4 present a complicate oscillating behavior with the position along the Kitaev chain with average charge density given byρ = (∑ L n=1 ρ n )/L. Our findings confirm that zero-bias peaks in the conductance curves are related to Majorana modes which are weakly hybridized with the electrodes quantum states, the degree of hybridization being related to the deviation of the internal mode charge from the charge neutrality condition. On the other hand, satellite peaks in the differential conductance curves (Figure 8a,d,g) are related to internal modes characterized by a relevant deviation from the charge neutrality requirement, the latter condition providing evidence of a strong contamination of the Majorana character of the internal mode.  Hybridized modes with prevalent electron-like (hole-like) character are defined byρ > 0 (ρ < 0).