Field-linked resonances of polar molecules

Scattering resonances are an essential tool for controlling the interactions of ultracold atoms and molecules. However, conventional Feshbach scattering resonances1, which have been extensively studied in various platforms1–7, are not expected to exist in most ultracold polar molecules because of the fast loss that occurs when two molecules approach at a close distance8–10. Here we demonstrate a new type of scattering resonance that is universal for a wide range of polar molecules. The so-called field-linked resonances11–14 occur in the scattering of microwave-dressed molecules because of stable macroscopic tetramer states in the intermolecular potential. We identify two resonances between ultracold ground-state sodium–potassium molecules and use the microwave frequencies and polarizations to tune the inelastic collision rate by three orders of magnitude, from the unitary limit to well below the universal regime. The field-linked resonance provides a tuning knob to independently control the elastic contact interaction and the dipole–dipole interaction, which we observe as a modification in the thermalization rate. Our result provides a general strategy for resonant scattering between ultracold polar molecules, which paves the way for realizing dipolar superfluids15 and molecular supersolids16, as well as assembling ultracold polyatomic molecules.

Scattering resonances are an essential tool for controlling the interactions of ultracold atoms and molecules. However, conventional Feshbach scattering resonances 1 , which have been extensively studied in various platforms [1][2][3][4][5][6][7] , are not expected to exist in most ultracold polar molecules because of the fast loss that occurs when two molecules approach at a close distance [8][9][10] . Here we demonstrate a new type of scattering resonance that is universal for a wide range of polar molecules. The so-called field-linked resonances [11][12][13][14] occur in the scattering of microwave-dressed molecules because of stable macroscopic tetramer states in the intermolecular potential. We identify two resonances between ultracold ground-state sodiumpotassium molecules and use the microwave frequencies and polarizations to tune the inelastic collision rate by three orders of magnitude, from the unitary limit to well below the universal regime. The field-linked resonance provides a tuning knob to independently control the elastic contact interaction and the dipole-dipole interaction, which we observe as a modification in the thermalization rate. Our result provides a general strategy for resonant scattering between ultracold polar molecules, which paves the way for realizing dipolar superfluids 15 and molecular supersolids 16 , as well as assembling ultracold polyatomic molecules.
Ultracold polar molecules with tunable dipole moments provide a powerful platform for quantum simulations 17,18 , quantum computation 19,20 and ultracold chemistry 21 . Scattering resonances are long-sought-after tools in these systems, which have been essential in ultracold-atom experiments to control the contact interaction and for creating strongly correlated quantum phases 22 , as well as for producing ultracold diatomic molecules 1 . Independent control over contact and long-range interactions in ultracold molecules has been predicted to enable the realization of new quantum phenomena such as exotic self-bound droplets and supersolid quantum phases 16 . Moreover, measurements of scattering resonances provide an accurate benchmark for calculations of the molecular potential energy surface 3,21 and open a new route in controlled quantum chemistry 5 .
A scattering resonance occurs when the scattering state strongly couples to a quasibound state. Based on whether the quasibound state is hosted by the same or a different channel than the scattering channel, the resonance is categorized as a shape resonance or a Feshbach resonance, respectively. Shape and Feshbach resonances have been observed in atom-molecule and molecule-molecule collisions by scanning the collision energy by using molecular beams at kelvin and subkelvin temperatures 21,[23][24][25][26] . In the ultracold (submicrokelvin) regime, scattering resonances are often induced by an external electromagnetic field that shifts the relative energy between the quasibound state and the scattering state 1 . Magnetically tunable Feshbach resonances have been observed in collisions between weakly bound Feshbach molecules 2,4 and recently between NaLi molecules in the spin-triplet ground state 6 . However, the magnetic tuning scheme essential to Feshbach resonances requires a non-zero electronic spin, and is, therefore, unlikely to find application for bialkali molecules in the spin-singlet ground state. The spin-singlet absolute ground state of bialkaline molecules is of special interest, as it is the only long-lived state in which the molecules feature strong electric dipole-dipole interactions (DDI). Moreover, Feshbach resonances are not expected to occur between ground-state molecules in the presence of nearly universal loss, owing to the high density of tetramer states near the collisional threshold and the loss mechanisms associated with collisional complexes [8][9][10] . A general method to realize collisional resonances of ultracold dipolar molecules, therefore, remains open.
Here we demonstrate a general approach to create such resonances in collisions between dipolar molecules by coupling them to so-called field-linked (FL) states 11,12 . These weakly bound states are induced by engineering an attractive well in the long-range intermolecular potential through microwave dressing 14,27 . Unlike in conventional resonances, in which an external field merely tunes an existing short-range quasibound state into resonance, the long-range FL states exist only in the presence of the microwave field. The sensitivity of the FL states to the microwave field leads to an unprecedented level of control over the intermolecular interaction. Here we demonstrate this tunability by observing two resonance branches in the inelastic scattering rate, the peak positions of which continuously shift with the microwave frequency and polarization. We further characterize the change of the thermalization rate caused by the diverging scattering length in a resonant collision channel.

