Mie scattering 2026-10-06
The electromagnetic scattering solution for a sphere describes the transition between Rayleigh scattering and geometrical optics. Grain-size-dependent efficiencies can produce structure in the radiation-pressure coefficient.
Choose a vibrational Raman shift magnitude , excitation vacuum wavelength and temperature . These are original illustrative inputs, not a measured 1928 optical spectrum or a claim about a particular liquid. They give
Take reciprocal lengths in consistent units before converting to nanometres. At , and . The predicted line-integrated radiant intensity ratio is , whereas the ideal photon count-rate ratio is . Thus the radiometric ratio is about percent, rather than the percent obtained by using the Boltzmann factor alone.
Figure 1.
Molecular energy transfer and an original one-mode Raman scattering line model
. Upper panels show initial and final vibrational levels with schematic virtual intermediates. Lower panel uses the stated , , example. Chosen equal Gaussian function linewidths have standard deviation ; the Rayleigh scattering line has an arbitrarily chosen integrated area of relative to a unit Stokes area. The anti-Stokes area follows the line-integrated radiant intensity formula. The vertical axis is logarithmic. No line shape, width, relative elastic strength or line position is a historical measurement. Source. Details.
The gross requirement is an anisotropic polarizability tensor, so that changing molecular orientation changes the induced response. A permanent electric dipole moment is not required; a homonuclear diatomic can therefore have a rotational Raman scattering spectrum.
For the inelastic lines of a linear rigid rotor, the detailed rule is , with . The contribution is unshifted Rayleigh scattering, not a rotational Stokes Raman scattering line. Stokes Raman scattering has because the molecule gains rotational energy.
The size parameter , complex refractive index and particle shape determine scattering strength and direction. In the small-sphere Rayleigh scattering regime, scattering cross-section scales as ; the Mie scattering regime is needed when particle size approaches the wavelength. Particle number density and particle size distribution, aerosol hygroscopicity and the angular fraction returned to space control an ensemble's radiative effect. High scattering does not imply strong cooling if most light is scattered forward or the particle also absorbs strongly; single-scattering albedo distinguishes scattering from absorption.
Cloud particles can scatter or absorb, so not every atmospheric haze makes a high-albedo planet. Weak molecular bands alone also admit low abundance or high mean molecular weight; a consistent spectral combination is stronger evidence than a single signature.
Three observable signatures are a blue optical transmission slope, muted molecular transmission features, and enhanced reflected optical light.
Small-particle atmospheric haze can increase the apparent transit radius toward short wavelengths, roughly --. For Rayleigh scattering-like extinction, and an isothermal atmospheric scale height model gives . Larger particles instead produce a flatter, approximately grey continuum.
An optically thick exoplanet cloud deck can truncate the slant path through gas and flatten molecular bands in a exoplanet transmission spectrum, for example over -- and the water band around . The observed quantity is reduced band-to-continuum radius contrast, not necessarily a low molecular abundance.
exoplanet cloud scattering can produce a high geometric albedo and distinctive optical secondary eclipses/phase curves around --. Longitudinal exoplanet cloud variation can shift a reflected-light phase maximum. A spectral slope or flat spectrum alone is not a unique exoplanet cloud identification: particle size, gas scattering, metallicity and stellar contamination must be considered together.
For a spherical grain of density , the radiation-pressure coefficient is
The dependence of radiation pressure cancels that of stellar gravity. Here radiation-pressure efficiency includes absorption and the momentum-transfer part of scattering, averaged over the stellar spectrum. In the geometrical optics regime, exceeds the important stellar wavelengths, is of order unity, and : cross-sectional area grows as , whereas mass grows as . Around the stellar wavelength, Mie scattering can produce a broad maximum and material-dependent structure. In the Rayleigh scattering regime, absorption gives , while scattering gives . Thus small absorbing grains approach an approximately constant , whereas nearly transparent grains have . A universal fall to zero at small is therefore inappropriate without specifying the optical properties. The sketch shows both possible small-grain limits; its vertical normalization is illustrative.
Figure 1.
Radiation-to-gravity ratios for absorbing and scattering grains, and the signed angular lag over five dust orbits for beta equal to 0.01 and 0.04
. Details.
Assume release with negligible velocity relative to the planet, negligible planetary gravity after release, and a constant radiation-pressure coefficient. Put and . The inherited speed and specific angular momentum are and . The new specific orbital energy is
For , the dust orbit released from a circular parent ring is an ellipse. Using and gives
For , release is at periapsis, since . At the radiation-pressure blowout threshold, , the orbit is parabolic. For it is hyperbolic: the signed semimajor axis is negative and . For the central force ceases to be attractive, so the elliptical interpretation and the positive-eccentricity formula cannot be extended unchanged.
With measured from release, it is the true anomaly. The polar equation of a Kepler orbit and conservation of specific angular momentum give
Define the signed dust-tail angular lag by , with and both angles unwrapped. For fixed , expansion to first order in gives
so . The dust trails the planet, hence the signed lag is negative. Its slope is ; it decreases monotonically, has horizontal tangents at successive release-direction passages, and has a sinusoidal ripple about . At the end of orbit , and . At five orbits this is for and for . The upper- sketch remains a first-order approximation; its accumulated error is .
For the subsequent positive-distance statements, set . The mean motion of the grain is
A full relative wrap, , therefore takes after averaging over the orbital ripple. During this time Poynting–Robertson drag changes the semimajor axis by
The signed change is negative. This estimate treats as fixed in the rate and computes the drag accumulated over the radiation-pressure wrap time. Since , drag causes little migration over that time. For an arbitrarily tiny , however, this small migration can itself alter the relative phase appreciably: neglecting that feedback on the wrap time additionally requires .
Rapid dust sublimation is a plausible sink near a hot planet: a grain can evaporate long before it completes a relative wrap. Destructive collisions can also remove visible grains. Poynting–Robertson drag alone, with the estimate above, does not explain immediate removal from a short tail. If the destruction time is , a short tail requires roughly .
For a steady dust-tail continuity equation, assume a constant injection rate , grains with the same fixed , negligible release-speed dispersion, and negligible destruction within the segment being calculated. A finite destruction lifetime can terminate that segment, or multiply its density by a survival probability. To first order,
To obtain the explicit trigonometric profile, make the additional secular phase approximation for a dust tail, , discarding the bounded term in the angle-to-age map while retaining the periodic instantaneous drift speed. Then
Since the square is even, the same written profile applies to the signed if number intervals are interpreted with positive orientation. The result is an approximate phase-substitution profile, not a uniformly valid consequence of . The consistent first-order description is instead parametric:
In particular, immediately after release, and , whereas the phase-substitution expression would give . Nor does by itself justify neglecting the orbital ripple. The formal enhancements where are dust-tail caustics; a spread of radiation-pressure coefficients, finite release velocities, and destruction smooth them. With constant lifetime , the parametric profile acquires . These qualifications state precisely the extra assumptions behind the explicit profile.
One approach uses exoplanet transmission spectra. Three aerosol signatures are:
A second approach uses the dayside spectrum and brightness obtained with exoplanet secondary eclipses, optionally supported by orbital phase variations. Its three useful aerosol signatures are:
Transmission measures aerosol extinction along the limb; eclipse/reflection measurements test dayside scattering and emission. No single flattened feature proves clouds: small atmospheric scale height, low molecular abundances, stellar heterogeneity and thermal-gradient degeneracies must be considered. Joint wavelength coverage and geometries make the inference stronger.
Assume that the optical continuum is dominated by Rayleigh scattering by particles much smaller than the wavelength, and that their column abundance is fixed. Their cross-section, hence the slant optical depth, scales as . Differentiating the annulus model for transmission spectroscopy exactly gives the logarithmic slope of total transit depth
For and , the slope lies between and zero, but generally depends on wavelength. In the optically thin limit, , so its slope is . It vanishes as the atmospheric signal becomes negligible compared with the opaque disc; it also approaches zero in the saturated optically thick limit.
The familiar value applies to the optically thin atmospheric excess,
or formally to a model in which is negligible. For an ordinary thin exoplanet atmosphere, however, , so neglecting the opaque planetary baseline in the total transit depth is not usually justified. A universal for the logarithm of the total transit depth does not follow from the stated model.
In a stratified isothermal atmosphere, the related scattering slope of a transmission spectrum is instead , obtained by moving the effective tangent level to keep its slant optical depth of order unity. Consequently there. The radius slope, atmospheric-excess slope and total-depth slope are different observables.
In the annulus model for transmission spectroscopy, an added aerosol optical depth changes the signal to . Large particles can supply nearly wavelength-independent extinction. An opaque high exoplanet cloud deck then masks deeper gas, flattens the optical exoplanet transmission spectrum, and weakens atomic or molecular features.
Small particles can instead produce a rising transit radius toward short wavelengths. In an isothermal atmosphere, the slant optical depth of an isothermal atmosphere is
Taking the effective radius near gives . For Rayleigh scattering, , so the scattering slope of a transmission spectrum is
A steeper-than-expected optical slope or suppressed gas features can therefore indicate atmospheric haze or exoplanet clouds. These signatures are not unique: high mean molecular mass reduces , gas itself can produce Rayleigh scattering, and stellar surface heterogeneity can mimic slopes. Consistent optical and infrared features help distinguish these explanations.
For a spherical grain of diameter , density , and radiation-pressure efficiency , comparison of the stellar momentum flux with Newton's law of universal gravitation gives the radiation-pressure coefficient
For much larger than the characteristic stellar wavelength, , so . When becomes comparable to the optical wavelength, diffraction and the grain's composition make size dependent and reaches a broad maximum. Deep in the Rayleigh scattering regime, absorption can give and hence nearly constant , while scattering alone gives and hence . A realistic curve therefore rises roughly as toward micron sizes, turns over near the stellar spectral peak, and falls or flattens for still smaller grains.
Sunlight is approximately unpolarized before its interaction with the environment. For Rayleigh scattering by small isotropic scatterers, the induced electric dipoles radiate only their components transverse to the viewing direction. For unpolarized incident light scattered through angle , the perpendicular and parallel intensities are proportional to 1 and . The degree of linear polarization is therefore
At single scattering is fully linearly polarized perpendicular to the scattering plane. This explains the polarized sky away from the Sun; multiple scattering and aerosols reduce the ideal contrast.
Reflection at a dielectric surface treats the two linear polarizations differently. At Brewster incidence, , the reflected p component vanishes for ideal transparent nonmagnetic media. The remaining reflected glare is s polarized, perpendicular to the incidence plane, approximately horizontal for a horizontal water surface. Away from this angle the reflection is generally partially polarized. The transmitted beam also has unequal s/p transmission.
Figure 1.
Original diagrams of polarization perpendicular to a single-scattering plane and of s-polarized reflection near Brewster incidence at a dielectric surface
. Details.
Refraction and internal reflection in water droplets likewise have polarization-dependent Fresnel coefficients, producing strong linear polarization in rainbow light. These are generic scattering and boundary-condition mechanisms; no biological detection mechanism is required.
Both Thomson scattering and Rayleigh scattering can be understood as radiation from a dipole driven by the incident electric field. Their different frequency dependences come from the response of a free charge compared with a bound charge.
For a free particle of charge and mass , a harmonic field gives , hence . The radiated field contains two time derivatives of , cancelling the response's . The Thomson cross section is consequently independent of frequency in the classical low-recoil regime:
For incident linear polarization with unit polarization , summing over outgoing polarizations gives . Radiation vanishes along the acceleration and is strongest perpendicular to it. Averaging an incident unpolarized light beam gives
where is the scattering angle. Scattering through is fully linearly polarized perpendicular to the scattering plane for this ideal isotropic dipole model. The electron case dominates over equally charged heavier particles because the cross section scales as . recoil corrections become necessary when photon energy is no longer small compared with .
For a bound charge with natural frequency and damping rate ,
Far below resonance, is approximately constant, so
This is Rayleigh scattering. It requires a scatterer small compared with the wavelength and, for the simple expression, an isotropic linear dipole response. Atoms, molecules and small dielectric particles provide examples. The stronger scattering of shorter visible wavelengths explains blue skylight; absorption and the nonuniform solar spectrum prevent treating that simple scaling as a complete color calculation.[ref]
For isotropic electric scatterers the normalized angular and polarization patterns are the same as in Thomson scattering. The frequency scaling distinguishes their limits: a free charge has a displacement proportional to , whereas a tightly bound charge has an approximately frequency-independent displacement. Near resonance, far outside the small-particle limit, or with anisotropic polarizability, the simple Rayleigh formula must be modified. In particular, the electric-only pattern must not silently be applied to a perfectly conducting sphere whose magnetic dipole response is also retained.
A scattering measurement should subtract stray-light and detector backgrounds, control scatterer density and optical thickness, and calibrate detector response and incident power when comparing wavelengths. A polarization analyzer at tests the dipole pattern; multiple scattering or unresolved nonisotropic contributions can reduce its polarization.
Spherical-particle frequency scaling. For a dielectric sphere of radius , the quasistatic electric polarizability is
The electric dipole radiation formula and give
Thus Rayleigh scattering scales as , or , when the polarizability varies slowly with frequency. Material resonances and absorption can change this simple scaling. For the incident polarization, the dipole is along : intensity is proportional to , where is the angle from . There is no radiation on the dipole axis and a maximum in the plane. In a direction the scattered electric field is along , hence linearly polarized wherever nonzero.[ref]
Figure 1.
A z-directed electric dipole radiates most strongly perpendicular to its axis and has axial zeros
. Source. Details.
Needle model used in parts (a)–(d). Let its axis be . A very slender, high-permittivity needle has dominant longitudinal response, so use the rank-one needle polarizability . Its transverse response is neglected at leading order. For observation along , define , and
For a unit incident Jones vector , the common far-field relation is , and . Here is irradiance at distance ; multiplying by gives power per solid angle. No absolute value of is given or needed for the requested relative intensities.
Diffuse random cloud, with the incident polarization of (a). The single-particle amplitudes through and analysers are proportional to and . The isotropic orientation fourth moment gives the needle-scattering analyzer moments
For example, ; symmetry and give . Since the particles scatter mutually incoherently, intensities add. For illuminated particles,
The cross correlation vanishes by reflection symmetry. The degree of linear polarization is : random anisotropic scatterers produce partially depolarized light, even though every individual needle produces a single linear polarization.
Figure 2.
Independent isotropically oriented needles give the analyser ratio one to three, while coherent mean amplitudes have no x component
. Source. Details.
Continuous homogeneous dielectric. One must sum amplitudes coherently and include their spatial phases, rather than retain the preceding sum of intensities. Random orientations give , so the averaged response is isotropic and its induced polarization is along . Its coherent radiation toward , when allowed by the sample geometry, is therefore purely polarized, with . The independent-needle prefactor is no longer a quantitative dense-medium model; local-field and collective effects determine the effective permittivity.
There is also coherent continuum phase matching: a uniform bulk driven by a plane wave has no arbitrary outgoing plane-wave component. In an ideal infinite homogeneous medium there is no diffuse side scattering at all; propagation, refraction and reflection are determined by the medium and its interfaces. A finite sample's coherent side signal depends on its shape and boundaries. Thus the old diffuse ratio is replaced by a coherent isotropic-medium response, not by a universal nonzero side-scattering intensity.
Molecules in clean air are much smaller than visible wavelengths. Their induced electric dipole moments undergo Rayleigh scattering. Away from molecular resonances, the induced dipole amplitude is approximately , with only weak wavelength variation in the polarizability . The radiation-zone electric field is proportional to , so the scattered intensity satisfies
This is the Rayleigh scattering wavelength dependence. For equal incident intensities, light scatters about times as strongly as light. Thus light scattered out of the direct solar beam and into our line of sight is enriched in blue wavelengths.
The dipole mechanism also determines the Rayleigh scattering angular pattern. For unpolarized light, its angular factor is ; near scattering, the scattered light has strong linear polarization. The sky's perceived colour depends additionally on the solar spectrum, atmospheric transmission and the eye's spectral response, so the law alone does not predict a violet sky. Over a long near-horizontal path the direct beam loses more short-wavelength light, explaining red sunsets. Large aerosol particles instead produce Mie scattering, which is less strongly wavelength-selective and can make haze look whitish.
Put , where is the binding spring constant rather than an optical wavenumber. The steady complex amplitude solution of the forced Lorentz oscillator model is
The acceleration complex amplitude is . Hence
This is the scattering cross-section of a harmonically bound charge. The displacement and induced electric dipole moment are parallel to the incident linear polarization, and integrating their dipole radiation over all directions gives the total cross section above.
For , binding dominates and
recovering Rayleigh scattering with static electric polarizability . For , inertia dominates: , so its second time derivative becomes frequency independent and
This is a classical Thomson scattering regime only while . Exactly at , an undamped driven oscillator has no finite steady-state amplitude; the divergent formula must not be used. Finite damping and radiation reaction regulate the physical resonant response. Coherent contributions from many charges in an extended medium can also require a structure calculation rather than multiplying this single-charge result by a particle count.
In Rayleigh scattering, a small bound system develops an induced electric dipole moment . Well below its binding resonances, the electric polarizability is approximately its static value , independent of frequency. The mean radiated power and incident light intensity are
Dividing gives the total scattering cross section
This Rayleigh scattering wavelength dependence comes from the two time derivatives of an oscillating electric dipole moment, whose radiated field is proportional to . It requires the same scatterer with approximately constant low-frequency electric polarizability; resonances and strong optical dispersion invalidate the simple fourth-power law.
Now take , with the same angle to the incident direction. The perpendicular direction is and the in-plane transverse direction is . The roles of the two electric dipole populations interchange: electric dipoles radiate at full strength, whereas the observed wave amplitude from a electric dipole is reduced by .
Thus and , giving
This is also obtained by the symmetry , in part (a). For , both planes give the familiar Rayleigh scattering degree of polarization : the diffuse signal is unpolarized forward and fully linearly polarized at right angles. At a dark direction with zero denominator, no degree of polarization can be measured.
The given static induced surface charge and the dipole integral from part (b) give
Here is the dimensionless relative permittivity, as required by the stated expression; is the SI electric polarizability. Take and . In the Rayleigh scattering regime the field is approximately uniform across the sphere. Its induced dipole radiates mean power
The incident time-averaged light intensity in vacuum is . Dividing power by intensity gives the total cross section,
The units are ; the comparison is with the geometric projected area , not with a length:
The inequality holds for an ordinary nonabsorbing dielectric, with finite contrast factor and . It explains why a sufficiently small transparent sphere scatters weakly and why its scattered power has the characteristic dependence. Near a material resonance where is small, this simple inequality and static-polarizability approximation need reconsideration, including loss and radiative corrections. For the stated low-frequency ordinary dielectric the supplied polarizability is the intended approximation.
The far fields obey . The radial component of their Poynting vector is . Averaging and integrating gives
For a small isotropically polarizable particle, and incident intensity is . Thus the total Rayleigh scattering cross-section is
Each incident linear polarization gives the same integrated power; averaging unpolarized light does not add another factor of two. The scalar real polarizability approximation neglects absorption and molecular anisotropy.
For Rayleigh attenuation from refractivity in dilute air, polarization density is . Consequently and, with and ,
The ideal-gas number density at the supplied pressure and temperature is , yielding . For an explicit representative visible wavelength ,
Direct-beam intensity follows . Equivalently . Its strong wavelength dependence gives roughly – across – if refractivity is held fixed. This is a uniform-density dry-air scattering estimate, not absorption or transmission through a vertically stratified atmosphere; aerosols and molecular depolarization corrections are outside the given model.