Connection between f-electron correlations and magnetic excitations in UTe2

The heavy-fermion paramagnet UTe₂ has been a subject of focused investigation within the strongly correlated electron community since the recent discovery of its highly nontrivial superconducting ground state¹. UTe₂ strongly violates the Pauli limit for BCS superconductors, which relates the temperature of the onset of superconductivity T_c ≈ 2K^2,3 to the critical magnetic field. The Pauli limit is exceeded by factors of approximately 2, 4, and 2.5 for fields along the crystalline (hat{a}), (hat{b}), and (hat{c}) directions respectively, and magnetization scaling strongly suggests proximity to ferromagnetism^1,4 despite no observation of long-ranged magnetic order, which was taken as initial supporting evidence of a spin-triplet superconducting ground state in UTe₂.

Nuclear magnetic resonance (NMR) studies indicate that the spin susceptibility remains constant upon cooling below T_c^5,6,7. This violates expectations for spin-singlet BCS superconductivity, where a Knight shift signifies is associated with the condensation of the superconducting quasiparticles. This, along with other experimental results such as polar Kerr effect measurements^8,9,10, exotic reentrant superconductivity at fields above 35 T^11,12, and fluctuations in μSR^13,14,15, are all taken to be supporting evidence for a potential spin-triplet superconducting ground state in UTe₂. Such a superconducting pairing is exceedingly rare in condensed matter systems, with the only conclusive example being that of the rigorously studied superfluid phase of ³He¹⁶. A spin-triplet superconductor is of particular interest in the context of applications to quantum devices.

Many open questions remain regarding the ground state superconductivity in UTe₂. The superconducting nodal gap function^17,18,19, the order parameters associated with each of the superconducting phases^5,8, time reversal symmetry breaking⁹, the role of magnetic interactions, and the electronic band structure and Fermi surface are all unresolved. As spin-triplet superconductivity can be mediated by ferromagnetism (although an antiferromagnetic mechanism^20,21 is also possible) unambiguous understanding of the spin interactions is required to describe the paramagnetic ground state of UTe₂. Additionally, band hybridization in heavy fermion systems involves the Kondo effect and RKKY interactions, and the delicate interplay of these effects with f-electron band has been linked to unconventional superconductivity and conduction electron mediated magnetic exchange interactions²². Their role in UTe₂ is currently not well understood.

Inelastic neutron scattering (INS) experiments are the most direct probe of spin fluctuations and magnetism in heavy fermion superconductors, many of which possess many similarities to UTe₂. In the cases of UCoGe^23,24 and UGe₂^24,25, ferromagnetic fluctuations are observed, whereas antiferromagnetic correlations associated with superconductivity have been observed in materials like URu₂Si₂^26,27 and the cuprates²⁸. However, the imaginary susceptibility measured by INS ({chi }^{{prime} }({boldsymbol{Q}},hslash omega )) in correlated f-electron systems often arises not only from magnetic excitations like magnons, but from band excitons as in the cases of CePd₃²⁹ and SmB₆³⁰. Further, it has recently been demonstrated that in f-electron systems even the magnons in ordered phases sensitively depend upon the underlying correlated band structure³¹. An applied magnetic field can help to distinguish between magnons and band excitons. For example, in the case of CeRu₂Si₂³², magnetic field strongly suppresses magnon scattering due to interionic correlations, as opposed to the gentle perturbations expected in the band exciton picture.

Earlier INS studies of UTe₂ report a condensation of excitations emerging at the Brillouin zone boundary upon cooling below the Kondo hybridization temperature of T_K ≈ 40K³³. Additional INS results consistently report the presence of a broad continuum of excitations in the superconducting state peaked in intensity at energies of ℏω = 4 meV^34,35, where ℏω denotes neutron energy transfer, with the scattering being confined to the (hk0) plane. Additionally, an excitation observed a ℏω=1 meV was suggested to be a resonant excitation which originates from a bound state within the particle-hole continuum gap^35,36. This was taken as evidence that the spin fluctuations were primarily antiferromagnetic in nature, as one would expect ferromagnetic fluctuations at integer Q vectors, where Q denotes momentum transfer.

