{"id":"26ba4362-72bb-4562-a6cc-564ef34d2fed","arxiv_id":"2507.02544","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":11,"one_line_summary":"Callisto's partially differentiated interior is more naturally explained by pebble accretion than by satellitesimal accretion, giving a potential fossil test of planet formation.","lead":"A new model argues that Jupiter's moon Callisto could only have kept its partially differentiated interior if it formed by pebble accretion rather than by accretion of larger satellitesimals. If the JUICE mission confirms the partial differentiation, this would offer a rare observational test of pebble accretion, a leading theory for planet formation.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The entire discrimination rests on the asserted, underived η dichotomy (0 for pebbles vs 0.1–0.3 for satellitesimals); a few percent of subsurface deposition for pebbles would shift Fig. 3 to Cold-like full melting, destroying the claim.","rationale":"I read the paper in good faith and find the central claim clearly conditional: if JUICE confirms partial differentiation with mmelt ∼ 1/3, then the paper argues this would be the first observed evidence for pebble accretion. The argument is orderly and the parameter exploration is broad (τform 0.5–20 Myr, tstart 0.5–10 Myr, six accretion cases, a Ganymede comparison, and a 26Al sensitivity check). The weak point is the separation between the two scenarios: it is entirely the η dichotomy. The reader's weakest_assumption identifies exactly this, and I agree. The paper does not supply a self-contained computation of η for either impactor population; it cites Monteux et al. (2014) for 0.1–0.3 and Bennacer et al. (2025) for η≈0 in pebbles, but the manuscript itself does not show how η scales with impactor size, velocity, porosity, or regolith structure. Because Large and Small pebbles coincide numerically with the η=0 Cold-limit satellitesimal case, the conclusion is a one-parameter sensitivity, not a robust multimodal contrast. The paper also explicitly acknowledges thermal blanketing (Section 4.3) as an omitted effect that could act in the same direction (more melting), but dismisses it without a quantitative bound; this is a secondary, independent concern, but the η issue is more load-bearing because it is the pivot of the entire discrimination. I therefore recommend keeping the reader's CONDITIONAL verdict: the paper is worth publishing as a hypothesis-generating work, but the 'only pebble accretion' claim should be accompanied by a physically grounded η or, at minimum, an η-sweep demonstrating the stability of the mmelt ≤ 1/3 region to the assumed values.","tokens_in":13852,"tokens_out":9535,"duration_ms":103773,"concrete_test":"Recompute the melt-mass-fraction grid of Fig. 3, keeping all equations and parameter ranges unchanged, for pebble accretion with η = 0.01, 0.02, 0.05, 0.08 and for satellitesimal accretion with η = 0.05, 0.08; report the area (in the τform–tstart plane) with mmelt ≤ 1/3. If the pebble η-sweep loses the mmelt ≤ 1/3 region for any η above ~0.03, the 'pebble accretion can maintain partial differentiation' claim fails; if the satellitesimal η-sweep gains such a region at η = 0.08, the exclusivity claim fails. Either outcome would falsify the current version of the central claim.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's two central assertions—'satellitesimal accretion inevitably leads to significant differentiation' and 'pebble accretion can maintain partial differentiation'—are separated by exactly one parameter, η, the fraction of impact energy deposited below the radiative surface. The results in Fig. 3 show that Large/Small pebbles are numerically identical to the η=0 Cold-limit satellitesimal case, so every difference between the scenarios is carried by η. But η is not derived in this manuscript. For satellitesimals it is taken from Monteux et al. (2014) as 0.1–0.3; for pebbles it is asserted to be 0, citing Bennacer et al. (2025) with no reproduction of the size/velocity/porosity dependence. If real pebble impacts deposit even η=0.02–0.05 of their energy below the surface (e.g., through crater ejection or a briefly insulated heated pocket), the pebble panels of Fig. 3 should interpolate toward the Cold-accretion panel, where mmelt>0.5 over the entire grid; the 'can be maintained' part collapses. Symmetrically, if the smallest satellitesimals deposit η<0.1, the 'cannot be achieved' part weakens. Section 4.3 acknowledges thermal blanketing but dismisses it without a bound; that is a second, independent threat to the pebble panel, but the η dichotomy is the more fundamental unvalidated pivot.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper argues that Callisto's partially differentiated interior, inferred from Galileo gravity data, can only be explained by pebble accretion rather than satellitesimal accretion. The author sets up a two-layer energy balance model in which a fraction η of impact kinetic energy is deposited below the surface, uses analytic accretion rates for the two scenarios, and computes melt mass fractions over broad ranges of formation timescales and start times. The central result is that satellitesimal accretion with η=0.1–0.3 leads to widespread melting (mmelt≳0.5), whereas pebble accretion with η=0 can keep mmelt≤1/3 for formation times longer than ~1.2 Myr, matching the observationally inferred partial differentiation.","tokens_in":14238,"tokens_out":3278,"duration_ms":38212,"significance":"If the scenario holds, the paper would provide a novel observational discriminator between accretion mechanisms in a satellite context, with implications for planet formation. The work is clearly structured, uses simple analytic models, and explores a wide parameter grid, including a comparison with Ganymede. It also transparently lists the assumptions and includes an appendix on 26Al abundance. However, the central conclusion rests on a small number of assumed parameters, especially the subsurface energy fraction η, which the paper does not derive from impact physics. The significance is therefore conditional on whether the η dichotomy between pebbles and satellitesimals is physically justified.","major_comments":[{"comment":"The entire discrimination between the accretion scenarios is carried by the assumed subsurface energy fraction η: pebble accretion sets η=0 while satellitesimal accretion sets η=0.1–0.3. As the paper itself notes in Section 3.1, the Large and Small pebble results are identical to the Cold-limit satellitesimal case (η=0) by construction. In Fig. 3, the difference between the 'Cold accretion' panel (η=0.1, mmelt>0.5) and the 'Large pebbles' panel (η=0, mmelt≤1/3) is entirely due to this assumed η contrast. The paper does not derive η from impactor size, velocity, porosity, or target properties, and it does not provide a sensitivity analysis. If pebble impacts deposit even a few percent of their energy below the radiative surface (e.g., η≈0.02–0.05), the pebble panels of Fig. 3 should interpolate toward the Cold accretion result and the central claim that pebble accretion 'can maintain' partial differentiation collapses. The authors should either derive η from microphysical impact models or show that the conclusions are robust to a small nonzero η for pebbles and to the uncertain range for satellitesimals.","section":""},{"comment":"The paper acknowledges that thermal blanketing by a water-vapor atmosphere, other radiogenic isotopes, hydration heat, and long-term convection are ignored, but then dismisses these effects with the statement that 'such additional heat sources do not alter the conclusions of this paper; only pebble accretion has a chance to explain the partially differentiated interior of Callisto.' This is not a valid dismissal for the positive claim that pebble accretion can maintain partial differentiation. If thermal blanketing is effective at Callisto's mass (which is just above the ~0.02 Earth-mass threshold cited from Johansen et al. 2023), the pebble case would experience additional melting and might no longer reproduce mmelt≤1/3. The paper does not provide a bound on the magnitude of these ignored effects for the specific Callisto parameters, so the 'can be maintained' conclusion is not robust to the paper's own acknowledged omissions.","section":""},{"comment":"The thermal model assumes instantaneous radiative cooling at the surface and ignores heat conduction, convection, and latent heat within the satellite. These processes can redistribute heat from the outer accreted layers toward the interior, potentially changing the melt fraction even for fixed η. The paper does not justify that the two-layer energy balance captures the relevant physics for a growing Callisto; it simply adopts the framework. A quantitative estimate of the importance of conduction/convection for the parameters considered, or a comparison with a more complete thermal model, would be needed to support the term 'robustly demonstrate' used in the abstract and conclusions.","section":""}],"minor_comments":[{"comment":"The symbol η is used both for the subsurface energy fraction (Eqs. 1–2) and for the headwind prefactor ηhw in Eq. (11), which is confusing; please use a distinct symbol such as η_hw in the relative velocity expression.","section":""},{"comment":"The text states 'More than 90 wt% of Callisto is differentiated' but mmelt is defined as a mass fraction normalized by the total satellite mass; please use consistent terminology (e.g., 'more than 90% of the mass' rather than 'wt%' if that is intended).","section":""},{"comment":"The panels for 'Large pebbles' and 'Small pebbles' appear identical, which is a consequence of the assumed impact velocities; consider merging or explicitly stating in the caption that they are indistinguishable in this model.","section":""},{"comment":"The citation to Bennacer et al. (2025) is used to justify η=0 for pebbles, but that paper may not have analyzed pebble-sized impactors in the same regime; please clarify the exact basis for this assumption and whether it holds for the pebble sizes considered here.","section":""}],"recommendation":"major_revision","confidential_remarks":"The paper addresses an interesting and potentially high-impact question, and the model is transparent. The main concern is that the central discrimination relies entirely on the assumed value of η for pebbles versus satellitesimals, which is not derived or tested for sensitivity. The dismissal of thermal blanketing and other ignored effects in Section 4.3 is logically insufficient for the positive claim. These issues are fixable with additional analysis (e.g., a sensitivity scan over η, a discussion of impact-energy deposition for pebble-sized projectiles, and a quantitative bound on blanketing effects), so I recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Read it. The paper extends Shibaike et al. (2019) with a broad parameter sweep of accretion heating for Callisto, and it is a useful reference for the JUICE era. The new elements are real: a full tform-tstart grid for six accretion scenarios, a Ganymede comparison, a fragmentation-limited pebble case, and a check on lower 26Al. The model is simple and transparent, and the author is honest about the ignored effects (conductive/convective heat transport, latent heat, vapor-blanketing) and frames the conclusion as a necessary condition. That framing is the right one, because the paper does not prove sufficiency.\n\nThe soft spot is not hidden; it is in Section 2.2. The entire discrimination between satellitesimal and pebble accretion is the subsurface energy fraction eta: 0.1-0.3 for satellitesimals, exactly 0 for pebbles. The Large and Small pebble cases are identical to the Cold-limit satellitesimal case, which also has eta=0. So the 'only pebble accretion can do this' claim reduces to the assertion that pebbles deposit none of their kinetic energy below the radiative surface. That assertion is cited to Bennacer et al. (2025) but not derived, and the Monteux et al. values for satellitesimals are taken as given. A few percent of subsurface energy from pebble impacts would interpolate Fig. 3 toward the Cold panel, and the claimed uniqueness would vanish. The author should either derive eta from impact physics (impactor size, velocity, porosity) or sweep it explicitly for pebbles.\n\nThe second soft spot is thermal blanketing. Section 4.3 acknowledges Johansen et al. (2023) but dismisses it with an argument that the paper only needs necessary conditions. That works for accretion heating, but a vapor atmosphere generated by the same accretion could increase subsurface temperatures and melt the interior, which would directly undermine the pebble result. The author should bound that effect, not just set it aside.\n\nSo: the paper is a solid, clearly written necessary-condition study. It deserves a serious referee, because the diagnostic idea—using Callisto's partial differentiation as a pebble-accretion discriminator—is timely and testable with JUICE. But the 'robustly demonstrate' language in the abstract and summary overstates the case. A referee should ask for an eta sensitivity study and a quantitative treatment of blanketing before the claim can stand.\n\nWho is it for: planetary formation and satellite-interior modelers; anyone prepping JUICE interpretations. I'd bring it to a reading group to discuss the eta issue, but I wouldn't cite it as a proof of pebble accretion. Recommend: send to peer review, but expect major revision.","headline":"Clean parameter study; main 'only pebble accretion' claim rides entirely on an underived subsurface-energy fraction.","tokens_in":14747,"tokens_out":2332,"would_cite":false,"duration_ms":25890,"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":"Callisto's partially differentiated interior can only be explained by pebble accretion, making the moon a possible first observational fingerprint of the mechanism.","keywords":["Callisto","pebble accretion","satellitesimal accretion","Galilean satellites","partial differentiation","accretion heating","planet formation","moment of inertia"],"falsifier":"Measure, by impact experiment or computer simulation, how much kinetic energy centimeter-to-meter ice-rock projectiles deposit beneath the impact point at speeds near Callisto's escape velocity ($\\sim2.4$ km/s); if the buried fraction exceeds a few percent, pebble accretion would heat the subsurface much like satellitesimal accretion and the paper's proposed distinction would collapse.","tokens_in":13598,"feed_emoji":"🧊","tokens_out":20197,"duration_ms":179672,"temperature":0.7,"pith_summary":"The paper argues that Callisto, Jupiter's second-largest moon, can decide between two competing pictures of how large satellites form: accretion of kilometer-sized building blocks called satellitesimals, or accretion of small gas-dragged particles called pebbles. Callisto's gravity data suggest its interior is only partially differentiated, with roughly a third of the mass melted and the rest a mixed ice-rock slab. The paper's thermal model shows that satellitesimal impacts bury 10 to 30 percent of their kinetic energy below the surface, where it cannot radiate away, so satellitesimal accretion differentiates more than half of Callisto for every plausible formation time and start time. Pebble impacts, by contrast, release their energy at the surface and can be slowed by the circumplanetary gas disk, so pebble accretion keeps Callisto partially differentiated over broad parameter ranges. If the upcoming JUICE mission confirms the partial differentiation, this would be the first observed evidence that pebble accretion assembled a large satellite, with consequences for planet formation theory as a whole.","feed_headline":"Callisto's partially mixed interior points to pebble accretion","feed_subtitle":"A Galileo gravity hint, if JUICE confirms it, becomes the first observed evidence for pebble-built moons","key_machinery":"The argument is carried by the subsurface energy-deposition fraction $\\eta$, defined in the paper's energy-balance equations: it is the share of an impactor's kinetic energy deposited beneath the surface, where radiative cooling does not act. Satellitesimal accretion is assigned $\\eta = 0.1$ to $0.3$ following earlier impact studies, while pebble accretion is assigned $\\eta = 0$ because centimeter-to-meter pebbles stop at the surface. The model combines this with an accretion rate $\\dot M \\propto M^{2/3}$ for both mechanisms, impact velocities that include gas-drag-limited settling $v_{\\rm set} = gSt/\\Omega_K$, and $^{26}$Al radiogenic heating, to compute the melt mass fraction $m_{\\rm melt}$ by comparing internal temperature with the pressure-dependent ice melting curve. The $\\eta$ contrast alone is what separates the two scenarios: the large- and small-pebble runs coincide exactly with the $\\eta=0$ 'Cold-limit' satellitesimal case, and only the fragmentation-limited pebbles receive additional protection from reduced impact velocity.","core_discovery":"The central claim is that Callisto's measured normalized moment of inertia, $C/(M_sR_s^2) = 0.3549 \\pm 0.0042$, can be reproduced only if the moon formed by accreting pebbles rather than satellitesimals. In the model, satellitesimal impactors deposit a fraction $\\eta = 0.1$ to $0.3$ of their kinetic energy below the surface, where radiative cooling is ineffective, and this produces melt mass fractions $m_{\\rm melt}$ exceeding 0.5 in the Cold-accretion case and 0.9 in the Hot-accretion case across formation periods from 0.5 to 20 Myr and start times from 0.5 to 10 Myr. Pebble impactors deposit essentially no energy below the surface ($\\eta = 0$), and fragmentation-limited pebbles with Stokes number $St = 0.001$ are further slowed by aerodynamic drag, keeping the melt fraction at or below the observed $m_{\\rm melt} \\sim 1/3$ for formation periods longer than about 1.2 Myr and for late-enough start times. Only the 'Cold-limit' satellitesimal case with $\\eta = 0$, which the author argues is unrealistic for kilometer-sized impactors, avoids differentiation; therefore, if Callisto's partial differentiation is confirmed by JUICE, satellitesimal accretion is ruled out and pebble accretion is the only viable formation mechanism.","pith_inferences":["If the central claim holds, Callisto's measured moment of inertia could be inverted into constraints on pebble size and formation timescale, since the predicted melt fraction depends on both.","The same $\\eta$-based reasoning may apply to other partially differentiated icy bodies, such as medium-sized moons of Saturn or large Kuiper belt objects; finding more examples would show whether pebble-built, surface-heated accretion is a common formation mode rather than a Callisto-specific outcome.","The paper stops at melting during accretion; an editorial extension is to model the long-term thermal evolution of the partially differentiated structure with long-lived radioisotopes and convection to see whether the observed state survives to the present day."],"forward_implications":["If JUICE confirms Callisto is only partially differentiated, satellitesimal accretion is excluded as the formation mechanism and Callisto becomes the first observed case of pebble accretion building a large body.","Ganymede's full differentiation is not produced by pebble accretion in this model, so Ganymede requires additional post-formation heating to explain its fully differentiated state.","For Callisto, the model requires formation lasting at least about 1.2 Myr and starting late enough (roughly $t_{\\rm start} \\gtrsim 2$ Myr) to avoid $^{26}$Al melting; these are concrete, testable conditions on the accretion history.","The critical Stokes number for Callisto is $St_{\\rm crit} \\simeq 0.0086$, so aerodynamic drag matters only for the small, fragmentation-limited pebbles; Ganymede's larger critical value, $St_{\\rm crit} \\simeq 0.02$, means small pebbles are decelerated there more easily, offering a partial handle on why the two moons differ."],"supporting_citations":[{"why":"Provides the Galileo gravity measurement that identifies Callisto as partially differentiated via its normalized moment of inertia.","marker":"J. D. Anderson et al. 2001"},{"why":"Interprets that moment of inertia as a roughly 300 km outer ice shell over a mixed interior, setting the melt-mass-fraction target near one-third.","marker":"G. Schubert et al. 2004"},{"why":"Supplies the eta = 0.1-0.3 range for the fraction of satellitesimal impact energy deposited below the surface.","marker":"J. Monteux et al. 2014"},{"why":"Establishes the prior accretion-heating and aluminum-26 heating baseline and motivates the eta = 0 limiting case.","marker":"A. C. Barr & R. M. Canup 2008"},{"why":"Shows that larger satellitesimal impactors increase interior differentiation and that small pebble-size impactors deposit energy near the surface.","marker":"Y. Bennacer et al. 2025"},{"why":"Provides the earlier pebble-accretion satellite formation model, including the formation-mass threshold, disk parameters, and Stokes numbers used here.","marker":"Y. Shibaike et al. 2019"},{"why":"Supplies the two-dimensional pebble-accretion rate scaling used in the model.","marker":"S. Ida et al. 2016"},{"why":"Gives the pebble relative-velocity formula used to compute impact velocities.","marker":"B. Liu & C. W. Ormel 2018"},{"why":"Describes how JUICE will test hydrostatic equilibrium and refine the differentiation estimate, the observational trigger for the paper's conclusion.","marker":"P. Cappuccio et al. 2022"},{"why":"Provides the fragmentation-limited Stokes number formula used to set the St = 0.001 pebble case.","marker":"S. Okuzumi et al. 2016"}],"fun_headline_variants":["Callisto's partial mix singles out pebble accretion","If confirmed, Callisto is first evidence for pebble moons","Pebble accretion can keep Callisto partially undifferentiated","Callisto's gravity data favors pebble over satellitesimal formation","JUICE could turn Callisto into pebble-accretion proof"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The conclusion rests on the assumption that pebble impacts deposit essentially no kinetic energy below Callisto's surface, while satellitesimal impacts bury 10 to 30 percent of theirs; if real pebbles bury even a few percent at depth, or real satellitesimals bury less than 10 percent, the two formation scenarios would no longer be cleanly separated.","fun_headline_variants_meta":{"raw":{"variants":["Callisto's partial mix singles out pebble accretion","If confirmed, Callisto is first evidence for pebble moons","Pebble accretion can keep Callisto partially undifferentiated","Callisto's gravity data favors pebble over satellitesimal formation","JUICE could turn Callisto into pebble-accretion proof"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000241,"raw_usage":{"total_tokens":1556,"prompt_tokens":1014,"completion_tokens":542,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":630,"completion_tokens_details":{"reasoning_tokens":473}},"tokens_in":630,"tokens_out":542,"duration_ms":5683,"temperature":1.0,"reasoning_tokens":473,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T20:27:33.893922+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure, by impact experiment or computer simulation, how much kinetic energy centimeter-to-meter ice-rock projectiles deposit beneath the impact point at speeds near Callisto's escape velocity ($\\sim2.4$ km/s); if the buried fraction exceeds a few percent, pebble accretion would heat the subsurface much like satellitesimal accretion and the paper's proposed distinction would collapse.","supporting_citations":[{"cited_title":"D., Jacobson , R","cited_arxiv_id":null,"evidence_quote":"Provides the Galileo gravity measurement that identifies Callisto as partially differentiated via its normalized moment of inertia."},{"cited_title":"D., Spohn , T., & McKinnon , W","cited_arxiv_id":null,"evidence_quote":"Interprets that moment of inertia as a roughly 300 km outer ice shell over a mixed interior, setting the melt-mass-fraction target near one-third."},{"cited_title":"2014, title Can large icy moons accrete undifferentiated? , , 237, 377, 10.1016/j.icarus.2014.04.041","cited_arxiv_id":null,"evidence_quote":"Supplies the eta = 0.1-0.3 range for the fraction of satellitesimal impact energy deposited below the surface."},{"cited_title":"C., & Canup , R","cited_arxiv_id":null,"evidence_quote":"Establishes the prior accretion-heating and aluminum-26 heating baseline and motivates the eta = 0 limiting case."},{"cited_title":"2025, title Conditions for Accretion Favoring an Unmelted Callisto and a Differentiated Ganymede , , 6, 138, 10.3847/PSJ/add719","cited_arxiv_id":null,"evidence_quote":"Shows that larger satellitesimal impactors increase interior differentiation and that small pebble-size impactors deposit energy near the surface."},{"cited_title":"2022, title Callisto and Europa Gravity Measurements from JUICE 3GM Experiment Simulation , , 3, 199, 10.3847/PSJ/ac83c4","cited_arxiv_id":null,"evidence_quote":"Describes how JUICE will test hydrostatic equilibrium and refine the differentiation estimate, the observational trigger for the paper's conclusion."},{"cited_title":"2016, title Sintering-induced Dust Ring Formation in Protoplanetary Disks: Application to the HL Tau Disk, The Astrophysical Journal, 821, 82","cited_arxiv_id":null,"evidence_quote":"Provides the fragmentation-limited Stokes number formula used to set the St = 0.001 pebble case."}],"review_version":1}