Interaction potential
Polar molecules possess a permanent dipole moment d 0 in their body-fixed frame. To induce a dipole moment in the laboratory frame, external fields need to be applied to mix different rotational states and to break the rotational symmetry. Here we use microwave dressing between the two lowest rotational states of the molecules to polarize them. The induced dipole moment follows the alternating current electric field as , where ξ describes the ellipticity of the microwave radiation and e ± (t) are the σ ± polarization basis vectors. The time-averaged dipole moment d d ΔΩ = / 6(1 + ( / ) ) 0 2 is tunable via the microwave detuning Δ and the Rabi frequency, Ω. Because the rotational frequency of NaK at 5.643 GHz is much faster than all other dynamical timescales in the system, we consider the time-averaged DDI at long range 15 where r = (r, θ, φ) is the relative position between the molecules in polar coordinates defined by the microwave wavevector. The microwave propagates in the z direction, and the microwave polarization ellipse has its major and minor axes in the y and x direction. Remarkably, the symmetry of the interaction can be manipulated by the microwave ellipticity, as illustrated in Fig. 1a-c. For ξ = 0° (circular polarization) and ξ = 45° (linear polarization), where U dd resembles the typical DDI up to a constant prefactor 28 . In between linear and circular polarization, the interaction breaks the rotational symmetry along all directions.
As the molecules approach each other, microwave dressing induces an anisotropic van der Waals interaction U vdW ∝ 1/r 6 (Methods) 15 . With blue-detuned microwave dressing, U vdW is repulsive in all directions, which protects the molecules from loss processes at short range and reduces the inelastic cross-sections 14,[29][30][31] . Such a shielding potential arises owing to an avoided crossing between the attractive and the repulsive branch of the DDI. In a semi-classical picture, this avoided crossing can be understood as the reorientation of colliding dipoles through DDI 31,32 . A similar flipping of the dipoles occurs between polar molecules in a direct current electric field 33,34 , between Rydberg atoms 35 and between ions and Rydberg atoms 36 .
The full interaction potential between two dressed molecules is the sum of the DDI potential and the van der Waals potential U(r) = U dd (r) + U vdW (r). We can shape the interaction potential and control the scattering process, as illustrated in Fig. 1. Along the y axis, U(r) resembles a Mie potential 15,37,38 with a characteristic length of about 10 3 a 0 . A deviation from circular polarization breaks the azimuthal symmetry of the DDI and enhances the depth of the potential well along the y axis, which becomes deep enough to support one or more bound states 13 . These bound states are the FL states, the properties of which strongly depend on the external fields. By tuning the binding energy of the FL state across the collisional threshold, for example, with microwave detuning as shown in Fig. 1d,e, FL resonances occur, which drastically alter the scattering properties between the molecules.
The low-energy scattering in U(r) can be described by the associated partial-wave phase shifts, which are given by 39 Here l and m l are the angular momentum and its projection along the quantization axis, k is the relative wavevector, a lm l and a dd are the characteristic lengths associated with the contact interaction and the DDI respectively, and c lm l denote partial-wave dependent prefactors for the dipolar scattering phase shifts. Note that the contact interaction is suppressed by the centrifugal barrier for l ≠ 0 as k → 0, whereas the phase shift from the long-range DDI scales linearly with k in all partial waves, and is proportional to the dipolar length a μd π = /4 ϵ

