Quantum interference between quasi-2D Fermi surface sheets in UTe₂

Abstract

UTe₂ is a spin-triplet superconductor candidate for which high quality samples with long mean free paths have recently become available, enabling quantum oscillation measurements to probe its Fermi surface and effective carrier masses. It has recently been reported that UTe₂ possesses a 3D Fermi surface component [Phys. Rev. Lett. 131, 036501 (2023)]. The distinction between 2D and 3D Fermi surface sections in triplet superconductors can have important implications regarding the topological properties of the superconductivity. Here we report the observation of oscillatory components in the magnetoconductance of UTe₂ at high magnetic fields. We find that these oscillations are well described by quantum interference between quasiparticles traversing semiclassical trajectories spanning magnetic breakdown networks. Our observations are consistent with a quasi-2D model of this material’s Fermi surface based on prior dHvA-effect measurements. Our results strongly indicate that UTe₂ – which exhibits a multitude of complex physical phenomena – possesses a remarkably simple Fermi surface consisting exclusively of two quasi-2D cylindrical sections.

Young’s double slit experiment represents a powerful example of the wave-particle duality of photons [1]. A century later Davisson and Germer observed a similar phenomenon involving the quantum mechanical interference of a beam of electrons incident on a crystalline target [2, 3]. In the solid state, superconducting quantum interference devices provide exceptionally accurate measurements of magnetic flux via diffraction-modulated interferometry [4, 5]. For the case of normal metals, the manifestation of quantum interference (QI) effects in the magnetoconductance was first predicted by Shiba and Fukuyama [6], and soon thereafter experimentally realized by Stark and Friedberg in their measurements of the magnetoresistance of magnesium [7]. The concept of the Stark interferometer is based on interference between semiclassical quasiparticle trajectories across magnetic breakdown networks connecting separate Fermi surface (FS) sections, yielding oscillations in the conductivity that are periodic in inverse magnetic field strength [8, 9, 10, 11].

Since the seminal experiments by Stark and coworkers, quantum interference oscillations (QIOs) have been observed in a variety of materials [12, 13, 14, 15, 16] including, in particular, a number of organic metals with quasi-2D (Q2D) FSs [17, 18, 19, 20, 21, 22, 23, 24]. Unlike quantum oscillations (QOs) from the dHvA- or SdH-effects, in which phase coherence and Landau quantization of quasiparticles traversing orbits corresponding to closed FS sections provide a direct measurement of the FS [11], QIOs only yield an indirect probe of the FS, as their frequencies correspond to $k$-space orbits spanning separate FS sections. Therefore, QIOs are only observed in materials in which the $k$-space separation of FS sections is sufficiently small for quasiparticles to tunnel between FS sheets in accessible magnetic field strengths [25, 11]. It is important to note that QI is exclusively a kinetic effect and is thus observable in the electrical transport – unlike the dHvA-effect, QIOs do not correspond to an oscillatory component of the free energy, therefore QI effects cannot be observed in bulk thermodynamic properties such as the magnetization [9, 10, 17, 26].

Here, we report the observation of QIOs at high magnetic fields in contactless resistivity measurements of the heavy fermion paramagnetic metal UTe₂. This material has recently shown promising signs of being a spin-triplet superconductor [27, 28, 29], similar to the analogous ferromagnetic compounds UGe₂, URhGe and UCoGe [30, 31, 32]. Evidence indicating triplet pairing in UTe₂ comes from a number of sources including small changes in the NMR Knight shift on cooling through the superconducting critical temperature $T_{\text{c}}$ [29, 33] along with anisotropic upper critical fields that far exceed the Pauli limit for singlet pairing [34, 35, 27, 36]. Recent advances in the growth procedure of single crystal UTe₂ specimens have led to a marked enhancement in crystalline quality, enabling the observation of QOs from the dHvA-effect [37, 38]. The angular profile of the dHvA data is indicative of a relatively simple Q2D FS, consisting of one electron-type and one hole-type cylinder, each hosting quasiparticles of heavy effective masses $\sim$ 40 $m_{\text{e}}$ [37, 38].

