Home | About | Journals | Submit | Contact Us | Français |

**|**Nat Commun**|**PMC5730620

Formats

Article sections

Authors

Related links

Nat Commun. 2017; 8: 976.

Published online 2017 October 17. doi: 10.1038/s41467-017-01209-9

PMCID: PMC5730620

K. Iwaya,^{}^{1} Y. Kohsaka,^{1} K. Okawa,^{2} T. Machida,^{1} M. S. Bahramy,^{1,}^{3} T. Hanaguri,^{}^{1} and T. Sasagawa^{2}

K. Iwaya, Email: pj.nekir@ayawi.

Received 2017 May 25; Accepted 2017 August 24.

Copyright © The Author(s) 2017

A bulk superconductor possessing a topological surface state at the Fermi level is a promising system to realise long-sought topological superconductivity. Although several candidate materials have been proposed, experimental demonstrations concurrently exploring spin textures and superconductivity at the surface have remained elusive. Here we perform spectroscopic-imaging scanning tunnelling microscopy on the centrosymmetric superconductor *β*-PdBi_{2} that hosts a topological surface state. By combining first-principles electronic-structure calculations and quasiparticle interference experiments, we determine the spin textures at the surface, and show not only the topological surface state but also all other surface bands exhibit spin polarisations parallel to the surface. We find that the superconducting gap fully opens in all the spin-polarised surface states. This behaviour is consistent with a possible spin-triplet order parameter expected for such in-plane spin textures, but the observed superconducting gap amplitude is comparable to that of the bulk, suggesting that the spin-singlet component is predominant in *β*-PdBi_{2}.

Superconductivity arising from non-degenerate spin-polarised Fermi surfaces is expected to consist of a mixing of spin-singlet and triplet order parameters as proposed in non-centrosymmetric superconductors^{1, 2}. If superconductivity is induced in the spin-polarised topological surface states (TSSs), such a situation may be realised and the induced spin-triplet order parameter could play a key role for the emergence of Majorana fermions at edge states or in vortex cores, offering potential applications such as topological quantum computing^{3–6}. To accomplish topological superconductivity, various candidate systems including carrier doped topological insulators^{7–9} and superconductor/topological insulator heterostructures^{10–14} have been investigated, in addition to other systems such as magnetic^{15} and semiconducting^{16} nanowires fabricated on superconductors.

Among these systems, stoichiometric superconductors possessing the TSSs at the Fermi level *E*
_{F} are promising candidates. The non-centrosymmetric superconductor PbTaSe_{2}
^{17} is one of such materials for which a recent spectroscopic-imaging scanning tunnelling microscopy (SI-STM) study has revealed the existence of both the TSS and a fully opened superconducting (SC) gap^{18}. At the cleaved surface of PbTaSe_{2}, the inversion symmetry is broken not only along the surface normal direction but also along the in-plane direction because of the non-centrosymmetric crystal structure. This causes intricate spin-orbit coupling that gives rise to out-of-plane spin components. The SC order parameter favoured in such a situation is interesting but complicated.

Here we study another stoichiometric superconductor *β*-PdBi_{2} with a SC transition temperature *T*_{c}~5.4K^{19}. *β*-PdBi_{2} crystallises into centrosymmetric tetragonal lattice with space group I4/mmm. Since in-plane inversion symmetry is always maintained, the spin texture induced by spin-orbit coupling is expected to be aligned in the in-plane direction at the surface. Angle-resolved photoemission spectroscopy (ARPES) has revealed the presence of the TSS as well as the topologically trivial surface state at *E*
_{F}
^{20}. The SC gap has been studied by scanning tunnelling spectroscopy experiments^{21, 22}. However, reported tunnelling spectra do not agree with each other; there is only one BCS-like gap at the surface of the bulk single crystal cleaved in the air^{21}, whereas two gaps are identified in the case of the in situ molecular-beam-epitaxy-grown thin film^{22}. Moreover, the simultaneous investigation of the TSS and the SC gap, which is critically important to discuss topological superconductivity, has not yet been done. We have performed SI-STM on clean surfaces of *β*-PdBi_{2} prepared by cleaving bulk single crystals in ultra high vacuum to study the SC gap opening in the TSS. Quasiparticle interference (QPI) imaging combined with numerical simulations has been utilised to explore the spin textures at the surface. We find that the observed QPI patterns are qualitatively consistent with the spin textures predicted by the first-principles calculations, which reveal that not only the TSS but also all other states at the surface are spin polarised. This suggests that the surface of *β*-PdBi_{2} can potentially harbour the mixed spin-singlet and triplet superconductivity. A fully opened single SC gap Δ~0.8meV has been identified in the tunnelling spectrum without any detectable residual spectral weight near *E*
_{F}, being consistent with the result of air-cleaved sample^{21}. These results indicate that the electronic states detected by SI-STM including the TSS are all spin-polarised and exhibit fully gapped superconductivity. We argue that the spin-triplet order parameter expected from the spin textures is consistent with the observed isotropic gap. The observed gap amplitude, however, is comparable to that of the bulk^{23}, suggesting that the spin-singlet component is dominant. Therefore, the topological nature of the observed superconductivity, if any, may manifest itself at lower temperatures.

