Zodiacal Emission Model (ZEM) Overview

• The Zodiacal Emission Model (ZEM) is a quantitative framework that models interplanetary dust emission in both Solar System and exoplanetary systems.
• It employs axisymmetric density profiles, thermal emission, and scattering physics to convert brightness measurements into meaningful dust level estimates.
• ZEM facilitates rigorous foreground subtraction and instrument calibration, crucial for detecting faint cosmic backgrounds and Earth-analogue exoplanets.

The Zodiacal Emission Model (ZEM) provides a quantitative framework for interpreting the infrared and optical emission associated with interplanetary dust (IPD) clouds, both in the Solar System (“zodi”) and in extrasolar planetary systems (“exo-zodi”). ZEMs enable rigorous prediction and subtraction of zodiacal foregrounds in both broadband photometry and high-resolution imaging, crucial for accurate measurements of faint astrophysical signals such as the cosmic infrared background and for the detection of Earth-analogue exoplanets. Current ZEMs encode the spatial, thermal, and scattering properties of the IPD population using physically-motivated density models, grain absorption/scattering physics, and survey-specific response and calibration procedures. This article synthesizes the governing concepts, core mathematical formalism, astrophysical context, calibration strategies, and limitations of state-of-the-art ZEMs, with specific attention to the axisymmetric parametric approach as exemplified in exozodi modeling for nulling interferometers (Kennedy et al., 2014), as well as leading Solar System implementations.

1. Physical and Mathematical Structure of the ZEM

The essential structure of the ZEM is an optically thin, axisymmetric belt of warm dust parameterized by its midplane surface density $\Sigma_m(r)$, representing the cross-sectional area per unit area as a function of stellocentric distance $r$: $$\Sigma_{m}(r) = z\,\Sigma_{m,0}\,\biggl(\frac{r}{r_{0}}\biggr)^{-\alpha}$$ where $z$ is the “zodi” scaling (unity for Solar System levels), $\Sigma_{m,0}$ is the normalization at reference radius $r_0$, and $\alpha$ is the radial power-law index. The model is truncated at inner and outer radii $r_\text{in}$ and $r_\text{out}$, typically chosen for Solar-analogue models as $r_\text{in} = 0.034~\text{AU}$ and $r_\text{out} = 10~\text{AU}$ with $r_0 = \sqrt{L_\star/L_\odot}~\text{AU}$ to ensure the temperature at $r_0$ matches the equilibrium blackbody temperature near $278$ K (Kennedy et al., 2014).

The ZEM computes the infrared surface brightness (specific intensity) at an observing wavelength $\lambda$ via: $$S_\text{disk}(\lambda, r) = 2.35\times10^{-11}~\Sigma_{m}(r)~B_\nu\bigl[\lambda,T_\text{BB}(r)\bigr]$$ where $B_\nu$ is the Planck function and the blackbody temperature profile for grains is

$$T_\text{BB}(r) = 278.3~L_\star^{1/4}~r^{-1/2}~\text{K}$$

This assumes $Q_\text{abs}\approx 1$ for sufficiently large grains in the mid-IR, consistent with optically efficient, astronomical silicate/carbonaceous populations. The model is inherently axisymmetric and neglects small-scale azimuthal asymmetries except as parameterized in specific survey configurations.

2. Nulling Interferometry and Model–Observation Mapping

For applications in exozodiacal dust detection with interferometric instruments such as the Large Binocular Telescope Interferometer (LBTI), the ZEM is coupled to the spatial transmission profile of the instrument. The LBTI implements a one-dimensional fringe pattern with sky-plane transmission: $$T_\text{null}(x) = \sin^2\!\left(\frac{\pi x}{2\phi_\text{null}}\right), \quad \phi_\text{null} = \frac{\lambda}{2B}$$ with $B$ the baseline and $x$ the sky-coordinate perpendicular to the fringes. The observable quantity is the “null depth,” representing the fraction of total disk flux transmitted relative to the stellar point source, computed as: $$\text{null}(z) = \frac{1}{F_\star} \int_{r_\text{in}}^{r_\text{out}} 2\pi r\,S_\text{disk}(r)\langle T_\text{null}\rangle_\theta \, dr$$ where $F_\star$ is the stellar flux and $\langle T_\text{null}\rangle_\theta$ is averaged over disk orientations. The dust level $z_\text{obs}$ inferred from an observation is

$$z_\text{obs} = \frac{\text{null}_\text{obs}}{\text{null}_\text{model}(z=1)}$$

allowing direct linkage between the observed interferometric signal and the underlying IPD column (Kennedy et al., 2014).

3. Scattered-Light Predictions and Empirical Albedo Calibration

