pith. sign in

arxiv: 2604.14358 · v1 · submitted 2026-04-15 · 🌌 astro-ph.EP

Simulating the interplay between the snowline pebble flux and ongoing planet formation and migration

Danila Astrakhantsev , Sebastiaan Krijt , Sofia Savvidou , Bertram Bitsch This is my paper

Pith reviewed 2026-05-10 11:44 UTC · model grok-4.3

classification 🌌 astro-ph.EP
keywords planet formationpebble accretionsnowlineprotoplanetary disksplanet migrationpebble fluxdisk turbulence
0
0 comments X

The pith

The snowline pebble flux at protoplanet insertion strongly correlates with the final planet mass formed by pebble accretion.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper runs planet formation simulations that include pebble accretion and orbital migration to test how the pebble mass flux crossing the snowline shapes the masses and orbits of the planets that grow. A clear link appears between the flux value at the moment a seed protoplanet is placed and the mass the planet reaches by the end of the run. In high-turbulence disks the link is smooth, while in low-turbulence disks it becomes a sharp threshold: only fluxes above roughly 100 Earth masses per million years produce giant planets. The simulations also show that giant planets disturb the snowline flux for only a short time before they migrate inward, and that the modeled fluxes line up with those inferred from roughly one-million-year-old disks, especially larger ones with modest fragmentation speeds.

Core claim

Across the explored range of disk parameters, the snowline pebble mass flux recorded at the instant of protoplanet insertion sets the final planet mass. The relation is continuous when turbulence is high (alpha equal to 10 to the minus 3) but becomes a step function at lower turbulence (alpha equal to 10 to the minus 4), with giant planets appearing only when the initial flux exceeds 100 Earth masses per million years. Giant planets in the high-turbulence case alter the snowline flux for only about 10^5 to 10^6 years before they grow and cross the snowline. The resulting planet-containing models match observed pebble fluxes in disks of age about 1 Myr, particularly for disk sizes of 40 au or

What carries the argument

Numerical simulations of pebble accretion and migration that track both the growth of inserted protoplanets and their back-reaction on the radial pebble flux at the snowline.

If this is right

  • Pebble fluxes measured in young disks can be used directly to test or rule out families of planet-formation models.
  • Pebble-accretion growth remains consistent with the disk-evolution constraints now available from infrared observations.
  • Only disks larger than about 40 au, with turbulence near 10 to the minus 3 and low fragmentation speeds, produce planets whose masses and orbits match the observed sample.
  • A giant planet perturbs the snowline pebble supply for only a brief interval before it migrates inward and the flux recovers.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Higher-resolution maps of pebble flux across many disks could reveal whether real systems sit above or below the reported step-function threshold.
  • The temporary nature of the giant-planet perturbation suggests that time-resolved observations might catch disks in the act of recovering from such an event.
  • Repeating the runs with seeds inserted at earlier or later times would test how sensitive the mass-flux correlation is to the assumed formation epoch.

Load-bearing premise

The reported correlation and step-function threshold rest on the specific choices made for turbulence level, fragmentation velocity, overall disk size, and the precise moment when the seed protoplanet is inserted.

What would settle it

Discovery of a giant planet inside a disk whose measured snowline pebble flux at an age of about 1 Myr lies well below 100 Earth masses per million years would contradict the low-turbulence threshold.

Figures

Figures reproduced from arXiv: 2604.14358 by Bertram Bitsch, Danila Astrakhantsev, Sebastiaan Krijt, Sofia Savvidou.

