Higher-Order Mie Resonances in Dielectrics

• Higher-order Mie resonances are defined as electromagnetic modes in dielectric spheres where multipolar orders beyond dipoles, such as quadrupoles and octupoles, produce enhanced internal field buildup and narrow linewidths.
• Analytical methods, including Padé approximants, reveal that resonance conditions exhibit order-dependent scaling laws—with quality factors rising sharply with refractive index and multipolar order.
• These resonances impact practical applications by enabling advanced functionalities in nonlinear optics, high-resolution microsphere imaging, and collective photonic metastructures.

Higher-order Mie resonances are electromagnetic resonances of dielectric or related spherical resonators whose multipolar order exceeds the dipole and whose fields are governed by the partial-wave structure of Mie theory. In the cited literature, the term covers at least two closely related regimes: resonances with larger multipolar order such as quadrupoles and octupoles in subwavelength high-index particles, and ultra-sharp high-order internal partial-wave modes in microspheres, described as optical super-resonances, that occur only at exceptionally specific size parameters and refractive indices (Zambrana-Puyalto et al., 2024, Wang et al., 2022). Across these regimes, higher-order resonances are associated with more angular nodes, stronger spatial structure, narrower linewidths, higher quality factors, enhanced internal field buildup, and a pronounced sensitivity to geometry, material contrast, and spectral matching.

1. Definitions and multipolar taxonomy

In standard Lorenz–Mie theory, the electromagnetic field of a sphere is decomposed into electric and magnetic multipolar channels. The external scattered field is described by the electric and magnetic Mie coefficients $a_n$ and $b_n$, while the internal field is described by $c_n$ and $d_n$ (Tzarouchis et al., 2017, Tribelsky et al., 2015). The multipolar order is indexed as dipole, quadrupole, octupole, and so on; one paper states this explicitly as $j=1$ for dipole, $j=2$ for quadrupole, and $j=3$ for octupole (Zambrana-Puyalto et al., 2024).

For subwavelength dielectric nanoparticles, higher-order Mie resonances are often identified with quadrupolar channels. In a coupled-multipole treatment of resonant assemblies, the dominant optical response is written in terms of electric dipole (ED), magnetic dipole (MD), electric quadrupole (EQ), and magnetic quadrupole (MQ) contributions, with still higher multipoles neglected (Ustimenko et al., 2021). In this usage, “higher-order” mainly means the quadrupolar resonances above the usual dipole response.

A broader usage appears for dielectric microspheres. There, the sphere is illuminated by a plane wave and the internal field coefficients become resonant when their denominators become very small. The resulting modes are high-order internal Mie resonances, described as optical super-resonances, because they are resonant internal partial-wave modes that can generate field enhancements far beyond the photonic nanojet regime (Wang et al., 2022). This terminological spread suggests that higher-order Mie resonances should be understood not as a single narrowly delimited family, but as a hierarchy of multipolar and internal-cavity resonances whose detailed manifestation depends on scale and excitation conditions.

2. Resonance conditions, linewidths, and quality factors

The central geometric control parameter is the size parameter, written either as $x=ka$ or $q=2\pi a/\lambda$, depending on notation (Tzarouchis et al., 2017, Wang et al., 2022). Resonance occurs when the denominator of the relevant Mie coefficient approaches zero, so the response acquires a pole-like structure. A recurrent theme in the analytical literature is that simple Taylor expansions obscure this structure, whereas Padé approximants preserve it and separate static, dynamic, and radiative contributions more transparently (Tzarouchis et al., 2016, Tzarouchis et al., 2017).

For small dielectric spheres, Padé-based pole formulas make the size dependence of higher-order resonances explicit. For magnetic resonances up to fifth order, one paper summarizes the pole locations as

$$\varepsilon_{b_n}= -\frac{2}{2n-1}+\left(\frac{p_n}{x}\right)^2 -i\frac{2}{\left[\left(2n-1\right)!!\right]^2}x^{2n-1}\left(1-t_nx^2\right),$$

while for dielectric electric resonances up to fourth order it gives

$b_n$0

The same analysis emphasizes that the imaginary part of the pole condition is the radiative damping term and directly controls the linewidth; higher-order resonances are progressively narrower and harder to observe (Tzarouchis et al., 2017).

A complementary asymptotic description is given for high-index dielectric spherical resonators. For the lowest-energy resonance of each multipole family in a lossless sphere,

$b_n$1

Thus the exponent of $b_n$2 increases by $b_n$3 with each step in multipolar order, and electric modes carry a higher power of $b_n$4 than magnetic modes of the same $b_n$5 (Zambrana-Puyalto et al., 2024). The first few asymptotic examples are

$b_n$6

$b_n$7

$b_n$8

The paper further notes that $b_n$9 is achievable for octupolar modes at $c_n$0, and that the high-$c_n$1 asymptotic scaling becomes accurate at lower refractive index as the multipolar order increases (Zambrana-Puyalto et al., 2024).

