Radial Structure of Debris Discs

• Radial structure of debris discs is defined by inner and outer edges, rings, gaps, and halos that serve as diagnostic markers of collisional cascades and dynamical perturbations.
• Observations using ALMA, Herschel, HST, and JWST combined with parametric and non-parametric models quantify features like edge sharpness and fractional width to infer planet formation history.
• Statistical analyses reveal that collisional evolution and planetary influences shape evolving disc architectures, with trends in belt radius, width, and halo brightness informing models of disc dynamics.

Debris discs are circumstellar structures composed of dust and planetesimals, analogous to the Kuiper Belt and Zodiacal light in the Solar System. The dust is constantly replenished by the collisional grinding of planetesimals, resulting in a complex radial structure defined by key features such as inner and outer edges, rings, gaps, and extended halos. The radial architecture of these discs encodes fundamental information about planet formation, collisional evolution, and dynamical sculpting by planets or stellar companions. With the advent of high-resolution imaging from ALMA, Herschel, HST, and JWST, as well as robust non-parametric and dynamical modeling frameworks, the detailed radial structure of debris discs is now accessible for statistical analysis across large samples.

1. Observational Diagnostics and Canonical Metrics

Radial structure in debris discs is characterized by several key parameters measurable via resolved imaging or reconstructed from spectral energy distributions (SEDs) (Hughes et al., 2018, Han et al., 12 Feb 2025, Han et al., 20 Jan 2026):

• Inner edge ($R_{\rm in}$): Sharp truncation, often reflecting dynamical clearing by planets or sublimation boundaries. Inner edges can be extremely steep—$\alpha_{\rm in} \gg 2$—in planet-truncated belts, or shallow (e.g. $\Sigma(r) \propto r^{2}$) when set by collisional evolution (Blanco et al., 2023, Han et al., 20 Jan 2026).
• Outer edge ($R_{\rm out}$): The maximum detected radius, often marking a transition to a low surface-brightness halo dominated by small grains.
• Characteristic radius ($R_{\rm belt}$): Typically identified with the peak in surface brightness, or the mean radius of a ring/belt (Hughes et al., 2018, Marshall et al., 2014).
• Fractional width ($\Delta r / r$): The full width at half maximum (FWHM) of the profile normalized by the centroid radius; can distinguish narrow rings ($\Delta r / r \lesssim 0.2$) from broad belts ($\Delta r / r \sim 0.5$ or higher) (Han et al., 20 Jan 2026, Han et al., 12 Feb 2025, Marshall et al., 2023).
• Surface density profile ($\Sigma(r)$): Commonly parameterized as a power-law or Gaussian, with broken/sloped regions modeling edges, gaps, or halos.

Resolved surface-brightness profiles in scattered light or thermal emission directly trace $\Sigma(r)$ convolved with the appropriate kernel (scattering phase function or Planck function), enabling power-law fits $I(r) \propto r^{-m}$, with $m$ varying systematically by regime (typically $m \simeq 3-4$ outside dense rings) (Hughes et al., 2018, Thebault et al., 2023).

2. Physical Origins of Radial Structure: Collisions, Radiation Pressure, and Dynamical Sculpting

Radial distributions in debris discs arise from a combination of collisional processing, radiation pressure, and gravitational perturbations by planets or binary companions.

Collisional Cascade: The steady-state collisional grinding yields a classical size distribution $n(s)\propto s^{-3.5}$ (Han et al., 20 Jan 2026, Blanco et al., 2023), causing spatial redistribution via radiation pressure and collisional destruction.

• Inner profiles: In pure collisional evolution, interior to a critical radius $r_c$, the dust surface density characteristically rises as $\Sigma_{\rm dust}(r)\propto r^2$ (Blanco et al., 2023). This slope directly arises from the balance between collisional timescales and disk age.
• Halos: Outside the main belt, radiation pressure launches small grains onto high-eccentricity orbits, producing a power-law halo with $\Sigma(r)\propto r^{-1.5}$ and corresponding scattered-light surface brightness $I(r)\propto r^{-3.5}$ in the idealized isotropic case (Thebault et al., 2023, Thebault et al., 2012).
• Flattening mechanisms: In dense discs, unbound grains (with $s<s_{\rm blow}$), size-dependent scattering phase functions, and finite vertical resolution can flatten the outer SB slope from $-3.5$ to as shallow as $-2.3$ (Thebault et al., 2023).