Here, we report INS results on UTe₂ in the (hk0) scattering plane with an applied field along the (hat{c}) axis. Our results resolve the low energy excitations across multiple Brillouin zones in the (hk0) scattering plane at energies below 6 meV, accounting for nearly all the expected spectral weight from the magnetic moment of U³⁺. The measured scattering intensity is independent to an applied magnetic field of μ₀H ≤11 T, suggesting that the magnetic response is dominated by interband excitations excited across the df Kondo hybridization gap, rather than ferromagnetic or antiferromagnetic exchange between local spins.

Inelastic Neutron Scattering Results

The INS measurements were performed on a mosaic of coaligned single crystalline samples of UTe₂ ranging in T_c ∈ {1.8, 2.0} K, with bulk properties summarized in Fig. 1a, b. Scattering integrated over the elastic window (ℏω ∈ {−0.15, 0.15} meV) shows no evidence of magnetic or field-dependent scattering. A mosaic spread of full width at half maximum (FWHM) 5.0(1)^o in the (hk0) plane is determined from the nuclear Bragg peaks. A sample environment background is determined using a momentum averaged scattering from the regions in Q, ω space where magnetic scattering from the sample is taken to be absent, which is described in detail in the Supplementary Information (SI).

a Specific heat from two samples grown by chemical vapor transport, showing a single sharp transition indicative of good sample homogeneity. b Transport measurements showing the superconducting transition in four samples. The solid black line represents the same sample in both C(T) and R(T). c, d Two representative copper plates with coaligned samples used in the mosaic.

Full size image

Two representative constant energy slices of the magnetic scattering from UTe₂ are shown in Fig. 2. The scattering in all cases has been symmetrized by the two-fold C_2m lattice symmetry to enhance statistics. The lowest energy contribution to the scattering is shown in Fig. 2a, b. As in previous studies, the scattering is constrained to the Brillouin Zone (BZ) edge, emanating from the corners (T=((frac{1}{2},frac{1}{2}),0)-like points) and edges (Y=(0, 0.6, 0)-type points). Scattering at higher Brillioun zones is also observed.

The energy integration windows are ℏω ∈ [1.7, 2.3] meV in (a, b), and ℏω ∈ [3.2, 3.8] meV in (c, d). All observable scattering at the lowest energies originates from the (Y,(0frac{1}{2}0)) and T ((frac{1}{2}frac{1}{2})0) points, as shown in (a, b), with the lowest energy scattering emerging from the points denoted as Y (edge) and T (corner) type points. The magnetic field in (b, d) is applied along the (hat{c})-axis.

Full size image

Constant energy slices in energy window of maximum scattering intensity ℏω ∈ {3.2, 3.8} meV are shown for the 0 T measurement in Fig. 2c and 11 T in Fig. 2d. Here the scattering remains on the BZ edges, but the excitations have dispersed such that they have smeared in the (h00) direction. A minimum in intensity has emerged along the (0k0) line and zero intensity is observed at all Γ points. As is made more quantitatively clear in Fig. 3, there is no qualitative difference in scattering intensity for the 0 T and 11 T settings apart from a small deviation shown in Fig. 3(a). This also true for the elastic scattering, as shown in the SI.

a, b Constant Q cuts along the energy dimension of UTe₂ scattering integrated in a circular window of diameter ΔQ = 0.2 Å⁻¹, centered at the Y₁ (0,0.6,0) and T₂ (0.4,1.4,0) points. Black filled circles represent 0 T measurements, red unfilled square makers 11 T. c Intensity integrated in the low-energy portion of the measurement of 1 meV to 2 meV. The cut is along the BZ edge, with an integration with in Q of 0.3 Å⁻¹. d Scattering intensity in the same energy window along the (0k0) direction. Both cuts are fit simultaneously using a form described in the text, with the result shown in blue. e Goodness of fit χ² value for fit of unsymmetrized constant energy scattering integrated from ℏω ∈ {1.7, 2.3} meV to Eq. (3) as a function of the magnetic moment direction, where θ_a is the tilt away from the a axis towards the b axis. All error bars represent one standard deviation and cuts are from unsymmetrized data.

