# Majorana fermions in an out-of-equilibrium topological superconducting wire: an exact microscopic transport analysis of a p-wave open chain coupled to normal leads

###### Abstract

Topological superconductors are prime candidates for the implementation of topological-quantum-computation ideas because they can support non-Abelian excitations like Majorana fermions. We go beyond the low-energy effective-model descriptions of Majorana bound states (MBSs), to derive non-equilibrium transport properties of wire geometries of these systems in the presence of arbitrarily large applied voltages. Our approach involves quantum Langevin equations and non-equilibrium Green’s functions. By virtue of a full microscopic calculation we are able to model the tunnel coupling between the superconducting wire and the metallic leads realistically; study the role of high-energy non-topological excitations; predict how the behavior compares for increasing number of odd vs. even number of sites; and study the evolution across the topological quantum phase transition (QPT). We find that the normalized spectral weight in the MBSs can be remarkably large and goes to zero continuously at the topological QPT. Our results have concrete implications for the experimental search and study of MBSs.

Topological quantum computing Nayak08 is interesting not only because of the technological goal but also because it requires the realization and manipulation of new types of quasiparticle excitations that have never been observed before and the question of their actual existence tests the limits of our understanding of the quantum world at a fundamental level. The simplest proposals along this line involve the realization of Majorana fermions (MFs).

One of the pressing theoretical problems in this context is to sort out the experimental signatures of MFs. Several proposals for their realization in (effective) one-dimensional (1D) topological (p-wave) superconductors have been put forward recently Lutchyn10 ; *Oreg10; *Nersesyan11; *Sau11. A natural way to probe superconducting (SC) nanowires is via transport measurements and STM spectroscopy. Such transport studies of mid-gap features in the spectrum of unconventional superconductors have consistently attracted considerable attention Kashiwaya00 and particularly so in the quasi-one-dimensional case of (possibly p-wave) organic superconductors Sengupta01 ; Ha03 . The renewed interest in MBSs has given rise to more detailed studies including some exciting recent experimental claims of MBS detection Sasaki11 .

Important insights have been gained from past theoretical works Bolech07 ; Nilsson08 ; Flensberg10 that study low-energy effective models consisting of the MBSs and an effective coupling to and between them as proposed in Ref. Bolech07, . These works have made specific predictions for detecting MBSs in experiments. For instance, non-equilibrium electrical transport between a metallic probe and a collection of coupled MBSs has been derived Bolech07 ; Nilsson08 ; Flensberg10 and the differential electrical conductance calculated in the limit where the applied voltage () is much smaller than the SC energy gap. It is crucial to go beyond low-energy approximations and calculate the full current-voltage characteristics including the contribution of the non-topological higher-energy excitations. It is also important to check the regime of validity of these low-energy models and even better have a full microscopic calculation to predict and interpret experimental results in realistic systems.

While different proposals for realizing MBSs in 1D differ in details Lutchyn10 ; *Oreg10; *Nersesyan11; *Sau11, all essential features of realistic systems are expected to be minimally described by the so called Kitaev chain (linked to the early theoretical proposals for topological quantum computation) Kitaev00 . We thus study non-equilibrium electrical transport in Kitaev’s tight-binding model of a single-channel 1D quantum wire with triplet (p-wave) SC pairing between electrons with the same spin orientation (equal-spin pairing). The model’s spectrum includes two MBSs at the two ends of the wire, which are expected to exhibit non-Abelian braiding statistics (in wire networks Alicea11 ; *Clarke11; *Halperin11), are separated from the non-topological complex excitations by an energy gap and have been proposed as a basis for a decoherence-free topological quantum memory. In view of theoretically capturing the physics of realistic experiments, we are interested in studying non-equilibrium transport through the 1D quantum wire in an open-quantum-system geometry, i.e., the wire is connected to external baths which incorporate both decoherence and dissipation. This is a minimal setup that can capture the essential features of recent proposals for realizing p-wave SC wires Lutchyn10 ; *Oreg10. We develop a novel extension of the quantum Langevin Equations and Green’s Function (LEGF) approach Dhar03 to SC wires with BCS-type pairing. LEGF provides us with an appropriate exact method to determine the electrical current in the SC wire with arbitrary values of as compared with the strengths of the contacts and the SC pairing amplitude. Although this method has proven quite useful to investigate electrical Dhar03 ; DharSen06 ; Roy07 as well as thermal DharRoy06 non-equilibrium transport in non-interacting, ordered and disordered quantum systems connected to baths, in this letter we will focus only on the electrical transport.

