# Chiral Topological Superconductors Enhanced by Long-Range Interactions

###### Abstract

We study the phase diagram and edge states of a two-dimensional p-wave superconductor with long-range hopping and pairing amplitudes. New topological phases and quasiparticles different from the usual short-range model are obtained. When both hopping and pairing terms decay with the same exponent, one of the topological chiral phases with propagating Majorana edge states gets significantly enhanced by long-range couplings. On the other hand, when the long-range pairing amplitude decays more slowly than the hopping, we discover new topological phases where propagating Majorana fermions at each edge pair nonlocally and become gapped even in the thermodynamic limit. Remarkably, these nonlocal edge states are still robust, remain separated from the bulk, and are localized at both edges at the same time. The inclusion of long-range effects is potentially applicable to recent experiments with magnetic impurities and islands in 2D superconductors.

###### pacs:

74.20.Mn,03.65.Vf,71.10.Pm,05.30.Rt1. Introduction.— Topological superconductors are novel quantum phases of matter rmp1 (); rmp2 (); LibroBernevig () with unconventional pairing symmetries that host gapless Majorana zero modes. These are exotic quasiparticles that can realize non-Abelian statistics Read_et_al00 () unlike bosons or fermions. Remarkably, their intrinsic topological protection and robustness against disorder make them ideal candidates as building blocks in quantum computation rmp3 (); rmp4 (); Alicea_et_al_11 (); Baranov_et_al_13 () and quantum information processing Kitaev01 (); Mazza_et_al13 ().

As first shown in Fu_Kane_2008 (), the surface state of a 3D topological insulator can host Majorana zero modes by proximity coupling a conventional s-wave superconductor. Topological superconductivity can also be induced by replacing the topological insulator with a semiconductor with strong spin-orbit interaction Sau_et_al_2010 (); Sau2_et_al_2010 (); Alicea_2010 (); Sau3_et_al_2010 (); Oreg_et_al_2010 (). Indeed, several experiments in proximity-induced superconducting nanowires Mourik_et_al12 (); Deng_et_al12 (); Das_et_al12 (); Albrecht_et_al16 () and topological insulator-superconductor heterostructures Wang_et_al12 (); He_et_al14 (); Sun_et_al16 (); He_et_al17 () have shown signatures for the presence of Majorana zero modes and chiral edge states. In all these platforms, the superconducting pairing amplitude is of purely short-range nature.

An alternative approach to realizing chiral p-wave Hamiltonians is to deposit magnetic atoms on top of a conventional s-wave superconductor substrate Nadj_et_al14 (); Pawlak_et_al15 (). For a periodic 1D array of these magnetic impurities Pascual_et_al16 (); Ruby_et_al17 (), a p-wave Hamiltonian with intrinsic long-range pairing Nadj_et_al13 (); Pientka2013 (); Klinovaja_et_al13 (); Pientka2014 (); Li_et_al_16 (); Kaladzhyan_et_al16 () is induced, provided that the length of the chain is small compared to the coherence length of the host superconductor Pientka2013 (). Novel topological phase diagrams Viyuela_et_al16 (); Pachos_16 (); Luca_16 (); Alecce_17 () have been analyzed for certain 1D Hamiltonians with long-range interactions Niui_12 (); DeGottardi_13 (); Vodola_et_al14 (); Gong2015_1 (); Gong2015_2 (); Tudela_15 (); Vodola_et_al16 (); Vodola_et_al17 (); Dutta_17 (). Interestingly, massless Majorana end modes were shown to pair up into a topological massive Dirac fermion Viyuela_et_al16 () localized simultaneously at the two ends of a 1D long-range p-wave superconductor. More recently, spatially extended Shiba states Yu_65 (); Shiba_68 (); Rusinov_69 () of magnetic impurities in 2D superconductors have been experimentally observed Menard_et_al15 (); Menard_et_al16 (); Heinrich_et_al17 (). When they conform to planar arrays Ronty_et_al_15 (); Li_et_al_16_2D (); Kaladzhyan_et_al16_2 (), long-range couplings between the impurity states emerge, raising the question of whether new topological phases in 2D long-range systems can be found.

