Hyperfine van der Waals Interaction

• Hyperfine van der Waals interaction is a class of forces where hyperfine structure determines whether intermolecular coupling occurs in first or second order, affecting observable shifts.
• Experimental studies reveal its impact through scalar nuclear spin-spin coupling in liquids and hyperfine-resolved spectral shifts in hydrogen atoms, highlighting distinct coupling pathways.
• In ultracold open-shell molecules, hyperfine selection rules can suppress resonant dipole exchange, leading to repulsive long-range potentials that minimize inelastic losses.

Searching arXiv for recent and foundational papers on hyperfine van der Waals interactions. Hyperfine van der Waals interaction denotes a class of intermolecular effects in which hyperfine structure, or a hyperfine-mediated coupling mechanism, qualitatively determines a van der Waals observable. In the literature represented here, the term is used in three closely related but technically distinct senses: a scalar nuclear spin-spin coupling transmitted across transient van der Waals contact in liquids; hyperfine-resolved van der Waals shifts in quasi-degenerate excited hydrogen; and a second-order $R^{-6}$ interaction between open-shell polar molecules that emerges when resonant dipole-dipole exchange is forbidden by hyperfine selection rules (Ledbetter et al., 2011, Jentschura et al., 2017, Adhikari et al., 2017, Walraven et al., 16 Jul 2025). Taken together, these works suggest that hyperfine van der Waals physics is best understood not as a single universal Hamiltonian term, but as a family of hyperfine-conditioned long-range or contact-dominated interactions.

1. Conceptual scope and distinguishing features

The common theme is that hyperfine structure does not merely add a perturbative correction to an otherwise fixed van der Waals potential. Instead, it determines whether an interaction is present in first order or only in second order, which virtual states dominate the effective coupling, and which observable carries the signature of the interaction. In one setting, the relevant coupling is the ordinary NMR scalar interaction,

$H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$

whose microscopic origin is a second-order Fermi-contact hyperfine mechanism mediated by electrons. In another, the relevant operator is the dipole-dipole interaction, but hyperfine resolution determines whether degenerate exchange channels survive, producing $1/R^3$ behavior, or are suppressed so that the leading effect is $1/R^6$. In hydrogen, the decisive small denominators are the Lamb shift and hyperfine splittings; in open-shell molecules, they are hyperfine splittings within a rotational manifold; in hyperpolarized xenon solutions, the key object is an ensemble-averaged scalar coupling over transient encounters.

These distinctions are essential because the phrase “van der Waals” can otherwise be misleading. The repulsive interaction predicted for open-shell molecules is explicitly not ordinary electronic van der Waals attraction, and the scalar coupling observed in xenon-pentane is not a resolved line splitting of a bound complex. Likewise, the hydrogen calculations show that quasi-degenerate manifolds invalidate a naive nondegenerate dispersion treatment. A recurrent misconception addressed in these works is therefore that van der Waals physics is always isotropic, always attractive, or always describable by a single species-independent $C_6$. The cited studies show instead that hyperfine structure can select the sign, magnitude, scaling law, and even the appropriate experimental observable.

2. Scalar nuclear spin-spin coupling across van der Waals contact

The most direct experimental realization of a hyperfine van der Waals interaction is the observation of scalar coupling between unbound nuclei in a solution of hyperpolarized $^{129}\mathrm{Xe}$ and pentane. The central question was whether the scalar nuclear spin-spin coupling familiar from conventional NMR spectroscopy can survive in systems that are not covalently bound, but instead interact only transiently through van der Waals contact. The reported result was the first observation of such couplings in van der Waals molecules, establishing direct experimental evidence for a hyperfine van der Waals interaction of the scalar, electron-mediated type (Ledbetter et al., 2011).

The mechanism is the same second-order hyperfine process responsible for through-bond $J$-couplings. Each nucleus couples to the electronic spin density at its location, and second-order mixing through the electronic structure generates an effective bilinear interaction between the two nuclear spins. In the rapidly exchanging liquid, one does not resolve line splittings from individual Xe$\cdots$H encounters. Instead, if one spin species is strongly hyperpolarized, the averaged scalar interaction produces a mean-field frequency shift. For proton spins,

$$H_Z = -\hbar \gamma_1\, \mathbf I_1 \cdot \mathbf B_0,$$

