{"id":"accca08c-7e3c-47a9-a7ee-56e168fae25d","arxiv_id":"2411.17408","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Post-processing of 3D hydrodynamic simulations with astrochemistry shows that circumplanetary gas is heated to about 800 K, releasing all volatile ices and producing sulfur-bearing molecules that may trace embedded protoplanets.","lead":"This paper simulates the chemistry of gas flowing from a protoplanetary disk into the gravitational grip of a Jupiter-sized planet and finds that intense heating releases frozen water, carbon dioxide, and hydrogen sulfide, producing sulfur molecules like SO, SO2, and CS. A generalist might care because those sulfur molecules could become bright observational beacons that reveal young planets still forming inside their birth disks.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Acknowledged omission of refractory carbon may reverse the C/O change that drives the SO/SO2 tracer prediction; at the claimed ~800 K, refractory carbon release should be included.","rationale":"Good-faith reading: the paper is a transparent post-processing study; the central claim is a qualitative chemical differentiation of CPD gas, with quantitative column densities offered as observational checks. The weakest point is not the assumed gas-dust coupling, which is physically well motivated for small grains at the high volume densities reached in the CPD, but the acknowledged omission of refractory carbon. The model uses fixed elemental abundances and computes sulfur chemistry at T up to 800 K. At those temperatures, refractory carbon, if present, can be released; this changes C/O, and SO/SO2 are strongly C/O-dependent. The authors mention the possibility in Section 4.1 but do not test it. Because the paper's headline tracers are oxidized sulfur molecules, a C/O increase could suppress exactly the species used for the observational comparison. This is an internal limitation, not a disagreement with consensus; it is flagged by the authors themselves. A single, inexpensive network/abundance test would settle it. The reader's verdict of CONDITIONAL remains appropriate; this concern does not overturn the paper but reinforces the conditionality.","tokens_in":24256,"tokens_out":13932,"duration_ms":146686,"concrete_test":"Re-run the disk-streamline ALCHEMIC calculation with an additional refractory carbon abundance C_ref/H between 5e-5 and 1e-4, released when T > 500 K following Bergin et al. (2015), and recompute the SO, SO2, CS, and H2CS column densities of Section 4.4. If SO/SO2 columns fall by more than a factor of 2 or the CS/SO ratio moves outside the observationally inferred range, the tracer prediction is not robust; if the changes are below 50%, the omission is not load-bearing.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that ice sublimation at T up to ~800 K makes the circumplanetary gas oxygen-rich and generates SO/SO2/CS/H2CS. Section 4.1 explicitly states that refractory carbon is ignored 'for simplicity,' even though the same section notes refractory carbon sublimes near ~500 K (Bergin et al. 2015). Since the disk streamline reaches ~800 K (Fig. 6), a refractory carbon reservoir would be at least partially released if it is carried into the CPD. Adding C/H of order 1e-4 (comparable to O/H) would raise C/O above unity, and sulfur chemistry is known to be sensitive to C/O: high C/O suppresses SO and SO2 and favors CS (Semenov et al. 2018). The modeled SO/SO2 columns in Section 4.4 could therefore be an artifact of the assumed elemental abundances rather than a robust prediction of circumplanetary chemistry. The authors flag this possibility but do not quantify it, so the advertised observational consistency is not established. The reader's thermal-coupling concern is less decisive because small grains in the dense CPD are plausibly collisionally coupled to the gas; the refractory-carbon omission is acknowledged, sits at exactly the temperatures central to the claim, and directly affects the proposed tracers.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper presents a post-processing chemical evolution study of gas flowing from a protoplanetary disk into a circumplanetary disk (CPD) around an embedded Jupiter-mass planet. Using 3D hydrodynamics simulations from Lega et al. (2024, Paper I) and the ALCHEMIC chemical network, the authors follow selected streamlines and compute time-dependent chemistry as gas enters the planet's Hill sphere. The main result is that gas heats to ~800 K near the planet, sublimating all volatile ices (H2O, CO2, H2S) and driving high-temperature chemistry that produces sulfur-bearing molecules such as CS, SO, SO2, and H2CS, which are predicted to be enhanced relative to the surrounding protoplanetary disk. The authors estimate column densities of these species over one Hill radius (SO: 2.0e14, SO2: 7.3e14, CS: 4.6e14, H2CS: 1.7e17 cm-2) and argue they are consistent with observations of planet-hosting disks. The paper also discusses implications for volatile delivery to giant planets, the absence of ice lines in warm CPDs, and the timing of moon formation.","tokens_in":24439,"tokens_out":7706,"duration_ms":66115,"significance":"If the predicted tracer species and their column densities are robust, the work provides a concrete path toward detecting and characterizing young embedded protoplanets through molecular line observations. The study combines a high-resolution 3D simulation with a comprehensive chemical network (6065 reactions, 655 species), uses literature-based initial conditions without fitting to the observational targets, and makes falsifiable predictions (column densities, azimuthal asymmetries in SO and SO2 emission). The authors also clearly enumerate caveats, including neglected diffusion, uncertain cooling timescales, and missing accretion UV flux. The main scientific value is in linking gas dynamics to chemistry around young planets and in proposing specific observational tracers. However, the robustness of the central SO/SO2 prediction is not fully established because of the acknowledged omission of refractory carbon, which can alter C/O and thus sulfur chemistry.","major_comments":[{"comment":"The central prediction of enhanced SO and SO2 columns rests on the assumption that the gas in the CPD becomes oxygen-rich through the sublimation of H2O, CO2, and H2S ices. However, Section 4.1 notes that the gas reaches ~800 K, exceeding the sublimation temperature of refractory carbon (~500 K, Bergin et al. 2015), and states that refractory carbon is ignored 'for simplicity.' If refractory carbon is present on the co-moving small grains, it would be at least partially released at these temperatures, potentially raising C/O above unity. Since SO and SO2 formation is suppressed at high C/O (Semenov et al. 2018), the computed column densities in Section 4.4 could be significantly overestimated. The authors flag this possibility but do not quantify it; a sensitivity test including refractory carbon is needed to support the claim that these species are robust tracers of circumplanetary material.","section":"Section 4.1 and Fig. 6"},{"comment":"The column densities of SO, SO2, CS, and H2CS are computed from a single representative streamline that remains in the CPD, with abundances azimuthally averaged over concentric circles. The authors note that the streamline does not visit all regions (hatched areas in Figs. 13 and 14), yet they do not assess how representative this streamline is of the full CPD, despite having computed 10,000 streamlines for the residence-time and flux analysis (Section 3.1). The residence-time map in Fig. 4 shows considerable structure, so a single trajectory may not capture the global temperature and density distribution that sets the chemical abundances. The resulting column densities are therefore uncertain, which weakens the quantitative comparison with observed values.","section":"Sections 2.3 and 4.4"},{"comment":"The authors estimate a gas cooling timescale of 0.1-0.2 orbits, comparable to the orbital timescale of the streamlines (~0.16 orbits), implying that the rapid temperature oscillations along the orbit may not be physical. The authors argue that the highest temperatures and associated chemistry are unaffected, but this reasoning assumes that the peak temperature and its duration are correctly captured. If the gas cannot cool as fast as assumed, the temperature history—and thus the time spent above sublimation and reaction thresholds—would differ. This should be tested by recomputing the chemistry with a thermal model that accounts for the finite cooling timescale.","section":"Section 4.5.2"}],"minor_comments":[{"comment":"The word 'sync' appears to be a typo for 'sink' in the sentence 'could act as a sync for the freed OH.'","section":"Section 4.1"},{"comment":"The word 'spacial' should be 'spatial' in 'spacial resolution chemical survey.'","section":"Section 4.4"},{"comment":"The word 'aziumthally' should be 'azimuthally' in the figure captions.","section":"Figure captions 13 and 14"},{"comment":"The word 'magntiude' should be 'magnitude' in the sentence about the increase in gas density.","section":"Appendix A.1"},{"comment":"The phrase 'We thus possibly overestimating' should be 'We thus possibly overestimate.'","section":"Section 2.2.4"}],"recommendation":"major_revision","confidential_remarks":"The paper is well written and addresses an important topic. The main concern is the robustness of the SO/SO2 tracer prediction given the acknowledged omission of refractory carbon and the single-streamline column density derivation. The stress-test concern about refractory carbon appears valid: the authors themselves note that the gas reaches ~800 K, above the ~500 K sublimation temperature, and they do not quantify the effect on C/O or the resulting sulfur chemistry. I recommend major revision requiring a sensitivity study with refractory carbon and a more representative column density estimate. The paper fits the scope of A&A."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Know this much: the qualitative claim — gas entering a young giant planet's well heats enough to release all volatile ices and drive sulfur chemistry that makes SO, SO2, CS, and H2CS plausible tracers of circumplanetary material — is solid and new in its specifics. The quantitative column densities, though, rest on an acknowledged but untested carbon-budget assumption, and that is the part I would push on in review.\n\nThe genuine contribution: the authors post-process their own 3D fargOCA simulations with ALCHEMIC, and unlike earlier work they resolve deep enough that gas reaches about 800 K near the planet. That sublimates H2O, CO2, and H2S, and the resulting sulfur chemistry (H2S toward CS, SO, SO2) is quantified for circumplanetary material for the first time. The paper converts this into a concrete observational strategy — look for azimuthal asymmetries in SO/SO2/CS emission — and the reported column densities are forward predictions, not fits. The comparisons with PDS 70, HD 169142, and HD 100546 come after the fact, so circularity is low. Credit is also due for the caveats in Section 4.5: diffusion, radiative heating, cooling timescales, and the single snapshot are all stated honestly.\n\nSoft spots, in decreasing seriousness. First, refractory carbon. Section 4.1 notes that the gas reaches the sublimation temperature of refractory carbon (about 500 K; Bergin et al. 2015) and then drops the subject \"for simplicity.\" Since the streamline actually hits about 800 K, a refractory carbon reservoir would be at least partially released; adding C/H near 1e-4 could push C/O above unity, which suppresses SO and SO2 in favor of CS (Semenov et al. 2018). The authors flag this possibility but run no sensitivity test, so the SO/SO2 tracer prediction is conditional rather than established. I rate this above the reader's main concern about gas-dust thermal coupling: small grains in the dense CPD are plausibly collisionally coupled, so that assumption is defensible.\n\nSecond, the column numbers come from a single representative streamline, azimuthally averaged (Section 2.3). They are order-of-magnitude estimates, and the authors' own cooling-timescale argument (0.1-0.2 orbits to cool against roughly 0.16-orbit periods) undercuts the gas-ice oscillation picture in Figure 8 more than they concede. Minor: the H2CS column (1.7e17 cm-2) sits several hundred times above SO, SO2, and CS without comment — probably real, but it deserves a sentence.\n\nThis paper is for anyone working on embedded-protoplanet detection or volatile delivery to giant planet atmospheres. It deserves a serious referee, with a requested refractory-carbon sensitivity test. I would cite the qualitative result regardless.","headline":"Qualitative sulfur-tracer claim is solid and new; the SO/SO2 column densities hinge on an acknowledged but untested refractory-carbon assumption.","tokens_in":25035,"tokens_out":5700,"would_cite":true,"duration_ms":51703,"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":"As gas falls into the gravitational well of an embedded Jupiter-mass planet it heats to roughly 800 K, vaporizes all volatile ices, and produces sulfur molecules that can reveal young planets.","keywords":["circumplanetary disk","planet formation","astrochemistry","sulfur-bearing molecules","volatile delivery","embedded protoplanets","chemical tracers","streamline chemistry"],"falsifier":"A direct test would measure gas and dust temperatures separately inside the Hill sphere of an embedded Jupiter-mass planet: if the small grains remain below the water sublimation temperature while the gas approaches 800 K, the coupled-grain premise fails and the predicted volatile release and sulfur chemistry would not occur. Equivalently, if deep SO and SO$_2$ mapping of a young embedded planet shows column densities more than an order of magnitude below the predicted $2.0\\times 10^{14}$ and $7.3\\times 10^{14}$ cm$^{-2}$, the tracer claim would be falsified.","tokens_in":45,"feed_emoji":"🪐","tokens_out":8243,"duration_ms":150260,"temperature":0.7,"pith_summary":"This paper argues that gas streaming from a protoplanetary disk into the gravitational well of an embedded Jupiter-mass planet is heated to about $\\sim 800$ K, hot enough to vaporize every volatile ice carried by small, co-moving dust grains. The released water, carbon dioxide, and hydrogen sulfide, combined with the warmth, drive a chemical shift that produces sulfur-bearing molecules such as CS, SO, and SO$_2$ that are nearly absent in the cooler disk. The authors compute column densities for these molecules in the circumplanetary disk and report values consistent with existing observations, which would make them practical tracers for locating young embedded protoplanets. If correct, the result also changes how volatiles are delivered to growing giant planets, since key elements would arrive as gas rather than as ice-coated grains.","feed_headline":"Gas falling onto a young Jupiter hits 800 K and frees all ices","feed_subtitle":"This model turns sulfur molecules SO, SO2, CS into signposts for spotting embedded protoplanets.","key_machinery":"The central mechanism is streamline-resolved chemistry: gas parcels are drawn from a snapshot of a 3D hydrodynamic simulation of a Jupiter-mass planet embedded in a disk at 5.2 AU, with a Hill sphere of 0.355 AU, and a zero-dimensional, time-dependent chemical reaction network is integrated along each trajectory using the local temperature and density. Streamlines are classified into families by how long they reside inside the circumplanetary disk, and the chemistry is initialized with elemental abundances and evolved for 1 million years in the outer disk before the gas is released. The load-bearing physical step is the near-800 K temperature peak reached at closest approach to the planet: because the small dust grains are assumed thermally coupled and co-moving, H$_2$O, CO$_2$, and H$_2$S ices sublimate, providing oxidants such as OH that drive the sulfur chemistry toward SO and SO$_2$. Column densities are estimated by averaging abundances on concentric circles around the planet over ten disk heights and integrating up the vertical axis.","core_discovery":"The paper's central claim is that the circumplanetary material formed by an embedded Jupiter-mass planet is chemically distinct from the surrounding protoplanetary disk because gas entering deep into the planet's gravitational well is heated up to $\\sim 800$ K. That heat releases all volatile content from the ice phase of the small grains assumed to ride with the gas, and the resulting high-temperature chemistry converts the dominant sulfur ice H$_2$S into gas-phase CS, SO, SO$_2$, and H$_2$CS. On a length scale of one Hill radius the model yields column densities of $2.0\\times 10^{14}$ cm$^{-2}$ for SO, $7.3\\times 10^{14}$ for SO$_2$, $4.6\\times 10^{14}$ for CS, and $1.7\\times 10^{17}$ for H$_2$CS, which the authors find consistent with previous observational estimates. The paper further claims these species are promising observational tracers of young embedded protoplanets because their column densities exceed those of the surrounding disk, and that the warm circumplanetary disk has no canonical ice lines along its midplane, delaying moon formation until the disk cools.","pith_inferences":["If the tracer claim is right, sulfur-line asymmetries should be a generic feature of actively accreting embedded planets in a few-million-year-old disks, not just the few systems already observed; a systematic survey of planet-hosting disks in SO and SO$_2$ emission would test this.","The weakest step, thermal coupling of gas and small dust, is testable in principle: separate gas and dust temperature maps in the Hill sphere would show whether the ice sublimation trigger actually occurs.","Applying the same streamline-chemistry method to lower-mass planets or wider orbits would predict weaker or absent sulfur tracers, giving a testable trend across planet mass and disk age.","The single representative streamline used for column densities could be replaced by a full 3D chemical map; if the azimuthally averaged columns remain within a factor of a few, the observational conclusions would be more secure."],"forward_implications":["Volatile elements (C, O, N, S) should be delivered to a growing giant planet in the gas phase once the flow enters the circumplanetary disk, rather than as ice mantles on grains.","SO, SO$_2$, CS, and H$_2$CS column densities in the circumplanetary disk are high enough that these molecules should appear as bright azimuthal asymmetries in face-on disks, providing a way to localize young embedded protoplanets.","No water or CO$_2$ ice line exists along the circumplanetary disk midplane in this model, so moon formation is pushed to a later, cooler stage of the system's evolution.","In a dense spherical-envelope accretion mode, SO$_2$ is generated much more abundantly than in the circumplanetary disk case, so sulfur-oxide chemistry may distinguish accretion geometries.","The warm chemistry also dissociates CH$_4$ and produces long-chain hydrocarbons and formic acid, changing the molecular inventory inherited by the planet."],"supporting_citations":[{"why":"It supplies the 3D hydrodynamic simulations and streamline trajectories whose temperatures and densities are post-processed here.","marker":"Lega et al. 2024"},{"why":"It provides the zero-dimensional time-dependent chemical network used to evolve abundances along the streamlines.","marker":"Semenov et al. 2010"},{"why":"It supplies the chemistry solver used to integrate the reaction network along each streamline.","marker":"Semenov 2017"},{"why":"It supplies the initial elemental abundances adopted for the first million years of chemical evolution.","marker":"Bosman et al. 2021"},{"why":"It provides the disk chemical model used for comparing circumplanetary and protoplanetary column densities.","marker":"Cridland et al. 2023"},{"why":"It gives the observed SO column density in an embedded-planet system against which the model's tracer predictions are compared.","marker":"Law et al. 2023"},{"why":"It reports detections of sulfur-bearing molecules in planet-hosting disks that motivate the tracer interpretation.","marker":"Booth et al. 2023a"},{"why":"It provides a weak SO detection in a mature planet-hosting disk that the authors use to bound when sulfur tracers should appear.","marker":"Rampinelli et al. 2024"}],"fun_headline_variants":["Jupiter's accretion heat frees ices, creates SO and CS","Young Jupiter's 800 K inflow yields sulfur fingerprints","Gas falling on forming planet reveals unique chemistry","Ices vaporize near baby Jupiter, forging tracer molecules"],"cache_read_input_tokens":27136,"weakest_assumption_plain":"The load-bearing premise is that the small dust grains are thermally coupled to the gas and move with it, so when the gas reaches about 800 K the ice mantles warm and sublimate; if the grains stay cold, the H$_2$O, CO$_2$, and H$_2$S ices remain frozen and the predicted SO and SO$_2$ chemistry collapses.","fun_headline_variants_meta":{"raw":{"variants":["Jupiter's accretion heat frees ices, creates SO and CS","Young Jupiter's 800 K inflow yields sulfur fingerprints","Gas falling on forming planet reveals unique chemistry","Ices vaporize near baby Jupiter, forging tracer molecules"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000626,"raw_usage":{"total_tokens":2982,"prompt_tokens":1113,"completion_tokens":1869,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":729,"completion_tokens_details":{"reasoning_tokens":1803}},"tokens_in":729,"tokens_out":1869,"duration_ms":21794,"temperature":1.0,"reasoning_tokens":1803,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T12:08:04.341777+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A direct test would measure gas and dust temperatures separately inside the Hill sphere of an embedded Jupiter-mass planet: if the small grains remain below the water sublimation temperature while the gas approaches 800 K, the coupled-grain premise fails and the predicted volatile release and sulfur chemistry would not occur. Equivalently, if deep SO and SO$_2$ mapping of a young embedded planet shows column densities more than an order of magnitude below the predicted $2.0\\times 10^{14}$ and $7.3\\times 10^{14}$ cm$^{-2}$, the tracer claim would be falsified.","supporting_citations":[{"cited_title":"2010, A&A, 522, A42 Semenov,D.A.2017,ALCHEMIC:Advancedtime-dependentchemicalkinetics, Astrophysics Source Code Library, record ascl:1708.008","cited_arxiv_id":null,"evidence_quote":"It provides the zero-dimensional time-dependent chemical network used to evolve abundances along the streamlines."},{"cited_title":"D., Alarcón, F., Bergin, E","cited_arxiv_id":null,"evidence_quote":"It supplies the initial elemental abundances adopted for the first million years of chemical evolution."},{"cited_title":"J., Facchini, S., van Dishoeck, E","cited_arxiv_id":null,"evidence_quote":"It provides the disk chemical model used for comparing circumplanetary and protoplanetary column densities."},{"cited_title":"Gas dynamics around a Jupiter-mass planet I. Influence of protoplanetary disk properties","cited_arxiv_id":"2408.12233","evidence_quote":"It gives the observed SO column density in an embedded-planet system against which the model's tracer predictions are compared."},{"cited_title":"2024, A&A, 689, A65","cited_arxiv_id":null,"evidence_quote":"It provides a weak SO detection in a mature planet-hosting disk that the authors use to bound when sulfur tracers should appear."}],"review_version":1}