Planetary Infrared Excess (PIE) Overview

• Planetary Infrared Excess (PIE) is defined as the surplus IR flux beyond stellar atmospheric predictions, indicating additional circumstellar sources.
• Detection relies on precise flux ratios and SED fitting, using high-fidelity photometry and robust calibration to identify debris disks and planetary companions.
• PIE analysis drives exoplanet discovery and atmospheric characterization, providing insights into system evolution and the prevalence of circumstellar material.

Planetary Infrared Excess (PIE) refers to the observed surplus of infrared (IR) flux in a spatially unresolved system above the level predicted for the primary star’s photosphere alone. Quantitatively, it is the residual thermal emission—typically at mid- or far-IR wavelengths—not accounted for by stellar atmospheric models, and is attributed to circumstellar material (e.g., debris disks, planetary companions) or to planetary thermal emission in the star+planet system. PIE provides a direct probe of extrasolar planetary systems, including debris from planetesimal collisions, young giant planets, white dwarf survivors, and, with advanced retrieval methods, the atmospheres of both transiting and non-transiting exoplanets.

1. Formal Definition and Physical Basis

PIE is defined as the excess photometric or spectroscopic signal, at specified IR wavelengths, above the stellar model atmosphere. This is evaluated by comparing either observed IR flux densities, flux ratios (e.g., $f_{22}/f_{12}$ for WISE bands), or the full spectral energy distribution (SED), to synthetic photospheric predictions that match the star's $T_{\mathrm{eff}}$, $\log g$, and metallicity:

$$\mathrm{PIE}(\lambda) = F_{\rm obs}(\lambda) - F_*(\lambda),$$

where $F_{\rm obs}$ is the measured flux and $F_*$ is the model photosphere. In the context of combined light from a star and a non-transiting planet, the total flux is:

$$F_{\rm tot}(\lambda) = F_*(\lambda) + F_p(\lambda),$$

with $F_p$ representing the planet's thermal emission. The PIE signal is thus $F_p(\lambda)$, or equivalently, the contrast ratio $F_p(\lambda) / F_*(\lambda)$ (2207.13727, Stevenson, 2020).

In multi-object systems or for population studies, PIE can also refer to the ensemble rate of detected IR excesses above a defined threshold, serving as a statistical measure of warm/cold debris or planetary emission occurrence (Maldonado et al., 2017, Wang et al., 13 Apr 2026).

2. Detection Methodologies and Statistical Significance

Detection of PIE relies on high-fidelity photometry and robust calibration against stellar models:

• Flux Ratio Approach: For surveys such as WISE, the observed flux ratio at two IR bands (e.g., 22μm/12μm) is compared to synthetic ratios from convolved ATLAS9 or PHOENIX spectra. The excess parameter,

$T_{\mathrm{eff}}$0

is evaluated, where $T_{\mathrm{eff}}$1 and $T_{\mathrm{eff}}$2 are observed and synthetic ratios, and $T_{\mathrm{eff}}$3 includes photometric, calibration, and model uncertainties. Significant excess is flagged when $T_{\mathrm{eff}}$4 (3σ) or higher (Maldonado et al., 2017).

• SED Fitting: For white dwarfs, fluxes from optical to mid-IR bands are simultaneously fit with white dwarf atmosphere models (e.g., Koester models), optionally adding a blackbody disk, brown-dwarf atmosphere, or additional component for excess. A source exhibits PIE if one or more IR points exceed model predictions by >3σ (Wang et al., 13 Apr 2026, Xu et al., 2020).
• Empirical Stellar Subtraction: For exoplanet systems observed in combined light (with JWST or similar), the stellar SED can be empirically anchored using in-eclipse or out-of-transit data, with the PIE defined as the normalized planet–star flux ratio:

$T_{\mathrm{eff}}$5

(Lustig-Yaeger et al., 7 Oct 2025).

3. Origins: Physical Scenarios for PIE Signals

The astrophysical origins of PIE encompass several physical processes:

• Debris Disks: Thermal emission from collisionally produced dust, often tracing planetesimal belts in planetary systems. Debris disks are diagnosed by their spectral slope and fractional luminosity $T_{\mathrm{eff}}$6, with dust temperatures from $T_{\mathrm{eff}}$7 K (cold, $T_{\mathrm{eff}}$8100 μm) to $T_{\mathrm{eff}}$91000–1800 K (warm; near- and mid-IR). Classical examples include both single and binary star debris environments (Maldonado et al., 2017, Matranga et al., 2010).
• Planetary Companions: For white dwarfs, brown dwarfs, or giant planet remnants at projected separations $\log g$0 few AU, the companion’s blackbody or atmospheric spectrum dominates beyond $\log g$1 μm, giving rise to a pronounced PIE. Unresolved mid-IR excess in WD systems is robustly associated with cold (200–300 K) giant planets at 0.1–10 AU (Limbach et al., 2024, Poulsen et al., 26 Jan 2026).
• Multiple Planets: In tightly packed multi-planet systems (e.g., TRAPPIST-1), the aggregate PIE signal is modeled as a sum of planetary emission components, each parametrized by radius, temperature/albedo, and orbital separation, with spectral separation or retrieval required to disentangle them (Mayorga et al., 2023).
• Phase-Dependent Thermal Emission: For transiting and non-transiting exoplanets, time-resolved PIE or phase curves can reveal the distribution of thermal emission (day/night contrast, atmospheric circulation regimes) (Johnson et al., 6 Feb 2026, Lustig-Yaeger et al., 7 Oct 2025).

4. Instrumentation, Calibration, and Retrieval Strategies

Detection and interpretation of PIE require stringent calibration and sophisticated modeling:

• Instrumental Calibration: The AllWISE photometry requires band-dependent linear corrections to match synthetic spectra, with offsets up to 10–40% in some bands, especially at low fluxes. Systematic errors in calibration dominate the error budget at high SNR (Maldonado et al., 2017).
• Spectral Coverage: Broad simultaneous wavelength coverage (1–18 μm or wider) is critical. JWST NIRISS+NIRSpec+MIRI combinations or future missions (e.g., MIRECLE) are designed to maximize PIE sensitivity for both warm and temperate planetary regimes (2207.13727, Hammond et al., 1 Apr 2025).
• Bayesian Retrievals: Joint fitting frameworks, often using MCMC or nested sampling, recover planetary and stellar parameters along with atmospheric compositions (CO$\log g$2, H$\log g$3O, O$\log g$4, CH$\log g$5, etc.). For non-transiting planets, radius–temperature degeneracy can be broken if both the Wien and Rayleigh–Jeans portions of the planetary SED are observed (Lustig-Yaeger et al., 2021, Stevenson, 2020).

| Instrument/Technique | Wavelength Range | Key Usage | |--------------------------|------------------|-----------------------------------| | WISE/AllWISE | 3.4–22 μm | Survey for debris/planet PIE | | JWST MIRI Imager/Spectra | 5–28 μm | Spectroscopy of PIE at high SNR | | ExoCAM+PSG/MIRECLE | 1–18 μm | GCM-driven PIE spectra/surveys |

• Empirical and Data-Driven Approaches: Variable PIE (VPIE) employs empirical stellar templates derived from the data, reconstructing the stellar SED across epochs and isolating the planet’s phase-dependent signal from the residuals (Johnson et al., 6 Feb 2026).

5. Quantitative Results and Prevalence

PIE occurrence rates and quantitative results vary by astrophysical context and wavelength:

• Main Sequence Stars: The detection rate of warm PIE ($\log g$622 μm) is extremely low: $\log g$7 for both planet-hosting and non-hosting solar-type stars; no statistically significant difference exists between hosts and non-hosts. Cold debris disks ($\log g$870–100 μm) are an order of magnitude more common ($\log g$910–20\%) (Maldonado et al., 2017).
• White Dwarfs: IR-excess occurrence rates among bright WDs are 6.6–8.4%, encompassing both dust disks and cool companions. Newly identified systems extend dust disk detection to cooling ages $$\mathrm{PIE}(\lambda) = F_{\rm obs}(\lambda) - F_*(\lambda),$$05 Gyr (Wang et al., 13 Apr 2026, Xu et al., 2020).
• Planets around White Dwarfs: PIE detections of unresolved, cold giant planets at $$\mathrm{PIE}(\lambda) = F_{\rm obs}(\lambda) - F_*(\lambda),$$1–$$\mathrm{PIE}(\lambda) = F_{\rm obs}(\lambda) - F_*(\lambda),$$2 AU around WDs (T$$\mathrm{PIE}(\lambda) = F_{\rm obs}(\lambda) - F_*(\lambda),$$3250–260 K, $$\mathrm{PIE}(\lambda) = F_{\rm obs}(\lambda) - F_*(\lambda),$$43–7 M$$\mathrm{PIE}(\lambda) = F_{\rm obs}(\lambda) - F_*(\lambda),$$5) are now possible with JWST/MIRI, ruling out simple dust disk origins via spectral fits and systematics arguments (Limbach et al., 2024, Poulsen et al., 26 Jan 2026).
• Phase-Resolved and Atmospheric Characterization: For nearby transiting and non-transiting exoplanets, PIE can deliver SNR$$\mathrm{PIE}(\lambda) = F_{\rm obs}(\lambda) - F_*(\lambda),$$6 for temperate rocky planets within $$\mathrm{PIE}(\lambda) = F_{\rm obs}(\lambda) - F_*(\lambda),$$75 pc in 1.5–2 m-class telescopes (Hammond et al., 1 Apr 2025, 2207.13727). Retrievals yield atmospheric molecular abundances to 0.5–1 dex, dayside/nightside temperatures to 10–20 K, and Bond albedos to $$\mathrm{PIE}(\lambda) = F_{\rm obs}(\lambda) - F_*(\lambda),$$80.1 (Lustig-Yaeger et al., 7 Oct 2025).

6. Astrophysical and Methodological Implications

• Planet Formation and System Evolution: PIE from debris disks traces ongoing planetesimal collisions and the late dynamical evolution of planetary systems, including post-main-sequence survival and the outcomes of binary interactions (Maldonado et al., 2017, Matranga et al., 2010).
• Atmospheric Science of Non-Transiting Planets: PIE enables retrieval of atmospheric properties (composition, temperature profile, cloud structure, circulation regimes) well beyond the transiting exoplanet sample, leveraging broad IR spectral coverage (Hammond et al., 1 Apr 2025, Lustig-Yaeger et al., 2021).
• Planet Discovery and Multiplicity: The PIE technique, when applied in multi-planet systems and with next-generation sensitivity, can identify previously unknown planets via their stacked IR excess signatures, though semi-major axis constraints can be degenerate without external priors (Mayorga et al., 2023).
• Strong Calibration/Model Dependency: The reliability of PIE results rests on absolute flux calibration, accurate stellar models, robust treatment of blended or background sources, and systematic error mitigation, especially for small planet/star contrast regimes (10⁻⁴–10⁻⁶) (Maldonado et al., 2017, Wang et al., 13 Apr 2026).

7. Open Questions and Future Prospects

• Pushing Sensitivity Boundaries: Future missions with $$\mathrm{PIE}(\lambda) = F_{\rm obs}(\lambda) - F_*(\lambda),$$91–5 ppm stability, broad 1–20 μm coverage, and cryogenic telescopes (e.g., MIRECLE) are required for PIE detection of true Earth analogs around M dwarfs and for robust non-transiting atmospheric spectroscopy (2207.13727, Hammond et al., 1 Apr 2025).
• Spectroscopic Discrimination in White Dwarfs: Confirming planetary origins of PIE in post-main-sequence systems demands mid-IR spectroscopy to distinguish molecular signatures of giant planets (H$F_{\rm obs}$0O, CH$F_{\rm obs}$1, CO$F_{\rm obs}$2) from silicate-rich debris spectra (Limbach et al., 2024, Poulsen et al., 26 Jan 2026).
• Temporal and Multi-Epoch Variability: PIE can map thermal phase curves and extract information on planet circulation and atmospheric dynamics via high-cadence observations, especially with empirical (VPIE) or phase-resolved techniques (Johnson et al., 6 Feb 2026, Lustig-Yaeger et al., 7 Oct 2025).
• Statistical and Demographic Constraints: As the catalog of PIE systems grows (e.g., from DESI, MEOW, MEAD), long-term monitoring and multi-wavelength follow-up will refine occurrence rates, correlations with planetary architectures, and the fate of planetary systems through stellar evolution (Wang et al., 13 Apr 2026, Limbach et al., 2024).

PIE thus remains a central diagnostic in both the census of circumstellar debris and emergent exoplanet atmospheric science, with its full potential contingent on advances in mid-IR instrumentation, modeling sophistication, and statistical characterization of both stellar and planetary populations.