Resonance map
We map out the resonances by measuring the inelastic rate coefficient β in of collisions between the dressed molecules. The optically trapped ground-state 23 Na 40 K molecules with nuclear spin projections (m i,Na , m i,K ) = (3/2, −4) are formed from an ultracold atomic mixture by means of magnetoassociation and subsequent stimulated Raman adiabatic passage (STIRAP) at a magnetic offset field of 72.35 G (ref. 40 ). For most measurements, the temperature T of the molecular ensemble is 230 nK and the initial average density n 0 is about 5 × 10 11 cm −3 . Next, the microwave is ramped on in 100 μs to dress the molecules. After a variable hold time, the remaining molecules are released from the optical trap and we determine the number of molecules and the temperature from the time-of-flight images.

= 19°
= 45°y x Intermolecular distance (a 0 ) 0 -5,000 5,000 Intermolecular distance (a 0 ) 0 -5,000 5,000 Intermolecular distance (a 0 ) 0 -5,000 5,000 As we tune the ellipticity of the microwave from circular to linear, up to two FL states emerge from the dressed potential. Figure 2a shows exemplary loss rate coefficients for three different ellipticities. At ξ = 6(2)°, the potential is too shallow to support a bound state, thus no resonance is observed. In this regime, the inelastic collision is suppressed by the shielding potential at small detunings 31 . For ξ = 19(2)°, the interaction potential supports a single bound state near zero microwave detuning, leading to enhanced inelastic scattering at Δ ≈ 2π × 10 MHz. For ξ = 37(2)°, the potential becomes deep enough to support two bound states, leading to two resonance peaks.
One special feature of the FL resonances is their sensitivity to external fields. We show that the resonance position continuously changes with the microwave parameters by mapping out the two resonance branches while varying the microwave detuning and polarization. Figure 2b shows two branches of FL resonances, starting at ξ ≈ 10° and ξ ≈ 32°. As the polarization ellipticity ξ increases, less DDI is needed to support the bound states and the resonances, therefore, shift to larger detuning. However, the global inelastic rate coefficient increases as the polarization becomes more elliptical, because of the increased coupling to the other dressed states 28 . Overall our measurements show good agreement with our theory predictions (Fig. 2c). We attribute the broadening and shift of the resonance peaks compared to the theory to an increase of the Rabi frequency as we scan the detuning (Methods). These systematic errors affect mostly the FL resonances near linear polarization, where the potential depth is more sensitive to the relative detuning.

Temperature dependence of the inelastic scattering
The temperature dependence of the inelastic scattering rate varies with the detuning. At large detuning where the DDI is reduced, the inelastic scattering rate is universal 41 and scales for identical fermions linearly with temperature. At small detuning, the scattering enters the semi-classical regime for which β in is independent of the temperature 42 . On the scattering resonance, the collision rate has a temperature dependence that is reminiscent of the unitarity limit 5,32 , whereas the loss remains substantially smaller than this limit because of shielding. Meanwhile, the width of the resonance feature is broadened by thermal averaging. As a consequence, for temperatures as high as 700 nK, the resonance becomes less visible as shown in Fig. 3. When the collision energy becomes lower than the centrifugal barrier of the interaction potential, the resonance peak would be further narrowed owing to the increased lifetime of the quasibound state. Therefore, reaching ultracold temperatures is crucial for the observation of FL resonances.