Full size image

In Fig. 3, representative intensity cuts along Q and energy transfer are shown. All data presented in cuts is not symmetrized. Figure 3(a) and (b) plot constant Q cuts along energy at the high symmetry points Y₁ = (0, 0.6, 0) in (a) and T₂ = (0.4, 1.4, 0) in (b). Both cuts are integrated in the Q dimension in a circle of radius δQ = 0.3 Å⁻¹ (h ± 0.2 r.l.u., k ± 0.29 r.l.u.). The energy dependence of the intensity in both plots shows a monotonic increase in intensity from ℏω = 0, with a peak of intensity at 3.2(1) meV at Y₁ and 4.0(1) meV T₂ meV respectively. There is no statistically significant field-dependence to the scattering at the T₂ point. At energies below ℏω = 3meV, the zero field scattering is slightly more intense at the Y₁ point, as shown in Fig. 3a. The integrated intensity within ℏω ∈ {0.3, 3.0} is 2.0(3) b/U for μ₀H = 0 T, and 1.5(3) b/U for μ₀H = 11 T. We are unable to resolve the proposed spin-resonance excitation reported in refs. ^34,35.

Figure 3c, d show constant energy cuts integrated in ℏω ∈ {1, 2} meV. This integration window accounts for most of the scattering before dispersion becomes significant. Consistent with previous studies^{33,34,35,36,37}, we find broad peaks in intensity centered at the Brillouin zone edges. Again, there is no difference between the 0 T and 11 T scattering intensities. In (c) the cut is along the (h, 1.5, 0) direction, i.e. along the BZ edge, and in (d) the cut is along the (0, k, 0) direction. In both cases the scattering was integrated in a perpendicular Q range of 0.3 Å⁻¹. The cuts are simultaneously fit to a model of nearest neighbor correlated spins described later in the text.

Finally, Fig. 4a, b present the full dispersion of the magnetic excitations in UTe₂ along the zone boundary, with Fig. 4c including the Γ₀ = (000), Γ₁ = (010), and Γ₂ = (020) points. The Q path in Fig. 4a is primarily along the (h00) direction and is depicted by the blue path in the schematic of the Brillouin zone in (c), and (b) is along the (0k0) direction and depicted by the red path in (f). The third path along (0k0) is shown in (c), revealing no dispersion at the Y-type points and zero intensity at the zone center. The scattering is broader in energy than the expected instrumental resolution (δE_FWHM = 0.11 meV for ℏω=0 meV), and at every point in Q the energy dependent intensity I(ℏω) may be fit to the form of a simple Lorentzian peak of FWHM width ΔE_FWHM=3.0(2) meV. For all values of Q with significant enough intensity to perform this fit, the extracted dispersion of magnetic excitations along the representative paths in Q is presented in Fig. 4d, e. Along both directions, the peak in scattering intensity has dispersion with clear periodicity. Due to the magnetic form factor of U^3+/4+ and a net spin polarization factor in the neutron scattering cross-section, these excitations become faint at the highest accessible values of Q. The dispersion is steeper along the (0k0) direction, as shown both in Fig. 4b, e, and in all cases peaks at the T points.

a, b Magnetic excitation in UTe₂ using paths along the BZ edge as depicted in (f). The path in Q for panel (a) is primarily along the (h,1.5,0) direction and is in blue, and the path for (b) is primarily along the (0.5,k,0) direction and is shown in red. At every point in Q, a Lorentzian form is used to fit the intensity maxima as a function of energy transfer, giving effective dispersions shown in (d, e). (d) corresponds to subplot (a), and (e) corresponds to (b). Values of Q with insufficient statistics to find a stable intensity maxima have been omitted from these plots. Cyan triangles in (d) are the dispersion observed in ref. ³³. c Magnetic excitation dispersion along the (0, k, 0) direction through the zone center, as depicted by the cyan path in (f). All error bars represent one standard deviation and all slices use symmetrization.

