p-wave triggered superconductivity in

Full text

Turn on search term navigation

ARTICLE

Received 9 Jun 2016 | Accepted 21 Nov 2016 | Published 19 Jan 2017

A. Di Bernardo1, O. Millo2, M. Barbone3, H. Alpern2, Y. Kalcheim2, U. Sassi3,w, A.K. Ott3, D. De Fazio3, D. Yoon3,M. Amado1, A.C. Ferrari3, J. Linder4 & J.W.A. Robinson1

Electron pairing in the vast majority of superconductors follows the Bardeen CooperSchrieffer theory of superconductivity, which describes the condensation of electrons into pairs with antiparallel spins in a singlet state with an s-wave symmetry. Unconventional superconductivity was predicted in single-layer graphene (SLG), with the electrons pairing with a p-wave or chiral d-wave symmetry, depending on the position of the Fermi energy with respect to the Dirac point. By placing SLG on an electron-doped (non-chiral) d-wave superconductor and performing local scanning tunnelling microscopy and spectroscopy, here we show evidence for a p-wave triggered superconducting density of states in SLG. The realization of unconventional superconductivity in SLG offers an exciting new route for the development of p-wave superconductivity using two-dimensional materials with transition temperatures above 4.2 K.

DOI: 10.1038/ncomms14024 OPEN

p-wave triggered superconductivity in single-layer graphene on an electron-doped oxide superconductor

1 Department of Materials Science and Metallurgy, University of Cambridge, 27 Charles Babbage Road, Cambridge CB3 0FS, UK. 2 Racah Institute of Physics and the Hebrew University Center for Nanoscience and Nanotechnology, The Hebrew University of Jerusalem, Jerusalem 91904, Israel. 3 Cambridge Graphene Centre, University of Cambridge, Cambridge CB3 0FA, UK. 4 Department of Physics, Norwegian University of Science and Technology, N-7491 Trondheim, Norway. w Present address: Nokia Bell Labs, Broers Building, Cambridge CB3 0FA, UK. Correspondence and requests for materials should be addressed to O.M.

(email: mailto:[email protected]

Web End [email protected] ) or to J.W.A.R. (email: mailto:[email protected]

Web End [email protected] ).

NATURE COMMUNICATIONS | 8:14024 | DOI: 10.1038/ncomms14024 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 1

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms14024

Electron pairing in the vast majority of superconductors follows the Bardeen, Cooper and Schrieffer (BCS) theory of superconductivity1, which describes the condensation of

electron pairs with opposite spins in a singlet state with an isotropic order parameter (s-wave symmetry). Superconductivity in a number of compounds is predicted to be unconventional in that the electrons pair in a triplet state with parallel spins and anisotropic p-wave symmetry, most notably in UPt3 (ref. 2)

and Sr2RuO4 (SRO)3,4. Although the investigation of compounds that are thought to be bulk p-wave superconductors has been intensive over the past two decades5,6 there is still much debate over the nature of the superconductivity in these compounds. Theoretically, a p-wave superconducting state is highly sensitive to electron scattering which has meant that, with exception to ref. 7, all experiments to date on potential p-wave superconductors have focused on single-crystal samples. A further complication relates to the typically low transition temperature (Tc) of

these superconductors (SRO has a Tc below 1.5 K)5, which, in conjunction with the requirement of extreme purity, makes the fabrication of p-wave devices challenging. More recently, tunnelling studies8,9 have suggested that p-wave superconductivity can emerge on the surface of topological insulators coupled to s-wave superconductors, but the manipulation of these states for device applications is still an open issue. Triplet (s-wave) states have also been reported in spectroscopic experiments involving magnetically inhomogeneous s-wave superconductor/ferromagnet hybrids1012.

Electrons in single-layer graphene (SLG) are predicted to condense to a superconducting state, either intrinsically by doping1320 or by placing SLG on a superconductor with a BCS or a non-BCS pairing symmetry21,22. The resulting symmetry depends on the position of the Fermi energy (EF) with respect to the Dirac point. In particular, for EF shifts up to 1 eV, a p-wave15,16 state is predicted. As the doping approaches the van Hove singularity (EFB3 eV; ref. 19), a singlet chiral d-wave and triplet f-wave symmetry are also predicted17,20, whereas ref.

18 found dominant chiral d-wave superconductivity near van Hove doping and argued that weak coupling superconductivity for doping levels between half-lling and the van Hove density is of KohnLuttinger type and likely to be f-wave pairing for disconnected Fermi pockets. Reference 16 predicted that a nonchiral p-wave symmetry is favoured for small nearest-neighbour repulsion V (o1.1 eV), small onsite interaction U (B8.4 eV) or large doping (above 10%), whereas the chiral p-wave state occurs as U or V are increased or the doping level diminishes with respect to the aforementioned values (in pure SLG at half-lling U is B9.3 eV and V is B5.5 eV; ref. 16). At low density (B20%)

and including next-nearest neighbour hopping, a chiral p-wave state can emerge15. Moreover, the possibility of spin-triplet s-wave pairing has been considered in bilayer graphene23.

Although intrinsic superconductivity in SLG has not been observed24, superconductivity has been induced by doping SLG with Li adatoms25, intercalating SLG sheets with Ca (ref. 26) or by placing SLG on a superconductor27. In the latter case, the intrinsic pairing potential for p- or chiral d-wave superconductivity can be enhanced theoretically21,22 to the point that a full transition to a superconducting state is triggered and manifested in the SLG superconducting density of states (DoS). Achieving p-wave or chiral d-wave superconductivity in SLG via a proximity effect would enable the fabrication of devices and allow the full investigation of how these symmetry states can be utilized in cryogenic technology related to spin transport with low dissipation. In addition, it would be attractive for applications to achieve p-wave superconductivity in SLG above 4.2 K.

The standard technique to probe electron pair symmetry involves measuring the voltage (V) bias dependence of the

non-linear differential conductance (dI/dV V) near the

superconducting gap edge, which is proportional to the quasiparticle DoS in the tunnelling limit28. For a conventional superconductor with s-wave (zero spin) pairing, dI/dV diminishes below the gap edge in the tunnelling limit. In contrast, unconventional p- and chiral d-wave symmetries result in zero-energy quasiparticle excitations and surface bound states29, which manifest themselves in the tunnelling spectra as subgap structures in the DoS, including V-shaped gaps or strong zero-bias conductance peaks (ZBCPs)3033.

Tonnoir et al.27 locally probed the superconducting DoS in SLG on the s-wave superconductor Re by scanning tunnelling microscopy (STM). They inferred induced superconductivity in SLG from the observation of a gapped DoS that matched the underlying layer of Re (s-wave). The absence of a subgap structure and, therefore, unconventional superconductivity, may indicate a modication of the SLG band structure34,35 due to the high carrier density of Re (neB4.5 1023 cm 3) resulting in

