pith. sign in

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

Icy Volatile Enhancements in Evolving Protoplanetary Disks

Elizabeth Yunerman , Ellen Price , Karin \"Oberg This is my paper

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

classification 🌌 astro-ph.EP
keywords protoplanetary disksice linesvolatile icesplanet formationC/O ratioN/O ratiodisk dynamicssublimation
0
0 comments X

The pith

In evolving protoplanetary disks, particle drift and advection enhance hypervolatile ices relative to water by up to 100 times beyond their ice lines, producing solid C/O and N/O ratios near 1.

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

The paper establishes that dynamic processes in protoplanetary disks, including particle drift and outward advection, cause volatile ices to become enriched relative to H2O ice compared to what static models predict. Early on, drift removes water-rich material from the outer disk, and later, advection and deposition build up the enhancements further. This matters because the resulting compositions of solids influence what planets and planetesimals form from, particularly their carbon and nitrogen content. The enhancements are robust for hypervolatiles but vary for mid-volatiles depending on disk parameters.

Core claim

Incorporating additional carbon, nitrogen, and oxygen species along with more particle sizes into models of viscous disks with drift, thermal evolution, and desorption shows that before about 0.5 million years, the outer disk is desiccated by drift, enhancing relative volatile ices, while at later times outward advection and volatile deposition increase these enhancements further. The combined effect produces solid C/O and N/O ratios of approximately 1 beyond the hypervolatile ice lines, with hypervolatiles like N2, CO, and CH4 increasing by about 100 times across the parameter space, and mid-volatiles like CO2 and NH3 by 2 to 50 times. This demonstrates the necessity of coupling disk dynam

What carries the argument

The relative enhancement of icy volatiles to H2O through the interplay of particle drift desiccating the outer disk early and outward advection with deposition later in the disk's evolution.

If this is right

  • Hypervolatile ices such as CO, N2, and CH4 show robust enhancements of approximately 100 times beyond their ice lines.
  • Mid-volatile ices like CO2 and NH3 exhibit enhancements between 2 and 50 times, sensitive to specific model parameters.
  • Solid C/O and N/O ratios reach values near 1 in the outer disk regions, higher than in static disk models.
  • Advection plays a central role in redistributing volatiles across different disk radii.
  • The compositions of grains and planetesimals must account for these dynamic effects to match observed planetary properties.

Where Pith is reading between the lines

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

  • These dynamic enhancements could help explain observed carbon and nitrogen abundances in comets and outer solar system bodies.
  • Models that ignore disk evolution might underestimate the volatile content available for planet formation in the outer regions.
  • Testing with additional chemical species or different disk viscosities could reveal further variations in enhancement factors.
  • Observations targeting the radial distribution of ices in young disks at early and late stages would provide direct tests of the predicted time-dependent enhancements.

Load-bearing premise

The drift, advection, and desorption or adsorption processes from simpler models continue to control the volatile distributions even after adding more molecular species and particle size bins without introducing new dominant processes.

What would settle it

Spectroscopic measurements of the relative abundances of CO and H2O ices at radii beyond the CO ice line in a protoplanetary disk at an age greater than 0.5 million years that show no significant enhancement compared to static predictions would falsify the central claim.

Figures

Figures reproduced from arXiv: 2604.14124 by Elizabeth Yunerman, Ellen Price, Karin \"Oberg.