The Hamiltonian of the full system is given by where

(1) | ||||

The term denotes the Hamiltonian of the 1D wire and that of the bath, while the coupling between them is . Here and denote standard fermion operators on the SC wire and on the bath, respectively. The p-wave BCS-pairing amplitude of the wire is and the parameter can be tuned to drive the SC wire in or out of a topological phase (cf. Refs. Kitaev00, ; Motrunich01, ).

Next we introduce a local re-parametrization for the electron operators on the SC wire in terms of MFs Kitaev00 , , with and . The wire Hamiltonian in the MF basis (see Fig. 1) is given by . The two MFs at any one end of the wire are connected to the end site of the bath with equal-strength coupling. (Note that the contact coupling to the inner MF of the end site of the wire is outside the scope of all previous works, as they deal either with effective MBS models Bolech07 ; Nilsson08 ; Flensberg10 or a continuum representation Sengupta01 .) We assume that for times , the wire is disconnected from the baths (each in equilibrium at a specified temperature and chemical potential ). At time we connect both baths to the 1D SC wire and we are interested in the steady-state properties, which can be derived by integrating out the bath operators using the quantum LEGF formalism and considering the limit . We get the following solution

(2) | |||

where and same for . The single-particle retarded matrix Green’s function of the bath is given by and . Its first diagonal element (in the MF ‘site’ basis) is used to define , which corresponds to dissipation due to coupling with the bath. The ‘noise’ in the leads,

(3) |

depends on the bath’s initial distribution and obeys DharSen06

(4) |

where is the Fermi distribution function. The above identity is a form of fluctuation-dissipation relation. We use the definition where is the local density of states at the first site on the bath. Finally, is the Green’s function of the full system (wire plus baths) DharSen06 .

The outward current between the wire and the baths is given by

(5) |

where denotes averaging over noise like in Eq. (4). For the sake of concreteness and simplicity of presentation, we will consider the ‘symmetric’ case for which the two baths to are similar () and connected with the SC wire by identical contacts ( or ) while and . In this case the steady-state currents are equal up to a sign and

(6) |

We have checked that the expression in Eq. (6) reduces to the known result for the case of a normal metallic wire () where is the transmission coefficient of a single electron through the wire. For a SC wire, that contains contributions to transport coming from both electrons and holes, can be simply interpreted in the symmetric case via the zero-temperature differential conductance (with , )

(7) |

Moreover, for , is an even function and thus proportional to .

At and , there exist two exactly zero-energy localized MBSs (decoupled from rest of the wire) and the corresponding zero-temperature has a peak at zero bias and two more peaks in at coming from all the degenerate high-energy excitations at these parameter sets. As can be seen from Fig. 2(a-f), the shows distinct behaviors for finite depending on whether is odd (panels e,f) or even (panels a-d) once . We have developed a physical understanding of different behaviors by closely looking at the behavior of a finite Kitaev chain with odd vs. even (see Fig. 1). For odd , there are pairs of degenerate modes including a pair of zero-energy modes thereby showing peaks (including one at ). For even , all modes, including the ones at zero energy are split and the zero-temperature has a dip at zero-bias. However, the splitting in the MBSs decreases with increasing wire length, and finally goes to zero in the long-wire limit as with for . (This is in contrast to the result of an effective model of a chain of coupled MBSs, for which the zero-bias conductance is zero for the case of two MBSs Flensberg10 .) Note that the splitting of the non-topological high-energy excitations survives for a long even-site wire as can be seen by comparing heights in Fig. 2(a-d).

