Resonant Ideal in Constrained Systems

• Resonant Ideal is a conceptual framework describing systems where resonance emerges under maximally constrained conditions, revealing otherwise hidden responses.
• Practical examples include radiationless anapole modes in silicon nanodisks, edge-bound states in photonic crystals, and ideal MHD responses in fusion plasmas.
• Applications in devices such as graphene sensors and cold-atom experiments show that even near ideal limits, non-ideal factors like disorder and dissipation are crucial.

“Resonant ideal” is not a single standardized term with one fixed definition; across the cited literature, it denotes a family of constructions in which resonance is examined in an ideal limit, or in which an ideally selective resonant response is isolated. In that usage, “ideal” may mean a disorder-free resonant photonic crystal, a smooth incompressible ideal-MHD evolution, a ballistic graphene nanoribbon without defects, a high-vacuum nanoresonator whose dissipation is largely extrinsic, or an ultracold gas whose background interactions are negligibly small until an optical resonance is applied. Resonance, correspondingly, may refer to a radiationless optical state, a rational-surface response, a quasi-bound transport channel, or a mechanically or spectroscopically shifted eigenfrequency (Skipetrov, 2019, Pfefferlé et al., 2019, Ihnatsenka, 2022, More et al., 2023, Blatt et al., 2011).

1. Conceptual scope

A recurrent structure across these works is the coexistence of two ingredients: a simplified or symmetry-constrained baseline system, and a resonant perturbation that reveals an otherwise hidden response. In optics, the canonical example is the anapole, where destructive interference cancels the electric-type dipolar far field while large internal fields remain. In tokamak or stellarator MHD, the corresponding ideal limit forces the resonant normal field to vanish at a rational surface, with singular or sharply localized shielding currents carrying the response. In cold-atom scattering, an almost ideal gas acquires a complex scattering length only when light couples the continuum to a near-threshold molecular level (Lu et al., 2021, Pharr et al., 18 Mar 2026, Blatt et al., 2011).

The literature also shows that “ideal” rarely means lossless or trivial. Instead it usually denotes a constraint class. In the ideal diamond photonic crystal, the central observable is the density of states inside a band gap of a finite, perfectly ordered lattice. In ideal MHD, frozen-in flux and force balance constrain accessible equilibria. In ideal graphene nanoribbons, ballistic transport and defect-free interfaces do not suppress resonant behavior; they make the interface-induced quasi-bound states more transparent. This suggests that “resonant ideal” is best understood as a cross-disciplinary descriptor for resonance under a maximally constrained background model, rather than as a universal formalism (Skipetrov, 2019, Pfefferlé et al., 2019, Ihnatsenka, 2022).

2. Optical and electromagnetic realizations

In open optical resonators, an anapole is a radiationless field distribution that can still be externally excited. Its lowest-order condition is the cancellation of electric dipole and toroidal dipole radiation,

$\mathbf{p} + i k \mathbf{T} = 0,$

with

$$\mathbf{p} = \frac{i}{\omega} \int \mathbf{J}(\mathbf{r})\, d^3 r,\qquad \mathbf{T} = \frac{1}{10 c} \int \big[\mathbf{r}(\mathbf{r}\cdot\mathbf{J}) - 2 r^2 \mathbf{J}\big]\, d^3 r.$$

The distinctive result of the silicon-nanoparticle study is that a radiationless anapole and a magnetic resonant state can coexist at the same wavelength and be selected purely by the polarization of a tightly focused cylindrical vector beam. For a Si nanodisk with radius $r=150$ nm and height $h=160$ nm on glass, radially polarized excitation yields an anapole near $\lambda \approx 720$ nm, while azimuthally polarized excitation at the same wavelength produces a resonant quadrupolar state. Experimentally, the AP-excited resonance appears near $735$ nm, and at that wavelength the RP-excited back-scattering is about one order of magnitude smaller than under AP excitation at equal power. The underlying selection rule is symmetry based: tightly focused RP beams satisfy $p_{Ml}^{m}=0$, whereas AP beams satisfy $p_{El}^{m}=0$ in the focal multipole expansion, so the illumination itself selects electric-type or magnetic-type channels (Lu et al., 2021).