and under rapid exchange the scalar term becomes

$$H_J = h\langle J\rangle I_{1z}\langle I_{2z}\rangle,$$

leading to

$H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$0

The shift can also be written in contact-shift form,

$H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$1

with

$H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$2

The experiment used hyperpolarized liquid $H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$3 and pentane in a spherical cell half-filled with pentane and then condensed xenon, giving an approximate xenon number density

$H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$4

Typical xenon polarizations were $H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$5–$H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$6. NMR was performed in a static field of about $H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$7 in a magnetically shielded environment, and high-$H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$8 SQUID magnetometers provided a magnetometric sensitivity of about

$H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$9

After a proton $1/R^3$0 pulse, the proton and xenon free induction decays were observed near $1/R^3$1 and $1/R^3$2, respectively. Proton $1/R^3$3 was typically about $1/R^3$4, xenon $1/R^3$5 about $1/R^3$6, and about $1/R^3$7 proton transients could be acquired from a single xenon batch.

A major issue was separation of the scalar shift from ordinary dipolar fields and SQUID-induced magnetic artifacts. The data were taken with $1/R^3$8 along three orthogonal directions $1/R^3$9, $1/R^6$0, and $1/R^6$1, and the results were combined so that orientation-dependent dipolar demagnetization fields were compensated. After correction for SQUID-induced shifts, the $1/R^6$2 frequency shift was consistent with zero, while the proton shift retained a nonzero component. Averaging over field directions and subtracting the SQUID contribution gave a net proton shift of $1/R^6$3 per unit xenon magnetization, corresponding to

$1/R^6$4

and

$1/R^6$5

The microscopic interpretation combined relativistic DFT and liquid-state averaging. The coupling was expressed as

$1/R^6$6

or, equivalently,

$1/R^6$7

Relativistic DFT with ZORA implemented in ADF at the scalar ZORA BP86/TZ2P level yielded $1/R^6$8 coupling constants on a cubic grid around pentane. The calculated $1/R^6$9 values were mostly negative and ranged between about $C_6$0 and $C_6$1, with a few positive values, and were fit over $C_6$2–$C_6$3 to

$C_6$4

Molecular dynamics with OPLS-AA for pentane and a Lennard-Jones model for xenon gave the RDF, and the resulting average,

$C_6$5

agreed well with the measured $C_6$6.

The significance of this result is specific. The measured quantity is an ensemble average over rapidly exchanging configurations, not a resolved pairwise coupling for a single van der Waals complex. The experiment therefore established that electron-mediated scalar coupling can persist across transient noncovalent contact and be read out as a magnetization-dependent NMR line shift, rather than as a conventional splitting pattern.

3. Hyperfine-resolved van der Waals structure in the hydrogen $C_6$7-$C_6$8 system

In the excited hydrogen problem, hyperfine van der Waals interaction appears in a different form. The system consists of two hydrogen atoms restricted to the $C_6$9, $^{129}\mathrm{Xe}$0 manifold, namely $^{129}\mathrm{Xe}$1 and $^{129}\mathrm{Xe}$2, each split into $^{129}\mathrm{Xe}$3 and $^{129}\mathrm{Xe}$4 hyperfine levels. There are therefore $^{129}\mathrm{Xe}$5 one-atom states and $^{129}\mathrm{Xe}$6 two-atom product states before symmetry reduction. The total Hamiltonian is

$^{129}\mathrm{Xe}$7

with the nonretarded dipole-dipole interaction

$^{129}\mathrm{Xe}$8

The key scales are

$^{129}\mathrm{Xe}$9

with the hierarchy $J$0 valid for $J$1 (Jentschura et al., 2017).

The paper’s main conceptual point is that hyperfine resolution determines whether a given long-range interaction is first-order or second-order. The conserved quantum number is

$J$2

which decomposes the $J$3-dimensional space into $J$4, $J$5, and $J$6 manifolds of dimensions $J$7, $J$8, and $J$9, respectively. Each further splits into two irreducible uncoupled subspaces, one containing $\cdots$0-$\cdots$1 and $\cdots$2-$\cdots$3 combinations and one containing $\cdots$4-$\cdots$5 and $\cdots$6-$\cdots$7 combinations. The most complicated irreducible block is $\cdots$8-dimensional.