Figure 1
Figure 1. Figure 1: The disk midplane temperature may evolve signif￾icantly over the lifetime of a disk, in particular when incorpo￾rating an evolving stellar luminosity and mass accretion rate. Our model uses the MIST stellar evolution track for stellar luminosity with time, L⋆(t), of a 0.5M⊙ M dwarf shown in panel (a), and uses the relation M˙ ⋆(t) ∝ t −1 normalized to an average T Tauri star (red point) for an evolving mas… view at source ↗
Figure 2
Figure 2. Figure 2: Bulk gas surface density (light blue) decreases and viscously spreads outward over 3 Myr, while bulk solid surface density (dark blue) significantly decreases, with a small portion in the outer disk that increases. The line styles correspond to different times in the evolution with solid, dot￾ted, dotted-dashed, and dashed taken at initial (1 yr), 0.03, 0.3, and 3 Myr, respectively. For an accretion disk w… view at source ↗
Figure 3
Figure 3. Figure 3: Panel (a) displays the normalized mass distri￾butions for particles in a collisional cascade, showing how the bulk solid disk mass is partitioned across five particle size bins from smin = 0.1 µm to smax = 10 cm in our model. Varying the fiducial slope αs = 11/6 (blue) by ±5% skews the bulk of the solid mass either towards larger (red) or smaller (yellow) particles. Panel (b) then shows the evolved solid s… view at source ↗
Figure 4
Figure 4. Figure 4: Compositional ice lines can be numerically approximated where the total solid surface density for [PITH_FULL_IMAGE:figures/full_fig_p009_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Volatile vapor (Σg) and ice (Σs) surface density evolution varies significantly across different species and particle sizes, indicative of interconnected thermal, physical, and molecular processes setting the distri￾bution of planet forming materials. Solid and dashed lines are times of 10 kyr and 3 Myr respectively, while in between timesteps are more transparent. Each color, and row, corresponds to a par… view at source ↗
Figure 6
Figure 6. Figure 6: Analytic viscous (Eq. 9), drift (Eq. 19), adsorption (inverse of Eq. 22), and desorption (in￾verse of Eq. 21) timescales calculated using the fidu￾cial model bulk gas surface density and temperature at disk times of 0.03, 0.3, and 3 Myr. The drift and ad￾sorption timescales assume a particle size of 0.1 µm, tracing the dominant solid surface area. The species labeled lines correspond to desorption timescal… view at source ↗
Figure 7
Figure 7. Figure 7: Hypervolatile and mid-volatile ices can over time become ∼ 10−100× enhanced beyond their ice lines as compared to H2O ice, where the emerging enhanced region is advected according to the evolv￾ing critical radius. The left column is the radial–tem￾poral relative enhancement evolution as defined by Eq. 25, bounded by the species ice line where the color first appears. Contours of 2, 5, 10, 20, 50, and 100× … view at source ↗
Figure 8
Figure 8. Figure 8: Relative enhancement at 3 Myr of various volatile ices to H2O ice as compared with initial abundances. All species become relatively enhanced beyond ∼20 au for the fiducial disk model of an actively evolving disk where all species are inherited in their solid phase. efficient (re)adsorption for species whose ice line is near or beyond the critical disk radius at any point in the disk’s lifetime. 4.3. Compa… view at source ↗
Figure 9
Figure 9. Figure 9: Relative icy enhancements occur regard￾less of variation to physical disk parameters. Relative enhancement of CO2 and CO ice to H2O ice as compared with initial abundances at 3 Myr assuming different efficien￾cies for the physical dependencies. The solid line is the fidu￾cial model, where the drift efficiency is fd = 0.1, dust-to– gas ratio is 10−3 , and turbulence parameter is α = 10−3 . The dashed line i… view at source ↗
Figure 10
Figure 10. Figure 10: displays the 3 Myr CO2 and CO icy rel￾ative enhancements compared to H2O ice for particle 1 10 100 disk radius, r [au] 1 10 100 relative icy enhancement CO2 CO fiducial higher drift efficiency higher d2g lower alpha [PITH_FULL_IMAGE:figures/full_fig_p015_10.png] view at source ↗
Figure 12
Figure 12. Figure 12: As long as H2O ice is inherited, all species become significantly relatively enhanced as compared to H2O ice. Relative enhancement of CO2 and CO ice to H2O ice as compared with initial abundances at 3 Myr are shown for a variety of different maximum disk temperatures reached during the disk formation stage. In the fiducial cold case (solid), all species are initially frozen out with everything starting in… view at source ↗
Figure 13
Figure 13. Figure 13: Relative icy enhancements occur when the species ice line is near, beyond, or evolves through the [PITH_FULL_IMAGE:figures/full_fig_p019_13.png] view at source ↗
Figure 14
Figure 14. Figure 14: Outer disk regions become increasingly H [PITH_FULL_IMAGE:figures/full_fig_p020_14.png] view at source ↗
Figure 15
Figure 15. Figure 15: compares the evolving fiducial gas and solid volatile C/O and N/O ratios across the disk at different time steps. The evolving species’ ice lines and critical 1 10 100 0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 C/O (volatile) rc H2O CO2 CH4 CO GAS SOLIDS 1 10 100 disk radius, r [au] 10 2 10 1 10 0 10 1 10 2 10 3 10 4 N/O (volatile) NH+ 4 NH3 N2 GAS SOLIDS [PITH_FULL_IMAGE:figures/full_fig_p021_15.png] view at source ↗
Figure 16
Figure 16. Figure 16: Expected C/N/O ratios are sensitive to model assumptions. [PITH_FULL_IMAGE:figures/full_fig_p022_16.png] view at source ↗
read the original abstract

Protoplanetary disk ice lines shape a multitude of planet formation processes, setting the environmental composition through evolution. Ice line locations depend on molecular sublimation and deposition properties, but in dynamic disks where temperature and density structures change, so do the expected compositions of planets and planetesimals. In turbulent viscous disks with particle drift, thermal evolution, and desorption/adsorption, Price et al. 2021 demonstrated that the CO/H$_2$O ice ratio beyond the CO ice line can become enhanced by $\sim10\times$. We expand on their work by incorporating additional carbon, nitrogen, and oxygen species, more particle sizes, and a broader disk parameter exploration. We find that before $\sim0.5$Myr, volatile ices are enhanced relative to H$_2$O as the outer disk is desiccated by drift, while at later disk times outward advection and volatile deposition further increase relative volatile icy enhancements beyond the evolving critical disk radius. The outcome of these combined relative icy enhancement to H$_2$O mechanisms is solid C/O $\sim$ N/O $\sim1$ beyond the hypervolatile ice lines, much higher than expected in static disks. Hypervolatiles (N$_2$, CO, and CH$_4$) robustly increase to $\sim100\times$ across the explored parameter space, while mid-volatiles (CO$_2$ and NH$_3$) are sensitive to model choices, with enhancements ranging from $\sim2-50\times$. Together these results demonstrate that coupling disk dynamics with simple sublimation and deposition chemistry is fundamental to predicting grain, planetesimal, and planetary compositions, particularly the role of advection in redistributing volatiles across disk radii.

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 manuscript extends the Price et al. (2021) framework for volatile transport in turbulent viscous protoplanetary disks by adding CO2, NH3, N2, and CH4, increasing the number of particle size bins, and exploring a wider range of disk viscosity, turbulence, and initial abundance parameters. Time-dependent simulations of drift, advection, thermal evolution, and desorption/adsorption show that the outer disk is desiccated early on, followed by outward advection and deposition that enhance icy volatiles relative to H2O beyond the evolving critical radius. The central claim is that this produces solid C/O ∼ N/O ∼ 1 beyond hypervolatile ice lines, with hypervolatiles (N2, CO, CH4) enhanced by ∼100× and mid-volatiles (CO2, NH3) by 2–50× relative to static-disk expectations.

Significance. If the numerical results hold, the work demonstrates that coupling disk dynamics to simple sublimation/deposition chemistry is essential for predicting grain and planetesimal compositions, offering a dynamical explanation for elevated C/O and N/O ratios in solar-system bodies and exoplanets. The reported robustness of the hypervolatile enhancement across the explored parameter space is a clear strength, as is the explicit demonstration that advection redistributes volatiles after the early desiccation phase.

major comments (2)
  1. [Model extension and results (likely §2–3)] The robustness claim for the ∼100× hypervolatile enhancement rests on the assumption that adding CO2, NH3, N2, CH4 and extra particle sizes does not alter the dominant drift/advection/desorption pathways identified in Price et al. (2021). The abstract already notes that mid-volatile enhancements vary from 2–50× depending on model choices; a quantitative demonstration (e.g., a direct comparison run with and without the new species) that the critical-radius evolution and hypervolatile transport remain unchanged is needed to support the “robustly” qualifier.
  2. [Results and parameter exploration (likely §4)] Enhancement factors are reported as approximate ranges without accompanying error bars, standard deviations across runs, or convergence tests on particle-size resolution. Because the central quantitative claims are the 100× and 2–50× factors, the absence of these diagnostics makes it impossible to judge whether the reported values are numerically stable or sensitive to binning choices.
minor comments (2)
  1. [Abstract and §3] The time boundary “before ∼0.5 Myr” and “later disk times” should be tied to specific figures or tables so readers can map the two enhancement regimes to the plotted snapshots.
  2. [Methods] Notation for the evolving critical radius and the desorption/adsorption rates should be defined once in the methods and used consistently; occasional re-use of symbols from Price et al. (2021) without re-definition can confuse readers.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their constructive review and positive assessment of the work's significance. We address each major comment below, committing to revisions where they strengthen the manuscript.

read point-by-point responses

  1. Referee: The robustness claim for the ∼100× hypervolatile enhancement rests on the assumption that adding CO2, NH3, N2, CH4 and extra particle sizes does not alter the dominant drift/advection/desorption pathways identified in Price et al. (2021). The abstract already notes that mid-volatile enhancements vary from 2–50× depending on model choices; a quantitative demonstration (e.g., a direct comparison run with and without the new species) that the critical-radius evolution and hypervolatile transport remain unchanged is needed to support the “robustly” qualifier.

    Authors: We agree that a direct comparison would provide stronger support for the robustness qualifier. The hypervolatile species (N2, CO, CH4) are governed by the same low-sublimation-temperature physics as in Price et al. (2021), with transport dominated by drift and advection relative to the evolving critical radius set by disk thermal structure. Mid-volatile variability is already highlighted in the abstract. In revision we will add a side-by-side comparison run excluding the new mid-volatile species, confirming that hypervolatile enhancement factors and critical-radius evolution change by less than 5%. This will be presented in a new panel or appendix to justify the claim. revision: yes

  2. Referee: Enhancement factors are reported as approximate ranges without accompanying error bars, standard deviations across runs, or convergence tests on particle-size resolution. Because the central quantitative claims are the 100× and 2–50× factors, the absence of these diagnostics makes it impossible to judge whether the reported values are numerically stable or sensitive to binning choices.

    Authors: We thank the referee for identifying this presentational gap. The reported ranges reflect systematic variation across the explored disk parameters (viscosity, turbulence, initial abundances). In the revised manuscript we will add error bars to the enhancement plots showing the standard deviation across the model ensemble. We will also include a particle-size convergence test (in the methods or an appendix) demonstrating that increasing the number of bins beyond our fiducial choice alters the reported factors by less than 10%. These additions will allow quantitative assessment of numerical stability. revision: yes

Circularity Check

0 steps flagged

No significant circularity in simulation-derived enhancements

full rationale

The paper reports outcomes from forward time-dependent numerical simulations of disk evolution, particle drift, advection, and volatile deposition/adsorption, extended from a prior base model by adding species and size bins. No quoted equations or steps show outputs reducing by construction to fitted inputs, self-defined ratios, or unverified self-citation chains. The Price et al. 2021 reference supplies the starting framework but the new parameter explorations and reported enhancements (e.g., hypervolatile factors) are independently generated results, not tautological renamings or forced predictions.

Axiom & Free-Parameter Ledger

2 free parameters · 1 axioms · 0 invented entities

The model rests on standard protoplanetary-disk assumptions rather than new postulates; free parameters are typical choices for viscosity and initial conditions in such simulations.

free parameters (2)
  • disk viscosity and turbulence strength
    Chosen to explore evolutionary timescales and particle drift rates in the viscous disk framework.
  • initial molecular abundances
    Set to standard interstellar or disk values before evolution begins.
axioms (1)
  • domain assumption Turbulent viscous disk evolution with particle drift, thermal evolution, and desorption/adsorption
    Invoked as the physical framework inherited from Price et al. 2021 and extended here.

pith-pipeline@v0.9.0 · 5610 in / 1337 out tokens · 33105 ms · 2026-05-10T11:48:09.077409+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

3 extracted references · 3 canonical work pages

  1. [1]

    F., Feaga, L

    A’Hearn, M. F., Feaga, L. M., Keller, H. U., et al. 2012, ApJ, 758, 29, doi: 10.1088/0004-637X/758/1/29 Altwegg, K., Balsiger, H., & Fuselier, S. A. 2019, Annu. Rev. Astron. Astrophys., 57, 113, doi: 10.1146/annurev-astro-091918-104409 Altwegg, K., Balsiger, H., Hänni, N., et al. 2020, Nat Astron, 4, 533, doi: 10.1038/s41550-019-0991-9 Andrews, S. M. 2020...

    work page doi:10.1088/0004-637x/758/1/29 2012

  2. [2]

    adsabs.harvard.edu/abs/1973A%26A....24..337S/abstract Siess, L., Dufour, E., & Forestini, M

    https://ui. adsabs.harvard.edu/abs/1973A%26A....24..337S/abstract Siess, L., Dufour, E., & Forestini, M. 2000, Astronomy and Astrophysics, 358,

    work page 2000

  3. [3]

    https://ui.adsabs.harvard.edu/abs/2000A&A...358..593S Simon, A., Rajappan, M., & Öberg, K. I. 2023, The Astrophysical Journal, 955, 5, doi: 10.3847/1538-4357/aceaf8 Sirono, S.-i., & Kudo, D. 2021, The Astrophysical Journal, 911, 114, doi: 10.3847/1538-4357/abec7c Smith, R. S., May, R. A., & Kay, B. D. 2016, J. Phys. Chem. B, 120, 1979, doi: 10.1021/acs.jp...

    work page doi:10.3847/1538-4357/aceaf8 2023