{"id":"a5a7d061-3da6-410b-ae1e-4b39eea90528","arxiv_id":"1908.02742","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Planets with well-mixed, vapor-polluted envelopes enter runaway gas accretion when core plus vapor mass exceeds a predicted critical metal mass, and later cooling is slowed by dilution of the heavy vapor.","lead":"This paper builds an analytical model showing that growing planets with hot, vapor-polluted atmospheres can suddenly start swallowing gas once the total mass of solid material, in the core plus vapor, passes a threshold called the critical metal mass. The result gives a new formation path for super-Earths and mini-Neptunes and explains why gas giant formation may be more likely farther out.","discovery_kind":"first_principles","skeptic_critique":{"model":"deepseek-v4-flash","headline":"At the crossover mass defining Mz,crit, the paper's own energy estimate (Eq. 42b) gives Eenv > 0 for heavily polluted envelopes, so the assumed quasi-static, well-mixed structure is unbound exactly where the new runaway criterion applies.","rationale":"The paper's strongest claim is the critical metal mass criterion and its scaling. The reader's weakest assumption (perfect mixing and dynamical stability) is indeed the right place to look. My stress-test sharpens it: the model's own virial/energy calculation implies that at the crossover mass where Mxy = Mz, heavily polluted envelopes have positive total energy and cannot be in the assumed equilibrium. This is more specific than a generic concern about unknown mixing efficiency: it is an internal tension between §4.2/§5.2.1 and §7. However, the paper is transparent about the ideal-gas and uniform-mixing limitations, explicitly flags the stability issue, and presents the analysis as trend-identifying rather than precise. A conditional acceptance requiring a numerical check of envelope stability, mixing, and the crossover approximation remains the appropriate disposition; the reader already reached that verdict. I agree with the reader's identification of the load-bearing assumption, but I add the internal energy contradiction and a concrete test that the reader did not specify, hence partial agreement.","tokens_in":29846,"tokens_out":15354,"duration_ms":168706,"concrete_test":"Evaluate the total envelope energy (Eq. 42b, without the rvap≫rc limit) at the critical metal mass for the fiducial heavily polluted case of Fig. 7 (0.1 AU, Mc = 1 M⊕, fz obtained from Eq. 29). If Eenv > 0, no bound hydrostatic solution exists at Mz,crit under the model's own assumptions and the runaway criterion loses its foundation; if Eenv < 0 once the full radial integral is performed, the concern is resolved. A hydrostatic solve with the paper's own bvp_solve setup at the same parameters would corroborate the sign of Eenv.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is Eq. 30: runaway begins when the total metal mass exceeds Mz,crit ≈ 4 M⊕ (...). The derivation assumes the envelope remains a quasi-static, uniformly mixed hydrostatic structure with γg < 4/3 up to that point. However, §5.2.1 shows (Eq. 42b) that the total energy of such a polluted envelope turns positive when Menv/Mc ≳ 0.43. At the crossover that defines the critical metal mass, Menv = Mz(1+fz) and Mc = Mz(1−fz), so Menv/Mc = (1+fz)/(1−fz) ≥ 1 for every fz ≥ 0. For the heavily polluted cases (small cores, large fz) in which the new criterion differs most from the classical critical core mass and in which the inner high-Z region dominates the envelope mass, this ratio is well above the threshold. The paper explicitly states in §5.2.1 that such a positive-energy configuration is 'clearly non-physical' and in §7 that dynamical instability 'would fundamentally change their evolution and likely invalidate our results,' but it does not connect this to the validity of the runaway criterion at Mz,crit. Thus the load-bearing condition—existence of a bound, well-mixed envelope up to the critical metal mass—is not established and is internally contradicted in the regime where the criterion is most distinctive.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents an analytical three-zone model for the interior structure of proto-planets whose envelopes are polluted by vaporized solids. The model assumes an ideal-gas equation of state, a uniformly mixed inner high-Z vapor region with constant adiabatic index γg < 4/3, a step-function compositional transition, and neglect of Ledoux stabilization. Using this model, the authors describe four evolutionary phases: direct core growth, polluted envelope growth, embedded cooling, and post-disk indirect core growth. The central result is a critical metal mass, Eq. (30), that marks the onset of runaway gas accretion when the total metal mass (core plus vapor) exceeds ~4 M⊕ for nominal parameters; this is claimed to supersede the traditional critical core mass. A second major result is that compositional dilution during embedded cooling slows gas accretion, and a third is that post-disk vapor sedimentation produces only slow indirect core growth without significant mass-loss. The derivation is internally consistent in a formal sense and is cross-checked against a numerical integration of the structure equations (Fig. 3) and against phase-I core masses from BVO18 (Fig. 5).","tokens_in":30182,"tokens_out":11763,"duration_ms":137638,"significance":"If the central claims hold, this would be the first analytical framework for the evolution of polluted planetary envelopes, with testable scalings for the critical metal mass as a function of orbital distance, opacity, accretion rate, volatility, and core mass. The model usefully identifies dilution as a self-limiting mechanism in embedded cooling and makes a falsifiable prediction that volatile-poor inner-disk planets require more solids to reach runaway than volatile-rich outer-disk planets. The paper is transparent about many of its simplifications, and the comparison with a numerical structure integrator in Fig. 3 is a genuine check. However, the main result is conditional on the perfect-mixing assumption, which the paper itself notes is in direct opposition to the no-mixing numerical models of Bodenheimer et al. (2018). More seriously, the paper's own energy calculation in §5.2.1 implies that the quasi-static well-mixed envelope is unbound at the very crossover that defines the critical metal mass, a contradiction that is not reconciled in the text. The manuscript deserves publication only after this load-bearing issue is addressed.","major_comments":[{"comment":"The runaway criterion Mz,crit is derived for the crossover Mz = Mxy, where the envelope mass is Menv = (1+fz)Mz and the core mass is Mc = (1−fz)Mz. For any fz ≥ 0, Menv/Mc = (1+fz)/(1−fz) ≥ 1. Equation (42b), with the default γg = 1.25, gives Eenv > 0 for Menv/Mc ≳ 0.43, and the authors state in §5.2.1 that this positive-energy configuration is “clearly non-physical” and in §7 that dynamical instability “would fundamentally change their evolution and likely invalidate our results.” The paper does not connect this statement to the validity of the runaway criterion. In the heavily polluted regime (small cores, large fz), where the inner high-Z region dominates the envelope mass, the quasi-static, well-mixed envelope assumed in the derivation of Eq. (30) is unbound at the crossover. The central claim therefore needs either a demonstration that the envelope reaches Mz,crit before becoming dynamically unstable, or a re-derivation of the criterion from the bound-energy limit. As written, the criterion is internally contradicted in exactly the regime where it differs most from the classical critical core mass.","section":"§4.2 (Eqs. 25–30) and §5.2.1 (Eq. 42b)"},{"comment":"The conclusion that pollution accelerates gas accretion and lowers the mass at runaway rests entirely on the assumption that the convective high-Z region is perfectly mixed with uniform composition. The paper itself notes that Bodenheimer et al. (2018), who assume no compositional mixing, find the opposite trend: in their model a saturated inner region slows gas accretion. Since the mixing efficiency in planetary envelopes is poorly constrained, the sign of the effect is not established by this model. I am not asking for a full parameter study, but the abstract's statement that the critical metal mass “supersedes the traditional critical core mass” should be qualified to make clear that this is a consequence of the perfect-mixing assumption. A concrete step would be to state the regime of validity explicitly and to identify what would falsify the prediction — for example, a numerical simulation with finite compositional diffusion that brackets the mixed and unmixed limits.","section":"§2.3.2 and §7"}],"minor_comments":[{"comment":"The normalization of κrcb is written as “0.01 g cm−2”; opacity κ has units of cm2 g−1 in Eq. (6), so the label should read “0.01 cm2 g−1”.","section":"Eq. (30)"},{"comment":"The superscript “1” before “Impacts and rainout” appears to be a misplaced footnote marker; the sentence reads awkwardly and should be reformatted.","section":"Fig. 5 caption"},{"comment":"The text says the sedimentation condition is Tcg = Tvap, but Eq. (46) is written in terms of Trcb and rrcb/rc; please clarify how the condition T(rc) = Tvap translates to the stated expression.","section":"§6.1, Eq. (46)"},{"comment":"The phrase “a negative scaling with MenvMc” is ambiguous; it should be “a negative scaling with Menv Mc” (the product of the envelope and core masses).","section":"§5.2.1"}],"recommendation":"major_revision","confidential_remarks":"The stress-test concern is genuine: Eq. (42b) directly undermines Eq. (30) unless the runaway is explicitly redefined as the dynamical instability itself. I believe this is fixable in revision, which is why I recommend major revision rather than rejection. The authors are clearly aware of the issue but have not connected it to their principal claim; a revision that either derives the critical mass from the bound-energy limit or carefully restricts the claim to the regime where Eenv < 0 would substantially strengthen the paper."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Read this if you care about super-Earth formation. It provides the first analytical model where the runaway criterion is a critical metal mass (core + vapor), replacing the classical critical core mass, and it identifies dilution of the mean molecular weight as a self-limiting mechanism during embedded cooling. The four-phase framework is a useful organizing picture, and the structure is transparent: they check the analytical structure against a numerical integration (Fig. 3) and against their earlier BVO18 core masses (Fig. 5, with corrected opacity). The authors also clearly flag their main assumptions.\n\nThe soft spot is not hidden, but it is more serious than they present it. In Sect. 5.2.1 they show that a uniformly mixed polluted envelope with gamma_g < 4/3 has positive total energy when Menv/Mc >~ 0.43, and they call such a configuration 'clearly non-physical.' Yet at the crossover mass that defines the critical metal mass, Menv = Mz(1+f_z) and Mc = Mz(1-f_z), so Menv/Mc >= 1 for every non-zero pollution fraction. For exactly the cases where their new criterion departs from the classical one—small cores, large f_z—the assumed quasi-static structure is unbound well before the criterion is reached. They mention in Sect. 7 that dynamical instability would invalidate their results, but they don't make the link to the derivation of the central criterion. That is a load-bearing internal tension, not a peripheral caveat.\n\nOther concerns run smaller. The perfect-mixing assumption sits at the opposite extreme from Bodenheimer et al. (2018), and mixing efficiency in hot envelopes really is unknown; both trends cannot hold. The ideal-gas EOS is a rough approximation, and Tvap is fitted to BVO18 in a mildly circular way since it sets the normalization of Eq. 30. These are flagged by the authors themselves, and they are acceptable if the stability tension is resolved.\n\nWho is this for? Planet formation theorists interested in the origin of the super-Earth/sub-Neptune population, and people working on envelope pollution in core accretion. It deserves a serious referee. I would send it out, but with a demand for either a stability analysis of the polluted envelope approaching runaway (e.g., along the lines of Wuchterl 1990) or a numerical simulation that tracks whether the supposedly well-mixed structure remains bound. As it stands, the central claim is a good idea whose supporting structure is not yet self-consistent.","headline":"A clean analytical framework for polluted-envelope runaway that is internally contradicted by its own energy calculation at the regime it targets.","tokens_in":30718,"tokens_out":3436,"would_cite":true,"duration_ms":36961,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Critical metal mass, not core mass, triggers runaway gas accretion","keywords":["pebble accretion","envelope pollution","critical metal mass","runaway gas accretion","super-Earth formation","sub-Neptune envelopes","rainout","planet formation"],"falsifier":"Compute the deep interior with a realistic equation of state for hydrogen-helium plus silicate vapor at the temperatures and pressures of a forming super-Earth and check whether the effective adiabatic index stays below 4/3 while the composition remains uniform; a mixed envelope that is dynamically unstable, or that develops a compositional gradient instead of mixing, would break the critical metal mass and dilution-limited cooling.","tokens_in":29567,"feed_emoji":"🪐","tokens_out":10887,"duration_ms":103950,"temperature":0.7,"pith_summary":"This paper argues that a planet growing by pebble accretion turns into a gas giant not when its core reaches a traditional critical core mass, but when the total mass of heavy elements it has collected, core plus vapor held in the envelope, exceeds a value the authors call the critical metal mass. The authors construct an analytical interior model with a hot, uniformly mixed layer of vaporized rock around the core and derive an approximate threshold, about $4.0\\,M_\\oplus$ for silicate pollution at nominal parameters, that scales weakly with opacity, disk distance, and solids accretion rate and more strongly with the pollutant's evaporation temperature and with the core mass. If this criterion is right, planets that form close to their star, accreting less volatile material, must swallow more solids before runaway accretion begins, while planets outside the ice line, accreting volatile ices, can become gas giants at smaller solid masses. The same model splits a polluted planet's evolution into four phases, ending with slow rainout of the vapor onto the core after the disk disperses, with no strong mass loss.","feed_headline":"Critical metal mass, not core mass, triggers runaway gas accretion","feed_subtitle":"Core plus vapor mass, not core mass alone, sets when a growing planet runs away into a gas giant.","key_machinery":"The load-bearing object is the three-zone analytical envelope: an outer isothermal radiative region, an intermediate metal-free convective region with adiabatic index $\\gamma_{xy} > 4/3$, and an inner high-Z vapor layer that is assumed perfectly mixed, with mean molecular weight $\\mu_g \\simeq \\mu_{xy}/(1-Z)$ and adiabatic index $\\gamma_g < 4/3$. Because $\\gamma_g$ lies below $4/3$, most of the vapor mass sits near the core, so the envelope mass is controlled by this inner layer, and runaway is approximated by the crossover $M_z = M_{xy}$; combining these pieces produces the critical metal mass formula (Eqs. 29-30). The same machinery generates the dilution-limited gas accretion scaling for embedded cooling (Eq. 43) and the rainout timescale estimate after disk dispersal (Eq. 51).","core_discovery":"The paper's central claim is that the classical critical core mass should be replaced, for planets with polluted envelopes, by a critical metal mass $M_{z,\\mathrm{crit}}$: runaway gas accretion begins when $M_z = M_c + M_{\\mathrm{vapor}}$ exceeds this threshold. For the nominal silicate case the model gives $M_{z,\\mathrm{crit}} \\approx 4.0\\,M_\\oplus\\,(\\kappa_{\\mathrm{rcb}}/0.01\\,\\mathrm{g\\,cm^{-2}})^{1/6}(d/\\mathrm{AU})^{7/108}(T_{\\mathrm{vap}}/2500\\,\\mathrm{K})^{8/27}(\\dot{M}_z/10^{-5}\\,M_\\oplus\\,\\mathrm{yr^{-1}})^{1/6}(M_c/M_\\oplus)^{1/2}$. Pollution lowers the solid mass at which runaway sets in relative to a metal-free envelope, because the heavy vapor raises the mean molecular weight and lowers the adiabatic index, pulling in nebular gas; the threshold nevertheless rises for larger cores, which dilute the pollution, and for less volatile pollutants, which produce a smaller vapor region. The same framework yields a four-phase evolutionary sequence, with direct core growth, vapor-dominated envelope growth, dilution-limited embedded cooling, and post-disk rainout as the final stage.","pith_inferences":["If the critical metal mass holds, the observed close-in pile-up of sub-Neptunes may partly reflect a formation barrier, the extra solids needed to trigger runaway inside the ice line, rather than only later photo-evaporation sculpting the distribution.","A natural next test is a realistic equation-of-state calculation for mixed hydrogen-helium and silicate vapor in the deep interior, to check whether the effective adiabatic index really stays below $4/3$ and whether such layers are dynamically stable; the paper itself flags this as its largest uncertainty.","The model suggests a correlation between a sub-Neptune's envelope heavy-element abundance and its formation location, with more volatile-rich interiors tracing formation farther out, a prediction that could be checked against atmospheric abundances and orbital architectures.","Because the competing non-mixed picture gives the opposite gas-accretion trend, measuring how efficiently hot planetary interiors mix, through rainout outcomes or through the observable dilution of envelope metals, would decide which scenario shapes super-Earth evolution."],"forward_implications":["Runaway gas accretion is set by the total metal mass, so planets that vaporize their accreted solids can become gas giants at lower solid mass than the classical critical core mass predicts.","Planets inside the ice line, accreting silicates with high evaporation temperature, need more solids to reach runaway, biasing the inner disk toward super-Earths and sub-Neptunes and the outer disk toward gas giants.","If solids accretion stops while the disk is still present, the inflow of hydrogen-helium dilutes the heavy-element vapor and slows further contraction, giving sub-Neptune-mass planets a longer runway before runaway.","After disk dissipation, the contracting envelope eventually condenses its vapor, adding mass to the core over several Gyr, and the energy released comes too late to power significant outflow.","Planets that keep their primordial envelopes can retain part of their heavy elements outside the core after billions of years, so a diffuse heavy-element region in a giant planet need not imply core erosion."],"supporting_citations":[{"why":"It supplies the impact-vaporization and rainout behavior that phase I reproduces and against which the model is calibrated.","marker":"Brouwers et al. (2018)"},{"why":"It is the non-mixed, saturated-interior counter-model whose gas-accretion trend is opposite to the paper's, motivating the uniform-mixing assumption.","marker":"Bodenheimer et al. (2018)"},{"why":"These numerical studies show that envelope pollution enhances nebular gas accretion, the effect the critical metal mass quantifies.","marker":"Venturini et al. (2015, 2016)"},{"why":"It gives the early identification that high mean molecular weight in a polluted envelope lowers the mass threshold for runaway, the mechanism generalized here.","marker":"Stevenson (1982)"},{"why":"It provides the analytical critical-mass and envelope-mass framework and the self-gravity caveat that the model adapts to polluted interiors.","marker":"Piso & Youdin (2014)"},{"why":"It is the source of the adiabatic convective-structure equations used for the intermediate metal-free envelope region.","marker":"Ginzburg et al. (2016)"},{"why":"It is the benchmark for embedded cooling and nebular gas accretion against which the dilution-limited cooling is compared.","marker":"Lee & Chiang (2015)"}],"fun_headline_variants":["Polluted envelopes rewrite the rule for runaway gas accretion","Runaway gas accretion hinges on metal mass, not core mass","Metal mass, not core mass, sets the gas-giant trigger","Vapor-laden envelopes shift the trigger for planet runaway"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole argument assumes that vaporized rock mixes evenly through the hot inner atmosphere and that this atmosphere remains dynamically stable during formation; if mixing is weak or the atmosphere can collapse, the critical metal mass and dilution mechanisms fail.","fun_headline_variants_meta":{"raw":{"variants":["Polluted envelopes rewrite the rule for runaway gas accretion","Runaway gas accretion hinges on metal mass, not core mass","Metal mass, not core mass, sets the gas-giant trigger","Vapor-laden envelopes shift the trigger for planet runaway"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000999,"raw_usage":{"total_tokens":4325,"prompt_tokens":1139,"completion_tokens":3186,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":755,"completion_tokens_details":{"reasoning_tokens":3117}},"tokens_in":755,"tokens_out":3186,"duration_ms":23543,"temperature":1.0,"reasoning_tokens":3117,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T14:36:35.901457+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the deep interior with a realistic equation of state for hydrogen-helium plus silicate vapor at the temperatures and pressures of a forming super-Earth and check whether the effective adiabatic index stays below 4/3 while the composition remains uniform; a mixed envelope that is dynamically unstable, or that develops a compositional gradient instead of mixing, would break the critical metal mass and dilution-limited cooling.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It gives the early identification that high mean molecular weight in a polluted envelope lowers the mass threshold for runaway, the mechanism generalized here."}],"review_version":1}