For a finite , the high-energy excitations are split for any even or odd number of sites ^{1}^{1}1We have checked that the isolated spikes are not numerical artifacts and that they can be traced back to the structure of .. Gating the wire () also splits any degenerate zero modes at for the finite but, remarkably, in the large limit, they regain their integrity and yield a zero-bias peak. (The splitting for finite- case can be understood by adding a matrix element between the two MFs within each site in Fig. 1.) Next, in the long-wire case, we closely study how evolves as is varied across the topological QPT [i.e. the quantum phase transition from the topological () to non-topological () SC phase].
We observe that as approaches , the gap shrinks and the zero-bias peak corresponding to the MBSs loses weight; at the gap closes completely; and for the gap reopens with no zero-bias anomaly [see Fig. 2(h)].

In Fig. 3(a) we compare the I-V characteristics of the wire for different parameters (varying ). In the topological phase, displays a sharp rise for small voltages due to the presence of the low-energy MBSs and then saturates with increasing voltage reflecting the gap in the energy-spectrum. The gap and hence the plateau is most sharply defined for . The current increases again for even higher voltages, when non-topological excitations in the SC wire start to contribute to the transport, and finally saturates again when the applied voltage is larger than the band-width of the wire. For comparison, we plot in Fig. 3(b) the I-V of the topologically trivial SC phase along with the corresponding metallic phase with the same .

We find that the normalized spectral weight of the MBSs initially decreases with increasing wire length [see Fig. 4(a)] but quickly saturates to a length-independent finite value (cf. with the ideas of electronic teleportation in these systems Semenoff07 ; *Fu10) which depends on and but is independent of the strength of the contacts. As can be seen in Fig. 4(b), the maximum spectral weight for long wires is (a remarkably large fraction carried by a single pair of states) and is achieved when and when the two MBSs have exactly zero energy, are completely localized (one at each end) and decoupled from the rest of the wire (see Fig. 1). The contribution of the MBSs (which can be interpreted as the order parameter of the topological phase) decreases for increasing and at the topological QPT goes to zero in the way shown in Fig. 4(b).

If we depart from the symmetric case considered above by setting , we find that gives an I-V similar to that of the symmetric case while becomes zero (in agreement with Ref. Akhmerov11, ). We have also verified that moderate asymmetry between the left and right baths or the respective contacts does not significantly affect the results for the current as compared to the symmetric case and is thus meaningful to relate these to existing or future measurements.

In conclusion, by carrying out an exact microscopic calculation of non-equilibrium electrical transport in a p-wave (Kitaev) chain connected to baths, we have been able to go beyond the previous works using either effective models of MBSs Bolech07 ; Nilsson08 ; Flensberg10 or a continuum representation such as the Bogoliubov-deGennes formalism Sengupta01 . First, we have captured the interplay and differences between Majorana and regular complex fermions (like the Bogoliubov excitations in the superconductor) and overcome the limitation of having voltages smaller than the SC gap. Second, we have modeled the contacts more realistically and with greater detail (going beyond the scope of scattering-formalism approaches; cf. Refs. Nilsson08, ; Sengupta01, ). Third, we have shown how the microscopic properties and correspondingly the differential conductance varies with the length of the SC quantum wire for different parameter regimes. Fourth, we have compared and contrasted the behaviors when the number of sites in the wire is even vs. odd. Fifth, we have calculated the normalized spectral-weight contribution of the MBSs to the total differential conductance and thereby provided a starting point to meaningfully compare the contributions to the transport coming from MBSs vs. other excitations in realistic experiments. And finally, we have predicted how the differential conductance changes across the topological QPT and shown the way in which the spectral-weight of the MBSs continuously goes to zero at the QPT.

As shown in this letter, all these points not only make concrete experimental predictions, but also have strong implications for the scope and validity of the low-energy effective models. Key features of the experimental curves Ha03 ; Sasaki11 are in agreement with our results, in particular the overall shape of the zero-bias conductance peak and the pronounced dips on its sides (which reflect the plateaux shown in Fig. 3). We suggest systematic measurements of the differential conductance based on our predictions and, in particular, the analysis of the spectral weight of the zero-bias conductance peak (such study would require a careful material-dependent protocol to take into account the effect of higher bands and dimensionality when defining the fractional spectral weight) as a way to strengthen and advance the experimental efforts to observe MBSs.