For states involving one $\cdots$9 atom and one $$H_Z = -\hbar \gamma_1\, \mathbf I_1 \cdot \mathbf B_0,$$0 atom, exact or near-exact degeneracy under interchange $$H_Z = -\hbar \gamma_1\, \mathbf I_1 \cdot \mathbf B_0,$$1 allows direct coupling by $$H_Z = -\hbar \gamma_1\, \mathbf I_1 \cdot \mathbf B_0,$$2, giving first-order shifts linear in

$$H_Z = -\hbar \gamma_1\, \mathbf I_1 \cdot \mathbf B_0,$$3

For example, certain $$H_Z = -\hbar \gamma_1\, \mathbf I_1 \cdot \mathbf B_0,$$4 states have energies

$$H_Z = -\hbar \gamma_1\, \mathbf I_1 \cdot \mathbf B_0,$$5

By contrast, when both atoms are metastable $$H_Z = -\hbar \gamma_1\, \mathbf I_1 \cdot \mathbf B_0,$$6 states, the hyperfine-resolved states are not directly coupled to energetically degenerate partners, so the leading interaction is second order,

$$H_Z = -\hbar \gamma_1\, \mathbf I_1 \cdot \mathbf B_0,$$7

The characteristic scale is

$$H_Z = -\hbar \gamma_1\, \mathbf I_1 \cdot \mathbf B_0,$$8

The hyperfine-resolved treatment makes the denominators state dependent. For $$H_Z = -\hbar \gamma_1\, \mathbf I_1 \cdot \mathbf B_0,$$9-$$H_J = h\langle J\rangle I_{1z}\langle I_{2z}\rangle,$$0 states, virtual transitions involve combinations such as $$H_J = h\langle J\rangle I_{1z}\langle I_{2z}\rangle,$$1, $$H_J = h\langle J\rangle I_{1z}\langle I_{2z}\rangle,$$2, and $$H_J = h\langle J\rangle I_{1z}\langle I_{2z}\rangle,$$3. In the $$H_J = h\langle J\rangle I_{1z}\langle I_{2z}\rangle,$$4 manifold, the effective Hamiltonian for a degenerate $$H_J = h\langle J\rangle I_{1z}\langle I_{2z}\rangle,$$5-$$H_J = h\langle J\rangle I_{1z}\langle I_{2z}\rangle,$$6 pair yields eigenvalue shifts

$$H_J = h\langle J\rangle I_{1z}\langle I_{2z}\rangle,$$7

and analogous expressions with $$H_J = h\langle J\rangle I_{1z}\langle I_{2z}\rangle,$$8 for the companion pair. In the $$H_J = h\langle J\rangle I_{1z}\langle I_{2z}\rangle,$$9 sector, even some exactly degenerate states remain unshifted in first order because the direct van der Waals matrix elements vanish; their splittings arise only through second-order effective interactions.

The physical consequence is that hyperfine structure reorganizes the long-range interaction landscape. A hyperfine-unresolved treatment would miss which channels possess resonant $H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$00 exchange and which remain purely dispersive at $H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$01. For spectroscopy, the practically important statement is explicit: all states with both atoms in metastable $H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$02 levels are shifted only in second order, while states with one atom in a $H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$03 state can acquire large first-order hyperfine-resolved van der Waals splittings.

4. Dirac-$H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$04 hyperfine modification of the hydrogen $H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$05-$H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$06 interaction

The $H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$07-$H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$08 problem isolates another aspect of hyperfine van der Waals physics: a hyperfine-induced modification of an already existing dispersion interaction. The reference configurations,

$H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$09

are degenerate under exchange. The leading interaction comes in second-order perturbation theory through virtual $H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$10-state excitations, governed by

$H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$11

The crucial feature is that the excited $H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$12 atom has quasi-degenerate $H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$13 and $H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$14 virtual states, separated only by the Lamb shift $H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$15 and fine-structure interval $H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$16. If $H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$17 and $H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$18 are set to zero too early, an important contribution to $H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$19 is lost; the paper resolves a discrepancy in the literature precisely on this point (Adhikari et al., 2017).

In the nonretarded van der Waals range,

$H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$20

the interaction has the symmetry-dependent form

$H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$21

The direct and exchange coefficients are

$H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$22

so

$H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$23

The large quasi-degenerate contribution

$H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$24

is the missing piece in earlier calculations that obtained only about $H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$25.

