Solving for the 2D Water Snowline with Hydrodynamic Simulations. Emergence of gas outflow, water cycle and temperature plateau

Chris W. Ormel; Shoji Mori; Xue-Ning Bai; Yu Wang

REVIEW 3 major objections 3 minor 115 references

Solving for the 2D Water Snowline with Hydrodynamic Simulations. Emergence of gas outflow, water cycle and temperature plateau

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

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

Pith's one-line read 2D simulations show that vapor injection at the snowline drives an outflow that piles up ice, and that active disks add a water cycle trapping still more.

desk verdict A genuine 2D advance on the water snowline, with the observable 40 K dip as the least secure part. read the letter →

arxiv 2502.08936 v1 pith:NKGB26HQ submitted 2025-02-13 astro-ph.EP

Yu Wang , Chris W. Ormel , Shoji Mori , Xue-Ning Bai This is my paper · ORCID

classification astro-ph.EP

keywords watersnowlineprotoplanetarydiskspebblepile-uplatentheattwo-streamradiationtransfermultifluidhydrodynamicscycleplanetesimalformation

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 the water snowline in a protoplanetary disk cannot be understood with 1D vertically averaged isothermal models: in 2D, vapor released by sublimating pebbles drives a radial gas outflow that piles ice up just outside the snowline, and in actively heated disks a vertical water cycle, sublimation at the hot midplane, recondensation in cooler upper layers, and settling back down, traps substantially more ice. It further shows that the latent heat of sublimation flattens the temperature profile across the snowline, broadening the pile-up, lowering its peak solid-to-gas ratio, and producing a temperature dip of about 40 K that should appear as a continuum intensity dip at millimeter wavelengths in outbursting disks. If correct, this changes where and how strongly planetesimal formation is expected at snowlines, favors young active disks as planetesimal factories, and offers a new observational channel for locating water snowlines. The paper establishes these results with 2D multifluid hydrodynamic simulations that self-consistently couple pebble dynamics, vapor transport, phase change, and a two-stream radiation transfer temperature solution.

What carries the argument

The carrying mechanism is the coupling between the phase-change module and the 2D gas flow: sublimating pebbles inject vapor that is viscously spread both ways, steepening the midplane density gradient and creating a radial outflow, while in vertically stratified active disks the closed loop of midplane sublimation, upward vapor diffusion, upper-layer recondensation, and settling of grown pebbles constitutes the water cycle. The temperature is set by a two-stream radiation transfer solver whose heating terms include stellar irradiation, midplane viscous dissipation, heat diffusion, and the latent-heat term $q_{\mathrm{latent}} = -d\rho_{\mathrm{vap}}/dt \cdot L_{\mathrm{latent}}$, so the phase change feeds back on the thermal structure that defines the snowline.

What would settle it

Rerun the active-disk simulation with Rosseland opacity tied to ice mass and with viscous or Joule heating deposited two to three scale heights above the midplane; if the midplane no longer stays warmer than the upper layers, the water cycle, the factor of about 1.5 ice-mass enhancement, and the about 40 K plateau should disappear, and observationally this can be tested by mapping continuum brightness temperature and spectral index across the snowline of an outbursting disk such as V883 Ori.

Watch

Extended reading notes

Core claim

At steady state, icy pebbles drifting inward release vapor at the snowline; the vapor's viscous spreading steepens the radial density gradient and drives a midplane gas outflow roughly ten times stronger than the background accretion flow, producing an advective outward vapor flux comparable to the diffusive one and an ice pile-up exterior to the snowline. In active disks, where the midplane is hotter than the upper layers, the paper identifies a water cycle: ice sublimes at the midplane, vapor diffuses upward and recondenses on pebbles in cooler layers, and the grown pebbles settle back to the midplane, suppressing the upper-layer vapor escape channel and increasing trapped ice mass by about a factor of 1.5 relative to passive or isothermal disks. Latent heat exchange flattens the temperature gradient across the sublimation front, widening the pile-up and reducing its peak solid-to-gas ratio, with the strongest effect in optically thick active disks; the associated roughly 40 K cooling plateau is proposed as the origin of intensity dips in the dust continuum of outbursting disks such as V883 Ori.

Load-bearing premise

The thermal structure, and with it the water cycle and the roughly 40 K latent-heat plateau, assumes constant Rosseland and visible opacities and that all viscous heating is deposited at the midplane; if heating is deposited a few scale heights above the midplane or ice grains dominate the opacity, the effects weaken or shift, as the paper concedes in Section 6.4.

Editorial extensions

If this is right

