{"id":"8409406a-6471-4477-aa68-d1928c121dcf","arxiv_id":"1909.00759","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":14,"one_line_summary":"The characteristic mass of super-Earths is set by the pebble isolation mass, giving roughly 1 Earth mass around a 0.08-solar-mass star and about 20 Earth masses around a solar-mass star.","lead":"This paper simulates how small planets grow by eating pebbles in disks around stars of different sizes, and finds their final masses follow the mass of their star. It suggests a simple rule: a growing planet stops growing when it opens a gap that blocks pebbles, and that stopping mass sets the typical super-Earth mass.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The observed Mp–Mstar trend is also produced by migration-limited growth in the paper's own Appendix A.2 and may be raised above Miso by giant impacts in multi-planet systems (Sect. 5.2.1), so the trend does not uniquely support pebble isolation as the mass-setting mechanism.","rationale":"The reader's weakest assumption was that the Bitsch et al. (2018) pebble isolation formula extrapolates unchanged to 0.08–1 Msun disks. My reading identifies a different but related load-bearing concern: even if that formula is perfectly valid, the paper's own appendix and caveats show that the observed Mp–Mstar trend is not a unique consequence of pebble isolation. Appendix A.2 demonstrates a migration-limited scenario without Miso that also yields a linear scaling, and Sect. 5.2.1 concedes that giant impacts in multi-planet systems can push final masses above Miso. Since the observational comparison is qualitative and the simulated systems are single-planet, the central claim that the characteristic mass 'may be set' by pebble isolation is plausible but underdetermined. This supports a CONDITIONAL verdict rather than ACCEPT, exactly as the reader concluded, so I do not propose changing the verdict. A secondary arithmetic issue reinforces the need for caution: with the paper's stated Mdot_g ∝ Mstar^1.8 scaling, Eqs. (10) and (18) give hg ∝ Mstar^(1/15) at the ice line, not Mstar^(1/9), so Eq. (39) should be Miso ∝ Mstar^(6/5), not Mstar^(4/3); this shows the quantitative prediction is sensitive to the adopted disk scaling, though the qualitative trend survives. The proposed test—disabling Miso in the fiducial model—would directly settle whether the reproduced trend actually requires pebble isolation.","tokens_in":40377,"tokens_out":12257,"duration_ms":127418,"concrete_test":"Run the fiducial two-component-disk Monte Carlo from Table 2 (both Scenario A and Scenario B) with the pebble-isolation truncation disabled, letting embryos continue pebble accretion until they reach the inner cavity, and compare the upper envelope of Mp versus Mstar against Eq. (39) and Fig. 7. If the no-isolation run still yields a similar linear or slightly super-linear Mp–Mstar envelope, then the observed trend cannot be uniquely attributed to pebble isolation; if the envelope flattens or lies far above the observed super-Earth masses, the pebble-isolation mechanism is more strongly supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim requires that pebble isolation is the mechanism that sets the characteristic planet mass and produces the Mp–Mstar correlation. Two pieces of evidence inside the paper weaken this causal attribution. First, Appendix A.2 shows that when pebble accretion is allowed to continue past Miso (in a pure stellar-irradiation disk), the resulting planet masses still scale linearly with stellar mass (Fig. A.2), and the text explicitly notes that 'even though there would be no Miso constraint, the masses of the growing planets would still be limited by migration.' Thus the linear Mp–Mstar scaling is not a unique diagnostic of pebble isolation; a migration-limited growth channel can produce the same qualitative trend. Second, Sect. 5.2.1 states that the simulated planets are single embryos, while most observed close-in super-Earths reside in multi-planet systems where post-disk giant impacts can raise final masses above Miso. The paper therefore concedes that Miso is only a lower limit for the final masses of planets in multiples. The observed trend reproduced by the model is thus consistent with pebble isolation, but it does not demonstrate that pebble isolation is the operative mass-setting process. The central claim is underdetermined by the presented single-embryo simulations and qualitative observational comparisons.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents a pebble-driven population synthesis model for planet formation around stars from 0.08 to 1 solar mass. The model includes viscous and irradiated disk structures, pebble accretion, gas accretion, type I/II migration, and Monte Carlo sampling of disk parameters. The authors report that the characteristic super-Earth mass is set by the pebble isolation mass and scales with stellar mass, from about 1 Earth mass around late M dwarfs to about 20 Earth masses around solar-mass stars. They also discuss water fractions, metallicity trends, and compare the simulated populations with observed exoplanet masses, orbital distances, metallicities, and water-content inferences. The paper's central claim is that the observed Mp–Mstar correlation is produced by pebble isolation, but the authors themselves show in an appendix that a migration-limited growth scenario can produce a similar scaling.","tokens_in":40733,"tokens_out":3543,"duration_ms":39334,"significance":"If the central claim is established, the paper would provide a physically motivated explanation for the characteristic masses of super-Earths and their stellar-mass dependence, connecting pebble isolation to a major observed demographic trend. The model is quite comprehensive for a population synthesis study: it includes a two-component disk model, a wide parameter study, and multiple robustness tests in Appendices A.1–A.7, covering alternatives such as pure irradiated disks, early embryo formation, different stellar luminosity relations, and low turbulence levels. The paper also makes concrete, falsifiable predictions for planets around very low-mass stars. However, the manuscript's main attribution of the mass trend to pebble isolation is not fully supported by its own tests, and the observational comparison is qualitative rather than statistical. The result is therefore significant conditional on a sharper causal test.","major_comments":[{"comment":"The central attribution of the Mp–Mstar scaling to pebble isolation is underdetermined by the paper's own tests. Section 5.2.1 concedes that the simulated planets are single embryos and that for observed multi-planet systems the pebble isolation mass is only a lower limit, since post-disk giant impacts can raise final masses. Appendix A.2 then shows that a model with no pebble isolation mass, in which growth is truncated by migration, also produces a linear Mp–Mstar scaling (Fig. A.2). The paper's key conclusion in Sect. 6 therefore requires a quantitative model comparison that can distinguish pebble isolation from migration-limited growth; the qualitative agreement shown in Figs. 1 and 7 does not provide that discrimination.","section":"§5.2.1 and Appendix A.2"},{"comment":"The predicted scaling is inherited from the adopted Bitsch et al. (2018) isolation mass formula rather than independently derived. Equation (26) already contains an explicit linear Mstar factor, and Eq. (39) evaluates that formula at the ice line using the assumed Mdot_g ∝ Mstar^1.8 relation, yielding Miso ∝ Mstar^(4/3) rather than a purely linear relation. The paper should state this inheritance explicitly and test the sensitivity of the central slope to the extrapolation of Eq. (26) below about 0.3 solar masses, because the validity of that hydrodynamic fit at very low stellar masses is a load-bearing assumption.","section":"Eqs. (26), (27), (39)"},{"comment":"The observational support for the central trend is currently qualitative and does not account for selection effects. Fig. 1 mixes radial-velocity minimum masses with true masses from transit/TTV measurements, and no detection completeness function is modeled. The statement that the observed Mp–Mstar trend is not due to observational bias is asserted rather than demonstrated. A quantitative comparison, even a forward-modeled occurrence rate with a simple detection probability, would be needed to support the claim that the simulated population 'agrees well' with observations.","section":"§5.2.1 and Fig. 1"}],"minor_comments":[{"comment":"The abstract and Sect. 1 describe the relation as linear, while Eq. (39) gives Miso ∝ Mstar^(4/3); this discrepancy should be clarified, for example by stating that the superlinear exponent is weak and approximately linear over the considered range.","section":"Abstract and Eq. (39)"},{"comment":"The unit for the pebble internal density in Eq. (25) is written as g cm^-1 but should be g cm^-3.","section":"Eq. (25)"},{"comment":"There are minor typographical errors such as 'mangetohydynamical' in Sect. 2.1 and 'sumarrize' in Sect. 2.1.2; these should be corrected.","section":"Sect. 2.1 and Sect. 5.2.2"},{"comment":"The description of water fraction ranges in the text is sometimes given as ranges like '& 10% to 1%' in Sect. 6, which appears to be an inverted or incomplete interval; the intended ordering should be checked.","section":"Sect. 4.2, Fig. 6"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Worth reading carefully. This is the first pebble-driven population synthesis that spans 0.08 to 1 solar mass, and it makes concrete predictions: super-Earth masses should scale roughly linearly with stellar mass, from about an Earth mass around late M dwarfs to tens of Earth masses around solar-type stars, with distinctive water-content signatures depending on where embryos start. That is a real contribution. The paper also does something I respect: it tests its own assumptions in Appendices A.1–A.7, including the pure-irradiation disk, early seed formation, different luminosity relations, and very low turbulence. The authors are not hiding the warts.\n\nThe soft spots are real but not fatal. The central mass scaling enters through the Bitsch et al. (2018) pebble isolation formula, Eq. (26), which already contains Miso ∝ Mstar. Equation (39) then repackages that dependence. So the headline prediction is not an independent derivation; it is an application of an existing result to a broader stellar mass range. That is fine for a population synthesis paper, but it means the paper is not a test of the isolation-mass mechanism itself. The stress-test note is right: Appendix A.2 shows that migration-limited growth in a pure-irradiation disk also produces a linear Mp–Mstar scaling, so the observed trend does not uniquely point to pebble isolation. The paper also concedes that planets in multi-planet systems can be pushed above Miso by giant impacts, making Miso a lower limit for final masses. These caveats are stated clearly, but they do reduce the force of the main claim.\n\nThe observational comparisons are qualitative. Figure 1 shows a nice visual correlation with a handful of low-mass-star systems, and the model reproduces the rough metallicity trends for super-Earths versus gas giants, but there is no statistical fit, no completeness correction, and no quantitative comparison between the simulated and observed mass distributions. The model also overproduces hot Jupiters relative to cold giants, which the authors acknowledge and attribute to missing multi-body dynamics. That is a minor blemish at this level of modeling, but worth noting.\n\nNo code or data artifacts are provided, so the Monte Carlo implementation cannot be checked independently. The equations are mostly internally consistent, and the parameter choices are reasonable and clearly stated.\n\nWho should read this? Anyone working on super-Earth formation, pebble accretion, or M-dwarf planet populations. It is a useful framework for making predictions that TESS and SPECULOOS can test. I would send it to a serious referee. My own recommendation would be conditional acceptance with two requests: a quantitative statistical comparison to the observed sample, and a section that honestly discusses whether the observed scaling can discriminate between pebble isolation and migration-limited growth.","headline":"First pebble-accretion population synthesis across 0.08–1 solar masses, with a clean super-Earth mass–stellar mass prediction that is largely inherited from the adopted isolation-mass formula; worth refereeing, but not decisive on its own.","tokens_in":41283,"tokens_out":1680,"would_cite":true,"duration_ms":22403,"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":"The paper argues that the pebble isolation mass, the point at which a growing planet halts pebble accretion, sets the characteristic super-Earth mass and predicts a near-linear scaling with host-star mass.","keywords":["planet formation","pebble accretion","pebble isolation mass","super-Earths","M dwarfs","population synthesis","exoplanet demographics","protoplanetary disk evolution"],"falsifier":"Run 3D hydrodynamic simulations of a growing planet in a disk around a 0.1 solar-mass star and measure the actual pebble isolation mass; if it differs substantially from the extrapolated value, the central scaling fails. Observationally, precise masses of close-in planets around late M dwarfs provide a test: the model predicts an upper envelope near 1 to 2 Earth masses around a 0.08 solar-mass star, so discovering a population of rocky planets at several Earth masses around such stars would falsify the isolation-mass truncation.","tokens_in":40170,"feed_emoji":"🪐","tokens_out":11552,"duration_ms":188892,"temperature":0.7,"pith_summary":"Planets forming in disks around stars from 0.08 to 1 solar mass stop growing by pebble accretion when they reach the pebble isolation mass, the mass at which a planet's gravity reverses the local disk pressure gradient and halts inward-drifting pebbles. The paper builds a pebble-driven population synthesis model and shows that this cutoff sets the characteristic super-Earth mass, giving roughly one Earth mass around a late M dwarf and about 20 Earth masses around a solar-mass star. The predicted scaling, approximately $M_{\\rm iso} \\propto M_\\star^{4/3}$, reproduces the observed trend that more massive stars host more massive super-Earths and that gas giants appear only above a threshold stellar mass. The same model also yields the observed weak metallicity dependence of super-Earth formation, the strong metallicity preference of gas giants, and water-content differences that depend on where embryos are born.","feed_headline":"Pebble isolation sets super-Earth mass, scaled by host star","feed_subtitle":"One disk process yields ~1 Earth-mass planets around late M dwarfs and ~20 Earth-mass planets around Sun-like stars.","key_machinery":"The key mechanism is the pebble isolation mass, defined as the planet mass at which the planet's gravitational perturbation opens a gap and reverses the disk's local pressure gradient, stopping the inward drift of pebbles and therefore ending pebble accretion. The paper adopts a 3D hydrodynamic fitting formula, Eq. (26), for this mass and, using the disk aspect ratio at the water ice line in its viscously heated inner disk, reduces it to the simple scaling $M_{\\rm iso} \\simeq 25\\,(M_\\star/M_\\odot)^{4/3}\\,M_\\oplus$ (Eq. 39). This scaling carries the argument: it converts the observed planet$-$stellar-mass correlation into a prediction of core-growth truncation, and it is insensitive to details such as embryo birth location or the turbulent $\\alpha_t$ value, which changes $M_{\\rm iso}$ by only about 25% across an order of magnitude in $\\alpha_t$.","core_discovery":"The central claim, stated in Sect. 6, is that the characteristic core-dominated planet mass may be set by the pebble isolation mass. During growth a planet eventually opens a shallow gap and reverses the local pressure gradient, so pebbles stop drifting inward and accretion terminates; the mass at which this happens is $M_{\\rm iso} = 25\\,(M_\\star/M_\\odot)\\,(h_g/0.05)^3\\,M_\\oplus \\simeq 25\\,(M_\\star/M_\\odot)^{4/3}\\,M_\\oplus$ for the adopted inner-disk structure. Using a Monte Carlo synthesis of embryos that grow by pebble accretion, with gas accretion and type I/II migration included, the paper finds that the upper mass envelope of super-Earths follows this scaling from about 1 $M_\\oplus$ around a 0.08 $M_\\odot$ star to about 20$-$25 $M_\\oplus$ around a solar-mass star. Excluding gas giants, the simulated population reproduces the observed planet-mass$-$stellar-mass trend, the rarity of giant planets around low-mass stars, and the stronger metallicity dependence of gas giants compared with super-Earths.","pith_inferences":["The authors note that giant impacts after disk dispersal are not modeled; if collisions add mass, final planet masses in multi-planet systems could exceed the isolation mass, so the cleanest test of the scaling is single-planet systems or the lower mass envelope of close-in super-Earths.","A direct observational discrimination could come from transit and radial-velocity surveys of ultracool dwarfs: if planets around 0.08$-$0.1 $M_\\odot$ stars are found at masses well above about 2 $M_\\oplus$, the isolation-mass truncation would need revision or an extra mass source.","Because the scaling is carried by the inner-disk aspect ratio, disks that are significantly hotter or colder than the adopted structure would shift the normalization and exponent of the planet$-$star mass relation, making the observed correlation a potential indirect probe of inner disk temperatures."],"forward_implications":["Super-Earth masses around low-mass stars should be capped near the pebble isolation mass, so systems around late M dwarfs are expected to host Earth- to few-Earth-mass planets and no gas giants.","Because $M_{\\rm iso}$ stays below roughly 10 $M_\\oplus$ for stars below about 0.3 $M_\\odot$, such systems should rarely or never form gas giants through this channel, matching the observed absence of massive planets around very low-mass stars.","Gas giant formation should be strongly favored around stars above roughly 0.3 $M_\\odot$ and in metal-rich disks, while super-Earth formation should be nearly metallicity-independent; both trends agree with current exoplanet populations.","The water content of a super-Earth depends on where its embryo formed: ice-line embryos in low-turbulence disks end up with about 10$-$15% water by mass, while embryos spread across the disk produce a bimodal mix of nearly dry and water-rich planets."],"supporting_citations":[{"why":"Supplies Eq. (26), the 3D hydrodynamic pebble isolation mass fitting formula that is the direct source of the stellar-mass scaling.","marker":"Bitsch et al. 2018"},{"why":"Introduced the pebble isolation mass concept, the stopping condition that sets the super-Earth mass scale.","marker":"Lambrechts et al. 2014"},{"why":"Established pebble accretion as the growth mechanism the model uses before isolation quenches it.","marker":"Lambrechts & Johansen 2012"},{"why":"Provides the type I migration torque prescription, including outward migration, that lets planets remain near the transition radius while growing.","marker":"Paardekooper et al. 2011"},{"why":"Reports the observed linear planet-mass$-$stellar-mass scaling among Kepler planets that the model aims to explain.","marker":"Wu 2019"},{"why":"Finds a universal truncation near planet-to-star mass ratio $3\\times10^{-5}$, the observational anchor for the mass correlation.","marker":"Pascucci et al. 2018"}],"fun_headline_variants":["Super-Earth mass scales with host star via pebble isolation","Pebble isolation sets super-Earth mass; star mass sets the scale","Super-Earth mass follows pebble isolation, increasing with star mass","Pebble isolation yields 1 Earth-mass planets around M dwarfs, 20 around Sun"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument stands on the assumption that the pebble isolation mass formula fitted to hydrodynamic simulations for solar-type disks remains valid for disks around stars from 0.08 to 1 solar mass, so that its dependence on stellar mass is as predicted.","fun_headline_variants_meta":{"raw":{"variants":["Super-Earth mass scales with host star via pebble isolation","Pebble isolation sets super-Earth mass; star mass sets the scale","Super-Earth mass follows pebble isolation, increasing with star mass","Pebble isolation yields 1 Earth-mass planets around M dwarfs, 20 around Sun"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001425,"raw_usage":{"total_tokens":5852,"prompt_tokens":1152,"completion_tokens":4700,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":768,"completion_tokens_details":{"reasoning_tokens":4618}},"tokens_in":768,"tokens_out":4700,"duration_ms":27763,"temperature":1.0,"reasoning_tokens":4618,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T05:37:31.187699+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run 3D hydrodynamic simulations of a growing planet in a disk around a 0.1 solar-mass star and measure the actual pebble isolation mass; if it differs substantially from the extrapolated value, the central scaling fails. Observationally, precise masses of close-in planets around late M dwarfs provide a test: the model predicts an upper envelope near 1 to 2 Earth masses around a 0.08 solar-mass star, so discovering a population of rocky planets at several Earth masses around such stars would falsify the isolation-mass truncation.","supporting_citations":[{"cited_title":"2011, MNRAS, 410, 293","cited_arxiv_id":null,"evidence_quote":"Provides the type I migration torque prescription, including outward migration, that lets planets remain near the transition radius while growing."},{"cited_title":"D., Gould, A., & Fernandes, R","cited_arxiv_id":null,"evidence_quote":"Finds a universal truncation near planet-to-star mass ratio $3\\times10^{-5}$, the observational anchor for the mass correlation."}],"review_version":1}