{"id":"0f34919f-309c-4de3-9a5a-774e3e8bb849","arxiv_id":"2411.13351","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":8,"one_line_summary":"A model of Jupiter's early circumplanetary disk shows that self-shadowing creates a ~100 K cold trap around 10 Jupiter radii that can condense volatile ices and may shape Galilean moon compositions.","lead":"This paper models the gas disk around young Jupiter with a 2D quasi-static model and finds a self-shadowed zone about 10 Jupiter radii out that is roughly 100 K colder than its surroundings. That cold trap could let ices like ammonia, carbon dioxide and hydrogen sulfide condense closer to Jupiter, potentially affecting the composition of the Galilean moons.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. 26 multiplies the vertical temperature gradient by an extra factor q_s; until corrected, the 100 K shadowed cold trap is not reproducible from the manuscript's equations.","rationale":"The reader's weakest-assumption pick was the accretion-rate prescription, and that concern is genuine: the paper itself shows the cold trap vanishes with the Sasaki et al. (2010) decay. I focus instead on Eq. 26 because it is a verifiable internal inconsistency: the normalization from q to q' introduces only one power of q_s, not two. If the printed equations are used literally, the vertical temperature gradient is wrong by a large dimensional factor, so the calculated photosurface and shadow geometry cannot be trusted. This is exactly the kind of missing or incorrect step that the review rules say to flag. I do not accuse the authors of anything; the error may be typographical, or the implemented code may differ from the text. But with no code or data, the manuscript alone does not allow the central claim to be checked. The appropriate disposition remains CONDITIONAL, which is why I leave the reader's verdict unchanged: the paper should be accepted only after the corrected vertical-temperature equation is shown to reproduce the shadowed cold trap, or after code is released that confirms the printed equations were not used literally.","tokens_in":20304,"tokens_out":14226,"duration_ms":166147,"concrete_test":"Independently derive Eq. 24 from the first two moments of the grey radiative transfer equation under the stated Eddington approximation; then correct the normalization to dT^4/dq' = q_s dT^4/dq and recompute the nominal (chi = 0.1, t = 150 kyr) snapshot. If the 9-15 RJ region still shows a ~100 K temperature deficit relative to its surroundings, the cold-trap result survives this internal check; if the feature shifts, weakens, or disappears, the headline claim is unsupported. As a secondary check, evaluate Eq. 33 directly as a line-of-sight shadow mask and compare it with the hatched regions in the figures.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim depends on the vertical temperature structure obtained by integrating Eq. 26. As printed, Eq. 26 is not a valid normalization of Eq. 24: with q' = q/q_s, dq' = dq/q_s, so dT^4/dq' = q_s dT^4/dq. The right-hand side of Eq. 26 instead contains q_s^2, an extra factor of a column density that can change the gradient by orders of magnitude. Because the photosurface altitude (Eq. 15), the shadow condition (Eq. 33), and the midplane temperatures all feed into the claimed 100 K drop at 9-15 RJ, this is not a cosmetic typo. No derivation of Eq. 24 is given and no code is released, so the published equations cannot reproduce Figures 3-6. The accretion-rate sensitivity in Sec. 4.1 is also real, but the normalization error is more fundamental: it puts the model's internal consistency in question before any external parameter choice is considered.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript develops a 2D, quasi-stationary model of Jupiter's circumplanetary disk (CPD) combining an alpha-viscosity surface density profile, Gaussian vertical structure, grey atmosphere radiative transfer, and a geometric self-shadowing prescription. Using analytic fits to the Mordasini (2013) post-runaway accretion and luminosity histories, the authors find that after about 100 kyr Jupiter's radiative heating dominates and the opaque photosurface casts a shadow between roughly 9 and 15 Jupiter radii, producing a ~100 K colder annulus that persists for tens of kyr. They argue that this shadowed zone acts as a cold trap for NH3, CO2, and H2S, with consequences for the composition of the Galilean moons' building blocks. The paper includes sensitivity tests over metallicity, alpha viscosity, and the accretion-rate decay timescale.","tokens_in":20441,"tokens_out":7416,"duration_ms":79736,"significance":"If the result holds, the paper identifies a new and physically plausible mechanism—self-shadowing of a rapidly depleting circumplanetary disk—that can create a persistent cold annulus at ~10 RJ during the epoch of Galilean moon formation. The claimed cold trap is a derived output of the model, not imposed by construction, and the authors test a reasonable range of metallicity and alpha values. They also honestly report that the cold trap disappears when a slower, 1 Myr accretion decay is used, which is an important external limitation but not an internal circularity. A particular strength is that the model makes falsifiable predictions (volatile ice condensation closer than canonical icelines, C/O and N/O enhancements, possible inner/outer moon compositional differences) that can be confronted with JUICE and Europa Clipper data. Reproducibility is currently limited, however, because the key vertical-gradient equation is dimensionally inconsistent as printed and no derivation or code is provided.","major_comments":[{"comment":"I checked the normalization step leading to Eq. (26): with q' = q/q_s, one has dT^4/dq' = q_s dT^4/dq, so the right-hand side of Eq. (26) should indeed contain q_s^2; the earlier reviewer concern about an 'extra factor' in Eq. (26) is therefore not valid. The real problem is Eq. (24) itself: dT^4/dq has units K^4 cm^2 g^-1, whereas the right-hand side ν(q) Σ_g^2 Ω_K^2 κ_R(q) q has units g^2 cm^-2 s^-3, so a Stefan-Boltzmann factor and possibly additional terms from the moment equations are missing. No derivation of Eq. (24) from the Eddington moments is supplied, and no code is released. Since the midplane temperature, the photosurface altitude, and the shadow location all follow from Eq. (24), the published equations as written cannot reproduce Figures 3–6. Please provide the full derivation, correct the equation, and make the numerical implementation available or tabulate enough vertical temperature profiles to allow independent reproduction.","section":"§2.3, Eq. (24)"},{"comment":"The shadow condition is not typeset correctly: 'zs(r) > RJ + zs(r′− RJ) r′ r, ∀r′ < r' is missing parentheses and an operator, and as printed it is not a well-formed inequality. The intended condition appears to be a comparison involving RJ + [zs(r') - RJ] r'/r, but this needs to be written in standard notation. Please also state explicitly that zs = 0 is used when the disk is optically thin, since that convention is invoked in the sentence following Eq. (33).","section":"§2.5, Eq. (33)"},{"comment":"The abstract and conclusion present the ~100 K cold trap as the main finding, but Section 4.1 shows that with the slower Sasaki et al. (2010) accretion decay the shadowed regions exhibit no significant temperature drop. This dependence is acknowledged in the text, which is commendable, but the central claim is therefore conditional on the rapidly decaying Mordasini (2013) prescription. I recommend that the abstract and Summary/Conclusion explicitly state 'for a CPD with a rapidly decaying accretion rate' whenever the cold trap is described, so that readers do not overgeneralize the result to all CPD evolutionary scenarios.","section":"§4.1, Fig. 8"}],"minor_comments":[{"comment":"The word 'axisymetric' should be 'axisymmetric'.","section":"§2.1"},{"comment":"Equation (13) uses κ_R = χ κ_0 T^β, but Table 1 quotes κ_0 in cm^2 g^-1 without specifying that its units depend on β; please state the units as cm^2 g^-1 K^-β or provide κ_0 values in a consistent unit system.","section":"§2.2, Table 1"},{"comment":"The exponential fits to the Mordasini (2013) accretion and luminosity histories are quoted without any comparison to the underlying model data or fit residuals; please provide a small figure or state the fit range and typical accuracy so readers can judge the quality of the fits.","section":"§2.6, Eqs. (34)–(35)"},{"comment":"The word 'metalicity' appears twice and should be 'metallicity'.","section":"§3.2"},{"comment":"The sentence 'which building blocks have to form before 1 Myr of CPD evolution, otherwise, the shadows’ influence are weaker' is grammatically unclear and should be rephrased.","section":"§4.4"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within the scope of A&A and addresses a timely question in Galilean satellite formation. The novelty is moderate but real: a self-consistent treatment of Jupiter's heating and disk self-shadowing in a quasi-stationary CPD model, with a concrete prediction for a cold trap near 10 RJ. I found no evidence of circularity and no citation anomalies. The main obstacle is internal: Eq. (24) is dimensionally inconsistent as printed and no derivation or code is given, so the central numerical result is not currently reproducible from the manuscript. If the authors can supply a correct derivation, fix the equation, and release or archive the code (or a detailed data table), the paper would be acceptable after a major revision. The accretion-rate sensitivity is honestly disclosed but should be foregrounded in the abstract."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"One thing to know: this paper applies self-shadowing to Jupiter's circumplanetary disk and gets a cold trap for ices at ~10 RJ—but the central result rests on a vertical temperature equation that doesn't check out as printed. The idea is worth taking seriously; the manuscript is not yet reproducible.\n\nWhat's new: the gas-starved CPD models of Canup & Ward, Makalkin & Dorofeeva, Heller & Pudritz don't include self-shadowing. Adding a geometric shadowing prescription to the photosurface is a natural and previously missing ingredient. The paper is also honest about the model's sensitivity: it varies metallicity and alpha, and it explicitly shows that switching to the slower Sasaki et al. accretion decay eliminates the cold trap. That transparency is real credit.\n\nWhere it falls down: Eq. 24, the vertical energy equation, is dimensionally inconsistent as printed and has no derivation. Eq. 26, the normalized version, carries an extra q_s factor: with q' = q/q_s, the right side should have q' in the denominator, not q' q_s^2. That extra column density changes the gradient by orders of magnitude, so the temperatures and shadow geometry in Figures 3–6 are not reproducible from the equations given. Eq. 33, the shadow condition, is ambiguous as written. No code or data are provided. The accretion-rate dependence is disclosed, but it means the 100 K drop is a model-dependent hypothesis under one particular formation scenario, not a robust prediction.\n\nFor whom: anyone working on Galilean moon formation or CPD thermal structure, especially those connecting to JUICE/Clipper. It deserves a serious referee, but the review should require the authors to fix or derive the vertical equations, correct the normalization, and release the code before the cold-trap claim can be evaluated. As it stands, I wouldn't cite the cold-trap result.","headline":"Self-shadowing is a sensible new idea for Jovian CPDs, but an extra q_s factor in the vertical temperature equations puts the 100 K cold trap on shaky ground.","tokens_in":21108,"tokens_out":5025,"would_cite":false,"duration_ms":51923,"reading_group":"maybe","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 circumplanetary disk cast shadows that chilled gas by 100 K, creating a cold trap that could set the Galilean moons' icy compositions.","keywords":["Galilean moons","circumplanetary disk","self-shadowing","cold trap","volatile ices","grey atmosphere radiative transfer","Jupiter formation","satellite formation"],"falsifier":"Observations of the Galilean moons' surface and subsurface ices by upcoming spacecraft: if Europa, Ganymede, and Callisto show no radial enrichment pattern for NH3, CO2, or H2S ices, the cold-trap scenario for moon building blocks is contradicted. Alternatively, a measurement or simulation establishing that Jupiter's post-formation CPD accretion rate decayed on a roughly million-year timescale, rather than the rapid decay, would remove the shadowed cold trap, as the paper itself demonstrates.","tokens_in":19974,"feed_emoji":"🪐","tokens_out":9802,"duration_ms":90417,"temperature":0.7,"pith_summary":"This paper argues that Jupiter's circumplanetary disk cast its own shadow: the opaque photosurface of the hot inner disk blocked the young planet's radiation from a ring of outer disk gas, lowering the temperature there by roughly 100 K. Using a two-dimensional, quasi-stationary disk model with grey atmosphere radiative transfer, the authors find the shadowed annulus around 10 Jupiter radii persists for tens of thousands of years during the disk's early depletion. That cold, higher-pressure zone could act as a cold trap, letting volatiles such as NH3, CO2, and H2S condense closer to Jupiter than their normal ice lines. The result matters because it offers a concrete physical route by which the building blocks of the Galilean moons acquired their icy compositions.","feed_headline":"Jupiter's disk shadow chilled gas by 100 K","feed_subtitle":"The cold ring near 10 Jupiter radii could have trapped NH3, CO2, and H2S ices that built the Galilean moons.","key_machinery":"The load-bearing object is the photosurface of the circumplanetary disk: the altitude at which the optical depth reaches tau = 2/3 and the disk becomes opaque to Jupiter's radiation. The model splits the disk into an inner, optically thick adiabatic zone and an outer, optically thin isothermal zone; the surface temperature at the photosurface is set by a balance of viscous heating, accretion heating, and Jupiter's radiative flux, and the vertical temperature profile is then propagated to the midplane with a grey atmosphere Eddington approximation. Self-shadowing is decided by a geometric condition (Eq. 33): a disk patch is in shadow if any inner patch's photosurface altitude projects an opaque wall over it. The whole structure is integrated quasi-statically in time, following an exponentially decaying accretion rate and luminosity fit to a Jupiter formation model, until the photosurface disappears.","core_discovery":"The central claim is that self-shadowing, not just radiative heating, sets the thermal structure of Jupiter's post-formation circumplanetary disk. When the disk becomes optically thin in its outer parts, the optically thick inner region's photosurface acts as an opaque wall that projects a shadow onto the annulus between roughly 9 and 15 Jupiter radii. In the nominal model, this shadowed zone is about 100 K colder than the surrounding gas for the interval from about 100 to 160 kyr after the disk formed, while pressure and density are locally enhanced. The authors show that at 150 kyr this region pulls the condensation fronts of H2O, NH3, CO2, and H2S inward, so these species can form ices at 9-12 RJ rather than at their normal, more distant ice lines. They identify the effect as a cold trap that could seed the Galilean moons with volatile-rich building blocks, with the inner moons potentially enriched relative to the outer ones.","pith_inferences":["A testable extension follows for the Galilean moons: if the cold trap operated, Europa and Ganymede's surface and subsurface ices should show a measurable enrichment of NH3, CO2, and H2S relative to Callisto, a pattern upcoming spacecraft observations could look for.","The same self-shadowing mechanism should operate in the circumplanetary disks of other young giant planets; their disks may show the same roughly 100 K annular cold spots, observable as brightness or molecular-abundance dips at a few planetary radii.","The cold, high-pressure annulus is a natural site for pebble concentration; coupling the thermal model to a dust evolution code would predict whether the shadowed ring seeds a preferred moon-formation radius near 10 Jupiter radii, linking the model to the present-day Galilean satellite architecture.","Because the cold trap's existence hinges on the accretion-rate decay, the model effectively constrains Jupiter's post-formation gas depletion history: if future observations of the moon system's volatile distribution require the cold trap, then the circumplanetary disk must have drained on the fast depletion timescale rather than the slower one."],"forward_implications":["During roughly 127-160 kyr after the circumplanetary disk formed, the annulus near 10 Jupiter radii was a cold trap that let NH3, CO2, and H2S condense closer to Jupiter than their normal ice lines.","The inner Galilean moons, forming in or near this shadowed region, could have accreted a larger share of volatile ices than the outer moons, opposite to a simple radial temperature gradient expectation.","With higher CPD metallicity the shadow lasts longer, up to about 120 kyr for a tenfold enrichment, so the cold trap's influence scales with how dusty the infalling gas was.","If the disk's gas depleted on a slower, roughly million-year timescale, the shadowed zone no longer produces a significant temperature drop, so the cold-trap scenario is specific to a rapidly draining disk.","Local rises in pressure and density inside the shadow could act as dust traps, concentrating solids where proto-moons could form via streaming instability."],"supporting_citations":[{"why":"Provides the gas-starved viscous accretion disk formulation used for the surface density and vertical structure.","marker":"Canup & Ward 2002"},{"why":"Establishes the gas-starved scenario and the centrifugal-radius constraint that the disk extend beyond Callisto's orbit.","marker":"Canup & Ward 2006"},{"why":"Supplies the grey-atmosphere Eddington temperature structure and the ks = 0.2 fraction of Jupiter's light absorbed at the photosurface.","marker":"Makalkin & Dorofeeva 1995"},{"why":"Gives the semi-analytical surface-temperature equation that balances viscous, accretion, and Jupiter radiative heating.","marker":"Makalkin & Dorofeeva 2014"},{"why":"Source of the opacity data used for the Rosseland mean opacity in the radiative transfer.","marker":"Pollack et al. 1994"},{"why":"Fits the Pollack opacity data into the piecewise power law used in the model.","marker":"Makalkin & Dorofeeva 2006"},{"why":"Supplies the rapidly decaying post-formation accretion rate and luminosity curves whose time evolution produces the shadowed cold trap.","marker":"Mordasini 2013"},{"why":"Provides the slower, million-year CPD depletion prescription that, when substituted, removes the cold traps and defines the result's dependence on accretion history.","marker":"Sasaki et al. 2010"}],"fun_headline_variants":["Jupiter's disk self-shadow traps volatiles in cold zone","Self-shadowing cools Jovian disk by 100 K, trapping ices","Shadowed annulus in Jupiter's disk may seed Galilean moons","Cold trap from disk shadow shapes moon-building ices"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The cold trap appears only if Jupiter's circumplanetary disk drained its gas on the rapid timescale the paper adopts from a Jupiter formation model; with the slower, million-year depletion prescription the paper also tests, the shadowed zone stops being cold enough to trap volatiles.","fun_headline_variants_meta":{"raw":{"variants":["Jupiter's disk self-shadow traps volatiles in cold zone","Self-shadowing cools Jovian disk by 100 K, trapping ices","Shadowed annulus in Jupiter's disk may seed Galilean moons","Cold trap from disk shadow shapes moon-building ices"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000286,"raw_usage":{"total_tokens":1722,"prompt_tokens":1023,"completion_tokens":699,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":639,"completion_tokens_details":{"reasoning_tokens":625}},"tokens_in":639,"tokens_out":699,"duration_ms":7180,"temperature":1.0,"reasoning_tokens":625,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T16:31:03.683490+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Observations of the Galilean moons' surface and subsurface ices by upcoming spacecraft: if Europa, Ganymede, and Callisto show no radial enrichment pattern for NH3, CO2, or H2S ices, the cold-trap scenario for moon building blocks is contradicted. Alternatively, a measurement or simulation establishing that Jupiter's post-formation CPD accretion rate decayed on a roughly million-year timescale, rather than the rapid decay, would remove the shadowed cold trap, as the paper itself demonstrates.","supporting_citations":[{"cited_title":"B., & Dorofeeva, V","cited_arxiv_id":null,"evidence_quote":"Supplies the grey-atmosphere Eddington temperature structure and the ks = 0.2 fraction of Jupiter's light absorbed at the photosurface."}],"review_version":1}