3. Internal-field localization and optical super-resonances

One of the most extreme realizations of higher-order Mie physics is the optical super-resonance of dielectric microspheres. The internal mode strength can be quantified by an internal scattering efficiency

$c_n$2

with

$c_n$3

These resonances are extremely sensitive to $c_n$4: coarse sampling with $c_n$5 shows only modest peaks, whereas refinement to $c_n$6–$c_n$7 reveals many narrow resonances and dramatically larger internal efficiency (Wang et al., 2022).

The field distribution is not uniform throughout the particle. Instead, it localizes into highly concentrated regions, often as two nearly symmetric hotspots along the propagation axis. For $c_n$8, the most intense electric-field resonance occurs near $c_n$9, with $d_n$0, and the strongest magnetic-field resonance occurs near $d_n$1, with $d_n$2 (Wang et al., 2022). As the refractive index increases, the reported peak values rise to $d_n$3 and $d_n$4 for $d_n$5, $d_n$6 and $d_n$7 for $d_n$8, and $d_n$9 and $j=1$0 for $j=1$1 (Wang et al., 2022).

A related analytical perspective comes from high-refractive-index spheres. There the outer scattering problem and the inner problem behave very differently as $j=1$2 increases. Outside, each partial scattered wave can be decomposed into a background term corresponding to a perfectly reflecting sphere and a resonant term, yielding asymmetric Fano profiles. Inside, by contrast, electric and magnetic Mie resonances of different orders overlap substantially and can produce giant in-particle field concentration (Tribelsky et al., 2015). The asymptotic electric and magnetic resonance ladders,

$j=1$3

show why different orders can crowd together in the Fraunhofer regime and promote overlap (Tribelsky et al., 2015).

The microsphere literature connects these super-resonances directly to focusing and nanoscopy. Conventional photonic nanojets are often near or only modestly below the diffraction limit, whereas super-resonant internal Mie modes are proposed as the missing mechanism behind the $j=1$4 super-resolution reported in microsphere-assisted imaging (Wang et al., 2022).

4. Collective, hybridized, and finite-size higher-order resonances

Higher-order Mie resonances are not restricted to isolated particles. In multiple-scattering problems, higher multipoles modify both the physical response and the convergence properties of reduced models. A coupled-multipole/Born-series formulation that retains ED, MD, EQ, and MQ shows that the critical separation needed for convergence is larger for quadrupole resonances than for dipole resonances. For nonabsorbing particles, the quoted thresholds are approximately

$j=1$5

with comparable but polarization-dependent values for MD and MQ. The same work notes that even the third-order Born approximation can keep the scattering-error below about $j=1$6 when the interaction parameter is small (Ustimenko et al., 2021).

In finite metastructures, the collective problem is not exhausted by Bloch-wave BIC terminology. One paper states that, unlike a common belief, the bound states in the continuum derived by the Bloch-wave theory do not directly determine the resonance with the highest $j=1$7 value in large but finite arrays. Higher $j=1$8 factors are associated instead with collective resonances formed by nominally guided modes below the light line, with strong effect of both electric and magnetic multipoles (Hoang et al., 2024). In its 1D example, the intrinsic MQ mode splits into bonding and antibonding collective states, with $j=1$9 and exponents $j=2$0 for MQ-B, $j=2$1 for MQ-A, and $j=2$2 for MO (Hoang et al., 2024).

Hybridization can also involve non-dielectric channels. In a nonreciprocal Tellegen sphere, the axion coupling $j=2$3 mixes electric and magnetic multipoles so that an electric source radiates a magnetic multipole outside the sphere and vice versa. The same coupling hybridizes the Mie resonances and produces characteristic double-peak structures, especially clearly for higher-order multipoles, whose resonances are narrower and therefore more sensitive to hybridization (Seidov et al., 6 Feb 2025). In hybrid WS$j=2$4-on-gold nanoantennas, a higher-order anapole mode (HOAM) couples strongly to a Fabry–Pérot-plasmonic mode, yielding a supercavity mode with experimental $j=2$5 and anti-crossing splitting of $j=2$6 meV (Randerson et al., 2023).

5. Functional consequences: magneto-optics, nonlinearity, topology, and temporal dynamics

Higher-order Mie resonances qualitatively reshape optomagnetic excitation landscapes. In a 400 nm diameter Bi-substituted iron-garnet sphere with $j=2$7, the inverse Faraday effect field

$j=2$8

becomes highly structured at resonance. Lower-order MD is the most extended mode, with $j=2$9 for $j=3$0, whereas higher-order modes are more confined, more nonuniform, more sign-changing, and richer in hotspots and nodes (Krichevsky et al., 2024). The electric resonance family shows a direct hotspot count progression: ED gives one $j=3$1 hotspot, EQ gives two, and EO gives four. The in-plane effective field forms a vortex at every Mie resonance and is described as a Neel-type optical skyrmion-like field; for higher-order modes, especially MQ and EQ, it develops a double-ring pattern with inner and outer regions of opposite helicity (Krichevsky et al., 2024).

