Detection of Majorana zero mode in the vortex

Majorana zero modes, which behave like Majorana fermions, are quasiparticle excitations in condensed matter systems. They obey non-Abelian statistics, and have been proposed as building blocks of topological quantum computers. They are predicted to exist in the vortex of topological superconductors. In 2012, such a topological superconductor was engineered by depositing topological insulator thin films on top of an s-wave superconductor. Thereafter, several evidences have been reported to prove the Majorana zero modes’ existence in the vortex. In this review, by putting all experimental and theoretical results together, we show that these experimental evidences are consistent and they are also strongly supported by the theories, so the existence of Majorana zero mode is firmly established. Moreover, the adjacent Majorana zero modes annihilate when two vortices are close enough, which demonstrate that they have the nature of Majorana fermions. Finally, their potential application in topological quantum computing is discussed.


INTRODUCTION
Majorana fermions (MFs), the elusive particles that are their own antiparticle, were firstly predicted by Ettore Majorana. 1 Now, their concept is finding ever wider relevance in modern physics, not only in neutrino physics, supersymmetry, and dark matter, but also on some exotic states of condensed matter. 2,3 The recent effort in searching such a genuine particles focuses on the neutrinoless double beta decay experiment. 4 Some quasiparticles, which behave like MFs, may emerge as a novel excitation in some condensed matter systems. These Majorana quasiparticles obey non-Abelian statistics and may be used as robust building blocks in quantum computing. 5 Due to the scientific importance and potential applications in quantum computing, Majorana quasiparticles have attracted lots of attentions recently.
Theoretically, Majorana quasiparticles are predicted to emerge in the vortex core of topological superconductors (TSC) as a zero energy excitation mode called Majorana zero mode (MZM). Unfortunately, there are at present no known materials with this property. In 2008, the search for MZMs took a big step forward when Fu and Kane proposed an ingenious idea to artificially engineer a topological superconductor based on a normal superconductor and the exotic properties of the electronic surface states of a topological insulator (TI). 6 After then, a race for searching MZMs started all over the world. [7][8][9][10][11] By superconducting proximity effect, 12 superconductivity can be introduced into TI. However, because of the interface interactions, the huge lattice mismatches and the instability of TIs, the interface of a superconductor (SC) and a TI is very difficult to be established. In 2012, our group successfully fabricated the high quality single crystal thin films of TI Bi 2 Se 3 on a conventional s-wave superconductor NbSe 2 by means of molecular beam epitaxy (MBE) and demonstrated the co-existence of topological surface states (TSSs) and superconductivity in such system. 13 By putting TI on the top of SC, this system has several significant advantages that promote the research on MZMs in this platform. First, MBE with the help of reflection high-energy electron diffraction, allows us to control the thickness of TI films as precise as one layer during the growing process. As a result, we can grow heterostructures with different thicknesses and study their superconductivity related electronic properties in detail, since the topological properties and proximity-induced superconductivity are directly affected by the TI thin film thickness. And more importantly, if the TI film is thin enough to make the up surface states superconducting, the resulted TSC can be studied by surface sensitive techniques directly. Second, state of the art MBE method has a reputation of growing thin films with extremely high quality. TI films with better quality compared with bulk samples will have Fermi level closer to the Dirac point. Also, Fermi level of the TI film can be tuned by doping. The Fermi level position plays a crucial role in searching MZM which will be discussed later. Third, in other systems, it is required to apply a sufficiently large magnetic field to make the system topological in the first place, in which case the bands are already spin polarized. It will then be difficult to detect the spin property of the MZM since both MZM (if present) and finite-energy midgap states would be spin polarized. In TI/SC system, a small magnetic field is sufficient to induce vortices and, hence, the MZM in the vortex core. Small magnetic field also makes this system more promising for further application.
This article mainly focuses on the recent experimental works on the 3D TI/SC heterostructure system, which is also referred as an artificial TSC. We review the structure and fabrication of the heterostructure in Section 'Topological insulator-superconductor heterostructure', and its superconducting property in Section 'Topological superconductor'. Then we discuss the detection of Majorana Zero Mode in such system in Section 'Experimental detection of Majorana zero mode' and finally give a brief summary and outlook in Section 'Summary and outlook'.

