# Majorana bound states and non-local spin correlations

in a quantum wire
on an unconventional superconductor

###### Abstract

We study theoretically the proximity effect of a one-dimensional metallic quantum wire (in the absence of spin-orbit interaction) lying on top of an unconventional superconductor. Three different material classes are considered as a substrate: (i) a chiral superconductor in class D with broken time-reversal symmetry; a class DIII superconductor (ii) with and (iii) without a nontrivial number. Interestingly, we find degenerate zero energy Majorana bound states at both ends of the wire for all three cases. They are unstable against spin-orbit interaction in case (i) while they are topologically protected by time-reversal symmetry in cases (ii) and (iii). Remarkably, we show that non-local spin correlations between the two ends of the wire can be simply controlled by a gate potential in our setup.

###### pacs:

74.45+c, 74.20.RpIntroduction.– Proximity effects of superconductors and normal metals as well as semiconductors have been a subject of continued interest due to the rich physical phenomena of these hybrid systems. Especially, the topological nature of the superconducting proximity effect is an issue of current research activities Fu and Kane (2008); Beenakker (2011); Tanaka et al. (2012); Alicea (2012) stimulated by the possible occurrence of Majorana fermions in solid state systems Majorana (1937); Wilczek (2009); Franz (2010); Read and Green (2000); Ivanov (2001); Das Sarma et al. (2005); Nayak et al. (2008). For example, the interface between an -wave superconductor and the surface state of a three-dimensional topological insulator (TI) is predicted to offer a platform to realize Majorana fermions Fu and Kane (2008). The underlying idea here is holographic principle that entails the reduction of electronic degrees of freedom, i.e., electron fractionalization, at the surface of a topologically non-trivial bulk state. Furthermore, the proximity effect of a TI to unconventional superconductors has also been studied theoretically Linder et al. (2010). Another (highly interesting) proposal is to use semiconductors with Rashba spin-orbit interaction in combination with -wave superconductors for this purpose. In the presence of a sufficiently strong magnetic field, an inverted gap opens at the -point in this system and topological superconductivity as well as Majorana fermions may appear when the Fermi energy lies inside the gap Sato et al. (2009); Sau et al. (2010a); Alicea (2010); Oreg et al. (2010). A recent experiment observed the zero-bias conductance peak at the edge of an InSb quantum wire (QW) coupled to an -wave superconducting substrate, which might be the first experimental observation of a Majorana fermion Mourik et al. (2012). It has also been shown theoretically that a one-dimensional Rashba quantum wire coupled to a -wave superconductor hosts doubly degenerate Majorana bound states (MBSs) at the edge Wong and Law (2012). Besides these distinct works, there are various other research activities on the basis of QWs (for example Refs. Sau et al. (2010b); Lutchyn et al. (2010); Potter and Lee (2011); Oreg et al. (2010); Tewari et al. (2011); Stanescu et al. (2011); Lutchyn et al. (2011); Kim et al. (2012)). In all the cases listed above, spin-orbit interaction and/or the Zeeman energy in the QW is essential. Therefore, an open question is if zero energy Majorana bound states are also realizable in a one-dimensional system in the absence of these interactions, which opens new opportunities as we will show below.

In this paper, we investigate electronic states caused by the proximity effect between a metallic QW and a two-dimensional (2D) unconventional superconductor to search for zero energy MBSs at the ends of the QW. We discover that the resulting MBSs are doubly degenerate, i.e., characterized by a spin degree of freedom. This leads to nonlocal spin-correlations between the two ends of the QW which can be manipulated by all-electric means. For a substrate in symmetry class D Schnyder et al. (2008); *NewJPhys.12.065010, i.e., in the absence of time reversal symmetry (TRS) the degenerate end states can be gapped out by switching on TRS preserving local imperfections that couple opposite spin. This is reflected in the fact these systems are characterized by a trivial topological invariant. In contrast, for a TRS preserving substrate in symmetry class DIII, we find helical MBS pairs that are topologically protected by TRS.

Model.– Our setup is shown schematically in Fig. 1 and the model Hamiltonian given by

(1) |

(2) | ||||

