Gas dynamics around a Jupiter mass planet: II. Chemical evolution of circumplanetary material

Alex J. Cridland; Elena Lega; Myriam Benisty

REVIEW 3 major objections 5 minor 61 references

Gas dynamics around a Jupiter mass planet: II. Chemical evolution of circumplanetary material

T0 review · 3 major / 5 minor · reviewed 2026-08-12 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read 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.

desk verdict Qualitative sulfur-tracer claim is solid and new; the SO/SO2 column densities hinge on an acknowledged but untested refractory-carbon assumption. read the letter →

arxiv 2411.17408 v1 pith:KVRC4VZL submitted 2024-11-26 astro-ph.EP

Alex J. Cridland , Elena Lega , Myriam Benisty This is my paper · ORCID

classification astro-ph.EP

keywords circumplanetarydiskplanetformationastrochemistrysulfur-bearingmoleculesvolatiledeliveryembeddedprotoplanetschemicaltracersstreamlinechemistry

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

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.

What carries the argument

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.

What would settle it

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.

Watch

Extended reading notes

Core claim

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.

Load-bearing premise

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.

Editorial extensions

If this is right

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.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

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.

Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Desk Editor's Note

read the 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.

The 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.

Soft 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.

Second, 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.

This 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.

Referee Report

3 major / 5 minor

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.

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 (3)

[Section 4.1 and Fig. 6] 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.
[Sections 2.3 and 4.4] 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 4.5.2] 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.

minor comments (5)

[Section 4.1] The word 'sync' appears to be a typo for 'sink' in the sentence 'could act as a sync for the freed OH.'
[Section 4.4] The word 'spacial' should be 'spatial' in 'spacial resolution chemical survey.'
[Figure captions 13 and 14] The word 'aziumthally' should be 'azimuthally' in the figure captions.
[Appendix A.1] The word 'magntiude' should be 'magnitude' in the sentence about the increase in gas density.
[Section 2.2.4] The phrase 'We thus possibly overestimating' should be 'We thus possibly overestimate.'

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the chemical predictions are a forward post-processing calculation with externally sourced physical and chemical inputs, compared to observations only after the fact.

full rationale

The paper does not fit any model parameter to the observations it compares against. The temperature, density, and velocity structure are taken from the companion hydrodynamic simulation of Paper I (Lega et al. 2024), which has overlapping authorship but is a separately published simulation used as a physical input rather than as a justification of the chemical result. The chemical initial abundances are taken from Bosman et al. (2021), the chemical network is the public ALCHEMIC osu_03_2008/KIDA-based network used in Semenov et al. (2018), and the CPD comparison model is the previously published PDS 70 model of Cridland et al. (2023). None of these inputs contains the target result, namely the predicted SO, SO2, CS, and H2CS column densities. The column densities reported in Section 4.4 are computed from the post-processed streamline chemistry and then compared, after the fact, with literature observations such as Law et al. (2023) and Booth et al. (2023a); the paper explicitly labels the values as maximum possible columns because optical depth effects are ignored. The acknowledged omissions — refractory carbon (Section 4.1), diffusion (Section 4.5.1), radiative heating and unresolved shocks (Section 4.5.2), and extra UV from accretion (Section 4.5.3) — are robustness limitations, not circular steps, because none of them is defined in terms of the predicted column densities. The use of a single 'cpd' streamline as a representative of the disk is an estimation procedure stated in Section 2.3, and the residence-time family classification is defined independently from the simulated dynamics and flux weighting; it is not a fit to the observed columns. Self-citations to Paper I and to Cridland et al. (2023) are normal use of prior published models and do not import a forbidden premise such as a uniqueness theorem or an ansatz that already contains the conclusion. Verdict: no significant circularity.

Assumptions & free parameters 4 free parameters · 6 assumptions · 0 invented entities

The model uses a standard astrochemical network and literature abundances; no new physical entities are introduced. The main inputs chosen by hand are the low dust-to-gas ratio for chemically relevant grains, the cosmic ray ionization rate, and the initial reset time, all of which affect the volatile ice budget and hence the sulfur chemistry. The assumption that gas and dust are thermally coupled is the most consequential domain assumption.

free parameters (4)

