Nodeless superconductivity in a quasi-two-dimensional superconductor AuTe2Se4/3*

Project supported by the Key Basic Research Program of China (Grant Nos. 2015CB921401 and 2016YFA0300503), the National Natural Science Foundation of China (Grant Nos. 11422429 and 11421404), China Postdoctoral Science Foundation (Grant No. 2016T90332), the Program for Professor of Special Appointment (Eastern Scholar) at Shanghai Institutions of Higher Learning, and STCSM of China (Grant No. 15XD1500200), and the Strategic Priority Research Program (B) of the Chinese Academy of Sciences (Grant No. XDB04040200).

Jia Xiao-Yu1, Yu Yun-Jie1, Chen Xu2, Guo Jian-Gang2, Ying Tian-Ping1, He Lan-Po1, Chen Xiao-Long2, 3, 4, Li Shi-Yan1, 5, †
State Key Laboratory of Surface Physics, Department of Physics, and Laboratory of Advanced Materials, Fudan University, Shanghai 200433, China
Beijing National Laboratory for Condensed Matter Physics, Institute of Physics, Chinese Academy of Sciences, Beijing 100190, China
School of Physical Sciences, University of Chinese Academy of Sciences, Beijing 100049, China
Collaborative Innovation Center of Quantum Matter, Beijing 100084, China
Collaborative Innovation Center of Advanced Microstructures, Nanjing 210093, China

 

† Corresponding author. E-mail: shiyan_li@fudan.edu.cn

Project supported by the Key Basic Research Program of China (Grant Nos. 2015CB921401 and 2016YFA0300503), the National Natural Science Foundation of China (Grant Nos. 11422429 and 11421404), China Postdoctoral Science Foundation (Grant No. 2016T90332), the Program for Professor of Special Appointment (Eastern Scholar) at Shanghai Institutions of Higher Learning, and STCSM of China (Grant No. 15XD1500200), and the Strategic Priority Research Program (B) of the Chinese Academy of Sciences (Grant No. XDB04040200).

Abstract

We performed ultra-low temperature thermal conductivity measurements on the single crystal of a new gold-based quasi-two-dimensional superconductor AuTe2Se4/3, which has a superconducting transition temperature Tc = 2.70 K. A negligible residual linear term κ0/T in zero magnetic field is observed, which suggests fully gapped superconducting state. Furthermore, the field dependence of κ0/T is similar to that of the multi-band s-wave superconductor BaFe1.9Ni0.1As2 at low field. These results reveal multiple nodeless superconducting gaps in this interesting quasi-two-dimensional superconductor with Berezinsky–Kosterlitz–Thouless topological transition.

1. Introduction

Two-dimensional (2D) superconductivity has been extensively studied in recent years to explore exotic quantum phenomena, such as high-temperature superconductivity and various quantum phase transitions.[13] There are several ways to realize the 2D superconductivity as described below.

(i) It can be realized in ultra-thin films grown by molecular beam epitaxy method, e.g., in monolayer FeSe/SrTiO3 (Tc = 65 K),[1] ultra-thin Ga films ( ),[2] and highly crystalline Pb thin films (Tc = 1.5 K).[4,5]

(ii) The 2D superconductivity can be found in some particular interfaces or surfaces, e.g., the LaAlO3/SrTiO3 interface (Tc = 0.2 K),[6] the interface of an insulator and a metal La2CuO4/La1.55Sr0.45CuO4 (Tc exceeds 50 K),[7] and the surface of MoS2 thin flake by ionic liquid gating (the maximum Tc = 10.8 K).[8,9]

(iii) By vertically putting two graphene sheets together with certain twisted angles, both a Mott-like insulating state and a superconducting state (Tc ≈ 1.7 K) arise from the electrons localized in these 2D superlattices.[10,11]