Elastic scattering
Scattering resonances are associated not only with enhanced losses of the molecules, but more importantly, the ability to control elastic scattering. With FL resonances, we can tune the elastic scattering rate while keeping the inelastic rate small.
We characterize the effect of the FL resonances on the elastic collision rate by measuring the thermalization rate of the samples. This is commonly done by quenching the trapping confinement in one dimension and observing the global cross-dimensional thermalization 43 . However, for small Δ/Ω our samples are in the hydrodynamic regime, in which the global thermalization rate is limited by the trapping frequencies 31 . Instead, we perturb the momentum distribution of the molecular cloud by pulsing on an optical-lattice beam for 300 μs. The lattice pulse diffracts some molecules and sends them to collide along the y axis, defined as the long-axis of the microwave field ( Fig. 1), along which the DDI is most attractive. Fast local thermalization smears out the diffraction pattern that is formed in momentum space during Article the lattice pulse (insets in Fig. 4a). From the contrast of the diffraction pattern we can estimate the thermalization rate Γ th (Methods). We observe that the measured thermalization rate follows a similar trend as the calculated value of the elastic collision rate coefficient, β el (Fig. 4a). Besides the contribution from the DDI, which scales with d 2 and decreases as Δ increases, a clear resonance feature is visible around the FL resonance. The shift of the resonance between the experiment data and theory is within the systematic uncertainty of the ellipticity. The average number of elastic collisions that is needed per particle to reach thermalization is so far unknown in the regime of elliptical microwave polarization. From the comparison between the measured Γ th /n and β el from our coupled-channel calculations, we find that this factor is close to 1 under the present experimental conditions.
The observed elastic scattering rate has contributions from both contact interaction and long-range DDI. As shown in equation (2), the DDI contributes in multiple channels, whereas the contact interaction only has a notable contribution in the resonant channel, owing to its unfavourable scaling with the wavevector. Although such an interplay between contact interaction and DDI limits the change in the total elastic scattering rate, the underlying scattering length, however, shows divergent behaviour in the resonant channel. Figure 4b shows the energy-dependent scattering length a k δ k k ( ) = −tan ( )/ lm lm l l ∼ (refs. 41,44 ) for the three p-wave channels at the average collision energy. The real (imaginary) part of the scattering length corresponds to the elastic (inelastic) scattering. The FL resonance occurs in the p y channel where the interaction is most attractive. The corresponding scattering length shows a resonance feature, in which the real part can be tuned to large positive or negative values, whereas the imaginary part remains small. The ratio of elastic-to-inelastic collisions is about ten on the resonance, and can be further enhanced at smaller ellipticity, in which the resonance shifts towards higher Rabi frequencies (Methods).
A simple analytic formula that describes the resonant elastic scattering length is given by where Δ 0 (k) and Δ * (k) denote the position and the width of the resonance. The width Δ * (k) ∝ k 2 follows the scaling of the p-wave contact interaction. For the collision energy considered here, we extract Δ 0 ≈ 2π × 10.99 MHz and Δ * ≈ 2π × 3.29 MHz from the fit to the coupled-channel calculations.
The resonance position Δ 0 (k) also has a weak energy dependence. As a consequence, the resonance position with thermal averaging is broadened and slightly shifted towards lower detuning. These thermal effects, however, will be greatly suppressed in a degenerate Fermi gas, in which the scattering predominantly occurs near the sharp Fermi energy 15 .

