{"id":"e81f3f89-5ff9-464e-83f2-99f24b3088c8","arxiv_id":"2507.14814","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":9,"one_line_summary":"Terrestrial planets form best in simulations with an inner planetesimal ring near 0.5 au plus an outer source near 1.5 to 2 au, with outward Type-I migration.","lead":"The paper runs 1000 simulations per model to argue that Mercury, Venus, and Earth grew from a narrow planetesimal ring near 0.5 au and later migrated outward, while Mars formed from a separate reservoir at 1.5 to 2 au. This two-source model is presented as the best statistical match to the real terrestrial planet system.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Isotopic claim is internally contradicted: model203 gives Mars ~90% outer-source material, yet §5 admits this 'would not have enough EC material to explain observations'; the abstract nevertheless claims the model explains Earth-Mars isotopic differences.","rationale":"The reader correctly identifies the parameterized MDW-like disk profile of Eq. 9 as a load-bearing assumption, and that concern is genuine: if real protoplanetary disks lacked the outward Type-I migration zone near 1 au, the inner-ring scenario for Venus and Earth would fail. The authors themselves defer self-consistent MDW disk models to future work. However, this is a physical assumption with some literature support and is not contradicted by the paper's own text. The isotopic mismatch is more damaging because it is an internal inconsistency: the paper's Section 5 explicitly admits that its preferred mass split produces a Mars composition inconsistent with the EC/OC constraints, while the abstract and conclusions present the isotopic difference as explained. The model may still be the best available dynamical match for masses and orbits, and the paper deserves credit for the large simulation campaign, resolution checks, and explicit caveats. But a headline claim that is contradicted by the authors' own quoted result cannot be accepted at face value. The reader's conditional verdict remains appropriate; no change in verdict is needed, but the paper should either revise the isotopic claim, adjust the reservoir mapping or mass split, or add a quantitative comparison to the Dauphas et al. (2024) posteriors before the central claim is treated as established.","tokens_in":48662,"tokens_out":7308,"duration_ms":93250,"concrete_test":"Use the already-recorded accretion histories behind Fig. 13 to compute, for every successful model203 simulation, the fraction of Earth's and Mars' mass accreted from the inner ring versus the outer source. Assign inner-ring material = EC and outer-source material = OC, the mapping the paper tentatively adopts, and compare the resulting bivariate distribution of (Earth EC fraction, Mars EC fraction) to the posterior constraints from Dauphas et al. (2024) including uncertainties. If fewer than ~5% of successful simulations fall within the observed 95% credible region, then the 'explains isotopic differences' claim fails as stated and should be removed or the reservoir assignment/mass split revised. This is a pure post-processing test requiring no new N-body runs.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The most load-bearing weakness is not the disk parameterization alone but the paper's admitted quantitative failure on one of its three headline explanations. In §5 the authors state that for the reference two-source model, 'Mars accretes ~90% of material from the outer source... This composition would not have enough EC material to explain observations.' Yet the abstract and §4.9 claim the model 'explains the isotopic differences of the Earth and Mars.' The isotopic constraints in §3.2.3 (Dauphas et al. 2024) are Mars ≈65% EC / ≈33% OC and Earth ≈95% EC. If, as the paper tentatively assumes, the inner ring is the EC reservoir and the outer source is the OC reservoir, model203 predicts Earth ≈70% EC / 30% OC and Mars ≈10% EC / 90% OC. The predicted Mars is OC-dominated, whereas the observed Mars is EC-dominated. The claimed isotopic explanation is therefore not demonstrated: the model can produce two isotopically distinct planets, but not the specific EC/OC mix inferred for Mars (or, strictly, the inferred Earth value either, unless an unsampled reservoir is invoked). This is not an outside-consensus disagreement; it is an internal inconsistency between a headline claim and the authors' own caveat. The paper's central claim should be treated as conditional on closing this gap.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents a large suite of N-body simulations of terrestrial planet formation starting from planetesimal disks with one or two radial reservoirs, a gas disk stage with Type-I migration and aerodynamic drag, a giant-planet instability, and a long late-stage integration. The authors explore roughly one hundred model configurations with one thousand simulations each and grade outcomes against objective criteria for Mercury, Venus/Earth, Mars, and the Venus-Earth separation. Their preferred two-source model (model203), with an inner ring at 0.6 au and an outer source at 1.7 au, yields success rates (V/E, Mars, Merc, V/E Sep) = (50%, 35%, 32%, 33%) and 6.1% of runs satisfying all four criteria. The paper claims that this model explains the masses and orbits of the terrestrial planets, the reduced/oxidized accretion history of Earth, and the isotopic differences between Earth and Mars, and that it provides a plausible accretion history for the Moon-forming impactor Theia.","tokens_in":49037,"tokens_out":2625,"duration_ms":34511,"significance":"If the central claims hold, this would be a substantial advance: the paper demonstrates that a two-reservoir initial condition plus outward Type-I migration in a wind-like disk can simultaneously match the radial mass concentration, the Mercury problem, the small-Mars problem, and the Earth's bimodal accretion chemistry. The statistical methodology is a clear strength: one thousand simulations per model, objective planet-matching criteria defined in Section 3.1, and systematic parameter variation with results tabulated in Tables 1 and 2. The fragmentation tests in Section 4.10 and the resolution checks with N=1000 and N=2000 are also valuable. However, the paper's third headline claim, that it explains the isotopic differences of Earth and Mars, is contradicted by the authors' own quantitative caveat in Section 5, and the Earth's 70/30 composition is largely built into the chosen mass split. The central mechanism also rests on a parameterized disk profile rather than a self-consistent disk model. These issues make the paper's central conclusion conditional rather than established.","major_comments":[{"comment":"There is an internal inconsistency between the abstract and the quantitative Mars composition. The abstract states that the model explains 'the isotopic differences of the Earth and Mars,' and Section 4.9 and Conclusion item 6 repeat this claim. However, Section 5 states that with the 2:1 mass split, 'Mars accretes ~90% of material from the outer source... This composition would not have enough EC material to explain observations.' Given the constraints quoted in Section 3.2.3 (Mars ~65% EC, Earth ~95% EC) and the tentative mapping of the inner ring to EC and outer source to OC, model203 predicts Mars to be ~90% OC and ~10% EC, which is the opposite of the inferred EC-dominated Mars. The model can produce two isotopically distinct planets, but it does not reproduce the specific EC/OC mixture inferred for Mars. This gap must be closed, either by an explicit treatment of an unsampled reservoir, by a different assignment of isotopic reservoirs to the two source rings, or by softening the claim in the abstract and conclusions.","section":"Section 5 and Section 3.2.3"},{"comment":"The Earth's ~70% inner and ~30% outer composition is largely a consequence of the chosen 2:1 mass split (w2 = 1/3 in Eq. (2)), not an independent prediction of the model. The authors state in Section 3.2.1 that the 70/30 chemical split is itself one of the geochemical constraints motivating a two-source model, and Section 5 notes that the 2:1 split 'would represent the best match to the chemical composition of the BSE.' Since model203 was selected as the best among roughly one hundred tested configurations, at least partly with this compositional prior in mind, the claimed match to the Earth's chemistry is partially circular. The paper should quantify the Earth composition for the other mass splits reported in Table 2 (model211 with w2=1/2 and model212 with w2=1/4), and demonstrate that the 70/30 outcome is not simply read off from the initial mass ratio.","section":"Section 4.8, Eq. (2), and Table 2"},{"comment":"The central dynamical mechanism relies on a parameterized gas disk profile, Eq. (9) with beta1 = -1, beta2 = 0, and Sigma0 = 3000 g/cm^2, which produces outward Type-I migration inside ~1.2 au. The authors acknowledge that this is a simple parameterization 'that roughly approximated the results of MHD simulations,' but it is load-bearing: Section 4.6 shows that without this outward migration the inner-ring model does not place Venus and Earth at 0.7-1 au, and Section 4.8 uses the same profile for the reference model. The height of the density bump and the zero-torque radius are not derived from a self-consistent disk model, and the paper does not explore time-dependent beta1 and beta2 even though the cited MHD simulations show the profile evolving. The authors should either provide a more direct connection to published MDW disk calculations, including a range of plausible profiles and evolution histories, or explicitly frame the results as conditional on this assumed disk structure.","section":"Section 2.2.2 and Eq. (9)"}],"minor_comments":[{"comment":"The sentence 'the best results were obtained with 1.4 MEarth in the inner ring and 0.7 MEarth in the inner ring' should read '0.7 MEarth in the outer ring;' otherwise it misstates the mass split that defines model203.","section":"Section 6, Conclusion item 5"},{"comment":"The Notes column describes a gas disk with 'alpha0 = -1,' but the relevant parameter appears to be the power-law index beta in Eq. (8), not alpha0; the notation should be corrected for consistency with Section 4.4.","section":"Table 1, model mmsn"},{"comment":"The affiliation for the Nice group contains a typo: 'Bouldervard de l'Observatoire' should be 'Boulevard de l'Observatoire.'","section":"Title page and affiliations"},{"comment":"The claim that late accretion estimates are unreliable because each initial planetesimal has mass ~0.5% of Earth, comparable to the HSE-inferred late veneer, is clearly stated, but it would be helpful to also state whether the N=1000 and N=2000 resolution runs reported later in the section were used to assess this specific quantity, since those runs are said to give nearly identical results.","section":"Section 4.9"}],"recommendation":"major_revision","confidential_remarks":"The paper is a serious and useful contribution with unusually large simulation statistics, and the authors are honest about many caveats. The main reason for major revision rather than rejection is that the isotopic claim is internally contradicted in Section 5 and the Earth composition claim is partly encoded in the chosen mass split; both are fixable by additional analysis or by reframing the claims. The disk profile concern is a legitimate physical caveat that should be made more prominent, but it is not by itself disqualifying given the paper's exploratory scope."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First, what you should know. This is the most extensive published parameter study of terrestrial planet formation I've seen: roughly 100 model variants at 1000 simulations each, with objective success criteria. The two-source architecture — an inner planetesimal ring near 0.5 au plus an outer source at 1.5–2 au, with outward Type-I migration from a disk-wind-like gas surface density profile — is genuinely new relative to the single-annulus, Grand Tack, and convergent-migration models. The model substantially improves the radial mass concentration and Mercury/Mars mass problems. That part is real and it is the paper's main contribution.\n\nThe soft spot is the paper's own isotope claim. The abstract says the model 'explains the isotopic differences of the Earth and Mars,' but Section 5 admits that in the reference model Mars accretes about 90% of its mass from the outer source, and that 'this composition would not have enough EC material to explain observations.' Geochemical constraints require Mars to be roughly 65% EC. If the inner ring is EC, model203 gives Mars roughly 10% EC. That is an internal contradiction between a headline claim and a caveat written by the same authors. It is not an outside-consensus disagreement. It is a load-bearing failure of one of the three named successes.\n\nTwo more moderate concerns. First, the Earth's 70/30 reduced/oxidized mix follows almost directly from the chosen 2:1 mass split between the reservoirs; it is not an independent prediction in any strong sense. Second, the disk profile in Eq. 9 is a parameterized ad hoc shape. The whole mechanism of leaving Mercury behind and moving Venus/Earth outward relies on it, and it is not derived from a self-consistent MHD disk model. A referee should ask how sensitive the results are to plausible profile variations.\n\nI would not desk-reject this. The dynamical study is careful, the statistics are handled honestly, and the two-source idea is worth testing seriously. The paper is for anyone working on terrestrial planet formation or the dynamical origin of inner solar system architecture. But a referee should demand that the authors either resolve the Mars EC/OC discrepancy with a concrete mechanism or stop claiming the isotopic difference is explained. They should also release code and data, and ideally quantify the selection effect from picking the best of about 100 models. I would send it to review with a request for major revision on the isotopic claim.","headline":"A serious and systematic simulation study with a genuinely new two-source architecture, but the Mars isotope claim is internally contradicted by the paper's own admitted composition numbers.","tokens_in":49571,"tokens_out":2897,"would_cite":true,"duration_ms":34226,"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":"Two separate planetesimal reservoirs — an inner ring near 0.5 au and an outer source at 1.5–2 au — can produce the terrestrial planets, with outward Type-I migration carrying Venus and Earth outward while Mars is the lone survivor of the…","keywords":["terrestrial planet formation","two-source reservoir model","Type-I migration","magnetically driven winds","silicate sublimation line","planetesimal rings","inner solar system","Mars formation"],"falsifier":"A concrete test is to measure or compute the radial gas-density profile of a protoplanetary disk inside 1.2 au and check whether its local slope $\\partial \\ln \\Sigma/\\partial \\ln r$ is positive there. If observations or detailed disk models that include magnetic winds show a monotonically decreasing density with no inner bump, the outward-migration regime disappears; embryos formed near 0.5 au would migrate inward by Lindblad torques and would not end up at the Venus-Earth separation seen in the reference model.","tokens_in":48456,"feed_emoji":"🪐","tokens_out":11419,"duration_ms":121495,"temperature":0.7,"pith_summary":"This paper claims that the terrestrial planets assembled from two distinct planetesimal reservoirs rather than a single radially extended disk. In the reference model (model203), an inner ring at about 0.6 au — where planetesimals would form near the silicate sublimation line — supplies Mercury, Venus, and Earth, while a second, broader source at 1.5–2 au supplies Mars and many Mars-sized bodies that mostly fall onto Earth. Outward Type-I migration driven by a magnetically-wind-like gas disk with a surface-density bump near 1 au carries Venus and Earth to their present orbits and leaves Mercury behind near the ring. The model matches all four terrestrial-planet criteria in 6.1% of 1,000 simulations, produces Earth with about 70% inner-ring (reduced) and 30% outer-source (oxidized) material, and gives Mars a distinct outer-source composition, explaining the Earth-Mars isotopic difference and suggesting Theia shared Earth's accretion history.","feed_headline":"Two planetesimal rings built the inner Solar System","feed_subtitle":"Match all four planets with an inner ring for Mercury-Venus-Earth and an outer source for Mars.","key_machinery":"The load-bearing mechanism is torque-driven migration in a gas disk with a surface-density bump: the paper writes the gas surface density as $\\Sigma_g = \\Sigma_0 (r/r_0)^{\\beta(r)} e^{-t/\\tau}$ with $\\beta(r)=\\beta_1 \\ln(r/r_0)+\\beta_2$, choosing $\\beta_1=-1$, $\\beta_2=0$, $\\Sigma_0=3000$ g cm$^{-2}$, $r_0=0.9$ au. This makes the local slope $\\partial \\ln \\Sigma/\\partial \\ln r$ positive inside about 1.2 au, so protoplanets of roughly 0.1–1 Earth mass experience an outward torque and drift toward the zero-torque radius near 1 au. The inner planetesimal ring at 0.5–0.6 au is the nursery: the most massive embryos grow, migrate out, and concentrate at 0.7–1 au, while lower-mass Mercury stays behind. A second reservoir at 1.5–2 au feeds Mars and, through inward migration, supplies about 30% of Earth's mass later. A giant-planet instability at about 11 Myr breaks the resonant chains formed in the disk and triggers the late giant-impact stage.","core_discovery":"The central discovery is that a two-source initial condition, not a single annulus or an extended disk, best explains the modern terrestrial system. With an inner ring at $r_1 = 0.6$ au ($\\sigma_1 = 0.05$ au), an outer reservoir at $r_2 = 1.7$ au ($\\sigma_2 = 0.1$ au), total initial mass 2.1 Earth masses split 2:1 between inner and outer sources, and a gas disk whose surface density flattens or rises with radius inside about 1.2 au ($\\Sigma_0 = 3000$ g cm$^{-2}$ at $r_0 = 0.9$ au, $\\beta_1 = -1$, $\\beta_2 = 0$), the simulations produce good Venus/Earth pairs in 50% of runs, good Mars in 35%, good Mercury in 32%, and tight Venus/Earth separation in 33%; 6.1% of runs satisfy all four criteria simultaneously. In those runs Earth accretes about 70% of its mass from the inner ring and about 30% from the outer source, while Mars accretes about 90% from the outer source. The paper argues that this accretion sequence matches the reduced-then-oxidized chemistry inferred for Earth's mantle, the distinct isotopic building blocks of Mars, and the Earth-like isotopic composition of the Moon, because Theia would have grown in the same inner reservoir and on a similar timescale as the proto-Earth.","pith_inferences":["A testable extension, not made in the paper: the same disk profile that pushes terrestrial embryos outward near 1 au is also invoked for close-in super-Earth formation, so a unified disk evolution model could tie the terrestrial ring to the super-Earth population and predict their relative formation timescales.","If the inner-ring reservoir is truly unsampled in meteorites, future samples from Mercury or the near-Sun small-body population should carry the reduced, refractory-enriched signature of the inner ring rather than ordinary-chondrite-like material.","The model's rapid assembly of Mars from the 1.5–2 au source implies measurable consequences for Mars' volatile and siderophile budgets; Mars sample-return measurements of tungsten and noble gases could test whether the last large impactor came from that outer reservoir.","The 6.1% all-good success rate is based on N = 400 super-particles; repeating the reference model with several thousand bodies would show whether the stated success probability is stable or partly a resolution artifact."],"forward_implications":["Mercury's small mass and orbit are a natural by-product: the inner ring is narrow and the most massive embryos migrate outward, leaving a low-mass planet behind near 0.5 au.","Earth's mantle composition and the Earth-Mars isotopic difference are explained by the same mass split: about 70% inner-ring (reduced) and 30% outer-source (oxidized) material for Earth, and about 90% outer-source material for Mars.","Theia, the Moon-forming impactor, plausibly grew from the same inner reservoir on a similar schedule as the proto-Earth, explaining the Earth-Moon isotopic similarity without special orbital or mixing assumptions.","Inner-ring planetesimals should be absent from the meteorite collection, because their implantation efficiency into the asteroid belt is estimated to be below about $3 \\times 10^{-7}$.","The best-fit model implies a gas disk with an inner density bump, most naturally produced by magnetically driven winds, and an early giant-planet instability at about 11 Myr; both are testable predictions for disk and planet-formation models."],"supporting_citations":[{"why":"Supplies the narrow-annulus baseline model whose Venus/Earth/Mars successes and Mercury/Venus-separation failures the paper reproduces and then surpasses.","marker":"Hansen (2009)"},{"why":"Defines the radial mass concentration problem and the success metric that the two-source model is built to solve.","marker":"Chambers (2001)"},{"why":"Motivates the inner planetesimal ring by placing early planetesimal formation near the silicate sublimation line.","marker":"Morbidelli et al. (2022)"},{"why":"Supplies the magnetically driven wind disk picture and the Type-I migration and damping prescriptions used in the gas-disk stage.","marker":"Ogihara et al. (2024)"},{"why":"Provides the convergent-migration ring model and giant-instability timing framework that the paper tests against and improves.","marker":"Woo et al. (2024)"},{"why":"Gives the Earth/Mars enstatite- versus ordinary-chondrite building-block constraints that define the two-reservoir interpretation.","marker":"Dauphas et al. (2024)"},{"why":"Establishes the heterogeneous-accretion constraint that Earth accreted about 70% reduced then 30% oxidized material.","marker":"Rubie et al. (2011)"},{"why":"Shows that single-reservoir narrow-ring models cannot match Earth's composition, motivating a second source.","marker":"Dale et al. (2025)"}],"fun_headline_variants":["Earth and Mars trace back to two planetesimal rings","Two rings of rock made Mercury, Venus, Earth, and Mars","Inner ring seeds Earth, outer ring seeds Mars","Two distinct reservoirs forged the inner planets","How two rings of planetesimals built the rocky worlds"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the protoplanetary gas disk had a density profile that pushed newborn planets outward inside about 1.2 au; if real disks instead pushed them inward, Venus and Earth would never move from the 0.5 au nursery out to 0.7–1 au and the inner-ring story fails.","fun_headline_variants_meta":{"raw":{"variants":["Earth and Mars trace back to two planetesimal rings","Two rings of rock made Mercury, Venus, Earth, and Mars","Inner ring seeds Earth, outer ring seeds Mars","Two distinct reservoirs forged the inner planets","How two rings of planetesimals built the rocky worlds"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000293,"raw_usage":{"total_tokens":1818,"prompt_tokens":1166,"completion_tokens":652,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":782,"completion_tokens_details":{"reasoning_tokens":576}},"tokens_in":782,"tokens_out":652,"duration_ms":7915,"temperature":1.0,"reasoning_tokens":576,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T15:47:32.280796+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A concrete test is to measure or compute the radial gas-density profile of a protoplanetary disk inside 1.2 au and check whether its local slope $\\partial \\ln \\Sigma/\\partial \\ln r$ is positive there. If observations or detailed disk models that include magnetic winds show a monotonically decreasing density with no inner bump, the outward-migration regime disappears; embryos formed near 0.5 au would migrate inward by Lindblad torques and would not end up at the Venus-Earth separation seen in the reference model.","supporting_citations":[{"cited_title":"doi:10.3847/1538-4357/ad65d5","cited_arxiv_id":null,"evidence_quote":"Supplies the magnetically driven wind disk picture and the Type-I migration and damping prescriptions used in the gas-disk stage."}],"review_version":1}