(3) |

(4) |

where describes the QW put along -direction on the surface of the superconductor described by which infinitely spreads in the -plane. These two are connected by a tunneling Hamiltonian with the spin-independent constant matrix element . We consider the following three cases for the order parameter matrix : (i) for the chiral superconductor of class D with the topological invariant , and (ii) for the time reversal symmetric class DIII superconductor with the topological invariant . Here, the denotes Pauli matrices in spin space. As the third case (iii), we consider a system where two layers of a 2D superconductor of case (ii) are coupled to each other by an interlayer transfer integral, which is a topologically trivial case with . The model (i) is relevant to SrRuO where the chiral superconductivity with broken time-reversal symmetry is believed to be realized Mackenzie and Maeno (2003). The model (ii) is the 2D analogue of the He B phase, and is continuously connected to the helical non-centrosymmetric topological superconductor with Rashba spin-orbit interaction Tanaka et al. (2009); Sato and Fujimoto (2009); Nakosai et al. (2012). We present the calculational steps and results mainly for case (i) below, but it is straightforward to apply them to the other two cases (ii) and (iii). Our main purpose is to investigate the electronic states in the one-dimensional superconductivity formed in this way and to discuss its topological nature as well as resulting Majorana bound states.

Effective Green’s function.– First, we derive the effective Green’s function of the QW by integrating over the superconductor

(5) | ||||

(6) |

where denotes Pauli matrices in the Nambu space. As long as there is neither spin-orbit interaction in the wire nor spin-dependent tunneling between the wire and the superconductor, theinduced superconducting order parameter is also spin-triplet and can be decoupled into two sectors . The integration can be analytically performed and we obtain

(7) | |||

(8) | |||

(9) | |||

(10) |

where is a momentum cut-off and , , likewise. In the following, we look at the sector, which is equivalent to a spinless -wave superconductor. From Eq. (8), we can then calculate the spectral function of the electron as

(11) |

where is the retarded Green’s function obtained from via the analytic continuation , where is a infinitesimally small positive number.

In Fig. 2, we plot Eq. (11) with , and several values of and . The insets show the schematic structure of the spectrum. The red lines represent states in the wire without and with an induced superconducting gap, in cases (a) and (b) respectively. The blue shaded area shows the continuum of the bulk states in the superconductor, which arises in the effective Green’s function of the QW as the self-energy (Eqs. (5) and (6)). The green lines are edge states in the 2D SC that will be further discussed below. Figures 2(a) and (b) correspond to the cases where the band of the QW is outside and inside of the superconducting gap of the substrate superconductor, respectively. We can see that in Fig. 2(b) a finite gap is induced by the proximity effect. An interesting situation arises when the magnitude of the transfer integral related to the tunneling is strong enough that bonding and anti-bonding states are formed at the interface of the wire and the superconductor. The energy of those states is pushed away from the low energy regime, and they act as a potential barrier to other states. Therefore, a boundary is effectively induced in the substrate superconductor. It is known that edge modes appear at the boundary of 2D topological superconductors. The energy of the midgap states around approach zero as the magnitude of becomes larger (see Fig. 2(c)). This shows the formation of edge channel that are gapped out because the edge modes located at both sides of the wire weakly interact with each other. In Fig. 2(d) we show the spectral function for an -wave superconductor substrate for comparison. There are no states corresponding to the edge channel in case (c).

Topological nature.– Now, we focus on the topological nature of the QW system. For this purpose, we use the method for the calculation of the topological number in terms of Green’s function at zero frequency Wang and Zhang (2012); Budich and Trauzettel (2012). We derive an effective Hamiltonian from Eq. (5) on the basis of for a given . This 22 Hamiltonian is similar to that of Kitaev’s original proposal of a 1D topological superconductor Kitaev (2001). It is known that the system is in the topological phase when connects antipodal points of the unit sphere as is varied from the center to the boundary of the Brillouin zone Alicea (2012). Therefore, it is easy to verify that when , the QW becomes two copies of Kitaev’s 1D topological superconductor with . The topological nature of the 1D superconductor manifests itself in bound states at the ends of the QW, which will benumerically studied in the following. However, in case (i), the two identical copies of the -nontrivial Kitaev model form a trivial composite system. Physically, this implies that the degenerate end states can be gapped out by switching on terms that couple the opposite spin sectors as we confirmed by numerical calculations.