(iv) An relatively easier way is to mechanically exfoliate thin flakes directly from bulk van der Waals (vdW) quasi-2D superconductors. Bi2Sr2CaCu2O8+x (Bi-2212) and NbSe2 have been successfully cleaved to atomic thin sheets with different thicknesses, in which superconductivity has been observed (Tc = 82 K for half-unit-cell Bi-2212 and Tc = 3.1 K for monolayer NbSe2, respectively).[1215] In this context, it is important to find more attractive vdW quasi-2D superconductors.

Very recently, a new gold-based quasi-2D superconductor AuTe2Se4/3 has been discovered.[16] Despite that metastable Au1−xTex (0.6 < x < 0.85) compounds are superconductors with Tc between 1.5 K and 3.0 K,[17,18] calaverite AuTe2 is a non-superconducting compound.[19] By introducing Se atoms into AuTe2, electrons from Te have been attracted by the Se anions, and superconductivity with was found in AuTe2Se4/3.[16] Furthermore, it manifests the 2D Berezinsky–Kosterlitz–Thouless (BKT) topological transition based on the measurements of standard IV curves around Tc, which indicates the superconductivity is mainly confined in the ab plane.[16] This is consistent with the fact that the interaction of interlayer is Van der Waals force in AuTe2Se4/3, and the AuTe2Se4/3 single crystal can be very easily exfoliated into thin flakes by Scotch tape.[16] Therefore, AuTe2Se4/3 is a new vdW quasi-2D superconductor.

To understand the superconducting pairing mechanism of a new superconductor, it is crucial to know the superconducting gap structure first. In this work, we measure the ultra-low-temperature thermal conductivity of AuTe2Se4/3 single crystal down to 100 mK, to determine its superconducting gap structure. A negligible residual linear term κ0/T in zero field and a slow field dependence of κ0/T at low field are revealed, which suggest multiple nodeless superconducting gaps in this interesting vdW quasi-2D superconductor. This agrees with our previous band structure calculations in Ref. [16].

2. Experiment

Single crystals of AuTe2Se4/3 were grown by the self-flux method.[16] The obtained single crystals have ribbon shape with shining mirror-like surfaces, as shown in the inset of Fig. 1(a). The largest surface was identified as the ab plane by x-ray diffraction (XRD) measurement, as shown in Fig. 1(a), in which only (00l) Bragg peaks show up. The sample for transport measurements was cut to a rectangular shape of 2.20 mm × 0.60 mm in the ab plane, with a thickness of 0.20 mm perpendicular to the ab plane. Four gold electrodes were fabricated on the sample by thermal evaporation to ensure a good contact at low temperature. Four silver wires were then attached to the electrodes with silver paint, which were used for both resistivity and thermal conductivity measurements, with electrical and heat currents in the ab plane. The contacts are metallic with typical resistance of 5 mΩ at 0.3 K. The thermal conductivity was measured in a dilution refrigerator, using two RuO2 chip thermometers, calibrated in situ against a reference RuO2 thermometer. All the magnetic fields for resistivity and thermal conductivity measurements were applied perpendicular to the ab plane. To ensure a homogeneous field distribution in the sample, all fields were applied at temperature above the temperature Tc.

Fig. 1. (color online) (a) X-ray diffraction pattern of AuTe2Se4/3 single crystal. Only (00l) Bragg peaks show up, demonstrating that the largest surface is ab plane. Inset: Optical image of our AuTe2Se4/3 single crystals. (b) The temperature dependence of resistivity at zero field for AuTe2Se4/3 single crystal. Top inset: The fit of the normal-state resistivity to the Fermi-liquid behavior ρ(T) = ρ0 + AT2 from 3 K to 15 K. Bottom inset: The resistive superconducting transition at low temperature.
3. Results and discussion

