{"id":"8df92a84-a248-414c-813e-d8af7a077a5b","arxiv_id":"2412.16278","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"Planet formation rate scales as gas surface density to a power between 4/3 and 2, depending on formation mechanism.","lead":"This paper derives a power-law relation between the rate of planet formation and the gas surface density of protoplanetary disks, with the exponent depending on how the planets form. It is an attempt to give planet formation a scaling law like the Kennicutt-Schmidt law that links star formation to gas density.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Gas-giant n=2 hinges on unstated Σ_core ∝ Σ_g; the definition ϵ_core = Σ_core/Σ_p is circular, so Eq. (7) is not derived from first principles.","rationale":"The reader's CONDITIONAL verdict is the right one, and this stress-test identifies the same soft spot with more specificity. The exponents in Eq. (7) are the paper's entire scientific output; they are what a future observer or modeler would take away. The growth-rate pieces (pebble accretion, Bondi, gap-opening) have sound density scalings and are properly cited. The fragile step is the conversion from growth rate to planet formation rate in Eq. (6): the seed density Σ_p must be known as a function of Σ_g. For pebble-fed terrestrial growth, the chain Σ_peb = f_peb Σ_g and Σ_pro = ϵ_pro Σ_peb provides that relation, so the 3/2 and 4/3 exponents are reasonably supported. For gas giants, the manuscript's own efficiency definitions are circular, and no analogous chain fixes Σ_core. The proposed population-synthesis run is a direct, feasible check: it measures q in Σ_core ∝ Σ_g^q and thereby decides whether the n=2 line in Eq. (7) is correct or should be n=3 (or some other value). This concern does not require rejecting the paper; it requires the authors to either derive the seed-density scaling or relabel the result as a conditional relation. The observational comparison with only TW Hya and PDS 70 is too sparse to settle the slope, so it does not mitigate this concern. With the requested clarification, the paper remains a useful scaling-relation proposal. Hence the verdict stays CONDITIONAL (UNCHANGED).","tokens_in":9968,"tokens_out":13549,"duration_ms":115091,"concrete_test":"Run a published pebble-accretion population synthesis (e.g., Lambrechts & Johansen 2014, or a modern N-body/pebble code) on a grid of disks that differ only in gas surface density Σ_g (spanned over 1–2 dex), with fixed α, temperature profile, and disk lifetime. At t = 1 Myr, measure the total core mass per unit area Σ_core and fit q in Σ_core ∝ Σ_g^q. Since Γ_grow ∝ Σ_g for Bondi/gap-opening growth, Eq. (6) predicts n_gas_giant = 1+q. If q = 1, the n=2 values in Eq. (7) are confirmed; if q = 2, the exponents become 3 and the unstated Σ_core ∝ Σ_g assumption is falsified.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is Eq. (7), a set of exponents. For the bottom-up channels, Eq. (6) gives PFR = ϵ_PFR Σ_p Γ_grow, and the quoted exponents follow only if the seed surface density Σ_p scales linearly with Σ_g. The paper never derives this. For terrestrial planets Eq. (2) gives Σ_peb = f_peb Σ_g, so with constant ϵ_pro the linear scaling is reasonable. For gas giants the text defines Σ_p as the surface density of cores and then sets ϵ_core = Σ_core/Σ_p; that is tautological (ϵ_core = 1) and contradicts the adopted value 0.2. If the intended definition is ϵ_core = Σ_core/Σ_peb, the n=2 result still requires the extra assumption that core surface density is a constant fraction of pebble surface density. In drift-limited pebble accretion, core mass is the time-integrated pebble flux, while core number density is set by embryo spacing; neither is guaranteed to be proportional to Σ_g. If n_core ∝ Σ_g and M_core ∝ Σ_g, then Σ_core ∝ Σ_g^2 and the gas-giant exponent becomes 3 (since Γ_grow ∝ Σ_g). Thus the n=2 value in Eq. (7) is one plausible scaling, not a first-principles prediction, and the abstract's 'derived from first principles' is an overstatement.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript proposes an analog of the Kennicutt-Schmidt law for planet formation, claiming a first-principles derivation of a power-law relation PFR ∝ Σ_g^n for different formation channels. It derives n = 3/2 for terrestrial planets in the headwind pebble-accretion regime, n = 4/3 in the shear regime, n = 2 for gas giants, n = 2 for gap-opening gas giants, and n = 3/2 for gravitational instability (Eq. 7). The paper then compares the predicted scaling with crude PFR estimates for the TW Hya and PDS 70 disks, concluding that the observed trend is consistent with the theory. The central result is a set of characteristic exponents that could serve as observational discriminants for formation mechanisms.","tokens_in":10378,"tokens_out":7158,"duration_ms":61228,"significance":"If the derivation were fully rigorous, this would be a useful conceptual bridge between planet formation and star formation, offering falsifiable predictions that could be tested with future ALMA/JWST observations of planet-forming disks. The paper is clearly written and the bottom-up exponents follow from textbook accretion rates, so the core idea is attractive. However, the gas-giant exponent rests on an unstated assumption, the observational test is only a two-point comparison with orders-of-magnitude uncertainties, and the gravitational-instability timescale is chosen without clear physical justification. The strengths are the compact formulation of a potentially testable scaling law and the explicit enumeration of channel-specific exponents.","major_comments":[{"comment":"The gas-giant exponent n=2 is not derived from first principles because the scaling of the core surface density Σ_p with gas surface density is never established. The text defines ε_core = Σ_core/Σ_p, which is tautological (it forces ε_core=1 and contradicts the adopted value 0.2), and even if the intended definition is ε_core = Σ_core/Σ_peb, the result PFR ∝ Σ_g^2 requires the additional unstated assumption that Σ_core ∝ Σ_g. In drift-limited pebble accretion, the core mass is set by the integrated pebble flux and the core spacing by isolation conditions, so a linear proportionality is not guaranteed; the authors should either derive this scaling or explicitly state it as an assumption and discuss its range of validity.","section":"§II.B, Eq. (7)"},{"comment":"The observational comparison is too weak to support the claim that the data show the predicted trend. Only two disks are used, the area A = π(100 AU)^2 used to convert a single-planet accretion rate into PFR_obs is arbitrary, and the ranges for Σ_g and PFR overlap substantially between the two objects. The statement that 'this trend could be seen also from the observational data' is an overstatement; the comparison should be framed as an illustrative order-of-magnitude check, not as a validation of the derived exponents.","section":"§III, Fig. 1"},{"comment":"The use of the free-fall time t_ff as the growth timescale for gravitational instability is not justified for a rotating disk, where the natural dynamical timescale is the orbital period Ω^{-1}. Since Ω is independent of Σ_g at fixed radius and stellar mass, adopting Ω would give PFR_GI ∝ Σ_g rather than n=3/2. The authors should justify the t_ff choice physically or acknowledge that the GI exponent is only one of several possible scalings.","section":"§II.C, Eq. (5)"}],"minor_comments":[{"comment":"There is a typo: 'terrestiral' should be 'terrestrial'.","section":"Eq. (7)"},{"comment":"The two-case equation for the pebble accretion rate is missing a closing brace or delimiter; the Ṁ_PA expression is not fully typeset, which hampers readability.","section":"Eq. (1)"},{"comment":"The free-fall time is written as sqrt(2π/32Gρ_g); the standard expression is sqrt(3π/(32Gρ_g)) or equivalent. The numerical factor does not affect the scaling, but the formula should be correct.","section":"Eq. (5)"},{"comment":"The fiducial value ε_PFR = 1 is not an efficiency but means no efficiency reduction; the table entry should be clarified. Some rows are also misaligned in the table formatting.","section":"Table A"},{"comment":"The text says PFR 'measures the number of planets' but the formula ε_PFR Σ_p Γ_grow has dimensions of mass per area per time (a surface density growth rate). The wording should be aligned with the dimensional meaning.","section":"§III, Eq. (6)"}],"recommendation":"major_revision","confidential_remarks":"The paper presents a simple and appealing scaling idea, but the gas-giant exponent—one of the two headline results—rests on a hidden proportionality assumption that is not stated or derived. This is fixable by rewriting the efficiency definitions and explicitly discussing the Σ_core–Σ_g relation. The observational section is currently too assertive for a two-point, arbitrarily normalized comparison. If the authors address the assumption and soften the observational claims, the letter could be publishable. The novelty relative to existing planet-formation scaling literature may also warrant a brief comparative discussion."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper proposes a Kennicutt-Schmidt analog for planet formation, PFR ∝ Σ_g^n, with n depending on the formation channel: 3/2 for terrestrial planets in the headwind regime, 4/3 in the shear regime, 2 for gas giants and gap-opening gas giants, and 3/2 for gravitational instability. The core idea is genuinely useful: packaging known mass-growth scalings into a compact set of exponents that could be compared with disk surveys. The terrestrial and GI derivations are clean and follow directly from published growth rates (pebble accretion in the two regimes, free-fall time for GI). The author is also honest that this is a first proposal, not a settled law.\n\nThe soft spots are in the gas-giant channel and the observational validation. For gas giants, the n=2 result requires the seed surface density Σ_core to scale linearly with Σ_g, but that assumption is never stated or derived. The definition ε_core = Σ_core/Σ_p is tautological when Σ_p = Σ_core, and the adopted value 0.2 is inconsistent with that definition. This needs to be cleaned up: either define Σ_p differently or state explicitly that core surface density is assumed proportional to pebble surface density, with the efficiency absorbing the coefficient. If the core surface density scales as Σ_g^2 instead, the exponent changes to 3. So the 'first principles' claim in the abstract is too strong.\n\nThe observational test is a two-point comparison (TW Hya and PDS 70) with error bars spanning orders of magnitude. It is suggestive but not a validation. The author acknowledges this, so it's not a fatal flaw, but the paper should not present Fig. 1 as strong support.\n\nOverall, this is a reasonable, short contribution that could be useful after revision. It deserves peer review, but a referee should push for a clearer treatment of the seed-density scaling and a more cautious interpretation of the data. I'd suggest major revision rather than rejection.","headline":"A useful scaling-law package, but the gas-giant exponent rests on an unstated seed-density assumption and the observational test is only two points; worth peer review with major revision.","tokens_in":10827,"tokens_out":4934,"would_cite":false,"duration_ms":41966,"reading_group":"maybe","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 derives a star-formation-law analog for planet formation: $PFR \\propto \\Sigma_g^n$, with $n$ from $4/3$ to $2$ depending on the formation mechanism.","keywords":["planet formation rate","protoplanetary disks","gas surface density","pebble accretion","gas giant formation","gravitational instability","power-law scaling","star-formation scaling law"],"falsifier":"Take a sample of protoplanetary disks with measured gas surface densities and with directly detectable embedded planets whose accretion rates give independent planet formation rates, extend the two-point comparison made here to tens of disks, and fit $\\log PFR$ versus $\\log \\Sigma_g$ separately for terrestrial-type and giant-type systems; the central claim fails if the giant-planet slope is not close to $2$ or the terrestrial slope is not close to $3/2$ or $4/3$. A numerical simulation that varies disk gas density while holding the seed population fixed should reproduce the same slopes.","tokens_in":9791,"feed_emoji":"🪐","tokens_out":10656,"duration_ms":91882,"temperature":0.7,"pith_summary":"The paper tries to establish that the rate at which planets form in a protoplanetary disk is a power-law function of the gas surface density, $PFR \\propto \\Sigma_g^n$, in direct analogy to the star-formation scaling law measured in galaxies. The exponent $n$ is derived separately for each major formation channel, and the values are $3/2$ for terrestrial planets in the headwind pebble-accretion regime, $4/3$ in the shear regime, $2$ for gas giants in both runaway and gap-opening accretion, and $3/2$ for planets formed by gravitational instability. If the relation is correct, the gas content of a disk directly sets the pace of planet formation, giving observers a simple diagnostic for which mechanism produced a given planet population.","feed_headline":"Planet formation obeys a gas-density power law","feed_subtitle":"Each formation channel has its own slope, from 4/3 for terrestrial planets to 2 for gas giants, testable with disk surveys.","key_machinery":"The load-bearing estimator is $PFR \\approx \\epsilon_{PFR}\\,\\Sigma_p/t_{grow}$, where $t_{grow} = M/(dM/dt)$ is the mass-growth timescale of a typical seed. The power-law transfer from growth rate to gas density is carried by four physical inputs: the pebble-to-gas fraction $f_{peb}=\\Sigma_{peb}/\\Sigma_g=0.01$ that ties pebble accretion (capture of small gas-coupled particles by a growing protoplanet) to the gas reservoir, the Bondi accretion rate $\\dot{M}_{RA}\\propto\\rho_g$ for runaway envelope growth, a gap-opening accretion rate proportional to $\\Sigma_g$, and the free-fall timescale $t_{ff}\\propto\\rho_g^{-1/2}$ for gravitational instability. Efficiencies $\\epsilon_{pro}$ and $\\epsilon_{core}$ are set to $0.2$, following the pebble-flux treatment of core formation, and they set the overall normalization rather than the exponents.","core_discovery":"Defining the planet formation rate as $PFR \\approx \\epsilon_{PFR}\\,\\Sigma_p/t_{grow}$, with $\\Sigma_p$ the surface density of growing seeds and $t_{grow}$ the growth timescale from each accretion mechanism, the paper finds that all channels collapse to $PFR \\propto \\Sigma_g^n$ with $n = 3/2$ (terrestrial, headwind), $4/3$ (terrestrial, shear), $2$ (gas giants and gap-opening gas giants), and $3/2$ (gravitational instability). The key steps are writing the pebble surface density as a fixed fraction $f_{peb}\\Sigma_g$, treating runaway gas accretion as Bondi accretion proportional to gas density, using a gap-opening accretion rate linear in $\\Sigma_g$, and using the free-fall time for gravitational instability. The paper further claims this local law yields a global disk-integrated relation, applies above a threshold gas density, and predicts self-regulation as the gas is depleted.","pith_inferences":["The paper leaves implicit that, if the relation holds, surveys measuring disk gas mass could be converted into expected planet formation rates, letting modelers predict where planets should appear before they are directly imaged.","A clean statistical test would compare exoplanet occurrence rates against host-disk gas mass: a break from $n\\approx 4/3$-$3/2$ at low masses to $n\\approx 2$ at high masses would support mechanism-dependent slopes, while a single global slope would indicate that seed densities decouple from gas density.","Migrating planets complicate the comparison: planets form at the local gas density but are later observed far from their birth site, so matching the predicted slopes to observed planet positions requires accounting for orbital migration and disk evolution."],"forward_implications":["In the headwind pebble-accretion regime, terrestrial planet formation should speed up as $\\Sigma_g^{3/2}$, while in the shear regime the slope is shallower, $\\Sigma_g^{4/3}$, so resolved disk density maps can pick out which accretion mode operates.","Gas giant formation should scale as $\\Sigma_g^2$, making giant planets strongly concentrated in the most gas-rich disks and rare in gas-poor ones.","A measured slope near $3/2$ in a planet population would point to gravitational instability or headwind pebble accretion, whereas a slope near $2$ points to gas-accretion-dominated channels.","Integrated over a disk, the local law becomes a global relation between total planet formation rate and disk gas content, and it predicts that planet formation self-regulates downward as gas is depleted."],"supporting_citations":[{"why":"Supplies the empirical star-formation scaling law that the planet-formation relation is designed to mirror.","marker":"[11]"},{"why":"Provides the original star-formation rate power law whose exponent is compared with the derived planet-formation exponents.","marker":"[12]"},{"why":"Gives the gas-drag radial-drift physics that fixes the pebble flux and its proportionality to gas surface density.","marker":"[3]"},{"why":"Defines the headwind and shear pebble-accretion regimes and their mass-growth rates used for terrestrial planets.","marker":"[22]"},{"why":"Relates pebble surface density to gas surface density and supplies the core and protoplanet formation efficiencies adopted in the PFR estimator.","marker":"[23]"},{"why":"Provides the gap-opening gas accretion rate whose linear dependence on $\\Sigma_g$ gives the $n=2$ exponent for gap-opening gas giants.","marker":"[24]"},{"why":"Supports gravitational instability as a planet-formation route whose collapse timescale produces the $\\Sigma_g^{3/2}$ scaling.","marker":"[29]"},{"why":"Provides the TW Hya disk gas densities, planet masses, and accretion rates used as observational comparison data.","marker":"[33]"},{"why":"Provides the PDS 70 protoplanet masses and accretion rates used as the second observational comparison point.","marker":"[35]"}],"fun_headline_variants":["Planet formation rate follows gas density power law","New Kennicutt-Schmidt analog for planet formation","From 4/3 to 2: planet formation power-law slopes","Gas density dictates planet formation rate","Plantary growth obeys gas-density scaling law"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the surface density of growing seeds, $\\Sigma_p$, stays proportional to the gas surface density through a fixed efficiency factor; if seed production saturates or decouples from the gas reservoir, the exponents and normalizations of the planet formation rate change.","fun_headline_variants_meta":{"raw":{"variants":["Planet formation rate follows gas density power law","New Kennicutt-Schmidt analog for planet formation","From 4/3 to 2: planet formation power-law slopes","Gas density dictates planet formation rate","Plantary growth obeys gas-density scaling law"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000115,"raw_usage":{"total_tokens":1073,"prompt_tokens":948,"completion_tokens":125,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":564,"completion_tokens_details":{"reasoning_tokens":61}},"tokens_in":564,"tokens_out":125,"duration_ms":1983,"temperature":1.0,"reasoning_tokens":61,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T10:44:14.990100+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take a sample of protoplanetary disks with measured gas surface densities and with directly detectable embedded planets whose accretion rates give independent planet formation rates, extend the two-point comparison made here to tens of disks, and fit $\\log PFR$ versus $\\log \\Sigma_g$ separately for terrestrial-type and giant-type systems; the central claim fails if the giant-planet slope is not close to $2$ or the terrestrial slope is not close to $3/2$ or $4/3$. A numerical simulation that varies disk gas density while holding the seed population fixed should reproduce the same slopes.","supporting_citations":[{"cited_title":"Kennicutt, Robert C., The Star Formation Law in Galactic Disks, ApJ 344, 685 (1989)","cited_arxiv_id":null,"evidence_quote":"Supplies the empirical star-formation scaling law that the planet-formation relation is designed to mirror."},{"cited_title":"Schmidt, The Rate of Star Formation., ApJ 129, 243 (1959)","cited_arxiv_id":null,"evidence_quote":"Provides the original star-formation rate power law whose exponent is compared with the derived planet-formation exponents."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the gas-drag radial-drift physics that fixes the pebble flux and its proportionality to gas surface density."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the headwind and shear pebble-accretion regimes and their mass-growth rates used for terrestrial planets."},{"cited_title":"Final Masses of Giant Planets II: Jupiter Formation in a Gas-Depleted Disk","cited_arxiv_id":"1510.06848","evidence_quote":"Provides the gap-opening gas accretion rate whose linear dependence on $\\Sigma_g$ gives the $n=2$ exponent for gap-opening gas giants."},{"cited_title":"Goldreich and W","cited_arxiv_id":null,"evidence_quote":"Supports gravitational instability as a planet-formation route whose collapse timescale produces the $\\Sigma_g^{3/2}$ scaling."},{"cited_title":"Outflow Driven by a Protoplanet Embedded in the TW Hya Disk","cited_arxiv_id":"2407.14395","evidence_quote":"Provides the TW Hya disk gas densities, planet masses, and accretion rates used as observational comparison data."},{"cited_title":"Keck/NIRC2 $L$'-Band Imaging of Jovian-Mass Accreting Protoplanets around PDS 70","cited_arxiv_id":"2004.09597","evidence_quote":"Provides the PDS 70 protoplanet masses and accretion rates used as the second observational comparison point."}],"review_version":1}