Observation of the spin-polarized surface state in a noncentrosymmetric superconductor BiPd

Abstract

Recently, noncentrosymmetric superconductor BiPd has attracted considerable research interest due to the possibility of hosting topological superconductivity. Here we report a systematic high-resolution angle-resolved photoemission spectroscopy (ARPES) and spin-resolved ARPES study of the normal state electronic and spin properties of BiPd. Our experimental results show the presence of a surface state at higher-binding energy with the location of Dirac point at around 700 meV below the Fermi level. The detailed photon energy, temperature-dependent and spin-resolved ARPES measurements complemented by our first-principles calculations demonstrate the existence of the spin-polarized surface states at high-binding energy. The absence of such spin-polarized surface states near the Fermi level negates the possibility of a topological superconducting behaviour on the surface. Our direct experimental observation of spin-polarized surface states in BiPd provides critical information that will guide the future search for topological superconductivity in noncentrosymmetric materials.

ARTICLE

Received 20 Apr 2015 | Accepted 22 Sep 2016 | Published 7 Nov 2016

Madhab Neupane1,*, Nasser Alidoust2,*, M. Mofazzel Hosen1,*, Jian-Xin Zhu3, Klauss Dimitri1, Su-Yang Xu2, Nagendra Dhakal1, Raman Sankar4, Ilya Belopolski2, Daniel S. Sanchez2, Tay-Rong Chang5, Horng-Tay Jeng5,6, Koji Miyamoto7, Taichi Okuda7, Hsin Lin8,9, Arun Bansil10, Dariusz Kaczorowski11, Fangcheng Chou4,M. Zahid Hasan2 & Tomasz Durakiewicz12,13

Recently, noncentrosymmetric superconductor BiPd has attracted considerable research interest due to the possibility of hosting topological superconductivity. Here we report a systematic high-resolution angle-resolved photoemission spectroscopy (ARPES) and spin-resolved ARPES study of the normal state electronic and spin properties of BiPd. Our experimental results show the presence of a surface state at higher-binding energy with the location of Dirac point at around 700 meV below the Fermi level. The detailed photon energy, temperature-dependent and spin-resolved ARPES measurements complemented by our rst-principles calculations demonstrate the existence of the spin-polarized surface states at high-binding energy. The absence of such spin-polarized surface states near the Fermi level negates the possibility of a topological superconducting behaviour on the surface. Our direct experimental observation of spin-polarized surface states in BiPd provides critical information that will guide the future search for topological superconductivity in noncentrosymmetric materials.

1 Department of Physics, University of Central Florida, Orlando, Florida 32816, USA. 2 Joseph Henry Laboratory and Department of Physics, Princeton University, Princeton, New Jersey 08544, USA. 3 Theoretical Division and Center for Integrated Nanotechnologies, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA. 4 Center for Condensed Matter Sciences, National Taiwan University, Taipei 10617, Taiwan. 5 Department of Physics, National Tsing Hua University, Hsinchu 30013, Taiwan. 6 Institute of Physics, Academia Sinica, Taipei 11529, Taiwan. 7 Hiroshima Synchrotron Radiation Center, Hiroshima University, 2-313 Kagamiyama, Higashi-Hiroshima 739-0046, Japan. 8 Centre for Advanced 2D Materials and Graphene Research Centre, National University of Singapore, Singapore 117546, Singapore. 9 Department of Physics, National University of Singapore, Singapore 117542, Singapore.

10 Department of Physics, Northeastern University, Boston, Massachusetts 02115, USA. 11 Institute of Low Temperature and Structure Research, Polish Academy of Sciences, 50-950 Wroclaw, Poland. 12 Condensed Matter and Magnet Science Group, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA. 13 Institute of Physics, Maria Curie - Sklodowska University, 20-031 Lublin, Poland. * These authors contributed equally to this work. Correspondence and requests for materials should be addressed to M.N. (email: mailto:[email protected]

Web End [email protected] ).

NATURE COMMUNICATIONS | 7:13315 | DOI: 10.1038/ncomms13315 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 1

DOI: 10.1038/ncomms13315 OPEN

Observation of the spin-polarized surface state in a noncentrosymmetric superconductor BiPd

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms13315