Full size image

The second difference between the observed scattering and previous reports is that the resonant excitation reported in refs. ^35,36 is not seen in our measurement. While this resonance was initially reported at the Y₁ = (0, 0.6, 0) point, independent measurements also observed the excitation at the Y₂ = (0, 1.4, 0) point³⁶. Here, our results do not show evidence of a peak in scattering at ℏω ≈ 1 meV at either Y₁ or Y₂. A more detailed examination is in the SI, but this does not preclude the existance of the resonant excitation UTe₂. There are many different details of the experiment that may be responsible for this, such as mosaic distribution, detector coverage, or sample-dependence, and the resonance is not the focus of this work.

Overall, the scattering presented is compatible with previous INS studies in the superconducting phase. There are two significant differences. First, these data allow for the full extraction of the dispersion of magnetic excitations in the (hk0) scattering plane. An important result that we extract from this is the quantity of the fluctuating moment for both the 0 T and 11 T configurations. Because the scattering has been normalized, the well-known total moment sum rule^38,39 of the form

The dynamical structure factor S(Q, ω) is related to the scattering intensity by

Here, the scattering has been normalized by the kinematic ratio of final and initial momenta (frac{{k}_{f}}{{k}_{i}}), the Lande g-factor g is taken to be 2, ∣F(Q)∣² is the modulus of the magnetic form factor associated with scattering from U³⁺ or U⁴⁺, which are indistinguishable in this measurement, r₀ = 5.391 ⋅ 10⁻¹³ cm is the characteristic magnetic neutron scattering length, and the summation is over the Cartesian x, y, z spin directions.

As the experiment cannot access to the (00l) scattering direction, to perform the integration over the full Brillouin zone we refer to the results of ref. ³⁷ which suggest no dispersion along the (00l) direction and intensity with a Gaussian distribution centered at l = 0, with a width of σ_FWHM ≈ 0.5 rlu. The scattering was integrated in energy over the full range of observed excitations, ℏω ∈ {0.3, 6} meV, and in two Brillouin zones h ∈ {0, 1} and k ∈ {0.5, 2.5}. From this, we find an overall fluctuating moment of μ_eff =1.7(5) μ_B, which may be compared to the itinerant moment of approximately 1.1 μ_B as measured by high field magnetization⁴⁰. It may also be compared with the moment found from the paramagnetic susceptibility of 3.6μ_B⁴¹ or free-ion U³⁺ of 3.62 μ_B. One may also compare the moment inferred from the total sum rule to the itinerant moment, which high-field magnetization finds to be near 1.1 μ_B⁴⁰. By comparing the zero field and 11 T integrated Bragg peak intensities at the (110), (1(bar{1})0), and (020) Bragg peaks we find an upper limit of the field-induced ordered moment of μ = 0.0(2) μ_B. This is slightly at odds with what one would expect from magnetization, which for fields along the b axis finds an induced moment on the order of μ = 0.15 μ_B¹ at 11 T. A more detailed diffraction study using a large single crystal rather than a mosaic of coaligned samples would be more appropriate to resolve this moment.

The lowest energy part of the scattering emerges from the Y and T points. The Q-dependent intensity of the scattering can be directly calculated by assuming a form of S(Q, ω), where all spins in the unit cell have the same preferential spin orientation. The scattering takes the form

The sum is over all of the Y and T points, where the intensity is assumed to have a Gaussian line shape of amplitude A_T for T-points and A_Y for Y-points. The peaks are of width σ_h and σ_k in the h and k directions, respectively. The prefactor ({mathcal{C}}) is defined by

Here, the preferential spin orientation (hat{M}) is constrained to be in the (hk0) plane, with the tilt away from the a axis towards the b axis being denoted as θ_a. No significant difference could be found between constraining the moment to the (hk0) plane and an entirely free orientation, apart from an upper limit of 35^o from the (hat{a})-axis towards the (hat{c})-axis. Fitting the constant energy slices integrated within ℏ ∈ {1.7, 2.3} meV finds best fit values of θ_a=θ_a=17(3) degrees as depicted in Fig. 3e. The minimum at θ_a=17(3) degrees is shallow, but the scattering is clearly incompatible with a spin orientation beyond 30 degrees. This result demonstrates that the spins contributing to the scattering prefer to orient along the crystalline (hat{a})-axis within the ab plane, which is consistent with previous neutron studies³⁴ and bulk magnetometry indicating that (hat{a}) is the magnetic easy axis¹.