Planetary Sculpting: Planets carve gaps and sharp edges, with the width of the gap related to planet mass and semimajor axis through the chaotic zone scaling $\Delta r/a_p \sim 1.3\mu^{2/7}$, where $\mu=M_p/M_*$ (Thebault et al., 2012, Han et al., 20 Jan 2026).

• Gaps induced by planets can be sharp in pure N-body models but are systematically filled in by collisions and the replenishment of small grains; the gap contrast decreases significantly for high optical depth disks ($\tau_0\gtrsim2\times10^{-3}$) (Thebault et al., 2012).
• Multi-ring structures and broad depleted annuli may be produced by secular resonances driven by two or more planets, naturally generating double-ring morphologies and ring offsets (Yelverton et al., 2018).

Stellar Companions: In binaries, the presence of a secondary star sets a critical semimajor axis $a_{\rm crit}$, beyond which orbits are unstable (Thebault et al., 2010, Thebault, 2011). However, steady collisional production and radiation pressure still populate the dynamically unstable zone with high-$\beta$ grains, and the “forbidden” region outside $a_{\rm crit}$ hosts a faint, collisionally maintained halo.

3. Parametric and Non-parametric Radial Profile Models

The functional fitting of resolved profiles relies on both parametric and non-parametric approaches:

Parametric Profiles (Han et al., 20 Jan 2026, Marshall et al., 2014, Hengst et al., 2020):

• Single/double power laws: $\Sigma(r) \propto r^{-p}$, broken at characteristic radii with sharpness exponent $\gamma$
• Gaussian and double-Gaussian rings: Used to describe both narrow and broad belts as well as multi-ring systems
• Composite models: Gap-carving or halo components added on top of power-law or Gaussian models.

Non-parametric Methods (Han et al., 12 Feb 2025, Han et al., 20 Jan 2026):

• Algorithms such as rave and frank allow for direct deprojection and deconvolution from resolved images or visibilities, yielding $\Sigma(r)$ and vertical thickness $H(r)$ with uncertainties independent of parametric model assumptions.
• Trends recovered from large samples show increasing belt radii and fractional widths with system age, and a strong positive correlation of outer-edge width with age, supporting evolutionary broadening via dynamical stirring and planet-induced scattering (Han et al., 12 Feb 2025).

Empirical Examples

System	$r_0$ (au)	$\Delta r / r_0$	Edge Model	Reference
HD 107146	$80\pm2$	$0.25\pm0.03$	Double power law, gaps	(Han et al., 12 Feb 2025)
HD 16743	$158$	$0.50\pm0.05$	Gaussian belt, broad	(Marshall et al., 2023)
HR 4796A	$76.4$	$<0.05$	Gaussian core, $r^{-3.5}$	(1908.10378)
Vega	$85$	$0.18$	Gaussian, exponential halo	(Sibthorpe et al., 2010)

4. Halos, Gaps, and Edge Phenomenology

Halos: Universal halos arise from the injection of small, radiation-pressure-affected grains onto eccentric orbits outside the main belt (Thebault et al., 2023, Thebault et al., 2012, Thebault et al., 2010). The halo surface density is set by a combination of the collisional production rate and the sink timescale (dynamical ejection in binaries or unbound orbits for $s<s_{\rm blow}$). Typical slopes for the scattered-light surface brightness range from $-3.5$ (ideal isotropic) to $-3.0$ or shallower once the effect of unbound grains, non-isotropic SPFs, or finite resolution is included (Thebault et al., 2023). In thermal emission, halos dominate only at $\lambda\leq90\mu$m; at longer wavelengths, their contribution to the flux declines to a few percent.

Gaps and Ring Multiplicity: High-resolution surveys find that a substantial fraction of discs ($\sim$50% in the ARKS ALMA survey) exhibit more structure than a single ring, including wide gaps, substructure, and faint broad halos (Han et al., 20 Jan 2026). Gap widths and depths are closely linked to planetary mass and disc optical depth, but collisions act to diminish the contrast of dynamically induced gaps.

Edge Slopes: Steep inner edges ($\alpha_{\rm in} \gtrsim 15$) are markers of recent or ongoing planet-disk sculpting; in contrast, shallow edges ($\Sigma \propto r^2$ interior rise) are indicative of pure collisional relaxation, with the inner knee radius $r_c$ marking where the largest planetesimals attain collisional equilibrium (Blanco et al., 2023).

5. Eccentricity, Asymmetries, and Resonances