###### Acknowledgements.

We acknowledge the financial support by the University of Cincinnati. CJB & NS are grateful for the hospitality of the Tata Institute of Fundamental Research (TIFR) where the last part of the research was done and a major part of the manuscript written.## References

- (1) C. Nayak, S. H. Simon, A. Stern, M. Freedman, and S. Das Sarma, Rev. Mod. Phys. 80, 1083 (Sep 2008)
- (2) R. M. Lutchyn, J. D. Sau, and S. Das Sarma, Phys. Rev. Lett. 105, 077001 (Aug 2010)
- (3) Y. Oreg, G. Refael, and F. von Oppen, Phys. Rev. Lett. 105, 177002 (Oct 2010)
- (4) A. A. Nersesyan and A. M. Tsvelik, Europhys. Lett. 96, 17002 (Aug 2011)
- (5) J. D. Sau and S. Das Sarma, ArXiv e-prints(Nov. 2011), arXiv:1111.6600 [cond-mat.supr-con]
- (6) S. Kashiwaya and Y. Tanaka, Rep. Prog. Phys. 63, 1641 (Oct 2000)
- (7) K. Sengupta, I. Žutić, H.-J. Kwon, V. M. Yakovenko, and S. Das Sarma, Phys. Rev. B 63, 144531 (Mar 2001)
- (8) H. Ha, J. Oh, J. Moser, and M. Naughton, Synthetic Metals 137, 1215 (Apr 2003)
- (9) S. Sasaki, M. Kriener, K. Segawa, K. Yada, Y. Tanaka, M. Sato, and Y. Ando, Phys. Rev. Lett. 107, 217001 (Nov 2011)
- (10) C. J. Bolech and E. Demler, Phys. Rev. Lett. 98, 237002 (Jun 2007)
- (11) J. Nilsson, A. R. Akhmerov, and C. W. J. Beenakker, Phys. Rev. Lett. 101, 120403 (Sep 2008)
- (12) K. Flensberg, Phys. Rev. B 82, 180516 (Nov 2010)
- (13) A. Y. Kitaev, Usp. Fiz. Nauk (Suppl.) 171, 131 (Oct 2000), English: Phys.-Usp. 44 131 (2001).
- (14) J. Alicea, Y. Oreg, G. Refael, F. von Oppen, and M. P. A. Fisher, Nature Phys. 7, 412 (Feb 2011)
- (15) D. J. Clarke, J. D. Sau, and S. Tewari, Phys. Rev. B 84, 035120 (Jul 2011)
- (16) B. I. Halperin, Y. Oreg, A. Stern, G. Refael, J. Alicea, and F. von Oppen, ArXiv e-prints(Dec. 2011), arXiv:1112.5333 [cond-mat.supr-con]
- (17) A. Dhar and B. Sriram Shastry, Phys. Rev. B 67, 195405 (May 2003)
- (18) A. Dhar and D. Sen, Phys. Rev. B 73, 085119 (Feb 2006)
- (19) D. Roy and A. Dhar, Phys. Rev. B 75, 195110 (May 2007)
- (20) A. Dhar and D. Roy, J. Stat. Phys. 125, 801 (Nov 2006)
- (21) O. Motrunich, K. Damle, and D. A. Huse, Phys. Rev. B 63, 224204 (May 2001)
- (22) We have checked that the isolated spikes are not numerical artifacts and that they can be traced back to the structure of .
- (23) G. W. Semenoff and P. Sodano, J. Phys. B 40, 1479 (Mar 2007)
- (24) L. Fu, Phys. Rev. Lett. 104, 056402 (Feb 2010)
- (25) A. R. Akhmerov, J. P. Dahlhaus, F. Hassler, M. Wimmer, and C. W. J. Beenakker, Phys. Rev. Lett. 106, 057001 (Jan 2011)