We first examine the spin textures at the surface of *β*-PdBi_{2} by first-principles calculations within the framework of the density functional theory. The surface exposed after cleaving a single crystal along (001) plane is expected to be composed of Bi atoms since the weakest bond in the crystal structure is the van der Waals bond between adjacent Bi layers (Fig. 1a inset). Therefore, we mainly focus on low-energy electronic states composed of Bi 6p orbitals at the top-most layer predominating the tunnelling current. (Details of the calculation is described in Methods section.)

Calculated surface electronic states of the top-most Bi layer of *β*-PdBi_{2}. **a** Energy dispersions of the surface (topologically trivial S1 and nontrivial S2) and the bulk states (*α*, *β* and *γ*) along the $\overline{\mathrm{\Gamma}}$− **...**

The calculated band dispersions and the constant energy contours (CECs) at energy *E*=+30meV are shown in Fig. 1a, b, respectively. Each of the surface and bulk states is labelled in the same manner as the previous report^{20}. The surface states are found in the close vicinity of the bulk bands. The trivial surface state S1 is located near the bulk *β* band, whereas the TSS S2 is located at the bottom of the bulk *γ* and *δ* bands. For comparison with the experimental results, we performed a Wannier transformation as written in Methods section. Figure 1c shows the spectral function at +30meV with the Wannier transformation. The spectral weights of bulk *γ* and *δ* bands are very small. For the TSS S2 surrounding $\overline{\mathrm{M}}$ points, the weight is negligibly small except for the section parallel to the $\overline{\mathrm{M}}$−$\overline{\mathrm{X}}$ direction.

Figure 1d–f show *x*, *y* and *z* components of the spin polarisations weighted by the spectral function. It is clear that the spin polarisations are predominantly aligned in-plane, as expected from the symmetry of the crystal structure. It should be noted that not only the TSS S2 but also S1, and even the bulk bands are all spin polarised. The spin orientations of S1 and S2 are consistent with the results of spin-resolved ARPES^{20}. The spin polarisations of the bulk bands are likely to be induced by local broken inversion symmetry at the Bi site in the bulk as well as at the surface^{24} (Fig. 1a inset). The spin textures of these bulk bands have never been studied before and play an important role in QPI as we shall show below.

In order to confirm the calculated spin textures experimentally, we utilise QPI imaging. The QPI effect is nothing but the formation of electronic standing waves. In the presence of elastic scatterers, two quantum states on the same CEC with a scattering wave vector **q** may interfere and generate a standing wave with a wavelength of 2*π*/|**q**|. This yields *E*-dependent spatial oscillations in the local density-of-states distribution that manifest themselves in the tunnelling conductance d*I*/d*V*|_{V=E/e} images, where *I* is the tunnelling current, *V* is the bias voltage and *e* is the elementary charge. Fourier analysis allows us to identify various QPI signals appearing at different **q** vectors. The pre-requisite to have QPI signals at **q** is that **q** must connect two states with high enough spectral weights. In addition, the spin polarisations of these two states should not be antiparallel because spin-antiparallel states are orthogonal and thereby cannot interfere with each other. In this way, the QPI imaging can provide information on the spin textures in the **k** space^{25, 26}.