A more singular electromagnetic realization appears in moving Fabry–Perot cavities built from metallic plates separated by a nanoscale vacuum gap. For evanescent waves, propagation across the gap contributes amplitude decay $e^{-2\kappa L}$ but no phase advance, so the resonance condition is carried entirely by the reflection phases and amplitudes. Relative motion Doppler-shifts the frequency seen by the moving plate to $\omega' \approx \omega - v k_x$, and at the phase-balance wavevector $$\mathbf{p} = \frac{i}{\omega} \int \mathbf{J}(\mathbf{r})\, d^3 r,\qquad \mathbf{T} = \frac{1}{10 c} \int \big[\mathbf{r}(\mathbf{r}\cdot\mathbf{J}) - 2 r^2 \mathbf{J}\big]\, d^3 r.$$0 one has $$\mathbf{p} = \frac{i}{\omega} \int \mathbf{J}(\mathbf{r})\, d^3 r,\qquad \mathbf{T} = \frac{1}{10 c} \int \big[\mathbf{r}(\mathbf{r}\cdot\mathbf{J}) - 2 r^2 \mathbf{J}\big]\, d^3 r.$$1, so the moving reflection coefficient becomes the complex conjugate of the stationary one. Near the surface-plasmon resonance, $$\mathbf{p} = \frac{i}{\omega} \int \mathbf{J}(\mathbf{r})\, d^3 r,\qquad \mathbf{T} = \frac{1}{10 c} \int \big[\mathbf{r}(\mathbf{r}\cdot\mathbf{J}) - 2 r^2 \mathbf{J}\big]\, d^3 r.$$2 for p-polarized evanescent waves, and motion supplies the gain needed to satisfy

$$\mathbf{p} = \frac{i}{\omega} \int \mathbf{J}(\mathbf{r})\, d^3 r,\qquad \mathbf{T} = \frac{1}{10 c} \int \big[\mathbf{r}(\mathbf{r}\cdot\mathbf{J}) - 2 r^2 \mathbf{J}\big]\, d^3 r.$$3

The resulting critical scales are

$$\mathbf{p} = \frac{i}{\omega} \int \mathbf{J}(\mathbf{r})\, d^3 r,\qquad \mathbf{T} = \frac{1}{10 c} \int \big[\mathbf{r}(\mathbf{r}\cdot\mathbf{J}) - 2 r^2 \mathbf{J}\big]\, d^3 r.$$4

In this regime the scattering pole reaches the real axis, the quality factor diverges formally, and near-field heat transfer is predicted to grow as $$\mathbf{p} = \frac{i}{\omega} \int \mathbf{J}(\mathbf{r})\, d^3 r,\qquad \mathbf{T} = \frac{1}{10 c} \int \big[\mathbf{r}(\mathbf{r}\cdot\mathbf{J}) - 2 r^2 \mathbf{J}\big]\, d^3 r.$$5. The paper is equally explicit that nonlocality and nonlinearities regularize the singularity (Guo et al., 2013).

3. Spectral density, point scatterers, and superdimensional resonators

In finite resonant photonic crystals, ideality enters through perfect order rather than through vanishing dissipation. The model studied in the diamond-lattice work consists of resonant point scatterers on two interpenetrating fcc lattices, treated with a $$\mathbf{p} = \frac{i}{\omega} \int \mathbf{J}(\mathbf{r})\, d^3 r,\qquad \mathbf{T} = \frac{1}{10 c} \int \big[\mathbf{r}(\mathbf{r}\cdot\mathbf{J}) - 2 r^2 \mathbf{J}\big]\, d^3 r.$$6 non-Hermitian Green’s matrix. For $$\mathbf{p} = \frac{i}{\omega} \int \mathbf{J}(\mathbf{r})\, d^3 r,\qquad \mathbf{T} = \frac{1}{10 c} \int \big[\mathbf{r}(\mathbf{r}\cdot\mathbf{J}) - 2 r^2 \mathbf{J}\big]\, d^3 r.$$7, the infinite crystal has a full band gap over

$$\mathbf{p} = \frac{i}{\omega} \int \mathbf{J}(\mathbf{r})\, d^3 r,\qquad \mathbf{T} = \frac{1}{10 c} \int \big[\mathbf{r}(\mathbf{r}\cdot\mathbf{J}) - 2 r^2 \mathbf{J}\big]\, d^3 r.$$8

For finite spherical fragments, the density of states is defined by

$$\mathbf{p} = \frac{i}{\omega} \int \mathbf{J}(\mathbf{r})\, d^3 r,\qquad \mathbf{T} = \frac{1}{10 c} \int \big[\mathbf{r}(\mathbf{r}\cdot\mathbf{J}) - 2 r^2 \mathbf{J}\big]\, d^3 r.$$9