Methods – UTe₂ single crystals were grown by a molten salt flux technique [39] using the methodology detailed in ref. [38]. This technique has been shown to yield high quality specimens of $T_{\text{c}}\approx$ 2.1 K with long mean free paths of the order of 100 nm [36, 37, 39, 38]. Details regarding sample characterization are given in the Supplemental Materials . Contactless resistivity measurements were performed in static fields to 41.5 T at the National High Magnetic Field Lab, Tallahassee, Florida, using the tunnel diode oscillator (TDO) technique [41]; similar measurements were obtained in pulsed fields to 70 T at the Hochfeld-Magnetlabor, HZDR, Dresden, using the proximity detector oscillator (PDO) technique [42, 43].

Results – Figure 1 shows the background-subtracted TDO signal ($\Delta f_{\text{TDO}}$) for magnetic field oriented 8$\degree$ away from the crystalline $c$-axis towards the $a$-axis ($\theta_{c}=8\degree$). The FFT of the $\Delta f_{\text{TDO}}$ data reveals four clear frequency branches, which we label as $\upalpha$-$\updelta$. Notably, the FFT spectra at $\theta_{c}=8\degree$ of the TDO signal are very different to the spectra we observed in our prior dHvA study at the same angle (ref. [38]), implying that these are not QOs stemming from the SdH-effect. Furthermore, the amplitude of dHvA QOs diminished by almost an order of magnitude between 19 mK and 200 mK – whereas here the signal is large and very well resolved at 400 mK. These observations indicate that the oscillations in $f_{\text{TDO}}$ are likely QIOs not QOs, as QIOs generally correspond to reciprocal space areas constructed from sums and differences between FS sections, and often exhibit effective masses much lower than those of dHvA and SdH QOs [17, 18, 13].

Using our FS model from ref. [38], we illustrate in Figs. 1 & 2 how the frequencies of the $\upalpha$-$\updelta$ FFT peaks correspond remarkably well to $k$-space areas between the cylindrical Fermi sheets, which are centred at the centre and corners of the first Brillouin zone (BZ). Each of these frequency components can thus be well understood as coming from QI between two quasiparticles – one making two orbits around a FS cylinder, and the other traversing a magnetic breakdown (MB) network between two cylinders of the same carrier type.

Figure 1: (a) QIOs in the contactless resistivity of UTe₂ and (b) the corresponding FFT spectra. (c) Our Fermi surface model for UTe₂ adapted from ref. [38], with the planes $k_{z}=0$ and $k_{z}=\frac{\pi}{2c}$ indicated. (d) An extended-zone view, with the $c$-axis into the page, at $k_{z}=0$ and (e) at $k_{z}=\frac{\pi}{2c}$. QI trajectories enclosing areas that correspond to the QIO frequencies $\upalpha$-$\updelta$ are indicated; note that each enclosed area $\upalpha$-$\updelta$ has two distinct MB networks corresponding to it (denoted as $\lambda_{1,2}$).

Figure 2: Schematic of two semiclassical quasiparticle trajectories that interfere to give the $f_{\upbeta}$ frequency component. The difference in area of the paths ACDEA and ABABA is equal to $\mathcal{A}_{\upbeta}$, as shown in the text. Note that the orbit CDEBC is not possible due to the direction of the Lorentz force (indicated with arrows). The frequency components $\upalpha$, $\upgamma$ and $\updelta$ come from QI between analogous networks, as traced out in Figure 1.

Figure 3: QIOs, FFTs, and apparent effective masses for (a-c) steady field measurements focussing on the $\upalpha$ and $\upbeta$ frequency components, and (b-e) higher temperature pulsed field measurements focussing on the $\upgamma$ and $\updelta$ components. The effective masses $m^{*}_{\upgamma}$ and $m^{*}_{\updelta}$ are markedly lower than those observed in dHvA-effect measurements for the same field orientation [38].

Figure 4: (a) Evolution of FFT spectra as a function of magnetic field tilt angle away from the $c$-axis. (b) Frequency versus angle for the $\upalpha$ and $\upbeta$ branches (in purple) and the $\upgamma$ and $\updelta$ branches (in green). (c) The expectation of the angular frequency profile corresponding to the areas $\mathcal{A}_{\upalpha,\upbeta,\upgamma,\updelta}$ computed from our FS model depicted in Fig. 1. Surprisingly good correspondence between model prediction and experimental measurement is observed, given the simplicity of the model assumptions (given in the Supplementary Materials).