The snowline solid pile-up in real disks is stronger and broader than 1D isothermal models predict: active disks store about 1.5 times more ice mass, and the resolved outflow boosts the peak solid-to-gas ratio by about 1.7 over the fixed-gas 1D run.
Active disks with hot midplanes are more conducive to planetesimal formation than passive or isothermal disks, because the water cycle traps water at the snowline.
Latent heat cooling produces a roughly 20 percent (40 K) temperature plateau extending about 50 percent of the snowline radius, which should appear as a dark annulus in millimeter continuum independent of dust opacity changes.
The reduced headwind near the snowline, by a factor of about 0.6 to 0.7, lowers the thresholds for streaming instability and pebble accretion, so planetesimals born there can start pebble accretion immediately and grow icy cores on timescales of about 0.05 Myr.
The vapor-injection outflow slows pebble drift and can even reverse it at high viscosity; at $\alpha = 0.03$ this produces a silicate pile-up with peak solid-to-gas ratio near 1.5, making the snowline a more general solid barrier than pure vapor-retrodiffusion predicts.

Reading between the lines

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

If the latent-heat dip is real, then dark annuli in FUor-class outbursting disks could be used to map the snowline radius and accretion luminosity without relying on dust-opacity modeling, a step the paper suggests but does not fully develop.
The water cycle implies that pile-up strength and composition depend on vertical thermal structure, so observable proxies such as dust scale height and accretion rate in young Class 0 and I disks could indicate where planetesimal formation is most favored.
A testable extension is to check whether the latent-heat-broadened pile-up produces a wider region of enhanced ice content than 1D models predict, using spatially resolved water isotopologue line observations across the snowline in outbursting systems.
Because the paper assumes constant opacities and midplane-deposited viscous heating, the water-cycle enhancement is likely an upper bound; a direct simulation with ice-dominated opacity and heating deposited several scale heights above the midplane would show whether the enhancement weakens to the passive or isothermal level.

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

The paper is worth reading. It does something genuinely new: it solves the water snowline in 2D (R-Z) with coupled pebble dynamics, vapor transport, and a self-consistent thermal structure including latent heat. The main mechanisms it identifies — a vapor-driven radial outflow near the snowline, a 'water cycle' in actively heated disks that traps ice, and a latent-heat temperature plateau — are not present in the 1D vertically averaged models they cite. The numerical work is careful: they verify steady state via flux conservation, compare against 1D runs, and decompose the radial velocity to show the outflow is viscous spreading of vapor. That is solid evidence for the qualitative picture.

The soft spots are the quantitative edges. The most striking prediction — a ~40 K temperature dip in active disks manifesting as an ALMA intensity dip, possibly explaining V883 Ori — is asserted, not computed. There is no synthetic emission calculation, so the connection from the two-stream temperature to an observed continuum depression is a plausible argument, not a demonstrated result. The stress-test on this point is fair and the authors themselves flag the reliance on midplane heating and constant ice-free opacity in Sect. 6.4. If MHD/Joule heating deposits energy higher up, or icy mantles change the opacity, the blanketing regime that makes latent heat important weakens. So the qualitative mechanisms are likely robust; the specific 40 K dip and its observability are the least secure part.

I also note the lack of code/data release, which lowers reproducibility, and the narrow parameter exploration (though they do show the water cycle survives a range of alpha). These are minor relative to the main claims. The paper is honest about its limitations and does not over-fit its conclusions to observations.

Who is this for? Anyone working on snowlines, pebble growth, and planetesimal formation. It deserves a serious referee. The right revision path is: keep the hydrodynamics and water cycle as the core contribution, and either add a radiative transfer post-processing of the intensity or reframe the observational section as a testable prediction rather than an explanation of V883 Ori.

My recommendation: send to peer review. It is a genuine advance, and the flaws are fixable without overturning the main results.

Referee Report

3 major / 3 minor

Summary. This paper presents 2D (R-Z) multifluid hydrodynamic simulations of the water snowline in protoplanetary disks, using Athena++ with a phase-change module and a two-stream radiation transfer solver. The authors report three main findings: (i) vapor injection at the snowline drives a radial gas outflow that augments the outward vapor flux; (ii) in actively heated disks, a 'water cycle'—sublimation in the hotter midplane, recondensation in cooler upper layers, and settling—traps ice and enhances the pile-up; and (iii) latent heat exchange produces a ~40 K temperature plateau/dip that they argue should appear as an intensity dip in dust continuum observations of outbursting disks such as V883 Ori. They conclude that 2D gas dynamics and realistic thermal structure qualitatively change snowline pile-up predictions and favor planetesimal formation in active disks.

Significance. If the main claims hold, this is a significant contribution to planet formation theory. It extends previous 1D, vertically-averaged snowline models to include the 2D flow pattern and vertical temperature structure, and it identifies a vapor-driven outflow and a water cycle that plausibly enhance solid pile-up at the snowline. The paper is careful in its numerics: it verifies steady state with flux conservation, compares 1D and 2D runs, and provides a transparent table of parameters and outputs. The latent-heat temperature dip is an intriguing, potentially observable signature. However, the observational claim is not yet backed by a synthetic emission calculation, and the active-disk results are conditional on the assumed midplane heating and constant ice-free opacity.

major comments (3)

