The History of the Solar Systems Debris Disc: Observable Properties of the Kuiper Belt

The Nice model of Gomes et al. (2005) suggests that the migration of the giant planets caused a planetesimal clearing event which led to the Late Heavy Bombardment (LHB) at 880 Myr. Here we investigate the IR emission from the Kuiper belt during the history of the Solar System as described by the Nice model. We describe a method for easily converting the results of n-body planetesimal simulations into observational properties (assuming black-body grains and a single size distribution) and further modify this method to improve its realism (using realistic grain properties and a three-phase size distribution). We compare our results with observed debris discs and evaluate the plausibility of detecting an LHB-like process in extrasolar systems. Recent surveys have shown that 4% of stars exhibit 24 um excess and 16% exhibit 70 um excess. We show that the Solar System would have been amongst the brightest of these systems before the LHB at both 24 and 70 um. We find a significant increase in 24 um emission during the LHB, which rapidly drops off and becomes undetectable within 30 Myr, whereas the 70 um emission remains detectable until 360 Myr after the LHB. Comparison with the statistics of debris disc evolution shows that such depletion events must be rare occurring around less than 12% of Sun-like stars and with this level of incidence we would expect approximately 1 of the 413 Sun-like, field stars so far detected to have a 24 um excess to be currently going through an LHB. We also find that collisional processes are important in the Solar System before the LHB and that parameters for weak Kuiper belt objects are inconsistent with the Nice model interpretation of the LHB.

💡 Research Summary

The paper investigates how the Kuiper Belt (KB) would have appeared as an infrared (IR) debris disc throughout the history of the Solar System, using the Nice model (Gomes et al. 2005) as a framework for the Late Heavy Bombardment (LHB) at ~880 Myr. The authors develop a two‑stage methodology for converting the output of N‑body planetesimal simulations into observable IR fluxes. In the first, simplified stage they assume black‑body grains and a single characteristic size for all bodies; orbital distances are used to compute equilibrium temperatures (T ∝ r⁻¹ᐟ²) and the total flux at any wavelength is obtained by summing the Planck emission of each particle weighted by its geometric cross‑section. This approach provides a rapid, order‑of‑magnitude estimate of the disc’s excess emission at 24 µm and 70 µm.

In the second, more realistic stage, the authors replace the black‑body assumption with wavelength‑dependent absorption and emission efficiencies (Q_abs, Q_emit) derived from Mie theory for astronomical silicates and water ice. They also introduce a three‑phase size distribution: (1) large planetesimals (km‑scale), (2) collisional fragments (10 µm–km), and (3) fine dust (sub‑10 µm). The evolution of this distribution is governed by a collisional cascade model that incorporates a fragmentation threshold Q_D* and a size‑dependent collision rate. Small grains are further subject to radiation pressure and Poynting‑Robertson drag, which set their removal timescales. By integrating these processes over the Nice‑model dynamical history, the authors generate time‑dependent synthetic spectral energy distributions (SEDs).

The results show that before the LHB the Solar System would have been among the brightest known debris discs. At 24 µm the excess would have been roughly 5–10 times the median of observed Sun‑like stars, while at 70 µm the excess would have been 3–5 times higher. The onset of the LHB, driven by the rapid migration of Jupiter and Saturn through mutual resonances, injects a large fraction of KB material into the inner Solar System. This raises the dust temperature, producing a sharp spike in 24 µm emission that reaches up to ~15 times the pre‑LHB level. However, the spike is short‑lived: collisional grinding and radiation‑pressure blow‑out remove the small grains within ~30 Myr, after which the 24 µm excess drops below current detection thresholds. In contrast, the 70 µm excess, dominated by larger, colder grains that remain at greater heliocentric distances, decays more slowly and stays above typical detection limits for roughly 360 Myr after the LHB.

To place these findings in an observational context, the authors compare the synthetic excess lifetimes with Spitzer and Herschel survey statistics: ~4 % of Sun‑like stars show a 24 µm excess, while ~16 % show a 70 µm excess. By matching the fraction of time a Solar‑type system would be observable in each band to the observed fractions, they infer that LHB‑like clearing events must be relatively rare, occurring in less than about 12 % of Sun‑like stars. Given the current sample of 413 Sun‑like field stars with measured 24 µm excesses, the model predicts that roughly one star should be caught in the brief 24 µm bright phase associated with an LHB.

A sensitivity analysis on the collisional parameters reveals that assuming weak Kuiper Belt objects (low Q_D*) leads to an overly rapid depletion of dust, producing IR excesses far below the observed levels both before and after the LHB. Conversely, adopting stronger bodies (higher Q_D*) yields excesses that match the data, implying that the real KB population must be relatively robust against fragmentation. This conclusion challenges versions of the Nice model that rely on a very fragile planetesimal belt to trigger the LHB.

Overall, the paper provides a robust framework for translating dynamical planetesimal simulations into observable debris‑disc signatures, demonstrates that the Solar System would have been a conspicuous IR source before and briefly during the LHB, and uses current debris‑disc statistics to argue that such cataclysmic clearing events are uncommon among Sun‑like stars. The work also highlights the importance of realistic grain physics and collisional strength in interpreting debris‑disc observations and in testing dynamical models of planetary system evolution.