Recently, noncentrosymmetric (NCS) superconductors (SCs) have attracted considerable research interest due to the possibility that they host several exotic states13. The lack

of inversion symmetry creates an asymmetric potential gradient, which may split the electron bands by lifting the spin degeneracy, allowing hybrid-pairing of spin-singlet and spin-triplet states within the same orbital channel36. Furthermore, with the advent of topological insulators718, it was recently proposed that NCS SCs with a strong spinorbit coupling (SOC) are potential candidates for realizing topological superconductivity1922, holding promise for hosting protected Majorana surface states15,23,24. But so far, evidence of topological surface states in NCS materials is still lacking. Such convincing evidence may be obtained by angle- and spin-resolved photoemission, which provides an energy- and momentum-resolved probe of the electronic and spin structure.

A NCS system provides a good platform to look for exotic states such as the Weyl semimetal in TaAs25 and the topological nodal-line phase in PbTaSe2(ref. 26). Another NCS SC BiPd can provide a platform to study the interplay of SOC effects with superconductivity. This compound undergoes a structural transition from b-BiPd (orthorhombic) to a-BiPd (monoclinic) at 210 C and then becomes superconducting at the transition temperature (Tc)B3.7 K (refs 4,5). In comparison to many other

NCS SCs, BiPd is a weakly correlated compound possessing a heavy atom, Bi. Measurements by point-contact spectroscopy27 and nuclear quadrupole resonance28 indicate a complex gap structure in BiPd, which might be caused by the lack of inversion symmetry4,5. Recent scanning tunnelling microscopy measurements reveal that the superconducting state of BiPd appears to be topologically trivial, consistent with BardeenCooperSchrieffer theory with an s-wave order parameter6. Further momentum and spin-resolved experimental evidence is highly desirable to establish the presence of topological Dirac surface states, and to study the underlying microscopic mechanism of superconductivity. A detailed high-resolution angle-resolved photoemission spectroscopy (ARPES) and spin-resolved ARPES study is necessary to prove or disprove the topological nature of the surface states in the normal state of BiPd. Such a characterization of the normal state is currently missing, but necessary and needed as a rst step towards opening a discussion on the relationship between superconductivity and the topological properties of NCS systems.

In this paper, we report the experimental observation of spin-polarized surface states in the NCS material BiPd using ARPES and spin-resolved ARPES. Our experimental results show the presence of a surface state at higher-binding energy, with the Dirac node located at B700 meV below the

Fermi level. Furthermore, the Dirac node located at B650 meV below the Fermi level is also observed for different surface cleaving, providing direct evidence for the noncentrosymmetic crystal structure presented in BiPd. Our measurements suggest that the observed states may be topological in nature. Our results are further supported by our rst-principles calculations. The absence of spin-polarized surface states near the Fermi level of BiPd negates the possibility that this system might host topological superconductivity on the surface.

ResultsCrystal structure and sample characterization. The crystal structure of BiPd at low temperatures (o210 C) has a monoclinic unit cell with a 5.63 , b 10.66 , c 5.68 ,

a g 90 and b 101 with the b axis being as its unique axis

(Fig. 1a; ref. 4). Detailed characterization of the single crystals

used in our study indicated their high quality (Supplementary Figs 17; Supplementary Notes 14; Supplementary References). They exhibit a simple metallic behaviour in the normal state and a sharp superconducting transition at TcB3.7 K. The low-temperature magnetic susceptibility data is shown in Fig. 1b (see also Supplementary Fig. 1 and Supplementary Note 1 for additional data and discussion). A schematic bulk Brillouin zone is shown in Fig. 1c, where the projected surface along (010) is also illustrated. To experimentally identify its electronic structure, we study the electronic structure of BiPd on the cleaved (010) surface. Figure 1d (a picture of the single crystal measured is displayed in the inset of Fig. 1d) shows momentum-integrated ARPES spectral intensity over a wide-energy window. Sharp ARPES intensity peaks at binding energies EBB23 and

26 eV corresponding to the bismuth 5d3/2 and 5d5/2 energy levels

are observed.

Electronic structure of BiPd. We study the overall electronic structure of BiPd using ARPES. Figure 2a shows an ARPES dispersion map in a 1.3 eV binding energy window, where several dispersive bands within the valence band are identied. Moreover, several crossing or metallic bands in the vicinity of the Fermi level are observed. Remarkably, a nearly linearly dispersive Dirac cone-like state is observed at the Brillouin zone centre, showing a Dirac node located at a binding energy of EBB700