signicant charge transfer.

An alternative approach to ref. 27 is to place SLG on an oxide cuprate superconductor as the lower carrier concentrations of these materials (B1020 cm 3; ref. 36) will reduce unwanted effects due to doping, such as band structure modications.

However, most oxide-based superconductors have a complex DoS due to their anisotropic superconducting order parameters3032, which will necessarily complicate the analysis of the SLG DoS spectra, unless care is taken.

Pr2 xCexCuO4 (PCCO) shares similarities with other

high-temperature oxide superconducting materials, such as YBa2Cu3O7 (YBCO). In particular, it has a d-wave symmetry, albeit with a much longer coherence length (xs) of B30 nm (ref. 37) compared with that of YBCO of B2 nm (ref. 37). Tunnelling studies on hole-doped (xo0.13) PCCO (ref. 38)

reveal anisotropic superconducting DoS consistent with d-wave symmetry with enhancements of the conductance at zero energy as observed on hole-doped YBCO (ref. 39). In contrast, for electron-doped PCCO (x40.13), which is the case considered in our experiment where x 0.15 in PCCO, the superconducting

DoS is found to be isotropic37,40 with a gapped DoS for all tunnelling directions. This is due to the electron-doped PCCO being in the dirty limit40 meaning xs4le, with le being the electron mean free path.

Here we report a scanning tunnelling spectroscopy (STS) investigation of SLG placed on the electron-doped cuprate superconductor PCCO to probe the possibility of inducing unconventional superconductivity in SLG such as p-wave15,16 via a superconductor proximity effect. Our results reveal modications of the superconducting DoS in SLG compared with the underlying PCCO at 4 K: a mixture of V-shaped gaps, strong ZBCPs and split ZBCPs are observed, depending on the position of the STM tip. A theoretical analysis indicates that the subgap structures match p-wave superconductivity with the different symmetry components of the p-wave state depending on the surface orientation of the underlying PCCO.

ResultsPCCO characterization and Raman spectroscopy of SLG/ PCCO. The PCCO (200-nm-thick) is grown by pulsed laser deposition on (001) oriented SrTiO3 (STO) (see Methods) and has a near-bulk superconducting transition of 20.5 K (see Fig. 1a and Supplementary Fig. 1a). High-angle X-ray diffraction performed on control PCCO lms grown in the same conditions as those measured in this study by STM (Fig. 1b and Supplementary Fig. 1b) reveals poor c-axis texturing with strong peaks from the (110) family of planes. In particular, the rocking curves of the

2 NATURE COMMUNICATIONS | 8:14024 | DOI: 10.1038/ncomms14024 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms14024 ARTICLE

c d

STO(001)

STO(002)

1011

1010

PCCO (008)

STO(003)

Intensity (Counts pers)

109

PCCO (010)

Intensity

(Counts per s)

108

PCCO(220)

Intensity

(Counts per s)

Intensity

(Countspers)

107

106

e f

105

(Degrees)

104

200 nm

103

102

101

10 20 30 40 50 60 70 80

Intensity (Counts per s)

2 (Degrees)

(Degrees) (Degrees)

Intensity (arb. units)

293 K G

150 K

4 K

PCCO/STO substrate

SLG/Cu

250

500 750 1,000

1,400

1,600

2,600

2,800

1,400

1,600

2,500

2,750

Raman shift (cm1)

Raman shift (cm1) Raman shift (cm1)

Figure 1 | Characterization of SLG on electron-doped PCCO. (a) Schematic of SLG (hexagonal lattice) on PCCO/STO. (b) High-angle X-ray diffraction data on four control lms of 5 mm 5 mm PCCO/STO deposited using identical growth conditions (curves are vertically offset for clarity). (cf) Rocking curves

(omega scans) of the (006) PCCO diffraction peaks from the same samples investigated in (b) with matching colours showing a FWHM of 0.85 (c), 0.42 (d),0.51 (e) and 0.61 (f). (g) STM topography map showing SLG on PCCO/STO measured at 4.2K (scale bar in g has a length of 0.5 nm). (h) Raman spectra at 514.5 nm of PCCO as-grown on STO (grey) and following transfer of SLG (red) at 293 K. (i) Raman spectra at 514.5 nm of SLG as-grown on Cu (grey) and following transfer onto PCCO/STO (red) measured at 293 K. The PCCO/STO background is subtracted to allow identication of D, G and 2D peaks. (j) Raman spectra at 514.5nm of SLG on PCCO/STO, after background subtraction, at 293 K (red), 150 K (blue) and 4.2 K (light blue).

(006) diffraction peak of the PCCO lms show a spread in the full width at half maximum (FWHM) in the range 0.420.85 (Fig. 1cf), which are large values for perfect (001) texturing. These values are comparable to the FWHMs of the diffraction peaks of the (110) family of planes (with differences typically of B0.3 for neighbouring peaks belonging to the two families of planes) as shown in Supplementary Fig. 1cf, which conrms that the PCCO surface is a mixture of (100) and (110) planes. This is in agreement with previous studies38,40 which have shown that secondary (110) phases would always form during the PCCO growth.

SLG is grown on Cu by chemical vapour deposition (CVD), as described in ref. 41, and transferred onto PCCO/STO, following the procedure in ref. 42. The analysis of the Raman spectrum acquired on as-grown graphene on Cu (Fig. 1i, grey curve) shows that the 2D peak can be tted with a single Lorentzian with a position, Pos (2D), B2681 cm 1 and FWHM(2D) B37 cm 1, indicating SLG (ref. 43). The G peak position, Pos(G), and

FWHM(G) are 1,587 cm 1 and 20 cm 1, respectively; and the intensity ratio, I(2D)/I(G), and area ratio, A(2D)/A(G), are 2.9 and5.6, respectively. These indicate a doping of B200 meV for SLG on Cu at room temperature4446. The spectrum shows a small D peak with I(D)/I(G) 0.15, corresponding to a defect concentration of

nDB3.8 1010 cm 2 (refs 45,47). A compressive biaxial strain of

0.07% is also estimated from the Raman analysis. To compare the SLG quality before and after transfer, the background signal of PCCO on STO is measured under identical conditions (see Fig. 1h and Methods section) and subtracted after normalization to the intensity of the Raman peak of STO at B300 cm 1. After transfer,

Pos(2D) and FWHM(2D) are B2681 cm 1 and 37 cm 1, Pos(G)

and FWHM(G) are 1,581 cm 1 and 17 cm 1, while I(2D)/I(G) and A(2D)/A(G) are 3.4 and 7.1, respectively (Fig. 1i, red curve).