The character of the scattering in Fig. 2a gives important information about the electronic band structure near the Fermi level ϵ_F. As shown by the fit in 3c, d, the U^3+/4+ form factor is sufficient to describe the intensity modulation of the scattering with increasing ∣Q∣, indicating that the spin density of electrons in the conducting band originates from the 5f orbital. This is consistent with recent RIXS reports⁴² suggesting a 5f² configuration. Secondly, we note that the scattering is found at Y = (0, 0.6, 0) and (T=(frac{1}{2}frac{1}{2}0)) type points, but is absent at the (M=(frac{1}{2}00)) type points. This cannot be explained exclusively through the effect of spin polarization, as there are many symmetry-equivalent M-type points available in the (hk0) scattering plane where no scattering is observed.

The scattering observed in this experiment confirms results from previous works^{33,34,35,36,37}, and has direct physical consequences. The first is that at temperatures below T_K ≈ 40 K, a dispersive magnetic excitation mode emerges which is peaked in intensity at 3.5(2) meV and is broad in energy (ΔE_FWHM ≈ 1.5 meV). While the scattering intensity may persist to zero energy transfer, at the lowest accessible energies the intensity becomes vanishingly small and the maxima of all excitations bands are fully gapped. There is no evidence for long-range magnetic order, but the fluctuating moment has preferential spin orientation along the (hat{a}) axis. The excitations themselves are constrained to the BZ edge, with no significant scattering at the Γ point at all measured energy transfers. Finally, the momentum dependence of the scattering is described well by ionic U^3+/4+ and a total moment sum rule suggests that the scattering captures the majority of the relevant magnetic excitations. Though these results are informative, the interpretation of the origin of the magnetic scattering is less straightforward.

Previous INS studies of UTe₂ have characterized these excitations at non-integer Q as antiferromagnetic spin fluctuations, suggesting a picture with dominant antiferromagnetic interactions between nearest neighbor U ions^35,37. Subsequent analysis of neutron scattering spectra suggested the presence of both ferro and antiferromagnetic interactions in a spin-ladder model, where spin exchange interactions are ferromagnetic between the U³⁺ chains along the (hat{a})-axis, ferromagnetic within the U-dimers along the (hat{c})-axis, and antiferromagnetic between ladders along the (hat{b})-axis³⁴. Additionally, a pressure-induced antiferromagnetic ordered phase in UTe₂ was found at an experimentally modest hydrostatic pressure of about P_c=1.5 GPa^43,44. These studies report an incommensurate AFM ordering wavevector of k = (0.07, 0.33, 1)⁴⁴, and the close proximity of the ground state to an AFM ordered phase was taken as evidence for strong AFM magnetic interactions. This contradicts initial suggestions that spin-triplet superconductivity in UTe₂ is promoted by proximity to a ferromagnetic critical point. While spin-triplet superconductivity has been suggested to be possible in the case of AFM interactions⁴⁵, the precise nature of the spin interactions is still an important open question.

A magnon-like picture where the dispersion of the magnetic excitations originates from AFM interactions between spins raises a number of questions. The first is how magnons can emerge from a ground state that is paramagnetic. In this scenario, the excitations would be classified as paramagnons, where spins are proximate to local magnetic order due to strong spin correlations. Such a scenario is compatible with the broad signal observed in UTe₂. However, our data show that the magnetic response has no dependence on the applied magnetic field up to a relatively large field scale of μ₀H=11 T, contrary to what one would expect in the case of spins coupled with a Heisenberg interaction on the scale of 4 meV. Additionally, in the generic Heisenberg picture using linear spin wave theory the magnon dispersion minima occur at the magnetic ordering wave vector k_m. Here, and in previous works^33,35,37, these minima are observed at the in equivalent Y and T points, requiring the further complexity of two different species of magnetic order or a multi-k structure. Because such a scenario quickly becomes quite complicated, we consider an alternative picture.