and the central result is the finite-size law

$r=150$0

throughout the remaining band gap. The explanation is geometric: in-gap states are boundary-dominated, so their total weight scales as surface-to-volume, $r=150$1. Moderate positional disorder narrows the gap and increases fluctuations near the shifted edges, but the disorder-averaged mid-gap DOS still obeys the same $r=150$2 scaling as long as the gap remains open (Skipetrov, 2019).

A closely related nearly ideal resonant medium is the ultracold $r=150$3Yb slab probed on the $r=150$4 transition at $r=150$5 nm. Here “ideal” refers to frozen atomic positions, negligible Doppler dephasing on the measurement timescale, low-intensity linear response, and an effectively two-level cycling transition. The complex transmission is measured interferometrically as $r=150$6, and the data agree with first-principles coupled-dipole simulations up to $r=150$7. The generalized Beer–Lambert model overestimates the minimum transmission and phase amplitude by about $r=150$8 at higher densities, whereas the coupled-dipole model captures both the amplitudes and the mild asymmetry of the line shape. The paper also resolves an experimental controversy: apparent saturation and broadening in earlier direct-transmission measurements are attributed to off-axis scattering collected by the imaging system and to imaging noise, not to a failure of microscopic theory (Vatré et al., 2024).

Superdimensional metamaterial resonators push this logic into spectral design. In the ideal two-dimensional construction, the Helmholtz equation is

$r=150$9

which, under separation $h=160$0, becomes

$h=160$1

The associated counting function grows as $h=160$2 for $h=160$3, so a planar resonator behaves spectrally like a medium of dimension $h=160$4. In three dimensions, the analogous bulk construction yields $h=160$5. The same designs exhibit giant focusing through the singular Green’s function and highly eccentric sub-Riemannian ray geometry near the degeneracy set $h=160$6 (Greenleaf et al., 2014).

4. Ideal MHD, rational surfaces, and singular response

In toroidal MHD, the ideal resonant condition is set by rational surfaces, $h=160$7, or equivalently by $h=160$8. In the ideal limit, islands cannot open, so the resonant normal field at the rational surface must vanish and is replaced by pitch-aligned shielding currents. The recent resistive study formalizes two complementary measures of this response. The shielding-current metric is

$h=160$9

while the penetrated resonant field is

$\lambda \approx 720$0

Within asymptotically matched resistive equilibria, $\lambda \approx 720$1 until low-$\lambda \approx 720$2 saturation, whereas $\lambda \approx 720$3 remains close to its ideal meaning at high $\lambda \approx 720$4 but is influenced by global kink structure at lower $\lambda \approx 720$5. In a low-rotation ITER equilibrium, both metrics shift the dominant coupling spectrum toward lower poloidal mode number $\lambda \approx 720$6, and they predict optimal relative coil phasings that differ substantially from ideal-MHD expectations (Pharr et al., 18 Mar 2026).

The older rigidity analysis shows how restrictive the ideal limit can become near a resonant layer. For the Hahm–Kulsrud–Taylor slab equilibrium $\lambda \approx 720$7, the neutral line $\lambda \approx 720$8 is a degenerate critical manifold of the flux function. Under smooth incompressible ideal motion, advection obeys

$\lambda \approx 720$9

while force balance requires $735$0. If mirror symmetry across $735$1 is imposed, the near-identity generator vanishes to all orders: every Taylor coefficient of $735$2 at $735$3 is forced to zero, and mirror-symmetric boundary displacements of the form $735$4 are formally unsupported by smooth incompressible ideal motion that preserves force balance. The paper identifies the underlying obstruction as a $735$5 resonant denominator in the generating potential (Pfefferlé et al., 2019).

Direct numerical resolution of the HKT singularity confirms the same picture in a different language. In the slab with resonant wall ripple, the ideal response requires a discontinuity in the tangential field across the resonant surface, equivalent to a $735$6-function current sheet with surface current density $735$7. The nonlinear inner solution gives

$735$8

with

$735$9

The Grad–Shafranov solver and SPEC, which implements multi-region relaxed MHD, converge to the same singular ideal equilibrium. In SPEC, the error in the sheet-amplitude diagnostic scales approximately as $p_{Ml}^{m}=0$0, while uniform-refinement flux-surface geometry errors scale as $p_{Ml}^{m}=0$1. The paper also corrects a longstanding interpretation: the “sine qua non” condition is a linear regularity condition, not a necessary existence criterion for the resonant equilibrium (Huang et al., 2021).