This implies a reduction in doping to o100 meV (refs 4446), with I(D)/I(G) B0.09, corresponding to nDB2.2 1010 cm 2.

Similar values of defect concentration and strain for as-grown SLG on Cu and after its transfer onto PCCO imply homogeneous SLG on PCCO.

To assess doping and sample quality in the same temperature range used for the superconductivity studies, Raman measurements are also performed between 4.2 K and room temperature (Fig. 1j and Supplementary Fig. 2a). The low-temperature Raman data indicate that doping remains below 100 meV and that there are no structural changes compared with room temperature.

Local density of states measurements on SLG/PCCO. Using STM we locally measure dI/dV V spectra of SLG on PCCO at

4.2 K and correlate these to surface topography (Fig. 1g). The most predominant superconducting-related spectra show either V-shaped gaps (Fig. 2a) or a subgap structure, including ZBPCs and split ZBCPs (Fig. 2b,c). Additional supporting data are shown in Supplementary Fig. 3. V-shaped gaps are observed in about 45% of the scans, while in all the other areas we observe either ZBCPs (B30%) or split peaks (B25%). In the normal state none of these spectral features is observed, which rules out the possibility that these are due to electronic inhomogeneity in the sample, since this should give rise to the same features independently of temperature. In some areas, in all of which the SLG local topography cannot be properly resolved, single-electron

NATURE COMMUNICATIONS | 8:14024 | DOI: 10.1038/ncomms14024 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 3

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms14024

1.8

3.0

2.0

1.5

1.8

2.5

1.6

dI/dV(arb. units)

1.2

2.0

1.4

0.9

1.5

1.2

0.6

1.0

0.3

Voltage (mV)

Figure 2 | STM differential conductance versus bias-voltage spectra on SLG/PCCO/STO at 4.2 K. (ac) Typical proximity-induced V-shaped gaps (a), ZBCPs (b) and split ZBCPs (c). Different colours in ac are used to distinguish between spectra recorded in different sample areas. Insets in ac show the typical topography of a sample area where the corresponding spectra in ac are obtained (the scale bars in the insets have a length of 0.5 nm).

tunnelling spectra are measured in the superconducting state (Supplementary Fig. 4); such features persist above Tc (at least up to 50 K), as also shown in Supplementary Fig. 4.

The evolution of superconducting-related spectral features is also studied as the sample is warmed up above its superconducting transition. All subgap features including ZBCPs and split ZBCPs are suppressed as the sample is warmed up (Supplementary Fig. 5) giving spectra that are either consistent with lightly doped SLG with EF within 100 meV from the Dirac point48, or to structureless (at) spectra, or single-electron tunnelling effects. The latter two appear in o10% of the total scanned area, where no superconducting-related features are observed below Tc and the STM topographic images do not exhibit a clear SLG structure. This shows that the proximity takes place only in regions that structurally and electronically conform to the SLG behaviour. Whilst proving that the V-shaped gaps, ZBCPs or split peaks are related to superconductivity, the absence of such structures (particularly conductance peaks) above the superconducting transition temperature, rules out spurious effects from defects (for example magnetic impurities) or structural inhomogeneity that may cause Kondo scattering and thus enhancements of the DoS, as such features would also be present in the normal state. Further support for this conclusion comes from the fact that ZBCPs are found (below Tc) only in regions where the STM images show clear SLG topography and not in defected regions. We also measure the evolution of ZBCPs in an applied out-of-plane magnetic eld (Supplementary Fig. 6), and their magnitude is always found to decrease, with no splitting ever observed, which is also inconsistent with a Kondo effect. This effect of magnetic eld is similar to those observed in tunnelling spectroscopy studies of proximity-induced p-wave superconductivity in Bi2Se3 (ref. 9) and odd-parity topological superconductivity in CuxBi2Se3 (ref. 49).

Control experiments on bare PCCO and Au/PCCO. Control samples of bare PCCO are also investigated, but no subgap structure or V-shaped gaps are observed (Fig. 3a), even at facets that expose the nodal ab-plane. The only DoS features we can observe are smeared BCS-like gaps, consistent with previous STS experiments on electron-doped PCCO (refs 37,40). We investigate the effect of substrate choice on the superconductor proximity effect by fabricating SLG/PCCO on (001) LaAlO3 (LAO) and observe similar results to STO (Supplementary Fig. 7).

The ZBCPs on SLG/PCCO might be related to the penetration of the anisotropic components of PCCO that are masked in bare

PCCO, but which may appear in SLG due to its long electron mean free path (B100 nm near the Dirac point)50 and spin diffusion length (41 mm) (ref. 51) being much longer than the coherence length in PCCO of B30 nm (ref. 37). These are important parameters to consider, since the superconducting condensate is induced in the entire graphene plane. To check this, we replace SLG with Au and fabricate Au/PCCO/STO lms (Fig. 3b), where a 10-nm-thick polycrystalline Au layer is deposited in-situ without breaking vacuum by pulsed laser deposition, and perform STS. Au is chosen since it has an electron mean free path B30 nm at room temperature52. The topography maps of Au/PCCO/STO (Fig. 3b) reveal a low surface roughness (B1 nm over a 1 mm2 area;

Supplementary Fig. 8), and the corresponding tunnelling spectra on Au show superconducting gaps with no subgap structure, thus supporting our claim that the subgap structures in SLG/PCCO are related to SLG and not the underlying PCCO. We note that the gapped spectra on Au/PCCO are shallower than those of bare PCCO (Fig. 3a,b), but qualitatively similar, suggesting that Au is fully superconducting.

Tuning the superconducting proximity effect in SLG/PCCO. To further conrm our claim that the spectral features on SLG/PCCO are related to a superconductor proximity effect in SLG and not the underlying PCCO, we also measure spectra on SLG covered with evaporated microislands (10 mm in diameter and 30 nm in height) of either Ag or Au (Fig. 3c,d). Unlike metals, such as Pd, which have a stronger binding interaction with SLG (ref. 34), the low binding energy of Au and Ag (0.04 eV per carbon atom34) results in a reduced modication of the SLG band structure34. Ag (Au) has a work function of 4.4 eV (5.2 eV) which is lower (higher) than the work function of SLG (B4.6 eV; ref. 34). Therefore, if a superconducting proximity effect occurs on SLG/PCCO, this should be stronger near a Ag microisland than a Au one, since Ag (Au) acts as a donor (acceptor) of electrons in SLG, with a doping prole that can extend up to hundreds nanometre from the microisland step edge53. The tens of dI/dV V spectra acquired on these