Here we outline expected QPI patterns from the calculated spin textures at the surface (Fig. 1d–f). In the case of non-magnetic scatterers, any intra-band back-scattering of spin-polarised bands is forbidden because the spin polarisations of the two states relevant for scattering are always antiparallel due to time-reversal symmetry. The inter-band back-scatterings from S1 to $\overline{\mathrm{\Gamma}}$-centred S2 and from *α* to $\overline{\mathrm{\Gamma}}$-centred S2 are also forbidden, because the spin orientations of S1, S2 and *α* are the same in the $\overline{\mathrm{\Gamma}}$−$\overline{\mathrm{M}}$ direction. The inter-band forward-scatterings from S1 to $\overline{\mathrm{\Gamma}}$-centred S2 and from S1 to *α* are allowed but the expected **q** vectors are small; it may be difficult to distinguish these QPI signals from extrinsic small **q** modulations associated with inhomogeneously distributed unavoidable defects in a real material. As a result, three scattering channels may govern the QPI patterns; inter-band forward-scattering from *α* to $\overline{\mathrm{\Gamma}}$-centred S2 (**q**
_{1}), inter-band back-scatterings from S1 to *β* band (**q**
_{2}) and from $\overline{\mathrm{\Gamma}}$-centred S2 to *β* band (**q**
_{3}) (Fig. 1c).

We start our experiments by checking the quality of the cleaved surface using constant-current STM imaging. An atomically flat area as large as 100×100nm^{2} is observed as shown in Fig. 2a. In a magnified image (Fig. 2a inset), at least two kinds of defects are observed as a local suppression and a subtle protrusion surrounded by a suppressed area. Although the exact nature of these defects is unknown, they may work as quasiparticle scatterers that cause QPI. Fourier-transformed STM image (Supplementary Fig. ^{1}) exhibits Bragg peaks that are consistent with the in-plane lattice constant *a*
_{0}=0.337nm^{19}. This allows us to identify the crystallographic directions in real and reciprocal spaces.

Spectroscopic imaging on the *β*-PdBi_{2} surface. **a** A constant-current STM image of the *β*-PdBi_{2} surface prepared by cleaving in ultra high vacuum. Set-up conditions for imaging were *V*=+200mV, *I*=200pA. **...**

QPI imaging has been performed in the same field of view of Fig. 2a. Figure 2b, c shows a d*I*/d*V* map at *E*=+30meV and its Fourier-transformed image, respectively. We identify three distinct line-like QPI signals in Fig. 2c, which are parallel with each other and normal to the $\overline{\mathrm{\Gamma}}$−$\overline{\mathrm{M}}$ direction. These three QPI signals show hole-like dispersions (d|**q**|/d*V*
_{s}<0) as seen in Fig. 2d. The line-like shapes and hole-like dispersions are consistent with the square-shaped CECs of S1, S2, *α* and *β* bands shown in Fig. 1b and the band dispersions shown in Fig. 1a, respectively. Namely, the observed QPI signals may correspond to **q**
_{1}, **q**
_{2} and **q**
_{3}.

We confirm this conjecture by comparing experimental data with the numerically simulated QPI patterns calculated by using the results of first-principles calculations with the standard *T*-matrix formalism. In addition to a full (standard) simulation, we have also performed a spinless simulation. The spin degree of freedom is included in the former and hypothetically suppressed in the latter. The contrast between the two simulations highlights the role of the spin textures. (See Methods section for details.)

The result of the full simulation at *E*=+30meV is shown in Fig. 3a. Three parallel line-like QPI signals normal to the $\overline{\mathrm{\Gamma}}$−$\overline{\mathrm{M}}$ direction are identified. These features correspond to **q**
_{1}, **q**
_{2} and **q**
_{3} indicated in Fig. 1c and their energy evolutions reasonably agree with the experimental observation (Fig. 3b), indicating that **q**
_{1}, **q**
_{2} and **q**
_{3} scatterings indeed dominate the QPI patterns.

Simulated QPI patterns of *β*-PdBi_{2} surface. **a** Simulated QPI pattern at *E*=+30meV with taking spin degrees of freedom into account (full simulation). The locations of **q**
_{1}, **q**
_{2} and **q**
_{3} are shown by arrows. A white square **...**