Experimental and wave-theoretic studies show that ideal resonant response is not merely a formal boundary-layer construction. In ASDEX Upgrade, a rigidly rotating $p_{Ml}^{m}=0$2 magnetic perturbation produces edge displacements larger than vacuum-field predictions, and ECE-I measurements give a poloidal mode number $p_{Ml}^{m}=0$3 at $p_{Ml}^{m}=0$4, close to the synthetic VMEC value $p_{Ml}^{m}=0$5. The key ideal-MHD observation is that the magnetic perturbation spectrum peaks at non-resonant $p_{Ml}^{m}=0$6, while the displacement spectrum is almost resonant, $p_{Ml}^{m}=0$7, as expected from $p_{Ml}^{m}=0$8 (Willensdorfer et al., 2016). In coronal flux tubes, ideal and resistive analyses yield the same kink-wave frequency and damping rate in the small-resistivity limit, but radically different eigenfunctions: ideal eigenfunctions remain global, whereas resistive eigenfunctions localize around the resonance. The thin-boundary damping law,

$p_{Ml}^{m}=0$9

is recovered in the appropriate limit, and a separate invariant-imbedding analysis shows that resonant absorption in cylindrical geometry obeys a universal scaling in the variable $p_{El}^{m}=0$0 for linear, parabolic, and sinusoidal density profiles (Soler et al., 2013, Yu et al., 2019).

5. State selection, transport, and nonlinear resonant exchange

In nonlinear plasma stability, the ideal limit again does not imply a unique or especially benign resonant outcome. Nonlinear MHD simulations of high-$p_{El}^{m}=0$1 Wendelstein 7-X plasmas show that increasing parallel thermal conductivity reduces the linear growth rate of ideal ballooning modes, but the saturated pressure profile is barely affected. More importantly, broad and peaked pressure profiles behave differently: peaked profiles at lower $p_{El}^{m}=0$2 can degrade more severely than broad profiles at higher $p_{El}^{m}=0$3 and larger linear growth rate. A rotational-transform scan further shows that, when the linear growth rate is matched, resonant configurations containing the low-order $p_{El}^{m}=0$4 surface and non-resonant configurations relax the pressure profile by similar amounts. The paper’s explicit conclusion is that benign saturation is not guaranteed and is not dictated by linear growth or by the mere presence of a low-order resonance (Zhou et al., 1 Mar 2026).

In synthetic spin-flop bilayers, the contrast between ideal and asymmetric resonant behavior is especially sharp. An ideal synthetic antiferromagnet has two degenerate antiparallel ground states, identical coupling to a uniform hard-axis drive, and a single optical resonance for both states. Under resonant excitation it therefore enters a dynamic running state rather than a deterministically selectable one. Introducing thickness asymmetry $p_{El}^{m}=0$5 or bias asymmetry $p_{El}^{m}=0$6 creates a synthetic ferrimagnet, splits the optical resonance of the two antiparallel states, and permits state selection by frequency or field. In the reported devices, the optical resonance is near $p_{El}^{m}=0$7 GHz; $p_{El}^{m}=0$8 ns microwave pulses at the split frequencies toggle the state deterministically, with reported error rates below $p_{El}^{m}=0$9 (Koop et al., 2013).

Mesoscopic transport papers use the same motif in a different setting. In ideal armchair graphene nanoribbons of width $e^{-2\kappa L}$0 nm attached to leads twice as wide, the ribbon–lead interfaces act as electronic mirrors. In the noninteracting model, the conductance oscillations in magnetic field follow Fabry–Perot quantization,

$e^{-2\kappa L}$1

with an effective cavity length $e^{-2\kappa L}$2 nm. In the Hartree model, the pattern changes qualitatively to Aharonov–Bohm oscillations with period

$e^{-2\kappa L}$3

and the calculations give $e^{-2\kappa L}$4 T, implying $e^{-2\kappa L}$5 nm$e^{-2\kappa L}$6. Both models attribute the oscillations to resonant backscattering through quasi-bound states, but the spatial structure of those states differs sharply between the Fabry–Perot and Aharonov–Bohm regimes (Ihnatsenka, 2022).