Figure 1
Figure 1. Figure 1: Snowline pebble mass fluxes over time for disks of varying 𝑅0, 𝛼 and 𝑣frag. The shaded region shows pebble mass fluxes assumed in the accretion simulations of Lambrechts et al. (2019). then in turn results in faster dust depletion and lower pebble fluxes at 𝑡 ≳ 1Myr. 2.4 Protoplanet injection For a given disc model, we insert (proto)planets (one planet per disc) at different times 𝑡0 = [0.1, 0.25, 0.5, 0.7… view at source ↗
Figure 2
Figure 2. Figure 2: Planet growth tracks, mass vs semi-major axis. Pink dots represent the final planet mass and positions, with the numbered red dots showing selected planets numbered 1-4. Crosses show where the planet crosses the water snowline, if it does so. The black dashed line represents the approximate pebble isolation mass (Eq. 2) for a disk with 𝑅0 = 40 au. The difference to other disks is minor, within a few 𝑀⊕. Pl… view at source ↗
Figure 3
Figure 3. Figure 3: Final planet mass vs pebble flux across the snowline at time of insertion. Plots A and B are coloured by the planet’s initial and final semi-major axes respectively. Plot C is coloured by the planet’s insertion time 𝑡0, and plot D is coloured by the characteristic disk radius 𝑅0. Plot E is coloured by the disk’s 𝛼 parameter, and plot F is coloured by the fragmentation velocity 𝑣frag. All plots contain outl… view at source ↗
Figure 4
Figure 4. Figure 4: Final planet mass vs pebble flux across the snowline at 1 Myr. Plots A and B are coloured by the planet’s final and initial semi-major axes respectively. Plot C is coloured by the planet’s insertion time 𝑡0, and plot D is coloured by the characteristic disk radius 𝑅0. Plot E is coloured by the disk’s 𝛼 parameter, and plot F is coloured by the fragmentation velocity 𝑣frag. All plots contain outlined numbere… view at source ↗
Figure 5
Figure 5. Figure 5: Pebble fluxes over time for disks of varying characteristic radius 𝑅0 (left to right), turbulence level 𝛼 (blue vs. green) and fragmentation velocity 𝑣frag (top and bottom). The thicker lines show the flux from [PITH_FULL_IMAGE:figures/full_fig_p009_5.png] view at source ↗
read the original abstract

Pebble drift plays a central role in modern planet formation models. In this work we carry out planet formation simulations (including pebble accretion and migration) for a range of disc parameters to investigate (a) the impact of the snowline pebble mass flux on final planet orbits and masses, and (b) the back-reaction of growing and migrating planets on the snowline pebble fluxes in their natal discs. We find a strong correlation between the snowline pebble flux (at the time of protoplanet insertion) and the final planet mass. The correlation is continuous in disks with high turbulence levels ($\alpha=10^{-3}$), but exhibits a step function at lower turbulence ($\alpha=10^{-4}$), with giant planet formation requiring (initial) snowline pebble mass fluxes exceeding $100~\mathrm{M_\oplus Myr^{-1}}$. We find qualitative agreement between pebble mass fluxes inferred for discs aged ${\sim}1~\mathrm{Myr}$ and our planet-containing models, especially for larger disks ($\geq$40 au), high $\alpha$ ($10^{-3}$), and low $v_\mathrm{frag}$ ($3\mathrm{~m~s}^{-1}$). Additionally, giant planets in high turbulence disks are found to perturb the snowline pebble flux only temporarily (for ${\approx}10^{5-6}\mathrm{~yr}$) due to them quickly growing and migrating across the snowline. Our simulations show that currently observed pebble fluxes can indeed be used to constrain planet formation simulations, emphasizing that planet formation via pebble accretion is broadly in agreement with the currently available constraints from disc evolution as provided by JWST.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The paper reports results from numerical simulations of planet formation that include pebble accretion, migration, and the back-reaction of growing planets on the snowline pebble flux. It identifies a strong correlation between the snowline pebble mass flux measured at the moment of protoplanet insertion and the final planet mass. This correlation is continuous for high turbulence (α = 10^{-3}) but shows a step-function behavior at lower turbulence (α = 10^{-4}), with giant planets forming only when the initial flux exceeds ~100 M⊕ Myr^{-1}. The work also examines how giant planets temporarily perturb the pebble flux and finds qualitative consistency between simulated fluxes in ~1 Myr disks and current observations, particularly for larger disks, high α, and low v_frag.

Significance. If the reported correlations and turbulence-dependent thresholds prove robust, the manuscript provides a useful observational constraint on pebble-accretion models by linking measurable snowline pebble fluxes to planet masses and orbits. The explicit treatment of back-reaction on the flux and the parameter survey across disk sizes and turbulence levels are strengths. The qualitative match to JWST-era disk observations adds relevance, though the central claims rest on forward simulations whose sensitivity to insertion timing and disk parameters requires further quantification.

major comments (2)
  1. [Results] Results section (description of the step-function threshold): The reported discontinuity in final planet mass versus snowline pebble flux at α = 10^{-4} is defined using the instantaneous flux value at the (fixed) protoplanet insertion epoch. Because insertion time simultaneously controls both the measured flux and the remaining growth/migration window, the step function could be an artifact of the chosen insertion timing rather than a general outcome. The manuscript does not report tests in which insertion time is varied while holding other parameters fixed, nor does it integrate the time-dependent flux history prior to insertion.
  2. [Methods] Methods section (numerical setup): The abstract and results present outcomes from a suite of simulations but provide no information on spatial resolution, time-stepping criteria, convergence tests, or how the pebble flux is computed at the snowline. Without these details it is impossible to assess whether the continuous-versus-step-function distinction is numerically converged or sensitive to the specific implementation of pebble drift and back-reaction.
minor comments (2)
  1. [Abstract] The abstract states 'qualitative agreement' with observed pebble fluxes but does not specify the quantitative metric or the exact observational sample used for the comparison.
  2. [Figures] Figure captions and text should explicitly state the range of disk outer radii explored and whether the reported fluxes are azimuthally averaged or taken at a single radial location.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful and constructive review. The comments have prompted us to strengthen the presentation of our numerical methods and to test the robustness of the reported step-function behavior. We respond to each major comment below and have revised the manuscript accordingly.

read point-by-point responses

  1. Referee: [Results] Results section (description of the step-function threshold): The reported discontinuity in final planet mass versus snowline pebble flux at α = 10^{-4} is defined using the instantaneous flux value at the (fixed) protoplanet insertion epoch. Because insertion time simultaneously controls both the measured flux and the remaining growth/migration window, the step function could be an artifact of the chosen insertion timing rather than a general outcome. The manuscript does not report tests in which insertion time is varied while holding other parameters fixed, nor does it integrate the time-dependent flux history prior to insertion.

    Authors: We acknowledge that fixing the insertion time at 0.5 Myr could in principle couple the measured flux to the available growth time. This epoch was selected because it corresponds to the typical onset of efficient pebble accretion in the inner disk in standard models. While the primary simulation suite uses a single insertion time, the step-function threshold and overall correlation are recovered across a broad range of disk sizes and fragmentation velocities, which themselves modulate the time evolution of the pebble flux. To directly address the concern, we have carried out additional runs with insertion times shifted to 0.3 Myr and 0.7 Myr for representative low-turbulence cases. The flux threshold for giant-planet formation remains near 100 M⊕ Myr^{-1} and the discontinuous behavior persists. These tests have been added to Section 3.2 together with a short discussion. We have also included a supplementary figure showing the cumulative pebble flux integrated up to insertion for the key models, as suggested. revision: yes

  2. Referee: [Methods] Methods section (numerical setup): The abstract and results present outcomes from a suite of simulations but provide no information on spatial resolution, time-stepping criteria, convergence tests, or how the pebble flux is computed at the snowline. Without these details it is impossible to assess whether the continuous-versus-step-function distinction is numerically converged or sensitive to the specific implementation of pebble drift and back-reaction.

    Authors: We apologize for the omission of these technical details in the original submission. The simulations employ a 1D radial grid with 250 logarithmically spaced cells extending from 0.1 au to 150 au. Time steps are adaptive and limited by the minimum of the viscous, migration, and accretion timescales (with a hard cap of 100 yr). Convergence was verified by repeating a subset of models at 500 cells; final planet masses and semi-major axes agree to within 10 %. The snowline pebble flux is obtained each time step by integrating the product of pebble surface density and inward radial velocity across the snowline location, with the gas and dust profiles updated to include the torque and gap-opening effects of the growing planet. A new subsection (2.3) has been added to the Methods that fully documents the grid, time-stepping algorithm, convergence tests, and the precise definition of the snowline flux. revision: yes

Circularity Check

0 steps flagged

No significant circularity; results emerge from forward simulations.

full rationale

The paper reports outcomes from numerical planet-formation simulations that include pebble accretion, migration, and back-reaction on the pebble flux. The claimed correlation between snowline pebble flux (measured at the fixed insertion epoch) and final planet mass, together with the turbulence-dependent step-function behavior, are computed results rather than quantities fitted to observations or defined in terms of themselves. No load-bearing step reduces by construction to a prior self-citation, an ansatz smuggled via citation, or a renaming of a known empirical pattern; the simulations are self-contained against the stated disk parameters and evolve the system forward from initial conditions.

Axiom & Free-Parameter Ledger

3 free parameters · 1 axioms · 0 invented entities

The work rests on standard assumptions of pebble accretion and type-I/II migration in viscous disks; specific free parameters include turbulence strength α, pebble fragmentation velocity, and disk radial extent, all chosen to explore parameter space rather than derived from first principles.

free parameters (3)
  • turbulence parameter α
    Set to 10^{-3} or 10^{-4} to explore high and low turbulence regimes; directly controls pebble diffusion and accretion efficiency.
  • pebble fragmentation velocity v_frag
    Set to 3 m/s in some runs; controls maximum pebble size and thus drift speed.
  • disk outer radius
    Varied around 40 au and larger; affects total pebble reservoir and flux at snowline.
axioms (1)
  • domain assumption Pebble accretion and migration prescriptions are valid in the chosen disk models
    Standard in the field but not re-derived here.

pith-pipeline@v0.9.0 · 5610 in / 1468 out tokens · 26658 ms · 2026-05-10T11:44:23.496889+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

1 extracted references · 1 canonical work pages

  1. [1]

    Andama G., Mah J., Bitsch B., 2024, A&A, 683, A118 Asplund M., Grevesse N., Sauval A. J., Scott P., 2009, ARA&A, 47, 481 Banzatti A., et al., 2020, ApJ, 903, 124 Banzatti A., et al., 2023, ApJ, 957, L22 Benítez-Llambay P., Masset F., Koenigsberger G., Szulágyi J., 2015, Nature, 520, 63 Birnstiel T., 2024, ARA&A, 62, 157 Birnstiel T., Klahr H., Ercolano B....

    work page doi:10.48550/arxiv.2203.09759 2024