Envelopes of embedded super-Earths I. Two-dimensional simulations
Pith reviewed 2026-05-25 09:44 UTC · model grok-4.3
The pith
Super-Earth envelopes acquire substantial rotational support once core mass exceeds the thermal mass, due to vortensity conservation in two-dimensional flow.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
In two-dimensional simulations of isothermal and adiabatic envelopes around protoplanetary cores, conservation of vortensity causes the gas to acquire substantial rotational support once the core mass exceeds the thermal mass. This rotational support is controlled by the finite size of the core across the entire envelope. Supersonic envelopes of higher-mass cores develop steady non-axisymmetric shocks that trigger mass accretion and turbulent mixing, while gas self-gravity has only weak effects on the envelope itself but alters the properties of the density waves in the disk.
What carries the argument
Conservation of vortensity in two-dimensional hydrodynamical flow, which enforces rotational support throughout the envelope once core mass surpasses the thermal mass.
If this is right
- One-dimensional models lose applicability for describing envelope structure around cores more massive than the thermal mass.
- The core radius directly determines the amount of rotational support present at all radii in the envelope.
- Supersonic envelopes develop non-axisymmetric shocks that enable ongoing mass accretion and internal turbulent mixing.
- Gas self-gravity produces only minor changes to envelope density and rotation but significantly affects the outer disk density waves.
Where Pith is reading between the lines
- Rotational support may slow the rate of envelope contraction and delay the runaway gas accretion phase that turns cores into giants.
- The two-dimensional shocks and mixing could leave observable signatures in the chemical composition or thermal structure of super-Earth atmospheres.
- Extending the same vortensity argument to three dimensions might reveal how vertical motions alter the rotational support and shock locations.
Load-bearing premise
Two-dimensional isothermal and adiabatic approximations with full resolution down to the core surface capture the essential envelope dynamics without needing three-dimensional effects or more realistic thermodynamics.
What would settle it
A three-dimensional simulation of an identical core and disk setup that shows negligible rotational support in the envelope above the thermal mass would falsify the central claim.
read the original abstract
Measurements of exoplanetary masses and radii have revealed a population of massive super-Earths --- planets sufficiently large that, according to one dimensional models, they should have turned into gas giants. To better understand the origin of these objects, we carry out hydrodynamical simulations of planetary cores embedded in a nascent protoplanetary disk. In this first paper of a series, to gain intuition as well as to develop useful diagnostics, we focus on two-dimensional simulations of the flow around protoplanetary cores. We use the pluto code to study isothermal and adiabatic envelopes around cores of sub- to super-thermal masses, fully resolving the envelope properties down to the core surface. Owing to the conservation of vortensity, envelopes acquire a substantial degree of rotational support when the core mass increases beyond the thermal mass, suggesting a limited applicability of one-dimensional models for describing the envelope structure. The finite size of the core (relatively large for super-Earths) also controls the amount of rotational support in the entire envelope. Steady non-axisymmetric shocks develop in the supersonic envelopes of high-mass cores, triggering mass accretion and turbulent mixing in their interiors. We also examine the influence of the gas self-gravity on the envelope structure. Although it only weakly alters the properties of the envelopes, the gas gravity has significant effect on the properties of the density waves triggered by the core in the protoplanetary disk.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript presents two-dimensional hydrodynamical simulations with the PLUTO code of isothermal and adiabatic envelopes around embedded planetary cores spanning sub- to super-thermal masses, fully resolved to the core surface. The central claim is that vortensity conservation produces substantial rotational support once core mass exceeds the thermal mass, limiting the applicability of one-dimensional models; additional findings include steady non-axisymmetric shocks that trigger accretion and mixing in high-mass cases, plus a weak effect of gas self-gravity on envelope structure but stronger effect on disk density waves.
Significance. If the vortensity-based mechanism holds under the stated 2D setup, the work supplies a concrete physical reason why 1D envelope models may fail for the observed super-Earth population and supplies useful diagnostics for future multidimensional studies. The direct integration of the hydro equations with explicit appeal to conservation properties, together with full resolution down to the core surface, constitutes a methodological strength.
major comments (3)
- [Abstract] Abstract: the claim that vortensity conservation produces rotational support is load-bearing, yet the same paragraph states that steady non-axisymmetric shocks develop in the supersonic envelopes of high-mass cores; shocks are known sources of vorticity, so the manuscript must quantify the fractional violation of vortensity conservation across the shocked regions (e.g., via a diagnostic plot or table) to substantiate the rotational-support conclusion.
- [Simulation setup] Simulation setup section: the abstract and introduction assert that the finite size of the core controls the amount of rotational support throughout the envelope, but no quantitative comparison (different core radii at fixed mass, or explicit scaling) is referenced; without this, the statement remains qualitative and cannot be assessed as load-bearing for the limited-applicability claim.
- [Results] Results section on envelope structure: the manuscript reports that envelopes acquire rotational support beyond the thermal mass, but provides no convergence tests or resolution study at the core surface; given that the central claim rests on numerical results whose details (resolution, boundary conditions, vortensity implementation) are not visible in the abstract, an explicit demonstration that the reported rotational support is insensitive to these choices is required.
minor comments (2)
- [Abstract] The abstract would be clearer if it stated the precise range of core masses (in Earth masses or thermal-mass units) and the disk parameters adopted.
- Figure captions should explicitly note whether the displayed quantities are time-averaged or instantaneous, especially for the shocked high-mass cases.
Simulated Author's Rebuttal
We thank the referee for the constructive and detailed report. We address each major comment below, agreeing where revisions are needed to strengthen the claims and providing the strongest honest defense of the existing analysis where appropriate.
read point-by-point responses
-
Referee: [Abstract] Abstract: the claim that vortensity conservation produces rotational support is load-bearing, yet the same paragraph states that steady non-axisymmetric shocks develop in the supersonic envelopes of high-mass cores; shocks are known sources of vorticity, so the manuscript must quantify the fractional violation of vortensity conservation across the shocked regions (e.g., via a diagnostic plot or table) to substantiate the rotational-support conclusion.
Authors: We agree that shocks are potential sources of vorticity and that explicit quantification is required to support the vortensity-conservation argument. The revised manuscript will include a new diagnostic figure (or table) computing the fractional vortensity change (Δζ/ζ) across the identified shock fronts for the high-mass cases. This diagnostic will demonstrate that violations remain localized and small enough (typically a few percent) that the global envelope retains the reported degree of rotational support from conservation along streamlines between shocks. We will also note that the shocks are steady and do not lead to cumulative vortensity generation that would invalidate the mechanism. revision: yes
-
Referee: [Simulation setup] Simulation setup section: the abstract and introduction assert that the finite size of the core controls the amount of rotational support throughout the envelope, but no quantitative comparison (different core radii at fixed mass, or explicit scaling) is referenced; without this, the statement remains qualitative and cannot be assessed as load-bearing for the limited-applicability claim.
Authors: The referee correctly notes the absence of an explicit quantitative comparison. While the manuscript already varies core radius implicitly through the thermal-mass scaling and reports the resulting envelope structure, we will add a dedicated panel or subsection in the revised version that directly compares rotational-support profiles at fixed core mass but different core radii (or an explicit scaling relation). This will make the controlling role of core size quantitative and directly testable. revision: yes
-
Referee: [Results] Results section on envelope structure: the manuscript reports that envelopes acquire rotational support beyond the thermal mass, but provides no convergence tests or resolution study at the core surface; given that the central claim rests on numerical results whose details (resolution, boundary conditions, vortensity implementation) are not visible in the abstract, an explicit demonstration that the reported rotational support is insensitive to these choices is required.
Authors: We acknowledge that a focused resolution study at the core surface is not presented. The revised manuscript will add an appendix or subsection containing resolution-convergence tests (varying grid spacing near the core while keeping the same boundary conditions and vortensity advection scheme) that demonstrate the measured rotational support and vortensity profiles remain insensitive above a minimum resolution threshold. This will confirm the robustness of the central numerical result. revision: yes
Circularity Check
No circularity: results follow from direct 2D hydro integration and standard vortensity conservation
full rationale
The paper reports outcomes of numerical simulations solving the 2D hydro equations with the PLUTO code, fully resolving down to the core surface under isothermal and adiabatic assumptions. The key statement attributes rotational support to vortensity conservation, which is an intrinsic property of the inviscid Euler equations being integrated rather than a fitted parameter, self-defined quantity, or result imported via self-citation. No derivation chain reduces to its own inputs by construction; the findings are presented as emergent from the simulations. The 2D-to-3D applicability concern is a limitation of scope, not evidence of circularity in the reported chain.
Axiom & Free-Parameter Ledger
axioms (2)
- standard math Vortensity is conserved in the flow around the core
- domain assumption 2D isothermal and adiabatic approximations suffice for envelope structure
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.