TOPOLOGICAL INSULATOR-SUPERCONDUCTOR HETEROSTRUCTURE
The original proposal of a TI-SC heterostructure by Fu and Kane requires the SC on top of the TI. 6 However the first reported heterostructure chose an opposite scenario. The well-known 3D TI Bi 2 Se 3 was grown on top of an s-wave SC NbSe 2 via MBE. 14 Soon another similar 3D TI Bi 2 Te 3 was also grown on the same substrate. 15,16 The successful fabrication of such heterostructure built a promising platform for the detection of MZMs. Here we review the material properties and growth procedures of the heterostructure and discuss some details of the heterostructure that are significant in the detection of the MZMs.
The SC substrate 2H-NbSe 2 is a transition-metal dichalcogenide material, a layered material that has very weak interaction between two adjacent Se-Nb-Se sheets. The in-plane lattice constant of NbSe 2 is 3.44 Å. In bulk NbSe 2 , a charge density wave phase emerges below 33 K and a superconducting phase emerges below 7.2 K. Bulk NbSe 2 crystals are prepared via chemical vapor transport method. The size of the flakes is typically larger than two by 2 mm 2 , suitable as a substrate for the further MBE growth and scanning tunneling microscopy (STM) measurements.
The 3D TIs grown on top of NbSe 2 are also layered materials Bi 2 Se 3 and Bi 2 Te 3 . They have a five-atom-layer sheet structure Se (Te)-Bi-Se (Te)-Bi-Se (Te). One five-atom-layer sheet is called a quintuple layer (QL). The in-plane lattice constants of Bi 2 Se 3 and Bi 2 Te 3 are 4.13 and 4.38 Å, respectively. The most significant properties of 3D TIs are their TSSs. [17][18][19][20][21] The TSSs lie in the bulk band gap of TIs and connect the valance band and the conduction band. Ideally, Fermi levels of the TIs lie inside the bulk band gap and cut through the TSSs only. However the Fermi levels of the grown bulk materials Bi 2 Se 3 and Bi 2 Te 3 always lie inside the conduction band while Fermi level of Sb 2 Te 3 , another 3D TI with similar structure, lies inside the valance band. This is due to the large amounts of defects inside the bulk materials.
The 3D TIs in the heterostructure are grown via MBE technique, a delicate thin film growing method that provide atomic layer thickness control ability and extreme high quality films with few defects. These two features of MBE method are crucial in the growth of TI/SC heterostructure: First, as mentioned before, defects in the samples will move the Fermi level away from the Dirac points of the TSSs towards conduction band of valance band. Electrons that form superconducting gap around Fermi level will then come from the bulk bands rather than from the TSSs. However, only when superconductivity is induced into the TSSs from a SC substrate via proximity effect, Cooper pairs in the TSSs will behave "spinless" and MZMs will be hosted inside an Abrikosov vortex on the TI surface. The high quality thin films grown by MBE can overcome this problem by making the Fermi level much closer to the Dirac point. Bi 2 Se 3 and Bi 2 Te 3 have been grown on several substrates. [22][23][24][25] The growth conditions of the two materials are similar. The substrate temperature is controlled to be lower than the evaporation temperature of Bi and higher than the evaporation temperature of Se or Te for the two kinds of TIs, which is typically range from 500 to 570 K. Under this temperature, Se or Te molecules will desorb from the substrate surface except that those form Bi 2 Se 3 and Bi 2 Te 3 with Bi atoms. In this case, the growing process can be performed under a Se (Te)-rich environment, where flux ratio of Se (Te) and Bi can be set to over 10. As a result, the extra amount of Se (Te) will suppress the formation of Se (Te) vacancies and antisite defects in the MBE grown films and the Fermi level of the thin films will be closer to the Dirac point.
Second, proximity effect that transfers Cooper pairs from the SC to the TI TSSs decays exponentially with the thickness of the TI film above. Once the TI films are too thick, the superconducting gap on its surface will be too small for the MZMs measurements via STM with current energy resolution (around 0.1 meV), since the excitations lie exactly at zero energy inside the superconducting gap. On the other hand, TSSs of 3D TIs is fully developed after a certain layers of the materials. If the films are too thin, i.e., less than 2 QLs on Si substrate, interactions between the top and bottom surfaces will destroy the TSSs. 23 Bi 2 Se 3 is grown on NbSe 2 at substrate temperature 500 K. It is grown in a layer-by-layer mode. 13 Bulk Bi 2 Se 3 crystal has a layer thickness of 0.94 nm while the first layer Bi 2 Se 3 on NbSe 2 is as thick as 1.5 nm. The larger thickness of the first QL suggests an interfacial layer exist between the first QL Bi 2 Se 3 and the NbSe 2 substrate. A detailed study of this interfacial layer reveals its structure in ref.
14. An NaCl-type BiSe structure forms on the NbSe 2 substrate before the formation of Bi 2 Se 3 to compensate the lattice mismatch between the two.
The growth of Bi 2 Te 3 on NbSe 2 is very similar with Bi 2 Se 3 on NbSe 2 . The structure of Bi 2 Te 3 on NbSe 2 is shown in Fig. 1. Combining MBE with RHEED technique, it is possible to control the exact QL number of the deposited thin films by monitoring RHEED oscillation pattern during the sample growth. Each intensity maximum in the RHEED oscillation pattern represents the formation of a complete QL of Bi 2 Te 3 . It is then convenient to control the thickness of the Bi 2 Te 3 on top of SC substrate precisely to obtain a heterostructure with a measureable superconducting gap and well developed TSSs. Typically 5 QL or 6 QL samples are suitable for the detection of MZMs. At such thickness, the Fermi energy just cut through the conduction band bottom of Bi 2 Te 3 , as shown in Fig. 1d, which provides enough contribution of the superconductivity from the TSSs.