Discussion
Field-linked resonances provide a new universal tool to control the collisions between ultracold polar molecules. These resonances occur as long as the Rabi frequency is sufficiently large, such that the interaction potential is deep enough to support the bound states. Quantitative predictions of FL resonances only require knowledge of the mass, the dipole moment and the rotational structure of the individual molecules as well as their loss rate at short range. This is in stark contrast to molecular collisions involving close contact between the molecules, in which a large number of collision channels are involved and the existing knowledge of the potential energy surface are too imprecise to predict the number of bound states, let alone their position. Precise knowledge on the FL states also makes them useful as intermediate states in photoassociation spectroscopy to probe the short-range potential 11 . The control over the scattering length opens up new possibilities to investigate many-body physics with both contact interaction and DDI. In a degenerate Fermi gas, the resonant interaction facilitates realization of dipolar superfluidity 27,45,46 . Specifically, pairing between molecules is enhanced because of the presence of the FL bound state. Therefore, the critical temperature for Bardeen-Cooper-Schrieffer superfluidity increases drastically near the FL resonance, that is, to about 14% of the Fermi temperature for NaK molecules 15 . The anisotropic nature of such a dipolar superfluid gives rise to new quantum phenomena such as gapless superfluidity 47 and topological p x + ip y symmetry 48 . In a Bose-Einstein condensate (BEC), independent control over the s-wave scattering length and dipolar length has led to the observation of self-bound droplets and the formation of supersolids in magnetic atoms 49 . Making use of FL resonance with bosonic polar molecules, the dipolar lengths of which are orders of magnitude larger than magnetic atoms, would enable the study of such exotic phenomena in entirely unexplored regimes 16 .
The observed resonances also demonstrate the existence of the FL states, a new set of exotic long-range polyatomic molecular states. Tetramer molecules with approximately twice the dipole moment of the individual diatomic molecules could potentially be created by adiabatically ramping the microwave field across a FL resonance or by radio-frequency association. Those composite bosonic tetramers are expected to be long lived at small binding energies 12 and could be collisionally stable owing to the shielding between the constituent molecules. Below the critical temperature, a tetramer gas could form a BEC 50 and may lead to a new crossover from a dipolar Bardeen-Cooper-Schrieffer superfluid to a BEC of tetramers.

Conclusion
We have observed a new type of universal scattering resonance between ultracold microwave-dressed polar molecules that is associated with FL tetramer bound states in the long-range potential well. The resonances are highly tunable via microwave power, frequency and polarization, which makes them a versatile tool for controlling molecular interactions. Because the FL states are insensitive to species-dependent short-range interactions, the FL resonance is applicable to a wide range of polar molecules. Our results provide a general route to strongly interacting molecular gases and open up new possibilities to investigate new quantum many-body phenomena and to produce long-lived dipolar tetramer molecules.

Online content
Any methods, additional references, Nature Portfolio reporting summaries, source data, extended data, supplementary information, acknowledgements, peer review information; details of author contributions and competing interests; and statements of data and code availability are available at https://doi.org/10.1038/s41586-022-05651-8.

Microwave field generation
We use a dual-feed square-waveguide antenna to generate a microwave field with tunable polarization. The waveguide is fabricated from copper-coated glass-fibre-reinforced epoxy laminate. The inner dimensions of the waveguide are 33 × 33 × 58 mm. The width of the waveguide is chosen such that the cut-off frequency of the transverse electric (TE 10 ) mode is below the rotational transition frequency of the molecules at 5.64 GHz and that the microwave field strength at the position of the molecules is optimal. To achieve impedance matching, the two feeds are 13 mm long, which is a quarter wavelength in free space. They are placed orthogonal to each other, and 22 mm away from the back plate of the waveguide, which is about a quarter wavelength in the waveguide. Each feed produces a close to linearly polarized electric field at the position of the molecules, which is about 25 mm away from the top of the waveguide. When the field strengths of these subfields are balanced, tuning the relative phase of the feeds allows straightforward tuning of the field polarization. A relative phase of approximately 90° (−90°) results in a σ + (σ − ) polarized field, whereas approximately 0° (180°) produces a linearly polarized field along the y (x) direction.
The microwave setup is sketched in Extended Data Fig. 1. A Rohde & Schwarz SMA100B signal generator (using a yttrium-iron-garnet oscillator) with the noise suppression option SMAB-B711 is used as microwave source. To independently control the two feeds of the antenna, the microwave is divided by a power splitter into two paths. Each path includes a voltage-controlled attenuator to balance the subfields, a 10-W amplifier (KUHNE electronic KU PA 510590-10 A) and a mechanical phase shifter (SHX BPS-S-6-120) for differential phase control.

