Contents

Mie CoefficientsEdit

Mie Coefficients refer to a pair of sets of complex numbers that entirely determine how a plane electromagnetic wave is scattered by a homogeneous sphere. Named after the German physicist Gustav Mie, who derived the solution in 1908, these coefficients—traditionally denoted a_n and b_n for each multipole order n—represent the electric and magnetic multipole contributions to the scattered field. Their calculation rests on the exact solution of Maxwell’s equations in a spherical geometry and relies on a careful handling of special functions, boundary conditions, and a convergent partial-wave expansion. From the a_n and b_n, one can obtain the angular distribution of scattered light, its polarization, and the various cross sections that quantify extinction, scattering, and backscattering.

The mathematical structure of the problem is an archetype in wave physics: a plane wave impinges on a sphere of radius a with a relative refractive index m (the index inside the sphere relative to the surrounding medium). The incident and scattered fields are expanded in vector spherical wave functions, and the boundary conditions at the sphere’s surface yield closed-form expressions for the coefficients a_n and b_n. The theory is equally valid for absorbing spheres and can accommodate complex refractive indices, which is essential for modeling realistic materials. In practical terms, the Mie coefficients provide a bridge from microscopic material properties to macroscopic optical observables.

Mathematical formulation

  • Size parameter and indices

    • The key dimensionless parameters are the size parameter x = k a, where k is the wavenumber of light in the surrounding medium, and the relative refractive index m = n_sphere / n_medium. The physics of scattering changes character as x and m vary, from Rayleigh-like behavior for small particles to complex resonant patterns for larger ones.
    • The fields are expressed in terms of Riccati-Bessel functions psi_n and xi_n, which are built from spherical Bessel and spherical Hankel functions. See Riccati-Bessel functions and spherical Bessel function for the mathematical background.

  • Coefficients a_n and b_n

    • For each integer n ≥ 1, the electric-type (a_n) and magnetic-type (b_n) coefficients are given by ratios of combinations of psi_n(mx), psi_n'(mx), psi_n(x), psi_n'(x), xi_n(x), and xi_n'(x). The derivatives are with respect to the argument.
    • A compact way to state them is:
    • a_n = [m psi_n(mx) psi_n'(x) − psi_n'(mx) psi_n(x)] / [m psi_n(mx) xi_n'(x) − psi_n'(mx) xi_n(x)]
    • b_n = [psi_n(mx) psi_n'(x) − m psi_n'(mx) psi_n(x)] / [psi_n(mx) xi_n'(x) − m psi_n'(mx) xi_n(x)]
    • Here psi_n(z) = z j_n(z) and xi_n(z) = z h_n^(1)(z), with j_n the spherical Bessel function and h_n^(1) the spherical Hankel function of the first kind. The primes denote derivatives with respect to the argument.
    • These coefficients fully determine the scattered field for a homogeneous sphere and underpin all derived observables. See vector spherical harmonics for the broader basis in which the fields are expanded.

  • Cross sections and angular observables

    • Extinction cross section (Q_ext) and scattering cross section (Q_sca) are obtained by summing over the multipole contributions of a_n and b_n:
    • Q_ext = (2/x^2) ∑_{n=1}^∞ (2n+1) Re(a_n + b_n)
    • Q_sca = (2/x^2) ∑_{n=1}^∞ (2n+1) (|a_n|^2 + |b_n|^2)
    • The backscattering cross section (Q_back) depends on the difference a_n − b_n and can be written as:
    • Q_back = (1/x^2) |∑_{n=1}^∞ (2n+1)(−1)^n (a_n − b_n)|^2
    • The differential scattering pattern, including polarization, can be expressed in terms of the same coefficients through standard Mie-scattering amplitudes S_1(θ) and S_2(θ), which are built from sums of a_n and b_n weighted by angular functions.

  • Special cases and extensions

    • Rayleigh limit (x ≪ 1): In this regime, only the lowest-order multipoles (primarily n = 1) contribute significantly. The scattering behavior mirrors the classical dipole response of a small particle, and the cross sections scale with x to the appropriate power (e.g., Q_sca ∝ x^4 in the simplest non-absorbing case). The particle’s polarizability α ≈ 4π a^3 (m^2 − 1)/(m^2 + 2) captures the leading order, linking Mie theory to simpler dipole models.
    • Geometric optics limit (x ≫ 1): For very large spheres relative to the wavelength, the scattering behavior becomes highly intricate with many contributing multipoles. The extinction cross section approaches the geometric-optics expectations, but a classic result known as the extinction paradox reveals that Q_ext tends toward 2 in certain regimes, effectively doubling the intuitive geometric cross section due to interference and shadowing effects.
    • Core-shell and layered spheres: Mie theory extends to coated spheres by applying the same boundary-condition machinery at multiple interfaces. The resulting coefficients incorporate the core and shell parameters (radii, refractive indices) and reduce to the homogeneous-sphere formulas when the shell vanishes. See core-shell particle for related concepts.
    • Non-spherical particles and alternatives: While Mie coefficients are exact for spheres, many real-world particles are non-spherical. In those cases, methods such as the T-matrix method or approximations based on spherical harmonics expansions around an effective shape are used, with trade-offs in accuracy and computational cost. For highly irregular particles, numerical techniques like finite-difference time-domain (FDTD) or discrete dipole approximation (DDA) may be employed.

Applications and impact

  • Atmospheric and environmental science: The Mie coefficients underpin models of how aerosols scatter sunlight, affecting radiative forcing, visibility, and climate predictions. Accurate retrievals of aerosol size distributions and compositions rely on robust Mie-based scattering calculations. See atmospheric scattering and aerosols.
  • Remote sensing and astronomy: Scattering by dust and droplets in planetary atmospheres and interstellar environments is analyzed using Mie theory, enabling interpretation of spectral and angular signatures. See astronomy and remote sensing.
  • Biomedical optics and nanophotonics: Metallic and dielectric nanoparticles exhibit resonant scattering governed by Mie-type coefficients, giving rise to tunable colorimetric responses, imaging contrast mechanisms, and engineered light-matter interactions at the nanoscale. See plasmonics and nanoparticles.
  • Materials design: Designing particles with tailored scattering properties—such as preferred forward scattering, minimal backward scattering, or specific color effects—draws directly on the behavior of a_n and b_n across wavelengths and size ranges. See optical engineering and materials science.

Computational aspects

  • Convergence and numerical stability: The infinite series in a_n and b_n must be truncated in practice. Convergence is rapid for small x but requires more terms as x grows. Efficient computation relies on stable recursion relations for the special functions involved (psi_n, xi_n) and careful handling near resonances.
  • Public codes and resources: Implementations of Mie theory range from dedicated packages to general-purpose physics libraries. Researchers often consult standard references and test their results against benchmark cases such as the Rayleigh limit or known coated-sphere solutions. See Mie theory for foundational material and related numerical methods.

See also