Envelopes of embedded super-Earths II. Three-dimensional isothermal simulations

Roman R. Rafikov; William B\'ethune

arxiv: 1907.02763 · v2 · pith:UUORFGLLnew · submitted 2019-07-05 · 🌌 astro-ph.EP

Envelopes of embedded super-Earths II. Three-dimensional isothermal simulations

William B\'ethune , Roman R. Rafikov This is my paper

Pith reviewed 2026-05-25 02:11 UTC · model grok-4.3

classification 🌌 astro-ph.EP

keywords super-Earthsprotoplanetary discsplanetary envelopeshydrodynamic simulationsBondi radiuscore accretionatmospheric recycling

0 comments

The pith

Envelope properties around planetary cores depend mainly on the ratio of the Bondi radius to the core's physical size.

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

The paper performs three-dimensional inviscid isothermal hydrodynamic simulations of gaseous envelopes forming around planetary cores embedded in protoplanetary discs. It examines how envelope properties change as core mass increases from sub-thermal to super-thermal values. The simulations reveal that global features such as rotational support and turbulent mixing depend primarily on the ratio of the Bondi radius to the core radius. For higher-mass cores, the flow pattern includes supersonic polar inflows that shock the envelope, equatorial outflows creating circulation, and a shocked shell separating regions of slow and fast recycling from disc gas. Inner envelope regions may remain bound despite outer mixing.

Core claim

Global envelope properties such as the amount of rotational support or turbulent mixing are mostly sensitive to the ratio of the Bondi radius of the core to its physical size. High-mass cores are fed by supersonic inflows arriving along the polar axis and shocking on the densest parts of the envelope, driving turbulence and mass accretion. Gas flows out of the core's Hill sphere in the equatorial plane, describing a global mass circulation through the envelope. The shell of shocked gas atop the core surface delimits regions of slow (inside) and fast (outside) material recycling by gas from the surrounding disc. While recycling hinders the runaway growth towards gas giants, the inner regions,

What carries the argument

The ratio of the Bondi radius of the core to its physical size, which controls the sensitivity of envelope properties such as rotational support and turbulent mixing.

If this is right

High-mass cores receive supersonic inflows along the polar axis that shock the envelope and drive turbulence and accretion.
Gas outflows in the equatorial plane create a global mass circulation through the envelope.
A shell of shocked gas separates inner slow-recycling regions from outer fast-recycling regions.
Recycling from the disc hinders runaway gas giant growth but inner envelope regions may stay bound.

Where Pith is reading between the lines

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

The same ratio sensitivity could hold in simulations with added viscosity if the Bondi-to-core scale remains dominant over viscous length scales.
Atmospheric observations of rotation or turbulence levels might indirectly constrain core radii for close-in super-Earths.
Extending the setup to radiative transfer could test whether the reported flow patterns persist when energy transport is included.

Load-bearing premise

The hydrodynamic equations are solved under the assumptions of inviscid flow and an isothermal equation of state throughout the domain.

What would settle it

A simulation including viscosity or a non-isothermal equation of state that shows envelope rotational support or mixing changing independently of the Bondi-to-core radius ratio.

read the original abstract

Massive planetary cores embedded in protoplanetary discs are believed to accrete extended atmospheres, providing a pathway to forming gas giants and gas-rich super-Earths. The properties of these atmospheres strongly depend on the nature of the coupling between the atmosphere and the surrounding disc. We examine the formation of gaseous envelopes around massive planetary cores via three-dimensional inviscid and isothermal hydrodynamic simulations. We focus the changes in the envelope properties as the core mass varies from low (sub-thermal) to high (super-thermal) values, a regime relevant to close-in super-Earths. We show that global envelope properties such as the amount of rotational support or turbulent mixing are mostly sensitive to the ratio of the Bondi radius of the core to its physical size. High-mass cores are fed by supersonic inflows arriving along the polar axis and shocking on the densest parts of the envelope, driving turbulence and mass accretion. Gas flows out of the core's Hill sphere in the equatorial plane, describing a global mass circulation through the envelope. The shell of shocked gas atop the core surface delimits regions of slow (inside) and fast (outside) material recycling by gas from the surrounding disc. While recycling hinders the runaway growth towards gas giants, the inner regions of protoplanetary atmospheres, more immune to mixing, may remain bound to the planet.

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.

Desk Editor's Note private letter to a colleague

The 3D runs map a clean geometric picture of polar inflows and equatorial outflows controlled by the Bondi-to-core radius ratio, but the inviscid isothermal setup leaves the robustness open.

read the letter

The main thing to know is that these simulations show global envelope properties like rotational support and mixing are mostly set by the ratio of Bondi radius to core size, with high-mass cores fed by supersonic polar inflows that shock and drive turbulence while gas leaves through the equator. A shocked shell then separates slow inner recycling from faster outer exchange with the disc. This organizes how retention versus recycling scales with core mass in a way that connects to population synthesis work on super-Earths and the gas giant transition. What is new is the three-dimensional treatment across sub-thermal to super-thermal cores that identifies that single ratio as the dominant control, going past the 2D or lower-resolution cases referenced in the abstract. The flow geometry is laid out clearly enough to be useful as a reference picture. The setup is stated plainly as inviscid Euler equations with an isothermal equation of state everywhere, which keeps the problem tractable and isolates the hydrodynamic circulation. That choice is standard for this class of problem and lets the authors focus on the mass flow patterns without added complications. The soft spot is the lack of any reported resolution studies, convergence tests, or direct comparison to analytic limits such as the expected Bondi accretion rate. Without those, the strength of the turbulence and the exact location of the shocked shell remain hard to judge numerically. The stress-test point about viscosity or non-isothermal thermodynamics is also on target: either addition could dissipate shocks or shift the effective sound speed, which might weaken or remove the claimed sensitivity to the radius ratio. The paper does not explore those extensions, so the result stays tied to the idealised case. This work is aimed at people who model envelope accretion in protoplanetary discs or who need hydrodynamic benchmarks for close-in planet formation. Readers working on recycling and retention in population synthesis codes will find the geometric description directly usable. It is coherent on its own terms and shows clear engagement with the relevant flow physics, so it deserves a serious referee even though the idealisations limit how far the conclusions can be taken without further checks.