TOPOLOGICAL SUPERCONDUCTOR
MZMs exist in the vortex of a spinless p x + ip y superconductor, [26][27][28][29] which is called a TSC. It is proposed that via the proximity effect from an s-wave SC, the two dimensional TSSs in a TI resemble a spinless p x + ip y like SC and thus host MZMs in vortices. 6 The purpose to build a TI/SC heterostructure is to induce superconductivity in the TI, or specifically in the TSSs, via proximity effect. In this section, we discuss the superconductivity in the TI/SC heterostructure. Superconductivity in the TI/SC heterostructure is analyzed via both scanning tunneling spectroscopy (STS) measurements and high energy resolution angle resolved photoemission spectroscopy (ARPES) measurements, 30 as shown in Fig. 2. The STS spectra clearly demonstrate gap features on the surface of the TI above SC. These gaps can be eliminated by raising temperature and applying magnetic field. Moreover, with increasing layer thickness of the TI, the size of the gap decays exponentially, which is a typical behavior of the proximity induced superconductivity. Undoubtfully, the spectra gaps on the surface of the TI are superconducting gaps. Since STS provides no momentum resolution, whether TSSs contribute to the gaps is unknown. ARPES measurements explicitly answer the question. From a series of measurements at k-space of the TSSs, the density of states clearly show gaps structures near the Fermi level. The spinresolved ARPES experiment reveals that superconductivity is induced into the TSSs of the TI and two dimensional topological superconductivity in a helical TSSs is distinctly different from that in an ordinary two-dimensional superconductor in terms of the spin degrees of freedom of electrons. These results from STS and APRES demonstrate that a TI/SC heterostructure successfully creates superconductivity in the TSSs which become a spinless px + ipy like TSC.
Applying magnetic field to a type-II superconductor will cause the emergence of Abrikosov vortices. 31 As mentioned before, a spinless p x + ip y SC will host MZMs inside these vortices excited by magnetic field. Those vortices can be imaged in real space via STS mapping. Low energy bound states caused by Andreev reflection at the center of the vortices will have a high density of state inside the superconducting gap while the surroundings have zero density of state at zero temperature and a low density of state at finite temperature. Thusly, positions of the vortices are revealed by the zero bias conductance map. As shown in Fig. 3, vortices represented by purple circles form hexagonal lattices on the surface of the SC substrate, as well as the surface of the TI in the TI/SC heterostructure. The lattice constant of the vortices, or the area owned by each vortex is determined by the magnitude of the field. Multiplying this area to the magnetic field, one should get the flux where A is the area of one vortex, as shown in the red diamond in Fig. 3a and e, and B is the strength of the field. As can be seen from Fig. 3b and c, decreasing the magnetic field will increase the distance between vortices, which has an interesting effect on the detection of the MZMs inside the vortices and will be discussed in the next section.
From the zero bias conductance map, we can get the coherence length of the SC in the Ginzburg-Landau theory, where σ is the zero bias conductance,σ 0 is the zero bias conductance at the vortex center, r is the distance from vortex center and ξ is the in-plane coherence length. 32 As shown in Fig. 3f, the in-plane coherence length of the heterostructure increases with the thickness of the TI thin film above. This is in consistence with the lowering of the upper critical field of the TI thin film H c2~2 .8 T in 3 QL Bi 2 Te 3 /NbSe 2 , compared with H c2 3.4 T in NbSe 2 , since H c2 = ϕ 0 /2πξ 2 . It is worth noting that the coherence length on the surface of the TI/SC heterostructure decreases with the increasing magnetic field, as shown in Fig. 3h. This phenomenon suggests that the superconductivity on the TI surface is not an ordinary single band s-wave superconductivity, since the coherence length of an NbSe 2 remains the same under different fields. 33,34 In the TI/SC heterostructure, superconductivity is induced into the TSSs via proximity effect. The above mentioned experimental evidence indicates that this induced superconductivity is the topological superconductivity that is different form the s-wave superconductivity of the substrate. Therefore, this heterostructure is truly a full gap artificial TSC.