Numerical study of a tight-binding model and Majorana bound states.- We study a tight-binding model by numerical calculations. We change our continuum model to a lattice model, i.e., and in Eq. (8), and construct the corresponding tight-binding Hamiltonian. Then, we calculate the energy spectrum in the system with open boundary conditions as depicted in Fig. 1. Figure 3(a) shows the topological phase transition as a function of the chemical potential in the QW. In the strong coupling regime, where , there are no states in the induced superconducting gap, while in the weak coupling regime, where there are zero energy states. It can be checked that four-fold degenerate zero energy states are formed, two at each end. Fig. 3(b) shows the probability distributions of one of them. Red solid circles represent the weight in the QW, and we can see the state is localized at the ends of the wire. In sharp contrast to the Kitaev model, we have two spin sectors in the QW which are degenerate. This fact attaches a spin degree of freedom to the MBSs. However, in the presence of spin-orbit interaction, which mixes the sectors, the degeneracy is lifted and the Majorana fermions will be pushed away from zero energy.

In the following, we briefly discuss the other two possibilities (ii) and (iii) for the SC substrate. We have confirmed that a similar analysis to case (i) outlined above applies to these cases, and we have found that the doubly degenerate zero energy Majorana bound states appear when the chemical potential of the QW is positive and large enough for the topological phase transition to occur. In cases (ii) and (iii), the spin degenerate pair of MBS stems from a non-trivial invariant for the composite system in symmetry class DIII. Therefore, these end states are topologically protected by TRS and cannot be gapped out by switching on spin-orbit interaction as long as the bulk gap is maintained. This is an essential difference to case (i) where the conservation of was needed in addition to the generic symmetries of the model to obtain zero energy MBSs, while TRS protects the doubly degenerate MBSs in cases (ii) and (iii). This is particularly remarkable for case (iii) where the substrate 2D superconductor is originally in a topologically trivial phase.

Local and non-local fermions.– There have been many theoretical proposals to use chiral MBSs for topological quantum computing Sau et al. (2011a, b); Halperin et al. (2012). Chiral MBSs are difficult to detect and to manipulate because they only couple in a limited way to their environment. This is different in our proposal. On the one hand, the doubly degenerate MBSs discussed above are more susceptible regarding their coupling to the environment. On the other hand, they can show novel responses through a strong coupling between the fermion parity and the spin degree of freedom. The idea behind is schematically shown in Fig. 4. The solid lines represent Kitaev chains with MBSs indicated by red circles. Let us elaborate a bit more on the two kinds of fermions that can be formed by our MBSs. Since Majorana fermions are real, normal fermions are formed from pairs of them. We can construct two kinds of normal (complex) fermions, non-local spin polarized ones and , and local ones and . The local pseudo spin operators are defined as . It is known that Majorana fermions have only Ising-type spin due to their anticommutation relation , then, and from the above definition. We would like to emphasize that these local variables are specific to the paired Majorana states. There will be various interesting physics as a consequence of the synthesis of the local and the non-local nature. Here, as an example, we focus on spin-spin correlation of two Ising spins formed at the ends. In preparation for that, we first review some aspects of the fermion parity. When we have two Majorana operators, and , there exist two orthogonal states and , which satisfy and , respectively. The former has an even fermion parity and the latter odd. They can be rewritten by normal fermion operators defined as where and . Note that if the wire is connected to ground by a capacitor, one can modify the fermion number via charging energy where is the charge, is the gate voltage, and . The importance of charging energy in the context of topological superconductors was first discussed in Ref. Fu (2010), where electron teleportation (mediated by the non-locality of the MBSs) has been proposed. We now apply this idea to our model. Our Hilbert space is 4-dimensional and there are two kinds of bases one can construct, i.e., from , , or from , and . Namely, one can define , etc., and , , …are local basis states while , , …are non-local basis states. The relation between these two kinds of basis states is given by