Subgap resonant transport in normal–superconductor quantum-dot devices provides an even cleaner “ideal resonant” limit. In carbon-nanotube quantum dots with optimized Pd/Pb/In contacts, the superconducting transport gap is large and hard, $e^{-2\kappa L}$7–$e^{-2\kappa L}$8 meV, and the level broadening is much smaller than the gap. Thermally populated quasiparticles in the superconducting lead then generate two resonant lines,

$e^{-2\kappa L}$9

denoted TL and TR. The current is described by a simple resonant-tunneling integral with BCS quasiparticle density of states, and the zero-bias peak amplitude scales as

$\omega' \approx \omega - v k_x$0

Fits give $\omega' \approx \omega - v k_x$1 $\omega' \approx \omega - v k_x$2eV and $\omega' \approx \omega - v k_x$3 $\omega' \approx \omega - v k_x$4eV, and the transport gaps close monotonically with both temperature and magnetic field (Gramich et al., 2016).

Ideal impact constraints can also define a resonant manifold in strongly nonlinear mechanics. For two linearly coupled oscillators with perfectly elastic bilateral impacts, the action–angle reduction yields a resonant-manifold Hamiltonian whose generalized limiting phase trajectory must be optimized because the coupling energy depends on the actions. Matching the LPT to the saddle of the reduced Hamiltonian gives the critical delocalization threshold. In the vanishing-foundation limit, the asymptotic critical coupling is

$\omega' \approx \omega - v k_x$5

which agrees closely with the previously known $\omega' \approx \omega - v k_x$6 result (Perchikov et al., 2018).

6. Sensing, interaction control, and limiting factors

The graphene resonant pressure sensor realizes an explicitly device-level version of the resonant ideal. A trilayer graphene beam of length $\omega' \approx \omega - v k_x$7 $\omega' \approx \omega - v k_x$8m and width $\omega' \approx \omega - v k_x$9 $$\mathbf{p} = \frac{i}{\omega} \int \mathbf{J}(\mathbf{r})\, d^3 r,\qquad \mathbf{T} = \frac{1}{10 c} \int \big[\mathbf{r}(\mathbf{r}\cdot\mathbf{J}) - 2 r^2 \mathbf{J}\big]\, d^3 r.$$00m is fabricated on a circular silicon diaphragm of radius $$\mathbf{p} = \frac{i}{\omega} \int \mathbf{J}(\mathbf{r})\, d^3 r,\qquad \mathbf{T} = \frac{1}{10 c} \int \big[\mathbf{r}(\mathbf{r}\cdot\mathbf{J}) - 2 r^2 \mathbf{J}\big]\, d^3 r.$$01 mm and thickness $$\mathbf{p} = \frac{i}{\omega} \int \mathbf{J}(\mathbf{r})\, d^3 r,\qquad \mathbf{T} = \frac{1}{10 c} \int \big[\mathbf{r}(\mathbf{r}\cdot\mathbf{J}) - 2 r^2 \mathbf{J}\big]\, d^3 r.$$02 $$\mathbf{p} = \frac{i}{\omega} \int \mathbf{J}(\mathbf{r})\, d^3 r,\qquad \mathbf{T} = \frac{1}{10 c} \int \big[\mathbf{r}(\mathbf{r}\cdot\mathbf{J}) - 2 r^2 \mathbf{J}\big]\, d^3 r.$$03m. The front side is maintained at $$\mathbf{p} = \frac{i}{\omega} \int \mathbf{J}(\mathbf{r})\, d^3 r,\qquad \mathbf{T} = \frac{1}{10 c} \int \big[\mathbf{r}(\mathbf{r}\cdot\mathbf{J}) - 2 r^2 \mathbf{J}\big]\, d^3 r.$$04 kPa, the backside pressure $$\mathbf{p} = \frac{i}{\omega} \int \mathbf{J}(\mathbf{r})\, d^3 r,\qquad \mathbf{T} = \frac{1}{10 c} \int \big[\mathbf{r}(\mathbf{r}\cdot\mathbf{J}) - 2 r^2 \mathbf{J}\big]\, d^3 r.$$05 is varied from $$\mathbf{p} = \frac{i}{\omega} \int \mathbf{J}(\mathbf{r})\, d^3 r,\qquad \mathbf{T} = \frac{1}{10 c} \int \big[\mathbf{r}(\mathbf{r}\cdot\mathbf{J}) - 2 r^2 \mathbf{J}\big]\, d^3 r.$$06 to $$\mathbf{p} = \frac{i}{\omega} \int \mathbf{J}(\mathbf{r})\, d^3 r,\qquad \mathbf{T} = \frac{1}{10 c} \int \big[\mathbf{r}(\mathbf{r}\cdot\mathbf{J}) - 2 r^2 \mathbf{J}\big]\, d^3 r.$$07 kPa, and the diaphragm strain shifts the graphene resonance frequency. The chain rule is explicit:

$$\mathbf{p} = \frac{i}{\omega} \int \mathbf{J}(\mathbf{r})\, d^3 r,\qquad \mathbf{T} = \frac{1}{10 c} \int \big[\mathbf{r}(\mathbf{r}\cdot\mathbf{J}) - 2 r^2 \mathbf{J}\big]\, d^3 r.$$08

With $$\mathbf{p} = \frac{i}{\omega} \int \mathbf{J}(\mathbf{r})\, d^3 r,\qquad \mathbf{T} = \frac{1}{10 c} \int \big[\mathbf{r}(\mathbf{r}\cdot\mathbf{J}) - 2 r^2 \mathbf{J}\big]\, d^3 r.$$09 and $$\mathbf{p} = \frac{i}{\omega} \int \mathbf{J}(\mathbf{r})\, d^3 r,\qquad \mathbf{T} = \frac{1}{10 c} \int \big[\mathbf{r}(\mathbf{r}\cdot\mathbf{J}) - 2 r^2 \mathbf{J}\big]\, d^3 r.$$10 GHz per unit strain for one device, the measured responsivity is $$\mathbf{p} = \frac{i}{\omega} \int \mathbf{J}(\mathbf{r})\, d^3 r,\qquad \mathbf{T} = \frac{1}{10 c} \int \big[\mathbf{r}(\mathbf{r}\cdot\mathbf{J}) - 2 r^2 \mathbf{J}\big]\, d^3 r.$$11 Hz/Pa; a second device yields $$\mathbf{p} = \frac{i}{\omega} \int \mathbf{J}(\mathbf{r})\, d^3 r,\qquad \mathbf{T} = \frac{1}{10 c} \int \big[\mathbf{r}(\mathbf{r}\cdot\mathbf{J}) - 2 r^2 \mathbf{J}\big]\, d^3 r.$$12 Hz/Pa. The reported pressure resolution is $$\mathbf{p} = \frac{i}{\omega} \int \mathbf{J}(\mathbf{r})\, d^3 r,\qquad \mathbf{T} = \frac{1}{10 c} \int \big[\mathbf{r}(\mathbf{r}\cdot\mathbf{J}) - 2 r^2 \mathbf{J}\big]\, d^3 r.$$13 Pa, corresponding to $$\mathbf{p} = \frac{i}{\omega} \int \mathbf{J}(\mathbf{r})\, d^3 r,\qquad \mathbf{T} = \frac{1}{10 c} \int \big[\mathbf{r}(\mathbf{r}\cdot\mathbf{J}) - 2 r^2 \mathbf{J}\big]\, d^3 r.$$14 of full scale, and the performance is attributed to two design choices stated explicitly in the paper: a high-quality vacuum environment for the nanoresonator and stimulus delivery through the thin silicon diaphragm (More et al., 2023).

Optical Feshbach resonance in bosonic $$\mathbf{p} = \frac{i}{\omega} \int \mathbf{J}(\mathbf{r})\, d^3 r,\qquad \mathbf{T} = \frac{1}{10 c} \int \big[\mathbf{r}(\mathbf{r}\cdot\mathbf{J}) - 2 r^2 \mathbf{J}\big]\, d^3 r.$$15Sr provides the cold-atom analogue. The background scattering length is $$\mathbf{p} = \frac{i}{\omega} \int \mathbf{J}(\mathbf{r})\, d^3 r,\qquad \mathbf{T} = \frac{1}{10 c} \int \big[\mathbf{r}(\mathbf{r}\cdot\mathbf{J}) - 2 r^2 \mathbf{J}\big]\, d^3 r.$$16, so the gas is nearly ideal until light near a photoassociation resonance induces a complex scattering length $$\mathbf{p} = \frac{i}{\omega} \int \mathbf{J}(\mathbf{r})\, d^3 r,\qquad \mathbf{T} = \frac{1}{10 c} \int \big[\mathbf{r}(\mathbf{r}\cdot\mathbf{J}) - 2 r^2 \mathbf{J}\big]\, d^3 r.$$17. The optical coupling is summarized by the optical length

