Platform of chiral Majorana edge modes and its quantum transport phenomena

TL;DR

This paper proposes a heterostructure-based method to realize two-dimensional topological superconductors hosting chiral Majorana edge modes, with potential experimental detection via Josephson junctions and implications for multiple edge modes with unconventional superconductors.

Contribution

It introduces a feasible heterostructure model for 2D topological superconductors with chiral Majorana modes, including experimental detection strategies and the possibility of multiple modes with unconventional superconductors.

Findings

Existence of a topological phase with a single chiral Majorana edge mode in realistic parameters.

Proposed experimental setup using a Josephson junction to detect Majorana modes.

Multiple chiral Majorana edge modes can emerge with unconventional superconductors.

Abstract

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.

Equations

$$H(\mathbf{k})=$$

$$+$$

$${\sigma }_{12}^{>}=\frac{eI}{{\mu }_1-{\mu }_2}=\frac{1}{3},\text{ when }I>I_c.$$

$${\sigma }_{12}^{<}=\frac{1-(I/I_c{)}^2}{2-(I/I_c{)}^2-{\sin }^2{\varphi }_0},\text{ when }I<I_c,$$

$$\gamma =-\text{Im}\sum _{E_n<0}\sum _{n^{\prime }\ne n}\int d^2\mathbf{k}\frac{⟨n|({\partial }_{k_x}{H}_{\mathbf{k}})|n^{\prime }⟩⟨n^{\prime }|({\partial }_{k_y}{H}_{\mathbf{k}})|n⟩-⟨n|({\partial }_{k_y}{H}_{\mathbf{k}})|n^{\prime }⟩⟨n^{\prime }|({\partial }_{k_x}{H}_{\mathbf{k}})|n⟩}{(E_n-E_{n^{\prime }}{)}^2}.$$

$$N_c=\gamma /2\pi .$$

$$T_{31}^{ee}=T_{31}^{he}=T_{61}^{ee}=T_{61}^{he}=1/4,$$

$$T_{34}^{ee}=T_{34}^{he}=T_{64}^{ee}=T_{64}^{he}=1/4,$$

$$T_{16}^{ee}=T_{43}^{ee}=1.$$

$$I_i=\sum _j{K}_{ij}({\mu }_j-{\mu }_s)/e,$$

$$K_{ij}\equiv {\delta }_{ij}-T_{ij}^{ee}+T_{ij}^{he},$$

$$K_{ii}=1,K_{16}=K_{43}=-1$$

$$\displaystyle e\left(\begin{array}[]{cccc}I_{1}\\ I_{6}\\ I_{3}\\ I_{4}\end{array}\right)=\left(\begin{array}[]{cccc}1&-1&&\\ 0&1&&\\ &&1&0\\ &&-1&1\end{array}\right)\left(\begin{array}[]{cccc}\mu_{1}-\mu_{s}\\ \mu_{6}-\mu_{s}\\ \mu_{3}-\mu_{s}\\ \mu_{4}-\mu_{s}\end{array}\right)\rightarrow\left(\begin{array}[]{cccc}\mu_{1}-\mu_{s}\\ \mu_{6}-\mu_{s}\\ \mu_{3}-\mu_{s}\\ \mu_{4}-\mu_{s}\end{array}\right)=e\left(\begin{array}[]{cccc}1&1&&\\ 0&1&&\\ &&1&0\\ &&1&1\end{array}\right)\left(\begin{array}[]{cccc}I_{1}\\ I_{6}\\ I_{3}\\ I_{4}\end{array}\right),$$

$${\mu }_1-{\mu }_s=-({\mu }_4-{\mu }_s)=eI,$$

$${\mu }_6-{\mu }_s={\mu }_3-{\mu }_s=0,$$

$${\mu }_3={\mu }_6={\mu }_s=({\mu }_1+{\mu }_4)/2,$$

$${\mu }_1-{\mu }_6={\mu }_3-{\mu }_4=eI.$$

$${\mu }_4={\mu }_5={\mu }_s^{\prime }=({\mu }_2+{\mu }_3)/2,$$

$${\mu }_5-{\mu }_2={\mu }_3-{\mu }_4=eI.$$

$$\tilde{\Delta }(\mathbf{k})=(\tilde{\Psi }+\tilde{\mathbf{d}}\cdot \mathbf{\sigma })i{\sigma }_y$$

$$\tilde{\Psi }=-i{d}_z{e}^{i\theta },$$

$${\tilde{d}}_x=[-\cos \alpha (d_x\cos \theta +d_y\sin \theta )+i\Psi \sin \alpha ]e^{i\theta },$$

$${\tilde{d}}_y=(d_y\cos \theta -d_x\sin \theta )e^{i\theta },$$

