A multifluid approach for polydisperse pebble accretion: From particles to fluids, establishing the multifluid framework
A new simulation framework models planet formation from a realistic mix of particle sizes, finding that pebble accretion rates differ meaningfully from single-size assumptions.
The Thesis
Planet formation theorists have long assumed that the pebbles spiraling inward to feed a growing protoplanet are all roughly the same size — a simplification that makes the math tractable but doesn't match reality. This paper proposes a 'multifluid' framework, treating each size class of pebble as its own fluid layer in a 2D disk simulation, and finds that accretion rates shift noticeably when the full range of particle sizes is included. The effect matters most for high-mass planets, where the planet's own gravity disturbs the gas disk and changes which particle sizes get captured most efficiently. The catch is that this is a computational methods paper: it validates a framework against prior single-size results, not a finished predictive model of any specific planetary system. Real payoff comes in future work that uses this framework to model specific disk observations or constrain formation timelines.
Catalyst
The FARGO3D hydrodynamics code (a widely used open-source disk simulation tool) has matured to the point where researchers can modify it to handle multiple coupled dust-gas fluid species at tractable computational cost. Observational data from ALMA (the Atacama Large Millimeter Array) now routinely resolves protoplanetary disk substructure, creating demand for simulations sophisticated enough to compare against real disks with real particle-size distributions. The combination of available code infrastructure and richer observational targets made this class of simulation feasible now.
What's New
Most prior pebble accretion studies — including widely cited analytic models by Ormel & Klahr and Lambrechts & Johansen — assume a single representative pebble size, described by a single Stokes number (a dimensionless ratio of a particle's stopping time to the local orbital period, which governs how strongly it couples to the gas). Some later work introduced polydisperse (multi-size) distributions but used simplified particle-trajectory methods rather than fluid simulations. This paper instead runs each size bin as a full fluid species inside FARGO3D, capturing back-reaction (the friction that particles exert on the gas, not just the gas on particles) and disk perturbations from the planet simultaneously. The authors claim this approach handles high planet masses and tightly gas-coupled particles more accurately than prior trajectory-based methods.
The Counter
This paper is primarily a methods validation, and the bar it clears — reproducing prior single-size results in a new code framework — is not a high one. The authors acknowledge that their polydisperse-to-monodisperse accretion rate ratio differs from earlier estimates, but it's unclear whether this reflects a genuine physical insight or artifacts of their specific setup, such as the choice of a Mathis-Rumpl-Nordsieck (MRN) power-law size distribution originally developed for interstellar dust, not protoplanetary disks. The simulations are 2D, which cannot capture the vertical settling and stratification of particles that strongly influences accretion efficiency in real 3D disks. The parameter space explored is narrow, and no comparison to observed disk systems is attempted. Until this framework is applied to make predictions that can be tested against ALMA observations, it remains a more sophisticated computational tool looking for a problem to solve — not a result that changes our understanding of planet formation.
Longs
- No direct public market plays — this is foundational astrophysics research with no near-term commercial application
Shorts
- Researchers whose analytic single-size pebble accretion models are embedded in planet population synthesis codes may find their efficiency estimates need revision, particularly in the high-Stokes-number regime above St~0.3
Enablers (Picks & Shovels)
- FARGO3D open-source hydrodynamics code (Universidad de Santiago de Chile / CNRS) — the simulation engine modified in this work
- ALMA observatory — provides disk observations that motivate and constrain models like this
- HPC (high-performance computing) allocations at European research computing centers — this class of multi-fluid simulation is computationally intensive
Private Watchlist
- No private companies have a specific connection to this research
Resources
The Paper
Pebble accretion offers an efficient pathway to form planets, driven by a constant supply of inward drifting mass and an accretion efficiency enhanced by gas drag. While most studies assume a single pebble size (monodisperse), real discs contain a range of sizes (polydisperse), that drift, interact, and accrete at different rates. We aim to model polydisperse pebble accretion with a fluid approach, validating the method and exploring how gas disc evolution, solid-to-gas back-reaction, and a polydisperse size distribution affect growth. We use FARGO3D, modified to allow pebble accretion, to run 2D hydrodynamic simulations in a global disc with multiple dust/pebble species representing an underlying continuous pebble size distribution. With our framework of a multifluid approach, we have found values for pebble accretion efficiency consistent with earlier studies for a static gas disc. This confirms that our approach gives an accurate representation of pebble accretion. Evolving the gas disc we find lower efficiencies compared to an unperturbed gas disc for high Stokes number ($\gtrsim 0.3$) and higher efficiency for smaller St ($\lesssim0.3$). This effect is increased for higher planet masses. The accretion rate is mostly dominated by the highest Stokes numbers in our parameter study ($\mathrm{St}\in[10^{-2},10^0]$). The ratio we find between polydisperse and monodisperse pebble accretion rate is higher than previous estimations. We have constructed a multifluid model framework capable of accurately simulating polydisperse pebble accretion consistent with previous studies. This framework has advantages for simulating higher planet masses, as well as modelling multiple pebble species which are coupled to the gas. We find the perturbation of the protoplanet on the gas-disc lowering accretion rate when assuming a MRN-distribution of solids.