samples show that pronounced ZBCPs are observed near and even on the Ag microislands, but no ZBCPs are ever observed near or on Au microislands, despite the well-dened SLG topography (some representative spectra are shown in Fig. 3c and in Supplementary Figs 9 and 10). This implies that the Au islands reduce the SLG electron density to the extent that the superconducting proximity effect for SLG on PCCO is suppressed. We also note that these data indicate that possible doping inhomogeneities can only affect the size of the unconventional superconducting features (even to their

4 NATURE COMMUNICATIONS | 8:14024 | DOI: 10.1038/ncomms14024 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms14024 ARTICLE

1.8

i iii

ii iv

1.5

1.2

dI/dV(arb. units)

0.9

0.6

0.3

10 8 6 4 2 0 2 4 6 8 10

Voltage (mV)

2.2

1.0

2.0

Ag(islands)/graphene/PCCO/STO Au(islands)/graphene/PCCO/STO Graphene/PCCO/STO

1.8

1.6

1.4

1.2

10 8 6 4 2 0 2 4 6 8 10

Voltage (mV)

Figure 3 | STM differential conductance versus bias-voltage spectra on control samples of PCCO/STO and Au/PCCO/STO at 4.2 K. (a) Typical STM spectra at 4.2 K for PCCO/STO with schematic of the sample structure (inset i) and STM topography (inset ii; the scale bar has a length of 0.5 nm). (b) Typical STM spectra at 4.2 K for Au/PCCO/STO (Au is 10-nm-thick) with schematic of the sample structure (inset iii) and STM topography (inset iv; the scale bar has a length of 0.5 nm). (c,d) Average ZBCP amplitudes for SLG/PCCO/STO (black shaded curve) and Au (green shaded curve) and Ag (red shaded curve) microislands on SLG/PCCO/STO (c) with corresponding schematics shown in (d) where the hexagonal lattice represents SLG.

suppression, along with the proximity effect altogether), but do not affect the unconventional nature of the proximity-induced superconductivity in SLG.

Calculation of the superconducting DoS in SLG/PCCO. We now consider the spatial variations in the superconducting DoS that we observe on SLG/PCCO (Fig. 2 and Supplementary Fig. 3). These indicate local variations in the proximity effect due to a combination of PCCO faceting and changes in the PCCO surface that is a mixture of different PCCO planes. Fig. 4 plots the theoretical DoS on SLG as a function of the STM orientation relative the graphene surface, which varies due to the PCCO faceting.

These are calculated by applying the model in refs 22,54, which predicts an effective p-wave pairing at the Dirac points in SLG on a d-wave superconductor and involves solving a tight-binding Hamiltonian for SLG assuming that the tunnelling between SLG

and the STM tip obeys an extended BlonderTinkhamKlapwijk (BTK) theory30,55. Following refs 22,54, we write the Hamiltonian in a band-basis where cq is a Fermion-operator in the conduction band and dq is a Fermion-operator in the valence band:

Xq;stEq m

qscqs tEq m d

qs dqs

h i

In equation (1) we dene the normal-state band dispersion as Eq

Pa eiqa where the sum is over nearest-neighbour vectors a, t is the hopping parameter, m is the chemical potential, while dq and uq are associated to the superconducting order parameter. The above Hamiltonian can be diagonalized and yields eigen

h i

i uq cq"d q# d q" c q#

h:c:

Xq;s

dq d q" d q# cq"c q#

NATURE COMMUNICATIONS | 8:14024 | DOI: 10.1038/ncomms14024 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 5

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms14024

1.4

b c

1.4 3.0

dI/dV(arb. units)

0.4

10 8 6 4 2 0 2 4 6 8 10

Voltage (mV) Voltage (mV) Voltage (mV)

[110]

[101]

[010]

[001]

[010]

[100]

(110) Phase

(110) Impurity Phase

Figure 4 | Theoretical proximity-induced superconducting DoS in SLG on PCCO. (ad) The plots show that the projections of the p-wave induced symmetries in SLG are dependent on the crystal orientation of the underlying PCCO. This is clearest in (d) which illustrates the relationship between the observed DoS on SLG and the plane normal to the underlying PCCO: on at regions oriented along the [001] crystallographic direction, the DoS corresponds to py-wave or antinodal dx y -wave with a strong tunnelling barrier (red curve in a); rough areas expose a mixture of normal planes oriented

between the [001] and [110] crystallographic directions, where the DoS corresponds to py-wave with a moderate tunnelling barrier (green curve in b); the DoS corresponding to the px-wave is projected on the phase oriented along the [110] crystallographic direction (blue curve in c).

values:

s 2

Since the shift in Fermi level (EF) in SLG (below 100 meV in SLG on PCCO at 4.2 K) is much larger than the superconducting order parameter (DB5 meV for PCCO), uq in equation (2) can be neglected and we are left with an effective superconductor with gap dq. It is this gap that acquires a p-wave symmetry near the

Dirac points when d-wave superconducting correlations are induced in SLG22,54. This is due to the fact that the transport properties of SLG are determined by its behaviour near the Dirac points for doping levels comparable to that in our study (o100 meV) and, as explained in refs 22,54, a d-wave symmetry in the full Brillouin zone corresponds to a p-wave symmetry in the vicinity of the Dirac points at K and K.

In this regime, where the shift in EF in SLG is much larger than D, a simplied model for tunnelling between SLG and the STM-tip based on BTK theory can be adopted to account for p-wave pairing. This is a commonly used modelling tool for tunnelling measurements with unconventional super-conductors5658. The procedure consists in setting up the wavefunctions on SLG and on the normal side of the tunnelling barrier, the strength of the latter being characterized by a dimensionless number Z. The Z parameter describes the strength of the tunnel barrier between SLG and the STM tip, that is, the barrier is effectively the vacuum in between. The Z value affects the DoS spectra as it is related to the ratio between Andreev

reections and quasiparticle tunnelling (Z 0 corresponds to a

perfect interface with zero resistance, whereas Z441 corresponds to a barrier with high tunnel resistance). The Z value can be controlled in STS via the current and bias-voltage settings (see below), which governs the tip-sample distance in the case of an ideally homogeneous sample. However, in our experiment, it is realistic to expect spatial variations in Z (for a given current-bias setting) since the STM tip approaches the graphene surface at slightly different distances at different locations, due to variations in the surface plane of the underlying PCCO and local variations in graphene-PCCO connectivity, which is hence accounted in our theoretical model.

The theoretical DoS in Fig. 4 are obtained with D 5 meV,