$${\tilde{d}}_z=[\sin \alpha (d_x\cos \theta +d_y\sin \theta )+i\Psi \cos \alpha ]e^{i\theta },$$

$$\alpha =\arccos \frac{V_z}{\sqrt{k^2+V_z^2}},\theta =\arg (k_x+i{k}_y).$$

$${\Delta }_{\pm }=i{\tilde{d}}_y\mp {\tilde{d}}_x=e^{i\theta }[\pm (d_x\cos \theta +d_y\sin \theta )\cos \alpha +i(d_y\cos \theta -d_x\sin \theta \mp \Psi \sin \alpha )].$$

Full text

Introduction —

Majorana fermions in condensed matter systems Wilczek ; Kitaev2001 , or Majorana quasi-particles, have been the subject of intense research due to their exotic properties such as fractional Josephson effect Kitaev2001 ; Tanaka2 ; Tanaka3 ; Kwon ; Fu1 ; Lutchyn ; Vic2011 , resonant Andreev reflection Vic2009 ; Wimmer ; Tanaka1 ; Bolech , enhanced cross Andreev reflection JamesNC , spin selection Simon ; JamesPRL ; Oreg ; XLiu ; Lunhui ; Noah , non-Abelian statistics RG ; Ivanov ; Fujimoto ; STF ; Alicea2011 , etc. Tremendous theoretical Fu2008 ; Sau2010 ; Lutchyn ; Alicea2010 ; ORV ; Potter ; TP ; Sho ; Vic2018 ; Ruixing and experimental Mourik ; Deng ; Das ; Yazdani2014 ; JiaPRL ; Marcus1 ; Marcus2 ; Yazdani2017 ; Kouwenhoven2018 progress has been achieved to create zero-dimensional Majornana bound states in one-dimensional topological superconductors (1D TSCs). Promising braiding methods have been proposed based on such 1D systems AliceaNP for the final proof of the non-Abelian particles, as well as for the creation of topological quantum computers Kitaev2003 ; Nayak .

On the other hand, 1D propagating Majorana modes at the edges of 2D TSCs have been realized recently in an experiment Qinglin2017 that combines quantum anomalous Hall insulators (QAHIs) with superconductors (SCs) Qi2010 . 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 $\Delta$. 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 Chung2011 ; JWang2015 ; CZ2017 .

However, controversy arises based on an observation that half-quantized conductance may also appear trivially without any Majorana mode Wen2018 ; Sau2018 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 the superconductivity order parameter, resulting in the requirement $M_{z}<\Delta$ Qi2010 . 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\leqslant\sigma_{12}\leqslant 1$ to $\sigma_{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,\uparrow}(\bm{k}),c_{t,\downarrow}(\bm{k}),$ $c_{b,\uparrow}(\bm{k}),c_{b,\downarrow}(\bm{k}),$ $c_{t,\uparrow}^{\dagger}(-\bm{k}),c_{t,\downarrow}^{\dagger}(-\bm{k}),$ $c_{b,\uparrow}^{\dagger}(-\bm{k}),c_{b,\downarrow}^{\dagger}(-\bm{k})\}$, the effective Hamiltonian of the system is

[TABLE]