(12) |

In our system, the total number of fermions can be controlled by the gate voltage . We consider correlation in the odd and even fermion number space, respectively. The density matrices are and . Then we obtain

(13) |

where represents the statistical average with the density matrices and , respectively.

These results show that by modulating charging energy we can control spin-spin correlation, which is independent of the distance between the two spins because they are connected through the superconducting wire in between. This effect is a consequence of both, the local and the non-local degrees of freedom of the system. One can also manipulate the two Ising spins by modulating the gate voltage as a function of time as . The equations of motion of and by the driving force yield

(14) | ||||

(15) |

for a constant during some time period . By tuning to satisfy one can reverse the pseudospins components when they are anti-parallel while they remain unchanged when they are parallel.

Conclusion.– In this paper, we have studied theoretically the proximity effect of a one-dimensional metallic quantum wire without spin-orbit interaction on the substrate of an unconventional superconductor. We have considered three different cases for the substrate: (i) a chiral superconductor; (ii) a class DIII superconductor with a non-trivial number, and (iii) a class DIII superconductor with a trivial number. Interestingly, we have found degenerate zero energy Majorana bound states at both ends of the quantum wire for all three cases. We have shown that these degenerate Majorana bound states for cases (ii) and (iii) can be nicely used for a combined spin/parity qubit. Furthermore, we have demonstrated that the resulting non-local spin correlations between the two ends of the wire can be controlled by the gate voltage potential acting on the wire. One of the candidate systems for the case (ii) substrate are thin films of SrRuO Tada et al. (2009) or bi-layer Rashba system Nakosai et al. (2012), where the class DIII superconducting state might be realized, and for the case (iii) super lattice structure of CeCoIn Mizukami et al. (2011).

###### Acknowledgements.

Acknowledgment.-SN was supported by Grant-in-Aid for JSPS Fellows and JCB by the Swedish Research Council. This work was further supported by Grant-in-Aid for Scientific Research (S) (Grants No. 24224009) from the Ministry of Education, Culture, Sports, Science and Technology of Japan, Strategic International Cooperative Program (Joint Research Type) from Japan Science and Technology Agency, Funding Program for World-Leading Innovative RD on Science and Technology (FIRST Program), the Deutsche Forschungsgemeinschaft, the European Science Foundation, and the Helmholtz Foundation.## References