0.1 meV and Z varying between 0.8 and 1.7, and they show that the surface orientation of PCCO sensitively controls the DoS, due to the projection of the px and py symmetry components of the p-wave state. Importantly, in this Z regime, enhanced conductance at zero bias, and therefore the ZBCPs, cannot be related to Andreev reections unless the order parameter is unconventional, meaning it does not have an s-wave symmetry. In addition, the sharp structure of the peaks that we observe experimentally (Fig. 2 and Supplementary Fig. 3) can be accounted for only by sign-changing order-parameter symmetry, such as p-wave.

DiscussionThe modications of the superconducting DoS on SLG/PCCO/ STO (Fig. 2 and Supplementary Fig. 3) compared with

tEq

u2q

d2q

6 NATURE COMMUNICATIONS | 8:14024 | DOI: 10.1038/ncomms14024 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms14024 ARTICLE

PCCO/STO (Fig. 3a) or Au/PCCO (Fig. 3b), in particular the ZBCPs and split peaks, indicate the induction of p-wave15,16 or chiral d-wave pairing in SLG20. However, chiral d-wave requires dominant electron-electron repulsion in SLG at much higher doping levels (up to 3 eV; ref. 20) than here (B100 meV). In our low doping regime, p-wave pairing should theoretically dominate15,16,22 with spectral features in the DoS that match our ndings in Fig. 2.

We also exclude other unconventional pairing symmetries that may be compatible with the spectral features we observe, such as f-wave pairing as suggested in ref. 59, since this would require a buckling in the honeycomb lattice, which we do not expect for SLG. Another possibility, however, would be to consider the emergent order parameter symmetry near the Dirac points in the presence of f-wave pairing17,18 in the Brillouin zone to determine if this could also provide zero-energy states.

Our theoretical model shows that an induced px-wave state is manifested in the SLG DoS as a ZBCP (Fig. 4c), while an induced py-wave state manifests as a split-ZBCP (Fig. 4b). If tunnelling occurs along the c-axis of the PCCO (along [001]), a V-shaped gapped DoS consistent with py or antinodal dx2 y

Finally, we point out that in ref. 61 subgap features in the DoS including ZBCPs and split zero-bias peaks could be obtained via Andreev bound states, formed under a small Pb slab coupled to SLG. Such a system effectively constituted a quantum dot coupled to superconducting leads. Here, the local and random doping inhomogeneities in SLG on PCCO are not likely to induce well-dened quantum dots, where connement would be further suppressed by Klein tunnelling, and thus cannot account for the subgap features seen in our STM data.

Although our analysis is consistent with the triggering of an unconventional p-wave order parameter in SLG, we cannot determine the exact symmetry of the p-wave state. Since the DoS spectra look qualitatively different from those reported for a chiral spin-triplet p-wave state in a topological superconductor8,9, we believe our results to be consistent with an effective p-wave order parameter emerging at the Dirac points in SLG as a result of Andreev bound states induced by the proximity coupling with the d-wave pairing potential in PCCO (ref. 62). The results therefore suggest the exciting possibility of creating p-wave superconducting devices on SLG in which quantum coherent electron states, gate tunability and spin/charge could be addressed in the superconducting state.

Measurements of quasiparticle interference modulations63 could provide further insights into the pairing symmetry induced in SLG on PCCO. Our results suggest that the superconducting pairing of the p-wave state in our system is spin-singlet, which means that the Andreev bound states providing the zero-energy peak seen in the STM data can be interpreted as an odd-frequency spin-singlet odd-parity state64,65. Interfacial spinorbit coupling due to broken inversion symmetry at the PCCO/SLG interface may also in principle induce a triplet component in the system66; however, the induction of a triplet component would necessary require signicant spinorbit coupling at the SLG/PCCO interface, unlikely here due to the low atomic numbers of the elements in PCCO and graphene.

Methods

PCCO lm growth. A 200-nm-thick PCCO lm is grown by pulsed laser deposition on (001) STO using a stoichiometric target fabricated by a solid-state reaction method from high purity Pr6O11 (99.9% purity), CeO2 (99.9% purity) and

CuO (99.99% purity) powders. The growth is carried out at 820 C by ring a Lambda Physik KrF excimer laser (l 248 nm) on a rotating PCCO target with a

pulse rate of 7 Hz and energy density of 1.5 J cm 2, after introducing 230 mTorr of N2O in a ultra-high vacuum chamber (base pressure of 10 8 mbar). To obtain an optimal superconducting transition (20.5 K; Supplementary Fig. 1a), the PCCO lms are annealed in situ67 by evacuating the chamber immediately after growth and holding the substrate at 720 C for 4 min. Films are then cooled to room temperature in a vacuum of 10 5 mbar within 2.5 h.

Growth and transfer of SLG on PCCO. SLG is grown by CVD on a 35-mm-thick Cu foil loaded into a hot wall tube furnace, which is subsequently evacuated to a pressure B1 mTorr (ref. 41). The Cu foil is annealed in a hydrogen atmosphere (H2, 20 sccm) at 1,000 C for 30 min to reduce the Cu oxide component and to increase the grain size41,68. SLG growth is then initiated by adding 5 sccm CH4 to the H2 ow. After 30 min, the substrate is cooled in vacuum (1 mTorr) to room temperature and then unloaded from the reactor.

CVD is chosen over micromechanical exfoliation42,69 so as to cover large areas (5 5 mm2) of PCCO with SLG. This is benecial for a number of reasons: the

localization and alignment of the sample is simplied (the scan size of our STM is limited to 1 1 mm2, and the optical contrast of SLG on PCCO does not allow for

an easy identication of SLG); the large area allows us to fabricate large arrays of metallic islands (covering 2 2 mm2 of the sample surface); the

measurements benet from larger statistics and the presence of local defects can be avoided, thus not disturbing tunnelling-current measurements, which only require nm-size areas to be performed.

The CVD-grown SLG is then transferred onto 200-nm-thick epitaxial lms of electron-doped PCCO on (001) STO. To do so, a 500-nm-thick layer of polymethyl methacrylate (PMMA) is spin coated onto the SLG/Cu sample. The PMMA/SLG/Cu stack is subsequently immersed in an aqueous solution of ammonium persulfate (APS) to etch the Cu foil41,42. The PMMA/SLG stack is then placed in de-ionized

2 symmetry is

observed (Fig. 4a).

The theoretical analysis is consistent with the distribution in the DoS spectral features that we measure on the SLG/PCCO (B45% V-shaped gaps, B30% ZBCPs and B25% split peaks).