Figure 1(b) shows the temperature dependence of the resistivity ρ(T) for AuTe2Se4/3 single crystal at zero field. From the bottom inset of Fig. 1(b), the width of the resistive superconducting transition (10%–90%) is less than 0.05 K. The Tc defined by ρ = 0 is 2.70 K, which is consistent with our previous magnetic susceptibility measurement.[16] From 3 K to 15 K, the normal-state resistivity can be well fitted by the Fermi-liquid behavior ρ(T) = ρ0 + AT2, as shown in the top inset of Fig. 1(b). The obtained residual resistivity ρ0 = 87.3 μΩ · cm gives the residual resistivity ratio (RRR) ρ (290 K)/ρ0 ≈ 3.69.

In order to estimate the upper critical field Hc2(0), the low-temperature resistivity of the AuTe2Se4/3 single crystal was measured in magnetic fields up to 0.7 T, as shown in Fig. 2(a). The temperature dependence of Hc2(T), defined by ρ = 0, is plotted in Fig. 2(b). One can see an apparently linear temperature dependence of Hc2. With a linear fit to the data, Hc2(0) ≈ 0.49 T is roughly estimated. Actually, the linear temperature dependence of Hc2 in Fig. 2(b) is very interesting. It may be explained by a two-band Fermi surface topology as in MgB2,[2022] or an unconventional superconducting state as in heavy-fermion compound UBe13.[23] Density function theory calculation shows there are four bands in AuTe2Se4/3,[16] therefore its linear behavior of Hc2(T) likely originates from multiple-band Fermi surface topology.

Fig. 2. (color online) (a) Low-temperature resistivity of AuTe2Se4/3 single crystal in magnetic fields up to 0.7 T. (b) Temperature dependence of the upper critical field Hc2(T), defined by ρ = 0. The dashed line is a linear fit to the data, which points to Hc2(0) ≈ 0.49 T.

Ultra-low temperature thermal conductivity measurement is an established bulk technique to probe the superconducting gap structure.[24] The temperature dependence of thermal conductivity for the AuTe2Se4/3 single crystal in zero and magnetic fields is plotted in Fig. 3. The measured thermal conductivity contains two contributions, κ = κe + κp, which come from electrons and phonons, respectively. In order to separate the two contributions, all the curves are fitted to κ/T = a + bTα−1 below 0.3 K.[25,26] The two terms aT and bTα represent contributions from electrons and phonons, respectively. The residual linear term κ0/Ta is obtained by extrapolating κ/T to T = 0 K. Due to the specular reflections of phonons at the sample surfaces, the power α in the second term is typically between 2 and 3.[25,26]

Fig. 3. (color online) Low-temperature thermal conductivity of AuTe2Se4/3 single crystal (a) in zero field, and (b) in magnetic fields applied perpendicular to the ab plane. The solid lines represent the fits to κ/T = a + bTα − 1 for the data in different H. The dashed line is the normal-state Wiedemann–Franz law expectation L0/ρ0(0.4 T), with the Lorenz number L0 = 2.45 × 10−8 W · Ω · K−2 and ρ0(0.4 T) = 87.3 μΩ · cm.

In zero field, the fitting gives κ0/T = −9.9 ± 11 μW · K−2 · cm−1 and α = 2.41, as seen in Fig. 3(a). Comparing with our experimental error bar ±5 μW · K−2 · cm−1, the κ0/T of AuTe2Se4/3 in zero field is negligible. For s-wave nodeless superconductors, there are no fermionic quasiparticles to conduct heat in the zero-temperature limit, since all electrons form Cooper pairs.[24,25] Therefore there is no residual linear term of κ0/T, for examples in V3Si and NbSe2.[25,27] However, for nodal superconductors, a substantial κ0/T in zero field contributed by the nodal quasiparticles has been found.[24] For example, κ0/T of the overdoped d-wave cuprate superconductor Tl2Ba2CuO6+δ (Tl-2201, Tc = 15 K) is 1.41 mW · K−2 · cm−1.[28] For the p-wave superconductor Sr2RuO4 (Tc = 1.5 K), κ0/T = 17 mW · K−2 · cm−1 was reported.[29] Therefore, the negligible κ0/T of AuTe2Se4/3 strongly suggests nodeless superconducting gap in it.