- Fu and Kane (2008) L. Fu and C. L. Kane, Phys. Rev. Lett. 100, 096407 (2008).
- Beenakker (2011) C. W. J. Beenakker, arXiv:1112.1950 (2011).
- Tanaka et al. (2012) Y. Tanaka, M. Sato, and N. Nagaosa, J. Phys. Soc. Jpn. 81, 011013 (2012).
- Alicea (2012) J. Alicea, Rep. Prog. Phys. 75, 076501 (2012).
- Majorana (1937) E. Majorana, Nuovo Cimento 14, 171 (1937).
- Wilczek (2009) F. Wilczek, Nature Phys. 5, 614 (2009).
- Franz (2010) M. Franz, Physics 3, 24 (2010).
- Read and Green (2000) N. Read and D. Green, Phys. Rev. B 61, 10267 (2000).
- Ivanov (2001) D. A. Ivanov, Phys. Rev. Lett. 86, 268 (2001).
- Das Sarma et al. (2005) S. Das Sarma, M. Freedman, and C. Nayak, Phys. Rev. Lett. 94, 166802 (2005).
- Nayak et al. (2008) C. Nayak, S. H. Simon, A. Stern, M. Freedman, and S. Das Sarma, Rev. Mod. Phys. 80, 1083 (2008).
- Linder et al. (2010) J. Linder, Y. Tanaka, T. Yokoyama, A. Sudbø, and N. Nagaosa, Phys. Rev. Lett. 104, 067001 (2010).
- Sato et al. (2009) M. Sato, Y. Takahashi, and S. Fujimoto, Phys. Rev. Lett. 103, 020401 (2009).
- Sau et al. (2010a) J. D. Sau, R. M. Lutchyn, S. Tewari, and S. Das Sarma, Phys. Rev. Lett. 104, 040502 (2010a).
- Alicea (2010) J. Alicea, Phys. Rev. B 81, 125318 (2010).
- Oreg et al. (2010) Y. Oreg, G. Refael, and F. von Oppen, Phys. Rev. Lett. 105, 177002 (2010).
- Mourik et al. (2012) V. Mourik, K. Zuo, S. M. Frolov, S. R. Plissard, E. P. A. M. Bakkers, and L. P. Kouwenhoven, Science 336, 1003 (2012).
- Wong and Law (2012) L. M. Wong and K. T. Law, arXiv:1211.0338 (2012).
- Sau et al. (2010b) J. D. Sau, S. Tewari, R. M. Lutchyn, T. D. Stanescu, and S. Das Sarma, Phys. Rev. B 82, 214509 (2010b).
- Lutchyn et al. (2010) R. M. Lutchyn, J. D. Sau, and S. Das Sarma, Phys. Rev. Lett. 105, 077001 (2010).
- Potter and Lee (2011) A. C. Potter and P. A. Lee, Phys. Rev. B 83, 094525 (2011).
- Tewari et al. (2011) S. Tewari, T. D. Stanescu, J. D. Sau, and S. D. Sarma, New J. Phys. 13, 065004 (2011).
- Stanescu et al. (2011) T. D. Stanescu, R. M. Lutchyn, and S. Das Sarma, Phys. Rev. B 84, 144522 (2011).
- Lutchyn et al. (2011) R. M. Lutchyn, T. D. Stanescu, and S. Das Sarma, Phys. Rev. Lett. 106, 127001 (2011).
- Kim et al. (2012) Y. Kim, J. Cano, and C. Nayak, arXiv:1208.3701 (2012).
- Schnyder et al. (2008) A. P. Schnyder, S. Ryu, A. Furusaki, and A. W. W. Ludwig, Phys. Rev. B 78, 195125 (2008).
- Ryu et al. (2010) S. Ryu, A. P. Schnyder, A. Furusaki, and A. W. W. Ludwig, New J. Phys. 12, 065010 (2010).
- Mackenzie and Maeno (2003) A. P. Mackenzie and Y. Maeno, Rev. Mod. Phys. 75, 657 (2003).
- Tanaka et al. (2009) Y. Tanaka, T. Yokoyama, A. V. Balatsky, and N. Nagaosa, Phys. Rev. B 79, 060505 (2009).
- Sato and Fujimoto (2009) M. Sato and S. Fujimoto, Phys. Rev. B 79, 094504 (2009).
- Nakosai et al. (2012) S. Nakosai, Y. Tanaka, and N. Nagaosa, Phys. Rev. Lett. 108, 147003 (2012).
- Wang and Zhang (2012) Z. Wang and S.-C. Zhang, Phys. Rev. X 2, 031008 (2012).
- Budich and Trauzettel (2012) J. C. Budich and B. Trauzettel, arXiv:1207.1104 (2012).
- Kitaev (2001) A. Y. Kitaev, Phys.-Usp. 44, 131 (2001).
- Sau et al. (2011a) J. D. Sau, D. J. Clarke, and S. Tewari, Phys. Rev. B 84, 094505 (2011a).
- Sau et al. (2011b) J. D. Sau, B. I. Halperin, K. Flensberg, and S. Das Sarma, Phys. Rev. B 84, 144509 (2011b).
- Halperin et al. (2012) B. I. Halperin, Y. Oreg, A. Stern, G. Refael, J. Alicea, and F. von Oppen, Phys. Rev. B 85, 144501 (2012).
- Fu (2010) L. Fu, Phys. Rev. Lett. 104, 056402 (2010).
- Tada et al. (2009) Y. Tada, N. Kawakami, and S. Fujimoto, New J. Phys. 11, 055070 (2009).
- Mizukami et al. (2011) Y. Mizukami, H. Shishido, T. Shibauchi, M. Shimozawa, S. Yasumoto, D. Watanabe, M. Yamashita, H. Ikeda, T. Terashima, H. Kontani, and Y. Matsuda, Nature Phys. 7, 849 (2011).