Although it is not possible to determine the orientation of the underlying PCCO below the STM tip during measurements, the X-ray data in Fig. 1b indicate that the surface of PCCO is a mixture of the (001) and (110) family of planes, and therefore the DoS according to our theory should be a mixture of V-shaped gaps (occurring for c-axis tunnelling), ZBCPs or split peaks (occurring for tunnelling from the ab-plane). The surface roughness of the PCCO will also expose different crystal orientations. Therefore, the spectroscopic features on SLG/PPCO will be related to a combination of PCCO crystal structure and surface roughness. Figure 4 also shows that the V-shaped DoS calculated for SLG/PCCO does not reach zero, which is consistent with our experimental observations (Fig. 2a and Supplementary Fig. 3), but also with previous tunnelling experiments on PCCO even at temperatures much lower (1.8 K) than our measurement set-up of 4.2 K (ref. 40).

Local variations in the Z parameter can determine the amplitude of subgap features (ZBCPs or split peaks) and the depth of V-shaped gaps, as shown in Supplementary Fig. 11. Although Z values can be changed intentionally by varying the set current and set bias-voltage values (before disabling the feedback loop and acquiring a dI/dV spectrum), care was taken to work with low set currents (B0.10.3 nA) and bias voltages just above the gap region (810 mV), yielding junction resistances B5 107 O, much larger than the quantum resistance

h e 2 (B25.8 kO) and thus well within the tunnelling regime. In this regime, the overall spectral features (ZBCPs, split ZBCPs or

V-shaped gaps) measured at a specic location on SLG/PCCO do not change upon varying the current and voltage set values, ruling out the possibility that these are related to single-electron tunnelling effects60. It should be noted, however, that the underlying orientation of the PCCO lm and the degree of local PCCO-SLG electrical connectivity can also inuence the tip-SLG distance and hence the Z value (for a given current-bias setting). We exclude that variations in the spectral features are due to defects in SLG, since STM images do not reveal any structural inhomogeneity in SLG and no damage after transfer (consistent with the Raman analysis), and, in particular, no defects are observed in regions where the spectra are acquired. Further, spectral anomalies disappear in the normal state.

NATURE COMMUNICATIONS | 8:14024 | DOI: 10.1038/ncomms14024 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 7

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms14024

water to rinse off acid residuals and shed off the de-ionized bath using the PCCO substrate. Finally, the PMMA/SLG/PCCO sample is dried at room temperature and placed in acetone to remove the PMMA layer, leaving SLG on PCCO.

Raman measurements. Raman spectroscopy is used to investigate the quality and uniformity of the as-grown SLG on Cu as well as to quantify the presence of defects and evaluate doping. Raman spectra are recorded at 514.5 nm using a Renishaw InVia spectrometer equipped with a Leica DM LM microscope and a 100 objective

(numerical aperture NA 0.85). Under these conditions the laser spot size is

B1 mm. The laser power is kept at B500 mW to avoid any possible sample damage.

To assess the doping and quality of SLG as a function of temperature, in particular below the superconducting critical temperature of PCCO of 20.5 K, low temperature Raman measurements are performed between room temperature and4.2 K (Supplementary Fig. 2a). These are carried out using an Oxford Instrument cryostat coupled with a Horiba LabRAM HR Evolution Raman spectrometer. An Olympus LUCPlanFL N 40 objective with an aberration correction ring is used.

This is adjusted to match the glass window thickness of the cryostat to enhance the signal-to-noise ratio.

Each sample is xed to the holder in the cryostat with vacuum grease (Apiezon N Grease) and the cryostat is pumped to 10 6 Torr before measurements, to reproduce the same vacuum conditions used in the STM. The sample is cooled to 4.2 K using liquid helium and the temperature is raised stepwise with a programmable temperature controller. For each temperature, a Raman spectrum of SLG on PCCO is recorded as well as a reference spectrum on a bare PCCO lm on (001) STO. To reveal the SLG contribution to the Raman signal, the spectrum measured on PCCO at a given temperature is subtracted from the corresponding one on SLG/PCCO.

Data availability. The data set generated and analysed during this study are available for access at http://dx.doi.org/10.17863/CAM.6228

Web End =http://dx.doi.org/10.17863/CAM.6228 .

References

1. Bardeen, J., Cooper, L. N. & Schrieffer, J. R. Theory of superconductivity. Phys. Rev. 108, 1175 (1957).

2. Tou, H. et al. Odd-parity superconductivity with parallel spin pairing in UPt3:

evidence from 195Pt Knight shift study. Phys. Rev. Lett. 77, 13741377 (1996).3. Maeno, Y. et al. Superconductivity in a layered perovskite without copper. Nature 372, 532534 (1994).

4. Ishida, K. et al. Spin-triplet superconductivity in Sr2RuO4 identied by 17O Knight shift. Nature 396, 658660 (1998).

5. Maeno, Y., Kittaka, S., Nomura, T., Yonezawa, S. & Ishida, K. I. Evaluation of spin-triplet superconductivity in Sr2RuO4. J. Phys. Soc. Jpn 81, 011009 (2012).

6. Linder, J. & Robinson, J. W. A. Superconducting spintronics. Nat. Phys 11, 307315 (2015).

7. Krockenberger, Y. et al. Growth of superconducting Sr2RuO4 thin lms. Appl. Phys. Lett. 97, 082502 (2010).

8. Xu, J. P. et al. Experimental detection of a majorana mode in the core of a magnetic vortex inside a topological insulator-superconductor Bi2Te3/NbSe2 heterostructure. Phys. Rev. Lett. 114, 017001 (2015).

9. Koren, G., Kirzhner, T., Kalcheim, Y. & Millo, O. Signature of proximity-induced px ipy triplet pairing in the doped topological insulator Bi2Se3 by the

s-wave superconductor NbN. Europhys. Lett. 103, 67010 (2013).10. Di Bernardo, A. et al. Signature of magnetic-dependent gapless odd frequency states at superconductor/ferromagnet interfaces. Nat. Commun. 6, 8053 (2015).

11. Kalcheim, Y., Millo, O., Di Bernardo, A., Pal, A. & Robinson, J. W. A. Inverse proximity effect at superconductor-ferromagnet interfaces: evidence for induced triplet pairing in the superconductor. Phys. Rev. B 92, 060501 (R) (2015).

12. Di Bernardo, A. et al. Inverse paramagnetic Meissner effect due to s-wave odd-frequency superconductivity. Phys. Rev. X 5, 041021 (2015).

13. Kopnin, N. B. & Sonin, E. B. BCS superconductivity of Dirac electrons in graphene layers. Phys. Rev. Lett. 100, 246808 (2008).

14. Uchoa, B. & Castro Neto, A. H. Superconducting states of pure and doped graphene. Phys. Rev. Lett. 98, 146801 (2007).

15. Ma, T., Yang, F., Yao, H. & Lin, H. Q. Possible triplet p ip superconductivity

