Platform of chiral Majorana edge modes and its quantum transport phenomena

We propose a method to create two-dimensional topological superconductors with a heterostructure of ferromagnet (FM), topological insulator (TI) thin film and superconductor, in which the two surfaces of the TI thin film are treated as a two-dimensional system. One of surfaces is superconducting due to proximity effect and the other feels an exchange field from the FM. We show that there is a topological phase with single chiral Majorana edge mode that exists in readily achievable parameter regions and does not require magnetization to be small. An experimental setup is proposed based on our model to uniquely determine the existence of Majorana chiral modes using a Josephson junction. Also, we show that multiple chiral Majorana edge modes may appear when unconventional superconductors are used.

On the other hand, 1D propagating Majorana modes at the edges of 2D TSCs have been realized recently in an experiment [46] that combines quantum anomalous Hall insulators (QAHIs) with superconductors (SCs) [47]. In such a system, the TSC with single chiral Majorana edge mode appears only when the out-of-plane magnetization, M z , is small compared to the superconducting gap ∆. This parameter regime may be achieved by tuning M z with external magnetic field. In this way, narrow regions of half-quantized conductance were observed, which are considered as a signature of the chiral Majorana mode [48][49][50].
However, controversy arises based on an observation that half-quantized conductance may also appear trivially without any Majorana mode [51,52] if the edge states of the QAHI go through a long-enough path in the SC so that the SC part behaves just like a metal connecting two QAHIs in series. Due to the requirement of very small magnetization and the fact that the QAHI was obtained by magnetic doping (which induces disorder), domains are likely to form and long paths for the QAHI edge states can exist. This difficulty of distinguishing the TSC explanation from the trivial one stems from the theoretical model where both surfaces of the TI have the same magnetization that competes directly against FIG. 1. Schematic FM/TI/SC heterostructure. The two surfaces of the TI thin film have Dirac dispersions. One of them is superconducting due to proximity effect and the other feels an out-of-plane exchange field from the FM. the superconductivity order parameter, resulting in the requirement M z < ∆ [47]. Thus, Majorana systems beyond this limitation are desired.
In this letter, we propose an alternative method to create 2D TSCs with a heterostructure of ferromagnet (FM), topological insulator (TI) thin film and superconductor (SC), in which the two surfaces of the TI thin film form a two-dimensional system (FIG.1). One surface is superconducting due to proximity effect and the other feels an exchange field from the FM. We show that there is a topological phase with single chiral Majorana edge mode that exists in readily achievable parameter regions and does not require magnetization to be smaller than the SC gap. An experimental setup containing a Josephson junction is proposed to uniquely determine the existence of Majorana chiral modes, in which a smoking-gun evidence is a change of conductance from 1/2 σ 12 1 to σ 12 = 1/3 (the unit of conductance is e 2 /h throughout this manuscript) as the current is tuned up across the Josephson critical current. We also show that multiple chiral Majorana edge modes may appear when unconventional superconductors are used.
Model -Assuming the FM to be insulating, the low-energy properties of the system is determined by the two surface states of the TI thin film. The two surfaces, however, experience different environments. The bottom surface is in good contact with the FM and feels a strong out-of-plane exchange field, whereas the top surface becomes superconducting due to proximity of the SC, as shown in FIG.1. In the Nambu basis {c t,↑ (k), c t,↓ (k),

the effective Hamiltonian of the system is
The Pauli matrices s x,y,z , σ x,y,z , τ x,y,z act on spin, layer and particle-hole spaces respectively. M z is the exchange field felt by the bottom layer, ∆(k) and d(k) are the singlet and triplet SC order parameters on the top layer. The constant t c is the hybridization energy between the two layers, which depends on the film thickness. µ is the chemical potential while δE denotes the energy shift between two surfaces. We assume the decay length of the exchange field to be smaller than the TI film thickness so that the top surface does not couple with the exchange field directly, although an indirect coupling can be conveyed by the hybridization t c .
The only discrete symmetry of this system is the redundant particle-hole symmetry and thus it belongs to D class [53,54] in which topological phases in two dimensions can be identified by Chern numbers. This model describes a currently accessible experimental system where the heterostructure can be fabricated by molecular beam epitaxy (MBE) method with controllable thickness [55]. Note that it becomes a QAHI if the SC on the top surface is replaced by another FM same as the bottom surface [56]. This property is useful when we discuss junctions of QAHIs and TSCs.
With s-wave SCs -In the limit t c = 0, the stacked TI surfaces may be unfolded and regarded as two sections of one surface aligned side by side. The edge of the stacked system then becomes the boundary between the sections. When they are gapped by superconductivity and magnetization respectively, a chiral Majorana mode appears at the boundary [58]. Similarly, such a mode is expected at the edge of the heterostructure in FIG.1. The phase diagram obtained by calculating the total Berry curvature γ is shown in FIG.2(a). In normal states, γ corresponds to Hall conductance. If a bulk gap exists, it is quantized so that N c = γ/2π with the integer N c being the Chern number [57]. In superconducting states, the Chern number can be defined in the same mathematical way using the Bogoliubov-de Gennes Hamiltonian (although it is no longer related to Hall conductance). For |µ| < |M z |, the system is gapped and N c = sign[M z ], corresponding to a single chiral Majorana edge mode whose chirality is determined by the magnetization direction. If |µ| > |M z |, the bottom surface becomes gapless and the total Berry curvature is not quantized. When t c = 0, it competes with M z trying to generate a trivial hybridization gap. Consequently, larger |M z | is required to obtain a non-zero Chern number at µ = 0, as shown in FIG.2(b). If M z > 0 is small and µ deviates from zero, the system first enters a non-trivial region with N c = sign[M z ] and then transitions into a trivial phase again for large |µ|, as shown by the dashed line. This is easily understood by looking at the normal state band structures as shown in the left panel of FIG.2(d). The bands are non-degenerate due to the inter-layer hybridization t c and the exchange field M z . The two subbands closest to zero energy have a trivial gap opened by t c and thus SC has no effect if µ = 0, giving N c = 0. As |µ| increases, the Fermi level cuts a single band and the resulting SC is topological with single chiral Majorana edge mode. As |µ| further increases, the Fermi level cuts two bands and the system becomes a trivial SC again.
This even-odd effect of the number of Fermi surfaces resembles that of Rashba systems [7,29]. The similarity becomes clearer in band structure when a difference in chemical potentials of the top and bottom surfaces, δE = 0.2, is considered. As shown in the right panel of FIG.2(d), the conduction (valence) band looks just like a Rashba band with positive (negative) mass and a Zeeman field which splits the degeneracy at k = 0. However, the splitting of the valence band is larger than that of the conduction band. This is because the states of the conduction (valence) band near k = 0 are from the top (bottom) surface of the TI film, and the top surface only couples with the FM indirectly through t c while the bottom surfaces feels the exchange field directly. Consequently, the topological region (where the number of Fermi surfaces is odd) is wider when µ < 0, as shown by the dashed line in FIG.2(c) where the whole phase diagram is obtained. The effect of δE has been discussed [59] and it happens in real experiments [60]. Note that no new phases emerge from this energy difference of two surfaces and its effect on the phase diagram is only quantitative.
Experimental detection -The chiral Majorana modes in the heterostructure discussed above may be detected by a QAHI-TSC-QAHI junction as shown by previous researches [46,48,49]. To avoid trivial explanations [51,52] and uniquely determine the existence of chiral Majorana modes, however, we propose an experimental setup including a Josephson junction [61,62]. As shown in FIG.3(a), the setup is achieved by adding SCs and FMs alternately on the top surface of the TI thin film while attaching a uniform FM to the bottom surface. Regions with both top and bottom surfaces coupled to FMs are QAHIs [56] while those with one surface coupled to FM and the other to SC are TSCs with single Majorana edge mode, as we have discussed. The two TSCs form a Josephson junction, which are connected to external electrodes 1 or 2 through another QAHI at each end. Note that the SCs are not grounded, contrary to that of reference [62]. Figure 3(b) schematically shows how the edge states propagate in Majorana basis in which each normal edge state of the QAHIs is regarded as two Majorana states. As we shall see in the following, a simple measurement of the current-dependent conductance from lead 1 to lead 2 provides smoking-gun evidence for the chiral Majorana edge modes.
Consider a current I < I bulk c (I bulk c being the critical current of the bulk SCs) flowing from lead 1 to lead 2. When I > I c , where I c is the Josephson critical current, the current through the junction is normal and a voltage difference between two SCs exists. The current across the junction is carried by the normal edge states of the QAHI and thus it is determined by the occupation numbers of the edge states. We can define the chemical potential µ i (i = 1, 2, ..., 6) for each QAHI edge state as shown in FIG.3(b). They must satisfy µ 1 − µ 6 = µ 3 − µ 4 = µ 5 − µ 2 = eI due to the quantized Hall conductance of the QAHIs. Additionally, the left half of the system is a QAHI-TSC-QAHI junction whose scattering matrix has been obtained [48]. A simple application of the scattering coefficients to multi-terminal measurement leads to µ 3 = µ 6 = (µ 1 +µ 4 )/2. Similarly, we get µ 4 = µ 5 = (µ 2 +µ 3 )/2 for the right half of the system. Combination of these relations leads to This result turns out to be the same as the case where the SCs are replaced by normal metals. However, it should be emphasized that the current I we consider here is smaller than the bulk critical current of the SCs and the TSCs with chiral Majorana edge states remain intact. When I < I c , the current flowing through the junction (the middle QAHI region) is a supercurrent carried by Cooper pairs and the aforementioned relations of µ i no longer hold. There is no voltage drop between the two TSCs and they can only differ by a phase of the order parameter φ which is related to the current by I = I c sin(φ/2) or φ = 2 arcsin(I/I c ). The existence of φ can drastically affect σ 12 by inducing interference between different Majorana paths and thus changing the tunneling amplitude [62]. Taking this into account, one obtains the conductance where φ 0 = k F L is the kinetic phase acquired by the edge states across the junction [62]. (k F is the Fermi wave vector and L is the length of the junction.) When I ≪ I c , φ is small and its effect on σ 12 is negligible. As I approaches I c , φ increases from zero to π and σ 12 decreases to zero, as shown in FIG.3(c). Note that in the special case with φ 0 = π/2 + nπ, σ < 12 is constantly unity. In summary, as the current increases, the two-terminal conductance σ 12 starts with a value between 1/2 and 1 and decreases until it exceeds the Josephson critical current, above which σ 12 is constantly 1/3. This is a unique consequence of the chiral Majorana edge states in the TSCs.
Multiple Majorana modes -With the setup in FIG.1, it is interesting to consider unconventional SCs rather than s-wave [63]. Particularly, previous studies show that multiple chiral Majorana edge modes may be achieved in two-dimensional Rashba system using p-wave and d-wave SCs with Zeeman field [64][65][66].
Consider a general pairing potential∆(k) = (Ψ +d · σ)iσ y that acts on TI surface states with the Hamiltonianĥ(k) = g · σ + V z σ z − µ 0 where V z is a perpendicular Zeeman field and g = (k y , −k x , 0) is the spin-orbit coupling vector. Let us, for clarity, assume the Fermi level to be above the band crossing point (µ 0 > 0), then the superconductivity gap function in the band basis be- If only singlet (s-wave and d-wave) pairings are considered, the gap function on the top surface state is simply Ψ cos α = (∆ s + ∆ d cos 2θ) cos α. When |∆ s | > |∆ d |, it is fully gapped and the topological property is the same as the pure s-wave case discussed previously. When  .4 (c)) but not in [010] direction (FIG.4 (a)). However, when M z is turned on, the change of the state is quite different from that of usual d-wave SCs. Particularly, for large M z , dispersive edge modes appear on [010] edges, as shown in FIG.4(b). When p-wave pairing is included, we consider three typical cases, d = ∆ p (sin θ, − cos θ, 0), d ⊥ = ∆ p (cos θ, sin θ, 0) and d ′ = (sin θ, cos θ, 0). In the first two cases, the d-vectors are either parallel or perpendicular to the spin-orbit vector g. For d , the gap function is (∆ p + Ψ cos α) so that p-wave behaves the same as s-wave. For d ⊥ , the gap function becomes ∆ p sin α − iΨ cos α, which is always fully gapped when ∆ p V z = 0. In the heterostructure, it can support single chiral Majorana mode. With the third choice, d ′ , the gap function is ∆ p sin 2θ sin α + i(∆ p − ∆ d cos α) cos 2θ. Note that, when ∆ p V z = 0, this is also fully gapped except at ∆ p = ∆ d cos α where the gap changes sign. A phase diagram with this pairing potential (∆ p = 0.3, ∆ p = 0.1) is obtained in FIG.5, where each phase is identified by a Chern number that is equal to the number of chiral Majorana modes. A maximum of four Majorana modes can exist in this system.

