{"id":"dbe31813-e867-4d16-a5f4-318b86f6895a","arxiv_id":"2411.14709","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"The Moon-forming giant impact is shown to be extremely unlikely and too thoroughly mixing under pebble accretion, making that model less favorable for the inner solar system.","lead":"This paper tests whether the leading 'pebble accretion' model for forming the inner planets can also explain the Moon. Using planet-formation and giant-impact simulations, it finds the required Moon-forming collision is rare and would erase observed Earth-Moon differences.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Unmodeled uniform prior over initial separation Δ makes the <1‰ collision probability a conditional statement, not a model prediction.","rationale":"The reader's weakest-assumption analysis identifies the same load-bearing weakness: the unmodeled prior over initial separation Δ. This is the most consequential issue because the paper's central quantitative claim—that the Moon-forming impact occurs with probability below 1‰ in the pebble-accretion scenario—is conditional on an arbitrary uniform distribution over Δ. Without deriving this distribution from the pebble-accretion model, the probability is not a robust prediction. The geochemical mixing leg is secondary: even if the probability were higher, the mixing argument would still need the separate assumption that proto-Earth and Theia are isotopically closer than Earth and aubrites (Section 3.2), which is also unmodeled. However, the probability leg is the headline and the basis for the paper's main conclusion. I credit the paper for using public codes (REBOUND, GIZMO), running a large number of simulations, and providing clear hydrodynamical evidence for thorough mixing. But because the prior issue is fixable by an additional derivation and does not invalidate the numerical methods, a CONDITIONAL verdict is appropriate—unchanged from the reader's assessment.","tokens_in":11768,"tokens_out":6989,"duration_ms":68357,"concrete_test":"Re-run the N-body simulations with initial separations sampled from the distribution produced by a full pebble-accretion simulation (including migration and growth), rather than a uniform prior over Δ. If the resulting probability of a delayed, correctly angled collision remains below 1‰ under this physically motivated distribution, the concern is resolved. A cheaper check: take the existing 2900 simulation outcomes and re-weight them by the physical Δ distribution, then recompute the quoted probability.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central probability claim rests on an arbitrary initial condition. In Section 2.1, proto-Earth–Theia pairs are initialized with separations Δ uniformly distributed in [1,30] mutual Hill radii (2300 runs over [1,21], 600 over (21,31]). The paper does not derive this distribution from the pebble-accretion model; it starts at t=5 Myr with mass growth turned off and simply samples Δ uniformly. The headline 'probability <1‰' (Section 3.1, Fig. 2a, abstract) is therefore the fraction of simulations that, under this hand-chosen prior, produce a delayed collision at the right time and angle. It is not a probability predicted by the pebble-accretion scenario. If the physical Δ distribution from a self-consistent pebble-accretion simulation (e.g., Johansen et al. 2021) is peaked at small separations or at the specific wide pair envisaged there, the quoted probability could shift by orders of magnitude. The conclusion even concedes 'uncertainty in how they get there.' This is the single most load-bearing weakness because the abstract's main quantitative claim and the 'disfavors pebble accretion' conclusion both depend on it.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper tests the pebble-accretion scenario for terrestrial-planet formation against the Moon-forming impact. It initializes 2900 N-body simulations of proto-Earth–Theia pairs on the Johansen et al. (2021) mass-growth track, samples initial separations uniformly in Hill radii, and counts collisions by timing and angle; it then runs GIZMO/MFM giant-impact simulations for representative collisions and quantifies mixing of target and impactor material in the post-impact planet and protolunar disk. The central claims are that a Moon-forming impact in this scenario has probability below 1 per mille of occurring at the right timing (70–120 Myr) and configuration, and that, if it does occur, perfect mixing leaves no room for primordial mantle heterogeneities or the observed small Earth–Moon isotopic difference. On this basis the authors conclude that pebble accretion is less favorable than chaotic collisional growth for the inner solar system.","tokens_in":12011,"tokens_out":7100,"duration_ms":70269,"significance":"The paper addresses a real and timely question and uses a plausible falsification strategy: the Moon is a sensitive, rarely considered constraint on pebble accretion. It deserves credit for combining N-body statistics with state-of-the-art MFM giant-impact simulations, for testing a range of disk, Jupiter, and migration variants, and for making an explicit falsifiable prediction about isotopic mixing in the protolunar disk. If the statistical claims were robust, this would be a significant argument against fast pebble-accretion formation of the terrestrial planets. However, the headline probability is conditional on an ad hoc uniform prior over initial separations and on a very small number of right-timing events, so the conclusion currently overstates the model-disfavoring power.","major_comments":[{"comment":"The quoted \"<1 per mille\" probability is computed by counting outcomes over a hand-chosen uniform sampling of initial separations Δ in [1,30] Hill radii (2300 runs over [1,21] and 600 over (21,31]). The pebble-accretion model is not used to derive this distribution; the simulation starts at t=5 Myr with mass growth turned off. The probability of a delayed collision is therefore a conditional frequency given the adopted prior, not a prediction of the pebble-accretion scenario. The authors' own conclusion in Section 4 (“Let alone the uncertainty in how they get there”) concedes this limitation. Please derive the Δ distribution from a self-consistent pebble-accretion calculation or clearly report all probabilities as conditional on Δ and test a range of physically motivated priors; without this, the headline probability is not a robust constraint.","section":"Section 2.1 and Section 3.1"},{"comment":"The right-timing count is 4 of 2900 simulations, yet the paper quotes \"probability <1 per mille\" in the abstract and Section 3.1 without statistical uncertainty. Under Poisson statistics, 4 events in 2900 trials gives a 95% confidence interval of roughly 1–10 events, i.e., about 0.3–3.5 per mille of trials, so the upper bound exceeds the quoted <1 per mille threshold. The small count also makes the subsequent \"1/3 angle selection\" factor fragile. Please report the count with a confidence interval and either increase the simulation count or soften the <1 per mille claim to match the statistical precision.","section":"Section 3.1"},{"comment":"The \"right configuration\" criterion is not fully specified. Table 1 lists many successful runs with initial or final angular momentum greater than 2 L_EM, and the authors themselves note that “it is controversial whether this excessive angular momentum can be removed (Rufu & Canup 2020).” If the Moon-forming requirement is tightened to avoid this uncertainty, the \"1/3 of delayed collisions\" factor and hence the <1 per mille estimate change; if it is not tightened, the mixing and disk-mass analysis may not correspond to the actual lunar-forming subset. Please specify the success criterion used to classify runs as Moon-forming in terms of disk mass, iron fraction, angular momentum, and angular-momentum-removal assumptions, and propagate its uncertainty into the probability estimate.","section":"Section 3.2 and Table 1"},{"comment":"The geochemical argument assumes, without model derivation, that proto-Earth and Theia in the pebble-accretion scenario have a smaller initial Δ17O difference than the Earth–aubrites difference of 22 ppm. The 695 ppm initial-difference calculation for run 20 is an illustrative estimate, but the claim that the observed 22 ppm Earth–Moon difference is incompatible depends on this unmodeled compositional prior. If proto-Earth and Theia can acquire a larger initial isotopic difference through heterogeneous pebble accretion, the mixing implied by the simulations may not be in tension with the observations. Please either derive the expected proto-Earth–Theia composition difference from the pebble-accretion model or present the geochemical conclusion as a conditional statement.","section":"Section 3.2"}],"minor_comments":[{"comment":"The five groups of test simulations and the three migrating-giant setups are described in prose; a summary table of the disk, eccentricity, and Jupiter parameters for each group would improve reproducibility.","section":"Section 2.1"},{"comment":"There is a typo in the phrase “ANEOSequation ofstate”; it should read “ANEOS equation of state.”","section":"Section 2.2"},{"comment":"The text states that impact angles are almost uniformly distributed, but the definition of impact angle and its relation to the impact parameter b used in Table 1 should be stated explicitly so the reader can connect Fig. 2b to the giant-impact runs.","section":"Figure 2b"},{"comment":"The statement that the collision probability for pairs with Δ>10 is “always minimal” is vague; please report the underlying counts for each test group so the reader can assess the strength of that conclusion.","section":"Section 3.1 and Figure 3"},{"comment":"The phrase “Let alone the uncertainty in how they get there” undermines the strength of the conclusion as written; consider moving this caveat into the methods or quantifying its effect on the probability estimate.","section":"Section 4"},{"comment":"The statement that data files “will be made available upon reasonable request” is weaker than the standard for a numerical study of this kind; depositing initial conditions and processed outputs in a public repository would increase confidence in the results.","section":"Data Availability"}],"recommendation":"major_revision","confidential_remarks":"The authors' own concession in Section 4 is honest, but it means the abstract's '<1 per mille' claim is not yet substantiated as a model prediction. I would encourage the editor to require either a self-consistent derivation of the initial-separation distribution from the pebble-accretion model or an explicit reframing of the results as conditional on that distribution. The paper is within MNRAS scope and has no apparent ethical concerns."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The thing to know about this paper is that it is the first concrete test of whether the pebble-accretion scenario for the terrestrial planets can actually produce a Moon-forming impact. The giant-impact mixing analysis is the strongest part: well-resolved MFM simulations show that half-Earth impacts at near-escape velocity thoroughly mix target and impactor, including the lower mantle, and put about 50% Theia material in the disk. That result is clean and credible. The N-body suite is also a real effort, with 2900 fiducial runs and several parameter variations, and it fills a gap that Johansen et al. explicitly left open.\n\nThe soft spot is where I would push the authors. The headline probability, <1‰, is a fraction of simulations drawn from a uniform prior over initial separations Δ. That prior is not derived from the pebble-accretion model; the paper starts at t=5 Myr, turns off mass growth, and simply samples Δ in [1,30] Hill radii. So the number is conditional on a hand-chosen distribution. The authors actually concede this in the conclusions with 'let alone the uncertainty in how they get there.' If the physical Δ distribution is peaked elsewhere, the probability could shift by orders of magnitude. That does not kill the paper, but it means the abstract's 'rendering it less favourable' is currently too strong. The fix is straightforward: derive the separation distribution from a self-consistent pebble-accretion simulation, or at least present the collision probability as a function of Δ so readers can apply their own prior. The other statistical issue is that the key count is 4 events out of 2900, with no Poisson error bars; that small number matters when you claim something is below 1‰.\n\nThe geochemical argument is plausible but more speculative. The mixing result itself is solid, but the claim that proto-Earth and Theia must be isotopically closer than Earth and aubrites is assumed rather than derived. It is a secondary point, not the load-bearing part.\n\nAll together, this is a serious paper with honest self-assessment and a genuinely new test. The central conclusion is plausible but not fully established, because the main probability claim depends on an unmodeled prior. I would accept it for peer review and would cite it, since it is the first real attempt to confront the pebble-accretion scenario with the Moon's origin. A referee should ask for a revised, more carefully conditioned probability statement and error bars on the 4/2900 count before this can be considered a strong disfavor of pebble accretion.","headline":"A credible and novel test of the pebble-accretion scenario via the Moon-forming impact, but the headline probability rests on an unmodeled prior and the strong conclusion needs to be reined in.","tokens_in":12551,"tokens_out":2682,"would_cite":true,"duration_ms":51925,"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 Moon-forming impact is so rare and so thoroughly mixing under pebble accretion that the model cannot explain the Earth-Moon system.","keywords":["Moon-forming impact","pebble accretion","terrestrial planet formation","giant impact","proto-Earth","Theia","primordial mantle heterogeneity","oxygen isotopes"],"falsifier":"Compute the orbital separation distribution of proto-Earth and Theia from a full pebble-accretion growth simulation; if a substantial fraction of such pairs evolve into the window for delayed collisions (2√3 to 16 Hill radii), the quoted sub-per-mille collision probability would not hold.","tokens_in":11572,"feed_emoji":"🌙","tokens_out":6349,"duration_ms":56553,"temperature":0.7,"pith_summary":"The paper argues that the Moon's origin, which has been thoroughly explored in the classical planetesimal-accretion picture, can discriminate between that picture and the alternative pebble-accretion scenario. It combines long-term N-body simulations of the proto-Earth and Theia pair with high-resolution giant-impact simulations to show that, under pebble accretion, a suitable Moon-forming impact is a rare event with probability below one per mille. Even when such an impact occurs, it mixes the two bodies so completely that it would erase the primordial mantle heterogeneities and the small Earth-Moon compositional difference that are observed. The paper concludes that the Earth-Moon system and the other terrestrial planets should preferably form via chaotic collisional growth among planetesimals, not via pebble accretion.","feed_headline":"Pebble accretion can't make the Moon, simulations show","feed_subtitle":"Earth-Theia collisions are too rare and too well mixed to match lunar geochemistry.","key_machinery":"The discriminator is the Moon-forming giant impact between proto-Earth (about 0.6 Earth masses) and Theia (about 0.4 Earth masses), which in the pebble accretion model is the only way to make the Moon. The quantitative machinery is a pair of simulations: N-body integrations that track whether the two bodies collide, when, and with what impact parameters, and global giant-impact simulations that follow the mixing of target and impactor material through the proto-lunar disk. A key intermediate is the Hill-radius separation Δ between the two bodies, which controls whether collisions occur quickly (Δ < 2√3) or are delayed (Δ between 2√3 and 16 Hill radii). The paper uses this separation to convert collision frequencies into probability statements about the Moon-forming event.","core_discovery":"The central claim is that the pebble accretion paradigm for terrestrial planet formation is strongly disfavored because it cannot reproduce the Moon-forming impact. In 2900 N-body simulations initialized with proto-Earth and Theia separated by 1 to 30 mutual Hill radii, collisions between the two half-Earths are common only when the separation is less than 2√3 Hill radii, in which case they occur within 1 Myr; for wider separations the collision probability drops to 1.5%, and only four simulations produce a giant impact within the 70–120 Myr window inferred for the Moon's formation. When such delayed impacts do occur, high-resolution hydrodynamic simulations show that the target and impactor mix almost completely, with 48–55% of the protolunar disk coming from Theia, erasing any primordial mantle reservoirs and forcing an initial oxygen isotope difference between proto-Earth and Theia far larger than expected for bodies with nearly identical accretion histories. The paper's conclusion is that the Moon's origin favours the classical collisional-growth scenario over pebble accretion.","pith_inferences":["A natural extension is that if pebble accretion is disfavored for the terrestrial planets, Earth's formation timescale would revert to tens of millions of years, bringing into question pebble-accretion-based interpretations of early isotopic chronometers.","The same half-Earth mixing argument could be applied to any proposed giant impact on Venus or Mars, predicting that such impacts, if they occurred, would also erase primordial heterogeneities in those bodies.","The statistical test could be sharpened by simulating the full pebble-accretion assembly of proto-Earth and Theia to derive a self-consistent initial orbital separation distribution instead of sampling it uniformly."],"forward_implications":["Under pebble accretion, the probability of a Moon-forming impact at the right timing and angle is below one per mille, so the model requires extreme fine-tuning to explain the Earth-Moon system.","A successful half-Earth impact produces a protolunar disk with 48–55% of its material from Theia and a thoroughly mixed mantle, so any primordial Earth mantle heterogeneity would not survive.","To match the observed 22 ppm oxygen isotope offset after such mixing, proto-Earth and Theia would need an initial isotopic difference of about 695 ppm, even larger than the Earth-Mars difference, which is implausible for close siblings formed by pebble accretion.","Delayed Moon-forming impacts in the 70–120 Myr window occur only for initial separations between 2√3 and 16 Hill radii, and the collision probability in that range is below three per mille."],"supporting_citations":[{"why":"Supplies the pebble-accretion model of terrestrial planet formation that this paper sets out to test.","marker":"Johansen et al. 2021"},{"why":"Provides the 70–120 Myr window for the Moon-forming impact used as the timing constraint.","marker":"Halliday & Canup 2023"},{"why":"Establishes the canonical giant impact scenario for the Moon that the classical model builds on.","marker":"Canup & Asphaug 2001"},{"why":"Gives the empirical relation between protolunar disk mass and final Moon mass used to classify successful impacts.","marker":"Canup 2012"},{"why":"Presents meshless-finite-mass giant impact simulations showing a layered post-impact proto-Earth, the basis of the mixing analysis.","marker":"Deng et al. 2019b"},{"why":"Defines the Hill-stability criterion Δ<2√3 used to separate fast from delayed collisions.","marker":"Gladman 1993"},{"why":"Measures the 22 ppm oxygen isotope offset between Earth and the Moon used for the mass balance calculation.","marker":"Cano et al. 2020"},{"why":"Reports primordial mantle heterogeneities that the paper's mixing argument says would be erased by the impact.","marker":"Mukhopadhyay 2012"},{"why":"Provides further isotopic evidence for reservoirs that predate the Moon-forming impact.","marker":"Rizo et al. 2016"}],"fun_headline_variants":["Pebble accretion can't make the Moon, new simulations show","Moon-forming impact is too rare and too mixed for pebble accretion","Pebble accretion fails to explain the Moon's geochemical fingerprint","Moon birth disfavors pebble accretion, favoring chaotic growth","Rare, perfectly mixed collisions make pebble accretion unlikely for Moon"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The results assume that the initial orbital separations between proto-Earth and Theia are distributed uniformly from 1 to 30 mutual Hill radii, a distribution that is adopted rather than derived from the pebble accretion model.","fun_headline_variants_meta":{"raw":{"variants":["Pebble accretion can't make the Moon, new simulations show","Moon-forming impact is too rare and too mixed for pebble accretion","Pebble accretion fails to explain the Moon's geochemical fingerprint","Moon birth disfavors pebble accretion, favoring chaotic growth","Rare, perfectly mixed collisions make pebble accretion unlikely for Moon"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001139,"raw_usage":{"total_tokens":4733,"prompt_tokens":955,"completion_tokens":3778,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":571,"completion_tokens_details":{"reasoning_tokens":3688}},"tokens_in":571,"tokens_out":3778,"duration_ms":26919,"temperature":1.0,"reasoning_tokens":3688,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T14:59:25.751008+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the orbital separation distribution of proto-Earth and Theia from a full pebble-accretion growth simulation; if a substantial fraction of such pairs evolve into the window for delayed collisions (2√3 to 16 Hill radii), the quoted sub-per-mille collision probability would not hold.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the pebble-accretion model of terrestrial planet formation that this paper sets out to test."},{"cited_title":"N., Canup R","cited_arxiv_id":null,"evidence_quote":"Provides the 70–120 Myr window for the Moon-forming impact used as the timing constraint."},{"cited_title":"M., 2012, Science, 338, 1052","cited_arxiv_id":null,"evidence_quote":"Gives the empirical relation between protolunar disk mass and final Moon mass used to classify successful impacts."},{"cited_title":"J., Sharp Z","cited_arxiv_id":null,"evidence_quote":"Measures the 22 ppm oxygen isotope offset between Earth and the Moon used for the mass balance calculation."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Reports primordial mantle heterogeneities that the paper's mixing argument says would be erased by the impact."},{"cited_title":"J., Carlson R","cited_arxiv_id":null,"evidence_quote":"Provides further isotopic evidence for reservoirs that predate the Moon-forming impact."}],"review_version":1}