EXPERIMENTAL DETECTION OF MZM
The fabrication of the artificial TSC was a big step towards the detection of the MZMs. Many theoretical and experimental works that focused on this platform gave a series of significant results on the existence of MZMs. In this section, we will first review some of the theoretical works of MZMs in the artificial TSC. Then we will discuss the experimental detection of MZMs in the artificial TSC which are strongly supported by those theories.
The features of MZM predicted by theory In the early theoretical work, the artificial TSC was designed to have a thin layer of SC on top of TI, 6 which is experimentally difficult to achieve. Growing TI on SC was proved to be much easier. 14,15 In this artificial TSC, two well-separated MZMs exist in a vortex core on both top and bottom surface of TI thin films. 35 As shown in Fig. 4a-d, the wave function distribution of the lowest energy bound state in the vortex evolves with the layer thickness n of the TI above the SC substrate. When n = 2, as shown in Fig. 4a, the bound state is connected across the top and bottom surface of the TI. In Fig. 4b where n = 4, the distribution of the state on the top and bottom surface is disconnected and localized, which represents the signature of two MZMs. For larger n in Fig. 4c and d, the coupling between the two surface layers is even weaker.  As a result, two MZMs become more disconnected in the z direction and more extended in the surfaces. It is worth noting that as the MZMs get wider, they become more fragile. Another effect caused by the increase of TI's thickness is the decay of the pairing gap induced via proximity effect. It is therefore important to control the TI's thickness of the heterostructure in order to get fully developed, as well as robust MZMs in the vortices of the surface. Since STM has both real space and energy space resolution, it is the best tool to detect such a MZM inside a vortex core on the TI's top surface. However there exist other low energy quasiparticle states in the vortex core that separated from the MZM by tiny energy gap, i.e., the Caroli-de Gennes-Matricon states. [36][37][38][39] The energy gap is estimated to be 0: , where Δ is the superconducting gap and E D is the Fermi energy relative to the Dirac point of the TSS band. 40 In this TI/SC system, the typical energy gap size is Δ∼1 meV and the Fermi level is E D ∼100 meV, which leads to the minigap with size ∼0.01 meV. This is much smaller than the present energy resolution (0.1 meV) in STS. There are two ways to increase the mini gap. One is to tune the Fermi level toward the Dirac point. However, in that case, the superconducting gap Δ becomes very small and the transition temperature becomes very low due to the weaker proximity effect. The other is using another SC material with larger gap, which is yet to be found. Hence, the direct observation of the Majorana mode is still difficult.
Fortunately, the MZM is not a point at the center of a vortex core, but extensively distributes around the core center, 41 which gives an opportunity to detect the MZM by investigating the spatial distribution of the bound states in the vortex core. This spatial distribution of MZM is different from other low energy bound states inside the vortex. Shifting the Fermi energy level from TSSs towards the bulk bands, the spatial distribution of the core state will change from a Y-shape when MZM exists to a Vshape when only other low energy excitation states exist in STS measurements with finite resolution 42 (see Fig. 4e-h).
Moreover, the MZM inside a vortex undergoes unique spinselective Andreev Reflection (SSAR) process. This process was first proposed in a semiconductor wire/SC structure in which MZMs exist at two ends of the wire. 43 The proposed SSAR of the MZM in 1D wire can be generalized to a 2D TSC. In the latter case, the spin of the MZM in the vortex core has a spatial distribution and is not conserved. However, the spin wave function at the center of the MZM at the vortex core is fully polarized along the external magnetic field. Therefore, the Andreev reflection at the vortex center is expected to be spin-selective. [42][43][44] Via an spin-polarized STM, 45,46 we can inject spin-polarized tunneling current in STS measurements and observe the zero bias conductance change caused by the SSAR from the MZM when flipping the spin of the tunneling electrons. 47,48 It is worth noting that only a small magnetic field is sufficient to excite vortices on the TI surface and hence the spin-polarized MZM in the vortex core. On the other hand, the spin polarization of other surface states and finite energy quasiparticle states is still negligibly small. This is a unique advantage to probe the SSAR in the TSC. In other systems, i.e., in Rashba wire systems, a sufficiently large magnetic field is required to lift the spin degeneracy and make the system topological in the first place, in which case the bands are already spin polarized. It will then be difficult to attribute spin-dependent zero-bias conductance to the MZM since both MZM (if present) and finiteenergy midgap states would be spin polarized.
Zero bias peak As mentioned above, energy resolution in most STS measurements is not sufficient to distinguish the MZM from other low energy bound states in the vortex core. All these excitations will merge together and form a single peak centered at zero energy with certain broadening width in an STS measurement. This peak is referred as the zero bias peak (ZBP). Though we cannot separate the MZM from other excitations, the conductance contributed from the MZM at zero energy reveals MZM's existence. In the artificial TSC formed by Bi 2 Te 3 /NbSe 2 heterostructure, a series of STS measurements focused on the ZBP at the center of the vortex on the TI's surface were performed and the results are given in Fig. 5. STS spectra taken at the vortex core under magnetic fields ranging from 0.025 T to 1.00 T on three different samples, 5 QL Bi 2 Te 3 /NbSe 2 , 2 QL Bi 2 Te 3 /NbSe 2 and bare NbSe 2 , are shown in Fig. 5a-c respectively. In Fig. 5a, the spectra taken on 5 QL Bi 2 Te 3 / NbSe 2 under magnetic field ranging from 0.25 T to 0.10 T are obviously different from others. The intensities of the ZBP in these spectra are dramatically higher than the others. In Fig. 5b, such ZBP intensity anomaly is absent on the 2 QL Bi 2 Te 3 /NbSe 2 . In Fig. 5c, the ZBP intensity anomaly is again absent on the NbSe 2 a b c substrate. This fact shows the ZBP intensity anomaly cannot be observed on normal s-wave superconductors. This ZBP anomalous behavior is therefore believed to come from the contribution of the MZM in the vortex core of the 5 QL Bi 2 Te 3 /NbSe 2 sample. In the 2 QL Bi 2 Te 3 /NbSe 2 , the Fermi level lies deep in the bulk bands of the TI and does not cut through the TSSs (see Fig. 1). Thusly, the bulk carriers dominate in the ZBPs, so they behavior just like those in s-wave superconductor. In this sample, the absence of the ZBP intensity anomaly also shows that the anomaly is not due to the proximity effect. In Fig. 5a, the ZBP anomaly only shows up at magnetic field lower than 0.18 T. At higher field, the intensities of the ZBPs are comparable to those in 2 QL samples. This phenomenon can be explained by the distance between two adjacent vortices and the spatial distribution of MZM's wavefunction. As mentioned before, two adjacent vortices will come closer when the excitation magnetic field is increased (see Fig. 3).
Meanwhile, MZM has a spatial distribution centered at vortex core (see Fig. 4a-d). The MZM spatial distribution is around 100 nm, which is close to the distance between two adjacent vortex at 0.18 T. When the field is large enough, the two MZMs at two adjacent vortices cores overlap and annihilate each other. 16 More evidence of this argument is discussed in the following experiment that focus of the spatial distribution of MZMs (see Figs. 6 and 7) and the SSAR effect (see ref. 48).