Discussion -
We have proposed a platform of Majorana edge channels by using superconductor/topological insulator/ferromagnet (SC/TI/FM) heterostructures. The topological phase is much easier to be realized compared with the setups studied thus far. The phase diagrams are revealed for s-wave, p-wave, and d-wave pairings for the SC. A smoking gun experiment is also proposed to confirm the Majorana edge channels which exclude the other possibilities. The heterostructures including TI have been already realized experimentally [55,60]. By this technique, the quantized anomalous Hall effect is realized at higher temperature due to the suppressed inhomogeneity of the exchange gap [55]. Also the different energy position of the Weyl point between the top and bottom surfaces enables the insulating phase which can support the topological magnetoelectric effect

Berry curvature and Chern number
In two-dimensional systems, the total berry curvature of the occupied electron states determines the topological property of the system. It is obtained as [S1] (S1) When the system is gapped, we can define the Chern number as Derivation of the conductance when I > Ic Consider the left half of FIG.3(b), which is a QAHI-TSC-QAHI junction. The electron-electron and electron-hole tunneling probabilities from edge channel j to edge channel i are [S2] T ee Other terms vanish. The current-voltage relation (at zero temperature) is given by [S3] with and other terms vanish. So we have e     Setting I 1 = −I 4 = I and I 3 = I 6 = 0, we obtain (S14) Similar analysis for the right QAHI-TSC-QAHI junction gives (S16) TI surface with general pairing Consider a general pairing potential∆(k) = (Ψ +d · σ)iσ y that acts on a single Dirac coneĥ(k) = g · σ + V z σ z − µ 0 . The constant V z is a Zeeman field along z-direction and g = (k y , −k x , 0) is the spin-orbit coupling vector. In the band basis, the pairing potential∆(k) is transformed tõ where we have defined two angles α and θ as α = arccos V z k 2 + V 2 z , θ = arg(k x + ik y ). (S22) When the Fermi level is not close to the band crossing point, the effect of inter-band pairing (Ψ andd z ) can be ignored. The remaining intra-band pairing is ∆ ± = id y ∓d x = e iθ [±(d x cos θ + d y sin θ) cos α + i(d y cos θ − d x sin θ ∓ Ψ sin α)].
Assuming the Fermi level to be above the band crossing point, then the gap function is just ∆ + . More about the phase diagram with unconventional SC As comparison, we obtained the phase diagrams for both d ⊥ and d ′ , as defined in the main text, with and without d-wave component as shown in FIG.S1. In the case of d ⊥ , the phase diagram is not essentially different from the s-wave case and adding a d-wave component does not change the topological property. For d ′ , the first observation is the inverted sign of the Chern number, which is due to the particular choice of d ′ and the basic topological property is not different from d ⊥ . However, adding a d-wave component to d ′ triplet pairing changes the topology dramatically and four Majorana modes can appear. This change is expected since increasing d-wave component will reverse the sign of the gap, as discussed in the main text.