Neutron scattering from localized U-f spins in heavy fermion compounds has previously been treated in the context of an extended Anderson lattice model^46,47. In the case of UTe₂, this would mean that the Kondo effect hybridizes the U-6d, U-5f, and Te-5p electronic bands near the chemical potential. In the extended Anderson lattice picture, the newly formed hybridized bands are constrained by symmetry to have saddle points with a significant density of states at high symmetry points at the Brillouin zone edge. This treatment provides an alternative description of the scattering, which is that of interband electronic transitions between regions with high density of state between the hybridized bands, which are at the zone edge and zone boundary⁴⁸ and naturally account for the field-independence of the excitations.

To explore this scenario, we compare to high resolution DFT+DMFT calculations of the band structure of UTe₂. These calculations were performed using the Rutgers code implemented in the WIEN2k package using the experimentally determined crystal structure⁴⁹, at temperatures of T=116 K, T=232 K, and T=580 K. A computational limitation of DFT+DMFT is the lower temperature limit, and thus the energy resolution of the band structure on the ≈ 10 meV scale. Thus, the detailed band structure near the chemical potential ϵ_f is not resolved in these calculations. Figure 5b, d plot the calculated quasiparticle band structure, which closely resembles our previously reported results in ref. ⁵⁰. We also note that these calculations are similar to those presented in ref. ³⁷, with the primary difference being that there is no infinite Hubbard U approximation. These results give insight into the hybridization associated with the f − band correlations, and while the precise energies within 10 meV of ϵ_f are unreliable due to the temperature constraint, the development of anticrossings of bands is evident. The anticrossings are denoted using arrows in Fig. 5c, and result in regions of high f-density of states at high symmetry points, between which interband excitations may occurr.

The high temperature calculation is shown in (a) and the low temperature in (b). The calculated low temperature electronic disperison throughout the full Brillouin zone is shown in (c), with the low energy scale most relevant to the quasiparticle excitations shown in (d). The green arrows in (c) highlight the anticrossings in the quasiparticle band structure that originate from f-band correlations.

Full size image

The calculated high temperature Fermi surface is shown in Fig. 5a, which finds two elliptical pockets similar to previous experimental and theoretical reports^51,52,53. At low temperature, Fig. 5(b), the effect of hybridization may be observed with the shape of the pockets changing dramatically and a new pocket emerging around the Z-point. This can be compared directly to photoemission experiments performed at T = 20 K, which report a heavy electron pocket at the Z point⁵³, and quantum oscillation measurements which suggest both 2D and 3D Fermi surface character, with conflicting conclusions as to the presence of such a pocket^54,55. In both plots, the blue/red surfaces originate from Te-5p band, and the blue/yellow surfaces originate from the U-6d band.

To examine the redistribution of f − weight at low temperature, we examine constant energy k_z=0 slices in Fig. 6. Figure. 6a, b show the k_z=0 Fermi surface for T=116 K and T=580 K respectively. We take the temperature difference of the spectral function A(Q, ω) to be a measure of the U-5f band weight, which accumulates in pockets centered about the X point at E=E_f, as shown as a function of energy transfer in Fig. 6c–f. The f-weight disperses towards the zone corner at higher energy transfers, which is suggestive of the origin of the dispersion of the excitations observed in INS. This dispersion, which originates from hybridization, is also anisotropic in the (hk0) plane, as most clearly demonstrated at the chemical potential in Fig. 6d.

Directly calculated Fermi surfaces plotted for T=116 K (a) and T=580 K (b) show the effect of hybridization and the accumulation of f-weight at low temperature. The difference between the two temperatures is shown at E − E_f = −10 meV (c), 0 meV (d), 10 meV (e), and 20 meV (f), which is taken as a measure of the distribution of f-electron character.

Full size image