In nonlinear nanocomposites, higher-order resonances enhance and even invert the effective Kerr response. For GaP spheres at $j=3$2 nm, the magnetic quadrupole and electric quadrupole resonances are explicitly identified near $j=3$3 nm and $j=3$4 nm, complementing the magnetic and electric dipoles near $j=3$5 nm and $j=3$6 nm. Near resonance, the effective nonlinear index reaches magnitudes on the order of $j=3$7 to $j=3$8, compared with the bulk GaP estimate $j=3$9, and the sign of the effective optical Kerr coefficient is inverted near the Mie resonances (Panov, 2018).

Higher-order resonances also support integrated photonic functionality. A TiO$x=ka$0-like dielectric building-block assembly uses a collective magnetic octupole resonance near $x=ka$1 nm for focusing and collective magnetic/electric dipole modes near $x=ka$2 nm for waveguiding. The reported electric-field enhancement at the focus is about $x=ka$3, the intensity enhancement is $x=ka$4, and the transfer efficiency along a 40-sphere chain reaches about $x=ka$5 at $x=ka$6 nm (Chattaraj et al., 2016). In confined Mie resonance photonic crystals, embedding PECs between dielectric rods suppresses the $x=ka$7 leakage of ordinary dielectric Mie states, yielding disentangled higher-orbital bands, complete band gaps, and third-order topology with bulk, surface, hinge, and corner states at $x=ka$8 GHz, $x=ka$9 GHz, $q=2\pi a/\lambda$0 GHz, and $q=2\pi a/\lambda$1 GHz, respectively (Li et al., 2023).

Temporal dynamics can also be decisive. In laser-driven plasma nanoshells, high-order Mie resonances produce about threefold electric-field enhancement at $q=2\pi a/\lambda$2 nm in optimized geometries, but the buildup time is roughly $q=2\pi a/\lambda$3 fs and the enhanced state persists for about $q=2\pi a/\lambda$4 fs before plasma expansion detunes the resonance. A 4-cycle pulse gives only about $q=2\pi a/\lambda$5 enhancement, whereas full enhancement typically needs tens of cycles (Gao, 30 Oct 2025). This suggests that, in driven nanoplasma systems, resonance establishment time is itself a design parameter.

6. Computation, approximations, and conceptual caveats

The numerical treatment of higher-order Mie resonances is unusually delicate because their spectral features can be exceptionally narrow and their field structure highly oscillatory. The microsphere super-resonance calculations explicitly show that sampling accuracy up to $q=2\pi a/\lambda$6 is needed to uncover resonance peaks hidden at coarser resolution (Wang et al., 2022). For general axially symmetric dielectric bodies, a Fourier–Nyström solver based on combined integral equations computes complex eigenwavenumbers and eigenfields with high accuracy even at very high wavenumbers, making it suitable for benchmarking high-$q=2\pi a/\lambda$7, high-order whispering-gallery and Mie-type resonances (Helsing et al., 2016).

Reduced-order models also have clear limits. A study of non-spherical Mie-resonant dielectric disks shows that induced dipole moments are defined not only by the field at the particle center but also by second-order spatial derivatives of the field. This intrinsic nonlocality is especially pronounced in the vicinity of the anapole minimum in the scattering cross-section and can reach up to $q=2\pi a/\lambda$8 of the local response (Bobylev et al., 2020). Accordingly, higher-order Mie physics can enter even when the observable response appears dipole-dominated.

Several recurrent misconceptions are corrected by the recent literature. Higher-order Mie resonances are not merely stronger versions of dipoles; they are a systematically more $q=2\pi a/\lambda$9-enhanced family with order-dependent scaling laws (Zambrana-Puyalto et al., 2024). They are not only external scattering peaks; the internal problem can display giant resonant buildup even when the outside field tends toward the perfectly reflecting-sphere limit (Tribelsky et al., 2015). Nor are the sharpest finite-array resonances necessarily the direct finite-size descendants of Bloch-wave BICs (Hoang et al., 2024). Finally, the phrase “higher-order” itself is context dependent: in one setting it denotes EQ and MQ channels above ED and MD, while in another it denotes extremely high-order internal partial waves in microspheres (Ustimenko et al., 2021, Wang et al., 2022).

Taken together, these results establish higher-order Mie resonances as a broad resonant hierarchy rather than a single phenomenon. They govern linewidth narrowing, giant internal field concentration, structured optomagnetic forcing, nonlinear-response enhancement, hybrid supercavity formation, tight-binding-like higher-orbital photonics, and transient nanoplasma dynamics. Their common foundation is the same: multipolar and internal-cavity solutions of Maxwell’s equations whose observability and functionality depend critically on size parameter, refractive index, radiative damping, loss, geometry, and collective coupling.