Spatial distribution
Besides from the ZBP anomaly, the spatial distribution of MZMs in the artificial TSC were also detected experimentally. The spatial distribution of a MZM has a wide extension centered at vortex core, with its energy being exactly at zero. The spatial distribution of other low energy quasiparticle excitations have energy dispersions correlated between their spatial positions of the maximal absolute value of their wave functions. As a general rule, a quasiparticle with higher energy locates at a position further away from the center of vortex (see Fig. 6a). In the Bi 2 Te 3 /NbSe 2 heterostructure, STS spectra were taken along a line starting from the center of the vortex to the edge of the vortex on samples with various layer thickness of TI above SC. Two kinds of spatial distributions of the quasiparticle excitations in the vortex are observed via STS, 16 as shown in Fig. 6. In samples from 1 QL to 3 QL layers of TI, the splitting of the excitations starts from the vortex core, and extend linearly towards the vortex edge, forming a V-shape pattern (see Fig. 6b-d). In samples from 4QL to 6 QL layers of TI, the splitting starting point shifts away from the vortex center, forming a Y-shape pattern. With increasing layer thickness, the splitting point is further away. The Y-shape splitting pattern originates from the MZM at the center of the vortex. The wide spatial distribution of MZM smears the splitting of other low energy excitations and causes a large ZBP around the vortex core. As shown in Fig. 4a-d, the spatial distribution of the MZM become wider when the TI become thicker. This explains the why splitting points shifts outside with the increasing TI layer thickness. More importantly, we found that when the ZBP anomaly appears, the Y-shape splitting pattern appears, and they also disappear together. As shown in Fig. 7a, b, d and e, the Y-shape splitting pattern can never be observed on s-wave superconductor NbSe 2 , or 1 QL Bi 2 Te 3 /NbSe 2 . It can only be observed on thick TI films, like 5QL Bi 2 Te 3 /NbSe 2 and at low magnetic field (below 0.18 T, Fig. 7c). When the magnetic field is larger than 0.18 T, the Yshape splitting pattern chances back to the V-shape splitting pattern (Fig. 7f). The difference between low field and heigh field which only exists in the sample with MZMs strongly suggests that this difference is caused by the destroying of the MZMs field due to the interaction between them when the field gets larger and the vortices get closer.
If we cut along zero bias in Fig. 7c and f, two curves are obtained. The one from Fig. 7c contains the contribution from MZM, and the one from Fig. 7f does not. So, we can extract the contribution of MZM alone by subtracting them. Those curves and their subtracted result are shown in Fig. 7g, with rotation along zaxis, from which the spatial distribution of a MZM can be deduced. The cone-like shape agrees quite well with the calculated one (Fig. 7h)   Spin-selective Andreev reflection Via a spin-polarized STM, we can observe the SSAR on the MZM in the artificial TSC. When an electron with the same spin as the MZM's inside is injected into the vortex core, SSAR will cause the reflection of a hole with the same spin. On the other hand, when the injected electron has the spin opposite to the MZM's spin, the SSAR will cause the reflection of the electron directly. The difference of the reflected particles in the SSAR process leads to a difference in the tunneling conductance in the STS. Injecting an electron and reflecting a hole will give a high tunneling conductance while injecting an electron and reflecting an electron as well will give a low tunneling conductance. The total tunneling conductance during the process consists of two parts, where the first term is the contribution from the normal tunneling, and the second term from the SSAR. E is the energy, r is the distance from the vortex core,BandMare the orientations of the external magnetic field B and spin polarization M, respectively. It is worth noting that the SSAR is most significant at the vortex center, where the superconducting order parameter vanishes and the SSAR is purely contributed from the MZM. The term dI/dV| n is proportional to the local density of states caused by normal tunneling, which is independent of spin polarization at r = 0 and E = 0 within the energy resolution of about 0.1 meV. Therefore the difference of spin-dependent conductance probes the SSAR at the vortex core.
In the experiment, the direction of the ferromagnetic Fe-coated W tip's polarization and the direction of the magnetic field that excited the vortex and the MZM within are controlled to be both parallel and antiparallel to each other. The STS measurements are conducted in all four tip-polarization and magnetic field configurations. In the STS measurements on a 5 QL Bi 2 Te 3 /NbSe 2 , the ZBP at vortex core clearly shows a conductance difference between the parallel and antiparallel situations, as shown in Fig. 8b. This result demonstrates the SSAR process caused by a spin polarized tunneling current reflected by the MZM. It is worth noting that this SSAR process can only happen at the center of the vortex. When probed at a distance from the vortex core, no such effect can be seen, as shown in Fig. 8c. The experimental results fit quite well with the theoretical calculations (See Fig. 8d and e).
The SSAR effect was also carefully examined in some other TI/SC samples without MZM, i.e., 3 QL Bi 2 Te 3 /NbSe 2 sample and bare NbSe 2 substrate. In those samples, the conductance in the ZBP of two different tip polarization-magnetic field direction configurations (parallel and antiparallel) exhibit no obvious difference. As mentioned before, the absence of the SSAR rules out MZM's origin from the proximity effect or from the substrate. Similarly, in the 5 QL Bi 2 Te 3 /NbSe 2 sample, when the excitation field is increased to a point where two adjacent vortices is too close to each other, the MZMs are destroyed and the SSAR effect disappears accordingly. Therefore, all evidences consistently prove that the MZMs exist in the vortex.