Using the distribution of the f-electron character presented in Fig. 6d, it is tempting to assign potential nesting vectors between regions within high spectral weight. This would lend support to the idea of an emergent charge density wave which has been observed in scanning tunneling microscopy (STM) studies, but has not been observed in bulk scattering measurements^{56,57,58,59,60}. As no such wavevector is observed in the elastic component of the magnetic scattering, such a phase would need to be a surface state rather than in the bulk.

Due to the energy resolution of the calculation, it is not possible to determine the precise scattering pathways that would lead to the measured S(Q, ℏω) in INS. However, the calculations do show that the hybridized band structure features parabolic-like bands with extrema at the zone boundary, as suggested in the Anderson picture. Two essential conditions are required for interband scattering, the first is the presence of a high density of states near the chemical potential at the high symmetry points at which scattering is observed, and the second being anisotropy in the band hybridization resulting in the lack of scattering at the X point, both of which are captured in the DFT+DMFT calculations.

Our work resolves the magnetic excitations in UTe₂ that peak at ℏω=3.5(2) meV, which we propose originate from electronic excitations near or across the Kondo hybridization gap. This conclusion is also supported by our previous work³³, in which the temperature dependence of the magnetic scatting intensity is linked to the Kondo coherence tmperature. Here, the excitations capture a total fluctuating magnetic moment of μ_eff =1.7(5) μ_B, which is approximately compatible with what one would expect for free-ion U³⁺. Our analysis shows that their origin is likely not from magnetic interactions, but no conclusion is made regarding the proximity of UTe₂ to any particular ordered state, or the nature of its magnetic interactions which have been reported as both ferromagnetic¹, and antiferromagnetic³⁷.

This is consistent with previous DFT+U studies, which suggest that ferromagnetic inter dimer interactions between U-ions may stabilize spin-triplet superconductivity^61,62,63. The scattering captures the full expected magnetic moment in the system, and we expect that this scenario could be supported by further DMFT and calculations of ({chi }^{{prime} }({boldsymbol{Q}},hslash omega )) similarly to the studies of refs. ^29,30. As we are unable to perform such calculations of ({chi }^{{prime} }({boldsymbol{Q}},hslash omega )) in this work, we cannot directly link the observed scattering to the calculated electronic band structure, and leave open the possibility of alternative scattering mechanisms. The DMFT results can be verified directly through photoemission high resolution experiments at low temperatures, which would be able to resolve the elliptical pockets predicted in our calculations. Additionally, we hope that this work will motivate the further development of DFT+DMFT methods to enable calculations of ({chi }^{{prime} }({boldsymbol{Q}},hslash omega )) in the meV scale energy regime, which is relevant not only in the case of UTe₂ but many other heavy fermion systems with Kondo hybridization.

The INS experiment was performed on a coaligned mosaic of 82 crystals of UTe₂ with a total mass of 1.1 g. Most crystals were grown by chemical vapor transport⁶⁴, with eight of total mass ≈ 80mg being grown by a salt flux method^2,3. All batches of crystals were screened by either transport measurements or susceptibility, with a range of critical temperatures of T_c ∈ {1.8, 2.0} K indicating high sample quality. The crystals were mounted on a oxygen-free copper sample holder using CYTOP (AGC Chemicals Company) and aligned using a x-ray Laue diffractometer.

The INS measurement was performed using the CAMEA spectrometer⁶⁵ at the Paul Scherrer Institut to measure the inelastic scattering in an energy transfer range of ℏω ∈ { − 1, 6} meV. The sample environment was the MB11 μ₀H=11 T vertical cryomagnet with a dilution refrigerator insert with a base temperature of T=45 mK throughout the experiment based on the sample thermometer. All measurements were performed in the (hk0) scattering plane, with the magnetic field oriented along the (hat{c}) direction. For the zero-field configuration, the total counting time was 61 hrs, and for the 11 T configuration 64 hrs. The scattering was normalized to absolute units of (b/meV/sr/U) using a vanadium standard of known mass. All analysis was performed using the MJOLNIR software^66,67.