meV. The Dirac-like state can be observed in the region of the blue rectangle in the Fig. 2a, at the binding energy region of 500 900 meV. At the Fermi level only the metallic bands,

but no other Dirac-like linearly dispersive bands are observed. On the other hand, the linearly dispersive Dirac-like bands are found to be in the region of higher-binding energies.

Figure 2be shows the Fermi surface and constant energy contour plots, respectively. In the vicinity of the Fermi level, many metallic bands are observed (Fig. 2b). We also study the ARPES measured constant energy contour maps (Fig. 2ce). At the Fermi level, the constant energy contour consists of many metallic pockets. With increasing binding energy at B600 meV, the circular pocket formed by the upper Dirac cone is observed (Fig. 2c). On increasing the binding energy, the size of the pocket decreases and eventually shrinks to a point (the Dirac point) near EBB700 meV (Fig. 2d). At even higher-binding energies, the nearly circular pocket formed by the lower Dirac cone is observed (for EBB750 meV in Fig. 2e).

To reveal the nature of the states observed in BiPd, we performed photon energy and temperature-dependent ARPES, as well as spin-resolved ARPES measurements that are further complemented by calculations. Figure 3 shows the energy-momentum cuts measured with varying photon energies from 30 to 58 eV with a 4 eV energy step (see Supplementary Figs 24 and Supplementary Note 2 for additional ARPES data and discussion). Clear Ek dispersion of the bulk bands is observed. Remarkably, the dispersion of the linearly dispersive states at high-binding energy (500900 meV) is found to be unchanged with respect to the varying photon energy, supporting the two-dimensional nature of this state. It is important to note that the Dirac-like two-dimensional states are not found to be perfectly linear in energy-momentum axis. Furthermore, it is important to recall that in real materials such as pure Bi, graphene or topological insulators, the Dirac cones are never perfectly linear over a large energy window yet they can be approximated as linear within a narrow energy window around the Dirac point. This linear part represents the massless dispersion, in contrast to the large effective mass of conventional band electrons in other materials. Moreover, we note that the weak asymmetry of the surface bands with respect to k and k is

2 NATURE COMMUNICATIONS | 7:13315 | DOI: 10.1038/ncomms13315 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms13315 ARTICLE

a b

0.00

FC

[afii9846] (e.m.u. g1 )

0.05

0.10

0.15

ZFC

0.20

30103

2.5 3.0 Temperature (K)

3.5 4.0

Bi

1.5 2.0

b

c

Pd

a

c

d

(010)

Bulk

Bi 5d3/2

Bi 5d5/2

Z

Z X

25

Intensity (a.u.)

20

15

10

(010)

5

0

30 25 20 15 Energy (eV)

10 5 0

Figure 1 | Crystal structure and sample characterization of BiPd. (a) Crystal structure of BiPd with the b axis shown as its unique axis. It crystallizes in a monoclinic structure at low temperatures. (b) The magnetic susceptibility as a function of temperature showing a sharp superconducting transition temperature at B3.7 K in the eld of 20 Gs. (c) Schematic drawing of surface and bulk Brillouin zones, corresponding to a primitive cell, are shown.

We note that all electronic structure calculations were performed on a conventional unit cell for the bulk and a supercell for the slab. High-symmetric points are also marked. (d) Core level spectroscopic measurement of BiPd showing sharp peaks of Bi 5d. The inset shows a photograph of theBiPd sample.

probably coming from the matrix element effects, and from the fact that these measured spectra were slightly off from the

X

X direction.

To test the robustness of the surface state observed in BiPd, we have performed a temperature-dependent measurement as shown in Fig. 4a. On raising the temperature, the Dirac-like surface states survive even at room temperature, which establishes that the Dirac-like surface states are robust to the rise of temperature and potential surface contamination, resulting from increased partial pressure of residual gases released upon heating, see Fig. 4a. Furthermore, the observed Dirac-like states are found to be robust against thermal cycling (2030020 K), since lowering the temperature back down to20 K results in the similar spectra with the strong presence of Dirac-like state features (see the panel of Fig. 4a with the note of Re_20K).

