Planet formation in chemically diverse and evolving discs II. Chemical fingerprints in planetary atmospheres
Reviewed by Pith T0 review T1 audit T2 compute T3 formal T4 kernel 2026-07-02 17:15 UTCgrok-4.3pith:4P42XQ2Rrecord.jsonopen to challenge →
The pith
Giant planet atmospheres divide into three classes set by the balance of gas versus solid accretion during formation and migration.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Giant-planet atmospheres fall into three classes tied to accretion regime. Gas-dominated cases show N/O* > C/O* > C/N* with near-stellar or substellar S/N* and C/S*. Planetesimal-dominated cases show N/O* < C/O* < C/N*, S/N* >= C/N*, and C/S* <= C/O*. Drift-enhanced cases show N/O* < C/O* < C/N* together with markedly superstellar volatile-to-refractory ratios. N/O*, C/N*, and S/N* change systematically with migration distance while C/O* stays largely insensitive; degeneracies appear for formation beyond the CO and N2 snowlines. Disc chemistry and grain size further modulate the volatile-ratio patterns within each class.
What carries the argument
Coupling of 1D viscously evolving disc models (with radial dust drift and volatile chemistry) to N-body simulations of planetesimals interacting with a growing and migrating giant planet.
If this is right
- N/O*, C/N*, and S/N* vary systematically with the extent of planetary migration.
- C/O* alone provides little information on the gas-to-solid accretion balance.
- Metallicity does not uniquely indicate the solid-to-gas ratio in drift-dominated regimes.
- Sulphur ratios supply an independent constraint that breaks some degeneracies between formation pathways.
- Changes in disc chemical state or grain size produce distinctive volatile-ratio signatures within each atmospheric class.
Where Pith is reading between the lines
- Multi-element atmospheric measurements from JWST or Ariel could be inverted to estimate both the migration distance and the dominant grain size in the birth disc.
- The framework implies that planets forming beyond the main snowlines may require additional tracers such as phosphorus or noble-gas abundances to resolve formation-location ambiguities.
- If the three classes are confirmed, population-level statistics on atmospheric ratios could constrain the typical timing of giant-planet migration relative to disc dispersal.
Load-bearing premise
The one-dimensional disc models with dust drift and chemistry, when combined with N-body planetesimal tracking, correctly reproduce the time-dependent chemical composition of gas and solids that a planet accretes.
What would settle it
A survey of measured C/O*, N/O*, C/N*, and S/N* in a statistically significant sample of giant-planet atmospheres would fail to recover the three predicted ratio clusters if the classification does not hold.
Figures
read the original abstract
Giant planets form in protoplanetary discs, where the coupled dynamical and chemical evolution of gas and solids determines the composition of the material they accrete. We investigate how planet formation and migration shape the primordial elemental makeup of giant-planet atmospheres. Our aim is to link atmospheric compositions to planets' formation pathways and the time-dependent chemical properties of their natal discs. We couple 1D models of viscously evolving discs - incorporating radial dust drift and volatile chemistry - with N-body simulations of planetesimals interacting with a growing and migrating giant planet. Four chemical scenarios and three representative grain sizes (0.1, 20, and 100 micron) are explored. We track the accretion of carbon, oxygen, nitrogen, and sulphur to derive atmospheric elemental ratios normalised to stellar values (* denotes stellar normalisation). We identify three atmospheric classes corresponding to distinct accretion regimes: gas-dominated, characterised by N/O* > C/O* > C/N* and unconstrained or substellar S/N* (near-stellar C/S*); planetesimal-dominated, showing N/O* < C/O* < C/N*, S/N* >= C/N*, and C/S* <= C/O*; and drift-enhanced, exhibiting N/O* < C/O* < C/N* and markedly superstellar volatile-to-refractory ratios. N/O*, C/N*, and S/N* vary systematically with migration extent, although degeneracies arise for planets forming beyond the CO and N2 snowlines; C/O* remains largely insensitive. Metallicity alone does not uniquely trace the solid-to-gas accretion balance in drift-dominated regimes. Variations in the disc's chemical state and dust size imprint distinctive volatile-ratio patterns across these classes, providing complementary constraints on disc properties. This multi-element framework establishes predictive trends to guide the interpretation of atmospheric spectra from facilities like JWST and Ariel.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript couples 1D viscously evolving protoplanetary disc models (including radial dust drift and volatile chemistry) with N-body simulations of planetesimal accretion onto a growing and migrating giant planet. Across four chemical scenarios and three grain sizes, it derives normalized elemental ratios (N/O*, C/O*, C/N*, S/N*, C/S*) and identifies three atmospheric classes tied to distinct accretion regimes: gas-dominated (N/O* > C/O* > C/N*, near-stellar C/S*), planetesimal-dominated (N/O* < C/O* < C/N*, S/N* >= C/N*, C/S* <= C/O*), and drift-enhanced (N/O* < C/O* < C/N* with superstellar volatile-to-refractory ratios). The work examines systematic variations with migration extent and disc parameters while noting degeneracies beyond certain snowlines.
Significance. If the reported class distinctions hold under the modeled conditions, the multi-element ratio framework supplies concrete, simulation-derived diagnostics that link atmospheric compositions to formation pathways and natal disc properties. This goes beyond single-ratio metrics such as C/O and supplies testable trends for JWST and Ariel spectra, with explicit accounting for time-dependent chemistry and parameter degeneracies.
major comments (1)
- [Methods and Results] The central claim that the three classes emerge from the simulations rests on the coupling of 1D disc models to N-body accretion; the manuscript should quantify the sensitivity of the reported ratio orderings to this coupling (e.g., by comparing runs with and without radial drift or with varied initial chemical abundances) to demonstrate that the class boundaries are not artifacts of the 1D approximation.
minor comments (2)
- [Abstract] The abstract states the class characteristics but supplies no example numerical ratio values or overlap metrics; adding a compact table of median ratios (with ranges) for each class across the four scenarios would allow readers to assess the distinctness of the patterns directly.
- [Introduction or Results] Notation for stellar-normalized ratios is introduced with an asterisk but the precise definition (e.g., whether it is log or linear) should be restated once in the main text near the first results figure for clarity.
Simulated Author's Rebuttal
We thank the referee for the constructive comment on robustness. We address the major point below and will incorporate clarifications in a revised manuscript.
read point-by-point responses
-
Referee: [Methods and Results] The central claim that the three classes emerge from the simulations rests on the coupling of 1D disc models to N-body accretion; the manuscript should quantify the sensitivity of the reported ratio orderings to this coupling (e.g., by comparing runs with and without radial drift or with varied initial chemical abundances) to demonstrate that the class boundaries are not artifacts of the 1D approximation.
Authors: The manuscript already varies initial chemical abundances across four scenarios and modulates radial drift via three grain sizes (0.1, 20, and 100 micron), with the latter directly controlling drift velocities. The reported ratio orderings and three atmospheric classes remain consistent across this full grid, providing evidence that the distinctions are not artifacts of any single parameter choice. We did not run cases with radial drift disabled, as the model framework incorporates it as a core process; however, the grain-size variations serve as a controlled test of drift strength. We will add a short discussion paragraph noting this robustness while acknowledging that full 2D/3D disc models could provide further validation. revision: partial
Circularity Check
No significant circularity in simulation-based classification
full rationale
The paper derives its three atmospheric classes directly as outputs from forward simulations that couple 1D viscously evolving disc models (with radial drift and chemistry) to N-body planetesimal accretion tracking. Elemental ratios such as N/O*, C/O*, C/N*, S/N* and C/S* are computed from the simulated accretion histories across four chemical scenarios and three grain sizes; the classes are then identified by inspecting the resulting patterns. No parameters are fitted to data, no equations reduce by construction to their own inputs, and no load-bearing self-citations or imported uniqueness theorems are invoked. The derivation chain is therefore self-contained within the described numerical framework.
Axiom & Free-Parameter Ledger
free parameters (2)
- grain size
- chemical scenario
axioms (3)
- domain assumption Viscous evolution governs the radial transport of gas and dust in the disc
- domain assumption Volatile chemistry sets the partitioning of C, O, N, and S between gas and solids at different disc radii
- domain assumption N-body dynamics accurately capture planetesimal accretion onto a migrating giant planet
Reference graph
Works this paper leans on
-
[1]
2001, Earth and Planetary Science Letters, 185, 49
Allègre, C., Manhès, G., & Lewin, É. 2001, Earth and Planetary Science Letters, 185, 49
work page 2001
- [2]
-
[3]
2023, in Astronomical Society of the Pacific Conference Series, V ol
Bae, J., Isella, A., Zhu, Z., et al. 2023, in Astronomical Society of the Pacific Conference Series, V ol. 534, Protostars and Planets VII, ed. S. Inutsuka, Y . Aikawa, T. Muto, K. Tomida, & M. Tamura, 423
work page 2023
-
[4]
Bardyn, A., Baklouti, D., Cottin, H., et al. 2017, MNRAS, 469, S712
work page 2017
-
[5]
Bate, M. R. 2022, MNRAS, 514, 2145
work page 2022
-
[6]
Bergin, E. A., Blake, G. A., Ciesla, F., Hirschmann, M. M., & Li, J. 2015, Pro- ceedings of the National Academy of Science, 112, 8965 Bernabò, L. M., Turrini, D., Testi, L., Marzari, F., & Polychroni, D. 2022, ApJ, 927, L22
work page 2015
- [7]
- [8]
-
[9]
Bitsch, B., Schneider, A. D., & Kreidberg, L. 2022, A&A, 665, A138
work page 2022
- [10]
- [11]
-
[12]
Booth, R. A., Clarke, C. J., Madhusudhan, N., & Ilee, J. D. 2017, MNRAS, 469, 3994
work page 2017
-
[13]
Booth, R. A. & Ilee, J. D. 2019, MNRAS, 487, 3998
work page 2019
-
[14]
Burn, R. & Mordasini, C. 2025, Planetary Population Synthesis, ed. H. J. Deeg & J. A. Belmonte (Cham: Springer Nature Switzerland), 1–60
work page 2025
-
[15]
A., Lothringer, J., & Blake, G
Chachan, Y ., Knutson, H. A., Lothringer, J., & Blake, G. A. 2023, ApJ, 943, 112
work page 2023
-
[16]
Chambers, J. E. 1999, MNRAS, 304, 793
work page 1999
-
[17]
Chiang, E. I. & Goldreich, P. 1997, ApJ, 490, 368
work page 1997
-
[18]
J., Eistrup, C., & van Dishoeck, E
Cridland, A. J., Eistrup, C., & van Dishoeck, E. F. 2019, A&A, 627, A127
work page 2019
-
[19]
Cridland, A. J., Pudritz, R. E., Birnstiel, T., Cleeves, L. I., & Bergin, E. A. 2017, MNRAS, 469, 3910
work page 2017
-
[20]
Cridland, A. J., van Dishoeck, E. F., Alessi, M., & Pudritz, R. E. 2020, A&A, 642, A229
work page 2020
-
[21]
Crossfield, I. J. M. 2023, ApJ, 952, L18
work page 2023
-
[22]
Cuzzi, J. N. & Zahnle, K. J. 2004, ApJ, 614, 490
work page 2004
- [23]
-
[24]
Dawson, R. I. & Johnson, J. A. 2018, ARA&A, 56, 175 Dr˛ a˙ zkowska, J., Alibert, Y ., & Moore, B. 2016, A&A, 594, A105
work page 2018
-
[25]
Dullemond, C. P., Dominik, C., & Natta, A. 2001, ApJ, 560, 957
work page 2001
-
[26]
Eistrup, C., Walsh, C., & van Dishoeck, E. F. 2016, A&A, 595, A83
work page 2016
-
[27]
Eistrup, C., Walsh, C., & van Dishoeck, E. F. 2018, A&A, 613, A14
work page 2018
-
[28]
Emsenhuber, A., Mordasini, C., Burn, R., et al. 2021, A&A, 656, A69
work page 2021
-
[29]
Feinstein, A. D., Booth, R. A., Bergner, J. B., et al. 2025, arXiv e-prints, arXiv:2506.00669
-
[30]
D., Radica, M., Welbanks, L., et al
Feinstein, A. D., Radica, M., Welbanks, L., et al. 2023, Nature, 614, 670
work page 2023
-
[31]
Fortier, A., Alibert, Y ., Carron, F., Benz, W., & Dittkrist, K. M. 2013, A&A, 549, A44
work page 2013
-
[32]
Hahn, J. M. & Malhotra, R. 2005, AJ, 130, 2392
work page 2005
-
[33]
1981, Progress of Theoretical Physics Supplement, 70, 35
Hayashi, C. 1981, Progress of Theoretical Physics Supplement, 70, 35
work page 1981
-
[34]
Isella, A., Guidi, G., Testi, L., et al. 2016, Phys. Rev. Lett., 117, 251101
work page 2016
-
[35]
Johansen, A. & Lambrechts, M. 2017, Annual Review of Earth and Planetary Sciences, 45, 359
work page 2017
- [36]
-
[37]
Klahr, H. & Schreiber, A. 2016, in IAU Symposium, V ol. 318, Asteroids: New
work page 2016
-
[38]
Lambrechts, M. & Johansen, A. 2012, A&A, 544, A32 Le Gal, R., Öberg, K. I., Teague, R., et al. 2021, ApJS, 257, 12
work page 2012
-
[39]
Lebreuilly, U., Vallucci-Goy, V ., Guillet, V ., Lombart, M., & Marchand, P. 2023, MNRAS, 518, 3326
work page 2023
-
[40]
K., Nakajima, M., & Fischer, R
Lichtenberg, T., Schaefer, L. K., Nakajima, M., & Fischer, R. A. 2023, in As- tronomical Society of the Pacific Conference Series, V ol. 534, Protostars and Planets VII, ed. S. Inutsuka, Y . Aikawa, T. Muto, K. Tomida, & M. Tamura, 907
work page 2023
-
[41]
J., Hubickyj, O., D’Angelo, G., & Bodenheimer, P
Lissauer, J. J., Hubickyj, O., D’Angelo, G., & Bodenheimer, P. 2009, Icarus, 199, 338
work page 2009
-
[42]
2010, in Astrophysics and Space Science Proceedings, V ol
Lodders, K. 2010, in Astrophysics and Space Science Proceedings, V ol. 16, Prin- ciples and Perspectives in Cosmochemistry, 379
work page 2010
- [43]
- [44]
-
[45]
MacDonald, R. J. & Madhusudhan, N. 2017, ApJ, 850, L15
work page 2017
- [46]
-
[47]
Mah, J., Bitsch, B., Pascucci, I., & Henning, T. 2023, A&A, 677, L7
work page 2023
-
[48]
Marchand, P., Lebreuilly, U., Mac Low, M. M., & Guillet, V . 2023, A&A, 670, A61 Mollière, P., Molyarova, T., Bitsch, B., et al. 2022, ApJ, 934, 74
work page 2023
- [49]
-
[50]
Mordasini, C., Mollière, P., Dittkrist, K. M., Jin, S., & Alibert, Y . 2015, Interna- tional Journal of Astrobiology, 14, 201
work page 2015
-
[51]
Mordasini, C., van Boekel, R., Mollière, P., Henning, T., & Benneke, B. 2016, ApJ, 832, 41
work page 2016
- [52]
- [53]
- [54]
-
[55]
Pacetti, E., Schisano, E., Turrini, D., et al. 2025, A&A, 701, A194
work page 2025
- [56]
-
[57]
2014, in Planets, Asteriods, Comets and The Solar System, ed
Palme, H., Lodders, K., & Jones, A. 2014, in Planets, Asteriods, Comets and The Solar System, ed. A. M. Davis, V ol. 2, 15–36
work page 2014
-
[58]
Pirani, S., Johansen, A., Bitsch, B., Mustill, A. J., & Turrini, D. 2019, A&A, 623, A169
work page 2019
-
[59]
B., Hubickyj, O., Bodenheimer, P., et al
Pollack, J. B., Hubickyj, O., Bodenheimer, P., et al. 1996, Icarus, 124, 62
work page 1996
-
[60]
Rustamkulov, Z., Sing, D. K., Mukherjee, S., et al. 2023, Nature, 614, 659
work page 2023
-
[61]
Shakura, N. I. & Sunyaev, R. A. 1973, A&A, 24, 337
work page 1973
- [62]
-
[63]
2018, Experimental Astronomy, 46, 135
Tinetti, G., Drossart, P., Eccleston, P., et al. 2018, Experimental Astronomy, 46, 135
work page 2018
-
[64]
Tsai, S.-M., Lee, E. K. H., Powell, D., et al. 2023, Nature, 617, 483
work page 2023
-
[65]
2023, The Compositional Dimension of Planet Formation, ed
Turrini, D. 2023, The Compositional Dimension of Planet Formation, ed. L. M. Lara & D. Jewitt (World Scientific Publishing Co Pte Ltd), 1–47
work page 2023
-
[66]
2022, Experimental Astronomy, 53, 225
Turrini, D., Codella, C., Danielski, C., et al. 2022, Experimental Astronomy, 53, 225
work page 2022
-
[67]
2026, Astronomy and Computing, 55, 101064
Turrini, D., Fonte, S., Politi, R., et al. 2026, Astronomy and Computing, 55, 101064
work page 2026
-
[68]
Turrini, D., Marzari, F., Polychroni, D., & Testi, L. 2019, ApJ, 877, 50
work page 2019
- [69]
-
[70]
2024, A&A, 690, A23 V orobyov, E
Vallucci-Goy, V ., Lebreuilly, U., & Hennebelle, P. 2024, A&A, 690, A23 V orobyov, E. I., Elbakyan, V . G., Johansen, A., et al. 2023, A&A, 670, A81
work page 2024
- [71]
-
[72]
J., Morbidelli, A., Raymond, S
Walsh, K. J., Morbidelli, A., Raymond, S. N., O’Brien, D. P., & Mandell, A. M. 2011, Nature, 475, 206
work page 2011
-
[73]
Winn, J. N. & Fabrycky, D. C. 2015, ARA&A, 53, 409 Article number, page 14 of 20 E. Pacetti et al.: Planet formation in chemically diverse and evolving discs Fig. A.1: Gas-phase elemental abundances of C (top row), O (middle row), and N (bottom row) in the natal disc as a function of radius and time. All abundances are normalised to their respective stell...
work page 2015
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.