Motivated by previous experiments with magnetic atoms Nadj_et_al14 (); Pascual_et_al16 (); Menard_et_al15 (); Heinrich_et_al17 () and islands Menard_et_al16 (), we compute the topological phase diagram (see Fig. 1) and the edge states of a 2D p-wave superconductor with long-range hopping and pairing couplings. We find topological phases and quasiparticles that were absent in the usual short-range model Read_et_al00 (); LibroBernevig (). When the power-law decay for hopping and for pairing amplitudes is the same, we observe different phenomena. (A) Enhanced Topological Chiral Phase (ETCP): a topological chiral phase with stable propagating Majorana edge states even for very strong long-range effects. (B) Suppressed Topological Chiral Phase (STCP): the topological phase with opposite chirality is slowly diminished and even disappears completely in the infinite-range limit. (C) Chiral Symmetry Inversion (CSI): a topological phase transition between the two chiral phases is driven when tuning the decaying exponent. On the other hand, when the pairing amplitude decays more slowly than the hopping, there is a critical value of the decaying exponent at which propagating edge states become gapped. Remarkably, these states remain chiral and localized at the edges, which leads to intrinsic protection against backscattering.

2. Long-range chiral topological superconductor.— We consider a 2D lattice of spinless fermions with long-range hopping and pairing amplitudes,

(1) |

where the sum runs over all integers and except . The function is the Euclidean distance between fermions and is the chemical potential, with a certain origin offset for convenience. The hopping amplitude decay with exponent of the distance , and the pairing amplitude

(2) |

decay with exponent . Without loss of generality, we fix in the rest of the paper. Eq. (2) is the most natural extension of a p-wave symmetry when considering long-range effects. Moreover, Hamiltonian (1) belongs to the D symmetry class of topological insulators and superconductors Ludwig (); Kitaev_2009 (), whatever and . Fluctuation effects over the mean field Hamiltonian have not been considered in the present work.

At large values of the decaying exponents (), Eq. (1) has a well-defined limit and we recover a short-range chiral p-wave Hamiltonian Volovik_99 (); Read_et_al00 (); Ivanov_01 (); LibroBernevig (). The system hosts vortices with non-Abelian anyonic statistics Ivanov_01 () that are of great relevance in proposals for topological quantum computation rmp3 (). In addition, propagating Majorana modes appear at the boundary, displaying a linear dispersion at low momentum that connects the particle and hole bands (see Fig. 2a).

Now we analyze the fate of these topologically protected edge states when we add long-range hopping and pairing amplitudes as in Eq. (1). To this end, we assume periodic boundary conditions in the direction, where the Fourier-transformed operators are

(3) |

with the number of sites along the direction. Thus, we can rewrite Hamiltonian (1) as , where

(4) | |||

with

(5) |

Since we assume periodic boundary conditions, the definition of the distance function and the pairing function need to be slightly modified by replacing whenever . We are now in a position to study the edge state structure of Hamiltonian (4). We consider two different cases: i) where both hopping and pairing amplitudes decay with the same exponent, and ii) where the pairing term is longer range than the hopping.

3. Chiral Majorana edge states.— Let us first consider the case where long-range hopping and pairing terms decay with the same exponent . The motivation comes from the experimental and theoretical results regarding magnetic atoms on top of s-wave topological superconductors both in one-dimensional chains Nadj_et_al13 (); Pientka2013 (); Pientka2014 (); Li_et_al_16 (); Nadj_et_al14 () and planar structures Ronty_et_al_15 (); Li_et_al_16_2D (); Menard_et_al15 (); Heinrich_et_al17 (), where long-range effects have been shown to play an even more important role than in 1D systems.

In what follows, we show the link between the edge state physics and the different topological phases. We demonstrate that the topological phase diagram gets modified by the effect of long-range couplings (see Fig. 1), creating an imbalance between different chiral phases. We assume periodic boundary conditions in both spatial directions, where the Hamiltonian (1) takes the form in the Nambu-spinor basis of paired fermions with , where are the Pauli matrices and