dust-to-gas mass ratio fdtg (chemically relevant small grains) = 1e-4
Assumed constant in space and time for the population of small, dynamically coupled grains (Section 2.1.2). Sets the grain surface area for freeze-out and surface reactions; directly controls the availability of volatile ices that drive the sulfur chemistry.
cosmic ray ionization rate chi_CR = 1e-17 s^-1
Assumed standard value (Table 2); sets ionization and destruction rates in the chemical network.
initial elemental abundances (C, O, S) = C/H=2.692e-4, O/H=2.880e-4, S/H=1.910e-8
Taken from Bosman et al. (2021), Table 2; sets the overall elemental budget available for chemistry. The authors ignore the CO depletion from that work, possibly overestimating carbon and oxygen.
initial chemical reset time = 1 Myr
The gas is evolved for 1 Myr at its starting position before being released into the streamline (Section 2.2.3), meant to mimic an old disk; the resulting initial molecular abundances depend on this choice.

assumptions (6)

domain assumption The hydrodynamic simulation snapshot from Paper 1 represents a quasi-steady state, and the velocity field is held static during chemical integration.
Section 2.1.1 states a single snapshot is taken to represent the physical structure; if the flow is time variable, the streamlines and temperatures used for chemistry would change.
domain assumption The ALCHEMIC network (osu_03_2008 with KIDA extensions) accurately describes the relevant gas-grain chemistry for C/O/S species.
Relied on for all molecular abundances; the network is standard but has known uncertainties, especially for sulfur chemistry (Section 2.2.2).
domain assumption Gas and dust are thermally coupled and co-moving, with a dust-to-gas ratio of 1e-4 for small grains.
Invoked in Sections 2.1.2 and 3.2.2; needed for volatile ices to desorb when gas warms. If dust remains colder, the sulfur chemistry would be suppressed.
domain assumption The ISM UV field is G0=1 at the disk surface and the host star is a standard T Tauri star with UV and X-ray luminosity 10^31 erg/s.
Used to compute photodissociation and photoionization rates; Section 2.2.1.
domain assumption Diffusion and mixing between the gas parcel and the surrounding gas are negligible.
Explicitly assumed in Section 4.5.1; diffusion could damp the desorption/adsorption oscillations and reduce SO production.
ad hoc to paper Refractory carbon is ignored in the chemistry.
Section 4.1: the authors note refractory carbon could be released at ~500 K and act as a sink for OH, possibly reducing SO and SO2 abundances; they exclude it 'for simplicity'.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Gas dynamics around a Jupiter mass planet: II. Chemical evolution of circumplanetary material." pith.science (2026). https://pith.science/paper/KVRC4VZL

@misc{pith2026241117408,
  author       = {Pith},
  title        = {Pith review of: Gas dynamics around a Jupiter mass planet: II. Chemical evolution of circumplanetary material},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/KVRC4VZL}},
  note         = {Machine review of arXiv:2411.17408}
}

abstract

In an ongoing effort to understand planet formation the link between the chemistry of the protoplanetary disk and the properties of resulting planets have long been a subject of interest. These connections have generally been made between mature planets and young protoplanetary disks through the carbon-to-oxygen (C/O) ratio. In a rare number of systems, young protoplanets have been found within their natal protoplanetary disks. These systems offer a unique opportunity to directly study the delivery of gas from the protoplanetary disk to the planet. In this work we post-process 3D numerical simulations of an embedded Jupiter-massed planet in its protoplanetary disk to explore the chemical evolution of gas as it flows from the disk to the planet. The relevant dust to this chemical evolution is assumed to be small, co-moving grains with a reduced dust-to-gas ratio indicative of the upper atmosphere of a protoplanetary disk. We find that as the gas enters deep into the planet's gravitational well, it warms significantly (up to $\sim 800$ K), releasing all of the volatile content from the ice phase. This change in phase can influence our understanding of the delivery of volatile species to the atmospheres of giant planets. The primary carbon, oxygen, and sulfur carrying ices: CO$_2$, H$_2$O, and H$_2$S are released into the gas phase and along with the warm gas temperatures near the embedded planets lead to the production of unique species like CS, SO, and SO$_2$ compared to the protoplanetary disk. We compute the column densities of SO, SO$_2$, CS, and H$_2$CS in our model and find that their values are consistent with previous observational studies.