$$\mathbf{p} = \frac{i}{\omega} \int \mathbf{J}(\mathbf{r})\, d^3 r,\qquad \mathbf{T} = \frac{1}{10 c} \int \big[\mathbf{r}(\mathbf{r}\cdot\mathbf{J}) - 2 r^2 \mathbf{J}\big]\, d^3 r.$$18

and near an isolated resonance the inelastic rate coefficient is

$$\mathbf{p} = \frac{i}{\omega} \int \mathbf{J}(\mathbf{r})\, d^3 r,\qquad \mathbf{T} = \frac{1}{10 c} \int \big[\mathbf{r}(\mathbf{r}\cdot\mathbf{J}) - 2 r^2 \mathbf{J}\big]\, d^3 r.$$19

with $$\mathbf{p} = \frac{i}{\omega} \int \mathbf{J}(\mathbf{r})\, d^3 r,\qquad \mathbf{T} = \frac{1}{10 c} \int \big[\mathbf{r}(\mathbf{r}\cdot\mathbf{J}) - 2 r^2 \mathbf{J}\big]\, d^3 r.$$20 when $$\mathbf{p} = \frac{i}{\omega} \int \mathbf{J}(\mathbf{r})\, d^3 r,\qquad \mathbf{T} = \frac{1}{10 c} \int \big[\mathbf{r}(\mathbf{r}\cdot\mathbf{J}) - 2 r^2 \mathbf{J}\big]\, d^3 r.$$21 is negligible. Using the narrow $$\mathbf{p} = \frac{i}{\omega} \int \mathbf{J}(\mathbf{r})\, d^3 r,\qquad \mathbf{T} = \frac{1}{10 c} \int \big[\mathbf{r}(\mathbf{r}\cdot\mathbf{J}) - 2 r^2 \mathbf{J}\big]\, d^3 r.$$22 nm intercombination line with $$\mathbf{p} = \frac{i}{\omega} \int \mathbf{J}(\mathbf{r})\, d^3 r,\qquad \mathbf{T} = \frac{1}{10 c} \int \big[\mathbf{r}(\mathbf{r}\cdot\mathbf{J}) - 2 r^2 \mathbf{J}\big]\, d^3 r.$$23 kHz, the experiment finds that OFR can drive cross-dimensional thermalization in a gas that would otherwise remain nonthermalizing on experimental timescales. A combined analysis gives an effective decay rate $$\mathbf{p} = \frac{i}{\omega} \int \mathbf{J}(\mathbf{r})\, d^3 r,\qquad \mathbf{T} = \frac{1}{10 c} \int \big[\mathbf{r}(\mathbf{r}\cdot\mathbf{J}) - 2 r^2 \mathbf{J}\big]\, d^3 r.$$24 kHz, and elastic effects become visible around $$\mathbf{p} = \frac{i}{\omega} \int \mathbf{J}(\mathbf{r})\, d^3 r,\qquad \mathbf{T} = \frac{1}{10 c} \int \big[\mathbf{r}(\mathbf{r}\cdot\mathbf{J}) - 2 r^2 \mathbf{J}\big]\, d^3 r.$$25 for a $$\mathbf{p} = \frac{i}{\omega} \int \mathbf{J}(\mathbf{r})\, d^3 r,\qquad \mathbf{T} = \frac{1}{10 c} \int \big[\mathbf{r}(\mathbf{r}\cdot\mathbf{J}) - 2 r^2 \mathbf{J}\big]\, d^3 r.$$26 ms pulse (Blatt et al., 2011).

Across all of these realizations, ideality is productive but fragile. In optics, substrate-induced redshifts, absorption, and beam-purity requirements limit perfect anapole suppression (Lu et al., 2021). In moving cavities, nonlocality and nonlinearities curtail the formal singularity (Guo et al., 2013). In resonant photonic crystals, positional disorder narrows the gap and enhances edge fluctuations (Skipetrov, 2019). In MHD, finite resistivity, reconnection, and global kink structure perturb ideal shielding (Pharr et al., 18 Mar 2026, Huang et al., 2021). In resonant devices, interface quality, lifetime broadening, and packaging determine whether the underlying resonance remains spectroscopically clean (Gramich et al., 2016, More et al., 2023). The term “resonant ideal” therefore marks not a perfect end state, but a technically useful asymptotic regime in which resonance exposes the clearest possible structure of a system before non-ideal corrections dominate.