Enhancement of superconductivity under pressure and the magnetic phase diagram of tantalum disulfide single crystals.

In low-dimensional electron systems, charge density waves (CDW) and superconductivity are two of the most fundamental collective quantum phenomena. For all known quasi-two-dimensional superconductors, the origin and exact boundary of the electronic orderings and superconductivity are still attractive problems. Through transport and thermodynamic measurements, we report on the field-temperature phase diagram in 2H-TaS2 single crystals. We show that the superconducting transition temperature (Tc) increases by one order of magnitude from temperatures at 0.98 K up to 9.15 K at 8.7 GPa when the Tc becomes very sharp. Additionally, the effects of 8.7 GPa illustrate a suppression of the CDW ground state, with critically small Fermi surfaces. Below the Tc the lattice of magnetic flux lines melts from a solid-like state to a broad vortex liquid phase region. Our measurements indicate an unconventional s-wave-like picture with two energy gaps evidencing its multi-band nature.

For more than four decades, one of the major subjects in condensed matter physics has been the coexistence of the charge density wave (CDW) order and superconductivity in transition metal dichalcogenides (TMDs)12. In CDW materials such a coupling between the electrons and the soft-phonon mode describes the phase transition from the CDW to a normal state3. The superconducting transition temperature (T) increases while the CDW lock-in temperature falls by doping4, critical thicknesses5, or by external pressure678. Recently, Klemm9 has shown that most of the pristine TMDs are highly unconventional in comparison with conventional superconductors. Amongst many TMD materials, 2H-TaS2 (H: hexagonal, see Methods, Extended Data Fig. 1) becomes superconducting at ambient pressure and without doping4. So far, this compound is one of the very few materials where a chiral and polar charge-ordered phase is suggested to exist1011. Based on scanning tunneling microscopy measurements, the nodal gap structure of a single-layer material has recently been proposed12. Moreover, the lack of agreement on the electronic properties of 2H-TaS2, the information on its magnetic properties and, the Abrikosov vortex dynamics, is also missing up to now13. Therefore, the appearance of superconductivity in 2H-TaS2 in the presence of a CDW is of great interest. This has motivated us to study the low temperature-field dependencies of both transport and thermodynamics in the normal and superconducting states of 2H-TaS2 single crystals to determine their superconducting properties.

Transport Measurements

The temperature dependencies of the in-plane and out-of-plane zero-field resistivity (ρ and ρ) are shown in Fig. 1(a). Both ρ and ρ exhibit a prominent CDW anomaly at 76 K (see The Methods, Extended Data Fig. 2). A parameter often used to characterize the interlayer coupling, is the anisotropy of the resistivity ρ/ρ. The largest anisotropy ratio found here is ρ/ρ ~ 16 just above the T. We noticed in particular that the anisotropy ratio is almost temperature independent. This anisotropy ratio behavior suggests that the in-plane and out-of-plane transport in 2H-TaS2 share the same scattering mechanism. Upon lowering the temperature below the CDW transition, the resistivity displays a drop to zero as shown in the inset of Fig. 1(a). The detailed magnetic field and temperature dependencies of ρ(H) at various temperatures ranging from, 60 mK to 3 K with the field direction parallel to the c-plane of the crystal, are presented in Fig. 2. At low temperatures, the curves are almost parallel to each other in the transition region. With increasing magnetic fields, the onset of superconductivity shifts to lower temperatures gradually. The suppression of superconductivity with a magnetic field applied along the c-direction is more obvious than that in the H||ab configuration, indicating a high anisotropy for a low T in 2H-TaS2.

Figure 1

Transport measurements 2H-TaS2 at ambient and high pressures.

(a) Temperature dependence of in-plane and out-of-plane resistivities at ambient pressure. The lower inset presents a zoom of the in-plane resistivity data around T. The upper inset shows an expanded layered structure of 2H-TaS2. The 2H form is based on edge sharing TaS6 trigonal prisms. Each layer of TaS2 has a strongly bonded 2D S-Ta-S layers, with Ta in either trigonal prismatic or octahedral coordination with S. The chemical bonding within the S-Ta-S layers are covalently bound. (b) Temperature dependence of the in-plane electrical resistivity in zero-field at 3.1 GPa and 8.7 GPa. The inset represents a zoom of the in-plane resistivity data with a very sharp superconducting transition and the T enhances up to 9.15 K at 8.7 GPa.

Figure 2

Low temperature and magnetic field dependencies of 2H-TaS2 resistivity.