Figures

Figures reproduced from arXiv: 2411.17408 by the authors.

**Figure 1.** Example of the output from fargOCA as seen in Paper 1. The figure axes are in units of AU while the colour range shows the surface density of the gas. The black circle shows the extent of the planet’s Hill radius. Model name Σ0/𝑀⊙𝑟 −2 𝑝 𝜈/𝑟 2 𝑝Ω𝑝 Nominal 6.67 × 10−4 10−5 LowMass 6.67 × 10−5 10−5 LowMassLowVis 6.67 × 10−5 10−6 [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗

**Figure 2.** Visualisation of the griding for the column number density calculation, viewing down along the z-axis. The left panel shows [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗

**Figure 3.** Example of some of the streamlines by the embedded planet. Three types of streamlines are shown: one misses the CPD [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗

Figures from the paper (13 more)

**Figure 4.** Figure 4: Map of the residence times for a set of initial conditions begin at a radial distance of 1 R [PITH_FULL_IMAGE:figures/full_fig_p007_4.png]

**Figure 5.** Figure 5: Families: missed (blue), escape (orange), and cpd (red) [PITH_FULL_IMAGE:figures/full_fig_p008_5.png]

**Figure 6.** Figure 6: Temperature evolution for test streamlines. The first col [PITH_FULL_IMAGE:figures/full_fig_p008_6.png]

**Figure 7.** Figure 7: The most abundant initial molecules in their gas and ice phases (prefaced with a ‘g’) for each of the streamline calculations. [PITH_FULL_IMAGE:figures/full_fig_p009_7.png]

**Figure 8.** Figure 8: Evolution of the most abundant oxygen bearing species for the 3 test streamlines. The solid lines represent the gas phase of [PITH_FULL_IMAGE:figures/full_fig_p010_8.png]

**Figure 9.** Figure 9: Same as in figure 8 but for the most abundant carbon bearing species. Again, note the change in y-axis range between the top [PITH_FULL_IMAGE:figures/full_fig_p011_9.png]

**Figure 10.** Figure 10: Same as in figure 8 but for relevant sulfur bearing species. Note here that the axes range on the y-axis has been reduced by [PITH_FULL_IMAGE:figures/full_fig_p012_10.png]

**Figure 11.** Figure 11: figure 11. SO is further enhanced in this high temperature gas [PITH_FULL_IMAGE:figures/full_fig_p012_11.png]

**Figure 11.** Figure 11: The evolution of the gas temperature of the streamline [PITH_FULL_IMAGE:figures/full_fig_p013_11.png]

**Figure 12.** Figure 12: A comparison between a streamline in the nominal model and the LowMass model, using the disk streamline. As before the [PITH_FULL_IMAGE:figures/full_fig_p014_12.png]

**Figure 13.** Figure 13: The aziumthally averaged H2O ice-to-gas mass ratio of the CPD and material surrounding the embedded planet. Recall that the CPD extends to ∼ 0.4 RHill. We base the calculation on the disk streamline and thus the hatched region shows where the streamline did not visit …

**Figure 14.** Figure 14: The aziumthally averaged CO2 ice-to-gas mass ratio of the CPD and material surrounding the embedded planet. The traditional ice line, where the abundance of gas and ice species are equal lies well above the disk midplane due to the warm nature of the circumplanetary m…

**Figure 15.** Figure 15: Left: A comparison between the C2H column density between the different chemical networks mentioned in the text. The networks are Miotello et al. (2019): C/H/O/N and Miotello et al. (2019) + Bruderer (2013): C/H/O/N/S. We tested that the resulting column density showe…

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

61 extracted references · 37 canonical work pages

[1]

M., Lim, P

Astropy Collaboration, Price-Whelan, A. M., Lim, P. L., et al. 2022, ApJ, 935, 167

2022
[2]

Ayliffe, B. A. & Bate, M. R. 2009, MNRAS, 397, 657

2009
[3]

& Morbidelli, A

Batygin, K. & Morbidelli, A. 2020, ApJ, 894, 143 Article number, page 17 of 20 A&A proofs: manuscript no. main

work page 2020
[4]

Bell, K. R. & Lin, D. N. C. 1994, ApJ, 427, 987

1994
[5]

2021, ApJL, 916, L2

Benisty, M., Bae, J., Facchini, S., et al. 2021, ApJL, 916, L2

2021
[6]

A., Blake, G

Bergin, E. A., Blake, G. A., Ciesla, F., Hirschmann, M. M., & Li, J. 2015, Proceedings of the National Academy of Science, 112, 8965

work page 2015
[7]

Bergin, E. A. & Snell, R. L. 2002, ApJL, 581, L105

work page 2002
[8]

P., & Brauer, F

Birnstiel, T., Dullemond, C. P., & Brauer, F. 2010, A&A, 513, A79

2010

Show all 61 references

[9]

S., Ilee, J

Booth, A. S., Ilee, J. D., Walsh, C., et al. 2023a, A&A, 669, A53 Booth,A.S.,Law,C.J.,Temmink,M.,Leemker,M.,&Macias,E.2023b,arXiv e-prints, arXiv:2308.07910

arXiv
[10]

S., Temmink, M., van Dishoeck, E

Booth, A. S., Temmink, M., van Dishoeck, E. F., et al. 2024, AJ, 167, 165

2024
[11]

S., van der Marel, N., Leemker, M., van Dishoeck, E

Booth, A. S., van der Marel, N., Leemker, M., van Dishoeck, E. F., & Ohashi, S. 2021, A&A, 651, L6

2021
[12]

D., Alarcón, F., Bergin, E

Bosman, A. D., Alarcón, F., Bergin, E. A., et al. 2021, ApJS, 257, 7

2021
[13]

D., Walsh, C., & van Dishoeck, E

Bosman, A. D., Walsh, C., & van Dishoeck, E. F. 2018, A&A, 618, A182

2018
[14]

2013, A&A, 559, A46 Bruderer,S.,vanDishoeck,E.F.,Doty,S.D.,&Herczeg,G.J.2012,A&A,541, A91

Bruderer, S. 2013, A&A, 559, A46 Bruderer,S.,vanDishoeck,E.F.,Doty,S.D.,&Herczeg,G.J.2012,A&A,541, A91

2013
[15]

Canup, R. M. & Ward, W. R. 2002, AJ, 124, 3404

2002
[16]

Canup, R. M. & Ward, W. R. 2006,Nature, 441, 834

2006
[17]

Canup, R. M. & Ward, W. R. 2009, in Europa, ed. R. T. Pappalardo, W. B. McKinnon, & K. K. Khurana, 59

2009
[18]

2024, MNRAS, 527, 2049

Chen, K., Kama, M., Pinilla, P., & Keyte, L. 2024, MNRAS, 527, 2049

2024
[19]

& Chiang, E

Choksi, N. & Chiang, E. 2024, arXiv e-prints, arXiv:2403.10057

2024 arXiv
[20]

2019, ApJL, 877, L33 Cleeves,L.I.,Bergin,E.A.,Qi,C.,Adams,F.C.,&Öberg,K.I.2015,ApJ,799, 204

Christiaens, V., Cantalloube, F., Casassus, S., et al. 2019, ApJL, 877, L33 Cleeves,L.I.,Bergin,E.A.,Qi,C.,Adams,F.C.,&Öberg,K.I.2015,ApJ,799, 204

2019
[21]

J., Bosman, A

Cridland, A. J., Bosman, A. D., & van Dishoeck, E. F. 2020, A&A, 635, A68

2020
[22]

J., Facchini, S., van Dishoeck, E

Cridland, A. J., Facchini, S., van Dishoeck, E. F., & Benisty, M. 2023, A&A, 674, A211

2023
[23]

J., van Dishoeck, E

Cridland, A. J., van Dishoeck, E. F., Alessi, M., & Pudritz, R. E. 2019b, A&A, 632, A63 D’Alessio, P., Calvet, N., Hartmann, L., Franco-Hernández, R., & Servín, H. 2006, ApJ, 638, 314 de Val-Borro, M., Edgar, R. G., Artymowicz, P., et al. 2006, MNRAS, 370, 529

2006
[24]

Draine, B. T. 1978, ApJS, 36, 595 Dr˛ ażkowska, J. & Dullemond, C. P. 2018, A&A, 614, A62

1978
[25]

Dullemond, C. P. & Dominik, C. 2005, A&A, 434, 971

2005
[26]

Eistrup, C., Walsh, C., & van Dishoeck, E. F. 2016, A&A, 595, A83

2016
[27]

Eistrup, C., Walsh, C., & van Dishoeck, E. F. 2018, A&A, 613, A14

2018
[28]

2019, ApJ, 887, 152

Fung, J., Zhu, Z., & Chiang, E. 2019, ApJ, 887, 152

2019
[29]

Gibb, E. L. & Horne, D. 2013, ApJL, 776, L28 Güver, T. & Özel, F. 2009, MNRAS, 400, 2050

2013
[30]

Y., Bohn, A

Haffert, S. Y., Bohn, A. J., de Boer, J., et al. 2019, Nature Astronomy, 3, 749

2019
[31]

1998, ApJ, 495, 385

Hartmann, L., Calvet, N., Gullbring, E., & D’Alessio, P. 1998, ApJ, 495, 385

1998
[32]

& Pudritz, R

Heller, R. & Pudritz, R. 2015, The Astrophysical Journal, 806, 181

2015
[33]

Hunter, J. D. 2007, Computing in Science and Engineering, 9, 90

2007
[34]

W., Krĳt, S., & Dong, R

Jiang, H., Wang, Y., Ormel, C. W., Krĳt, S., & Dong, R. 2023, A&A, 678, A33

2023
[35]

Kaufman, M. J. & Neufeld, D. A. 1996, ApJ, 456, 611

1996
[36]

2018, A&A, 617, A44

Keppler, M., Benisty, M., Müller, A., et al. 2018, A&A, 617, A44

2018
[37]

S., et al

Keyte, L., Kama, M., Booth, A. S., et al. 2023, Nature Astronomy, 7, 684

2023
[38]

& Kley, W

Klahr, H. & Kley, W. 2006, A&A, 445, 747

2006
[39]

M., Youdin, A

Krapp, L., Kratter, K. M., Youdin, A. N., et al. 2024, arXiv e-prints, arXiv:2402.14638 Krĳt, S., Bosman, A. D., Zhang, K., et al. 2020, ApJ, 899, 134

2024 arXiv
[40]

& Lega, E

Lambrechts, M. & Lega, E. 2017, A&A, 606, A146 Lambrechts,M.,Lega,E.,Nelson,R.P.,Crida,A.,&Morbidelli,A.2019,A&A, 630, A82

2017
[41]

J., Booth, A

Law, C. J., Booth, A. S., & Öberg, K. I. 2023, ApJL, 952, L19 Lega,E.,Benisty,M.,Cridland,A.,etal.2024,arXive-prints,arXiv:2408.12233

2023 arXiv
[42]

2014, MNRAS, 440, 683

Lega, E., Crida, A., Bitsch, B., & Morbidelli, A. 2014, MNRAS, 440, 683

2014
[43]

2023, ApJ, 952, 89

Marleau, G.-D., Kuiper, R., Béthune, W., & Mordasini, C. 2023, ApJ, 952, 89

2023
[44]

S., Rumpl, W., & Nordsieck, K

Mathis, J. S., Rumpl, W., & Nordsieck, K. H. 1977, ApJ, 217, 425

1977
[45]

K., Rocha, W

McClure, M. K., Rocha, W. R. M., Pontoppidan, K. M., et al. 2023, Nature Astronomy, 7, 431 Meĳerink,R.,Pontoppidan,K.M.,Blake,G.A.,Poelman,D.R.,&Dullemond, C. P. 2009, ApJ, 704, 1471

2023
[46]

F., et al

Miotello, A., Facchini, S., van Dishoeck, E. F., et al. 2019, A&A, 631, A69

2019
[47]

2014, Icarus, 232, 266

Morbidelli, A., Szulágyi, J., Crida, A., et al. 2014, Icarus, 232, 266

2014
[48]

2013, A&A, 558, A113

Mordasini, C. 2013, A&A, 558, A113

2013
[49]

2024, A&A, 689, A65

Rampinelli, L., Facchini, S., Leemker, M., et al. 2024, A&A, 689, A65

2024
[50]

2018, A&A, 617, A28

Semenov, D., Favre, C., Fedele, D., et al. 2018, A&A, 617, A28

2018
[51]

2010, A&A, 522, A42 Semenov,D.A.2017,ALCHEMIC:Advancedtime-dependentchemicalkinetics, Astrophysics Source Code Library, record ascl:1708.008

Semenov, D., Hersant, F., Wakelam, V., et al. 2010, A&A, 522, A42 Semenov,D.A.2017,ALCHEMIC:Advancedtime-dependentchemicalkinetics, Astrophysics Source Code Library, record ascl:1708.008

2010
[52]

& Mori, S

Shibaike, Y. & Mori, S. 2023, MNRAS, 518, 5444

2023
[53]

D., Cugno, G., et al

Stolker, T., Marleau, G. D., Cugno, G., et al. 2020, A&A, 644, A13 Szulágyi, J., Masset, F., Lega, E., et al. 2016, MNRAS, 460, 2853 Szulágyi, J., Morbidelli, A., Crida, A., & Masset, F. 2014, ApJ, 782, 65

2020
[54]

F., et al

Tabone, B., Bettoni, G., van Dishoeck, E. F., et al. 2023, Nature Astronomy, 7, 805 van der Walt, S., Colbert, S. C., & Varoquaux, G. 2011, Computing in Science and Engineering, 13, 22 van Dishoeck, E. F. & Black, J. H. 1982, ApJ, 258, 533 van Dishoeck, E. F., Kristensen, L. E...

2023
[55]

2020, scipy/scipy: SciPy 1.5.3, Zenodo

Virtanen, P., Gommers, R., Burovski, E., et al. 2020, scipy/scipy: SciPy 1.5.3, Zenodo

2020
[56]

2018, A&A, 615, A75

Visser, R., Bruderer, S., Cazzoletti, P., et al. 2018, A&A, 615, A75

2018
[57]

C., et al

Wakelam, V., Herbst, E., Loison, J. C., et al. 2012, ApJS, 199, 21

2012
[58]

2015, A&A, 582, A88

Walsh, C., Nomura, H., & van Dishoeck, E. 2015, A&A, 582, A88

2015
[59]

Wang, Y.-H., Leigh, N. W. C., Perna, R., & Shara, M. M. 2020, ApJ, 905, 136

2020
[60]

G., Ménard, F., Caceres, C., et al

Wolff, S. G., Ménard, F., Caceres, C., et al. 2017, AJ, 154, 26

2017
[61]

Zuckerman, B., Kuiper, T. B. H., & Rodriguez Kuiper, E. N. 1976, ApJL, 209, L137 Article number, page 18 of 20

1976

Pith tools

Reviewed August 12, 2026 · model on record in the stance chip above.

[1] [1]

M., Lim, P

Astropy Collaboration, Price-Whelan, A. M., Lim, P. L., et al. 2022, ApJ, 935, 167

2022

[2] [2]

Ayliffe, B. A. & Bate, M. R. 2009, MNRAS, 397, 657

2009

[3] [3]

& Morbidelli, A

Batygin, K. & Morbidelli, A. 2020, ApJ, 894, 143 Article number, page 17 of 20 A&A proofs: manuscript no. main

work page 2020

[4] [4]

Bell, K. R. & Lin, D. N. C. 1994, ApJ, 427, 987

1994

[5] [5]

2021, ApJL, 916, L2

Benisty, M., Bae, J., Facchini, S., et al. 2021, ApJL, 916, L2

2021

[6] [6]

A., Blake, G

Bergin, E. A., Blake, G. A., Ciesla, F., Hirschmann, M. M., & Li, J. 2015, Proceedings of the National Academy of Science, 112, 8965

work page 2015

[7] [7]

Bergin, E. A. & Snell, R. L. 2002, ApJL, 581, L105

work page 2002

[8] [8]

P., & Brauer, F

Birnstiel, T., Dullemond, C. P., & Brauer, F. 2010, A&A, 513, A79

2010