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.

Fig. 1: Characterisation of high purity UTe2.
figure 1

a Specific heat capacity (Cp) divided by temperature (T) of a UTe2 single crystal measured on warming to 3 K. A single, sharp bulk superconducting transition is exhibited. (Inset) The crystal structure of UTe2. b Resistivity (ρ) versus temperature squared up to 4 K for current sourced along the \(\overrightarrow{a}\) direction. A superconducting transition temperature of 2.1 K is observed (defined by zero resistivity, i.e. below the detection limit of 0.01 μΩ cm). A residual resistivity ratio (RRR, defined in the text) of 900 is found, with a residual resistivity ρ0 ≲ 0.5 μΩ cm, indicative of very high sample purity36,38. (Inset) The same dataset as the main panel extended up to 300 K.

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.

Fig. 2: Angular evolution of quantum oscillatory frequencies and amplitude.
figure 2

a Oscillatory component of magnetic torque (Δτ) at various angles and b, their corresponding Fourier frequency spectra. Angles were calibrated to within 2∘ of uncertainty. 0∘ corresponds to magnetic field, \(\overrightarrow{H}\), applied along the \(\overrightarrow{c}\) direction; 90∘ corresponds to field applied along the \(\overrightarrow{a}\) direction. At 0∘ a singular frequency peak of high amplitude is observed at f = 3.5 kT (with a second harmonic peak at 2f = 7.0 kT). Upon rotating away from \(\overrightarrow{c}\) towards \(\overrightarrow{a}\) this peak splits and the oscillatory amplitude diminishes markedly. c Raw magnetic torque signal (without background subtraction), τ, at 0∘ and 74∘. Both curves have been translated to have the same value of τ at 26 T for comparison. A clear oscillatory component is visible in the raw torque signal of the 0∘ data. In comparison, despite the 74∘ torque being of greater overall magnitude, it is markedly smoother than that at 0∘. d, Δτ at θ = 74∘ and e, the corresponding fast Fourier transform (FFT). Despite the very small oscillatory amplitude, a sharp high frequency peak of f = 18.5 kT is clearly resolved on top of the background noise profile.

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

Fig. 3: Singular quantum oscillatory frequency and heavy effective cyclotron mass along the \(\overrightarrow{c}\) direction.
figure 3

a Oscillatory component of magnetic torque for field applied very close (within 2∘) to the \(\overrightarrow{c}\) direction and b corresponding FFT amplitudes for incremental temperatures as indicated. A dominant frequency component of 3.5 kT (and second harmonic at 7.0 kT) is observed. c FFT amplitude as a function of temperature (coloured points). The solid line is a fit to the Lifshitz-Kosevich formula35, which fits the data well and yields a heavy effective cyclotron mass, m*, of 41(2) me, where me is the bare electron mass.

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

Fig. 4: The Fermi surface of UTe2.
figure 4

Angular dependence of the dHvA frequencies for a the \(\overrightarrow{c}\) to \(\overrightarrow{a}\) rotation plane, and b the \(\overrightarrow{c}\) to \(\overrightarrow{b}\) rotation plane. Blue and gold symbols are simulated frequencies (see Supplementary Information) for extremal orbit areas of our calculated electron-type (e−) and hole-type (h+) Fermi surface sections, respectively; red and purple symbols represent dHvA data from this study and ref. 33. c Side- and top-view of our simulated Fermi surface cylinders, with high symmetry points indicated. Again, blue (gold) represents electron- (hole-) type sections. d Extended-zone view of the UTe2 Fermi surface.

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

Fig. 5: Dependency of electrical conductivity on the warping of the Fermi surface cylinders.
figure 5

a Diagonal components of the electrical conductivity tensor determined from the Fermi surface fitting (see Supplementary Information for calculation details). σxx,yy,zz respectively refer to conductivity in the \(\overrightarrow{a}\), \(\overrightarrow{b}\), \(\overrightarrow{c}\) directions. 100% warping corresponds to our Fermi surface model of UTe2, which we find best fits the dHvA data (Fig. 4). b The same conductivity components as panel a plotted here on a logarithmic scale. For 0% warping (no undulation) we would expect σzz to be orders of magnitude lower than σxx and σyy. However, for our Fermi surface model (at 100% warping) each of σxx, σyy, and σzz are within an order of magnitude. c, Fermi surface simulations for warping factors between 0 and 200%.

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:

$${R}_{T} \,=\, \frac{X}{\sinh X}$$
(1)

where

$$X=\frac{2{\pi }^{2}{k}_{{{{{{{{\rm{B}}}}}}}}}T{m}^{*}}{e\hslash B},$$
(2)

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

$$g({E}_{{{{{{{{\rm{F}}}}}}}}})=\frac{3{\gamma }_{{{{{{{{\rm{N}}}}}}}}}}{{\pi }^{2}{k}_{{{{{{{{\rm{B}}}}}}}}}^{2}}.$$
(3)

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

$$g({E}_{{{{{{{{\rm{F}}}}}}}}})=\frac{1}{4{\pi }^{3}\hslash }\int\frac{d{{{{{{{\mathcal{S}}}}}}}}}{{v}_{{{{{{{{\rm{F}}}}}}}}}}.$$
(4)

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

$${\gamma }_{{{{{{{{\rm{N}}}}}}}}}=\frac{{k}_{{{{{{{{\rm{B}}}}}}}}}^{2}V{m}^{*}{k}_{z}}{6{\hslash }^{2}},$$
(5)

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.