(a) The temperature dependence of ρ at different magnetic fields for H||ab and H||c. (b) The inset shows the criterion for determining the T at 0.005 T. (c) The magnetic field dependence of the in-plane resistivity ρ for H||c.

It is worth mentioning that the T in resistivity at ambient pressure is anomalously wide. It is about 0.62 K from the onset T value to the zero resistivity value of the T at 0.98 K, i.e. about 50% of the T. This anomalous could have several sources: chemical or electronic inhomogeneity, fluctuations, or vortex effects. Inhomogeneity is indeed expected to widen the transition of this compound because studies show that even a small concentration of dopants enhance the T dramatically414. This is why superconductivity above 1 K in nominally pure 2H-TaS2 is explained by a small Ta excess or by the presence of sub 1% quantities of impurity atoms4. An intrinsic electronic inhomogeneity related to inhomogeneous CDW is quite possible as the chiral CDW reported for this compound supposes a domain structure4 which is inline with the observed narrowing of the T with CDW destruction, for example with doping by Ni14. On the other hand, both kinds of inhomogeneity could also affect the width of the transition in heat capacity; however that is only about 0.2 K, much narrower than the . This suggests that the effects of these inhomogeneities are limited by 0.2 K, while the rest of the is related to fluctuations and/or vortices. The dissipative vortex motion could either be due to the flux flow through low pinning centers or due to the free motion of individual vortices in the vortex liquid state. Since in our resistivity experiments we used the lowest current 0.1μA and the was not sensitive to its small enhancement, we may consider the vortex liquid state as the most probable mechanism of widening the transition, similarly to CuTiSe215. The vortex liquid can be considered a result of fluctuations in the vortex lattice below the T, while fluctuations in the superconducting order parameter lead to the appearance of preformed pairs above the T. The measurement of both fluctuation regions is the Ginzburg number Gi = δT/T, which is usually extremely small for the low-temperature superconductors Gi ~ (T/E)4 ~ 10−12–10−14, even for two-dimensional ones, for which Gi ~ T/E or τ−1/E for the clean and dirty limits respectively16. Here δT is the range of temperatures in which fluctuation corrections are relevant, and τ−1 is the quasiparticle scattering rate at the Fermi energy (E). However, in the CDW state, the reconstructed FS may have small and very shallow pockets for which the E could not be much larger than T or 1/τ17. Therefore, the broadening of the T due to the interplay with CDW is further supported by the sharp T after the suppression of the CDW upon compression.

Enhancement of Superconductivity Upon Compression

In low-dimensional electron systems, CDW and superconductivity are two of the most fundamental collective quantum phenomena12. Unconventional superconductivity is nearly always found in the vicinity of another ordered state, such as antiferromagnetism, CDW, or stripe order. This suggests a fundamental connection between superconductivity and fluctuations in some other order parameter18. To better understand this connection, we used high-pressure resistivity to directly study the CDW order in 2H-TaS2. The effect of pressure on 2H-TaS2 is presented in Fig. 1(b). Upon 3.1 GPa, the CDW slightly shifts to 69 K. The effects of 8.7 GPa illustrate a suppression of the T. In addition, a very sharp drop in resistivity indicates the onset of superconductivity and dramatically enhances the modest T to ~9.15 K upon 8.7 GPa. Similarly to recently reported data19, our resistance measurements show that the T increases from temperatures below 1 K up to 8.5 K at 9.5 GPa. Additionally, the authors observed a kink in the pressure dependence of T at about 4 GPa that they attributed to the lock-in transition from an incommensurate CDW to a commensurate CDW. Above this pressure, the commensurate T slowly decreases, coexisting with superconductivity within our full pressure range. These observations show that the enhancement in superconductivity is due to the consequent changes of Fermi surface (FS) upon compression. However, this is not direct evidence that confirms where such features act on superconductivity independently of the CDW. In the CDW state, a gap opens up over part of the FS in the direction of the q vectors of the CDW8. This reduces the average density of states at the FS. Upon compression, T is suppressed. The amplitude of the CDW lattice distortion also suppresses, thus gradually restoring the FS and increasing the T. Therefore, one can see that both superconductivity and the CDW involve widely different parts of the FS associated with the absence of or small interband correlations. It is worth noting that superconductivity in 2H-NbSe2 is only moderately affected by pressure2021 and the CDW already disappears at 5 GPa2022. The weak pressure dependence of the T at higher pressures indicates that the CDW in this pressure range is remarkably robust to a reduction in the lattice parameters19. Very recently23, in 2H-NbSe2 the rapid destruction of the CDW under pressure was found to be related to the quantum fluctuations of the lattice renormalized by the anharmonic part of the lattice potential. In addition, the connection between CDWs and superconductivity arises from the fact that high-energy optical phonon modes have a strong contribution to the Eliashberg function, whereas the low-energy longitudinal acoustic mode that drives the CDW transition barely contributes to superconductivity