NCS electronic signature. To better understand the electronic structure observed with ARPES, we perform rst-principles calculations of the bulk band structure (Supplementary Fig. 5; Supplementary Note 3) and slab calculations of BiPd using the generalized gradient approximation plus SOC method

(see Fig. 4b,c for the calculated electronic structure for the top and bottom surfaces, respectively). We note that the band structure calculations were performed on a conventional unit cell for the bulk system and on a supercell for the slab structure, which leads to folding of the states at the S point to the G point.

Our slab calculations show that the surface states are predicted to be at the G point with the location of the Dirac nodes at B0.5 and0.6 eV below the Fermi level, one of which comes from the top surface and the other from the bottom surface. The Dirac node of the surface state coming from the top and bottom surfaces are located at different binding energies (not degenerate), which comes from the fact that BiPd lacks inversion symmetry.

Since the calculations were made for slab geometry, the obtained band structure always involves bands from two surfaces, one on the top and the other on the bottom of the slab. Experimentally, one can only measure one surface (corresponding to either top or bottom surface of the slab geometry) at a time. In principle, it should be possible to measure either the top or bottom surface on the two sample pieces based on the cleaving. Therefore, depending on whether the cleaved surface is a top or bottom one, the location of the Dirac point energy is predicted to be different. Experimentally, the location of the Dirac point energy is found to be different based on the

NATURE COMMUNICATIONS | 7:13315 | DOI: 10.1038/ncomms13315 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 3

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms13315

a b c

0.4

600 meV

50 eV

0.5

0.0

0 meV

1 )

0.0

0.4

0.8

y(

k 0.5

1.0

0.5

0.0

0.5

Energy (eV)

d

e

700 meV 750 meV

1 )