The hyperfine modification is modeled by a local Dirac-$H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$26 perturbation. For $H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$27-states, the relevant term in the hyperfine Hamiltonian is the Fermi contact interaction,

$H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$28

To treat it generically, the paper introduces

$H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$29

with

$H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$30

and then substitutes

$H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$31

This gives the same physics in two equivalent representations: a correction to the intermolecular potential, and the shift of the $H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$32 hyperfine frequency caused by nearby $H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$33 atoms.

In the van der Waals regime, the correction renormalizes the $H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$34 coefficient: $H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$35 so

$H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$36

In the intermediate regime,

$H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$37

the hyperfine correction changes character. The energy-type contribution becomes more singular than the unperturbed leading term and scales as $H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$38: $H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$39 The paper notes that logarithmic terms generated by individual retardation contributions cancel in the final result. In the extreme Lamb-shift range,

$H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$40

the residual interaction is negligible, below $H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$41.

The experimental connection is explicit. The $H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$42 hyperfine splitting was measured as

$H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$43

and the theory gives shifts of the $H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$44 hyperfine interval due to a neighboring $H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$45 atom of

$H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$46

$H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$47

and

$H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$48

The paper therefore frames the hyperfine van der Waals effect both as a modification of the interatomic potential and as a measurable hyperfine-frequency shift in precision spectroscopy.

5. Hyperfine van der Waals repulsion in open-shell polar molecules

A newer use of the term refers to a distinct long-range interaction mechanism between ultracold open-shell polar molecules in rotational states differing by one quantum. The ordinary expectation for a pair such as $H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$49 is resonant dipole-dipole exchange, because $H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$50 is degenerate with $H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$51, giving a first-order interaction

$H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$52

For open-shell molecules, however, the rotational levels are split into hyperfine sublevels by rotation, electron spin, and nuclear spin couplings. For certain hyperfine choices, dipole selection rules forbid the exchange process

$H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$53

The paper states the relevant dipole selection rule as

$H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$54

When the exchange matrix element vanishes, the first-order resonant $H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$55 interaction is removed, and the leading term becomes second order in the dipole-dipole operator,

$H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$56

Because the denominators are only hyperfine splittings, typically tens to hundreds of MHz, the resulting $H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$57 interaction can be unusually strong and can be either attractive or repulsive (Walraven et al., 16 Jul 2025).

The paper focuses on laser-coolable $H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$58 molecules, especially CaF, while also showing the effect for MgF, SrF, BaF, and YO. In CaF, the $H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$59 manifold has $H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$60, and $H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$61 has $H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$62, with splittings on the order of $H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$63–$H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$64 MHz in the low rotational states shown. The dimer Hamiltonian contains radial kinetic energy, end-over-end rotation, the two single-molecule Hamiltonians, and the intermolecular dipole-dipole operator

$H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$65

equivalently written in spherical-tensor form.

CaF provides the clearest example of hyperfine van der Waals repulsion. In the $H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$66 manifold, the pair states $H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$67 and $H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$68 do not undergo resonant dipole-dipole exchange. For $H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$69, the transition

$H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$70

is forbidden because $H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$71, so the exchange process

$H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$72

does not occur. The dominant second-order coupling is then to a lower-lying $H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$73 pair state, so the induced $H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$74 interaction is repulsive. The lowest $H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$75-wave adiabat for the $H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$76 channel develops a long-range repulsive tail, then becomes strongly attractive at shorter distance, producing a barrier of order mK in height. Collisions below the barrier are reflected, and loss can occur only by tunneling to short range or by inelastic transitions.

The scattering calculations are full coupled-channels calculations with absorbing boundary conditions at short range. Most hyperfine combinations in CaF are highly lossy because resonant dipole-dipole interactions dominate. The attractive $H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$77 channel approaches universal loss at low temperature. By contrast, the repulsive $H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$78 channel can have loss rates as low as

$H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$79

whereas the universal loss benchmark for CaF is

$H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$80

The suppression therefore reaches about five orders of magnitude, and the paper emphasizes that merely changing the fluorine nuclear spin state can switch between strongly dipolar and strongly shielded behavior.

