{"id":"fa2a94de-b265-496f-bce7-62bc2d02b146","arxiv_id":"2505.12652","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Jupiter's radius at the end of the solar nebula was 2.0 to 2.5 times its present value, implying a warm start with entropy about 10.6 to 11 k_B per baryon.","lead":"By combining the orbits of Jupiter's inner moons with the planet's spin history, the authors infer that Jupiter was 2 to 2.5 times wider right after the solar nebula dispersed than it is today. The result gives an independent anchor for when and how Jupiter formed, and what its interior was like.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Equation (1) as printed cancels the ξ/ζ factors and would give R† ≈ 1.4 R_J; the quoted 2.0–2.5 R_J only follows if Ω_E denotes the present-day spin, not the equilibrium spin defined in the text.","rationale":"The most load-bearing step in the paper is not the astrophysical input but the mathematical relation that converts that input into R†. Eq. (1) is displayed as the result of rearranging angular-momentum conservation, and all subsequent claims (entropy, field, accretion rate) are evaluated over the radius range it produces. As printed, the equation is dimensionally consistent but internally cancels the observational factors ξ and ζ, reducing to a value near 1.4 R_J and contradicting the stated 2.0–2.5 R_J. The rest of the paper strongly indicates the intended formula uses the present-day spin rather than the equilibrium spin: the text invokes the constraint that J today is well established, and the 'self-consistent' recomputation gives 2.0–2.56 R_J. This makes a simple typo the most likely explanation, which is why the appropriate verdict remains conditional rather than rejection. However, a reader cannot reproduce the headline number from the printed central equation, and the independent ξ constraint is itself thinly sourced. The proposed check settles the typo question directly and, if the correction is confirmed, leaves the substantive scientific claim intact but in need of a corrected equation and a stronger reference for ξ.","tokens_in":30642,"tokens_out":8521,"duration_ms":89948,"concrete_test":"Re-derive Eq. (1) from J† = J_present with Ω† = χ Ω_br (ζ/ξ)^(3/2). First evaluate with the printed Ω_E (equilibrium spin): R† = sqrt(I_J/I†) R_J ≈ 1.4 R_J. Then evaluate with Ω_E = Ω_J (present-day spin): R† ≈ 2.0–2.5 R_J. Check which version reproduces the refined range 2.0–2.56 R_J obtained by the self-consistent calculation in §3 and the MESA/I† values in Table 1.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central derivation is the angular-momentum conservation relation, Eq. (1). The text defines Ω_E as the disk-locked equilibrium spin, Ω_E = χ (GM_J/R_t^3)^(1/2), and R_t = (ξ/ζ) R_J, so Ω_E/(χ Ω_br) = (ζ/ξ)^(3/2). Multiplying by the printed factor sqrt(ξ^3/ζ^3) gives unity: the ξ and ζ dependence disappears from Eq. (1) and the quoted range is not obtained. Using the same definitions but replacing Ω_E with the present-day spin Ω_J gives R†/R_J = sqrt[(I_J/I†)(Ω_J/(χ Ω_br))(ξ/ζ)^(3/2)], which evaluates to ≈2.2 for the central parameters. Thus the headline result depends on a variable substitution that is not what Eq. (1) states. The §3 self-consistent calculation uses the correct present-day spin and reproduces 2.0–2.56 R_J, so the most plausible reading is a typographical/subscript error in the printed equation; but because Eq. (1) is the paper's central quantitative claim, the manuscript as written does not support the inference. A secondary but related fragility is that the ξ constraint rests on a DPS meeting abstract, with the analytic model in Methods 4.1 requiring a 40% empirical correction; R† scales as ξ^(3/4), so a 10% error in ξ shifts the radius by ~7%.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript proposes an indirect method to determine Jupiter's radius and interior state at the time of proto-solar nebula dissipation. It combines constraints on Io's primordial semi-major axis from the resonant excitation of Amalthea and Thebe's inclinations (ξ = a†_Io/R_J = 4.02–4.98), a disk-locking equilibrium spin for Jupiter (Ω_E = χ√(GM_J/R_t^3)), and angular momentum conservation with MESA interior models that give the primordial moment-of-inertia factor I† as a function of R†. The authors derive R† = 2.02–2.59 R_J, infer a convective-envelope entropy S† ≈ 10.6–11 k_B per baryon, a surface magnetic field B† ≈ 21 mT, and an accretion rate Ṁ = 1.2–2.4 M_J/Myr, and place this state at roughly 3.8 Myr after CAI formation. A self-consistent numerical version of the calculation gives R† = 2.0–2.56 R_J, which the authors adopt as the refined range.","tokens_in":30927,"tokens_out":13716,"duration_ms":130963,"significance":"If the central inference survives scrutiny, the paper is significant: it offers a largely independent, observationally anchored constraint on Jupiter's primordial radius and entropy that bypasses detailed formation modeling, and it produces falsifiable predictions for the primordial field strength and accretion rate. The analytic methods are transparent, the MESA output is publicly available, and the comparison with the Stevenson et al. (2022) diffuse-core model is a useful internal consistency check. The main scientific payoff—a warm-start entropy of roughly 10.6–11 k_B per baryon and a radius of 2–2.5 R_J at disk dissipation—would help discriminate among giant-planet formation scenarios.","major_comments":[{"comment":"The printed equation is algebraically inconsistent with the definitions given in the text. With Ω_E = χ√(GM_J/R_t³) and R_t = (ξ/ζ) R_J, one obtains Ω_E/(χ Ω_br) = (ζ/ξ)^{3/2}, which exactly cancels the printed factor √(ξ³/ζ³). Equation (1) therefore reduces to R†/R_J = √(I_J/I†) ≈ 1.4, not the quoted 2.0–2.6 R_J. The quoted range is recovered only if the ratio Ω_E/(χ Ω_br) is replaced by Ω_J/(χ Ω_br), i.e., if the present-day spin Ω_J enters the angular-momentum conservation statement. The Section 3 self-consistent calculation appears to use the present-day spin and reproduces the range, so the error is likely typographical; nevertheless, because Eq. (1) is the central quantitative claim of the paper, it must be corrected and the variable definitions made explicit and unambiguous.","section":"Section 2, Eq. (1)"},{"comment":"The constraint ξ = a†_Io/R_J = 4.02–4.98 is the linchpin of the analysis, yet it rests on a DPS meeting abstract (Hamilton et al. 2001) and on the analytic model developed here, which requires a 40% empirical reduction to match Amalthea's inclination and gives Thebe's inclination within 7% assuming a particular sequence of resonance crossings. Because R† scales as ξ^{3/4}, a 10% uncertainty in ξ shifts R† by approximately 7%, which is comparable to the width of the quoted radius range. The manuscript should present a quantitative sensitivity analysis of the final radius, entropy, and field predictions to uncertainties in ξ, and should clarify how the 40% correction is propagated through the analytic resonance model.","section":"Section 2 and Methods 4.1"}],"minor_comments":[{"comment":"The sentence 'substituting equation (18) for the gas density' should refer to equation (17). In addition, the expression in Eq. (18) contains the ambiguous term √(2/R_t − √(4/r)), which appears dimensionally inconsistent and should be checked.","section":"Methods 4.5, Eq. (18)"},{"comment":"The primordial spin is denoted inconsistently: Ω_E appears in Section 2 and Eq. (1), while Ω† appears in Section 3 and Methods 4.3. Please use a single symbol for the primordial spin and define it once.","section":"Notation throughout"},{"comment":"The text states that B† ≈ 21 mT should be interpreted as an effective lower bound because the adopted f_ohm and γ are near the lower limits of dynamo models. The abstract and conclusion should carry this qualification, for example by saying 'at least ~21 mT', to avoid over-interpreting the quoted number.","section":"Abstract and Section 3"},{"comment":"The statement that any value of R† exceeding the orbital radius of Amalthea is 'unlikely to be physically meaningful' is not self-evident, because Amalthea's orbit at the epoch of disk dissipation need not equal its present-day orbit; please clarify the reasoning.","section":"Section 2"},{"comment":"The simultaneity of solar-nebula and circum-Jovian-disk dispersal is justified by a one-line energy argument. A brief discussion of possible asynchrony (e.g., different photoevaporation timescales for the solar nebula and a circumplanetary disk) would strengthen the epoch assignment of 3.8 Myr after CAI formation.","section":"Section 2"}],"recommendation":"major_revision","confidential_remarks":"The central equation error is likely fixable and the Section 3 calculation appears to support the quoted radius, so I do not recommend rejection. However, the dependence of the entire result on the unpublished Hamilton et al. (2001) abstract and on a 40% empirical correction in the analytic resonance model should be addressed squarely. The authors should be encouraged to provide a sensitivity analysis for ξ and to correct Eq. (1) with explicit variable definitions before publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague — the headline is this: the paper has a real, publishable idea, and the printed central equation has a typo that makes the main derivation look wrong. The result survives once you replace Ω_E with the present-day spin, and the paper's own §3 iteration confirms it. But as it stands, the equation as written cancels ξ and ζ and gives ~1.4 R_J, not the quoted 2.0–2.5 R_J. That needs to be fixed before anyone cites the formula.\n\nWhat's new: they chain Io's inferred primordial semi-major axis (ξ = 4.02–4.98 R_J from resonance kicks to Amalthea and Thebe) to the magnetospheric truncation radius R_t, then to the disk-locked spin, then to angular momentum conservation. That's a genuinely new diagnostic for the primordial state of a giant planet, and it gives a concrete radius, entropy (~10.6–11 k_B/baryon), field (~21 mT), and accretion rate (1.2–2.4 M_J/Myr). The MESA structure models are standard but the link to the dynamo scaling is sensible. This is a significant within-field result, not a paradigm shift.\n\nThe soft spots are real but narrow. First, Eq. (1): the text defines Ω_E as the equilibrium spin, which already contains the (ζ/ξ)^{3/2} factor; multiplying by √(ξ^3/ζ^3) kills the dependence. The quoted range only appears if Ω_E means the present-day spin. The §3 self-consistent calculation uses the correct present-day value and reproduces 2.0–2.56 R_J, so the right fix is a subscript or typo rather than a conceptual error, but it makes the central printed equation unsupported as written. In a paper whose headline is a number, that's a serious editing failure, not a nit.\n\nSecond, the ξ constraint rests on a DPS meeting abstract (Hamilton et al. 2001). The analytic resonance model in Methods 4.1 matches observations only after a ~40% correction calibrated to those simulations. The paper is upfront about that, but it means the anchor of the whole chain is weaker than it looks. Since R† scales like ξ^{3/4} (or like (ξ/ζ)^{3/2} with the corrected equation), a 10% shift in ξ changes the radius by ~7–8%, comparable to the quoted range's width. So the robustness is modest.\n\nBottom line: the idea is worth refereeing. Send it out, but with a firm instruction to fix Eq. (1) and either get the Hamilton et al. results into published form or present the resonance calculation with enough detail to stand alone. If those two things are handled, this becomes a citable constraint on Jupiter's formation. I'd bring it to reading group; I'd cite the corrected version.","headline":"A genuinely new diagnostic for Jupiter's primordial state, but the printed central equation has a typo that makes the derivation look wrong; the result survives once the present-day spin is used where the equilibrium spin is printed.","tokens_in":31465,"tokens_out":2797,"would_cite":true,"duration_ms":26567,"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":"Jupiter's primordial radius, entropy, magnetic field, and accretion rate are recovered from inner-moon resonances and spin conservation.","keywords":["Jupiter formation","primordial radius","circumplanetary disk","satellite resonances","magnetospheric truncation","angular momentum evolution","warm start entropy","solar nebula lifetime"],"falsifier":"High-precision astrometry of Thebe could decide whether its 1.09-degree inclination requires sequential crossings of the 6:4, 5:3, and 4:2 resonances with Io or could arise from fewer crossings. If the 4:2 crossing is ruled out, the upper bound on Io's primordial orbit disappears and the inferred radius range shifts; if all three crossings are confirmed, the radius range tightens toward its lower end.","tokens_in":30380,"feed_emoji":"🪐","tokens_out":7580,"duration_ms":74246,"temperature":0.7,"pith_summary":"Jupiter's formation is usually modeled from the inside out, with outcomes depending on uncertain details of gas accretion, shock physics, and hydrogen equations of state. This paper sidesteps those uncertainties by reading Jupiter's primordial state from two conserved or recorded quantities: the planet's rotational angular momentum and the orbital resonance histories of its inner moons. It concludes that when the proto-solar nebula dissipated, roughly 3.8 million years after the first solids formed, Jupiter was 2.0 to 2.5 times its present radius, had a convective-envelope entropy of about 10.6 to 11 $k_{\\rm B}$ per baryon, a surface magnetic field near 21 mT, and was accreting circumplanetary disk gas at 1.2 to 2.4 Jupiter masses per million years. These numbers are consistent with core accretion and supply a concrete snapshot of a giant planet at the end of its formation.","feed_headline":"Jupiter was 2–2.5 times current size when the solar nebula vanished","feed_subtitle":"Moon resonances plus spin conservation reveal a warm, fast-accreting Jupiter 3.8 Myr after the first solids.","key_machinery":"The central identity is conservation of rotational angular momentum, $J = I M R^2 \\Omega$, applied between the disk-bearing epoch and the present. The primordial spin is set by disk-locking, an equilibrium between accretion spin-up and magnetic braking that gives $\\Omega^{\\dagger} = \\chi \\sqrt{G M_{\\rm J}/R_{\\rm t}^3}$ with $\\chi \\approx 0.88$. The truncation radius is tied to Io's orbit by $R_{\\rm t} = a_{\\rm Io}^{\\dagger}/\\zeta$, with $\\zeta \\approx 1.13$ from numerical simulations of a three-satellite resonant chain. An integrable Hamiltonian model of second-order mean-motion resonances supplies the constraint on $a_{\\rm Io}^{\\dagger}$ from the inclination kicks to Amalthea and Thebe, and hydrostatic interior models supply the primordial moment-of-inertia factor as a monotonic function of radius. Equation (1) combines these pieces to solve for $R_{\\rm J}^{\\dagger}$.","core_discovery":"The central claim is that Jupiter's radius at the epoch of nebular dissipation can be inferred without invoking a specific accretion model. The paper uses the resonance-excited inclinations of the inner moons Amalthea and Thebe to fix Io's primordial semi-major axis at $a_{\\rm Io}^{\\dagger} = 4.02$ to $4.98\\,R_{\\rm J}$, then converts this to a magnetospheric truncation radius $R_{\\rm t} \\approx 3.6$ to $4.4\\,R_{\\rm J}$ using the established offset of a resonant satellite chain from the disk edge. Because magnetic coupling between Jupiter and the disk locks the primordial spin to a value determined by $R_{\\rm t}$, and because Jupiter's post-nebular evolution conserves rotational angular momentum, the present-day spin and moment of inertia can be run backward to solve for the primordial radius. The result, $R_{\\rm J}^{\\dagger} = 2.02$ to $2.59\\,R_{\\rm J}$, is then used with hydrostatic interior models to infer a 'warm start' entropy of roughly $10.6$ to $11\\,k_{\\rm B}$ per baryon, a dynamo-scaled surface field $B_{\\rm J}^{\\dagger} \\approx 21$ mT, and an accretion rate $\\dot{M} = 1.2$ to $2.4$ Jupiter masses per million years.","pith_inferences":["The same chain of reasoning could be applied to Saturn if its inner moon system preserves a similar resonance footprint, yielding an independent primordial radius and entropy for a second giant planet.","The predicted surface field of about 21 mT puts young Jupiter in a regime where cyclotron maser radio emission might be detectable from analogous young giant planets around nearby stars, offering an observational test of the scaling.","If the inferred accretion rate is representative, the circumplanetary disk must have been supplying matter vigorously right up to the end of the nebular epoch, favoring formation models in which the disk appears late and is short-lived."],"forward_implications":["Jupiter was roughly twice to 2.5 times its present radius at about 3.8 million years after the first solids, which is a direct, datable check on core-accretion formation chronologies.","The inferred envelope entropy of about 10.6 to 11 $k_{\\rm B}$ per baryon places young Jupiter in the 'warm start' regime rather than an extreme hot or cold start.","A primordial surface field near 21 mT, about 50 times today's value, combined with an accretion rate of 1.2 to 2.4 Jupiter masses per million years, constrains the dynamo and mass-feeding history of the circum-Jovian disk.","Io's orbit at disk dissipation was between 4.02 and 4.98 Jupiter radii, so the inner moons' resonance record becomes a usable chronometer for the late stages of giant-planet formation.","The method offers an observationally grounded alternative to accretion-model predictions for the terminal state of a giant planet's formation."],"supporting_citations":[{"why":"Supplies the constraint on Io's primordial semi-major axis from the resonance-excited inclinations of Amalthea and Thebe, which anchors the entire geometric reconstruction.","marker":"[26]"},{"why":"Gives the factor $\\zeta \\approx 1.13$ by which the Io-Europa-Ganymede resonant chain sits outside the disk truncation radius, converting Io's orbit into $R_{\\rm t}$.","marker":"[31]"},{"why":"Provides the magnetic braking and disk-locking framework for the terminal rotation rate of giant planets.","marker":"[32]"},{"why":"Dates solar nebula dissipation to roughly 3.8 million years after CAI formation, fixing the epoch of the inferred Jovian state.","marker":"[47]"},{"why":"Supplies the magnetospheric truncation formula used to relate disk radius, magnetic field strength, and accretion rate.","marker":"[48]"},{"why":"Provides the energy-flux dynamo scaling law used to translate interior heat flux into the primordial surface magnetic field.","marker":"[51]"},{"why":"Gives the dipole-field constant $\\alpha = 1.96$ used in the truncation and spin-equilibrium expressions.","marker":"[53]"},{"why":"Supplies the hydrostatic interior models from which the primordial moment-of-inertia factor and structural coefficient are computed.","marker":"[71]"}],"fun_headline_variants":["Jupiter was 2–2.5 times wider when the nebula vanished","Moon resonances infer Jupiter's swollen primordial state","Young Jupiter had a 50x stronger magnetic field than today","Jupiter's moons reveal its inflated radius and hot start"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"Everything rests on Io having begun its post-nebular tidal migration from a semi-major axis between 4.02 and 4.98 Jupiter radii, as inferred from the resonance excitation of Amalthea's and Thebe's orbital tilts; if that inference is wrong, the derived radius shifts by tens of percent.","fun_headline_variants_meta":{"raw":{"variants":["Jupiter was 2–2.5 times wider when the nebula vanished","Moon resonances infer Jupiter's swollen primordial state","Young Jupiter had a 50x stronger magnetic field than today","Jupiter's moons reveal its inflated radius and hot start"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000853,"raw_usage":{"total_tokens":3785,"prompt_tokens":1100,"completion_tokens":2685,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":716,"completion_tokens_details":{"reasoning_tokens":2626}},"tokens_in":716,"tokens_out":2685,"duration_ms":25381,"temperature":1.0,"reasoning_tokens":2626,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T20:30:16.932006+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"High-precision astrometry of Thebe could decide whether its 1.09-degree inclination requires sequential crossings of the 6:4, 5:3, and 4:2 resonances with Io or could arise from fewer crossings. If the 4:2 crossing is ruled out, the upper bound on Io's primordial orbit disappears and the inferred radius range shifts; if all three crossings are confirmed, the radius range tightens toward its lower end.","supporting_citations":[],"review_version":1}