[Sect. 6.3 and Conclusions item 6] The claim that the ~40 K latent-heat temperature dip "manifests as an intensity dip in the dust continuum" and can explain the V883 Ori depression is asserted rather than demonstrated; no radiative transfer post-processing of the emergent intensity is presented, and the argument implicitly assumes that the continuum is optically thick at the sublimation front, that the emission layer lies at or below the vertically extended sublimation surface, and that dust opacity is unchanged across the snowline. The authors should either compute synthetic emission maps (e.g., with a ray-tracing or Monte Carlo radiative transfer code) or, at minimum, explicitly reframe this as a qualitative suggestion; as written, the abstract and conclusions present it as a result.
[Sect. 2.5 (Eq. 29) and Sect. 6.4] The active-disk results—the water cycle, the enhanced pile-up, and the 40 K plateau—depend on the assumptions that viscous heating is deposited at the midplane (Eq. 25) and that the Rosseland opacity is constant and ice-free. The authors themselves concede in Sect. 6.4 that MHD/Joule heating at several scale heights above the midplane or ice-dominated opacity would weaken or shift these effects. Because the conclusion that "active disks are more conducive to planetesimal formation" rests on these assumptions, the manuscript should state this conditionality more prominently and ideally test the sensitivity with a variant in which the heating is deposited at higher altitude.
[Sect. 2.2 and Sect. 6.1] The backreaction of solids on gas is neglected, even though peak midplane solid-to-gas ratios reach ~0.5–0.8 in the 2D runs (Table 1). Since the paper's arguments about streaming-instability thresholds and the strength of the pile-up depend on ξpk, and the authors cite Schoonenberg & Ormel (2017) showing that backreaction can double the enhancement, the reported ξpk values should be treated as lower bounds; the paper should discuss whether the outflow and water-cycle conclusions are robust to including dust backreaction, or clearly state that they are not.

minor comments (3)

[Throughout] Throughout the manuscript, ligature spacing errors produce "e ffect", "di ffusion", "di fferent" in many places (e.g., Abstract, Sect. 1, Sect. 2.4); these should be corrected.
[Fig. 2 caption and Sect. 3.2] The quoted value "ρ0≈ 3.4×10−11 cm g−3" appears to have the units inverted; it should read "g cm−3".
[Page 1 header] The header on page 1 lists "Received September 15, 1996; accepted March 16, 1997", which are clearly incorrect dates; the manuscript should be updated.

Circularity Check

0 steps flagged · score 2.0 of 10

No significant circularity: the snowline pile-up, outflow, water cycle, and latent-heat effects emerge from the solved multifluid equations; the few self-citations are to previously published methods and are not load-bearing as unverified premises.

full rationale

The central derivations are self-contained numerical integrations. Inputs are the gas/pebble/vapor conservation equations (Eqs. 3-8), the phase-change rate law (Eq. 16) with Clausius-Clapeyron saturation (Eq. 17), and the two-stream radiative equilibrium temperature (Eqs. 21-29). The reported outputs, including the vapor-driven outflow, the water cycle, the pile-up morphology, xi_pk, FWHM, Mice, and the ~40 K latent-heat temperature plateau, are read off the steady-state solutions rather than fitted to those conclusions. Anchoring the initial midplane snowline at 2.1 au by choosing kappa_R and L_star is a calibration of the model setup, not a prediction of the absolute snowline radius; the subsequently reported rsnow,mid shifts and all pile-up indicators are emergent. The two most prominent self-citations, Wang et al. (2023) for the phase-change module and Mori et al. (2019) for the two-stream transfer formula, supply previously published, independently usable numerical and analytic methods; neither is an unverified premise whose conclusion is the present paper's result. The one extrapolative step, that the 40 K dip "manifests as an intensity dip" in V883 Ori (Sect. 6.3), is asserted without a synthetic emission calculation, and the authors themselves note in Sect. 6.4 that Joule-heating altitude and ice opacity could weaken the effect. That is a limitation or an unsupported observational inference, not a circular reduction: the temperature dip is computed, not assumed, and no quantity is renamed from its own definition. Therefore no step satisfies the criteria for circularity; score 2 reflects only the presence of minor, non-load-bearing self-citations in the methods.

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

The model relies on a set of standard disk and pebble assumptions plus several hand-chosen parameters. The opacity values and stellar luminosities are selected to anchor the initial snowline at 2.1 au and to produce the three thermal morphologies; these choices are physically motivated but are not fitted to the output claims. The most load-bearing assumptions are the midplane-deposited viscous heating and constant opacities, which directly determine the strength of the water cycle and latent-heat plateau.

free parameters (10)