The key role of the spin textures becomes evident if we compare the result of full simulation with that of spinless simulation (Fig. 3c), and the experimental data. Figure 3d compares line profiles along the $\overline{\mathrm{\Gamma}}$−$\overline{\mathrm{M}}$ direction from experimental and simulated QPI patterns. In the spinless simulation, there are three sharp peaks but their |**q**| values do not agree with the experimental observations. These sharp peaks correspond to the intra-surface-state back-scatterings in S1 (**q**
_{S1−S1}) and in $\overline{\mathrm{\Gamma}}$-centred S2 (**q**
_{S2−S2}), and inter-band back-scattering from S1 to $\overline{\mathrm{\Gamma}}$-centred S2 (**q**
_{S1−S2}). Because the spin orientations are antiparallel between initial and final states of these scatterings, they are suppressed in the full simulation that well captures the experimental observations.

The above comprehensive approach combining experiments and calculations provide the following important implications. First, our SI-STM data primarily reflect the electronic state of the top-most Bi layer as we assumed. Second, the spin degrees of freedom is crucial for QPI. Finally, the calculated spin textures shown in Fig. 1d–f adequately capture the real spin textures at the surface. Notably, the fourfold QPI pattern is similarly identified at *E*
_{F} in the normal state (*T*=1.5K, *B*=12T) as shown in Supplementary Fig. ^{2}. This indicates that the spin textures discussed here indeed exist at *E*
_{F} in the normal state.

Given the surface sensitivity of our measurement evidenced in the QPI patterns, we are able to argue the nature of superconductivity at the surface. Figure 4a shows the temperature *T* evolution of the tunnelling spectrum in the SC state. At *T*=0.4K, the SC gap fully opens and there is no residual spectral weight inside the gap, indicating that all of the states at the surface are gapped. Each spectrum can be fitted well with the Dynes function for the single isotropic gap, meaning that the SC gap is **k** independent. We also investigate QPI patterns near the SC gap energy and find no clear superconductivity-induced QPI signals except at **q**=0 (Supplementary Fig. ^{2}), which means that the SC gap is spatially uniform. The SC gap amplitude Δ(*T*) is estimated to be 0.8meV at *T*=0.4K, being reasonably consistent with the result of air-cleaved sample^{21}. We note that any noticeable spectroscopic features are not observed near *E*
_{F} at step edges (Supplementary Fig. ^{3}). The temperature dependence of the SC gap well follows the BCS behaviour (Fig. 4b).

Superconducting states at the surface. **a** d*I*/d*V* spectrum as a function of temperature. Set-up conditions were *V*=+10mV, *I*=200pA. Bias modulation amplitude *V*
_{r.m.s.} was 18μV. Solid lines **...**

We next examine the electronic state of vortices that may host Majorana fermions^{6, 27, 28}. Isotropic vortices are clearly imaged in the d*I*/d*V* map at *E*=0meV (Fig. 4c inset), as being similar to the previous result obtained at the surface of the air-cleaved bulk crystal^{21}. The line profile of the vortex core is also reasonably consistent with the previous result^{29} (Supplementary Fig. ^{4}). To estimate the in-plane coherence length and compare it with that obtained from *H*
_{c2} measurements, detailed magnetic field dependence of the vortex core size is required as proposed recently^{29}. Considering the good agreement with our result and those in ref. ^{29}, we believe that a similar coherence length as discussed in refs. ^{21, 29} would be expected in our samples. The spectrum taken in the vortex core is almost flat (Fig. 4c). This result indicates that the core is in the dirty limit (mean free path *l*<*ξ*) where vortex bound states as well as the Majorana zero mode are not well defined. Although the zero-energy local density-of-states peak in the vortex core cannot be an evidence of the Majorana zero mode by itself^{18}, it would be possible to investigate the spin structures that are unique to the Majorana state^{27, 28}, for example. We anticipate that *β*-PdBi_{2} will be a good touchstone for these theories in future because it is a stoichiometric material that can be made cleaner in principle.

