pith. machine review for the scientific record. sign in

arxiv: 2604.05178 · v1 · submitted 2026-04-06 · 🌌 astro-ph.EP · astro-ph.SR

Recognition: 2 theorem links

· Lean Theorem

Gas Giants Formed by Gravitational Instability May Accrete Atmospheres with Super-Stellar Carbon to Oxygen Ratios

Alan P. Boss
Authors on Pith no claims yet

Pith reviewed 2026-05-10 18:51 UTC · model grok-4.3

classification 🌌 astro-ph.EP astro-ph.SR
keywords exoplanet atmospheresgravitational instabilityprotoplanetary diskscarbon-to-oxygen ratiogas giant formationice evaporationsolid migration
0
0 comments X

The pith

Gas giants formed by disk gravitational instability can accrete atmospheres with super-stellar carbon-to-oxygen ratios.

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

Standard models assume that planets formed by gravitational instability in a gas disk will have atmospheric carbon-to-oxygen ratios identical to their host stars. This paper demonstrates that the process can instead produce elevated C/O ratios when the disk midplane sits in the temperature range where water ice has evaporated but carbon dioxide ice has not. Oxygen locked in solids can then be removed from the planet-forming region faster than the gas disperses, leaving behind carbon-enriched material for the protoplanet to accrete. The outcome holds only when the final orbit lies outside roughly 7 astronomical units from a solar-mass star; closer orbits yield near-stellar C/O. This finding revises how atmospheric composition can be used to diagnose formation pathways for gas giants.

Core claim

Analysis of fully three-dimensional radiative hydrodynamics models shows that GDGI protoplanets can form and accrete disk gas with super-stellar C/O ratios when midplane temperatures lie between the evaporation temperatures of water ice (~135 K) and CO2 ice (~47 K). Solids that avoid fragmentation and grow rapidly to centimeter sizes can be transported inward or outward faster than the disk gas dissipates, resulting in preferential accretion of carbon-rich gas over oxygen-rich ices, provided the protoplanet's orbit remains outside ~7 au from a solar-mass protostar. Orbits inside ~7 au lead to accretion of gas with nearly stellar C/O.

What carries the argument

The disk midplane temperature window between water-ice and CO2-ice evaporation temperatures, which enables differential removal of oxygen-rich solids while carbon-bearing gas remains available for accretion.

If this is right

  • GDGI planets that remain outside ~7 au are predicted to exhibit super-stellar atmospheric C/O ratios.
  • GDGI planets that end up inside ~7 au are predicted to exhibit near-stellar atmospheric C/O ratios.
  • Atmospheric C/O measurements can distinguish GDGI formation from core-accretion pathways in certain orbital regimes.
  • The long-standing assumption that all GDGI planets mirror host-star C/O ratios must be qualified by orbital distance.

Where Pith is reading between the lines

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

  • This selective-accretion process could produce high C/O signatures in wide-orbit, directly imaged planets without requiring core accretion or migration.
  • Similar temperature-dependent removal of solids might affect other volatile ratios if additional condensation fronts are present in the disk.
  • Atmospheric observations of young gas giants at 10-50 au could provide a direct test of whether gravitational instability is a common formation channel.

Load-bearing premise

Centimeter-sized solids must grow rapidly without fragmenting and must be carried away from the planet's feeding zone much faster than the disk gas dissipates.

What would settle it

A spectroscopic measurement of stellar C/O in the atmosphere of a confirmed gas giant planet orbiting beyond 7 au would contradict the predicted super-stellar enrichment.

Figures

Figures reproduced from arXiv: 2604.05178 by Alan P. Boss.