k y(

1.2

1.0

0.8 0.4 0.0

Momentum (1)

0.4

1.0 0.5 0.0 0.5

kx(1)

1.0 0.5 0.0 0.5

kx(1)

Figure 2 | Electronic structure of BiPd. (a) Dispersion map of BiPd along the zone centre obtained by using incident photon energy of 50 eV at a temperature of 10 K. The blue rectangle in the binding energy range of B500900 meV shows the linearly dispersive states. (b) Fermi surface map.

(c) Constant energy contour at binding energy of 600 meV shows the intensity map in the region above the Dirac point. (d,e) The constant energy contours at 700 and 750 meV show the intensity map at around and below the Dirac point, respectively. The values of binding energies are noted on the plots. Data were collected at ALS BL 10.0.1.

0.4 30 eV 34 eV 38 eV 42 eV

Energy (eV)

0.0

0.4

0.8

0.4

0.0

0.4

0.8

1.2

1.0 0.0

0.0

0.4

0.8

1.2

1.0 0.0 1.0 0.0

0.4 0.4

0.0

0.4

0.8

1.2

1.2

0.4

0.0

0.4

0.8

0.8 0.0

0.4

0.0

0.4

0.8

1.2

0.4

0.0

0.4

0.8

1.2

Energy (eV)

0.4

0.0

0.4

0.8

1.2

46 eV 50 eV 54 eV 58 eV

1.2

1.0 0.0 1.0 0.0 1.0 0.0

0.8 0.0

Momentum (1)

Figure 3 | Photon energy-dependent ARPES dispersion maps. The measured photon energies are noted on the plots. The linearly dispersive states at a binding energy of B700 meV do not show any dispersion with photon energy, which indicate its two-dimensional nature. These data were collected along

G

X high-symmetry direction at ALS BL 10.0.1 at a temperature of 10 K. Additional photon energy data are presented in Supplementary Information.

cleaving either the top or the bottom surface, in agreement with our calculations. Actually, we experimentally observed the different Dirac point energy for top and bottom surface cleaving. For the top surface the Dirac point is B700 meV, while for the

bottom surface cleaving Dirac point is B650 meV from the Fermi level (Supplementary Fig. 4; Supplementary Note 3). We note that such two terminations of the crystal surface in BiPd are directly observed by STM topographic image6.

4 NATURE COMMUNICATIONS | 7:13315 | DOI: 10.1038/ncomms13315 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms13315 ARTICLE

a

0.0

Energy (eV)

Energy (eV)

1.2

0.4

0.8

0.4

20 K 200 K 300 K Re_20 K

0.0 0.4 0.4 0.0 0.4

Momentum (1)

0.4 0.0 0.4 0.4 0.0 0.4

b c

0.3

0.2

0.1

0

0.1

0.2

0.3

0.4

0.5

0.6

0.7

0.8

0.3

0.2

0.1

0

0.1

0.2

0.3

0.4

0.5

0.6

0.7

0.8X Z X Z

Figure 4 | Temperature-dependent spectra and rst-principles calculations. (a) ARPES energy-momentum dispersion maps measured using a photon energy of 30 eV along the

X

G

X momentum space cut-direction with varying temperatures. The measured values of temperature are noted on the plots. The panel with the note of Re_20K is the spectrum measured after thermal cycling (20-300-20 K). These data were collected at SSRL BL5-4 with a photon energy of 30 eV. (b) Slab calculations of BiPd for the top surface and (c) the bottom surface along the high-symmetry lines. The blue (red) dots represent surface states for the top (bottom) surface. Details of the calculation method is given in the Method section and additional calculation plots are shown in Supplementary Information.

20.8

a b c

0.6

0.6

1 2

20.0 20.2 20.4

0.4

0.4

0.2

0.2

Kinetic energy (eV)

20.6 1 2

20.4

20.2

20.0

Intensity (a.u.)

0.0

Polarization

Intensity (a.u.)

0.0

0.2

0.2

Polarization

19.8

19.6

0.4

0.4

0.6

0.6

0.2 0.4 19.8

0.4 0.2 0.0

Momentum (1)

19.8 20.0 20.2 20.4

Kinetic energy (eV)

Kinetic energy (eV)

Figure 5 | Spin-resolved ARPES results. (a) Spin-integrated ARPES spectrum of surface states. The blue dashed lines and number indicate the momentum positions at which SR-ARPES measurements were carried out. (b) Spin-resolved energy-dispersive curvess (blue and red circles for up and down spin) and the in-plane spin polarization (red circle with error bar -right axis) at momentum position #1. (c) Same as b for momentum position #2. Helical spin-texture is observed in BiPd. Error bars represent the experimental uncertainties (s.d.) in determining the spin polarization. Spin-resolved measurements were carried out with photon energy of 25 eV and temperature of 20 K at Beamline-9B of the Hiroshima Synchrotron Radiation Center (HiSOR), Hiroshima, Japan.

To further understand the nature of the observed Dirac-like band, we measured the in-plane spin polarization or spin-texture properties of BiPd. Spin-resolved ARPES measurements were performed on the top-cleaved surface of BiPd. Two spin-resolved energy-dispersive curves are shown at momenta of B0.1 on the opposite sides of the Dirac-like dispersion. The blue dash lines

on the ARPES dispersion map shown in Fig. 5 mark the approximate momentum positions for the spin-resolved energy distribution curves. The spin-resolved energy-dispersive curves and corresponding net spin polarization are shown in Fig. 5b,c for momentum positions #1 and #2. The obtained spin data presented in Fig. 5b,c show observable net spin polarization.

NATURE COMMUNICATIONS | 7:13315 | DOI: 10.1038/ncomms13315 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 5

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms13315

Importantly, the observed spin polarization at one branch of the cone is opposite with that of another cone, thus conrming the spin-momentum locking behaviour. Spin-resolved ARPES data reveals a characteristic helical spin-texture, suggesting the possible topological origin of the states we observed. However, we caution that two-band crossings at a time-reversal invariant momentum should have an opposite value of the spin polarization independent of any topological invariant.

DiscussionThe experimental realization of the topological insulator phase in a NCS crystal structure is an object of intense research. Such a system may be utilized in testing several proposed exotic phenomena, such as crystalline-surface-dependent topological electronic states, pyroelectricity and natural topological pn junctions29. Recently, the rst-principles calculations predicted III-Bi to be an inversion asymmetric topological insulator with large bandgap possessing intrinsic topologically protected edge states and forming quantum spin Hall systems29, but these have not yet been realized experimentally. At the same time, the proposal of topological insulating nature in an inversion asymmetric compound BiTeCl still remains under debate30,31. Furthermore, the small bulk bandgaps of the realized inversion asymmetric topological insulators severely limit the manipulation and control of the topological surface states.

We note that BiPd consists of many bands near the Dirac-like states, which makes this system very complicated with no global bandgap in the normal state. The topology in BiPd is not well dened such as in Bi2Se3. In ARPES experiments, the samples were cleaved and measured in an ultrahigh vacuum environment that kept its surface clean during the measurements. The surface cleaving generates surface potential due to the surface charges, which may develop a large effective pressure along the b (perpendicular to the surface, see Fig. 1a) direction and thus drive the crystal (possibly the several top layers) into the topological insulator phase. Our calculations on strained BiPd show the bulk bandgap opening, (Supplementary Fig. 7; Supplementary Note 4). This study will enable further discussion on the theoretical origin of the topological order in NCS materials.

BiPd is a SC below TcB3.7 K. Since the spin-polarized surface state is located at high-binding energy, it negates the topological superconductivity behaviour on the surface at its native Fermi level. However, by electrical gating or surface deposition, the Fermi level can be tuned near the Dirac surface state, which provides an opportunity to realize the topological superconductivity in this noncentrosymmetic material.

Methods

Crystal growth and characterization. Single crystals of BiPd were grown by a modied Bridgman method as described elsewhere4. The crystals were characterized by means of X-ray diffraction, energy-dispersive X-ray spectroscopy, magnetic susceptibility, electrical resistivity and heat capacity measurements, using standard commercial equipment.

Electronic structure measurements. Synchrotron-based ARPES measurements of the electronic structure were performed at the Advanced Light Source (ALS), Berkeley at Beamline 10.0.1 and Stanford Synchrotron Radiation Lightsource (SSRL) at Beamline 54 both equipped with a high efciency R4000 electron analyser. The energy resolution was set to be better than 20 meV for the measurements with the synchrotron beamline. The angular resolution was set to be better than 0.2 for all synchrotron measurements. Samples were cleaved in situ and measured at 1080 K in a vacuum better than 10 10 torr. They were found to be very stable and without degradation for the typical measurement period of 20 h.

Spin-resolved ARPES measurements. Spin-resolved ARPES measurements were performed at the ESPRESSO end station installed at Beamline-9B of the Hiroshima Synchrotron Radiation Center (HiSOR), Hiroshima, Japan, equipped with state-of-the-art very low-energy electron diffraction spin detectors utilizing

preoxidized Fe(001)-p(1 1)-O targets32,33. The two spin detectors are placed at

an angle of 90 and are directly attached to a VG-Scienta R4000 hemispheric analyser, enabling simultaneous spin-resolved ARPES measurements for all three spin components. We also measured the x and z component of the spin that shows negligible spin polarization.

First-principles calculations. The rst-principles calculations were based on the generalized gradient approximation 34 using the projector augmented-wave method35, as implemented in the VASP package36,37. The experimental crystallographic structure was used38 for the calculations. The SOC was included self-consistently in the electronic structure calculations with a 6 4 5

Monkhorst-Pack k-mesh. To simulate surface effects, we used 1 5 1 supercell

for the (010) surface, with a vacuum thickness larger than 20 .

Data availability. All relevant data are available from the corresponding author upon request.

References

1. Bauer, E. & Sigrist, M. Non-Centrosymmetric Superconductors: Introduction and Overview, vol. 847 Lecture Notes in Physics (Springer, 2012).

2. Bauer, E. et al. Heavy fermion superconductivity and magnetic order in noncentrosymmetric CePt3Si. Phys. Rev. Lett. 92, 027003 (2004).

3. Frigeri, P. A., Agterberg, D. F., Koga, A. & Sigrist, M. Superconductivity without Inversion Symmetry: mnSi versus CePt3Si. Phys. Rev. Lett. 92, 097001 (2004).

4. Joshi, B., Thamizhavel,, A. & Ramakrishnan, S. Superconductivity in noncentrosymmetric BiPd. Phys. Rev. B 84, 064518 (2011).

5. Jiao, L. et al. Anisotropic superconductivity in noncentrosymmetric BiPd. Phys. Rev. B 89, 060507 (R) (2014).

6. Sun, Z. et al. Dirac surface states and nature of superconductivity in noncentrosymmetric BiPd. Nat. Commun. 6, 6633 (2015).

7. Hasan, M. Z. & Kane, C. L. Colloquium: topological insulators. Rev. Mod. Phys.

82, 30453067 (2010).

8. Qi, X.-L. & Zhang, S.-C. Topological insulators and superconductors. Rev. Mod. Phys. 83, 10571110 (2011).

9. Hasan, M. Z., Xu, S.-Y. & Neupane, M. in Topological Insulators: Fundamentals and Perspectives (eds Ortmann, F., Roche, S. & Valenzuela, S. O.) (John Wiley & Sons, 2015).

10. Neupane, M. et al. Observation of quantum-tunnelling-modulated spin texture in ultrathin topological insulator Bi2Se3 lms. Nat. Commun. 5, 3841 (2014).

11. Hasan, M. Z., Xu, S.-Y. & Bian, G. Topological insulators and superconductors. Phys. Scr. T164, 014001 (2015).

12. Neupane, M. et al. Observation of a topological 3D Dirac semimetal phase in high-mobility Cd2As3. Nat. Commun. 5, 3786 (2014).

13. Bansil, A., Lin, S. & Das, T. Colloquium: topological band theory. Rev. Mod. Phys. 88, 021004 (2016).

14. Neupane, M. et al. Non-KIondo-like electronic structure in the correlated rare-earth hexaboride YbB6. Phys. Rev. Lett. 114, 016403 (2015).

15. Xu, S.-Y. et al. Momentum-space imaging of Cooper pairing in a half-Dirac-gas topological superconductor. Nat. Phys. 10, 943950 (2014).

16. Neupane, M. et al. Surface electronic structure of the topological Kondo-insulator candidate correlated electron system SmB6. Nat. Commun. 4, 2991 (2013).

17. Xia, Y. et al. Observation of a large-gap topological-insulator class with a single Dirac cone on the surface. Nat. Phys. 5, 398402 (2009).

18. Neupane, M. et al. Gigantic surface lifetime of an intrinsic topological insulator. Phys. Rev. Lett. 115, 116801 (2015).

19. Chadov, S. et al. Tunable multifunctional topological insulators in ternary Heusler compounds. Nat. Mater. 9, 541545 (2010).

20. Lin, H. et al. Half-Heusler ternary compounds as new multifunctional experimental platforms for topological quantum phenomena. Nat. Mater. 9, 546549 (2010).

21. Schnyder, A. P., Ryu, S., Furusaki, A. & Ludwig, A. W. W. Classication of topological insulators and superconductors in three spatial dimensions. Phys. Rev. B 78, 195125 (2008).

22. Lu, C.-K. & Yip, S. Spin current in topologically trivial and nontrivial noncentrosymmetric superconductors. Phys. Rev. B 82, 104501 (2010).

23. Sato, M. & Fujimoto, S. Topological phases of noncentrosymmetric superconductors: Edge states, majorana fermions, and non-abelian statistics. Phys. Rev. B 79, 094504 (2009).

24. Hosur, P., Ghaemi, P., Mong, R. S. K. & Vishwanath, A. Majorana modes at the ends of superconductor vortices in doped topological insulators. Phys. Rev. Lett. 107, 097001 (2011).

25. Xu, S.-Y. et al. Discovery of a Weyl Fermion semimetal and topological Fermi arcs. Science 349, 613617 (2015).

26. Bian, G. et al. Topological nodal-line fermions in spin-orbit metal PbTaSe2.

Nat. Commun. 7, 10556 (2016).

6 NATURE COMMUNICATIONS | 7:13315 | DOI: 10.1038/ncomms13315 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms13315 ARTICLE

27. Mondal, M. et al. Andreev bound state and multiple energy gaps in the noncentrosymmetric superconductor BiPd. Phys. Rev. B 86, 094520 (2012).

28. Matano, K. et al. NMR and NQR studies on non-centrosymmetric superconductors Re7B3, LaBiPt, and BiPd. J. Phys. Soc. Jpn. 82, 084711 (2013).

29. Ma, Y. et al. Two-dimensional inversion asymmetric topological insulators in functionalized III-Bi bilayers. Phys. Rev. B 91, 235306 (2015).

30. Landolt, G. et al. Bulk and surface Rashba splitting in single termination BiTeCl. New J. Phys. 15, 085022 (2013).

31. Chen, Y. L. et al. Discovery of a single topological Dirac fermion inthe strong inversion asymmetric compound BiTeCl. Nat. Phys. 9, 704708 (2013).

32. Okuda, T. et al. Efcient spin resolved spectroscopy observation machine at Hiroshima Synchrotron Radiation Center. Rev. Sci. Instrum 82, 103302 (2011).

33. Okuda, T. et al. A double VLEED spin detector for high-resolution three dimensional spin vectorial analysis of anisotropic Rashba spin splitting. J. Ele. Spect. Rel. Phen. 201, 2329 (2015).

34. Perdew, P., Burke, K. & Ernzerhof, M. Generalized gradient approximation made simple. Phys. Rev. Lett. 77, 3865 (1996).

35. Blochl, P. E. Projector augmented-wave method. Phys. Rev. B 50, 17953 (1994).36. Kresse, G. & Hafner, J. Ab initio molecular dynamics for open-shell transition metals. Phys. Rev. B 48, 13115 (1993).

37. Kress, G. & Furthmuller, J. Efcient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Comput. Mater. Sci. 6, 1550 (1996).

38. Bhatt, Y. C. et al. Kristallstruktur von PdBi.r. J. Less-Common. Met. 64, 1724 (1979).

Acknowledgements

M.N. is supported by the start-up fund from University of Central Florida and LANL LDRD Program. T.D. was supported by NSF IR/D program. D.K. was supported by the National Science Centre (Poland) under research grant 2015/18/A/ST3/00057. J.-X. Z. is supported by the Center for Integrated Nanotechnologies, a U.S. DOE Ofce of Basic Energy Sciences user facility, in partnership with the LANL Institutional Computing Program for computational resources. The work at Princeton and synchrotron X-ray-based measurements are supported by the Ofce of Basic Energy Sciences, US Department of Energy (DOE) grant no. DE-FG-02-40105ER46200. The work at Northeastern University is supported by the DOE, Ofce of Science, Basic Energy Sciences grant number DE-FG02-07ER46352, and beneted from Northeastern Universitys Advanced Scientic Computation Center (ASCC) and the NERSC supercomputing center through DOE grant number DE-AC02-05CH11231. H.L. acknowledges the Singapore National Research Foundation for the support under NRF

award no. NRF-NRFF2013-03. T.R.C. and H.T.J. were supported by the National Science Council, Taiwan. We also thank NCHC, CINC-NTU and NCTS, Taiwan, for technical support. The measurements at HiSOR were performed with the approval of the Proposal Assessing Committee of HSRC (Proposal No. 15-A-66). We thank Sung-Kwan Mo and Makoto Hashimoto for beamline assistance at the LBNL and the SSRL.

Authors contributions

M.N. designed the study; D.K., R.S. and F.C. synthesized the samples and performed the transport characterization; M.N. performed the spectroscopy measurements with the help of N.A., S.-Y.X., I.B., D.S., M.M.H., K.D., N.D., K.M., T.O., M.Z.H., T.D.; J.-X.Z.,T.-R.C., H.-T.J and A.B. performed the calculations; M.N. wrote the manuscript and discussed with all the co-authors. M.N. was responsible for the overall research direction, planning and integration among different research units.

Additional information

Supplementary Information accompanies this paper at http://www.nature.com/naturecommunications

Web End =http://www.nature.com/ http://www.nature.com/naturecommunications

Web End =naturecommunications

Competing nancial interests: The authors declare no competing nancial interests.

Reprints and permission information is available online at http://npg.nature.com/reprintsandpermissions/

Web End =http://npg.nature.com/ http://npg.nature.com/reprintsandpermissions/

Web End =reprintsandpermissions/

How to cite this article: Neupane, M. et al. Observation of the spin-polarized surface state in a noncentrosymmetric superconductor BiPd. Nat. Commun. 7, 13315 doi: 10.1038/ncomms13315 (2016).

Publishers note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional afliations.

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the articles Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/

Web End =http://creativecommons.org/licenses/by/4.0/

r The Author(s) 2016

NATURE COMMUNICATIONS | 7:13315 | DOI: 10.1038/ncomms13315 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 7