Our concurrent investigations of the surface electronic states and superconductivity reveal that the SC gap fully opens in all of the spin-polarised surface states. Since the mixing of spin-singlet and triplet order parameters is generally expected in the presence of non-degenerate spin-polarised surface states^{1, 2} and the triplet component can possess nodes, it is intriguing to argue the possible nodal structure of the SC states of *β*-PdBi_{2}. It has been shown that the spin-triplet order parameter **d**(**k**) favours to be aligned along the spin direction when space inversion symmetry is broken^{30}. If this is the case at the surface of *β*-PdBi_{2} where the point group symmetry is lowered from *D*
_{4h} to *C*
_{4v}, the triplet pairing is of *p*-wave type with nodes along the out-of-plane direction^{30}. Since the surface states are two dimensional in nature, the out-of-plane nodes would not be active and the SC gap should look like isotropic as observed.

Next we discuss the amplitude of the SC gap. The SC gap amplitude in the presence of the parity mixing is given by |Ψ_{s}(**k**)|±|**d**(**k**)| where Ψ_{s} denotes *s*-wave order parameter that represents the bulk gap^{2, 30–32}. Thus the difference between the SC gap amplitudes of bulk and surface gives us an estimate of the amplitude of the *p*-wave component. The bulk SC gap has been estimated to be about 0.9meV by the specific heat measurement^{23}, which is close to the value of 0.8meV that we observed at the surface. According to the theory^{32}, the difference between the surface SC gap and the bulk counterpart depends on how far the surface state is separated from the bulk band at *E*
_{F} in **k** space. We estimate the quantity $\delta \equiv \left({k}_{\mathrm{F}}^{\mathrm{D}}-{k}_{\mathrm{F}}\right)\mathrm{\u2215}{k}_{\mathrm{F}}$, which is introduced in ref. ^{32}, to be about 0.05 (${k}_{\mathrm{F}}^{\mathrm{D}}$ and *k*
_{F} denote the Fermi momenta of the TSS and the bulk band, respectively). This is reasonably small to explain the similar SC gap amplitudes between the surface and the bulk. Namely, even though it is allowed, the *p*-wave component at the surface of *β*-PdBi_{2} is small. The absence of edge states shown in Supplementary Fig. ^{3} might be due to the small *p*-wave component. In order to detect this small *p*-wave component, if any, we need to avoid comparing the results of two different experiments, STM and specific heat, as well as to improve the energy resolution. As the bulk bands appear at the surface, too, it would be possible to detect the surface and bulk SC gaps simultaneously by STM alone, provided higher energy resolution could be achieved. To this end, the use of a SC tip at lower temperatures is important.

The *β*-PdBi_{2} single crystals were grown by a melt growth method^{20}, and characterised by x-ray diffraction (XRD) and transport measurements. Pd (3N5) and Bi (5N) at a molar ratio of 1:2 were encapsulated in an evacuated quartz tube, pre-reacted at temperature above 1000°C until it completely melted and mixed. Then, it was kept at 900°C for 20h, cooled down at a rate of 3°C/h down to 500°C, and finally rapidly quenched into cold water. PdBi_{2} has two different crystallographic and superconducting phases: *α*-phase with the space group C2/m (*T*_{c}~1.7K^{33}) and *β*-phase with I4/mmm (*T*_{c}~5.4K). The last quenching procedure is important to selectively obtain the *β*-phase. Any trace amount of *α*-phase was not detected by XRD. For STM measurements, we used single crystals with residual resistivity ratio of 15, larger than those of the previous reports^{19, 21, 23}. The crystals were cleaved along the (001) plane at room temperature in ultra high vacuum conditions to obtain clean and flat surfaces needed for STM. A commercial ^{3}He-based STM system (UNISOKU USM-1300) modified by ourselves^{34} was employed in this study. We used electrochemically etched tungsten tips, which were cleaned and sharpened by field ion microscopy. The tips were subsequently treated and calibrated on clean Au(100) surfaces before used for *β*-PdBi_{2}. We applied bias voltages to the sample whereas the tip was virtually grounded at the current-voltage converter (Femto LCA-1K-5G). Tunnelling spectra were measured by the software-based lock-in detector included in the commercial STM control system (Nanonis).