Figure 10
Figure 10. Figure 10 [PITH_FULL_IMAGE:figures/full_fig_p012_10.png] view at source ↗
Figure 1
Figure 1. Figure 1: — Log midplane density (g cm−3 ) contours for model fldAC after 196 yr of evolution. Numerous dense spiral arms and the four virtual protoplanets (VPs) are evident as dark red features in the inner disk. The disk has an inner radius of 4 au and an outer radius of 20 au [PITH_FULL_IMAGE:figures/full_fig_p017_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: — Log midplane gas temperature (K) contours for model fldAC [PITH_FULL_IMAGE:figures/full_fig_p018_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: — Midplane temperature profile along an azimuth at 9 o’clock in Figu [PITH_FULL_IMAGE:figures/full_fig_p019_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: — Vertical optical depth from the disk surface to the disk midp [PITH_FULL_IMAGE:figures/full_fig_p020_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: — Temperatures of the local disk gas accreted by VP2 during t [PITH_FULL_IMAGE:figures/full_fig_p021_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: — Radial profile of the Stokes number of 1-cm-size particles or [PITH_FULL_IMAGE:figures/full_fig_p022_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: — Same as Figure 6 except at the later time of 196 yrs of evolutio [PITH_FULL_IMAGE:figures/full_fig_p023_7.png] view at source ↗
read the original abstract

Characterizing the atmospheric compositions of exoplanets, along with determining properties such as their mass, mean density, and orbital configuration, is thought to be an effective means for differentiating between various formation and evolution scenarios. Exoplanet atmospheric C/O ratios, when compared to host star C/O ratios, have been advanced as discriminators of gas giant formation and evolution scenarios in the context of the core accretion mechanism. Gas giants formed by gas disk gravitational instability (GDGI), on the other hand, are thought to have atmospheres with C/O ratios identical to their host stars. We examine this assumption through analysis of fully three dimensional radiative hydrodynamics models of the GDGI in the flux-limited diffusion approximation. We show here that GDGI protoplanets may be able to form and accrete disk gas with super-stellar C/O ratios, as a result of their formation and orbital evolution in a disk with midplane temperatures in the range of the evaporation temperatures of water ($\sim$ 135 K) and CO$_2$ ($\sim$ 47 K) ices. Solids that avoid fragmentation and grow rapidly to cm-size could be transported inward to the central protostar or outward to the edge of the disk considerably faster than the disk gas is dissipated, leading to the preferential accretion of C-rich disk gas compared to the O-rich ices, provided that the protoplanet's orbit remains outside $\sim$ 7 au from a solar-mass protostar. Orbits inside $\sim$ 7 au, however, could result in the accretion of disk gas with nearly stellar C/O.

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 claims that gas giants formed by gas disk gravitational instability (GDGI) can accrete atmospheres with super-stellar C/O ratios. This follows from analysis of 3D radiative hydrodynamics simulations in the flux-limited diffusion approximation, which show protoplanet formation and orbital evolution in disks with midplane temperatures between the evaporation temperatures of water (~135 K) and CO2 (~47 K) ices. The key mechanism is preferential removal of O-rich solids (via rapid growth to cm sizes and radial transport) faster than gas dissipation, leaving C-rich gas for accretion if the orbit remains outside ~7 au from a solar-mass star; inside this radius, near-stellar C/O is expected.

Significance. If the result holds, it would undermine the standard assumption that GDGI-formed gas giants must have stellar C/O ratios, providing a potential observational discriminator between core accretion and gravitational instability pathways. The use of fully 3D RHD models to capture self-consistent gas dynamics and orbital evolution is a positive feature, as is the focus on a specific temperature range tied to ice evaporation.

major comments (2)
  1. [Methods and model analysis] The 3D flux-limited diffusion RHD models (as described in the methods and model analysis sections) address only gas dynamics, clump formation, and orbital evolution. They contain no dust growth, radial drift, fragmentation, or volatile chemistry. The central claim that GDGI protoplanets accrete super-stellar C/O gas rests on the untested premise that cm-sized solids (primarily O-rich water ice) undergo rapid radial drift or outward transport faster than disk gas dissipates while the protoplanet remains outside ~7 au. This assumption is load-bearing and not demonstrated by the simulations.
  2. [Discussion of orbital evolution and results] The ~7 au orbital threshold separating super-stellar from stellar C/O accretion is stated without derivation, sensitivity tests, or reference to specific outputs from the RHD simulations showing how this radius is determined.
minor comments (2)
  1. [Abstract] The abstract and main text provide no quantitative results from the simulations, such as achieved C/O values, error estimates, or details on the number and parameter space of runs, limiting assessment of the effect's magnitude.
  2. [General presentation] The manuscript would benefit from an explicit statement distinguishing results directly obtained from the RHD simulations versus the additional assumptions about solid transport and ice depletion.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their constructive and detailed comments on our manuscript. We address each major comment below and indicate planned revisions to strengthen the presentation of our results.

read point-by-point responses

  1. Referee: [Methods and model analysis] The 3D flux-limited diffusion RHD models (as described in the methods and model analysis sections) address only gas dynamics, clump formation, and orbital evolution. They contain no dust growth, radial drift, fragmentation, or volatile chemistry. The central claim that GDGI protoplanets accrete super-stellar C/O gas rests on the untested premise that cm-sized solids (primarily O-rich water ice) undergo rapid radial drift or outward transport faster than disk gas dissipates while the protoplanet remains outside ~7 au. This assumption is load-bearing and not demonstrated by the simulations.

    Authors: We agree that the simulations are limited to gas dynamics and do not include dust growth, radial drift, fragmentation, or volatile chemistry. The central claim therefore relies on the assumption that O-rich cm-sized solids are removed faster than the gas dissipates, which is not directly simulated here. This assumption draws from established literature on particle growth and radial drift timescales in protoplanetary disks. In the revised manuscript we will expand the discussion to include explicit timescale comparisons between dust radial transport and gas dissipation, supported by additional citations to relevant disk evolution studies, and will clarify that the preferential removal of O-rich solids is a theoretical expectation based on prior work rather than a direct result of the RHD runs. revision: partial

  2. Referee: [Discussion of orbital evolution and results] The ~7 au orbital threshold separating super-stellar from stellar C/O accretion is stated without derivation, sensitivity tests, or reference to specific outputs from the RHD simulations showing how this radius is determined.

    Authors: The ~7 au value is an approximate radius at which protoplanet orbits in our simulations remain outside the region where midplane temperatures would allow near-stellar C/O accretion for a solar-mass star. We will add an explicit derivation in the revised manuscript, referencing specific simulation outputs that track orbital evolution and linking the threshold to the adopted disk temperature profile between the water and CO2 ice lines. We will also include a qualitative discussion of sensitivity to this radius. revision: yes

Circularity Check

0 steps flagged

No circularity; simulations and qualitative argument are independent of target C/O

full rationale

The paper's derivation rests on 3D flux-limited diffusion RHD simulations that compute gas dynamics, clump formation, and orbital evolution for disks whose midplane temperatures fall between the evaporation points of water and CO2 ices. These simulations contain no dust, no chemistry, and no C/O tracking, so the reported super-stellar C/O outcome is not produced by any equation inside the models. The additional claim that cm-sized O-rich solids can drift or be transported faster than the gas dissipates is presented as a separate physical consideration drawn from standard radial-drift timescales; it is not fitted to the simulation output, not defined in terms of the target C/O ratio, and not justified by a self-citation chain that would render the result tautological. No uniqueness theorem, ansatz, or parameter fit is smuggled in via prior work by the same author. The logic therefore does not reduce to its inputs by construction.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The paper builds on standard protoplanetary disk physics and hydrodynamics without introducing new free parameters or entities beyond the model assumptions.

axioms (2)
  • domain assumption Midplane temperatures of ~135 K for water ice and ~47 K for CO2 ice evaporation in protoplanetary disks
    These are standard values from astrophysical literature on ice sublimation.
  • domain assumption Solids can grow to cm sizes without fragmenting and be transported faster than gas dissipates
    This is a key assumption for the preferential accretion mechanism.

pith-pipeline@v0.9.0 · 5592 in / 1366 out tokens · 52995 ms · 2026-05-10T18:51:36.776029+00:00 · methodology

discussion (0)

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

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Forward citations

Cited by 2 Pith papers

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. JWST high-contrast spectroscopy with speckle modelling: Atmospheric retrievals of the T dwarf companion HD 19467 B

    astro-ph.EP 2026-05 unverdicted novelty 7.0

    Joint speckle-atmosphere retrievals on JWST spectra of HD 19467 B and a field T dwarf detect H2O, CH4, CO, CO2 and NH3, yield carbon isotopic ratios of 154 and 85 respectively, and show near-solar metallicity with sub...

  2. HAT-P-70b through the Eyes of MAROON-X: Constraining Elemental Abundances of Metals and Insights on Atmosphere Dynamics

    astro-ph.EP 2026-05 conditional novelty 6.0

    New MAROON-X observations of HAT-P-70b detect multiple neutral and ionized metals with day-to-night wind signatures and demonstrate that ionization-aware retrievals yield abundance ratios closer to solar values except...

Reference graph

Works this paper leans on

3 extracted references · cited by 2 Pith papers

  1. [1]

    20 25, A&A, 700, A105 Birnstiel, T., Ormel, C

    Bachmann, N., Kreidberg, L., Molli` ere, P., Deming, D., & Tsai, S.-M. 20 25, A&A, 700, A105 Birnstiel, T., Ormel, C. W., & Dullemond, C. P. 2011, A&A, 525, A11 Boley, A. C., Helled, R., & Payne, M. J. 2011, ApJ, 735, 30 Boss, A. P. 1981, ApJ, 244, 40 Boss, A. P. 1984, ApJ, 277, 768 Boss, A. P. 1997, Science, 276, 1836 Boss, A. P. 2005, ApJ, 629, 535 Boss...

    2011

  2. [2]

    Gas giant protoplanet properties at the termination of mod el fldAC (a continuation of model fldA of A. P. Boss 2021): VP number, mass (J upiter masses), semimajor axis (au), orbital eccentricity, local VP disk temperatu re (K). VP # Mp a e T V P VP1 0.48 5.9 0.28 250 VP2 0.64 6.6 0.087 133 VP3 1.3 4.8 0.16 275 VP4 0.27 5.7 0.32 150 – 17 – Fig. 1.— Log midp...

    2021

  3. [3]

    3.— Midplane temperature profile along an azimuth at 9 o’clock in Figu res 1 and 2 for model fldAC

    – 19 – Fig. 3.— Midplane temperature profile along an azimuth at 9 o’clock in Figu res 1 and 2 for model fldAC. Snowlines for C, H 2O, CO 2 and CO are shown. The profile passes through regions of different gas C/O ratios based on the disk gas temperatu res (as in K. I. ¨Oberg et al. 2011): C/O ≈

    2011