In Fig. 3(b), between H = 0 T and 0.5 T, we fit all the curves to obtain the κ0/T for each magnetic field. Since κ/T starts to saturate at 0.40 T, we take H = 0.40 T as the bulk Hc2(0) of AuTe2Se4/3. This value is slightly lower than that determined by resistivity measurements in Fig. 2. Note that the choice of a slightly different Hc2(0) does not affect our discussion on the field dependence of κ0/T below. The obtained κ0/T = 0.28 ± 0.004 mW · K−2 · cm−1 in Hc2(0) = 0.40 T meets the normal-state Wiedemann–Franz law expectation L0/ρ0 = 0.28 mW · K−2 · cm−1 very well, with the Lorenz number L0 = 2.45 × 10−8 W · Ω · K−2 and ρ0(0.4 T) = 87.3 μΩ · cm. The verification of the Wiedemann–Franz law in the normal state demonstrates that our thermal conductivity measurements are reliable.

The field dependence of κ0/T can provide further information on the superconducting gap structure.[24] In Fig. 4, we plot the normalized κ0(H)/T as a function of H/Hc2 for AuTe2Se4/3. For comparison, similar data of the clean s-wave superconductor Nb,[30] the dirty s-wave superconducting alloy InBi,[31] the multiband s-wave superconductor NbSe2,[27] and iron-based superconductor BaFe1.9Ni0.1As2,[32] and an overdoped d-wave cuprate superconductor Tl-2201[28] are also plotted. For a single band s-wave superconductor Nb,[30] the κ(H)/T changes little even up to 40% Hc2, while for the nodal superconductor Tl-2201, a small field can yield a quick growth in the quasiparticle density of states (DOS) due to the Volovik effect, and the low-field κ0(H)/T shows a roughly dependence of .[28] In the case of NbSe2, the distinct κ0(H)/T behavior was well explained by multiple superconducting gaps with different magnitudes.[27]

Fig. 4. (color online) Normalized residual linear term κ0/T of AuTe2Se4/3 as a function of H/Hc2. For comparison, similar data are shown for the clean s-wave superconductor Nb,[30] the dirty s-wave superconducting alloy InBi,[31] the multiband s-wave superconductor NbSe2,[27] and iron-based superconductor BaFe1.9Ni0.1As2,[32] and an overdoped d-wave cuprate superconductor Tl-2201.[28]

From Fig. 4, the field dependence of κ0(H)/T for AuTe2Se4/3 grows slightly faster than that of the dirty s-wave superconductor InBi and lower than that of the multiband s-wave superconductor NbSe2. In fact, it almost overlaps that of BaFe1.9Ni0.1As2. In NbSe2, the gap on the Γ band is approximately one-third of the gap on the other two Fermi surfaces, and the magnetic field first suppresses the superconductivity on the Fermi surface with smaller gap,[27] thus κ0/T(H) increases much rapid at low field. The optimally electron-doped BaFe1.9Ni0.1As2 has a slower field dependence of κ0(H)/T than NbSe2.[32] Previously, it was explained by using the multiple superconducting gaps of its sister compound, the optimally electron-doped BaFe1.85Co0.15As2, which has the average gap values Δ(0) = 6.6 meV and 5.0 meV for hole and electron pockets, respectively.[33] This ratio of different gap magnitudes R = 6.6/5.0 = 1.3 is much smaller than that of NbSe2 (R ≈ 3). Due to the similarity between the κ0(H)/T of AuTe2Se4/3 and BaFe1.9Ni0.1As2, AuTe2Se4/3 is likely also a multiband s-wave superconductor. Indeed, electronic structure calculations of AuTe2Se4/3 show there are four bands crossing the Fermi level along different directions in Brillouin zone, including three 3D bands and one quasi-1D band.[16] In this sense, the ratio of the magnitudes for different superconducting gaps in AuTe2Se4/3 should be around 1.3 as well.