Specific Heat Measurements

To further elucidate the bulk superconductivity in 2H-TaSe2, we performed heat capacity studies down to 70 mK. Figure 3 summarizes the T-dependence of the specific heat data in various magnetic fields applied parallel and perpendicular to the ab plane. We observed a clear sharp anomaly at T = 1.4 K, close to that determined by our resistivity measurements. The specific heat jump systematically shifted to lower temperatures upon the application of magnetic fields. Our data of small fields close to the T shows the evolution of a small fluctuation, peak, overlapped with the specific-heat jump [see inset of Fig. 3(b)]. On the other hand, both kinds of chemical or electronic inhomogeneity should also affect the width of the transition in heat capacity, however that is only about 0.2 K, much narrower than the . This suggests that the effects of inhomogeneities are limited by 0.2 K, while the rest of could be related to fluctuations and/or vortices. A clear maximum of specific heat data at 76 K, typically found in 2H-TaS2 which is weakly first-order, is an indication of the CDW transition [see the inset of Fig. 3(a)]. Note that there is no upturn (Schottky nuclear contribution) in the specific heat data measured to temperatures as low as 70 mK, thus, the zero-field specific heat above T can be well fitted to C/T = γ + βT2, where γ and β are the electronic and lattice coefficients, respectively [see the dashed line in Fig. 3(b)]. The γ value is found to be around 8.8 mJ/mol K2, indicating that 2H-TaS2 in the CDW state is characterized by a modest density of states. This value agrees with the γ value found by refs 4 and 24 in which C was just measured between 1.8 and 10 K. The phononic coefficient β is found to be 0.35 mJ/mol K4. Using the relation θ = (12π4RN/5β)1/3, we obtained the Debye temperature θ = 249(2) K, which is comparable with values reported by DiSalvo et al.2. From the determined γn value, we found that ΔCel/γnTc = 0.72. This value is smaller than the prediction of the weak coupling BCS theory (ΔCel/γnTc = 1.43) and comparable to that in the intercalated compound24. This indicates that the specific-heat data cannot be described by a simple BCS gap (see Methods, Extended Data Fig. 3). However, in a clean situation with negligible pair-breaking effects, the reduced jump in the specific heat compared to that of a single-band s-wave superconductor might be related to unconventional superconductivity with nodes and/or a pronounced multiband character with rather different partial densities of states and gap values25. In addition, evidence of coupling effects arises from the normalized discontinuity value of the specific-heat slopes at the T, (T/ΔC)(dC/dT. In the single-band weak coupling BCS theory this ratio is 2.64, whereas a value of 3.35 can be deduced in the two-band superconductor MgB226. From our data, we obtained a value of (T/ΔC)(dC/dT ~ 3.54, which is very close to MgB2.

Figure 3

Temperature dependence of 2H-TaS2 specific heat.

T-dependence of the specific heat in various applied magnetic fields parallel to the ab axis (a) and parallel to the c plane (b). The dashed line in (b) is the fitting below 2.5 K by using Cp = γT + βT3. The inset in (b) shows a close-up of the superconducting state while the inset of (a) presents the CDW state.

H - T Phase Diagram

The H provides a valuable insight into the nature of the interaction responsible for the formation of Cooper pairs252728. The temperature dependencies of H and H obtained from C(T, H) and ρ(T, H) with both H||c and H||ab are plotted in Fig. 4 for both orientations. Specific heat T(H) values were deduced from the classical entropy conservation construction. The criteria of the normal state in resistivity was used to extract the T at each magnetic field. The irreversible magnetic field H was obtained from the zero value of T in ρ curves. However, the width of the resistive transition is shown in the inset of Fig. 4(b) and is proportional to μ0H2/3. This is inline with Tinkham’s theoretical prediction29 of the ΔT ∝ μ0H2/3. The large area between the H and H curves suggests that the vortex dissipation level is still low in this region. Moreover, the possible existence of a distinct H(T) far below H is due to the fact that the vortex lattices are soft and easily melted into vortex liquid by the magnetic field or thermal fluctuations30. The zero-temperature values for and are estimated to be approximately 0.31 and 1.38 T, respectively. From those we estimated the anisotropic coherence length  = 32.6 nm, and ξ = 7.3 nm. One can also estimate the coherence length from the uncertainty principle and BCS model. From the Faber-Pippard ratio, ξ = 0.18ħv/kT = 260 nm, for T = 1.4 K and an average Fermi velocity v ≈ 1.5 eVÅ17, which is considered to be similar for 2H-TaSe2, 2H-NbSe2, and 2H-TaS2. This shows that both anisotropy and CDW effects on electronic structure should be taken into account. Furthermore, it has been reported31 that the field-induced antiferromagnetism can extend outside the effective vortex core region where the superconducting order parameter is finite. Such an extended magnetic order is expected to suppress the superconducting order parameter around vortices. This effect will enlarge the vortex core size, which in turn will suppress the H. The effective core size has been found to be around three times that of the coherence length in 2H-NbSe232. From the behavior of Hc2 vs. T for the different field orientations, we have calculated the anisotropy as  = ξ/ξ. The anisotropy Γ increases upon approaching the Tc and reaches about 4(1). This indicates that the orbital pair breaking also accounts for the suppression of superconductivity close to T in 2H-TaS2.

Figure 4

H - T phase diagram of 2H-TaS2.

(a) The upper critical field μ0H and the irreversible field, μ0H for H||ab (a) and H||c (b). Open symbols in (b) are taken from ρ(H). The inset illustrates the transition width (ΔT vs. μ0H2/3). The dashed line is the linear fit.

In the case of multi-band superconductivity33343536 the low-temperature H-curve may exceed the single-band Werthamer-Helfand-Hohenberg predictions37. However, a noticeable upward curvature in the H(T) observed in some compounds has been attributed to multiband effects38. Using typical renormalized Fermi velocities derived from preliminary ARPES-data17 and T = 1.4 K, one also estimates, that in principle by a two-band approach adopting s-symmetry38, the slope-value is: Hc2, = , where c1 → c2 → 1/2 and in the case of a dominant interband pairing results in - = 0.14 T/K near the T which is already very close to our experimentally determined value. By fitting it using the two-band theory333940, one can obtain the band diffusivities D1, D2 and the intraband and interband coupling constants λ11, λ12, and λ21. The exact relations can be found in ref. 38. Using the band diffusivity ratio η = D2/D1 = 800, λ11 = 0.5, and λ12 = λ21 = 0.25, we fitted our data for 2H-TaS2. The obtained two-band fitting agrees well with the experimental data. To add more insight to the pairing symmetry for the 2H-TaS2 superconductor, we investigated the temperature dependence of the specific heat. The detailed electronic specific heat data and analysis are given in the Extended Data Fig. 3.

Summarizing, we have reported the first superconducting fluctuations investigation across the effect of pressure on the CDW state in 2H-TaS2. From an extensive thermodynamic study, we found a considerable broadening of the T at ambient pressure and its sharp transition at high pressures together with an unexpectedly broad region of vortex liquid phase in the vortex phase diagram. These results suggest the presence of the the superconducting fluctuations in the CDW state. Besides of a clear fundamental interest in our system, this finding can be used to control the fluctuations in quantum devices.

Methods Summary

Low-temperature transport (down to 60 mK) and specific heat (down to 70 mK) measurements were performed using a dilution refrigerator. The conductance anisotropy in layered material single crystals is large therefore using traditional four-terminal methods to determine the resistivity along the c axis, ρ, and in the ab plane, ρ, may be unreliable41. We used six terminals to determine each principal component of resistivity. In the latter method, the current was injected through the outermost contacts on one surface. Voltages were measured across the innermost contacts of each surface. The Laplace equation was then solved and inverted to find ρ and ρ42. In addition, this method allowed testing the sample homogeneity by permuting the electrodes which were used for the current and voltage4142. Four contacts were used to measure the high-pressure in-plane resistivity. The investigated 2H-TaS2 single crystals were synthesized at hq graphene and were of high purity (>99.995%). The resistivity and specific heat measurements down to 0.4 K were measured in a Physical Property Measurement System (Quantum Design) with an adiabatic thermal relaxation technique.

Online Content

Any additional Methods, Extended Data display items and Source Data are available in the online version of the paper; references unique to these sections appear only in the online paper.

Additional Information

How to cite this article: Abdel-Hafiez, M. et al. Enhancement of superconductivity under pressure and the magnetic phase diagram of tantalum disulfide single crystals. Sci. Rep. 6, 31824; doi: 10.1038/srep31824 (2016).