The same mechanism appears broadly across species, but the optimal channel is species dependent. For MgF, CaF, SrF, and BaF, the paper mainly examines the $H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$81 channel. All show suppression at low temperature, but CaF gives the lowest loss among that series because it best balances tunneling to short range and inelastic hyperfine-changing loss. YO is qualitatively different because its hyperfine constants have opposite signs, giving an inverted hyperfine structure relative to CaF. In YO the repulsive channel is $H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$82, not $H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$83, and predicted rates can be as low as

$H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$84

about eight orders of magnitude below universal.

The paper also gives a scaling law in terms of the hyperfine energy scale relative to the dipolar energy scale, with

$H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$85

In the $H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$86 limit, the loss rate becomes a universal function of $H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$87. This explains why the effect is strong in open-shell molecules with hyperfine splittings of tens to hundreds of MHz, but not in typical closed-shell bialkali molecules with hyperfine splittings of only tens of kHz. Magnetic fields up to about $H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$88 G preserve the suppression in CaF for the repulsive $H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$89 channel, with the lowest Zeeman component $H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$90 optimal below $H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$91 G. The proposed experimental test is a merged-optical-tweezer collision experiment in which one prepares one molecule in $H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$92, the other in $H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$93 in different hyperfine states, merges the tweezers, and measures survival or loss.

6. Unifying structure, observables, and limitations

Taken together, these studies indicate a unifying structural principle: hyperfine degrees of freedom decide whether the dominant intermolecular effect is first-order or second-order, and therefore whether the observable scales as $H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$94, $H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$95, a contact-shift-like mean field, or a hyperfine-frequency correction. In xenon-pentane, fast exchange converts an underlying scalar coupling into a xenon-magnetization-dependent proton frequency shift. In the hydrogen $H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$96-$H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$97 problem, hyperfine resolution decides which channels are resonant and which remain dispersive. In the hydrogen $H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$98-$H_J = h J\, \mathbf I_1 \cdot \mathbf I_2,$99 problem, the hyperfine interaction appears as a Dirac-$1/R^3$00 modification of the dispersion potential and of the $1/R^3$01 hyperfine interval. In open-shell molecules, hyperfine selection rules can eliminate first-order resonant dipolar exchange altogether, leaving a second-order hyperfine van der Waals interaction that may be repulsive (Ledbetter et al., 2011, Jentschura et al., 2017, Adhikari et al., 2017, Walraven et al., 16 Jul 2025).

The observables are correspondingly different. The liquid-state scalar-coupling experiment measures a line shift linear in xenon longitudinal magnetization, not a resolved splitting. The hydrogen calculations predict state-resolved level shifts and hyperfine-interval shifts in beam or gas environments. The ultracold-molecule analysis predicts adiabatic potential barriers, tunneling suppression, inelastic loss, and hyperfine-state-dependent collision rates. This suggests that “hyperfine van der Waals interaction” is an umbrella for several Hamiltonian reductions rather than a single measurement protocol.

The limitations are equally system specific. In the xenon-pentane experiment, the measured $1/R^3$02 is an ensemble average over rapidly exchanging configurations, and the analysis assumes fast exchange and near-parallel xenon polarization. The structural averaging in

$1/R^3$03

neglects explicit angular structure, even though the underlying DFT results show orientation dependence. In the $1/R^3$04-$1/R^3$05 hydrogen treatment, the neglect of the $1/R^3$06 manifold requires $1/R^3$07, and the interaction is intrinsically quasi-degenerate rather than an ordinary nondegenerate van der Waals calculation. In the $1/R^3$08-$1/R^3$09 problem, the most distinctive hyperfine correction, the $1/R^3$10 term, belongs only to the intermediate regime $1/R^3$11; in the extreme Lamb-shift range the effect is negligible. In open-shell molecules, the shielding depends on species-specific hyperfine constants and on the ratio $1/R^3$12, and residual loss may be dominated either by tunneling to short range or by inelastic coupling to other hyperfine states.

A final point of interpretation follows directly from the cited works. Hyperfine van der Waals interactions are not confined to covalent frameworks, not necessarily attractive, and not necessarily small. They can arise from second-order Fermi-contact physics across transient van der Waals contact, from hyperfine-resolved quasi-degenerate dipole-dipole coupling in atoms, or from hyperfine-selected second-order dipolar interactions in ultracold molecules. What unifies these cases is that hyperfine structure controls the effective intermolecular coupling channel, and therefore the observable long-range physics.