The Role of Self-Gravity in Debris Disk Warp Formation: The Case of HD 110058

TL;DR

The paper investigates how self-gravity in a massive debris disk influences warp formation under perturbation from an inclined inner planet, using GPU-accelerated $N$-body simulations for the HD 110058 system. It demonstrates that disk self-gravity enforces a semi-rigid, coherent precession that can rapidly generate a global warp and drive the system toward a quasi-equilibrium warp state, consistent with observations after projecting into scattered-light images. An analytic Laplace-Lagrange framework reproduces the main warp modes and yields an empirical relation linking the equilibrium warp angle to planetary and disk parameters, allowing a dynamical constraint of $M_{disk} \lesssim 1000\,M_\oplus$ for HD 110058. The work provides a translatable method for constraining unseen planets and disk masses from warp morphologies and highlights the importance of self-gravity in interpreting debris-disk structures and their masses.

Abstract

We investigate the crucial role of self-gravity in the formation of warps in debris disks, focusing on the HD 110058 system as an example. Using advanced, GPU-accelerated $N$-body simulations, we model the gravitational dynamics of a massive planetesimal disk perturbed by an inclined, inner planet. Our simulations reveal that self-gravity fundamentally alters the disk's evolution compared to massless models. It enforces a coherent, semi-rigid precession of the disk and enables the rapid formation of a global warp structure within 0.5 Myr. The warp angle undergoes a damped oscillation, eventually settling into a quasi-equilibrium state. By generating synthetic scattered-light images, we demonstrate that our model successfully reproduces the observed S-shaped warp morphology of the debris disk in HD 110058, supporting the existence of an unseen planet. Furthermore, we derive an empirical relationship that connects the equilibrium warp angle to the physical parameters of the disk and the planet. Applying this relation to HD 110058, we constrain its disk mass to be likely less than 1,000 $M_\oplus$, offering a new dynamical perspective on the debris disk mass problem.

Figures (7)

• Figure 1: Example snapshots of the debris disk for the standard case ($M_\mathrm{disk}=100M_\oplus$)at different epochs, including 0 Myr (initial setup), 5 Myr, 10 Myr, 15 Myr, and 17 Myr (the estimated age of the HD 110058). The black dots indicate the positions of the particles. The red curve and dot represent the instantaneous orbit and positions of the planet, and the black plus symbols indicate the central star. The upper panels (a1)-(e1) display the face-on perspectives of the disk, and the lower panels (a2)-(e2) display a nearly edge-on perspective.
• Figure 2: The evolution of the unit angular momentum vectors of the concentric rings. Red and blue solid lines represent the $x$- and $y$-components of the ring's angular momentum. Red and blue dashed lines show those of the planet. The left five panels (a1-e1) display the case with disk mass of $100~M_\oplus$. The right five panels (a2-e2) show the case with a massless disk formed by test particles.
• Figure 3: Evolution of relative inclinations of the rings. Panel (a): Solid lines show the evolution of relative inclination of the ring with respect to the planetary orbital plane. Dashed lines indicate the linear fitting results. Panel (b): Evolution of the relative longitude of the ascending node, $\Delta \Omega_j$, of the rings with respect to the planetary orbital plane.. Panel (c): Green dot-dashed line with points indicates the inclination profile of the disk at 17 Myr, Blue line with points shows the initial equilibrium inclination profile from linear fitting of $N$-body simulation results. Panel (d): Solid line displays the evolution of the warp angle (relative angle between the innermost and outermost rings. Dashed line shows linear fitting result).
• Figure 4: Comparison between N-body simulation and theoretical analysis. Panels (a), (b), and (c): Evolution of the warp angle for disk masses of $50~M_\oplus$, $100~M_\oplus$, and $200 ~M_\oplus$, respectively. Solid lines represent the simulation results, while dashed lines show those from Laplace-Lagrange secular theory. Panels (d), (e), and (f): Initial equilibrium inclination profile for the same three disk masses. Solid lines are from simulations, and dashed lines are from the analytical model. The dot-dashed lines with points indicate the inclination profile of the disk at 17 Myr. Panel (g): The relationship between the equilibrium warp angle and system parameters. The x-axis shows the warp angle predicted by Equation \ref{['eq:trends']}. Blue dots represent results from the Laplace-Lagrange theory. The initial warp angles and those at 17 Myr of $N$-body simulation are presented using orange dots and circles.
• Figure 5: Comparison between observed and simulated images of the HD 110058 debris disk. Panel (a): The observed image from stasevic23, identical to the lowest panel of their Figure 5. The red solid line is the fitted spine of the disk. Panel (b): The synthetic image from our $N$-body simulation. Panel (c): Spatial distribution of particles in the $N$-body simulation. The colored solid lines mark the five rings in Figure \ref{['fig:fig2_lxly']} and \ref{['fig:warp_evo']}. Due to the forward scattering of the dust grains, we only plot the front side of the rings.