To calculate the surface electronic structure and its corresponding QPI pattern, we first performed a DFT calculation using the Perdew−Burke−Ernzerhof exchange-correlation functional^{35} as implemented in the WIEN2K program^{36}. Relativistic effects including spin-orbit coupling were fully taken into account. For all atoms, the muffin-tin radius *R*
_{MT} was chosen such that its product with the maximum modulus of reciprocal vectors *K*
_{max} become *R*
_{MT}
*K*
_{max}=7.0. Considering the tetragonal phase of *β*-PdBi_{2}, the corresponding Brillouin zone was sampled using a 20×20×5 *k*-mesh. For the surface calculations, a 100 unit tight-binding supercell was constructed using maximally localised Wannier functions^{37–39}. The 6p orbitals of Bi atoms were chosen as the projection centres.

We simulated QPI patterns by incorporating the eigenvalues and eigenvectors obtained from our tight-binding supercell calculations into the standard *T*-matrix formalism^{40}. The top-most Bi 6p orbitals were considered in the calculations. For all the simulations, the broadening factor was chosen to be 5meV, and a localised, spin-preserving and orbital-preserving scatterer of strength 0.1eV were employed. For comparison with the experimental results, we performed a basis transformation from the lattice model to the continuum model with the Wannier function^{41} constructed by projecting the Bi 6p orbitals. We calculated the local density of state at 0.5nm above the top-most Bi atoms. The basis transformation is also applied to the spectral function shown in Fig. 1c to present the electronic structure used for the QPI calculation. For the spinless simulation, the sum of the Green’s functions constructed from the spin-up and spin-down subspaces were used for the *T*-matrix calculations.

All relevant data are available on request, which should be addressed to K.I.

The authors thank M. Sakano, K. Ishizaka, T. Mizushima, Y. Tanaka and Y. Yanase for valuable discussions and Y. Okada for experimental advices. This work was supported by a CREST project (JPMJCR16F2) from Japan Science and Technology Agency (JST) and JSPS KAKENHI Grant Nos.16K05465, 16H06012, and 16H03847. M.S.B. was supported by CREST, JST (Grant No. JPMJCR16F1).

Author contributions

K.I. carried out the experiments and the data analyses with assistance from Y.K., T.M. and T.H. Y.K., T.S., and M.S.B. carried out the first-principles calculations and numerical simulations. K.O. and T.S. synthesised and characterised the single crystals. T.H. supervised the project. K.I., Y.K. and T.H. wrote the manuscript.

The authors declare no competing financial interest.

A correction to this article is available online at https://doi.org/10.1038/s41467-017-02137-4.

**Electronic supplementary material**

**Supplementary Information** accompanies this paper at 10.1038/s41467-017-01209-9.

**Change history:** A correction to this article has been published and is linked from the HTML version of this paper.

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

1. Gor’kov LP, Rashba EI. Superconducting 2D system with lifted spin degeneracy: mixed singlet-triplet state. Phys. Rev. Lett. 2001;87:037004. doi: 10.1103/PhysRevLett.87.037004. [PubMed] [Cross Ref]

2. Sigrist M, et al. Superconductivity in non-centrosymmetric materials. J. Magn. Magn. Mater. 2007;310:536–540. doi: 10.1016/j.jmmm.2006.10.141. [Cross Ref]

3. Kitaev AY. Fault-tolerant quantum computation by anyons. Ann. Phys. 2003;303:2–30. doi: 10.1016/S0003-4916(02)00018-0. [Cross Ref]

4. Fu L, Kane CL. Superconducting proximity effect and Majorana Fermions at the surface of a topological insulator. Phys. Rev. Lett. 2008;100:096407. doi: 10.1103/PhysRevLett.100.096407. [PubMed] [Cross Ref]

5. Qi XL, Zhang SC. Topological insulators and superconductors. Rev. Mod. Phys. 2011;83:1057–1110. doi: 10.1103/RevModPhys.83.1057. [Cross Ref]

6. Hosur P, Ghaemi P, Mong RSK, Vishwanath A. Majorana modes at the ends of superconductor vortices in doped topological insulators. Phys. Rev. Lett. 2011;107:097001. doi: 10.1103/PhysRevLett.107.097001. [PubMed] [Cross Ref]

7. Wray LA, et al. Observation of topological order in a superconducting doped topological insulator. Nat. Phys. 2010;6:855–859. doi: 10.1038/nphys1762. [Cross Ref]