The ZEM further enables prediction of dust-scattering brightness at optical/near-IR wavelengths by invoking an empirically motivated single-scattering albedo $\omega$ (typically $\sim0.1$ for Solar System dust analogues) and approximating isotropic phase function behavior. The modeled scattered-light surface brightness at radius $r$ is given by: $$S_\text{sca}(r) = F_{\nu,\star}~4\pi\Bigl(\frac{d}{r}\Bigr)^2 \frac{\omega}{1-\omega}~\Sigma_m(r)$$ with $d$ the system distance and $F_{\nu,\star}$ the stellar flux density at that wavelength, assuming $Q_\text{sca}/(Q_\text{abs}+Q_\text{sca}) = \omega$ and $\Sigma_m\approx \Sigma_\text{true}(1-\omega)$. This formulation connects thermal-emission in the mid-IR to optical scattered-light limits, essential for direct-imaging exoplanet mission yield estimates and for interpreting upper limits or detections in multi-band surveys (Kennedy et al., 2014).

4. Parameterization and Sensitivity Regimes

The reference Solar System ZEM employs five principal parameters: $$(r_\text{in},~r_\text{out},~r_0,~\alpha,~\Sigma_{m,0})$$. For $L_\star=L_\odot$, canonical values are $r_\text{in}=0.034~\text{AU}$, $r_\text{out}=10~\text{AU}$, $r_0=1~\text{AU}$, $\alpha=0.34$, $\Sigma_{m,0}=7.12\times10^{-8}$. These selections reproduce the COBE/DIRBE-determined optical depth at $1~\text{AU}$. Survey sensitivity is governed by the single-measurement calibrated null uncertainty $\sigma_\text{null}$ (e.g., $10^{-4}$ for LBTI), translating into dust level limits of $z\sim3$–$10$ zodi for Solar analogues and $z\lesssim1$–$3$ zodi for more luminous stars (Kennedy et al., 2014).

Tabulated sensitivities:

Host Type	Null Depth ($\sigma_{\rm null}$)	Zodi Level Sensitivity	Scattered-Light Limit (mag arcsec⁻²)
Sun-like ($T_\star\sim5800$ K)	$10^{-4}$	$z\sim4$	$20$–$21$
Early-type (A)	$10^{-4}$	$z\sim1$	$>$21

For specific targets, e.g., $\eta$~Crv measured at null $=4.4\%\pm0.35\%$, ZEM yields $z\simeq1376^{+497}_{-140}$ and $S_\text{sca}\sim4.0^{+1.4}_{-0.4}$~mJy arcsec⁻² ($\sim15$~mag arcsec⁻²) at $r_0=2.3$ AU (Kennedy et al., 2014).

5. Physical Motivation and Survey Interpretation

The ZEM formalism is anchored in several physical assumptions:

• The dust belt is optically thin and axisymmetric;
• Cross-sectional area is a power-law in $r$;
• Large astrosilicate/carbonaceous grains dominate the mid-IR thermally active population, yielding $Q_\text{abs}\simeq1$ and $T\propto r^{-1/2}$;
• Scattered-light predictions rely on a measured Solar System albedo and, for first-order calculations, an isotropic (i.e., phase-independent) scattering phase function.

This simplicity allows both robust forward modeling for instrument design and rapid inversion of observed null depths into astrophysically meaningful dust levels, which is crucial for placing upper bounds on habitable-zone dust around nearby stars and thus for assessing the feasibility of future direct imaging surveys (Kennedy et al., 2014).

6. Connections to Solar System ZEMs and Broader Context

The exozodi ZEM formalism mirrors, in dimensionless terms, the Solar System-focused ZEMs employed for COBE/DIRBE, AKARI, and Planck analyses. All such models share a structure comprising dominant smooth (fan) clouds, discrete dust bands (from asteroidal families), circumsolar ring and trailing features, and typically posit axisymmetry but allow for mean offsets and vertical tilts. In Solar System models, explicit fits for grain size, composition, emissivity modifications, and small-scale components are implemented to match all-sky datasets and to calibrate the dominant sources of error. ZEM parameterizations also facilitate extension to exozodiacal environments by appropriate normalization of $z$ and rescaling of disk radii with stellar luminosity (Kennedy et al., 2014).

7. Implications for Exoplanet Observations and Future Work

The ZEM enables translation of direct instrument measurements—whether null depths or scattered-light upper limits—into physically interpretable dust masses, optical depths, and observability constraints for terrestrial planets. Accurate ZEMs are essential both for robust exoplanet yield forecasts and for minimizing spurious signals in transit photometry or direct imaging regimes. Ongoing model refinements include improved treatment of dust grain composition and size distributions, validation against scattered-light imaging, and detailed morphology of resonant structures. Increasingly, joint Bayesian frameworks and multi-instrument datasets (e.g., coupling LBTI with AKARI or JWST sky brightness monitoring) are required for optimal parameter inference and for resolving degeneracies in ZEM component fits.

In summary, the ZEM, as described in the context of exozodiacal detection and Solar System IR backgrounds, provides a tractable, physically-motivated template that enables accurate, survey-wide statements about habitable-zone dust, the dominant noise source for future Earth analogue imaging (Kennedy et al., 2014).