{"id":"6fee5693-0ab0-4724-835b-ba778d723ab8","arxiv_id":"1908.09840","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":8,"one_line_summary":"The Jovian sweeping secular resonance can excite 50-2000 km planetesimals to eccentricities above 0.6, producing bow shocks that convert 4-9% of dust in the 1.5-3 AU region into chondrules.","lead":"This paper proposes that a shifting resonance with young Jupiter stirred up mid-sized planetesimals as the early gas disk thinned, and the resulting bow shocks melted dust into chondrules. It links Jupiter's formation time and the disk's lifetime to a classic meteorite puzzle.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 4–9% chondrule yield is largely set by assumed planetesimal population inputs: Pc scales with the chosen surface-density normalization (1% of gas) and size distribution, and these are not varied or propagated as uncertainties.","rationale":"The reader's weakest_assumption focuses on the exponential disk-depletion law (Equations 6–7), which is indeed idealized. However, the paper's own parameter study (DEP0p5, DEP2, TDL5, EJ0p1, ST, STR) shows that Pc,tot changes by only a factor of about two across those variations, and the qualitative sweeping-resonance mechanism requires only a monotonic decline of the disk surface density. The more directly load-bearing issue for the quantitative central claim is the planetesimal population normalization: the headline 4–9% is computed after imposing a 1%-of-gas surface density and a particular size distribution, and the paper does not vary or propagate the uncertainty in these choices. The factor-of-two spread between WL2019 and MRN in Table 2 demonstrates that the quoted range is substantially an input choice rather than a derived prediction. This does not undermine the proposed eccentricity-excitation mechanism, which the N-body simulations support, so the verdict remains CONDITIONAL rather than moving to ACCEPT or REJECT. The recommended check would settle whether the quantitative yield survives a plausible downward revision of the normalization and cross-section.","tokens_in":19149,"tokens_out":9933,"duration_ms":99918,"concrete_test":"Recompute Pc,tot in the fiducial N-body model with (i) the planetesimal surface-density normalization reduced from 1% to the MMSN solid/gas ratio (~0.42%) and (ii) σc = πrp^2 for rp < 500 km while keeping π(1.6rp)^2 for larger bodies. If Pc,tot falls below about 2–3%, the 4–9% headline is not robust to these input choices and should be presented as conditional on the planetesimal population model.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The dynamical mechanism — 50–2000 km planetesimals reaching e>0.6 via the sweeping Jovian secular resonance — is supported by both the semi-analytic model and the N-body integrations, and I do not dispute it. The load-bearing weakness is in the quantitative headline. In Section 3.3, after computing the orbital evolution, the simulated planetesimals are scaled to a surface density of 1% of the initial gas disk and to the WL2019 size distribution (Figure 4); these choices enter Equation (32) essentially linearly through the total geometric cross-section of the planetesimal population. Table 2 already shows that switching from WL2019 to MRN size distributions changes Pc,tot by about a factor of two (8.1% vs 4.2%), so the quoted 4–9% range is largely spanned by the choice of size distribution. The 1% normalization itself is never varied; the MMSN solid-to-gas ratio is roughly 0.4%, and the shock cross-section is assumed to be π(1.6rp)^2 for all planetesimal sizes, although the Morris et al. (2012) calibration was for rp > 500 km and the paper acknowledges that smaller bodies likely approach σc ≈ πrp^2. A factor of 2–3 downward in either the normalization or the cross-section would push Pc below the inferred meteoritic range. These uncertainties are not propagated, so the agreement with the observed 1–10% chondrule fraction is not yet a strong test of the sweeping-resonance mechanism itself.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes that, as the protoplanetary disk depletes exponentially, the Jovian sweeping secular resonance excites 50-2000 km planetesimals in the 1.5-3.5 AU region to eccentricities larger than 0.6. The authors argue that bow shocks driven by these highly eccentric planetesimals can heat chondrule precursors, and they compute an average chondrule formation probability of about 4-9% between 1.5 and 3.0 AU. The dynamical mechanism is developed with a semi-analytic secular model and N-body simulations (HERMIT4), including parameter variations in disk depletion timescale, tidal damping strength, Jupiter's eccentricity, and the presence of Saturn. From the model, the authors further infer that the disk depletion timescale was approximately 1 Myr and that Jupiter formed no more than 0.7 Myr after CAIs.","tokens_in":19472,"tokens_out":11491,"duration_ms":106258,"significance":"If the dynamical result holds, this is an attractive and physically motivated solution to the long-standing difficulty of exciting planetesimals to the high eccentricities required for bow-shock chondrule formation. The agreement between the semi-analytic secular treatment and the N-body simulations is a genuine strength, and the paper explicitly models size-dependent aerodynamic and tidal gas drag. The model also produces falsifiable predictions: the disk depletion timescale, the timing of Jupiter's formation, and the size range of the planetesimals responsible for chondrule heating. These are concrete, testable outcomes that go beyond a purely qualitative scenario. The principal weakness is that the headline 4-9% chondrule fraction rests on several adopted scalings that are not varied or propagated as uncertainties, so the quantitative claim is less robust than the dynamical mechanism itself.","major_comments":[{"comment":"The chondrule formation probability Pc is directly proportional to the assumed planetesimal surface density (fixed at 1% of the initial gas surface density in Section 3.3.1) and to the assumed shock cross-section σc = π(1.6rp)^2 applied to all planetesimal sizes. The paper itself notes that the Morris et al. (2012) calibration is for rp > 500 km and that σc should approach πrp^2 for smaller bodies. The WL2019 vs. MRN comparison in Table 2 already shows a factor-of-two sensitivity to the size distribution; combining a more standard MMSN solid-to-gas normalization of ~0.4% with a geometric cross-section for the dominant rp < 500 km population would reduce the fiducial WL2019 value from 8.1% to roughly 2%. This places the quoted range at the edge of, if not below, the observational constraints discussed in Section 3.3.3. The authors should either justify these choices with additional evidence or explicitly propagate their uncertainty into the quoted 4-9% range.","section":"§2.3, Eq. (32), Table 2"},{"comment":"The adopted WL2019 planetesimal size distribution is evaluated at 1 Myr, but in the fiducial model the chondrule-forming resonance sweep occurs at t/τdep ≈ 2.5-4.5 (Figure 7), i.e., at 2.5-4.5 Myr after t=0 for τdep=1 Myr. The planetesimal population can evolve through growth, collisional fragmentation, and dynamical depletion over this interval, so the size distribution at the actual epoch of chondrule formation may differ from the 1 Myr distribution used in Figure 4. The paper does not test this sensitivity; the WL2019/MRN comparison brackets only two static choices, not an evolutionary sequence. This uncertainty directly affects the headline yield and should be addressed with a time-dependent size distribution or a sensitivity test.","section":"§3.3, Figure 4, Figure 7"},{"comment":"The inferred constraints τdep ≈ 1 Myr and 'Jupiter formed no more than 0.7 Myr after CAIs' rely on the assumption that the disk depletes as a global exponential decay with a fixed gap around Jupiter (Eqs. 6-7). While the paper appropriately labels this as an idealized prescription, it does not explore alternative depletion histories, such as inside-out clearing by photoevaporation or a time-variable gap width. Because the timing of the ν5 resonance sweep and the gas density at the time of heating are both consequences of this single depletion law, the quoted timing constraints are conditional on that law. The authors should either show that the main conclusions are robust to a broader class of depletion profiles or soften the statements in Section 4 and the abstract.","section":"§2.1.1, Eqs. (6)-(7), §4"}],"minor_comments":[{"comment":"The DEP0p5 model yields Pc,tot = 3.9% (WL2019) and 1.5% (MRN), both outside the 'about 4-9%' range quoted in the abstract and Section 5. Please qualify the range, for example '4-9% for τdep = 1-2 Myr,' or revise the abstract accordingly.","section":"Table 2 and abstract"},{"comment":"The comparison between the model's Pc (fraction of precursors that are heated) and the observational constraints on the current abundance of chondrules in the asteroid belt and IDPs is not direct, because the observed abundance also depends on survival, transport, and dilution processes. Please clarify this relationship or soften the claim of consistency.","section":"§3.3.3"},{"comment":"The phrase 'where the numerator is evaluated at time t = 0' is ambiguous, since dϖp,J/dt does not depend on the disk density and is constant in this formulation. Consider writing the equation with explicit time dependence to avoid confusion.","section":"Eq. (13)"},{"comment":"The curve labels 'with gap' and 'without gap' are clear, but the y-axis extends to 4 AU while the bulk of the asteroid belt lies inside 3.5 AU; a brief note that the resonance only reaches the belt at late times would aid the reader.","section":"Figure 2"},{"comment":"The condition 'ep,η ≪ 1' should be written as 'ep ≪ 1 and η ≪ 1' to avoid the implication that ep and η appear as a product.","section":"§2.1.2"},{"comment":"The statement that 'the sweeping secular resonance seems to be hardly avoidable' is stronger than the idealized disk model warrants; consider softening it to reflect the dependence on the assumed depletion law.","section":"§4.2"}],"recommendation":"major_revision","confidential_remarks":"The dynamical mechanism is novel and appears sound, and the paper is well within the scope of the journal. The main concern is that the quantitative chondrule fraction is not robust to plausible variations in the planetesimal surface-density normalization and shock cross-section, and the timing constraints depend on a single idealized disk-depletion law. These issues are addressable with additional sensitivity runs and a more careful propagation of uncertainties; I would be willing to review a revised version. The paper's acceptance status at another journal is not relevant to my assessment."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper is worth reading for the dynamical mechanism, not for the headline numbers. The novel piece is not the sweeping secular resonance itself—Zheng et al. (2017) already used it for asteroid belt structure—but the explicit size-dependent excitation in a depleting disk and the attempt to convert that into a chondrule formation probability.\n\nWhat the paper does well: the semi-analytic and N-body models agree on the key outcome, that 50–2000 km planetesimals reach eccentricities above 0.6 when the Jovian secular resonance sweeps inward. The physics of the \"weakly coupled\" size range is clearly laid out, and the parameter variations in Table 2 change Pc by only a factor of about two. The paper is also honest about the idealized disk prescription and about the breakdown of linear theory at high eccentricity. That central dynamical result looks solid.\n\nThe soft spots are where the stress-test note lands. The 4–9% probability is not a robust prediction. It scales linearly with the assumed planetesimal surface density (1% of the initial gas), with the size distribution (WL2019 gives twice MRN), and with the shock cross-section σc ≈ π(1.6rp)^2, which is applied to all sizes even though the Morris et al. (2012) calibration was for rp > 500 km; the paper acknowledges smaller bodies likely approach πrp^2, which would reduce Pc by another factor of 2.6. The MMSN solid-to-gas ratio is roughly 0.4%, so a downward normalization alone could push Pc below the meteoritic range. These uncertainties are not propagated, and the agreement with observed 1–10% chondrule fractions is therefore a model-conditioned consistency check, not an independent test of the mechanism. The inferred τdep ≈ 1 Myr in Section 4.1 is also largely reading back the fiducial input value rather than a new constraint. Minor: the modified HERMIT code and simulation data are not released, which makes it harder to check the N-body details.\n\nThe citation pattern is fine; Nagasawa et al. (2014) is clearly distinguished, and the self-citations are to the prior resonance work that this paper builds on. The central excitation mechanism holds up, and the limitations are addressable in revision.\n\nWho is this for? Planet formation theorists and meteoritics people who care about chondrule heating mechanisms. It deserves a serious referee—the dynamical result will be useful even if the probability estimate needs to be re-framed as conditional on uncertain planetesimal population inputs. I would send it to peer review and would cite the eccentricity-excitation result, with the caveat about the probability calibraation.","headline":"The sweeping-secular-resonance excitation of 50–2000 km planetesimals is a real dynamical result and the paper's best contribution; the 4–9% chondrule yield is a model-conditioned estimate, not a measurement.","tokens_in":20063,"tokens_out":1462,"would_cite":true,"duration_ms":16788,"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 Jovian sweeping secular resonance, moving inward as the disk depleted, drove 50-2000 km planetesimals to eccentricities above 0.6 and converted 4-9% of chondrule precursors between 1.5 and 3 AU into chondrules.","keywords":["chondrules","Jovian sweeping secular resonance","planetesimal eccentricity excitation","protoplanetary disk depletion","bow shock heating","asteroid belt","Jupiter formation","gas drag"],"falsifier":"Measure isotopic ages of the oldest chondrules and compare with the time Jupiter opened its gap: the model requires Jupiter to already be massive before the first chondrules, no later than about 0.7 Myr after the CAIs, and requires chondrule production to be spread over roughly 1-2 Myr as the resonance sweeps the belt. A chondrule dated earlier than Jupiter's gap opening, or an age spread much shorter than about 1 Myr, would rule the scenario out.","tokens_in":18888,"feed_emoji":"☄️","tokens_out":6434,"duration_ms":62683,"temperature":0.7,"pith_summary":"This paper argues that the same event that drained the solar nebula also made chondrules. As the gas disk's surface density decayed, Jupiter's secular resonance swept inward through the asteroid belt, pumping the eccentricities of intermediate-size planetesimals and turning them into bow-shock heating sources. In the authors' fiducial model, planetesimals 50-2000 km across reach eccentricities above 0.6, and their shocks convert 4-9% of chondrule precursors between 1.5 and 3 AU into chondrules. If correct, this ties chondrule ages to disk dispersal: the spread of chondrule formation times implies a disk depletion timescale near 1 Myr, and Jupiter must already have formed, within about 0.7 Myr after the calcium-aluminum inclusions.","feed_headline":"Jupiter's sweeping resonance melted 4-9% of chondrule precursors","feed_subtitle":"As the gas disk drained, 50-2000 km planetesimals were flung into orbits that drove chondrule-forming bow shocks.","key_machinery":"The load-bearing object is the Jovian sweeping secular resonance, the location in the disk where a planetesimal's apsidal precession rate matches Jupiter's, so that Jupiter's gravitational pull coherently stretches the orbit over many cycles; that location migrates inward as the disk surface density decays exponentially with timescale τdep. On either side of the resonance, gas drag acts as a size filter: aerodynamic drag dominates for planetesimals below about 50 km, tidal (Lindblad) drag dominates above about 2000 km, and the weakly coupled intermediate sizes are free to reach high eccentricity. The same drag laws set the speed of the sweep and the gas density during heating, and the probability calculation combines the bow-shock cross-section with a settled dust scale height to turn each high-eccentricity encounter into a chondrule-forming event.","core_discovery":"The central discovery is a single causal chain. After Jupiter opens a gap, the location where Jupiter's apsidal precession rate matches a planetesimal's precession rate, the secular resonance called ν5, moves inward as the disk mass falls. Planetesimals in the 50-2000 km size range are weakly coupled to the gas: small ones are damped by aerodynamic drag, and larger ones migrate inward too fast under tidal drag, but this middle group can be excited to eccentricities above 0.6 before gas drag stops them. Their bow shocks then heat chondrule precursors; the model finds most heating occurs in high-velocity shocks (vrel ≈ 12-18 km/s), in gas depleted to 1-10% of the minimum-mass solar nebula, and inside about 2.5 AU, giving a 4-9% probability that a precursor between 1.5 and 3 AU becomes a chondrule. The same sweep implies the disk depleted on a timescale of about 1 Myr and that Jupiter formed before the chondrules, no more than 0.7 Myr after the CAIs.","pith_inferences":["If the mechanism is robust, the same resonance sweep that made chondrules should have also shaped the asteroid belt's dynamical structure and mass deficit, so a combined model of chondrule yield and belt clearing could sharpen the timing predictions.","A straightforward test is to rerun the sweep with photoevaporation-driven depletion instead of pure exponential decay; if the resonance then passes through the belt faster, the predicted chondrule fraction and age spread would change in ways that isotopic chronologies could check.","The model's concentration of heating inside about 2.5 AU predicts a radial gradient in chondrule flash-heating signatures, so chondrules from parent bodies that formed farther out should show weaker or absent chondrule textures."],"forward_implications":["The 4-9% chondrule fraction matches the observationally estimated fraction of the present asteroid belt, so no separate heating mechanism is required once Jupiter formed.","Chondrule ages should trace the resonance sweep: production is spread over roughly 1-2 disk depletion timescales, with most chondrules forming 1-3 Myr after the CAIs.","Jupiter must have been massive enough to open a gap before the first chondrules, placing its formation no later than about 0.7 Myr after the CAIs.","The solar nebula's depletion timescale is about 1 Myr, consistent with inferred disk lifetimes around young stars.","Chondrule-producing planetesimals are confined to roughly 50-2000 km in radius; smaller and larger bodies are damped or migrate away before they can heat precursors."],"supporting_citations":[{"why":"Introduced the sweeping secular resonance concept that the paper adapts to a depleting protoplanetary disk.","marker":"Ward et al. (1976)"},{"why":"Supplies the disk-precession formulation, the planetesimal density prescription, and the late-stage application the paper builds upon.","marker":"Zheng et al. (2017)"},{"why":"Previous secular-resonance chondrule model with a fixed disk mass, the immediate precursor this paper generalizes to a depleting disk.","marker":"Nagasawa et al. (2014)"},{"why":"Provides the gas-density-dependent velocity window for chondrule melting used to define which shocks count as chondrule-forming.","marker":"Iida et al. (2001)"},{"why":"Gives the bow-shock cross-section used in the probability integral for chondrule heating.","marker":"Morris et al. (2012)"},{"why":"Supplies the MMSN surface density and temperature profile adopted for the disk model.","marker":"Hayashi (1981)"},{"why":"Provides the lowest-order aerodynamic drag formulae describing eccentricity and semi-major axis damping.","marker":"Adachi et al. (1976)"},{"why":"Gives the planetesimal size distribution used in the main chondrule probability estimates.","marker":"Walsh & Levison (2019)"}],"fun_headline_variants":["Jupiter's resonance melted 4-9% of chondrules","Sweeping resonance shocks formed meteorite beads","Jupiter's gap sweep heated planetesimals to melt","Jupiter's secular resonance drove chondrule bow shocks","How Jupiter's sweep melted early meteorite spheroids"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The results depend on the disk losing mass as a single global exponential decay of an MMSN-like surface density with a fixed gap around Jupiter, with no photoevaporation, radial transport, or gap evolution included.","fun_headline_variants_meta":{"raw":{"variants":["Jupiter's resonance melted 4-9% of chondrules","Sweeping resonance shocks formed meteorite beads","Jupiter's gap sweep heated planetesimals to melt","Jupiter's secular resonance drove chondrule bow shocks","How Jupiter's sweep melted early meteorite spheroids"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000279,"raw_usage":{"total_tokens":1731,"prompt_tokens":1090,"completion_tokens":641,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":706,"completion_tokens_details":{"reasoning_tokens":561}},"tokens_in":706,"tokens_out":641,"duration_ms":6436,"temperature":1.0,"reasoning_tokens":561,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T11:00:44.494653+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure isotopic ages of the oldest chondrules and compare with the time Jupiter opened its gap: the model requires Jupiter to already be massive before the first chondrules, no later than about 0.7 Myr after the CAIs, and requires chondrule production to be spread over roughly 1-2 Myr as the resonance sweeps the belt. A chondrule dated earlier than Jupiter's gap opening, or an age spread much shorter than about 1 Myr, would rule the scenario out.","supporting_citations":[{"cited_title":"R., Colombo, G., & Franklin, F","cited_arxiv_id":null,"evidence_quote":"Introduced the sweeping secular resonance concept that the paper adapts to a depleting protoplanetary disk."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the disk-precession formulation, the planetesimal density prescription, and the late-stage application the paper builds upon."},{"cited_title":"K., Tanaka, H., et al","cited_arxiv_id":null,"evidence_quote":"Previous secular-resonance chondrule model with a fixed disk mass, the immediate precursor this paper generalizes to a depleting disk."},{"cited_title":"2001, Icarus, 153, 430","cited_arxiv_id":null,"evidence_quote":"Provides the gas-density-dependent velocity window for chondrule melting used to define which shocks count as chondrule-forming."},{"cited_title":"A., Boley, A","cited_arxiv_id":null,"evidence_quote":"Gives the bow-shock cross-section used in the probability integral for chondrule heating."},{"cited_title":"J., & Levison, H","cited_arxiv_id":null,"evidence_quote":"Gives the planetesimal size distribution used in the main chondrule probability estimates."}],"review_version":1}