Radiation pressure mixing of large dust grains in protoplanetary disks

Dusty disks around young stars are formed out of interstellar dust that consists of amorphous, submicrometre grains. Yet the grains found in comets and meteorites, and traced in the spectra of young stars, include large crystalline grains that must have undergone annealing or condensation at temperatures in excess of 1,000 K, even though they are mixed with surrounding material that never experienced temperatures as high as that. This prompted theories of large-scale mixing capable of transporting thermally altered grains from the inner, hot part of accretion disks to outer, colder disk region, but all have assumptions that may be problematic. Here I report that infrared radiation arising from the dusty disk can loft grains bigger than one micrometre out of the inner disk, whereupon they are pushed outwards by stellar radiation pressure while gliding above the disk. Grains re-enter the disk at radii where it is too cold to produce sufficient infrared radiation pressure support for a given grain size and solid density. Properties of the observed disks suggest that this process might be active in almost all young stellar objects and young brown dwarfs.

💡 Research Summary

The paper addresses a long‑standing puzzle in protoplanetary disk (PPD) research: how large (≥1 µm) crystalline grains, which require temperatures above 1000 K to form, are found in the cold outer regions of disks where such temperatures never occur. Traditional large‑scale mixing mechanisms—turbulent diffusion, meridional circulation, and disk winds—each rely on assumptions that limit their applicability to these grains. Turbulent diffusion efficiently transports sub‑micron particles but loses effectiveness for larger grains because their stopping times become comparable to the orbital period. Meridional flows demand high viscosity (α ≈ 10⁻²) and steep temperature gradients that are not universally present. Disk winds can lift small dust but lack the momentum to carry millimetre‑scale particles outward.

The authors propose a two‑stage radiation‑pressure mechanism that operates entirely within the disk’s own radiation field and the stellar photon flux. First, infrared (IR) radiation emitted by the warm inner disk exerts an upward pressure on dust grains. By solving the vertical force balance that includes stellar gravity, gas pressure, and IR radiation pressure, they derive a critical grain size–density regime for which the IR pressure exceeds the combined downward forces. For typical silicate grains with densities ≤3 g cm⁻³, particles larger than about 1 µm become buoyant at radii of 0.1–1 AU. The IR pressure is proportional to the local IR luminosity (L_IR) and inversely to the square of the distance from the star, so the effect is strongest where the disk is hottest and most luminous.

Once lofted above the dense midplane, the grains encounter a low‑density environment where stellar radiation pressure dominates. The stellar flux provides a radially outward force characterized by the radiation‑pressure efficiency factor Q_pr. Because the grains are now essentially free‑floating above the gas, they are accelerated outward along the disk surface, reaching velocities of 10–30 km s⁻¹ in the simulations. The trajectory is largely a glide above the disk, minimizing gas drag and allowing the grains to travel several astronomical units before the IR radiation pressure that initially supported them diminishes.

The second crucial element is the “re‑entry radius.” As grains move outward, the disk temperature drops, reducing the IR radiation pressure. When the upward IR force falls below the sum of gravity and any residual gas drag, the grains lose buoyancy and fall back onto the disk surface. The authors compute this critical radius r_crit as a function of grain size, solid density, disk mass, and stellar luminosity. For a 2 µm silicate grain, r_crit is roughly 4 AU in a typical T Tauri disk, where the ambient temperature is below 150 K. At such temperatures, the crystalline structure formed in the hot inner disk is preserved, providing a natural explanation for the presence of high‑temperature minerals in cometary and meteoritic samples.

The paper validates the model against observational constraints. Infrared spectra of many young stellar objects display simultaneous signatures of amorphous sub‑micron dust and crystalline grains of several microns, a combination that is difficult to reconcile with pure turbulent mixing. The proposed mechanism reproduces the observed radial distribution of crystalline features without invoking unrealistically high turbulence levels. Moreover, the authors extend the analysis to disks around young brown dwarfs, whose lower luminosities still generate sufficient IR pressure to loft sub‑micron to micron‑scale grains, suggesting that the process may be universal among low‑mass objects.

Limitations are acknowledged. Grain–grain collisions could fragment lofted particles, altering the size distribution; gas drag, though reduced above the midplane, may still decelerate grains over longer timescales; and the radiation‑pressure efficiencies (absorption coefficients κ and Q_pr) depend on grain composition and morphology, which are not fully constrained experimentally. The authors propose future work incorporating full 3‑D radiative transfer coupled with magnetohydrodynamic (MHD) simulations to assess the interplay between radiation pressure and magnetic forces.

In summary, the study introduces a robust, physically grounded “buoyancy‑and‑radiation‑pressure” mixing pathway: infrared radiation from the inner disk lifts large grains, stellar radiation pressure drives them outward while they glide above the disk, and a temperature‑dependent loss of IR support causes them to re‑enter the disk at colder radii. This mechanism elegantly accounts for the widespread presence of large crystalline grains in protoplanetary disks and offers a compelling alternative to previously proposed mixing scenarios.