8. Levy N, et al. Experimental evidence for *s*-wave pairing symmetry in superconducting Cu_{x}Bi_{2}Se_{3} single crystals using a scanning tunneling microscope. Phys. Rev. Lett. 2013;110:117001. doi: 10.1103/PhysRevLett.110.117001. [PubMed] [Cross Ref]

9. Du G, et al. Drive the Dirac electrons into cooper pairs in Sr_{x}Bi_{2}Se_{3}. Nat. Commun. 2017;8:14466. doi: 10.1038/ncomms14466. [PMC free article] [PubMed] [Cross Ref]

10. Wang E, et al. Fully gapped topological surface states in Bi_{2}Se_{3} films induced by a d-wave high-temperature superconductor. Nat. Phys. 2013;9:621–625. doi: 10.1038/nphys2744. [Cross Ref]

11. Wang MX, et al. The coexistence of superconductivity and topological order in the Bi_{2}Se_{3} thin films. Science. 2012;336:52–55. doi: 10.1126/science.1216466. [PubMed] [Cross Ref]

12. Xu SY, et al. Momentum-space imaging of cooper pairing in a half-Dirac-gas topological superconductor. Nat. Phys. 2014;10:943–950. doi: 10.1038/nphys3139. [Cross Ref]

13. Xu JP, et al. Artificial topological superconductor by the proximity effect. Phys. Rev. Lett. 2014;112:217001. doi: 10.1103/PhysRevLett.112.217001. [Cross Ref]

14. Xu JP, et al. Experimental detection of a Majorana mode in the core of a magnetic vortex inside a topological insulator-superconductor Bi_{2}Te_{3}/NbSe_{2} Heterostructure. Phys. Rev. Lett. 2015;114:017001. doi: 10.1103/PhysRevLett.114.017001. [PubMed] [Cross Ref]

15. Nadj-Perge S, et al. Observation of Majorana fermions in ferromagnetic atomic chains on a superconductor. Science. 2014;346:602–607. doi: 10.1126/science.1259327. [PubMed] [Cross Ref]

16. Mourik S, et al. Signatures of Majorana fermions in hybrid superconductor-semiconductor nanowire devices. Science. 2012;336:1003–1007. doi: 10.1126/science.1222360. [PubMed] [Cross Ref]

17. Ali MN, Gibson QD, Klimczuk T, Cava RJ. Noncentrosymmetric superconductor with a bulk three-dimensional Dirac cone gapped by strong spin-orbit coupling. Phys. Rev. B. 2014;89:020505(R). doi: 10.1103/PhysRevB.89.020505. [Cross Ref]

18. Guan SY, et al. Superconducting topological surface states in the noncentrosymmetric bulk superconductor PbTaSe_{2}. Sci. Adv. 2016;2:e1600894. doi: 10.1126/sciadv.1600894. [PMC free article] [PubMed] [Cross Ref]

19. Imai Y, et al. Superconductivity at 5.4 K in *β*-Bi_{2}Pd. J. Phys. Soc. Jpn. 2012;81:113708. doi: 10.1143/JPSJ.81.113708. [Cross Ref]

20. Sakano M, et al. Topologically protected surface states in a centrosymmetric superconductor *β*-PdBi_{2}. Nat. Commun. 2015;6:8595. doi: 10.1038/ncomms9595. [PMC free article] [PubMed] [Cross Ref]

21. Herrera E, et al. Magnetic field dependence of the density of states in the multiband superconductor *β*-Bi_{2}Pd. Phys. Rev. B. 2015;92:054507. doi: 10.1103/PhysRevB.92.054507. [Cross Ref]

22. Lv YF, et al. Experimental signature of topological superconductivity and Majorana zero modes on *β*-Bi_{2}Pd thin films. Sci. Bull. 2017;62:852–856. doi: 10.1016/j.scib.2017.05.008. [Cross Ref]

23. Kačmarčik J, et al. Single-gap superconductivity in *β*-Bi_{2}Pd. Phys. Rev. B. 2016;93:144502. doi: 10.1103/PhysRevB.93.144502. [Cross Ref]

24. Zhang X, Liu Q, Luo JW, Freeman AJ, Zunger A. Hidden spin polarization in inversion-symmetric bulk crystals. Nat. Phys. 2014;10:387–393. doi: 10.1038/nphys2933. [Cross Ref]