in graphene at low lling. Phys. Rev. B 90, 245114 (2014).16. Faye, J. P. L., Sahebsara, P. & Snchal, D. Chiral superconductivity on the graphene lattice. Phys. Rev. B 92, 085121 (2015).

17. Kiesel, M. L., Platt, C., Hanke, W., Abanin, D. A. & Thomale, R. Competing many-body instabilities and unconventional superconductivity in graphene. Phys. Rev. B 86, 020507 (R) (2012).

18. Nandkishore, R., Thomale, R. & Chubukov, A. V. Superconductivity from weak repulsions in hexagonal lattice systems. Phys. Rev. B 89, 144501 (2014).

19. Gonzlez, J. Kohn-Luttinger superconductivity in graphene. Phys. Rev. B 78,

205431 (2008).

20. Nandkishore, R., Levitov, L. S. & Chubukov, A. V. Chiral superconductivity from repulsive interactions in doped graphene. Nat. Phys 8, 158163 (2012).

21. Black-Schaffer, A. M. & Doniach, S. Possibility of measuring intrinsic electronic correlations in graphene using a d-wave contact Josephson junction. Phys. Rev. B 81, 014517 (2010).

22. Linder, J., Black-Schaffer, A. M., Yokoyama, T., Doniach, S. & Sudb, A. Josephson current in graphene: Role of unconventional pairing symmetries. Phys. Rev. B 80, 094522 (2009).

23. Hosseini, M. V. & Zareyan, M. Model of an exotic chiral superconducting phase in a graphene bilayer. Phys. Rev. Lett. 108, 147001 (2012).

24. Kotov, V. N., Uchoa, B., Pereira, V. M., Guinea, F. & Castro Neto, A. H. Electron-electron interactions in graphene: current status and perspectives. Rev. Mod. Phys. 84, 10671125 (2012).

25. Ludbrook, B. M. et al. Evidence for superconductivity in Li-decorated monolayer graphene. Proc. Natl Acad. Sci. USA 112, 1179511799 (2015).

26. Chapman, J. et al. Superconductivity in Ca-doped graphene laminates. Sci. Rep. 6, 23254 (2016).

27. Tonnoir, C. et al. Induced superconductivity in graphene grown on Rhenium. Phys. Rev. Lett. 111, 246805 (2013).

28. Wolf, L. E. Principles of Electron Tunnelling Spectroscopy (Oxford University Press, 2012).

29. Hu, C. R. Midgap surface states as a novel signature of dxa2xb2 wave

superconductivity. Phys. Rev. Lett. 72, 15261529 (1994).30. Tanaka, Y. & Kashiwaya, S. Theory of tunnelling spectroscopy of d-wave superconductors. Phys. Rev. Lett. 74, 34513454 (1995).

31. Kashiwaya, S. & Tanaka, Y. Tunnelling effects on surface bound states in unconventional superconductors. Rep. Prog. Phys. 63, 16411724 (2000).

32. Yamashiro, M. & Tanaka, Y. Theory of tunnelling spectroscopy in superconducting Sr2RuO4. Phys. Rev. B 56, 78477850 (1997).

33. Li, X. W. Tunnelling conductance in normal metal/insulator/triplet superconductor junctions. Commun. Theor. Phys. 44, 381384 (2005).

34. Khomyakov, P. A. et al. First-principles study of the interaction and charge transfer between graphene and metals. Phys. Rev. B 79, 195425 (2009).35. Wintterlin, J. & Bocquet, M. L. Graphene on metal surfaces. Surf. Sci. 603, 18411852 (2009).

36. Li, P., Balakirev, F. F. & Greene, R. L. High-eld Hall resistivity and magnetoresistance of electron-doped Pr2-xCexCuO4-d. Phys. Rev. Lett. 99, 047003 (2007).

37. Kalcheim, Y., Millo, O., Egilmez, M., Robinson, J. W. A. & Blamire, M. G. Evidence for anisotropic triplet superconductor order parameter in half-metallic ferromagnetic La0.7Ca0.3Mn3O proximity coupled to superconducting

Pr1.85Ce0.15CuO4. Phys. Rev. B 85, 104504 (2012).

38. Qazilbash, M. M., Biswas, A., Dagan, Y., Ott, R. A. & Greene, R. L. Point-contact spectroscopy of the electron-doped cuprate superconductor Pr2 xCexCuO4: the dependence of conductance-voltage spectra on cerium

doping, barrier strength, and magnetic eld. Phys. Rev. B 68, 024502 (2003).39. Sharoni, A. et al. Scanning tunnelling spectroscopy of a-axis YBa2Cu3O7 d

lms: k-selectivity and the shape of the superconductor gap. Europhys. Lett. 62, 883889 (2003).40. Dagan, Y., Beck, R. & Greene, R. L. Dirty superconductivity in the electron-doped cuprate Pr2 xCexCuO4 d: tunnelling study. Phys. Rev. Lett. 99, 147004

(2007).41. Bae, S. et al. Roll-to-roll production of 30-inch graphene lms for transparent electrodes. Nat. Nanotechnol. 5, 574578 (2010).

42. Bonaccorso, F. et al. Production and processing of graphene and 2d crystals. Mater. Today 15, 564589 (2012).

43. Ferrari, A. C. et al. Raman spectrum of graphene and graphene layers. Phys. Rev. Lett. 97, 187401 (2006).

44. Das, A. et al. Monitoring dopants by Raman scattering in an electrochemically top-gated graphene transistor. Nat. Nanotechnol. 3, 201215 (2008).45. Bruna, M. et al. Doping dependence of the Raman spectrum of defected graphene. ACS Nano 8, 74327441 (2014).

46. Pisana, S. et al. Breakdown of the adiabatic Born-Oppenheimer approximation in graphene. Nat. Mater. 6, 198201 (2007).

47. Cancado, L. G. et al. Quantifying defects in graphene via Raman spectroscopy at different excitation energies. Nano Lett. 11, 31903196 (2011).

48. Xue, J. et al. Scanning tunnelling microscopy and spectroscopy of ultra-at graphene on hexagonal boron nitride. Nat. Mat. 10, 282285 (2011).

49. Sasaki, S. et al. Topological superconductivity in CuxBi2Se3. Phys. Rev. Lett. 107, 217001.

50. Giannazzo, F., Sonde, S., Lo Nigro, R., Rimini, E. & Ranieri, V. Mapping of the density of scattering centres limiting the electron mean free path in graphene. Nano Lett. 11, 46124618 (2011).

51. Yan, W. et al. Long-spin diffusion length in few-layer graphene akes. Phys. Rev. Lett. 117, 147201 (2016).

52. Jen, S. U., Yu, C. C., Liu, C. H. & Lee, G. Y. Piezoresistance and electrical resistivity of Pd, Au, and Cu lms. Thin Solid Films 434, 316 (2003).