Figure 5: (a) Raw TDO signal for magnetic field oriented along the $c$-axis ($\theta_{c}=0\degree$) and along the $a$-axis ($\theta_{c}=90\degree$). (b) QIOs at $\theta_{c}$ angles as indicated over 20-41.5 T and (c) their corresponding FFTs. (d) High frequency oscillations at $\theta_{c}=90\degree$ are quickly suppressed by rotation away from the $a$-axis. (e) FFTs of the curves in panel d. Only the data at $\theta_{c}=90\degree$ has a resolvable frequency component above the noise floor of the measurement. (f) A side-view of the cylindrical sheets of our Q2D FS model (compared to the axial-view given in Fig. 1). The trajectories $\upzeta$ and $\upeta$ are identified, which enclose areas with very good correspondence to the oscillatory frequencies observed for magnetic field orientated along the $a$-axis.

To show this, we consider the generalized theory of MB orbits given by Kaganov and Slutskin [26]. In a magnetic field $B$ the oscillatory component of a kinetic coefficient, such as the electrical conductivity, is composed of harmonics of the form

$$\sum_{\lambda,\lambda^{{}^{\prime}}}\exp[i\left(\phi_{\lambda}-\phi_{\lambda^{{}^{\prime}}}\right)]=\sum_{\lambda,\lambda^{{}^{\prime}}}\exp\left(\frac{i\hbar}{eB}\mathcal{A}_{\lambda,\lambda^{{}^{\prime}}}\right)$$ (1)

for phase $\phi$ where $\lambda$ and $\lambda^{{}^{\prime}}$ are the two semiclassical quasiparticle paths that share a common start and end point, enclosing between them an area in reciprocal space of $\mathcal{A}_{\lambda,\lambda^{{}^{\prime}}}$ [18].

Take for example the area $\mathcal{A_{\upbeta}}$ shaded in Fig. 2, which sits at the corner of the first BZ for $k_{z}=\frac{\pi}{2c}$. Writing the area of the hole-type FS cylinder as $\mathcal{A}_{h^{+}}$, we can see that the area $\mathcal{A_{\upbeta}}$ is equal to the difference of the areas enclosed by the paths $\lambda$ = ACDEA and $\lambda^{{}^{\prime}}$ = ABABA as $\mathcal{A_{\text{ACDEA}}}-\mathcal{A}_{\text{ABABA}}=(2\mathcal{A}_{h^{+}}+\mathcal{A}_{\upbeta})-2\mathcal{A}_{h^{+}}=\mathcal{A_{\upbeta}}$. Similarly, areas corresponding to the $\upalpha$, $\upgamma$ and $\updelta$ frequency components are formed by QI between the quasiparticle trajectories traced in Fig. 1 . The probability of a quasiparticle traversing a path depends on the number of MB tunnelling events (each of probability amplitude $p$) and Bragg reflections (of probability amplitude $q$) that are contained within the path, where $\absolutevalue{p}^{2}=P=\exp(\nicefrac{{-B_{0}}}{{B}})$ for breakdown field $B_{0}$ and $q=i\sqrt{\left(1-P\right)}$ [11, 25, 17, 47, 18]. Therefore, the probability amplitudes for quasiparticles to traverse the paths $\lambda$ = ACDEA and $\lambda^{{}^{\prime}}$ = ABABA, corresponding to the $\upbeta$ frequency in Fig. 2, are $q^{4}p^{4}\exp(i\phi_{\lambda})$ and $q^{8}\exp(i\phi_{\lambda^{{}^{\prime}}})$, respectively. Due to this exponentially suppressed tunnelling probability – which necessitates the application of high magnetic fields – we limit our discussion just to the lowest order relevant networks as depicted in Fig. 1, each of which requires only 4 instances of MB.

By Eqn. 1, the probability of quasiparticles traversing the paths in Fig. 2 will involve oscillating terms including some proportional to $\cos\left[{\phi_{\lambda}-\phi_{\lambda^{{}^{\prime}}}}\right]=\cos\left[{2\pi(2f_{h^{+}}+f_{\upbeta}-2f_{h^{+}})/B}\right]=\cos\left[{2\pi f_{\upbeta}/B}\right]$, which will contribute to the (real part of the) conductivity. Furthermore, in the low temperature limit the temperature dependence of QIOs simply follows the Lifshitz-Kosevich theory [48, 11, 13] with an apparent effective mass $m^{*}_{\lambda,\lambda^{\prime}}$, which is proportional to the dependence of the phase on the electron energy, $E_{k}$:

$$m^{*}_{\lambda,\lambda^{\prime}}=\frac{e\hbar B}{2\pi}\left|\frac{\partial\left(\phi_{\lambda}-\phi_{\lambda^{{}^{\prime}}}\right)}{\partial E_{k}}\right|=|m_{\lambda}^{*}-m_{\lambda^{{}^{\prime}}}^{*}|$$ (2)

where $m_{\lambda}^{*}$ ($m_{\lambda^{{}^{\prime}}}^{*}$) denotes the effective mass of path $\lambda$ ($\lambda^{{}^{\prime}}$) [18, 13]. Note that it is the difference in the effective masses of the two interfering paths that determines the apparent effective mass of QIOs – thus enabling QIOs to be observed to much higher temperatures than QOs from the dHvA- and SdH-effects [17, 18, 13, 19].

Figure 3 shows that the $\upgamma$ and $\updelta$ frequencies in the QIO spectra of UTe₂ exhibit apparent effective masses ($\approx$ 5 $m_{\text{e}}$) almost an order of magnitude lower than the quasiparticle effective masses reported for dHvA QOs ($\sim$ 40 $m_{\text{e}}$) [37, 38], showing that the subtraction of masses between the two trajectories in Eqn. 2 has almost cancelled out. By contrast, the $\upalpha$ and $\upbeta$ frequencies are much heavier with masses in the region of 20–35 $m_{e}$ (Fig. 3 and Fig. S3 in the Supplementary Materials). This implies that these MB networks span FS sections with a highly anisotropic distribution of the Fermi velocity, $v_{\text{F}}$. This is consistent with several experimental [49, 37, 38] and theoretical [50, 51] studies that indicate the hybridization between U $f$-electrons with the U $d$-bands and Te $p$-bands, which provides the dominant contribution to the Q2D FS, can result in significant variations in the effective quasiparticle masses at points around the cylindrical sheets. We note that our uncertainty in $m^{*}_{\upalpha}$ and $m^{*}_{\upbeta}$ is considerably larger than for $m^{*}_{\upgamma}$ and $m^{*}_{\updelta}$ due to these frequencies only being observable near the base temperature of the ³He cryostat used for this measurement, with the uncertainty in temperature dominating the uncertainty in $m^{*}_{\upalpha,\upbeta}$. Further measurements in the experimentally challenging temperature–field regime of $\leq$ 200 mK and $\geq$ 40 T are required to carefully probe the anisotropy of $v_{\text{F}}$ around the FS of UTe₂, and thus to better understand the hybridization of the $f$, $d$ and $p$ bands.

In principle, an infinite number of MB networks could give rise to QIOs. Thus, it is expected that orbits of the type $\mathcal{A_{\text{ACDEA}}}-\mathcal{A}_{\text{ABA}}=(2\mathcal{A}_{h^{+}}+\mathcal{A}_{\upbeta})-\mathcal{A}_{h^{+}}=\mathcal{A_{\upbeta}}+\mathcal{A}_{h^{+}}$ should occur. However, the effective mass associated with these orbits would be greater than the masses of the hole and electron orbits from which they arise. If in the most simple case we assume that the breakdown orbits of type $\mathcal{A_{\text{ACDEA}}}$ have masses of $2m^{*}_{h^{+}/e^{-}}+\epsilon_{m}$, where $\epsilon_{m}$ is a small difference to account for the fact that quasiparticles are in fact not traversing full FS sheets, then by Eqn. 2 these breakdown orbits need to interfere with two full FS sheet orbits to produce oscillations of $m^{*}=\epsilon_{m}$. By comparison, orbits of the type $\mathcal{A_{\text{ACDEA}}}-\mathcal{A}_{\text{ABA}}$ would instead have masses of $m^{*}=m^{*}_{h^{+}/e^{-}}+\epsilon_{m}$ and as such would be too heavy to observe at ³He temperatures.

Figure 4 shows the evolution of QIO frequency with magnetic field tilt angle, and compares with the prediction from our Q2D FS model (in panel c). Although this is only a crude approximation of the expected QIO frequency profile, we find remarkably good agreement between our FS model adapted from ref. [38] and the QIOs we observe in TDO measurements. This result gives strong confidence that the FS of UTe₂ is very well described by our Q2D model.