(6) | |||||

(7) |

Here, is the energy dispersion relation and is the so-called winding vector, expressed in terms of series of trigonometric functions

(8) | |||||

(9) | |||||

(10) |

where the sum runs over all integers and except , and is the Euclidean distance on the -torus.

We define a winding number in terms of the winding vector , which reads as

(11) |

where BZ stands for Brillouin zone. The winding number is an integer topological invariant that characterizes a continuous mapping between the crystalline momentum and the winding vector . This is also called the Pontryagin index. The topological transition points can be identified out of the gap-closing points of given in Eq. (6),

(12) | |||||

As shown in Fig. 1, for , the topological phase diagram coincides with the short-range chiral p-wave Hamiltonian Read_et_al00 (); LibroBernevig (), with , , . Thus, for , the Majorana modes have a definite chirality (), whereas for they have the opposite (). We now discuss the consequences of adding long-range effects that significantly change the boundaries and the phases of the short-range case. If the chemical potential or , then and the system is in a trivial phase with no edge states. However, when , then and there are propagating edge states with definite chirality. As shown in Fig. 2a the edge state at each boundary cross at when connecting the particle and hole bands. Interestingly enough, this topological phase gets significantly enlarged for , since the divergence of the function at the point shifts the location of the transition point towards as the system becomes long-ranged, leading to the ETCP mechanism. This is of great relevance for experimental realizations of p-wave superconductors and Majorana edge states without the need of a precise fine-tuning of the chemical potential. Physically, the effect is explained due to a broadening of the quasiparticle kinetic energy with respect to the short-range case, as shown in Fig. 2b. As we decrease the long-range exponent , the energy peak at the point is more pronounced; thus, the range of chemical potential at which topological superconductivity occurs gets enlarged.

A distinct topological phase can be found if . This phase is characterized by a winding number (see Fig. 1) and the presence of propagating edge states with opposite chirality. As shown in Fig. 3a the edge states at each boundary cross at . Surprisingly, this topological phase gets progressively suppressed when the long-range terms become important (STCP). This creates an effective asymmetry between the two chiral sectors, whose mathematical explanation is two-fold. First, in the short-range limit (), the two sectors and of Hamiltonian (1) are connected through the unitary transformation that inverts the chirality. However, the long-range terms break this symmetry leading to an asymmetry between the two sectors. The second reason is that is an alternating and convergent series at the high-symmetry points and ; thus, according to Eq. (12), the critical points and always remain finite. In fact, in the infinite range limit , they take on the same value and the topological phase involving the edge crossing at ends disappearing, as shown in Fig. 1. Finally, we can induce topological phase transitions between the two different chiral phases by tuning the decaying exponent at a fixed chemical potential. This mechanism, dubbed chiral symmetry inversion (CSI) is a distinct feature of these long-range topological systems remark (). We notice that for the pairing amplitude decays faster than the hopping, thus, the edge states never become gapped. Moreover, the phase boundaries in Fig. 1 only depend on the exponent (see Eq. (12)), and therefore would remain unaltered.

4. Non-local gapped edge states.— We now consider the situation where the long-range pairing amplitude decays slower than the hopping, i.e., . This situation could be achieved by periodic driving of a topological superconductor, as shown in Benito_et_al14 () for the 1D Kitaev chain Kitaev01 (). We demonstrate that depending on the pairing exponent and the chirality of the topological phase, the Majorana edge states become gapped despite remaining localized and protected at the boundary.

Without loss of generality, we assume . First, we focus on the topological phase where and Chern number . We have previously shown in Fig. 2a that there are propagating edge states crossing at , provided . But the situation is different for . By means of finite-size-scaling, we numerically show that when the chiral edge states crossing at become gapped due to purely long-range effects as depicted in Fig. 3b. This result coincides with a divergent group velocity at the point (with periodic boundary conditions), signaling strong long-range effects within this phase. In Fig. 3c, we plot the spatial distribution of the edge-states wave function for the mode proving edge localization.