SUMMARY AND OUTLOOK
Initially proposed by Fu and Kane in 2008, the artificial TSC system made of the 3D TI/SC heterostructure inspired vast fascinating theoretical works and significant experimental discovery in less than a decade. MZM has been proven to exist in TSC. Though a pwave TSC is yet out-of-reach in nature, topological superconductivity can be achieved by inducing superconductivity into the TSSs of TIs from normal s-wave SCs. The combination of state-of-the-art MBE and STM/STS technique plays a crucial role in the discovery of MZMs. Extreme high quality TI/SC heterostructures can be fabricated and transferred to an STM immediately for in-situ measurements. Low temperature condition (300 to 30 mK) and vector magnetic field help to observe the genuine SSAR properties of the mysterious MZMs. Although the detection of a single MZM lies in the minigap of other low energy quasiparticle states is still beyond current STS energy resolution, the nice agreements between theories and experimental observations and the selfconsistency in different independent experiments explicitly prove the existence of MZM in the artificial TSC.
The MZM in solid-state physics attracts great attention not only because it is a new quasiparticle in the diverse family formed by different topology phases of matter, it has serious applications in quantum computation as well. Its topological protected wavefunction can be used as qubits to proceed fault-tolerant quantum computation. 5,27 The non-Abelian braiding of MZMs that is needed in the fault-tolerant quantum computation is a challenging task. Recently, another approach called surface code provides an alternative way to apply MZMs in quantum computation. [49][50][51] This method does not require the non-Abelian statistics of MZMs, so it is not a fault-tolerant quantum computation. However, it still has remarkable ability for error detection and subsequent correction during qubit readout, as the nucleation of anyons through the action of a random operator can be reliably tracked through stabilizer measurements and it is much easier to achieve experimentally. 50 To realize the theoretical proposal that applies MZMs into quantum computation, much effort is needed based on the current 3D TI/SC platform. First, it requires novel fabrication methods to build a TI/SC heterostructure with fine structures to form a device that can be used to measure the interactions between multiple MZMs. Second, a robust way of manipulating MZMs is needed to perform the necessary gate operations for universal quantum computation. Third, it is important to explore new materials that host MZMs, i.e., a better SC substrate that has larger gap protecting the MZMs. Considering the fast development in the first decade in this field, these tasks will hopefully be accomplished in the near future.