So far, among all multi-band s-wave superconductors, only the iron-based superconductors are unconventional with the possible s±-wave pairing.[34] Therefore, the multiple s-wave gap we observe in AuTe2Se4/3 indicates conventional superconductivity in it. This is reasonable, since unconventional superconductivity usually involves magnetism,[35] while AuTe2Se4/3 is nonmagnetic. As a new vdW quasi-2D superconductor, it will be interesting to fabricate AuTe2Se4/3 thin flakes to investigate whether its superconductivity can be enhanced by gating techniques, like in other quasi-2D compounds.[8,3641]

4. Conclusion

In summary, the superconducting gap structure of the new gold-based quasi-2D superconductor AuTe2Se4/3 has been studied by the ultra-low temperature thermal conductivity measurements. A negligible κ0/T in zero field and the field dependence of κ0/T suggest multiple nodeless superconducting gaps. This is consistent with the band structure calculations which shows four Fermi surfaces in AuTe2Se4/3. To confirm those multiple gaps, further experiments like angle-resolved photoemission spectroscopy (ARPES) and scanning tunneling microscope (STM) are needed.

Reference
[1] Wang Q Y Li Z Zhang W H Zhang Z C Zhang J S Li W Ding H Ou Y B Deng P Chang K Wen J Song C L He K Jia J F Ji S H Wang Y Y Wang L L Chen X Ma X C Xue Q K 2012 Chin. Phys. Lett. 29 037402
[2] Xing Y Zhang H M Fu H L Liu H Sun Y Peng J P Wang F Lin X Ma X C Xue Q K Wang J Xie X C 2015 Science 350 542
[3] Saito Y Nojima T Iwasa Y 2016 Nat. Rev. Mater. 2 16094
[4] Qin S Kim J Niu Q Shih C K 2009 Science 324 1314
[5] Zhang T Cheng P Li W J Sun Y J Wang G Zhu X G He K Wang L L Ma X C Chen X Wang Y Y Liu Y Lin H Q Jia J F Xue Q K 2010 Nat. Phys. 6 104
[6] Reyren N Thiel S Caviglia A D Kourkoutis L F Hammerl G Richter C Schneider C W Kopp T Rütschi A S Jaccard D Gabay M Muller D A Triscone J M Mannhart J 2007 Science 317 1196
[7] Gozar A Logvenov G Kourkoutis L F Bollinger A T Giannuzzi L A Muller D A Bozovic I 2008 Nature 455 782
[8] Ye J T Zhang Y J Akashi R Bahramy M S Arita R Iwasa Y 2012 Science 338 1193
[9] Saito Y Nakamura Y Bahramy M S Kohama Y Ye J Kasahara Y Nakagawa Y Onga M Tokunaga M Nojima T Yanase Y Iwasa Y 2015 Nat. Phys. 12 144
[10] Cao Y Fatemi V Fang S Watanabe K Taniguchi T Kaxiras E Jarillo-Herrero P 2018 Nature 556 43
[11] Cao Y Fatemi V Demir A Fang S Tomarken S L Luo J Y Sanchez-Yamagishi J D Watanabe K Taniguchi T Kaxiras E Ashoori R C Jarillo-Herrero P 2018 Nature 556 80
[12] Jiang D Hu T You L Li Q Li A Wang H Mu G Chen Z Zhang H Yu G Zhu J Sun Q Lin C Xiao H Xie X Jiang M 2014 Nat. Commun. 5 5708
[13] Xi X Zhao L Wang Z Berger H Forro L Shan J Mak K F 2015 Nat. Nanotechnol. 10 765
[14] Staley N E Wu J Eklund P Liu Y Li L Xu Z 2009 Phys. Rev. 80 184505
[15] Xi X Wang Z Zhao W Park J H Law K T Berger H Forro L Shan J Mak K F 2016 Nat. Phys. 12 139
[16] Guo J G Chen X Jia X Y Zhang Q H Liu N Lei H C Li S Y Gu L Jin S F Chen X L 2017 Nat. Commun. 8 871
[17] Duwez P Willens R H Klement W 1960 J. Appl. Phys. 31 1136
[18] Luo H Klement W 1962 J. Chem. Phys. 36 1870
[19] Tsuei C C Newkirk L R 1969 Phys. Rev. 183 619
[20] Gurevich A 2003 Phys. Rev. 67 184515
[21] Gurevich A 2007 Physica C Supercond. 456 160
[22] Buzea C Yamashita T 2001 Supercond. Sci. Technol. 14 R115
[23] Maple M B Chen J W Lambert S E Fisk Z Smith J L Ott H R Brooks J S Naughton M J 1985 Phys. Rev. Lett. 54 477
[24] Shakeripour H Petrovic C Taillefer L 2009 New J. Phys. 11 055065
[25] Sutherland M Hawthorn D G Hill R W Ronning F Wakimoto S Zhang H Proust C Boaknin E Lupien C Taillefer L Liang R Bonn D A Hardy W N Gagnon R Hussey N E Kimura T Nohara M Takagi H 2003 Phys. Rev. 67 174520
[26] Li S Y Bonnemaison J B Payeur A Fournier P Wang C H Chen X H Taillefer L 2008 Phys. Rev. 77 134501
[27] Boaknin E Tanatar M A Paglione J Hawthorn D Ronning F Hill R W Sutherland M Taillefer L Sonier J Hayden S M Brill J W 2003 Phys. Rev. Lett. 90 117003
[28] Proust C Boaknin E Hill R W Taillefer L Mackenzie A P 2002 Phys. Rev. Lett. 89 147003
[29] Suzuki M Tanatar M A Kikugawa N Mao Z Q Maeno Y Ishiguro T 2002 Phys. Rev. Lett. 88 227004
[30] Lowell J Sousa J B 1970 J. Low Temp. Phys. 3 65
[31] Willis J O Ginsberg D M 1976 Phys. Rev. 14 1916
[32] Ding L Dong J K Zhou S Y Guan T Y Qiu X Zhang C Li L J Lin X Cao G H Xu Z A Li S Y 2009 New J. Phys. 11 093018
[33] Terashima K Sekiba Y Bowen J H Nakayama K Kawahara T Sato T Richard P Xu Y M Li L J Cao G H Xu Z A Ding H Takahashi T 2009 Proc. Natl. Acad. Sci. 106 7330
[34] Chen X H Dai P C Feng D L Xiang T Zhang F C 2014 Nat. Sci. Rev. 1 371
[35] Norman M R 2011 Science 332 196
[36] Ueno K Shimotani H Yuan H Ye J Kawasaki M Iwasa Y 2014 J. Phys. Soc. Jpn. 83 032001
[37] Saito Y Kasahara Y Ye J Iwasa Y Nojima T 2015 Science 125 9440
[38] Shi W Ye J Zhang Y Suzuki R Yoshida M Miyazaki J Inoue N Saito Y Iwasa Y 2015 Sci. Rep. 5 12534
[39] Lei B Cui J H Xiang Z J Shang C Wang N Z Ye G J Luo X G Wu T Sun Z Chen X H 2016 Phys. Rev. Lett. 116 077002
[40] Lei B Wang N Z Shang C Meng F B Ma L K Luo X G Wu T Sun Z Wang Y Jiang Z Mao B H Liu Z Yu Y J Zhang Y B Chen X H 2017 Phys. Rev. 95 020503
[41] Ying T P Wang M X Zhao Z Y Zhang Z Z Jia X Y Li Y C Lei B Li Q Yu Y Cheng E J An Z H Zhang Y Yang W Chen X H Li S Y 2018 arXiv: 1802.01484 [cond-mat]