The edge localization in the direction, and the bulk-edge gap depicted in Fig. 3 are robust against random fluctuations of the chemical potential in the direction. Disordered impurities in the direction could lead to elastic scattering between states with momentum and , respectively, which implies backscattering and localization of the propagating edge states. However, the chiral nature of the massive edge states protects them against backscattering since the modes at and are localized at different edges. It is only at very low momentum that the two modes at the edges hybridize, as shown in Fig. 3c of the manuscript. This is a difference with respect to the short-range case. However, this hybridization happens only very close to , thus, the density of states that are subject to disorder effects is extremely small. Moreover, at large momentum , the situation is quite the opposite. For strong long-range effects, the edge states remain localized at each edge and do not hybridize, whereas for the short-range model the edge states merge into bulk states. Hence, at large momentum our new edge states are more robust against disorder and back-scattering than the propagating Majorana modes in the short-range model. In addition, these states are topologically protected by a non-trivial Chern number .

On the other hand, if the Majorna edge states are stable. Second, we analyze the phase with and Chern number . We find propagating edge states crossing at , no matter the exponent . We note that at the group velocity is not divergent even for very strong long-range; thus, the edge-crossing point is protected as shown in Fig. 3a. In summary, we have proven that the long-range pairing induces a gap in the edge states depending on the chirality of the topological phase. It is also interesting to analyze the trivial region , where in the short-range limit there are no edge states. However, we numerically find that below protected gapped edge states are formed out of originally trivial bulk states. These are now localized at the two edges at the same time, being non-local. The energy dispersion relation and the edge localization are similar to Fig. 3b and Fig. 3c respectively.

5. Conclusions.— We found that long-range hopping and pairing amplitudes deeply modify the topological phase diagram of a two-dimensional p-wave topological superconductor. We have shown that when the decay exponents for the pairing and the hopping amplitudes are equal , one of the topological phases is greatly enhanced (ECTP), while the topological phase with opposite chirality is slowly suppressed (SCTP) for very strong long-range effects (see Fig 1). By tuning the decaying exponent , we can induce topological phase transitions between the two different chiral phases (CSI) at fixed chemical potential . On the other hand, if , i.e., the pairing decays slower than the hopping amplitude, rich phenomena occur depending on the exponent and the chemical potential : i) propagating edge states are topologically stable whatever the exponent ; ii) Majorana modes become gapped, though remaining localized at the edges and protected against backscattering; iii) previously trivial bulk states give rise to protected gapped edge states. The long-range couplings are motivated by recent experiments Nadj_et_al14 (); Pascual_et_al16 (); Menard_et_al15 (); Heinrich_et_al17 (); Menard_et_al16 () and theoretical proposals Ronty_et_al_15 (); Li_et_al_16_2D () of lattices of magnetic impurities and islands in 2D superconductors. In addition, the imbalance between pairing and hopping decaying exponents could be achieved through Floquet driving fields Benito_et_al14 (). The enhancement of topological phases by long-range couplings can be of importance for the observation and manipulation of Majorana modes without the need of a precise fine-tuning of the chemical potential.

###### Acknowledgements.

We thank Leonid Glazman and Tharnier Oliveira for valuable conversations. M.A.MD. and O.V. thank the Spanish MINECO grant FIS2012-33152, the CAM research consortium QUITEMAD+ S2013/ICE-2801. O.V. thanks Fundación Rafael del Pino, Fundación Ramón Areces and RCC Harvard. The research of M.A.M.-D. has been supported in part by the U.S. Army Research Office through Grant No. W911N F-14-1-0103. The work at MIT is supported by the DOE Office of Basic Energy Sciences, Division of Materials Sciences and Engineering under Award DE-SC0010526. LF is partly supported by the David and Lucile Packard Foundation.## References