Eccentric Belts: When significant free or forced eccentricity is present, overlapping orbits in narrow rings create distinctive radial features:

• For constant eccentricity, narrow rings display two radial brightness maxima, at pericenter and apocenter; for broader belts or lower spatial resolution, these merge, resulting in apocenter or pericenter glow depending on geometric and observational parameters (Lovell et al., 2023).
• Profiles are sensitive to width-to-radius ($w=\Delta a/a_0$) and eccentricity ($e$); the critical values $w_{\rm crit,a}=2e/(1+e)$ and $w_{\rm crit,p}=2e/(1-e)$ set the regime where the twin peaks merge (Lovell et al., 2023).
• Azimuthal asymmetries (e.g., pericenter glow) can arise from both eccentricity and collisionally enhanced dust production near pericenter, as seen in HR 4796A (1908.10378), or from dynamically forced pericenter alignment by planets (Löhne et al., 2017).

Secular Resonances and Multi-planet Gaps: Gapped, double-ring, or spiral features can be imprinted by secular resonances of interior multi-planet systems. The characteristic location and width of the resonance-induced gap depends on planet masses, semi-major axes, and eccentricities; typically, one resonance is broad and rapidly depletes surface density, while the other remains narrow (Yelverton et al., 2018).

6. Statistical Trends, Evolution, and Theoretical Implications

Large-sample, non-parametric studies find that the characteristic radius of debris belts scales roughly as $r_0\propto t^{0.26}$ (with age $t$), consistent with "self-stirring" models where the collisional cascade propagates outward (Han et al., 12 Feb 2025). Both fractional width and outer-edge width increase with age, with the broadening likely tracing a combination of collisional evolution and planetary perturbation-induced scattering.

Sharp edges, deep gaps, and high ring multiplicity can often be attributed to planetary sculpting, with inferred planet masses for ring-truncating bodies typically $0.1$–$5\,M_{\rm J}$ for belts at tens of au (Han et al., 20 Jan 2026). In contrast, broad, shallow-edged, or evolving belts require a combination of collisional physics (governing the r$^{2}$ inner rise and the $r^{-1.5}$–$r^{-3.5}$ outer halo) and dynamical stirring.

Comparison to protoplanetary discs indicates that the width distribution of debris rings substantially overlaps with substructure in young, gas-rich systems, but also that a population of especially broad debris belts (with $\Delta r / r > 0.8$) must be broadened post-gas phase by planet migration, scattering, or dynamical heating (Han et al., 20 Jan 2026).

7. Modeling Methodologies and Limitations

A suite of complementary modeling approaches are employed:

• Analytical and semi-analytical models provide fast prediction and inversion of azimuthally averaged profiles in terms of underlying semi-major axis and eccentricity distributions, enabling direct comparison with imaging data (Rafikov, 2022).
• Radiative transfer codes (e.g., MCFOST, RADMC-3D, HYPERION, GRaTeR) self-consistently fit the multi-wavelength SED and resolved images, accounting for grain size, composition, and stochastic heating (Hengst et al., 2020, Marshall et al., 2014).
• Dynamical and collisional simulations detail the interplay between gravitational perturbations and collisional cascades, revealing the persistence (or smearing) of rings, gaps, and halos (Thebault et al., 2012, Löhne et al., 2017, Thebault, 2011).
• Non-parametric inversion (e.g., rave, frank) robustly recovers radial profiles without bias from assumed functional forms, with empirical uncertainties derived via Monte Carlo trials or Gaussian-process regularization (Han et al., 12 Feb 2025, Han et al., 20 Jan 2026).

Limiting factors include beam convolution/resolution, sensitivity to low surface-brightness features, and degeneracies in SED-only fitting. Separating belt and halo components, as well as accounting for unbound grain populations and inclination effects, is essential for accurate characterization (Thebault et al., 2023).

---

In summary, the radial structure of debris discs reflects a dynamic interplay among collisional processes, radiation pressure, and gravitational sculpting by both planets and binary companions. Canonical inner and outer edges, ring widths, and halo slopes can be traced to specific physical processes, while observed diversity—e.g., double rings, gaps, broad belts, asymmetric features—hints at the rich dynamical histories and evolutionary pathways of planetary systems (Hughes et al., 2018, Han et al., 12 Feb 2025, Han et al., 20 Jan 2026, Blanco et al., 2023, Lovell et al., 2023, Thebault et al., 2023).