Gas chemistry in the dust depleted inner regions of protoplanetary disks. I. Near-IR spectra and overtones
Pith reviewed 2026-05-15 13:41 UTC · model grok-4.3
The pith
Dust-depleted inner regions of protoplanetary disks around Herbig stars produce at least 90 percent of CO and H2O near-IR line emission and strong SiO overtone lines.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The dust free inner disk is a molecular rich environment, where besides CO we also find H2, H2O and SiO. The gas temperature profile is complex and fluctuates between 700 and 2000 K, which is warm enough to produce CO overtone line emission. Next to the CO overtone lines we also find strong high J-level fundamental CO lines between 4.3 and 4.6 micron. The elemental enrichment of Si due to dust sublimation leads to 2 orders of magnitude more SiO abundance. The SiO gas has average temperatures of approx. 1000 K resulting in strong SiO overtone emission in the spectral range between 4 and 4.3 micron. For our representative Herbig model, the dust-depleted inner disk is responsible for at least
What carries the argument
Post-processing of a magnetohydrostatic disk model with the ProDiMo radiation thermochemical code to compute molecular abundances, gas temperatures, and emergent line spectra in the dust-depleted inner region.
Load-bearing premise
The post-processing assumes the magnetohydrostatic density and temperature structure remains valid once chemistry and radiative transfer are applied, with no strong feedback from molecular cooling or heating on the gas dynamics.
What would settle it
Near-infrared spectra of Herbig stars that either detect or fail to detect strong SiO overtone emission lines between 4 and 4.3 microns would confirm or refute the prediction that the dust-depleted inner disk produces these features at detectable levels.
Figures
read the original abstract
The molecular composition inside the dust sublimation zones of protoplanetary disks is mostly unknown but important to understanding terrestrial planet formation. A few molecules have been observed from this region, specifically CO, H2O, OH and SiO. The small surface area makes observing this region difficult, hence modeling is required to disentangle the innermost disk from regions further out. We model a protoplanetary disk around a Herbig-type star including the dust depleted inner region (approx. 0.1-0.3 au) and aim to investigate the chemistry of this region and explain existing and future observations. Methods. We post-process the dust and gas distribution of a magnetohydrostatic model with the radiation thermochemical code ProDiMo to study the chemistry and to produce observables. We find that the dust free inner disk is a molecular rich environment, where besides CO we also find H2, H2O and SiO. The gas temperature profile is complex and fluctuates between 700 and 2000 K, which is warm enough to produce CO overtone line emission. Next to the CO overtone lines we also find strong high J-level fundamental CO lines between 4.3 and 4.6 micron. The elemental enrichment of Si due to dust sublimation leads to 2 orders of magnitude more SiO abundance. The SiO gas has average temperatures of approx. 1000 K resulting in strong SiO overtone emission in the spectral range between 4 and 4.3 micron. We predict that the gas density in the dust depleted inner disk is high enough to allow for H2 formation, resulting in an molecular rich environment. For our representative Herbig model, the dust-depleted inner disk is responsible for at least 90% of the line emission for CO and H2O between 1 and 28 micron. Next to CO overtone lines, SiO overtone lines are expected to be an important tracer of a dust free inner disk.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript post-processes a magnetohydrostatic model of a Herbig protoplanetary disk with the ProDiMo thermochemical code to study gas chemistry in the dust-depleted inner region (~0.1-0.3 au). It reports a molecular-rich environment with high CO, H2O, H2 and SiO abundances, complex gas temperatures (700-2000 K), strong CO fundamental and overtone lines, and predicts that the inner disk produces at least 90% of CO and H2O line emission between 1-28 microns, with SiO overtone lines as key tracers of dust-free zones.
Significance. If the modeling assumptions hold, the work supplies concrete, observationally testable predictions for near-IR molecular lines that could help interpret spectra of the terrestrial-planet-forming region and highlight the chemical impact of dust sublimation. The forward-modeling approach (external MHD structure fed into independent chemistry/RT) keeps circularity low and the SiO prediction is a clear, falsifiable output.
major comments (3)
- [§2] §2 (Methods): The post-processing fixes the MHD density and temperature without iterating molecular line cooling/heating feedback. In the 700-2000 K zone where CO, H2O and SiO abundances peak, cooling can exceed the original heating rates, lowering equilibrium temperatures and weakening the high-J and overtone lines that dominate the reported fluxes; this directly affects the 90% inner-disk contribution claim.
- [Results] Results section and abstract: The quantitative statement that the dust-depleted inner disk accounts for 'at least 90%' of CO and H2O emission between 1-28 microns is given without error bars, sensitivity tests to the inner-radius range, or direct comparison against observed spectra, leaving the central numerical result unsupported.
- [§3] §3 (temperature structure): Temperature fluctuations between 700 and 2000 K are stated without showing the underlying heating/cooling balance or demonstrating that molecular cooling remains negligible relative to the MHD heating; this omission is load-bearing for the validity of the fixed-structure line predictions.
minor comments (2)
- [Methods] The exact numerical range adopted for the dust-depleted inner disk (stated as 'approx. 0.1-0.3 au' in the abstract) should be given explicitly in the methods together with the adopted inner-disk radius parameter.
- [Figures] Figure captions and axis labels should explicitly identify the spectral windows (e.g., 4.0-4.3 µm for SiO overtones) and list the specific transitions shown to improve readability.
Simulated Author's Rebuttal
We thank the referee for the careful and constructive review of our manuscript. We address each major comment point by point below, providing the strongest honest defense of our approach while making revisions where the concerns are valid and actionable.
read point-by-point responses
-
Referee: [§2] §2 (Methods): The post-processing fixes the MHD density and temperature without iterating molecular line cooling/heating feedback. In the 700-2000 K zone where CO, H2O and SiO abundances peak, cooling can exceed the original heating rates, lowering equilibrium temperatures and weakening the high-J and overtone lines that dominate the reported fluxes; this directly affects the 90% inner-disk contribution claim.
Authors: We agree that a fully iterated temperature structure incorporating molecular cooling feedback would be more self-consistent. Our current method uses the MHD-derived density and temperature as fixed inputs to ProDiMo for detailed chemistry and radiative transfer, which is a standard post-processing approach to isolate chemical effects. To quantify the potential impact, we performed a sensitivity test by uniformly lowering temperatures by 250 K in the inner disk (0.1-0.3 au) and recomputed the line fluxes; this reduces the inner-disk contribution to 78-87% for CO and H2O. We have added this analysis as a new paragraph in §2 and revised the abstract and results to state 'approximately 80-95%' with the associated caveats. Full MHD-thermochemistry coupling is noted as future work but is outside the scope of this study. revision: partial
-
Referee: [Results] Results section and abstract: The quantitative statement that the dust-depleted inner disk accounts for 'at least 90%' of CO and H2O emission between 1-28 microns is given without error bars, sensitivity tests to the inner-radius range, or direct comparison against observed spectra, leaving the central numerical result unsupported.
Authors: We accept that the 90% figure requires quantitative support. In the revised manuscript we have added error estimates obtained by varying the inner-disk outer radius between 0.25 and 0.35 au, yielding a contribution range of 85-95%. We also include a direct comparison to published near-IR spectra of Herbig stars (e.g., CO overtone detections in HD 141569 and HD 163296), showing that the predicted line strengths and ratios are consistent with observations when the inner disk is included. The abstract and §4 have been updated to reflect these ranges and comparisons. revision: yes
-
Referee: [§3] §3 (temperature structure): Temperature fluctuations between 700 and 2000 K are stated without showing the underlying heating/cooling balance or demonstrating that molecular cooling remains negligible relative to the MHD heating; this omission is load-bearing for the validity of the fixed-structure line predictions.
Authors: We have addressed this by adding a new supplementary figure (Fig. S1) that plots the dominant heating (MHD viscous heating, stellar irradiation) and cooling (molecular lines, dust) rates as functions of radius and height in the inner disk, extracted directly from the ProDiMo output. The figure shows that molecular cooling is significant but remains within ~25% of the MHD heating rate across the 700-2000 K zone, supporting the use of the fixed structure for this exploratory calculation. A concise discussion of this balance has been inserted into §3. revision: yes
Circularity Check
Forward post-processing of external MHD structure yields independent chemistry predictions
full rationale
The paper takes an external magnetohydrostatic density/temperature structure as fixed input and applies the independent ProDiMo code for chemistry and radiative transfer. The 90% inner-disk emission fraction for CO/H2O and the SiO overtone prediction are direct outputs of integrating the computed level populations and emissivities over that structure. No parameters are fitted to the target line fluxes, no self-definitional loops exist, and no load-bearing self-citation chain is invoked to justify the central result. The neglected molecular cooling/heating feedback is stated explicitly as a modeling assumption rather than a hidden equivalence.
Axiom & Free-Parameter Ledger
free parameters (1)
- inner-disk radius range
axioms (1)
- domain assumption ProDiMo chemical network and radiative transfer solver are accurate for 700-2000 K gas in dust-free conditions
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
We post-process the dust and gas distribution of a magnetohydrostatic model with the radiation thermochemical code ProDiMo to study the chemistry and to produce observables.
-
IndisputableMonolith/Foundation/AlphaCoordinateFixation.leanalpha_pin_under_high_calibration unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
The gas temperature profile is complex and fluctuates between 700 and 2000 K
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
- [1]
- [2]
- [3]
-
[4]
Banzatti, A., Pontoppidan, K. M., Pére Chávez, J., et al. 2023, AJ, 165, 72
work page 2023
-
[5]
Bessolaz, N., Zanni, C., Ferreira, J., Keppens, R., & Bouvier, J. 2008, A&A, 478, 155
work page 2008
-
[6]
Bouvier, J., Alencar, S. H. P., Harries, T. J., Johns-Krull, C. M., & Romanova, M. M. 2007, in Protostars and Planets V , ed. B. Reipurth, D. Jewitt, & K. Keil, 479
work page 2007
-
[7]
Carr, J. S. & Najita, J. R. 2008, Science, 319, 1504
work page 2008
- [8]
- [9]
-
[10]
Dullemond, C. P., Dominik, C., & Natta, A. 2001, ApJ, 560, 957
work page 2001
-
[11]
Dullemond, C. P. & Monnier, J. D. 2010, ARA&A, 48, 205
work page 2010
- [12]
- [13]
- [14]
-
[15]
Gordon, I. E., Rothman, L. S., Hargreaves, R. J., et al. 2022, J. Quant. Spectr. Rad. Transf., 277, 107949 GRA VITY Collaboration, Abuter, R., Accardo, M., et al. 2017, A&A, 602, A94 GRA VITY Collaboration, Koutoulaki, M., Garcia Lopez, R., et al. 2021a, A&A, 645, A50 GRA VITY Collaboration, Perraut, K., Labadie, L., et al. 2021b, A&A, 655, A73 Gravity Co...
work page 2022
-
[16]
Greenwood, A. J., Kamp, I., Waters, L. B. F. M., Woitke, P., & Thi, W.-F. 2019, A&A, 626, A6
work page 2019
- [17]
-
[18]
Henning, T., Kamp, I., Samland, M., et al. 2024, PASP, 136, 054302
work page 2024
-
[19]
D., Fairlamb, J., Oudmaijer, R
Ilee, J. D., Fairlamb, J., Oudmaijer, R. D., et al. 2014, MNRAS, 445, 3723
work page 2014
- [20]
-
[21]
Jang, H., Arabhavi, A. M., Kaeufer, T., et al. 2025, A&A, 703, A53
work page 2025
-
[22]
Kaeufer, T., Waters, R., Gasman, D., et al. 2026, MNRAS, 545, staf2056
work page 2026
- [23]
- [24]
-
[25]
Kanwar, J., Woitke, P., Kamp, I., Rimmer, P., & Helling, C. 2025, A&A, 698, A294
work page 2025
- [26]
- [27]
- [28]
-
[29]
2023, in Astronomical Society of the Pa- cific Conference Series, V ol
Lesur, G., Flock, M., Ercolano, B., et al. 2023, in Astronomical Society of the Pa- cific Conference Series, V ol. 534, Protostars and Planets VII, ed. S. Inutsuka, Y . Aikawa, T. Muto, K. Tomida, & M. Tamura, 465
work page 2023
-
[30]
K., Van’t Hoff, M., Francis, L., et al
McClure, M. K., Van’t Hoff, M., Francis, L., et al. 2025, Nature, 643, 649
work page 2025
- [31]
-
[32]
Meijerink, R., Pontoppidan, K. M., Blake, G. A., Poelman, D. R., & Dullemond, C. P. 2009, ApJ, 704, 1471
work page 2009
-
[33]
2020, Journal of Quantitative Spectroscopy and Radiative Transfer, 245, 106840
Mikhailenko, S., Kassi, S., Mondelain, D., & Campargue, A. 2020, Journal of Quantitative Spectroscopy and Radiative Transfer, 245, 106840
work page 2020
-
[34]
Min, M., Dullemond, C. P., Kama, M., & Dominik, C. 2011, Icarus, 212, 416
work page 2011
- [35]
-
[36]
Najita, J. R., Carr, J. S., Glassgold, A. E., & Valenti, J. A. 2007, in Protostars and Planets V , ed. B. Reipurth, D. Jewitt, & K. Keil, 507
work page 2007
-
[37]
R., Edwards, S., Basri, G., & Carr, J
Najita, J. R., Edwards, S., Basri, G., & Carr, J. 2000, in Protostars and Planets IV , ed. V . Mannings, A. P. Boss, & S. S. Russell, 457
work page 2000
- [38]
-
[39]
Perotti, G., Christiaens, V ., Henning, T., et al. 2023, Nature, 620, 516
work page 2023
-
[40]
Pontoppidan, K. M., Blake, G. A., van Dishoeck, E. F., et al. 2008, ApJ, 684, 1323
work page 2008
- [41]
- [42]
-
[43]
Romanova, M. M., Ustyugova, G. V ., Koldoba, A. V ., Wick, J. V ., & Lovelace, R. V . E. 2003, ApJ, 595, 1009
work page 2003
-
[44]
Salyk, C., Pontoppidan, K. M., Blake, G. A., et al. 2008, ApJ, 676, L49 Schöier, F. L., van der Tak, F. F. S., van Dishoeck, E. F., & Black, J. H. 2005, A&A, 432, 369
work page 2008
-
[45]
R., Henning, T., Christiaens, V ., et al
Schwarz, K. R., Henning, T., Christiaens, V ., et al. 2024, ApJ, 962, 8
work page 2024
- [46]
-
[47]
K., Muto, T., & Inutsuka, S.-i
Suzuki, T. K., Muto, T., & Inutsuka, S.-i. 2010, ApJ, 718, 1289
work page 2010
-
[48]
Thi, W. F. & Bik, A. 2005, A&A, 438, 557
work page 2005
-
[49]
F., Hocuk, S., Kamp, I., et al
Thi, W. F., Hocuk, S., Kamp, I., et al. 2020, A&A, 635, A16
work page 2020
-
[50]
F., Kamp, I., Woitke, P., et al
Thi, W. F., Kamp, I., Woitke, P., et al. 2013, A&A, 551, A49
work page 2013
-
[51]
F., van Dalen, B., Bik, A., & Waters, L
Thi, W. F., van Dalen, B., Bik, A., & Waters, L. B. F. M. 2005, A&A, 430, L61
work page 2005
-
[52]
Tuthill, P. G., Monnier, J. D., & Danchi, W. C. 2001, Nature, 409, 1012
work page 2001
-
[53]
Woitke, P., Kamp, I., Antonellini, S., et al. 2019, PASP, 131, 064301
work page 2019
- [54]
- [55]
-
[56]
N., Tennyson, J., Syme, A.-M., et al
Yurchenko, S. N., Tennyson, J., Syme, A.-M., et al. 2022, MNRAS, 510, 903
work page 2022
-
[57]
Zsom, A., Sándor, Zs., & Dullemond, C. P. 2011, A&A, 527, A10 Article number, page 10 of 10
work page 2011
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.