- (1) M. Z. Hasan and C. L. Kane, Rev. Mod. Phys. 82, 3045 (2010).
- (2) X.-L. Qi and S.-C. Zhang, Rev. Mod. Phys. 83, 1057 (2011).
- (3) B. A. Bernevig and T. L. Hughes, Topological Insulators and Topological Superconductors (Princeton University Press, New Jersey, 2013).
- (4) N. Read and D. Green, Phys. Rev. B 61, 10267 (2000).
- (5) C. Nayak, S. H. Simon, A. Stern, M. Freedman and S. Das Sarma, Rev. Mod. Phys. 80, 1083 (2008).
- (6) B. M. Terhal, Rev. Mod. Phys. 87, 307 (2015).
- (7) J. Alicea,Y. Oreg, G. Refael, F. von Oppen and M. P. A. Fisher, Nat. Phys. 7, 412-417 (2011).
- (8) C. V. Kraus, P. Zoller, M. A. Baranov, Phys. Rev. Lett. 111, 203001 (2013).
- (9) A. Y. Kitaev, Phys.-Usp. 44, 131 (2001).
- (10) L. Mazza, M. Rizzi, M. D. Lukin, J. I. Cirac, Phys. Rev. B 88, 205142 (2013).
- (11) L. Fu and C. L. Kane, Phys. Rev. Lett. 100, 096407 (2008).
- (12) J. D. Sau, R. M. Lutchyn, S. Tewari, and S. Das Sarma, Phys. Rev. Lett. 104, 040502 (2010).
- (13) J. Alicea, Phys. Rev. B 81, 125318 (2010).
- (14) J. D. Sau, S. Tewari, R. M. Lutchyn, T. D. Stanescu and S. Das Sarma, Phys. Rev. B 82, 214509 (2010).
- (15) R. M. Lutchyn, J. D. Sau, and S. Das Sarma, Phys. Rev. Lett. 105, 077001 (2010).
- (16) Y. Oreg, G. Refael, and F. von Oppen, Phys. Rev. Lett. 105, 077002 (2010).
- (17) V. Mourik, K. Zuo, S. M. Frolov, S. R. Plissard, E. P. A. M. Bakkers, L. P. Kouwenhoven, Science 336, 6084 (2012).
- (18) M. T. Deng, C. L. Yu, G. Y. Huang, M. Larsson, P. Caroff, H. Q. Xu, Nano Lett. 12, 6414-6419 (2012).
- (19) A. Das, Y. Ronen, Y. Most, Y. Oreg, M. Heiblum and H. Shtrikman, Nat. Phys. 8, 887-895 (2012).
- (20) S. M. Albrecht, A. P. Higginbotham, M. Madsen, F. Kuemmeth, T. S. Jespersen, J. Nygard, P. Krogstrup and C. M. Marcus, Nature 531, 206-209 (2016).
- (21) M-X. Wang et al., Science 336, 6077 (2012).
- (22) Q. L. He et al., Nat. Comm. 5, 4247 (2014).
- (23) H.-H. Sun et al., Phys. Rev. Lett 116, 257003 (2016).
- (24) Q. L. He et al., Science 357, 6348 (2017).
- (25) S. Nadj-Perge, I. K. Drozdov, J. Li, H. Chen, S. Jeon, J. Seo, A. H. MacDonald, B. A. Bernevig, A. Yazdani, Science 346, 6209 (2014).
- (26) R. Pawlak et al, npj Quantun Information 2, 16035 (2016).
- (27) D-J. Choi, C. Rubio-Verdú, J. de Bruijckere, M. Moreno Ugeda, N. Lorente, J. I. Pascual, Nat. Commun. 8, 15175 (2017).
- (28) M. Ruby, B. W. Heinrich, Y. Peng, F. von Oppen, K. J. Franke, Nano Lett. 17, 4473 (2017).
- (29) S. Nadj-Perge, I. K. Drozdov, B. A. Bernevig, A. Yazdani, Phys. Rev. B 88, 020407 (2013).
- (30) F. Pientka, L.I. Glazman and F. von Oppen, Phys. Rev. B 88 155420 (2013).
- (31) J. Klinovaja, P. Stano, A. Yazdani, and D. Loss, Phys. Rev. Lett. 111, 186805 (2013).
- (32) F. Pientka, L.I. Glazman and F. von Oppen, Phys. Rev. B 89 180505(R) (2014).
- (33) T. Neupert, A. Yazdani, B. A. Bernevig, Phys. Rev. B 93, 094508 (2016).
- (34) V. Kaladzhyan, C. Bena, P. Simon, J. Phys.: Condens. Matter 28 485701 (2016).
- (35) O. Viyuela, D. Vodola, G. Pupillo, M.A. Martin-Delgado, Phys. Rev. B 94, 125121 (2016).
- (36) K. Patrick, T. Neupert, J. K. Pachos, Phys. Rev. Lett. 118, 267002 (2017).
- (37) L. Lepori, L. Dell’Anna, New J. Phys. 19, 103030 (2017).
- (38) A. Alecce, L. Dell’Anna, Phys. Rev. B 95, 195160 (2017).
- (39) Y. Niu, S. B. Chung, C-H. Hsu, I. Mandal, S. Raghu, S. Chakravarty, Phys. Rev. B 85, 035110 (2012).
- (40) W. DeGottardi, M. Thakurathi, S. Vishveshwara, D. Sen, Phys. Rev. B 88, 165111 (2013).
- (41) D. Vodola, L. Lepori, E. Ercolessi, A. V. Gorshkov, G. Pupillo, Phys. Rev. Lett. 113, 156402 (2014).
- (42) Z.-X. Gong, M. F. Maghrebi, A. Hu, M. L. Wall, M. Foss-Feig, A. V. Gorshkov, Phys. Rev. B 93, 041102 (2016).
- (43) Z.-X. Gong, M. F. Maghrebi, A. Hu, M. Foss-Feig, P. Richerme, C. Monroe, A. V. Gorshkov, Phys. Rev. B 93, 205115 (2016).
- (44) A. González-Tudela, C.-L. Hung, D. E. Chang, J. I. Cirac and H. J. Kimble, Nature Photonics 9, 320-325 (2015).
- (45) D. Vodola, L. Lepori, E. Ercolessi, G. Pupillo, New J. Phys. 18 015001 (2016).
- (46) L. Lepori, A. Trombettoni, D. Vodola, J. Stat. Mech. 033102 (2017).
- (47) A. Dutta, A. Dutta, Phys. Rev. B 96, 125113 (2017).
- (48) L. Yu, Acta Phys. Sin. 21, 75 (1965).
- (49) H. Shiba, Prog. Theor. Phys. 40 435 (1968).
- (50) A.I. Rusinov, JETP Lett. 9, 85 (1969).
- (51) G. C. Menard et al, Nat. Phys. 11, 1013 (2015).
- (52) G. C. Menard et al, arXiv: 1607.06353 (2016).
- (53) B. W. Heinrich, J. I. Pascual, K. J. Franke, arXiv: 1705.03672 (2017).
- (54) J. Röntynen, T. Ojanen, Phys. Rev. Lett. 114, 236803 (2015).
- (55) J. Li, T. Neupert, Z. Wang, A.H. MacDonald, A. Yazdani, B. A. Bernevig, Nat. Comm. 7, 12297 (2016).
- (56) V. Kaladzhyan, J. Röntynen, P. Simon, and T. Ojanen, Phys. Rev. B 94, 060505 (2016).
- (57) A. P. Schnyder, S. Ryu, A. Furusaki and A. W. W. Ludwig, Phys. Rev. B 78, 195125 (2008).
- (58) A. Kitaev, AIP Conf. Proc. 1134, 22 (2009).
- (59) G. E. Volovik, JETP Lett. 70, 609 (1999).
- (60) D. A. Ivanov, Phys. Rev. Lett. 86, 268 (2001).
- (61) This is in contrast to the standard topological transition between different chiral phases tuning the chemical potential.
- (62) M. Benito, A. Gómez-León, V. Bastidas, T. Brandes, G. Platero, Phys. Rev. B 90, 205127 (2014).