Calibration of the field polarization
To characterize the polarization of the microwave field, we initially measure its π, σ + and σ − component in the frame of the magnetic offset field. We do this at low microwave power by measuring the Rabi coupling to excited rotational states with different projections m J of the rotational quantum number J on the magnetic field axis 31 . A direct measurement of the σ + and σ − component at large microwave power is impracticable, as the Rabi frequency is then orders of magnitude larger than the Zeeman splitting of the m J states. When we change the relative phase ϕ between the feeds, each field component, and thus each Rabi frequency, individually oscillates with a period of 360° owing to the interference between the imperfect subfields, as shown in Extended Data Fig. 2a. We fit these oscillations with the function with the fit parameters Ω 1 , Ω 2 and ϕ 0 , which define the contributions of the individual feed. The offset phases ϕ 0 have an uncertainty of 2.9°, which we attribute to the hysteresis and imperfect tuning of the mechanical phase shifters. Note, for the calibration measurements presented in Extended Data Fig. 2a, the power balance between the feeds was tuned to minimize the σ − component around ϕ = 90°. The finite ellipticity of the subfields causes the field strengths at other angles of ϕ to be unbalanced, so that we do not get pure σ + polarization at ϕ = −90° and the phase values that provide linear polarization shift away from 0° and 180°. We calculate the ellipticity ξ from the fitted Rabi frequencies for each individual relative phase, as shown in Extended Data Fig. 2b. At high microwave power, the electric field of the microwave replaces the magnetic offset field as the quantization axis, so that ξ is defined in the frame of the microwave field. Because the relative phase between the measured field components is unknown, we can only deduce the limits for the tilt angle of the microwave wavevector with respect to the magnetic offset field, which causes a systematic uncertainty of ξ. Close to σ + polarization the tilt angle is only around 10°, so that the uncertainty of ξ is dominated by the above-mentioned uncertainties of ϕ 0 .
When switching from the calibration measurements at low microwave power to the measurements at high power, we have to consider the non-linearity of the amplifiers close to their saturation power. This shifts the power balance between the two antenna feeds when we increase the microwave power and thereby changes ξ. We reestablish the power balance by tuning one voltage-controlled attenuator while optimizing the shielding at ϕ = 90°, where we initially minimized ξ. We find that we need to compensate the relative power by 10%. In addition, we also observe a small variation of the relative phase and power as we scan the microwave detuning. The variation in the relative phase leads to a systematic uncertainty of the ellipticity on the same order as the contribution from the uncertainty of ϕ 0 . The variation in Rabi frequency could broaden the resonance feature, especially for ellipticities close to linear polarization, as observed in Fig. 2.

Inelastic collision rate coefficient
The inelastic collision rate coefficients β in are experimentally determined from the time evolution of the measured molecule number N and the temperature T by numerically solving the differential equations 51 where ω is the geometric mean trapping frequency and Γ 1 is the one-body loss rate. We assume both the heating rate and the two-body loss rate coefficient to be constant during the loss process. As the overall heating is no more than 50%, the fitted values of β in do not significantly change when we instead assume a linear temperature dependence of the rate coefficient.
The comparison between loss measurements on the scattering resonance and away from resonance are shown in Extended Data Fig. 3a. To limit the number of free fit parameters, we determine Γ 1 = 0.53(2) Hz in independent measurements at low densities, as shown in Extended Data Fig. 3b.

