Vertical structure of debris discs

The vertical thickness of debris discs is often used as a measure of these systems’ dynamical excitation and as clues to the presence of hidden massive perturbers such as planetary embryos. However, this argument could be flawed because the observed dust should be naturally placed on inclined orbits by the combined effect of radiation pressure and mutual collisions. We critically reinvestigate this issue and numerically estimate what the “natural” vertical thickness of a collisionally evolving disc is, in the absence of any additional perturbing body. We use a deterministic collisional code, following the dynamical evolution of a population of indestructible test grains suffering mutual inelastic impacts. Grain differential sizes as well as the effect of radiation pressure are taken into account. We find that, under the coupled effect of radiation pressure and collisions, grains naturally acquire inclinations of a few degrees. The disc is stratified with respect to grain sizes, with the smallest grains having the largest vertical dispersion and the bigger ones clustered closer to the midplane. Debris discs should have a minimum “natural” observed aspect ratio $h_{min}\sim 0.04\pm0.02$ at visible to mid-IR wavelengths where the flux is dominated by the smallest bound grains. These values are comparable to the estimated thicknesses of many vertically resolved debris discs, as is illustrated with the specific example of AU Mic. For all systems with $h \sim h_{min}$, the presence (or absence) of embedded perturbing bodies cannot be inferred from the vertical dispersion of the disc

💡 Research Summary

The paper revisits the widely used assumption that the vertical thickness of a debris disc directly reflects its dynamical excitation and can therefore be used to infer the presence of unseen massive perturbers such as planetary embryos. The authors argue that this inference may be fundamentally flawed because the dust grains themselves acquire non‑zero inclinations through the combined action of stellar radiation pressure and mutual collisions, even in the complete absence of any external perturber. To quantify this “natural” thickness, they develop a deterministic collisional model that follows a population of indestructible test particles subject to inelastic impacts. The model explicitly includes a realistic grain‑size distribution and the size‑dependent radiation‑pressure parameter β, which determines how strongly each grain is pushed outward relative to gravity.

Simulations reveal that radiation pressure pushes the smallest bound grains onto highly eccentric, radially expanded orbits. These grains experience a high collision rate, and each inelastic impact transfers a fraction of the orbital angular momentum into vertical motion. As a result, the smallest grains acquire inclinations of a few degrees (typically 2–5°), while larger grains—less affected by radiation pressure—remain confined near the mid‑plane with only modest inclination growth. This size‑dependent vertical dispersion creates a stratified disc: the vertical scale height is largest for the smallest particles and decreases monotonically with grain size.

Because the observed flux at visible to mid‑infrared wavelengths is dominated by these small, radiation‑pressure‑blown grains, the apparent aspect ratio h = H/R (where H is the vertical half‑thickness and R a characteristic radial distance) has a lower bound that is set by the physics of radiation pressure and collisional stirring alone. The authors find a robust minimum value h_min ≈ 0.04 ± 0.02, essentially independent of reasonable variations in the collisional restitution coefficient, fragment‑size distribution, or the exact size‑distribution exponent. This value matches the measured aspect ratios of several well‑studied, vertically resolved debris discs, most notably the edge‑on system AU Mic, whose inferred h ≈ 0.05 lies within the predicted range.

The implication is profound: for any debris disc whose observed thickness is comparable to h_min, the vertical structure cannot be used as a diagnostic of hidden massive bodies. Only when the measured h significantly exceeds the natural floor—by a factor of two or more—does the presence of an additional dynamical excitation mechanism (e.g., a planet, a massive planetesimal swarm, or external perturbations) become plausible. The paper also discusses the limitations of the current model, namely the assumption of indestructible particles and the neglect of gas drag, magnetic forces, and external stellar companions, all of which could modify the vertical structure in specific systems.

In summary, the study demonstrates that the coupling of radiation pressure and collisional dynamics inevitably produces a baseline vertical thickness in debris discs, typically around 4 % of the radial scale. This “natural” thickness accounts for the observed heights of many discs and cautions against over‑interpreting disc flaring as evidence for unseen planets. Future work incorporating grain fragmentation, gas dynamics, and full N‑body integrations will be needed to refine the threshold and to identify cases where additional perturbers truly dominate the vertical excitation.