The Pauli matrices $s_{x,y,z},\sigma_{x,y,z},\tau_{x,y,z}$ act on spin, layer and particle-hole spaces respectively. $M_{z}$ is the exchange field felt by the bottom layer, $\Delta(\bm{k})$ and $\bm{d}(\bm{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. $\mu$ is the chemical potential while $\delta 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 Schnyder ; Teo 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 Mogi . Note that it becomes a QAHI if the SC on the top surface is replaced by another FM same as the bottom surface Yu . 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 Fu2009 . 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 $\gamma$ is shown in FIG.2(a). In normal states, $\gamma$ corresponds to Hall conductance. If a bulk gap exists, it is quantized so that $N_{c}=\gamma/2\pi$ with the integer $N_{c}$ being the Chern number TKNN . 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 $|\mu|<|M_{z}|$, the system is gapped and $N_{c}=\text{sign}[M_{z}]$, corresponding to a single chiral Majorana edge mode whose chirality is determined by the magnetization direction. If $|\mu|>|M_{z}|$, the bottom surface becomes gapless and the total Berry curvature is not quantized.

When $t_{c}\neq 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 $\mu=0$, as shown in FIG.2(b). If $M_{z}>0$ is small and $\mu$ deviates from zero, the system first enters a non-trivial region with $N_{c}=\text{sign}[M_{z}]$ and then transitions into a trivial phase again for large $|\mu|$, 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 $\mu=0$, giving $N_{c}=0$. As $|\mu|$ increases, the Fermi level cuts a single band and the resulting SC is topological with single chiral Majorana edge mode. As $|\mu|$ 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 Lutchyn ; Potter . The similarity becomes clearer in band structure when a difference in chemical potentials of the top and bottom surfaces, $\delta 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 $\mu<0$, as shown by the dashed line in FIG.2(c) where the whole phase diagram is obtained. The effect of $\delta E$ has been discussed Morimoto and it happens in real experiments Yoshimi . 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 Chung2011 ; JWang2015 ; Qinglin2017 . To avoid trivial explanations Wen2018 ; Sau2018 and uniquely determine the existence of chiral Majorana modes, however, we propose an experimental setup including a Josephson junction CZ2018 ; Shen2018 . 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 Yu 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 Shen2018 . 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_{c}^{bulk}$ ($I_{c}^{bulk}$ 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 $\mu_{i}$ ($i=1,2,...,6$) for each QAHI edge state as shown in FIG.3(b). They must satisfy $\mu_{1}-\mu_{6}=\mu_{3}-\mu_{4}=\mu_{5}-\mu_{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 Chung2011 . A simple application of the scattering coefficients to multi-terminal measurement leads to $\mu_{3}=\mu_{6}=(\mu_{1}+\mu_{4})/2$. Similarly, we get $\mu_{4}=\mu_{5}=(\mu_{2}+\mu_{3})/2$ for the right half of the system. Combination of these relations leads to

[TABLE]

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 $\mu_{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 $\phi$ which is related to the current by $I=I_{c}\sin(\phi/2)$ or $\phi=2\arcsin(I/I_{c})$. The existence of $\phi$ can drastically affect $\sigma_{12}$ by inducing interference between different Majorana paths and thus changing the tunneling amplitude Shen2018 . Taking this into account, one obtains the conductance

[TABLE]

where $\phi_{0}=k_{F}L$ is the kinetic phase acquired by the edge states across the junction Shen2018 . ($k_{F}$ is the Fermi wave vector and $L$ is the length of the junction.) When $I\ll I_{c}$, $\phi$ is small and its effect on $\sigma_{12}$ is negligible. As $I$ approaches $I_{c}$, $\phi$ increases from zero to $\pi$ and $\sigma_{12}$ decreases to zero, as shown in FIG.3(c). Note that in the special case with $\phi_{0}=\pi/2+n\pi$, $\sigma_{12}^{<}$ is constantly unity.

In summary, as the current increases, the two-terminal conductance $\sigma_{12}$ starts with a value between $1/2$ and $1$ and decreases until it exceeds the Josephson critical current, above which $\sigma_{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 Linder . 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 Yanase1 ; Yanase2 ; Yada .

Consider a general pairing potential $\hat{\Delta}(\bm{k})=(\Psi+\hat{\bm{d}}\cdot\bm{\sigma})i\sigma_{y}$ that acts on TI surface states with the Hamiltonian $\hat{h}(\bm{k})=\bm{g}\cdot\bm{\sigma}+V_{z}\sigma_{z}-\mu_{0}$ where $V_{z}$ is a perpendicular Zeeman field and $\bm{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 ($\mu_{0}>0$), then the superconductivity gap function in the band basis becomes $\Delta_{+}=e^{i\theta}[(d_{x}\cos\theta+d_{y}\sin\theta)\sin\alpha+i(d_{y}\cos\theta-d_{x}\sin\theta-\Psi\cos\alpha)]$ with $\alpha=\arcsin(V_{z}/\sqrt{k^{2}+V_{z}^{2}})$ and $\theta=\arg(k_{x}+ik_{y})$.

If only singlet (s-wave and d-wave) pairings are considered, the gap function on the top surface state is simply $\Psi\cos\alpha=(\Delta_{s}+\Delta_{d}\cos 2\theta)\cos\alpha$. When $|\Delta_{s}|>|\Delta_{d}|$, it is fully gapped and the topological property is the same as the pure s-wave case discussed previously. When $|\Delta_{s}|<|\Delta_{d}|$, it becomes nodal and there may be edge states depending on the direction of the open boundaries, similar to well-known usual d-wave superconductors Hu ; Ryu ; Tanaka2 ; Tanaka5 ; Sato ; FaWang ; Noah2 . For the heterostructure in FIG.1, the energy dispersion on the [010] and [110] edges are shown in FIG.4 for pure d-wave. (The s-wave component only shifts the positions of the nodal points.) If $M_{z}=0$, this system behaves just as usual d-wave SCs in which flat bands appear in [110] direction (FIG.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, $\bm{d}_{\parallel}=\Delta_{p}(\sin\theta,-\cos\theta,0)$, $\bm{d}_{\perp}=\Delta_{p}(\cos\theta,\sin\theta,0)$ and $\bm{d}^{\prime}=(\sin\theta,\cos\theta,0)$. In the first two cases, the d-vectors are either parallel or perpendicular to the spin-orbit vector $\bm{g}$. For $\bm{d}_{\parallel}$, the gap function is $(\Delta_{p}+\Psi\cos\alpha)$ so that p-wave behaves the same as s-wave. For $\bm{d}_{\perp}$, the gap function becomes $\Delta_{p}\sin\alpha-i\Psi\cos\alpha$, which is always fully gapped when $\Delta_{p}V_{z}\neq 0$. In the heterostructure, it can support single chiral Majorana mode. With the third choice, $\bm{d}^{\prime}$, the gap function is $\Delta_{p}\sin 2\theta\sin\alpha+i(\Delta_{p}-\Delta_{d}\cos\alpha)\cos 2\theta$. Note that, when $\Delta_{p}V_{z}\neq 0$, this is also fully gapped except at $\Delta_{p}=\Delta_{d}\cos\alpha$ where the gap changes sign. A phase diagram with this pairing potential ($\Delta_{p}=0.3,\Delta_{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 Mogi ; Yoshimi . By this technique, the quantized anomalous Hall effect is realized at higher temperature due to the suppressed inhomogeneity of the exchange gap Mogi . 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 Yoshimi . With these artificial structures, one can design various Majorana edge channels to realize the circuits with dissipationless current and even quantum computation SCZhangPNAS .

Acknowledgements.

Acknowledgment —

J.J.H. is very grateful to Chao-Xing Liu for discussions. N.N. was supported by Ministry of Education, Culture, Sports, Science, and Technology Nos. JP24224009 and JP26103006, the Impulsing Paradigm Change through Disruptive Technologies Program of Council for Science, Technology and Innovation (Cabinet Office, Government of Japan), and Core Research for Evolutionary Science and Technology (CREST) No. JPMJCR16F1 and No. JPMJCR1874, Japan. Y.T. was supported by Grant-in-Aid for Scientic Research on Innovative Areas, Topological Material Science (Grants No. No. JP15H05851, No. JP15H05853, and No. JP15K21717) and Grant-in-Aid for Scientic Research B (Grant No. JP18H01176) from the Ministry of Education, Culture, Sports, Science, and Technology, Japan (MEXT).

.1 **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 SupBernevigBook

[TABLE]

When the system is gapped, we can define the Chern number as

[TABLE]

.2 **Derivation of the conductance when $I>I_{c}$

**

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 SupChung2011

[TABLE]

Other terms vanish. The current-voltage relation (at zero temperature) is given by SupDatta

[TABLE]

with

[TABLE]

and other terms vanish. So we have

[TABLE]

Setting $I_{1}=-I_{4}=I$ and $I_{3}=I_{6}=0$, we obtain

[TABLE]

so that

[TABLE]

Similar analysis for the right QAHI-TSC-QAHI junction gives

[TABLE]

.3 **TI surface with general pairing

**

Consider a general pairing potential $\hat{\Delta}(\bm{k})=(\Psi+\hat{\bm{d}}\cdot\bm{\sigma})i\sigma_{y}$ that acts on a single Dirac cone $\hat{h}(\bm{k})=\bm{g}\cdot\bm{\sigma}+V_{z}\sigma_{z}-\mu_{0}$. The constant $V_{z}$ is a Zeeman field along z-direction and $\bm{g}=(k_{y},-k_{x},0)$ is the spin-orbit coupling vector. In the band basis, the pairing potential $\hat{\Delta}(\bm{k})$ is transformed to

[TABLE]

with

[TABLE]

where we have defined two angles $\alpha$ and $\theta$ as

[TABLE]

When the Fermi level is not close to the band crossing point, the effect of inter-band pairing ($\tilde{\Psi}$ and $\tilde{d}_{z}$) can be ignored. The remaining intra-band pairing is

[TABLE]

Assuming the Fermi level to be above the band crossing point, then the gap function is just $\Delta_{+}$.

.4 **More about the phase diagram with unconventional SC

**

As comparison, we obtained the phase diagrams for both $\bm{d}_{\perp}$ and $\bm{d}^{\prime}$, as defined in the main text, with and without d-wave component as shown in FIG.S1. In the case of $\bm{d}_{\perp}$, 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 $\bm{d}^{\prime}$, the first observation is the inverted sign of the Chern number, which is due to the particular choice of $\bm{d}^{\prime}$ and the basic topological property is not different from $\bm{d}_{\perp}$. However, adding a d-wave component to $\bm{d}^{\prime}$ 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.