alpha (viscosity parameter) = 3e-3
Sets gas accretion speed, diffusion, and flow pattern; chosen following Schoonenberg & Ormel (2017). Varying it to 0.03 or 0.001 changes pile-up strength (Appendix C).
Mdot_acc (mass accretion rate) = 1e-8 Msun/yr
Fixes gas surface density via Eq. (2); chosen following Schoonenberg & Ormel (2017).
fi/g (ice-to-gas flux ratio) = 0.4
Controls pebble injection and vapor abundance; assumed constant. Lower/higher values would change pile-up and latent-heat strength.
tau_s,0 (initial Stokes number at r0) = 0.03
Sets initial pebble size (~2 cm); chosen following Schoonenberg & Ormel (2017); affects drift and settling, and thus pile-up and water cycle.
kappa_R (Rosseland mean opacity) = 4.0, 0.01, 1.2e-5 cm2/g for active, passive, iso runs
Chosen by hand to place the initial midplane snowline at 2.1 au and to select the thermal morphology (Fig. 1). Directly sets the vertical temperature structure and the blanketing effect.
L* (stellar luminosity) = 1.0, 22.0, 1.0 Lsun
Chosen with kappa_R to anchor the snowline at 2.1 au; passive disk uses 22 Lsun, active and iso use 1 Lsun. Controls irradiation heating and thus snowline morphology.
beta (thermal relaxation parameter) = 50
Thermal relaxation timescale in Eq. (32); chosen to prevent numerical oscillations. In steady state the solution is independent of beta.
kappa_vi (visible opacity) = 10 cm2/g
Fixed for all runs, approximating a DSHARP-like grain size distribution; affects irradiation heating depth.
Hp (photosphere height) = 4 Hg
Assumed for the irradiation heating geometry in Eq. (24).
b (viscosity power-law index) = -1
Viscosity profile exponent; following Schoonenberg & Ormel (2017).

assumptions (6)

domain assumption Alpha-disk viscosity with viscous heating deposited at the midplane (Shakura & Sunyaev 1973)
Used in Sect. 2.1 and 2.5; the water cycle and latent-heat plateau require a hot midplane. The authors note in Sect. 6.4 that MHD (Joule) heating would deposit energy at several scale heights, weakening these effects.
domain assumption Constant Rosseland and visible opacities, independent of ice abundance
Assumed in Sect. 2.6; ice mantle opacity is excluded, which could shift the snowline and alter the blanketing effect (acknowledged in Sect. 6.4).
domain assumption No back-reaction of pebbles on the gas
Stated in Sect. 2.2 and discussed in Sect. 6.1; including back-reaction could boost the solid-to-gas ratio by a factor of a few and affect the outflow.
domain assumption Single-seed pebble model where the silicate core remains intact during sublimation
Sect. 2.3; if pebbles disaggregate, the pile-up would be stronger (traffic-jam scenario), changing the quantitative outcomes.
domain assumption Two-stream approximation separating stellar irradiation and disk thermal emission
Sect. 2.5; this standard radiative transfer assumption underlies the computed temperature structure and the latent-heat plateau.
standard math Ideal gas equation of state with variable mean molecular weight
Eq. (19)-(20) in Sect. 2.4; standard equation of state for the gas and vapor mixture.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Solving for the 2D Water Snowline with Hydrodynamic Simulations. Emergence of gas outflow, water cycle and temperature plateau." pith.science (2026). https://pith.science/paper/NKGB26HQ

@misc{pith2026250208936,
  author       = {Pith},
  title        = {Pith review of: Solving for the 2D Water Snowline with Hydrodynamic Simulations. Emergence of gas outflow, water cycle and temperature plateau},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/NKGB26HQ}},
  note         = {Machine review of arXiv:2502.08936}
}

read the original abstract

In protoplanetary disks, the water snowline marks the location where ice-rich pebbles sublimate, releasing silicate grains and water vapor. These processes can trigger pile-ups of solids, making the water snowline a promising site for forming planetesimals. However, previous studies exploring the pile-up conditions typically employ 1D, vertically-averaged and isothermal assumptions. In this work, we investigate how a 2D flow pattern and realistic temperature structure affect the pile-up of pebbles at the snowline and how latent heat effects can leave observational imprints. We perform 2D (R-Z) multifluid hydrodynamic simulations, tracking chemically heterogeneous pebbles and the released vapor. With a recent-developed phase change module, the mass transfer and latent heat exchange during ice sublimation are calculated self-consistently. The temperature is calculated by a two-stream radiation transfer method under various opacities and stellar luminosity. We find that vapor injection at the snowline drives a previously unrecognized outflow, leading to a pile-up of ice outside the snowline. Vapor injection also decreases the headwind velocity in the pile-up, promoting planetesimal formation and pebble accretion. In active disks, we identify a water-cycle: after ice sublimates in the hotter midplane, vapor recondenses onto pebbles in the upper, cooler layers, which settle back to the midplane. This cycle promotes ice-trapping at snowline. Latent heat exchange flattens the temperature gradient across the snowline, broadening the width while reducing the peak solid-to-gas ratio of pile-ups. Due to the water cycle, active disks are more conducive to planetesimal formation than passive disks. The significant temperature dip (~ 40K) caused by latent heat cooling manifests as an intensity dip in the dust continuum, presenting a new channel to identify the water snowline in outbursting systems.

