Abstract
The heavy fermion paramagnet UTe2 exhibits numerous characteristics of spin-triplet superconductivity. Efforts to understand the microscopic details of this exotic superconductivity have been impeded by uncertainty regarding the underlying electronic structure. Here we directly probe the Fermi surface of UTe2 by measuring magnetic quantum oscillations in pristine quality crystals. We find an angular profile of quantum oscillatory frequency and amplitude that is characteristic of a quasi-2D Fermi surface, which we find is well described by two cylindrical Fermi sheets of electron- and hole-type respectively. Additionally, we find that both cylindrical Fermi sheets possess considerable undulation but negligible small-scale corrugation, which may allow for their near-nesting and therefore promote magnetic fluctuations that enhance the triplet pairing mechanism. Importantly, we find no evidence for the presence of any 3D Fermi surface sections. Our results place strong constraints on the possible symmetry of the superconducting order parameter in UTe2.
Similar content being viewed by others
Introduction
Conventional phonon-mediated superconductivity involves the pairing of two fermions in a spin-singlet configuration1 forming a bosonic quasiparticle of total spin Sâ=â0. Unconventional superconductors may replace the attractive role played by phonons with a magnetically-mediated pairing interaction, often yielding a d-wave symmetry of the orbital wavefunction, but still with overall Sâ=â0 for the bound pair2. By contrast, the formation of superfluidity in 3He involves a triplet pairing configuration with Sâ=â1 and an odd-parity p-wave symmetry3. To date no bulk solid-state analogue of this exotic state of matter has been unequivocally identified, although actinide metals such as UPt3 and UGe2 are among the promising candidates4,5. The technological realisation of devices incorporating p-wave superconductivity is highly desirable, due to their expected ability to effect coherent quantum information processing6. For several years the layered perovskite Sr2RuO4 appeared a likely host of spin-triplet superconductivity7; however, recent experimental observations have cast considerable doubt on this interpretation8.
Since the discovery of unconventional superconductivity in UTe2 in 20199, characteristic features of a p-wave superconducting state in this material have been reported across numerous physical properties. These include a negligible change in the Knight shift upon cooling through Tc as probed by nuclear magnetic resonance (NMR)10, high upper critical fields far in excess of the Pauli paramagnetic limit11, high magnetic field re-entrant superconductivity12, chiral in-gap states measured by scanning tunneling microscopy (STM)13, time-reversal symmetry breaking inferred from the development of a finite polar Kerr rotation angle below Tc14,15, multiple point nodes detected from penetration depth measurements indicative of a chiral triplet pairing symmetry16, anomalous normal fluid properties consistent with Majorana surface arcs17, and ferromagnetic fluctuations coexisting with superconductivity measured by muon spin relaxation measurements18.
Several theoretical studies have sought to provide a microscopic description of how the exotic magnetic and superconducting features manifest in this material11. However, an outstanding challenge concerns the determination of the underlying electronic structure, with the question of the geometry and topology of the materialâs Fermi surface having been the subject of recent debate and speculation11,19,20,21,22,23,24,25,26,27,28,29,30. Two angle-resolved photoemission spectroscopy (ARPES) studies have given contrasting interpretations, with one inferring the presence of multiple small 3D Fermi surface pockets31, while the other identified spectral features characteristic of a large cylindrical quasi-2D Fermi surface section, along with possible heavy 3D section(s)32. A recent de Haas-van Alphen (dHvA) effect study33 also resolved features representative of a quasi-2D Fermi surface, but with a spectrally dominant low-frequency branch not captured by density functional theory (DFT) calculations, which could be indicative of a 3D Fermi surface pocket. Discerning the dimensionality of the UTe2 Fermi surface is important in order to determine the symmetry of the superconducting order parameter34.
Results
Here, we report direct measurements of the UTe2 Fermi surface, probed by magnetic torque measurements of the dHvA effect up to magnetic field strengths of 28 T at various temperatures down to 19 mK. Through measurements of magnetic quantum oscillations as a function of temperature and magnetic field tilt angle, we observe Fermi surface sections with heavy cyclotron effective masses up to 78(2) me, where me is the bare electron mass. We investigated the angular dependence of the quantum oscillations and used the resulting data to perform Fermi surface simulations. Our results indicate that the UTe2 Fermi surface is very well described by two cylindrical Fermi surface sheets of equal volumes and super-elliptical cross sections. In combination, the evolution of quantum oscillatory amplitude and frequency with the tilt angle of the magnetic field, our quantitative analysis of the oscillatory waveform and the contributions from separate Fermi sheets, and the correspondence between the density of states implied from specific heat measurements and our dHvA observations, makes the presence of any 3D Fermi surface pocket(s) very unlikely. We also measured Shubnikov-de Haas (SdH) oscillations in the contactless resistivity (see Supplementary Information). The SdH effect is generally more sensitive than the magnetic torque technique to symmetrical 3D Fermi pockets35. However, we find that the SdH response of UTe2 contains no additional frequency components than those probed by dHvA (up to field strengths of 28 T, see Supplementary Information). Furthermore, we find that our quasi-2D Fermi surface model very well describes the relatively small anisotropy of this materialâs electrical conductivity tensor. In combination these measurements and analysis provide no evidence indicating the presence of 3D sections in the Fermi surface of UTe2, which we instead find to be quasi-2D in nature and composed of two undulating cylindrical sections of hole- and electron-type, respectively.
Samples were grown by the molten salt flux (MSF) technique36 in excess uranium, to minimise the formation of uranium vacancies (see Methods for details). The MSF technique has been found to produce crystals of exceptionally high quality36, as demonstrated by specific heat capacity, Cp, and electrical resistivity, Ï, measurements in Fig. 1. For this batch of crystals on which quantum oscillation studies were performed, we observe a superconducting transition temperature (Tc) of 2.1 K and residual resistivity ratios (RRR) of up to 900. The RRR is defined as Ï (300 K)/Ï0, where Ï0 is the residual 0 K resistivity expected for the normal state in the absence of superconductivity, fitted by the dashed line (linear in T2) in Fig. 1b. By comparison, early UTe2 sample generations grown by the chemical vapour transport method exhibited RRR values ofâââ40 with Tc values ofâââ1.6 K (refs. 9,11,37). Optimization of the CVT growth parameters has yielded higher quality crystals with Tc = 2.0 K at a maximal reported RRR of 88 (ref. 38) â however, even these optimized samples exhibit uranium site vacancies ofâââ0.7%, compared to neglibile vacancies for MSF samples36,39. Furthermore, in this study we resolve quantum oscillatory components with frequencies up to 18.5 kT, implying a mean free path of itinerant quasiparticles of at least 1900 Ã (see Supplementary Information for calculation), further underlining the pristine quality of this new generation of UTe2 samples. With the advent of these higher quality MSF crystals, it may be necessary to revisit many prior experiments reported on CVT samples â for example, while this manuscript has been undergoing peer review, measurements of the NMR Knight shift40, Kerr effect41, thermal transport42 and muon spectroscopy43 have been reported on MSF specimens, each of which contrasts with prior interpretations drawn from CVT studies.
Figure 2 shows quantum oscillations measured in the magnetic torque of UTe2. The oscillatory component of the signal, ÎÏ, was isolated from the background magnetic torque by subtracting a smooth monotonic local polynomial regression fit (see Methods). In panel (a) we find that when a magnetic field, \(\overrightarrow{H}\), is applied along the \(\overrightarrow{c}\) direction (defined here as 0°), a monofrequency oscillatory waveform of large amplitude is observed. Upon rotating 8° away from \(\overrightarrow{c}\) towards the \(\overrightarrow{a}\) direction, we find here that ÎÏ incorporates a subtle beat pattern, with the corresponding Fourier spectra indicating the presence of at least one other frequency component in addition to the (now diminished) dominant frequency peak. As the field is tilted further away from the \(\overrightarrow{c}\) direction we find that this trend progresses: by 20° the oscillatory amplitude has diminished considerably, accompanied by further splitting and broadening of the FFT spectra, while at 74°ÎÏ appears almost featureless when plotted on this scale.
Figure 2 c compares the raw torque signal, Ï, at 0° and 74°. Both curves have been translated (without rescaling) to the same value at 26 T, for ease of comparison. A clear oscillatory component is visible in the 0Â°Ï curve, whereas the 74° data appear very smooth. Figure 2d, e give a zoomed-in view of ÎÏ at 74°, in which a very small amplitude (over an order of magnitude smaller than at 0°), high-frequency component of 18.5 kT is clearly present. In the Supplementary Information we perform a quantitative analysis of the 0° waveform, which reveals the oscillatory contribution of two distinct Fermi surface sections of identical cross-sectional area.
This angular evolution of the dHvA effect â of a singular, relatively low-frequency oscillatory component along a high symmetry direction subsequently evolving under rotation to yield much higher frequencies of smaller amplitude â is characteristic of a cylindrical quasi-2D Fermi surface. This is due to there being negligible smearing of the quantum oscillatory phase when averaging over kz along the cylindrical axis. Hence, one observes a large quantum oscillatory amplitude for magnetic field oriented in this direction (in this case the \(\overrightarrow{c}\) direction). Then, as the field is tilted away from the axis of the cylinder, the oscillatory amplitude falls considerably due to phase smearing that increases with the rate of change of the frequency with angle. Furthermore, the oscillatory frequency increases at higher angles because the cross-sectional area of the Fermi surface normal to the magnetic field grows as \(\frac{1}{\cos \theta }\) (ref. 35). Thus, our observation of the progression from large amplitude, low-frequency quantum oscillations with field oriented along \(\overrightarrow{c}\) evolving to small amplitude, high-frequency oscillations for field applied close to \(\overrightarrow{a}\) is strongly indicative of UTe2 possessing cylindrical quasi-2D Fermi surface sections, axially collinear with the \(\overrightarrow{c}\) direction.
Figure 3 shows the evolution in temperature of ÎÏ for magnetic field oriented along the \(\overrightarrow{c}\) direction. The quantum oscillation amplitude diminishes rapidly at elevated temperatures, with the signal at 200 mK being an order of magnitude smaller than at 19 mK. Figure 3c fits the quantum oscillatory amplitude to the temperature dependence of the Lifshitz-Kosevich formula35 (see Methods for details), yielding a heavy effective cyclotron mass, m*, of 41(2) me, consistent with observations reported in ref. 33. For an inclined angle of the magnetic field we observe heavier effective masses up to 78(2) me (see Supplementary Information).
We plot the angular evolution of quantum oscillatory frequency with magnetic field tilt angle in Fig. 4 for both the \(\overrightarrow{c}\) to \(\overrightarrow{a}\) and \(\overrightarrow{c}\) to \(\overrightarrow{b}\) rotation planes. Due to experimental limitations, our \(\overrightarrow{c}\) to \(\overrightarrow{b}\) measurements were constrained to within 45° of rotation. We also plot the dHvA frequency simulation for our calculated Fermi surface sections (see Supplementary Information for simulation details), and find remarkably good correspondence between measurement and theory for all frequency branches. Thus, we find that the dHvA profile of UTe2 is excellently described by two quasi-2D Fermi sheets, of âsquircularâ (super-elliptical) cross-section, with the hole- (electron-) type sheet centred at the X (Y) point of the Brillouin zone (Fig. 4c).
We note that in performing DFT calculations (see Supplementary Information) we were unable to capture the angular profile of the low, spectrally dominant frequency branch in the \(\overrightarrow{c}\) to \(\overrightarrow{a}\) rotation plane. This branch initially decreases in frequency as field is titled away from \(\overrightarrow{c}\), reaching a minimum at around 25°, before then increasing rapidly close to \(\overrightarrow{a}\). This is in sharp contrast to the expectation for a cylindrical Fermi surface section of circular cross-section with no undulations, for which the corresponding frequencies should increase monotonically under rotation away from the cylindrical axis. This property of the quantum oscillatory spectra thus places strong constraints on the possible Fermi surface geometry, and motivated our computation of data-led Fermi surface simulations (see Supplementary Information for simulation details). These simulations yielded cylindrical Fermi sheets with super-elliptical rather than circular cross-sections, with significant undulation along their lengths, but with a singular cross-sectional area present at extremal points normal to \(\overrightarrow{c}\). We find that this simulation excellently describes the evolution of the three frequency branches observed at intermediate angles (Fig. 4a).
We can compare the density of states at the Fermi energy inferred from these quantum oscillation experiments with that determined by measurements of the linear specific heat coefficient. Our specific heat measurements (in ambient magnetic field) give a residual normal state Sommerfeld coefficient γNâ=â121(1) mJ molâ1 Kâ2, consistent with prior reports11. Assuming that the quasiparticles of both the hole- and electron-sheets have m* = 41(2) me (Fig. 3c), we numerically calculated the contribution to the density of states from both Fermi surface sections (see Methods for calculation). We find that our simulated electron-type sheet should contribute â61.8 mJ molâ1 Kâ2, while the hole-type section should giveâââ58.7 mJ molâ1 Kâ2. Together the two cylinders thus comprise a total density of states at the Fermi energy corresponding to a Sommerfeld ratio of â120.5 mJ molâ1 Kâ2, in excellent agreement with the heat capacity measurements of γN. We therefore conclude that these two sheets are very likely the only Fermi surface sections present in UTe2.
A number of studies have sought to reconcile anomalous aspects of the superconducting and normal state properties of UTe2 by evoking models that assume the presence of one or more 3D Fermi surface sections22,23,24,25,26,27,28,29,30. The distinction between having a 3D or quasi-2D Fermi surface is important, as this sets strong constraints on the possible irreducible representations of the point group symmetry for the superconducting order parameter34. In the absence of direct observations clarifying the Fermi surface dimensionality these models appeared promising due to their apparent ability to explain physical properties, such as the slight anisotropy of the electrical conductivity tensor27. However, these studies did not consider the potential effects of pronounced undulations along the \(\overrightarrow{c}\) direction of quasi-2D cylindrical Fermi surface sheets. In the Supplementary Information we model the expected components of the electrical conductivity tensor for our simulated Fermi surface44,45. Due to the pronounced undulations giving a sizeable component in the \(\overrightarrow{c}\) direction to the unit vectors normal to the Fermi surface, we find that the low anisotropy of the electrical conductivity of UTe2 can be naturally explained by this pronounced warping of the quasi-2D cylindrical Fermi surface sheets (Fig. 5), without requiring any 3D Fermi surface pocket(s).
Furthermore, our angle-dependent dHvA measurements and corresponding Fermi surface simulation clearly resolve that the Fermi surface of UTe2 comprises two cylindrical sections, possessing quasiparticle effective masses that fully account for the linear specific heat coefficient. Therefore, any potential 3D Fermi surface sections must be very small in size such that their quantum oscillatory frequencies would be very low, or have very high effective masses to necessitate measurement temperatures considerably lower than 19 mK, or possess very high curvature around their entire surface so as to minimise the phase coherence of intersections with successive Landau tubes. However, as the contribution to the density of states per Fermi surface section is directly proportional to the effective mass of the quasiparticles hosted by that section (see Methods), it therefore seems unlikely that UTe2 possesses any 3D Fermi pockets with markedly heavier effective masses than the cylindrical sections we observe in our temperature-dependent measurements.
We add further confidence to our interpretation of the dHvA data by a quasi-2D Fermi surface model by fitting the 0° 19 mK ÎÏ curve from Fig. 3 as the sum of two distinct oscillatory components representing the hole- and electron-type sections, respectively (see Supplementary Information). The fit yields two sinusoids of identical frequencies (within error). Notably, there is only a very slight phase-smearing contribution, indicative of negligible small-scale corrugations along the lengths of the cylinders. Performing a Fourier analysis of the residual curve obtained by subtracting the fit from the data (Supplementary Fig. S15) reveals no additional frequency components that may have been obscured by the large peak in the 0° FFT in Fig. 2b. Therefore, we conclude that our dHvA and SdH data are very well interpreted based on UTe2 possessing a heavy, quasi-2D, charge-compensated Fermi surface.
Discussion
Our finding of a quasi-2D Fermi surface in UTe2 has important implications for determining the symmetry of the superconducting gap structure. Evidence indicating the presence of point nodes has been reported from thermal conductivity measurements46, with subsequent NMR47 and scanning SQUID48 studies interpreting the gap structure within the D2h point group as being of B3u character, with point nodes along the \(\overrightarrow{a}\) direction. However, recent NMR40 and thermal conductivity42 measurements performed on high-quality MSF crystals argue strongly in favour of a highly anisotropic full gap, of the Au representation. The quasi-2D cylinders of our Fermi surface model would be consistent with either a nodal B3u or fully gapped Au scenario. Further experimental and theoretical studies to distinguish between these symmetries â or the possibility of a nonunitary combination of both16,47,49,50 â are urgently called for to provide a complete microscopic understanding of the superconducting order parameter in UTe2.
It is interesting to consider how different the Fermi surface of UTe2 is from the complex multisheet Fermi surfaces found in the ferromagnetic superconductors URhGe, UGe2, and UCoGe, the latter two of which have small 3D pockets51. Contrastingly, the possession of a relatively simple quasi-2D Fermi surface comprising charge-compensated cylindrical components is remarkably similar to other unconventional superconductors, including the Fe-pnictides52,53, underdoped high-Tc cuprates54,55, and Pu-based superconductors56. It has been suggested19 that the near nesting between quasi-2D Fermi surface sections favours spin fluctuations in UTe2 and may thereby strengthen the spin-triplet pairing mechanism. Thus, given the pronounced undulation but negligible degree of small-scale corrugation of the UTe2 Fermi sheets, this is likely the reason why UTe2 exhibits such a markedly higher Tc than its ferromagnetic U-based cousins. Given the multitude of theoretical study into the effects of d-wave pairing symmetry hosted by such a Fermi surface in the case of e.g. the cuprates57, it is therefore interesting to consider what similarities and differences may be found when considering instead a p-wave symmetry.
At the first-order metamagnetic transition obtained at μ0Hâââ35 T for \(\overrightarrow{H}\parallel \overrightarrow{b}\), a Fermi surface reconstruction has been proposed to occur11 due to the observation of a change of sign of the Seebeck coefficient and a reported discontinuity in the carrier density as interpreted from Hall effect measurements58. This raises the interesting possibility that the high field re-entrant superconducting phase12, which is acutely angle-dependent and appears to persist to at least 70 T (ref. 59), may be markedly different in character compared to the superconductivity found below 35 T. We note that our Fermi surface simulations predict the occurrence of a Yamaji angle60 at the magnetic field orientation where the high field re-entrant superconducting phase is most pronounced (see Supplementary Fig. S22), similar to prior reports59,61. This implies that a sharp peak in the field-dependent density of states may underpin the microscopic mechanism driving this exotic superconducting state. A further quantum oscillation study beyond the scope of this work, in the experimentally challenging temperature-field regime of Tâ<â100 mK and μ0Hâ>â35 T, would therefore be of great interest in comparing the underlying fermiology of the magnetic field re-entrant superconducting phase with that of the lower field Fermi surface we uncover here.
We note that while our manuscript was undergoing peer-review, we became aware of a study of the oscillatory magnetoconductance of UTe2 in which the authors interpreted their measurements as being indicative of the presence of small, light, 3D Fermi surface sections62. This was surprising, as we observed no signatures of such pockets in our dHvA or SdH measurements, nor did an independent group in their detailed and careful dHvA measurements33,63. However, we have since been able to reproduce the magnetoconductance oscillations reported by ref. 62. We found that such observations could not actually be reconciled with the presence of 3D Fermi sections, but instead can be well understood as a result of high magnetic field-induced quasiparticle tunnelling between the cylindrical sheets, analogous to prior studies of quasi-2D organic metals in which such observations were also present in kinetic properties such as the conductance but, notably, were not observed in derivatives of the free energy such as the magnetization64,65. These findings will be reported in detail elsewhere66.
In conclusion, our quantum oscillation study on pristine quality crystals has revealed the quasi-2D nature of the Fermi surface in UTe2. Performing dHvA and SdH measurements in magnetic field strengths well above the \(\overrightarrow{c}\)-axis upper critical field have enabled us to compute the Fermi surface geometry, which consists of two cylindrical sheets of super-elliptical cross-section with a significant degree of undulation along their lengths. We find that this pronounced undulation can account for the relatively modest anisotropy of the electrical conductivity tensor, which had previously been cited as evidence favouring the presence of 3D Fermi surface pockets â however, we find no evidence for the presence of any 3D sections. We numerically calculated the contribution to the density of states for our computed Fermi surface, and find that it fully accounts for the normal state Sommerfeld ratio determined from specific heat measurements. Our findings indicate that despite having pronounced axial undulations the Fermi surface of UTe2 possesses a negligible degree of small-scale corrugation, implying that the Fermi sheets may nest very closely together, thereby favouring magnetic fluctuations that enhance the spin-triplet superconducting pairing mechanism.
Methods
Sample preparation
High-quality single crystals of UTe2 were grown using the molten salt flux technique adapted from ref. 36. We used an equimolar mixture of powdered NaCl (99.99%) and KCl (99.999%) salts as a flux, which had been dried at 200 °C for 24 hours. Natural uranium metal with an initial purity of 99.9% was further refined using the solid-state electrotransport (SSE) method67 under ultra-high vacuum (â~â10â10 mbar); by passing a high electrical current of 400 A through the initial uranium metal, impurities can be removed extremely effectively.
Following SSE treatment, a piece of purified uranium of typical mass â0.35 g was etched using nitric acid to remove surface oxides. It was subsequently placed in a carbon crucible of inner diameter 13 mm together with pieces of tellurium (99.9999%) with the molar ratio of 1:1.71; subsequently, the equimolar mixture of NaCl and KCl was added. The molar ratio of uranium to NaCl,KCl mixture was 1:60. The whole process was performed under a protective argon atmosphere in a glovebox. The carbon crucible was plugged by quartz wool, placed in a quartz tube, and heated up to 200â°C under dynamic high vacuum (~10â6 mbar) for 12 hours. Then it was sealed and placed in a furnace. It was initially heated to 450â°C in 24 hours, and left there for a further 24 hours. Then it was heated to 950°C at a rate of 0.35°C/min and kept there for an additional 24 hours. Afterwards the temperature was slowly decreased at a rate of 0.03°C/min down to 650°C, maintained there for 24 hours, and then cooled down to room temperature during the following 24 hours.
After the growth process, the ampoules were crushed and the contents of the carbon crucibles were immersed in water where the salts rapidly dissolved. Bar-shaped crystals were manually removed from the solution, rinsed with acetone, and stored under an argon atmosphere prior to characterisation and quantum oscillation studies. The longest edge of the produced single crystals was typically 3-12 mm (along the \(\overrightarrow{a}\) direction), with widths 0.5-1.2 mm (along \(\overrightarrow{b}\)) and thicknesses around 0.2-1 mm (along \(\overrightarrow{c}\)).
Capacitive torque magnetometry measurements
Torque magnetometry measurements were performed at the National High Magnetic Field Laboratory, Tallahassee, Florida, USA. Measurements were taken in SCM4 fitted with a dilution refrigerator sample environment. Single crystal samples were oriented with a Laue X-ray diffractometer. We note that our angular data obtained in the \(\overrightarrow{c}\)â\(\overrightarrow{a}\) plane is calibrated to withinâââ2° of experimental uncertainty; however, a possible azimuthal offset in the \(\overrightarrow{c}\)â\(\overrightarrow{b}\) angles means that these data should only be taken to be accurate to within â5°. Samples were mounted on beryllium copper cantilevers suspended above a copper base plate, thereby forming a capacitive circuit component. The capacitance of the cantileverâbase plate system was measured as a function of applied magnetic field strength by a General Radio analogue capacitance bridge in conjunction with a phase-sensitive detector. This configuration of cantilever and base plate was mounted on a custom-built rotatable housing unit, allowing for the angular dependence of the dHvA effect to be studied.
In analysing the measured torque data, the oscillatory component was isolated from the background magnetic torque by subtracting a smooth monotonic polynomial fit by the local regression technique68. The main benefit of this technique over simply subtracting a polynomial fitted over the whole field range is that the LOESS window over which the averaging occurs can be modified; for oscillations of faster (slower) frequency, smaller (larger) LOESS windows will achieve a better isolation of the dHvA effect signal. This averaging window then slides along the entire curve to extract the oscillatory component from the background magnetic torque.
Lifshitz-Kosevich temperature study
Figure 3 a shows quantum oscillations measured at various temperatures, with the quantum oscillatory amplitude being strongly diminished at elevated temperature. We extract an effective cyclotron mass, m*, by fitting the temperature dependence of the FFT amplitudes to the Lifshitz-Kosevich formula for temperature damping; this fit is plotted in Fig. 3c.
The temperature damping coefficient, RT, may be written as35:
where
in which e is the elementary charge, â is the reduced Planck constant, kB is the Boltzmann constant, T is the temperature, and B is the average magnetic field strength of the inverse field range used to compute the FFTs. Thus m* can be found by fitting the quantum oscillatory amplitude to Eqn. (1) as a function of temperature.
Evaluation of the density of states at the Fermi level
The quasiparticle density of states at the Fermi level, g(EF), may be expressed69 in terms of the linear specific heat coefficient extrapolated from the normal state, γN, as
We can compare this with the density of states predicted for a given Fermi surface geometry as measured by the dHvA effect35. For a Fermi surface section with surface element d\({{{{{{{\mathcal{S}}}}}}}}\), which hosts quasiparticles of effective mass m* that have Fermi velocity vF, we can write
For the simple geometric case of a cylindrical Fermi surface, of radius \({k}_{{{{{{{{\rm{F}}}}}}}}}=\sqrt{{k}_{x}^{2}+{k}_{y}^{2}}\) and height kz, combining these two expressions gives a contribution (per cylinder) to the linear specific heat coefficient of
for a metal of molar volume V. Therefore, comparing this simple case with our dHvA data, for the 3.5 kT quantum oscillatory frequency observed for magnetic field applied along the \(\overrightarrow{c}\) direction (Figs. 2, 3), assuming that both cylinders have the same (single) effective mass we found in our Lifshitz-Kosevich temperature study (Fig. 3) we can estimate a contribution per (circular) cylinder to γN ofâââ51.5 mJ molâ1 Kâ2.
Taking this result for the simple case of ideal, circularly cross-sectional cylindrical Fermi surface sections with no kz warping, we then numerically calculated the actual surface area of the squircular, warped Fermi sheets we generated in our Fermi surface simulations (Fig. 4). We found that the electron-type section possesses a surface area 1.20 times bigger than the case of the simple circular cross-sectional cylinder of the same cross-sectional area normal to \(\overrightarrow{c}\) (corresponding to a dHvA frequency of 3.5 kT). For the hole-type sheet, we found that its surface area is bigger by a factor of 1.14. Therefore, we obtained values ofâââ61.8 mJ molâ1 Kâ2 for the electron sheet andâââ58.7 mJ molâ1 Kâ2 for the hole sheet, giving a total contribution to γN ofâââ120.5 mJ molâ1 Kâ2.
We note that this treatment is only approximate, as we assume a constant vF along the entire surface of the Fermi sheets, using the value obtained from our measurement of m* for magnetic field applied along the \(\overrightarrow{c}\) direction. At inclined angles of magnetic field tilt angle a range of effective masses is observed, from as low as 32 me reported in ref. 33, up to the mass of 78 me we find in Supplementary Fig. S19. Therefore, this comparison between the density of states implied by specific heat capacity measurements and inferred from observations of the dHvA effect is only approximate, in the absence of a full determination of the profile of vF along the entire surface of the Fermi surface sections. Further measurements, to carefully determine the evolution in m* as a function of rotation angle in the \(\overrightarrow{c}\)-\(\overrightarrow{a}\) and \(\overrightarrow{c}\)-\(\overrightarrow{b}\) rotation planes, would thus be of great value in minimising the uncertainty of this calculation.
In our calculations we used the cyclotron effective mass measured for magnetic field parallel to the axis of the cylinders, as at this orientation the quantum oscillatory amplitude is largest and thus we are sampling a large proportion of the variation of vF along the cylindersâ surfaces. Notably, a special aspect of our Fermi surface model is that all cross-sectional orbits are extremal for field oriented along \(\overrightarrow{c}\), which likely explains why the signal amplitude is so much larger here. Given the close correspondence between the values of γN measured by specific heat experiments and calculated from our dHvA data and Fermi surface simulations, this adds strong confidence to our Fermi surface simulations and interpretation of the dHvA data that these two quasi-2D sections likely comprise the only Fermi surface sheets present in UTe2.
Data availability
The datasets supporting the findings of this study are available from the University of Cambridge Apollo Repository70.
References
Bardeen, J., Cooper, L. N. & Schrieffer, J. R. Theory of Superconductivity. Phys. Rev. 108, 1175â1204 (1957).
Monthoux, P., Pines, D. & Lonzarich, G. G. Superconductivity without phonons. Nature 450, 1177â1183 (2007).
Leggett, A. J. A theoretical description of the new phases of liquid 3He. Rev. Mod. Phys. 47, 331â414 (1975).
Schemm, E., Gannon, W., Wishne, C., Halperin, W. & Kapitulnik, A. Observation of broken time-reversal symmetry in the heavy-fermion superconductor UPt3. Science 345, 190â193 (2014).
Saxena, S. S. et al. Superconductivity on the border of itinerant-electron ferromagnetism in UGe2. Nature 406, 587â592 (2000).
Nayak, C., Simon, S. H., Stern, A., Freedman, M. & Das Sarma, S. Non-abelian anyons and topological quantum computation. Rev. Mod. Phys. 80, 1083â1159 (2008).
Mackenzie, A. P. & Maeno, Y. The superconductivity of Sr2RuO4 and the physics of spin-triplet pairing. Rev. Mod. Phys. 75, 657â712 (2003).
Pustogow, A. et al. Constraints on the superconducting order parameter in Sr2RuO4 from oxygen-17 nuclear magnetic resonance. Nature 574, 72â75 (2019).
Ran, S. et al. Nearly ferromagnetic spin-triplet superconductivity. Science 365, 684â687 (2019).
Nakamine, G. et al. Anisotropic response of spin susceptibility in the superconducting state of UTe2 probed with 125TeâââNMR measurement. Phys. Rev. B 103, L100503 (2021).
Aoki, D. et al. Unconventional superconductivity in UTe2. J. Phys. Condens. Matter 34, 243002 (2022).
Ran, S. et al. Extreme magnetic field-boosted superconductivity. Nat. Phys. 15, 1250â1254 (2019).
Jiao, L. et al. Chiral superconductivity in heavy-fermion metal UTe2. Nature 579, 523â527 (2020).
Hayes, I. M. et al. Multicomponent superconducting order parameter in UTe2. Science 373, 797â801 (2021).
Wei, D. S. et al. Interplay between magnetism and superconductivity in UTe2. Phys. Rev. B 105, 024521 (2022).
Ishihara, K. et al. Chiral superconductivity in UTe2 probed by anisotropic low-energy excitations. arXiv:2105.13721. (2021).
Bae, S. et al. Anomalous normal fluid response in a chiral superconductor UTe2. Nat. Commun. 12, 2644 (2021).
Sundar, S. et al. Coexistence of ferromagnetic fluctuations and superconductivity in the actinide superconductor UTe2. Phys. Rev. B 100, 140502 (2019).
Xu, Y., Sheng, Y. & Yang, Y.-f Quasi-Two-Dimensional Fermi Surfaces and Unitary Spin-Triplet Pairing in the Heavy Fermion Superconductor UTe2. Phys. Rev. Lett. 123, 217002 (2019).
Ishizuka, J., Sumita, S., Daido, A. & Yanase, Y. Insulator-Metal Transition and Topological Superconductivity in UTe2 from a First-Principles Calculation. Phys. Rev. Lett. 123, 217001 (2019).
Shaffer, D. & Chichinadze, D. V. Chiral superconductivity in UTe2 via emergent C4 symmetry and spin-orbit coupling. Phys. Rev. B 106, 014502 (2022).
Machida, K. Theory of Spin-polarized Superconductors â An Analogue of Superfluid 3He A-phase. J. Phys. Soc. Jpn. 89, 033702 (2020).
Yarzhemsky, V. & Teplyakov, E. Time reversal symmetry and the structure of Cooper pair wavefunction in topological superconductor UTe2. Phys. Lett. A 384, 126724 (2020).
Kittaka, S. et al. Orientation of point nodes and nonunitary triplet pairing tuned by the easy-axis magnetization in UTe2. Phys. Rev. Research 2, 032014 (2020).
Shishidou, T., Suh, H. G., Brydon, P. M. R., Weinert, M. & Agterberg, D. F. Topological band and superconductivity in UTe2. Phys. Rev. B 103, 104504 (2021).
Machida, K. Nonunitary triplet superconductivity tuned by field-controlled magnetization: URhGe, UCoGe, and UTe2. Phys. Rev. B 104, 014514 (2021).
Eo, Y. S. et al. c-axis transport in UTe2: Evidence of three-dimensional conductivity component. Phys. Rev. B 106, L060505 (2022).
Shick, A. B., Fujimori, S.-i & Pickett, W. E. UTe2: A nearly insulating half-filled \(j=\frac{5}{2}5{f}^{3}\) heavy-fermion metal. Phys. Rev. B 103, 125136 (2021).
Moriya, Y., Matsushita, T., Yamada, M. G., Mizushima, T. & Fujimoto, S. Intrinsic Anomalous Thermal Hall Effect in the Unconventional Superconductor UTe2. J. Phys. Soc. Jpn. 91, 094710 (2022).
Choi, H. C., Lee, S. H. & Yang, B.-J. Correlated normal state fermiology and topological superconductivity in UTe2. arXiv:2206.04876. (2022).
Fujimori, S.-i et al. Electronic Structure of UTe2 Studied by Photoelectron Spectroscopy. J. Phys. Soc. Jpn. 88, 103701 (2019).
Miao, L. et al. Low Energy Band Structure and Symmetries of UTe2 from Angle-Resolved Photoemission Spectroscopy. Phys. Rev. Lett. 124, 076401 (2020).
Aoki, D. et al. First Observation of the de Haas-van Alphen Effect and Fermi Surfaces in the Unconventional Superconductor UTe2. J. Phys. Soc. Jpn. 91, 083704 (2022).
Sato, M. & Ando, Y. Topological superconductors: a review. Rep. Prog. Phys. 80, 076501 (2017).
Shoenberg, D. Magnetic Oscillations in Metals. (Cambridge University Press, Cambridge, UK, Cambridge, UK, 1984).
Sakai, H. et al. Single crystal growth of superconducting UTe2 by molten salt flux method. Phys. Rev. Materials 6, 073401 (2022).
Aoki, D. et al. Unconventional Superconductivity in Heavy Fermion UTe2. J. Phys. Soc. Jpn. 88, 43702 (2019).
Rosa, P. F. S. et al. Single thermodynamic transition at 2 K in superconducting UTe2 single crystals. Commun. Mater. 3, 33 (2022).
Weiland, A., Thomas, S. & Rosa, P. Investigating the limits of superconductivity in UTe2. J. Phys. Mater. 5, 044001 (2022).
Matsumura, H. et al. Large Reduction in the a-axis Knight Shift on UTe2 with Tc = 2.1 K. J. Phys. Soc. Jpn. 92, 063701 (2023).
Ajeesh, M. et al. The fate of time-reversal symmetry breaking in UTe2. arXiv:2305.00589. (2023).
Suetsugu, S. et al. Fully gapped pairing state in spin-triplet superconductor UTe2. arXiv:2306.17549. (2023).
Azari, N. et al. Absence of Spontaneous Magnetic Fields Due to Time-Reversal Symmetry Breaking in Bulk Superconducting UTe2. arXiv:2308.09773. (2023).
Scheidemantel, T. J., Ambrosch-Draxl, C., Thonhauser, T., Badding, J. V. & Sofo, J. O. Transport coefficients from first-principles calculations. Phys. Rev. B 68, 125210 (2003).
Ziman, J. M. Principles of the Theory of Solids (Cambridge University Press, 1972).
Metz, T. et al. Point-node gap structure of the spin-triplet superconductor UTe2. Phys. Rev. B 100, 220504 (2019).
Fujibayashi, H. et al. Superconducting Order Parameter in UTe2 Determined by Knight Shift Measurement. J. Phys. Soc. Jpn. 91, 043705 (2022).
Iguchi, Y. et al. Microscopic Imaging Homogeneous and Single Phase Superfluid Density in UTe2. Phys. Rev. Lett. 130, 196003 (2023).
Anderson, P. W. & Morel, P. Generalized Bardeen-Cooper-Schrieffer States and the Proposed Low-Temperature Phase of Liquid He3. Phys. Rev. 123, 1911â1934 (1961).
Miyake, K. Theory of Pairing Assisted Spin Polarization in Spin-Triplet Equal Spin Pairing: Origin of Extra Magnetization in Sr2RuO4 in Superconducting State. J. Phys. Soc. Jpn. 83, 053701 (2014).
Aoki, D., Ishida, K. & Flouquet, J. Review of U-based Ferromagnetic Superconductors: Comparison between UGe2, URhGe, and UCoGe. J. Phys. Soc. Jpn. 88, 022001 (2019).
Dong, J. et al. Competing orders and spin-density-wave instability in La(O1âxFx)FeAs. EPL 83, 27006 (2008).
Carrington, A. Quantum oscillation studies of the Fermi surface of iron-pnictide superconductors. Rep. Prog. Phys. 74, 124507 (2011).
Sebastian, S. E. et al. Compensated electron and hole pockets in an underdoped high-Tc superconductor. Phys. Rev. B 81, 214524 (2010).
Proust, C. & Taillefer, L. The remarkable underlying ground states of cuprate superconductors. Annu. Rev. Condens. Matter Phys. 10, 409â429 (2019).
Maehira, T., Hotta, T., Ueda, K. & Hasegawa, A. Electronic Structure and the Fermi Surface of PuCoGa5 and NpCoGa5. Phys. Rev. Lett. 90, 207007 (2003).
Lee, P. A., Nagaosa, N. & Wen, X.-G. Doping a Mott insulator: Physics of high-temperature superconductivity. Rev. Mod. Phys. 78, 17â85 (2006).
Niu, Q. et al. Evidence of Fermi surface reconstruction at the metamagnetic transition of the strongly correlated superconductor UTe2. Phys. Rev. Res. 2, 033179 (2020).
Helm, T. et al. Suppressed magnetic scattering sets conditions for the emergence of 40 T high-field superconductivity in UTe2 arXiv:2207.08261. (2022).
Yamaji, K. On the Angle Dependence of the Magnetoresistance in Quasi-Two-Dimensional Organic Superconductors. J. Phys. Soc. Jpn. 58, 1520â1523 (1989).
Honda, F. et al. Pressure-induced Structural Phase Transition and New Superconducting Phase in UTe2. J. Phys. Soc. Jpn. 92, 044702 (2023).
Broyles, C. et al. Revealing a 3D Fermi Surface Pocket and Electron-Hole Tunneling in UTe2 with Quantum Oscillations. Phys. Rev. Lett. 131, 036501 (2023).
Aoki, D. Invited talk. SCES 2023, Incheon, South Korea (2023).
Kartsovnik, M. V. et al. Direct Observation of the Magnetic-Breakdown Induced Quantum Interference in the Quasi-Two-Dimensional Organic Metal κâââ(BEDTâTTF)2CU(NCS)2. Phys. Rev. Lett. 77, 2530â2533 (1996).
Harrison, N. et al. Magnetic breakdown and quantum interference in the quasi-two-dimensional superconductor κ-(BEDT-TTF)2Cu(NCS)2 in high magnetic fields. J. Phys. Condens. Matter 8, 5415 (1996).
Weinberger, T. I. et al. Quantum interference between quasi-2D Fermi surface sheets in UTe2 arXiv:2307.00568. (2023).
Haga, Y. et al. Purification of Uranium Metal using the Solid State Electrotransport Method under Ultrahigh Vacuum. Jpn. J. Appl. Phys. 37, 3604 (1998).
Cleveland, W. S. & Devlin, S. J. Locally weighted regression: An approach to regression analysis by local fitting. J. Am. Stat. Assoc. 83, 596â610 (1988).
Ashcroft, N. W. & Mermin, D. N. Solid State Physics. (Harcourt College Publishers, San Diego, CA, 1976).
Eaton, A. G. et al. Supporting data for: Quasi-2D Fermi surface in the anomalous superconductor UTe2 [https://doi.org/10.17863/CAM.96845] (2023).
Weinberger, T. Empirical_FS, GitHub repository [https://github.com/TheoWeinberger/Empirical_FS] (2023).
Weinberger, T. py_FS, GitHub repository [https://github.com/TheoWeinberger/py_FS] (2023).
Acknowledgements
We are grateful to J. Chen, D. Shaffer, D.V. Chichinadze, S.S. Saxena, C.K. de Podesta, P. Coleman, T. Helm, A.B. Shick, and W. Luo for fruitful discussions. We thank T.J. Brumm, T.P. Murphy, A.F. Bangura, D.E. Graf, S.T. Hannahs, S.W. Tozer, E.S. Choi, and L. Jiao for technical advice and assistance. This project was supported by the EPSRC of the UK (grant no. EP/X011992/1). A portion of this work was performed at the National High Magnetic Field Laboratory, which is supported by National Science Foundation Cooperative Agreement No. DMR-1644779 and the State of Florida. Crystal growth and characterization were performed in MGML (mgml.eu), which is supported within the program of Czech Research Infrastructures (project no. LM2023065). We acknowledge financial support by the Czech Science Foundation (GACR), project No. 22-22322S. The EMFL supported dual-access to facilities at MGML, Charles University, Prague, under the European Unionâs Horizon 2020 research and innovation programme through the ISABEL project (No. 871106). A part of this work was also supported by the JAEA REIMEI Research Program. DFT calculations were performed using resources provided by the Cambridge Service for Data Driven Discovery (CSD3) operated by the University of Cambridge Research Computing Service (https://www.csd3.cam.ac.uk), provided by Dell EMC and Intel using Tier-2 funding from the Engineering and Physical Sciences Research Council (capital grant EP/T022159/1), and DiRAC funding from the Science and Technology Facilities Council (https://www.dirac.ac.uk). This project also used the ARCHER2 UK National Supercomputing Service (https://www.archer2.ac.uk). A.G.E. acknowledges support from a QuantEmX grant from ICAM and the Gordon and Betty Moore Foundation through Grant GBMF5305; from the Henry Royce Institute for Advanced Materials through the Equipment Access Scheme enabling access to the Advanced Materials Characterisation Suite at Cambridge, grant numbers EP/P024947/1, EP/M000524/1 & EP/R00661X/1; and from Sidney Sussex College (University of Cambridge). T.I.W. and A.J.H. acknowledge support from EPSRC studentships (EP/R513180/1 & EP/M506485/1).
Author information
Authors and Affiliations
Contributions
A.G.E., N.J.M.P., A.J.H., R.N. & S.M.B. performed quantum oscillation measurements. P.O., H.S., Y.H. & M.V. optimised the molten salt flux growth technique of UTe2. High-quality single crystal specimens were grown and characterised by A.C., JiÅÃ P., Jan P., T.H., G.B., V.S. & M.V. Quantum oscillation data were analysed by A.G.E., N.J.M.P. & A.J.H., and subsequently interpreted by A.G.E., T.I.W., N.J.M.P., Z.W., A.J.H., G.G.L. & F.M.G. Fermi surface modelling and DFT calculations were performed by T.I.W., who also modelled the electrical conductivity tensor with guidance from G.G.L. The manuscript was written by A.G.E. & T.I.W., with input from all co-authors. A.G.E. conceived and oversaw the project.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Peer review
Peer review information
Nature Communications thanks the anonymous, reviewer(s) for their contribution to the peer review of this work. A peer review file is available.
Additional information
Publisherâs note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the articleâs Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the articleâs Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Eaton, A.G., Weinberger, T.I., Popiel, N.J.M. et al. Quasi-2D Fermi surface in the anomalous superconductor UTe2. Nat Commun 15, 223 (2024). https://doi.org/10.1038/s41467-023-44110-4
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41467-023-44110-4
This article is cited by
-
Orphan high field superconductivity in non-superconducting uranium ditelluride
Nature Communications (2024)
-
Field-induced compensation of magnetic exchange as the possible origin of reentrant superconductivity in UTe2
Nature Communications (2024)
-
Single-component superconductivity in UTe2 at ambient pressure
Nature Physics (2024)
-
Correlation-driven topological Kondo superconductors
Communications Physics (2024)
-
Correlated normal state fermiology and topological superconductivity in UTe2
Communications Physics (2024)