Referee Report

2 major / 1 minor

Summary. The paper presents three-dimensional inviscid isothermal hydrodynamic simulations of gaseous envelope formation around planetary cores embedded in protoplanetary discs, varying core mass across sub-thermal to super-thermal regimes relevant to close-in super-Earths. It claims that global envelope properties such as rotational support and turbulent mixing are mostly sensitive to the ratio of the Bondi radius to the core's physical size, with high-mass cores exhibiting supersonic polar inflows that shock and drive turbulence, equatorial outflows, and a global mass circulation pattern in which recycling hinders runaway growth while inner regions may remain bound.

Significance. If the reported sensitivity to the Bondi-to-core radius ratio holds within the simulated regime, the work identifies a controlling parameter for envelope-disc coupling and distinguishes regions of fast versus slow recycling, offering a dynamical explanation for why some super-Earths retain bound atmospheres. The 3D global simulations capture circulation features inaccessible to lower-dimensional models.

major comments (2)

[Abstract] Abstract: the claim that global properties are 'mostly sensitive' to the Bondi/core radius ratio rests on qualitative outcomes from the parameter exploration; no quantitative diagnostics (e.g., fractional variance explained by the ratio versus other parameters, or explicit sensitivity coefficients) are supplied to support the adverb 'mostly'.
[Abstract] Abstract (numerical setup): the simulations solve the inviscid Euler equations with a strictly isothermal EOS; because the headline sensitivity is diagnosed from shock-driven inflows and shear-driven turbulence, the absence of any resolution study, convergence test, or comparison to analytic limits leaves the robustness of the diagnosed global properties (rotational support, turbulent mixing) unquantified and load-bearing for the central claim.

minor comments (1)

[Abstract] Abstract: 'We focus the changes' is grammatically incomplete and should read 'We focus on the changes'.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their thoughtful review and constructive suggestions. We address each major comment below and indicate the revisions we will make to the manuscript.

read point-by-point responses

Referee: [Abstract] Abstract: the claim that global properties are 'mostly sensitive' to the Bondi/core radius ratio rests on qualitative outcomes from the parameter exploration; no quantitative diagnostics (e.g., fractional variance explained by the ratio versus other parameters, or explicit sensitivity coefficients) are supplied to support the adverb 'mostly'.

Authors: The parameter exploration in our study varies the core mass while keeping other parameters fixed, making the Bondi-to-core radius ratio the dominant varying parameter. The observed trends in envelope properties are indeed qualitative, as no formal variance analysis was performed. We will revise the abstract to replace 'mostly sensitive' with 'primarily depend on' to more accurately reflect the nature of our findings. revision: partial
Referee: [Abstract] Abstract (numerical setup): the simulations solve the inviscid Euler equations with a strictly isothermal EOS; because the headline sensitivity is diagnosed from shock-driven inflows and shear-driven turbulence, the absence of any resolution study, convergence test, or comparison to analytic limits leaves the robustness of the diagnosed global properties (rotational support, turbulent mixing) unquantified and load-bearing for the central claim.

Authors: We recognize the importance of demonstrating numerical robustness, particularly for properties influenced by shocks and turbulence in inviscid simulations. In the revised version, we will add a resolution study comparing results at different grid resolutions to confirm convergence of the key global properties. revision: yes

Circularity Check

0 steps flagged

No circularity: results from direct numerical simulation parameter sweep

full rationale

The paper performs 3D inviscid isothermal hydrodynamic simulations across a range of core masses, diagnosing global envelope properties (rotational support, turbulent mixing, recycling) directly from the flow fields. The central claim that these properties are mostly sensitive to the Bondi-to-core radius ratio is an empirical outcome of the simulation suite, not an analytic derivation, fitted parameter, or self-referential definition. No equations reduce to inputs by construction, no self-citations are load-bearing for the result, and the work contains no ansatz or uniqueness theorem imported from prior author work. The derivation chain is self-contained as a numerical exploration under stated assumptions.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claims rest on the numerical solution of the Euler equations under an isothermal closure and zero viscosity; these are standard domain assumptions rather than new postulates.

axioms (1)

domain assumption The gas obeys an isothermal equation of state and experiences no viscosity.
Explicitly stated in the abstract as the simulation framework.

pith-pipeline@v0.9.0 · 5770 in / 1276 out tokens · 20934 ms · 2026-05-25T02:11:12.499831+00:00 · methodology

Envelopes of embedded super-Earths II. Three-dimensional isothermal simulations

Core claim

What carries the argument

If this is right

Where Pith is reading between the lines

Load-bearing premise

What would settle it

discussion (0)