Our discussion so far has focussed on field aligned coaxially to the FS cylinders (along $c$), and at inclination angles close to $c$. Figure 5 shows that for field oriented along the $a$-axis, two additional frequencies $f_{\upzeta}$ = 220 T and $f_{\upeta}$ = 4.5 kT are observed. Again, the enclosed areas of these MB networks correspond very well to our Q2D FS model (Fig. 5f). The low frequency $\upzeta$ oscillations for field along $a$ are of considerable amplitude and are clearly observable in the raw TDO signal without background subtraction (Fig. 5a). Along the $a$-axis $\upzeta$ again corresponds to a QIO, whereas $\upeta$ is consistent with a conventional MB orbit, which may explain its small amplitude as well as its observation for $B$ only directly along the $a$-axis.

We note that a similar study of oscillations in the TDO signal of UTe₂ at high fields was recently reported [52]. For $H\parallel a$ ref. [52] reports an oscillatory frequency of 223 T, in very good agreement with the 220 T $\upzeta$ orbit we observe at this field orientation (Fig. 5). However, rather than being of a QI origin, the authors of ref. [52] interpreted the observed oscillatory waveform to comprise QOs from the SdH-effect caused by the presence of light 3D FS pocket(s). The distinction between Q2D and 3D FS dimensionality in the case of UTe₂ is important, as any 3D pockets could have significant implications regarding the topological properties of the putatively spin-triplet superconductivity [53, 54].

However, in our measurements we do not observe any indication of the presence of a 3D FS pocket. Fig. 5b shows the evolution of $\Delta f_{\text{TDO}}$ as the field is tilted away from $a$ towards $c$. For magnetic field oriented along the $a$-axis we observe low frequency large amplitude oscillations, in good agreement with the raw data presented in ref. [52]. A large oscillatory component is still visible 9$\degree$ away from $a$; however, after a rotation of 20$\degree$ (to $\theta_{c}=70\degree$) no oscillations are observed within the resolution of the measurement. This is inconsistent with this frequency branch coming from SdH-effect QOs due to a 3D pocket; however, this behavior is consistent with a QI interpretation of the oscillatory origin, as the $\upzeta$ trajectory is only possible for $B$ close to $a$. In the Supplementary Materials we show a similar evolution for rotating away from $a$ towards $b$. Furthermore, no slow oscillations at these tilt angles have been reported in prior dHvA measurements by the field modulation [37] or torque magnetometry [38] techniques – they appear only to be observed in the electrical conductivity, again consistent with a QI origin.

The stark difference in the effective masses of the $\upalpha,\upbeta$ and $\upgamma,\updelta$ components implies a strong anisotropy of $v_{\text{F}}(\textbf{k})$. In our recent study of dHvA QOs in UTe₂ we observed two-fold effective mass variations along the measured frequency branches under rotation away from the $c$-axis [38]. In order to attain such a variation, this implies a significant anisotropy of $v_{\text{F}}(k_{z})$, which in turn could account for the large difference in effective masses of the QIOs. Such a variation in effective mass likely stems from substantial hybridization between U $d$-bands and Te $p$-bands, which are the main contributors to the Q2D FS sheets [51], and an $f$-electron band sitting just above the Fermi level. This band has been detected in ARPES measurements, in which a significant spectral weight was observed at the Z-point of the BZ [49]. Models of UTe₂ that include the presence of such a band [50, 51] show that the effect of the U $f$-electrons hybridizing with U $d$-bands is to compress them in energy, effectively increasing their band mass . It is therefore likely that $v_{\text{F}}$ is lowest (and thus $m^{*}$ is highest) at the regions of the FS cylinders that are closest to the Z point, as here the spectral contribution of the $f$-electrons is largest and thus the hybridization with them will be the greatest.

In summary, we measured the contactless resistivity of UTe₂ to high applied magnetic field strengths. We observed oscillatory components that are well explained by quantum interference between semiclassical quasiparticle trajectories spanning magnetic breakdown networks. We find that the quantum interference frequencies correspond very well to a quasi-2D model of the UTe₂ Fermi surface. Our observations give no indication of the presence of any 3D Fermi surface pockets in this material.