53. Giovannetti, G. et al. Doping graphene with metal contacts. Phys. Rev. Lett. 101, 026803 (2008).

54. Jiang, Y., Yao, D.-X., Carlson, E. W., Chen, H.-D. & Hu, J. P. Andreev conductance in the d id wave superconducting states of graphene. Phys. Rev.

B 77, 235420 (2008).

8 NATURE COMMUNICATIONS | 8:14024 | DOI: 10.1038/ncomms14024 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms14024 ARTICLE

55. Blonder, G. E., Tinkham, M. & Klapwijk, T. Transition from metallic to tunnelling regimes in superconducting microconstrictions: excess current, charge imbalance, and supercurrent conversion. Phys. Rev. B 25, 4515 (1982).

56. Wei, J. Y. T., Yeh, N. C., Garrigus, D. F. & Strasik, M. Directional tunnelling and Andreev reection on YBa2Cu3O7-d single crystals: predominance of d-wave pairing symmetry veried with the generalized Blonder, Tinkham and

Klapwijk theory. Phys. Rev. Lett. 81, 2542 (1998).57. Daghero, D. et al. Strong-coupling d-wave superconducting in PuCoGa5 probed by point-contact spectroscopy. Nat. Commun. 3, 786 (2012).

58. Aggarwal, L. et al. Unconventional superconductivity at mesoscopic point contacts on the 3D Dirac semimetal Cd3As2. Nat. Mater. 15, 3237 (2016).

59. Zhang, L. D., Yang, F. & Yao, Y. Possible electric-eld-induced superconducting states in doped silicene. Sci. Rep. 5, 8203 (2015).

60. Hanna, A. E. & Tinkham, M. Variation of the Coulomb staircase in a two-junction system by fractional electron charge. Phys. Rev. B 44, 5919 (R) (1991).

61. Dirks, T. et al. Transport through Andreev bound states in a graphene quantum dot. Nat. Phys. 7, 386390 (2011).

62. Tanaka, Y., Sato, M. & Nagaosa, N. Symmetry and topology in superconductors odd-frequency pairing and edge states. J. Phys. Soc. Jpn 81, 011013 (2012).

63. Hamidian, M. H. et al. Atomic-scale electronic structure of the cuprate d-symmetry form factor density wave state. Nat. Phys. 12, 150156 (2016).

64. Tanaka, Y., Golubov, A., Kashiwaya, S. & Ueda, M. Anomalous Josephson effect between even- and odd-frequency superconductors. Phys. Rev. Lett. 99, 037005 (2007).

65. Tanaka, Y., Tanuma, Y. & Golubov, A. Odd-frequency pairing in normal-metal/superconductor junctions. Phys. Rev. B 76, 054522 (2007).

66. Eldenstein, V. M. Triplet superconductivity and magnetoelectric effect near the s-wave-superconductor-normal-metal interface caused by local breaking of mirror symmetry. Phys. Rev. B 67, 020505 (2003).

67. Maiser, E. et al. Pulsed-laser deposition of Pr2-xCexCuO4-y thin lms and the

effect of high-temperature post-annealing. Phys. C Supercond 297, 1522 (1998).68. Li, X. S. et al. Large-area synthesis of high-quality and uniform graphene on copper foils. Science 324, 13121314 (2009).

69. Novoselov, K. S. et al. Two-dimensional atomic crystals. PNAS 102, 1045110453 (2005).

Acknowledgements

The work was funded by the following agencies: Royal Society (Superconducting Spintronics), Leverhulme Trust (IN-2013-033), Schiff Foundation, the EPSRC (EP/N017242/1, EP/G037221/1, EP/K01711X/1, EP/K017144/1, EP/N010345/1, EP/

M507799/1, EP/L016087/1), ERC Grant Hetero2D, EU Graphene Flagship, COST Action MP-1201, MSCA-IFEF-ST No. 656485-Spin3, Outstanding Academic Fellows programme at NTNU, Research Council of Norway (205591, 216700 and 24080).

Author contributions

J.W.A.R. supervised the project and with A.D.B conceived it. A.D.B. prepared the PCCO samples and performed structural and transport measurements with help from M.A. Single-layer graphene was grown and transferred onto PCCO by U.S., D.D.F. and M.B. STM measurements were performed by O.M., H.A. and Y.K. Raman analysis was carried out by A.K.O., D.Y., M.B., A.C.F. M.B. fabricated the Ag and Au nanoislands on single-layer graphene on PCCO. The STM data were analysed by J.W.A.R., A.D.B. andO. M. with support from J.L. J.L. developed the theoretical model. The paper was written by A.D.B., O.M., A.C.F., J.L., J.W.A.R. All authors commented on the paper.

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: Di Bernardo, A. et al. p-wave triggered superconductivityin single-layer graphene on an electron-doped oxide superconductor. Nat. Commun. 8, 14024 doi: 10.1038/ncomms14024 (2017).

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) 2017

NATURE COMMUNICATIONS | 8:14024 | DOI: 10.1038/ncomms14024 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 9

Word count: 8266

Show less

Abstract

Translate

Electron pairing in the vast majority of superconductors follows the Bardeen-Cooper-Schrieffer theory of superconductivity, which describes the condensation of electrons into pairs with antiparallel spins in a singlet state with an s-wave symmetry. Unconventional superconductivity was predicted in single-layer graphene (SLG), with the electrons pairing with a p-wave or chiral d-wave symmetry, depending on the position of the Fermi energy with respect to the Dirac point. By placing SLG on an electron-doped (non-chiral) d-wave superconductor and performing local scanning tunnelling microscopy and spectroscopy, here we show evidence for a p-wave triggered superconducting density of states in SLG. The realization of unconventional superconductivity in SLG offers an exciting new route for the development of p-wave superconductivity using two-dimensional materials with transition temperatures above 4.2 K.

Details

Title

p-wave triggered superconductivity in single-layer graphene on an electron-doped oxide superconductor

Author

Di Bernardo, A; Millo, O; Barbone, M; Alpern, H; Kalcheim, Y; Sassi, U; Ott, A K; De Fazio, D; Yoon, D; Amado, M; Ferrari, A C; Linder, J; Robinson, J W A

Pages

14024

Publication year

2017

Publication date

Jan 2017

Publisher

Nature Publishing Group

e-ISSN

20411723

Source type

Scholarly Journal

Language of publication

English

DOI

https://doi.org/10.1038/ncomms14024

ProQuest document ID

1860303315

p-wave triggered superconductivity in single-layer graphene on an electron-doped oxide superconductor

Jump to:

Full text

Abstract

Details

Suggested sources