Fast thermalization measurements
The momentum distribution of the molecular cloud is disturbed by pulsing on a one-dimensional optical lattice for t lat = 300 μs. Subsequently the microwave power is ramped down in 100 μs and the trapping confinement is turned off. After 10 ms time of flight, the momentum distribution is imaged. The lattice is approximately aligned along the y direction. Its lattice constant is a lat = 532 nm and the lattice depth is 88E r , where E h ma = /(8 ) r 2 lat 2 is the lattice recoil energy. The pulse duration was chosen to optimize the contrast of the resulting interference pattern. Note that the pulse is short compared to the trap frequencies 2π × (82, 58, 188) Hz of the background confinement, so that cross-talk between the momentum distribution and the real-space density is small. Also two-body loss, even on resonance, is negligible on this time scale.
To estimate the degree of thermalization c th we make the simplifying assumption that the momentum distribution after the lattice pulse can be described by is a Gaussian (that is, thermal) distribution and n v ( ) 0 ∼ is the undisturbed interference pattern (Extended Data Fig. 4a). We determine ∼ n v ( ) 0 by averaging over 16 images of the interference pattern in the absence of microwave-induced interactions. We further assume that the interference pattern decays exponentially, such that the thermalization rate is given by Given the average number of elastic collisions that is required to reach thermal equilibrium N col and the average relative velocity v rel of the colliding molecules, the elastic scattering cross-section is given by Assuming that the thermalization is mainly driven by elastic collisions between molecules in the side peaks and molecules in the main peak of the diffraction pattern, the average relative velocity can be approximated as is the thermally averaged collision velocity in the undisturbed sample. The expression of N col is so far unknown in the regime of elliptical microwave polarization. We can, however, set an upper limit by assuming ξ = 0. In that case, N col ≈ 4 for a tilt angle of the microwave field of 10° (refs. 34,43 ). A comparison of this simplified model with our predictions of β vσ = el el is shown in Extended Data Fig. 4b.

Born approximation
The elastic scattering rate (Fig. 4a)  The scattering length given in equation (3) is a y

,Born
∼ plus the contribution from the FL resonance.

Resonance near circular microwave polarization
The elastic-to-inelastic collision ratio near the FL resonances can be improved with better shielding near the circular microwave polarization and with lower temperature. The ratio is then about 900 at the maximum elastic scattering rate and about 130 at the scattering resonance, as shown in Extended Data Fig. 5. Under these conditions, the FL resonance occurs at a much higher Rabi frequency compared to the observed resonances at more elliptical polarizations. However, this is still realistic to achieve by using an improved antenna design with increased microwave power. For Ω = 2π × 10 MHz and zero detuning, C 6i = 1.1 × 10 −72 kg m 8 s −1 is orders of magnitude larger than the direct isotropic van der Waals coefficient C 6 = −4.9 × 10 −74 kg m 8 s −1 (ref. 52 ). The induced C 6i reduces as the detuning increases and can be comparable to C 6 for Δ/Ω > 2.7.

Coupled-channel calculations
We perform coupled-channel scattering calculations by using the framework described in detail in refs. 28,29 . Here we summarize the numerical details of these calculations.
The NaK molecules are described as rigid rotors with the rotational states J = 0, 1. Interactions of the molecules with elliptically polarized microwaves are included as described in ref. 53 . The wavefunction for the relative motion is expanded in partial waves L = 1, 3, 5. Hyperfine interactions are not included as these were previously found to have a negligible effect 31 when operating at sufficiently large magnetic fields 29,53 . The colliding molecules interact with one another through the DDI, and undergo short-range loss, which is modelled as a capture boundary condition imposed at r = 20a 0 . We propagate the coupled-channel equations outwards to r = 10 6 a 0 and match the solution to the scattering boundary conditions. This yields the scattering matrix, from which collision cross-sections and rate coefficients are determined and thermally averaged by using an energy grid of 21 energies that are logarithmically spaced between 0.03k B T and 32k B T.
Binding energies of the FL bound states are calculated as follows. First, we compute adiabatic potentials by diagonalizing the Hamiltonian described above excluding the radial kinetic energy for fixed values of the molecule-molecule separation r. On each adiabatic potential curve, we compute bound states by using a sinc-function discrete variable representation 54 . We find that the position of the zero-energy bound states computed in this approximation agree well with the resonance positions, indicating that one can think of the FL bound states as living on a single adiabatic potential curve. Both resonances found here are supported by the lowest adiabatic potential, that is, the second resonance corresponds to a radial vibrational excitation, rather than an angular excitation. Interaction potentials that are shown in Fig. 1 are similarly computed as adiabatic potential curves, except that these are computed for a fixed orientation of the intermolecular axis relative to the microwave polarization, rather than by using a partial wave expansion.
Note, the coupled-channel calculations have no free parameters.

Data availability
The experimental data that support the findings of this study are available from the corresponding author upon reasonable request. Source data are provided with this paper.