AC16 · Dust as a dynamical tracer

Debris Disk Collisions, Radiation and Planet Sculpting

What can a bright dust ring tell us about an unseen planet? A planet can shift dust orbits, open gaps and organize resonant populations. Collisions continually change grain sizes, and starlight sends smaller grains onto different paths. A bright arc also depends on temperature, emissivity and resolution, so its shape alone does not uniquely establish a planet.

Move a planet and follow an evolving multi-size dust release. The same weighted grains produce its thermal image, spectrum and mass ledger.

Elapsed model time $t=374\,\mathrm{yr}$

Resolved mass retained $M/M_0=0.896$

Resolved dust flux at 16.65 pc $F_\nu=1.68\times10^{-4}\,\mathrm{Jy}$

Mass ledger closure error $\delta_M=0$

Let dust leave different orbits

Planet semimajor axis

13.0\,\mathrm{AU}

a_p\;[\mathrm{AU}]

Planet eccentricity

0

e_p

Thermal observing wavelength

850\,\mu\mathrm m

\lambda\;[\mu\mathrm m]

Drag the gold planet grip to change its orbit. The existing dust then responds to the changed external force.

Saved synthetic release state, ready for playback.

Views, model limits and numerical ledger

A finite released population, not a steady disk. The 2:1 preset starts already occupied. Collisions, missing small fragments and open boundaries change mass; resampling changes higher moments.

Cumulative mapped kinetic change $-4.19\times10^{21}\,\mathrm{J}$

Largest raw compression jump $3.16\%$

Maximum normal optical depth $0.00194$

Read the derivations and rerun the real SMA/SED exercise

Model assumptions

Synthetic gas-poor, optically thin initial-value experiment. A finite prescribed release spectrum is not a derived steady collisional cascade. Weighted parcels carry physical grain counts in eight radius bins. Catastrophic collisions use a stated constant strength and fragment kernel; unresolved grains, inner and outer boundaries and numerical stops have separate mass reservoirs. Local velocity interpolation and capped resampling introduce measured distribution errors. The resonant preset starts with an occupied parent family; it does not simulate capture. A prescribed Kepler planet drives test grains, with the heliocentric indirect acceleration, absorption-only radiation pressure and physical-c Poynting–Robertson drag. External driving, collisions and PR do not conserve total grain mechanical energy. Grain absorption and equilibrium temperature generate both the SED and sky map. Body and grain marker sizes are schematic; stage colors indicate grain size, not measured optical scattering. The scene displays resolved grains above 1.5 micrometres; the complete thermal map includes every retained bin. Pixel maps carry Jy/pixel and the illustrative Gaussian smoothing does not convert them to Jy/beam. Actual SMA data and its measured beam appear only in the separate public tutorial. Gas drag, cratering, self-gravity, realistic mineral optical constants and unique planet detection are outside this model.

One grain record, four different questions

Synthetic thermal sky map, spectral flux and mass by size. Every panel uses the same retained physical population.

Thermal sky map · Jy/pixel

East left, north up; 45 degrees inclination. Each band uses its own peak scale. Image flux plus the stated field loss equals the same SED sum.

Spectral energy distribution

Resolved dust only at 16.65 pc; sink emission omitted. Logarithmic axes. White marker: current observing band. No observed flux is mixed into this synthetic record.

Where the resolved mass remains

Initial-mass fractions: outer domain

0.0609

;unresolved

0.0428

;inner boundary

0

;numerical stop

0

.Map field loss

0

Cumulative collision mass flow

Rows: destroyed radius. Columns: daughter radius, then unresolved sink. The same event depletes both partners; coarse bins can retain a diagonal contribution. Color shows relative cumulative mass on a logarithmic scale; exports retain kilograms.

Start with the essentials

Focus question: Can one grain record explain both the sky image and its spectrum?
One-sentence intuition: The physical grain populations drive the force, collisions and thermal sum. During a paused comparison, a wavelength change alters emission weights while orbital positions remain fixed.

Core mathematical model

Pressure competes with gravity

\beta=\frac{3L_\star\langle Q_{\rm pr}\rangle}{16\pi cGM_\star\rho s}

Absorption-only optical efficiency and a blackbody star define the model beta.

The same grains form every thermal band

F_\nu=\sum_i\frac{N_i\pi s_i^2Q_{\rm abs}(s_i,\nu)B_\nu(T_i)}{D^2}

Each term uses its actual retained position through radiative equilibrium.

Common difficulties

Libration is stronger evidence than a period ratio

Typical misconceptionA 2:1 period ratio proves resonance.

Better mental modelFollow the critical angle over many orbits and distinguish bounded libration from circulation. Initial occupation does not demonstrate capture.

A normalized beam does not change units

Typical misconceptionSmoothing Jy/pixel gives Jy/beam.

Better mental modelA unit conversion additionally needs the ratio of actual beam and pixel solid angles.

Run the experiment

01
Compare grain sizes
Load the resonant family and compare 70 micrometres with 850 micrometres. Follow the changed thermal weighting without replacing the physical record.
What to observe: The same positions receive different thermal weights.
02
Move the driver
Drag the gold planet marker. Observe the dust respond with time; its past positions do not change instantaneously.
What to observe: The future trajectory changes continuously under the new force.
03
Audit the mass
Select the denser release, then export the live record. Sum retained mass and every sink before comparing its image with the reference.
What to observe: Resolved mass can fall while retained plus all reservoirs still closes.