Figures

Figures reproduced from arXiv: 2502.08936 by the authors.

**Figure 1.** Initial midplane snowline location resulting from the two-stream radiation transport model (equation (29)) as function of Rosseland mean opacity κR and stellar luminosity L⋆. Here the optical opacity is fixed to κvi = 10 cm2g −1 . The permitted parameter space is divided into three regimes (white dashed lines, for illustration purpose only): “passive disk” when L⋆ is high, “active disk” when κR is high and “vertical… view at source ↗

**Figure 2.** Steady state density structure of all simulated disk (as listed in [PITH_FULL_IMAGE:figures/full_fig_p008_2.png] view at source ↗

**Figure 3.** Steady state temperature structure of all the simulation runs ( [PITH_FULL_IMAGE:figures/full_fig_p009_3.png] view at source ↗

Figures from the paper (6 more)

**Figure 4.** Figure 4: Radial profiles of the passive run. (a): Vertically-integrated surface density of total gas, ice, silicate and vapor. The grey line shows the expected vapor surface density from the unperturbed viscous solution (equation (2)), where fi/g is the ice-to-gas flux ratio. …

**Figure 5.** Figure 5: Upper panels: Vertically-integrated radial mass flux of vapor. For three runs (passive, active and active-nqL), the diffusive, advective and total vapor flux are plotted. The dotted lines denote zero-flux levels. The blue arrows indicate the position of the midplane sn…

**Figure 6.** Figure 6: Azimuthal velocity deviation at the midplane, expressed in terms of the radial pressure gradient parameter η (equation (35)). The black dashed and solid lines show η before vapor injection and after the steady state is reached, respectively. The red lines show the vapo…

**Figure 7.** Figure 7: Comparisons of radial profiles of the active disks. For all three panels, the solid line represents the active disk while the dashed line represents active-nqL disk. (a): The surface density profiles. As in [PITH_FULL_IMAGE:figures/full_fig_p013_7.png]

**Figure 8.** Figure 8: Upper panels: Summary of the Lagrangian trajectory integration in the active-nqL and passive-nqL disks. Yellow-green color shading denotes the total time spent in each grid cell by the water particles, which are released from the outer boundary. The black dashed lines …

**Figure 9.** Figure 9: Summary of the pile-up properties of 1D and 2D runs as listed in [PITH_FULL_IMAGE:figures/full_fig_p015_9.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

115 extracted references · 48 canonical work pages

[1]

, " * write output.state after.block = add.period write newline

ENTRY address archiveprefix author booktitle chapter edition editor howpublished institution eprint journal key month note number organization pages publisher school series title type volume year label extra.label sort.label short.list INTEGERS output.state before.all mid.sentence after.sentence after.block FUNCTION init.state.consts #0 'before.all := #1 ...
[2]

write newline

" write newline "" before.all 'output.state := FUNCTION n.dashify 't := "" t empty not t #1 #1 substring "-" = t #1 #2 substring "--" = not "--" * t #2 global.max substring 't := t #1 #1 substring "-" = "-" * t #2 global.max substring 't := while if t #1 #1 substring * t #2 global.max substring 't := if while FUNCTION word.in bbl.in " " * FUNCTION format....
[3]

P., Simon , J

Abod , C. P., Simon , J. B., Li , R., et al. 2019, , 883, 192

2019
[4]

Armitage , P. J. 2020, Astrophysics of planet formation, Second Edition

2020
[5]

& Wurm , G

Aumatell , G. & Wurm , G. 2011, , 418, L1

2011
[6]

2017, , 845, 75

Bai , X.-N. 2017, , 845, 75

2017
[7]

& Stone , J

Bai , X.-N. & Stone , J. M. 2010 a , , 722, 1437

2010
[8]

& Stone , J

Bai , X.-N. & Stone , J. M. 2010 b , , 722, L220

2010

Show all 115 references

[9]

Balbus , S. A. & Papaloizou , J. C. B. 1999, , 521, 650

1999
[10]

A., Yang , C.-C., & Zhu , Z

Baronett , S. A., Yang , C.-C., & Zhu , Z. 2024, , 529, 275

2024
[11]

2024, , 62, 157

Birnstiel , T. 2024, , 62, 157

2024
[12]

P., & Brauer , F

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

2010
[13]

P., Zhu , Z., et al

Birnstiel , T., Dullemond , C. P., Zhu , Z., et al. 2018, , 869, L45

2018
[14]

2012, , 539, A148

Birnstiel , T., Klahr , H., & Ercolano , B. 2012, , 539, A148

2012
[15]

& Wurm , G

Blum , J. & Wurm , G. 2000, , 143, 138

2000
[16]

& Wurm , G

Blum , J. & Wurm , G. 2008, , 46, 21

2008
[17]

C., & D'Alessio , P

Calvet , N., Patino , A., Magris , G. C., & D'Alessio , P. 1991, , 380, 617

1991
[18]

Carrera , D., Johansen , A., & Davies , M. B. 2015, , 579, A43

2015
[19]

1935, , 96, 21

Chandrasekhar , S. 1935, , 96, 21

1935
[20]

& Cowling , T

Chapman , S. & Cowling , T. G. 1991, The Mathematical Theory of Non-uniform Gases

1991
[21]

Chiang , E. I. & Goldreich , P. 1997, , 490, 368

1997
[22]

Ciesla , F. J. 2009, , 200, 655

2009
[23]

Ciesla , F. J. 2010, , 723, 514

2010
[24]

A., Casassus , S., Tobin , J., et al

Cieza , L. A., Casassus , S., Tobin , J., et al. 2016, , 535, 258

2016
[25]

J., Lambrechts , M., van Kooten , E., & Johansen , A

Colmenares , M. J., Lambrechts , M., van Kooten , E., & Johansen , A. 2024, , 685, A114

2024
[26]

Cuzzi , J. N. & Zahnle , K. J. 2004, , 614, 490

2004
[27]

& Tielens , A

Dominik , C. & Tielens , A. G. G. M. 1997, , 480, 647

1997
[28]

& Alibert , Y

Dr a \.z kowska , J. & Alibert , Y. 2017, , 608, A92

2017
[29]

2023, in Astronomical Society of the Pacific Conference Series, Vol

Dr a \.z kowska , J., Bitsch , B., Lambrechts , M., et al. 2023, in Astronomical Society of the Pacific Conference Series, Vol. 534, Protostars and Planets VII, ed. S. Inutsuka , Y. Aikawa , T. Muto , K. Tomida , & M. Tamura , 717

2023
[30]

J., Hillenbrand , L

Fischer , W. J., Hillenbrand , L. A., Herczeg , G. J., et al. 2023, in Astronomical Society of the Pacific Conference Series, Vol. 534, Protostars and Planets VII, ed. S. Inutsuka , Y. Aikawa , T. Muto , K. Tomida , & M. Tamura , 355

2023
[31]

& Schmitt , B

Fray , N. & Schmitt , B. 2009, , 57, 2053

2009
[32]

Gammie , C. F. 1996, , 457, 355

1996
[33]

M., et al

Gillon , M., Jehin , E., Lederer , S. M., et al. 2016, , 533, 221

2016
[34]

Gillon , M., Triaud , A. H. M. J., Demory , B.-O., et al. 2017, , 542, 456

2017
[35]

M., Carrasco-Gonz \'a lez , C., Mac \' as , E., et al

Guerra-Alvarado , O. M., Carrasco-Gonz \'a lez , C., Mac \' as , E., et al. 2024 a , , 686, A298

2024
[36]

M., van der Marel , N., Di Francesco , J., et al

Guerra-Alvarado , O. M., van der Marel , N., Di Francesco , J., et al. 2024 b , , 681, A82

2024
[37]

M., S \'a ndor , Z., Ronco , M

Guilera , O. M., S \'a ndor , Z., Ronco , M. P., Venturini , J., & Miller Bertolami , M. M. 2020, , 642, A140

2020
[38]

& Blum , J

Gundlach , B. & Blum , J. 2015, , 798, 34

2015
[39]

P., Kreuzig , C., et al

Gundlach , B., Schmidt , K. P., Kreuzig , C., et al. 2018, , 479, 1273

2018
[40]

& Krijt , S

Houge , A. & Krijt , S. 2023, , 521, 5826

2023
[41]

2024, , 527, 9668

Houge , A., Mac \' as , E., & Krijt , S. 2024, , 527, 9668

2024
[42]

& Bai , X.-N

Hu , Z. & Bai , X.-N. 2021, , 503, 162

2021
[43]

M., P \'e rez , L

Huang , J., Andrews , S. M., P \'e rez , L. M., et al. 2018, , 869, L43

2018
[44]

& Bai , X.-N

Huang , P. & Bai , X.-N. 2022, , 262, 11

2022
[45]

1990, , 351, 632

Hubeny , I. 1990, , 351, 632

1990
[46]

Hyodo , R., Guillot , T., Ida , S., Okuzumi , S., & Youdin , A. N. 2021, , 646, A14

2021
[47]

2019, , 629, A90

Hyodo , R., Ida , S., & Charnoz , S. 2019, , 629, A90

2019
[48]

2016, , 591, A72

Ida , S., Guillot , T., & Morbidelli , A. 2016, , 591, A72

2016
[49]

2013, , 551, A75

Jacquet , E. 2013, , 551, A75

2013
[50]

2022, , 664, A86

Jang , H., Liu , B., & Johansen , A. 2022, , 664, A86

2022
[51]

M., & Carrasco-Gonz \'a lez , C

Jiang , H., Mac \' as , E., Guerra-Alvarado , O. M., & Carrasco-Gonz \'a lez , C. 2024, , 682, A32

2024
[52]

& Ormel , C

Jiang , H. & Ormel , C. W. 2023, , 518, 3877

2023
[53]

W., Krijt , S., & Dong , R

Jiang , H., Wang , Y., Ormel , C. W., Krijt , S., & Dong , R. 2023, , 678, A33

2023
[54]

S., Mac Low , M.-M., et al

Johansen , A., Oishi , J. S., Mac Low , M.-M., et al. 2007, , 448, 1022

2007
[55]

2021, Science Advances, 7, eabc0444

Johansen , A., Ronnet , T., Bizzarro , M., et al. 2021, Science Advances, 7, eabc0444

2021
[56]

& Youdin , A

Johansen , A. & Youdin , A. 2007, , 662, 627

2007
[57]

2009, , 704, L75

Johansen , A., Youdin , A., & Mac Low , M.-M. 2009, , 704, L75

2009
[58]

M., Alibert , Y., & San Sebasti \'a n , I

Kaufmann , N., Guilera , O. M., Alibert , Y., & San Sebasti \'a n , I. L. 2025, arXiv e-prints, arXiv:2502.02124

2025 arXiv
[59]

2013, Stellar Structure and Evolution

Kippenhahn , R., Weigert , A., & Weiss , A. 2013, Stellar Structure and Evolution

2013
[60]

A., et al

K \'o sp \'a l , \'A ., Cruz-S \'a enz de Miera , F., White , J. A., et al. 2021, , 256, 30

2021
[61]

W., Dominik , C., & Tielens , A

Krijt , S., Ormel , C. W., Dominik , C., & Tielens , A. G. G. M. 2016, , 586, A20

2016
[62]

& Johansen , A

Lambrechts , M. & Johansen , A. 2012, , 544, A32

2012
[63]

& Johansen , A

Lambrechts , M. & Johansen , A. 2014, , 572, A107

2014
[64]

J., Fuentes , J

Lee , E. J., Fuentes , J. R., & Hopkins , P. F. 2022, , 937, 95

2022
[65]

& Youdin , A

Li , R. & Youdin , A. N. 2021, , 919, 107

2021
[66]

N., & Simon , J

Li , R., Youdin , A. N., & Simon , J. B. 2019, , 885, 69

2019
[67]

Lichtenegger , H. I. M. & Komle , N. I. 1991, , 90, 319

1991
[68]

B., Li , R., et al

Lim , J., Simon , J. B., Li , R., et al. 2024, arXiv e-prints, arXiv:2410.17319

2024 arXiv
[69]

D., Li , Z.-Y., Tobin , J

Lin , Z.-Y. D., Li , Z.-Y., Tobin , J. J., et al. 2023, , 951, 9

2023
[70]

& Ji , J

Liu , B. & Ji , J. 2020, Research in Astronomy and Astrophysics, 20, 164

2020
[71]

2020, , 638, A88

Liu , B., Lambrechts , M., Johansen , A., Pascucci , I., & Henning , T. 2020, , 638, A88

2020
[72]

& Ormel , C

Liu , B. & Ormel , C. W. 2018, , 615, A138

2018
[73]

W., & Johansen , A

Liu , B., Ormel , C. W., & Johansen , A. 2019, , 624, A114

2019
[74]

& Pringle , J

Lynden-Bell , D. & Pringle , J. E. 1974, , 168, 603

1974
[75]

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

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

1977
[76]

2020, , 638, A1

Morbidelli , A. 2020, , 638, A1

2020
[77]

2019, , 872, 98

Mori , S., Bai , X.-N., & Okuzumi , S. 2019, , 872, 98

2019
[78]

2021, , 916, 72

Mori , S., Okuzumi , S., Kunitomo , M., & Bai , X.-N. 2021, , 916, 72

2021
[79]

& Wurm , G

Musiolik , G. & Wurm , G. 2019, , 873, 58

2019
[80]

1986, , 67, 375

Nakagawa , Y., Sekiya , M., & Hayashi , C. 1986, , 67, 375

1986
[81]

I., Facchini , S., & Anderson , D

\"O berg , K. I., Facchini , S., & Anderson , D. E. 2023, , 61, 287

2023
[82]

2011, , 738, 141

Oka , A., Nakamoto , T., & Ida , S. 2011, , 738, 141

2011
[83]

Ormel , C. W. & Klahr , H. H. 2010, , 520, A43

2010
[84]

Ormel , C. W. & Liu , B. 2018, , 615, A178

2018
[85]

W., Liu , B., & Schoonenberg , D

Ormel , C. W., Liu , B., & Schoonenberg , D. 2017, , 604, A1

2017
[86]

W., Shi , J.-M., & Kuiper , R

Ormel , C. W., Shi , J.-M., & Kuiper , R. 2015, , 447, 3512

2015
[87]

W., Spaans , M., & Tielens , A

Ormel , C. W., Spaans , M., & Tielens , A. G. G. M. 2007, , 461, 215

2007
[88]

Owen , J. E. 2020, , 495, 3160

2020
[89]

& Johansen , A

Ros , K. & Johansen , A. 2013, , 552, A137

2013
[90]

& Johansen , A

Ros , K. & Johansen , A. 2024, , 686, A237

2024
[91]

& Sirono , S.-i

Saito , E. & Sirono , S.-i. 2011, , 728, 20

2011
[92]

2017, , 597, A69

Sch \"a fer , U., Yang , C.-C., & Johansen , A. 2017, , 597, A69

2017
[93]

W., & Dorn , C

Schoonenberg , D., Liu , B., Ormel , C. W., & Dorn , C. 2019, , 627, A149

2019
[94]

Schoonenberg , D., Okuzumi , S., & Ormel , C. W. 2017, , 605, L2

2017
[95]

& Ormel , C

Schoonenberg , D. & Ormel , C. W. 2017, , 602, A21

2017
[96]

W., & Krijt , S

Schoonenberg , D., Ormel , C. W., & Krijt , S. 2018, , 620, A134

2018
[97]

Shakura , N. I. & Sunyaev , R. A. 1973, , 24, 337

1973
[98]

B., Armitage , P

Simon , J. B., Armitage , P. J., Li , R., & Youdin , A. N. 2016, , 822, 55

2016
[99]

L., Pommerol , A., et al

Spadaccia , S., Capelo , H. L., Pommerol , A., et al. 2022, , 509, 2825

2022
[100]

M., Tomida , K., White , C

Stone , J. M., Tomida , K., White , C. J., & Felker , K. G. 2020, , 249, 4

2020
[101]

& Lin , D

Takeuchi , T. & Lin , D. N. C. 2002, , 581, 1344

2002
[102]

J., van't Hoff , M

Tobin , J. J., van't Hoff , M. L. R., Leemker , M., et al. 2023, , 615, 227

2023
[103]

van 't Hoff , M. L. R., Tobin , J. J., Trapman , L., et al. 2018, , 864, L23

2018
[104]

2023, , 946, 70

Villenave , M., Podio , L., Duch \^e ne , G., et al. 2023, , 946, 70

2023
[105]

Visser , R. G. & Ormel , C. W. 2016, , 586, A66

2016
[106]

2013, , 559, A62

Wada , K., Tanaka , H., Okuzumi , S., et al. 2013, , 559, A62

2013
[107]

W., Huang , P., & Kuiper , R

Wang , Y., Ormel , C. W., Huang , P., & Kuiper , R. 2023, , 523, 6186

2023
[108]

Weidenschilling , S. J. 1977, , 180, 57

1977
[109]

& Kunz , M

Xu , W. & Kunz , M. W. 2021, , 508, 2142

2021
[110]

Xu , Z., Bai , X.-N., & Murray-Clay , R. A. 2017, , 847, 52

2017
[111]

2018, , 868, 27

Yang , C.-C., Mac Low , M.-M., & Johansen , A. 2018, , 868, 27

2018
[112]

Youdin , A. N. & Goodman , J. 2005, , 620, 459

2005
[113]

Youdin , A. N. & Lithwick , Y. 2007, , 192, 588

2007
[114]

2024, Reviews in Mineralogy and Geochemistry, 90, 27

Zhang , K. 2024, Reviews in Mineralogy and Geochemistry, 90, 27

2024
[115]

M., & Bai , X.-N

Zhu , Z., Stone , J. M., & Bai , X.-N. 2015, , 801, 81

2015

Pith tools

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

[1] [1]

, " * write output.state after.block = add.period write newline

ENTRY address archiveprefix author booktitle chapter edition editor howpublished institution eprint journal key month note number organization pages publisher school series title type volume year label extra.label sort.label short.list INTEGERS output.state before.all mid.sentence after.sentence after.block FUNCTION init.state.consts #0 'before.all := #1 ...

[2] [2]

write newline

" write newline "" before.all 'output.state := FUNCTION n.dashify 't := "" t empty not t #1 #1 substring "-" = t #1 #2 substring "--" = not "--" * t #2 global.max substring 't := t #1 #1 substring "-" = "-" * t #2 global.max substring 't := while if t #1 #1 substring * t #2 global.max substring 't := if while FUNCTION word.in bbl.in " " * FUNCTION format....

[3] [3]

P., Simon , J

Abod , C. P., Simon , J. B., Li , R., et al. 2019, , 883, 192

2019

[4] [4]

Armitage , P. J. 2020, Astrophysics of planet formation, Second Edition

2020

[5] [5]

& Wurm , G

Aumatell , G. & Wurm , G. 2011, , 418, L1

2011

[6] [6]

2017, , 845, 75

Bai , X.-N. 2017, , 845, 75

2017

[7] [7]

& Stone , J

Bai , X.-N. & Stone , J. M. 2010 a , , 722, 1437

2010

[8] [8]

& Stone , J

Bai , X.-N. & Stone , J. M. 2010 b , , 722, L220

2010