25. Pascual JI, et al. Role of spin in quasiparticle interference. Phys. Rev. Lett. 2004;93:196802. doi: 10.1103/PhysRevLett.93.196802. [PubMed] [Cross Ref]

26. Roushan P, et al. Topological surface states protected from backscattering by chiral spin texture. Nature. 2009;460:1106–1109. doi: 10.1038/nature08308. [PubMed] [Cross Ref]

27. Nagai Y, Nakamura H, Machida M. Spin-polarized Majorana bound states inside a vortex core in topological superconductors. J. Phys. Soc. Jpn. 2014;83:064703. doi: 10.7566/JPSJ.83.064703. [Cross Ref]

28. Kawakami T, Hu X. Evolution of density of states and a spin-resolved checkerboard-type pattern associated with the Majorana bound state. Phys. Rev. Lett. 2015;115:177001. doi: 10.1103/PhysRevLett.115.177001. [PubMed] [Cross Ref]

29. Fente A, et al. Field dependence of the vortex core size probed by scanning tunneling microscopy. Phys. Rev. B. 2016;94:014517. doi: 10.1103/PhysRevB.94.014517. [Cross Ref]

30. Frigeri PA, Agterberg DF, Koga A, Sigrist M. Superconductivity without inversion symmetry: MnSi versus CePt_{3}Si. Phys. Rev. Lett. 2004;92:097001. doi: 10.1103/PhysRevLett.92.097001. [PubMed] [Cross Ref]

31. Hayashi N, Wakabayashi K, Frigeri PA, Sigrist M. Temperature dependence of the superfluid density in a noncentrosymmetric superconductor. Phys. Rev. B. 2006;73:024504. doi: 10.1103/PhysRevB.73.024504. [Cross Ref]

32. Mizushima T, Yamakage A, Sato M, Tanaka Y. Dirac-fermion-induced parity mixing in superconducting topological insulators. Phys. Rev. B. 2014;90:184516. doi: 10.1103/PhysRevB.90.184516. [Cross Ref]

33. Mitra S, et al. Probing the superconducting gap symmetry of *α*-PdBi_{2}: a penetration depth study. Phys. Rev. B. 2017;95:134519. doi: 10.1103/PhysRevB.95.134519. [Cross Ref]

34. Hanaguri T. Development of high-field STM and its application to the study on magnetically-tuned criticality in Sr_{3}Ru_{2}O_{7}. J. Phys.: Conf. Ser. 2006;51:514–521.

35. Perdew JP, Burke K, Ernzerhof M. Generalized gradient approximation made simple. Phys. Rev. Lett. 1996;77:3865–3868. doi: 10.1103/PhysRevLett.77.3865. [PubMed] [Cross Ref]

36. Balha, P., Schwarz, K., Madsen, G., Kvasnicka, D. & Luitz, J. WIEN2K package, Version 13.1 (2013).

37. Souza I, Marzari N, Vanderbilt D. Maximally localized Wannier functions for entangled energy bands. Phys. Rev. B. 2001;65:035109. doi: 10.1103/PhysRevB.65.035109. [Cross Ref]

38. Mostofi AA, et al. Wannier90: a tool for obtaining maximally localised Wannier functions. Comp. Phys. Commun. 2008;178:685–699. doi: 10.1016/j.cpc.2007.11.016. [Cross Ref]

39. Kunes J, et al. Wien2wannier: from linearized augmented plane waves to maximally localized Wannier functions. Comp. Phys. Commun. 2010;181:1888–1895. doi: 10.1016/j.cpc.2010.08.005. [Cross Ref]

40. Kohsaka Y, et al. Spin-orbit scattering visualized in quasiparticle interference. Phys. Rev. B. 2017;95:115307. doi: 10.1103/PhysRevB.95.115307. [Cross Ref]

Articles from Nature Communications are provided here courtesy of **Nature Publishing Group**

PubMed Central Canada is a service of the Canadian Institutes of Health Research (CIHR) working in partnership with the National Research Council's national science library in cooperation with the National Center for Biotechnology Information at the U.S. National Library of Medicine(NCBI/NLM). It includes content provided to the PubMed Central International archive by participating publishers. |