REVIEW 2 major objections 4 minor 1 cited by
Water Ice in the Edge-On Orion Silhouette Disk 114--426 from JWST NIRCam Images
T0 review · 2 major / 4 minor · reviewed 2026-08-11 · deepseek-v4-flash
Pith's one-line read JWST images of the edge-on Orion disk 114–426 reveal a 3-micron water-ice absorption dip in both scattered light and silhouette, with ice-to-refractory ratios up to ~0.2 and 0.46 Earth masses of solids in the outer ansae.
desk verdict The 3 micron water ice detection is solid and the paper deserves a real referee; the quantitative abundance and mass numbers, however, rest on an unquantified pure-absorption assumption and should be treated as model-dependent until scattering is properly handled. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing identity is the pure-absorption law $F_\nu/F_{\nu,\mathrm{bg}} = e^{-\kappa_{\mathrm{abs}}\Sigma_{\mathrm{dust}}}$ (Equation 1), which turns each pixel's measured dimming of the nebular background into a dust column density. The absorption opacity $\kappa_{\mathrm{abs}}$ is generated with a dust model of refractory cores bearing water-ice mantles, parametrized by the ice-to-refractory mass ratio, the maximum grain size (with a power-law size distribution of index $-3.5$ and a fixed minimum size of 0.005 μm), and the dust surface density $\Sigma_{\mathrm{dust}}$; a grid search minimizes $\chi^2$ over the eight uncontaminated bands. Two sets of bands are deliberately excluded from the fit: the Paschen α bands F182M and F187N, and the narrow bands F212N and F470N centered on molecular hydrogen transitions. The second mechanism in the paper is the survivability calculation: an equilibrium between photodesorption of ice by ultraviolet photons and re-adsorption of water vapor, solved for the ice fraction as a function of the external FUV field, gas density, and water and dust abundances, which the authors use to argue that the detected ice is consistent with the disk's radiation environment.
What would settle it
A spatially resolved spectrum at $R \gtrsim 1000$ across 2.5–3.6 μm would measure the true profile and depth of the 3 μm feature on the silhouette; if the water-ice column required to fit that spectrum disagreed with the band-derived abundances beyond the quoted uncertainties, or if the feature shape showed no water-ice structure, the photometric detection and its quantified ratios would fail. A separate test targets the survival argument: mapping the outer-disk gas density with resolved molecular-line observations, since the adsorption–photodesorption balance predicts that below $n_H \approx 10^8$ cm$^{-3}$ most water should be in vapor; a firmly low gas density with a persistent 3 μm dip would mean the icy grains are not steady-state residents of this disk.
Extended reading notes
Core claim
The paper's central claim is that the 3 μm dip in the disk's spectrum, present both in the scattered-light lobes and in the silhouette of the outer ansae against the Orion Nebula, is the vibrational absorption feature of water ice, and that the silhouette version can be quantified in terms of ice abundance. Each pixel's brightness relative to the nebular background is taken to obey $F_\nu/F_{\nu,\mathrm{bg}} = e^{-\kappa_{\mathrm{abs}}\Sigma_{\mathrm{dust}}}$, with scattering neglected because it is preferentially forward-directed. The opacity $\kappa_{\mathrm{abs}}$ comes from models of ice-mantled refractory grains (87% amorphous pyroxene and 13% amorphous carbon by mass), and a grid search over three free parameters, the ice-to-refractory mass ratio, the maximum grain size, and the dust surface density, fits eight of the twelve NIRCam bands at each pixel of the silhouette. The fits return ice-to-refractory mass ratios of 0 to 0.18, maximum grain sizes of 0.25 to 5 μm, and a total solid mass of 0.46 $M_\oplus$ in the two ansae, with the ice fraction and grain size declining with radius. The authors further report excess absorption in the two bands containing the Paschen α line, which they attribute to excited atomic hydrogen, and a wavelength-dependent brightness asymmetry between the scattered-light lobes that they interpret as evidence of a tilted inner disk. A separate equilibrium calculation between ultraviolet photodesorption and re-adsorption of water leads them to conclude that water ice can survive the external radiation field of this intermediate cluster environment.
Load-bearing premise
The load-bearing premise is that the dark outer edges of the disk dim the background nebula by absorption alone, with scattering so strongly forward that it does not push background photons out of our line of sight and with no light emitted by the disk itself or by material in front of it, so that all of the derived ice fractions, grain sizes, and masses rest on that single radiative assumption.
Editorial extensions
If this is right
- If the 3 μm dip is water ice as modeled, ice makes up roughly a fifth of the solid mass relative to refractory material in the outer ansae, a direct compositional constraint for models of where planets get their water.
- Grain sizes of 0.25–5 μm in the silhouette, together with the absence of millimeter grains there in ALMA data, show that outer-disk dust has grown past interstellar sizes without forming pebbles, implying stalled growth or inward drift of the largest grains.
- The excess absorption in the Paschen α bands points to a reservoir of hydrogen in the n = 3 state, likely pumped by Lyβ photons, making the silhouette a probe of gas excitation as well as dust composition.
- The wavelength-flipping lobe asymmetry and lateral offset of the scattered-light lobes imply a tilted inner disk, so the outer silhouette and the inner disk are not coplanar in this system.
- Water ice that survives external ultraviolet fluxes near $10^2$–$10^4$ G$_0$ would establish that ice can persist in the typical clustered, intermediate-radiation environment of star formation, not only in quiescent regions.
Reading between the lines
- Because the F300M filter cannot capture the full depth of the 3 μm feature, the derived ice-to-refractory ratios are arguably lower limits; a deeper or broader feature would raise them and would also explain the model's systematic tendency to overpredict the F335M band.
- The same silhouette technique should transfer to other edge-on disks that project onto bright H II region backgrounds, where NIRCam's multi-band coverage could map ice abundance in absorption without the scattering complications of scattered-light measurements.
- If the Paschen α absorption is confirmed spectroscopically, the disk's outer edge becomes a natural screen for mapping n = 3 hydrogen across a protoplanetary disk, a new tracer of ultraviolet-driven excitation that requires only narrowband imaging to exploit.
- A radiative-transfer treatment that includes forward scattering on these same lines of sight would quantify the systematic bias in the derived grain sizes and ice fractions, and would show whether the pure-absorption assumption pushes the ice-to-refractory ratio toward its ceiling of 0.18.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper presents JWST/NIRCam observations of the edge-on silhouette disk 114–426 in 12 bands spanning 1–5 μm. The authors report a dip at 3 μm in both the scattered-light lobes and the outer silhouette ansae, which they attribute to water ice. They fit the silhouette absorption spectra with a pure-absorption model of icy grains (Eq. 1) to derive ice-to-refractory mass ratios up to 0.18, maximum grain sizes of 0.25–5 μm, and a total solid mass of 0.46 M⊕ in the silhouette. They further report excess absorption in the Paα bands, wavelength-dependent lobe asymmetries that suggest a tilted inner disk, and a photodesorption analysis arguing that water ice can survive in the external radiation field of the Orion Nebula.
Significance. If the quantitative results hold, this is a valuable addition to the small sample of edge-on disks with direct ice detections, and the silhouette geometry offers a rare opportunity to measure ice abundances in absorption rather than through scattering models. The detection of the 3 μm dip is well supported by consistent F300M decrements across multiple pixels in both the scattered-light and silhouette regions, and the paper is commendably transparent about the limitations of its model, including the exclusion of Paα-contaminated bands and the imperfect fit to the red edge of the ice feature. The use of public JWST data, a documented reduction pipeline, and a standard dust-opacity code (optool) makes the analysis reproducible in principle.
major comments (2)
- [§3.3, Eq. (1)] The quantitative results—ice-to-refractory ratios, grain sizes, and the integrated mass—rest entirely on the pure-absorption approximation F/F_bg = exp(−κ_abs Σ_dust), with scattering opacity set to zero. The justification that scattering is preferentially forward-directed is qualitative and, as argued in the text, most directly applies to a point source behind a slab. Here the background is an extended nebula filling the field, so forward-scattered photons from neighboring background directions can replace photons removed from a given line of sight; whether net attenuation survives depends on the grain albedo, the phase function, and the pixel geometry, all of which vary with wavelength and ice mantle thickness. The paper offers no quantitative estimate of the resulting bias in the derived ice-to-refractory ratio, grain size, or surface density. A scattering-inclusive test—even a simple slab calculation with a Henyey–Greenstein phase function—is needed before the quoted abundance and mass values can be taken at face value. The detection of the 3 μm dip itself is on firmer ground, but the model-dependent numbers are not.
- [§3.3 and §4.1] The model systematically overpredicts the F335M band, which is included in the fit (only F182M, F187N, F212N, and F470N are excluded). The paper acknowledges in §4.1 that the red edge of the water-ice feature is not reproduced by the adopted ice optical constants and grain shape, and suggests that more sophisticated modeling (additional ice species, discrete-dipole approximation, RADMC-3D) would be needed. Given that F335M is one of the few bands sampling the long-wavelength wing of the feature, the derived ice-to-refractory ratios and maximum grain sizes are partly constrained by a band that the model fails to fit. The authors should either restrict the fit to the bands the model describes (e.g., F277W, F300M, F360M, F444W) and state the resulting parameter ranges, or explicitly present the current best-fit values as provisional pending a better red-edge treatment.
minor comments (4)
- [§3.3] In Eq. (1), Σ_dust is described as the dust surface density, but the model includes the ice mantle as part of the solid mass; renaming this quantity or clarifying that it is the total solid (dust plus ice) surface density would avoid confusion, especially when the total mass is summed in §3.3.
- [Figure 6] The three panels of Figure 6 lack visible color-bar labels or units in the manuscript text; if they appear in the published figure, please ensure the units (e.g., g cm^−2 for surface density) are legible.
- [§4.4, Eq. (3)] The photodesorption balance assumes a single-angle incident FUV field, but the paper later notes that photons entering from varying angles are neglected (citing Cleeves et al. 2013). A brief sentence acknowledging that Eq. (3) shares this simplification would improve consistency.
- [General typos] There are a few minor typographical issues, such as a missing space in '1.875µm' (p. 4) and a missing period after 'v_H2O = 216 m s^−1' (p. 12); these should be corrected in the final version.
Circularity Check
No significant circularity: the ice detection and fitted abundances rest on independent NIRCam photometry and externally tabulated ice opacities.
full rationale
The central claims are (1) a 3 micron dip in scattered-light and silhouette spectra, and (2) a model fit to that dip yielding ice-to-refractory ratios, grain sizes, and surface density. The dip is measured directly from the ratio of on-disk pixel brightness to a fitted nebular background in Section 3.2.1, not defined by the model. The fitted parameters in Section 3.3 enter through the absorption opacity in Equation (1), with optical constants from Warren & Brandt (2008) and other external laboratory sources; no fitted parameter is used to construct the target it is then said to predict. The pure-absorption treatment in Equation (1) is an approximation justified by forward scattering; whether that approximation biases the derived abundances is a modeling-validity concern, not circularity, because scattering opacity is not fitted from the same data and then renamed as a prediction. The self-citations (Ballering et al. 2021; McCaughrean & Pearson 2023) are used for pipeline description and for consistency with prior radiative-transfer predictions; the detection and abundance derivation do not reduce to those citations. No uniqueness theorem or ansatz is imported from the authors' prior work to force the water-ice interpretation. I therefore find no circular step.
Assumptions & free parameters
free parameters (4)
- ice_to_refractory_mass_ratio =
0 to 0.18
- maximum_grain_size =
0.25 to 5 micron
- dust_surface_density =
0.13 to 3.7 x 10^-4 g cm^-2
- nebular_background_plane =
fitted per image to an annular ellipse
assumptions (5)
- domain assumption Dust grains in the silhouette are ice-mantled refractory particles with a power-law size distribution n(a) proportional to a^-3.5 from 0.005 micron to a fitted maximum, and core composition 87% amorphous pyroxene plus 13% amorphous carbon.
- domain assumption Scattering opacity is negligible because scattering is preferentially forward-directed, so the transmitted fraction is exp(-kappa_abs Sigma_dust).
- domain assumption The disk lies in the foreground of the Orion Nebula and no significant emission arises from the silhouette material or intervening material.
- domain assumption The 3 micron dip is attributed primarily to water ice; NH3 and CH3OH ices can also contribute to the 3 micron feature.
- domain assumption The laboratory optical constants for water ice, pyroxene, and carbon used in the model are representative of the ice and dust in the disk.
Cite this review
Pith. "Pith review of Water Ice in the Edge-On Orion Silhouette Disk 114--426 from JWST NIRCam Images." pith.science (2026). https://pith.science/paper/ISPXHY6N
@misc{pith2026241204356,
author = {Pith},
title = {Pith review of: Water Ice in the Edge-On Orion Silhouette Disk 114--426 from JWST NIRCam Images},
year = {2026},
howpublished = {\url{https://pith.science/paper/ISPXHY6N}},
note = {Machine review of arXiv:2412.04356}
}
abstract
We examine images of the protoplanetary disk 114--426 with JWST/NIRCam in 12 bands. This large disk is oriented edge-on with a dark midplane flanked by lobes of scattered light. The outer edges of the midplane are seen in silhouette against the Orion Nebula, providing a unique opportunity to study planet-forming material in absorption. We discover a dip in the scattered light of the disk at 3\,$\micron$ -- compelling evidence for the presence of water ice. The 3\,$\micron$ dip is also seen in the silhouette of the disk, where we quantify the ice abundance with models of pure absorption and avoid the complications of disk scattering effects. We find grain ice-to-refractory mass ratios of up to $\sim$0.2, maximum grain sizes of 0.25 to 5\,$\micron$, and a total dust plus ice mass of 0.46\,$M_\oplus$ in the silhouette region. We also discover excess absorption in the NIRCam bands that include the Paschen $\alpha$ line, suggesting there may be excited atomic hydrogen in the disk. Examining the morphology of the scattered light lobes reveals that they are laterally offset from each other and exhibit a brightness asymmetry that flips with wavelength -- both evidence for a tilted inner disk in this system.
Figures
Figures from the paper (9 more)
Forward citations
Cited by 1 Pith paper
-
Spatial distribution of water ice in the protoplanetary silhouette disk d216-0939
A radiative transfer model of the edge-on disk d216-0939 indicates 5.4% crystalline water ice in the cold outer layers, implying outward transport of ice.
Reference graph
Works this paper leans on
-
[1]
Ballering, N. P., Cleeves, L. I., & Anderson, D. E. 2021, ApJ, 920, 115, doi: 10.3847/1538-4357/ac17ed
-
[2]
Bally, J., O’Dell, C. R., & McCaughrean, M. J. 2000, AJ, 119, 2919, doi: 10.1086/301385
doi:10.1086/301385 2000
-
[3]
Bally, J., Mann, R. K., Eisner, J., et al. 2015, ApJ, 808, 69, doi: 10.1088/0004-637X/808/1/69
-
[4]
Boogert, A. C. A., Huard, T. L., Cook, A. M., et al. 2011, ApJ, 729, 92, doi: 10.1088/0004-637X/729/2/92
-
[5]
Cleeves, L. I., Adams, F. C., & Bergin, E. A. 2013, ApJ, 772, 5, doi: 10.1088/0004-637X/772/1/5
-
[6]
Dartois, E., d’Hendecourt, L., Thi, W., Pontoppidan, K. M., & van Dishoeck, E. F. 2002, A&A, 394, 1057, doi: 10.1051/0004-6361:20021228
-
[7]
Dartois, E., Noble, J. A., Caselli, P., et al. 2024, Nature Astronomy, 8, 359, doi: 10.1038/s41550-023-02155-x
-
[8]
2021, OpTool: Command-line driven tool for creating complex dust opacities, Astrophysics Source Code Library, record ascl:2104.010
Dominik, C., Min, M., & Tazaki, R. 2021, OpTool: Command-line driven tool for creating complex dust opacities, Astrophysics Source Code Library, record ascl:2104.010. http://ascl.net/2104.010
2021
Show all 46 references
-
[9]
1995, A&A, 300, 503 Duchˆ ene, G., M´ enard, F., Stapelfeldt, K
Mutschke, H. 1995, A&A, 300, 503 Duchˆ ene, G., M´ enard, F., Stapelfeldt, K. R., et al. 2024, AJ, 167, 77, doi: 10.3847/1538-3881/acf9a7
1995 doi
-
[10]
P., Juhasz, A., Pohl, A., et al
Dullemond, C. P., Juhasz, A., Pohl, A., et al. 2012, RADMC-3D: A multi-purpose radiative transfer tool, Astrophysics Source Code Library. http://ascl.net/1202.015
2012
-
[11]
Fatuzzo, M., & Adams, F. C. 2008, ApJ, 675, 1361, doi: 10.1086/527469 W ater Ice in 114–426 15
2008 doi
-
[12]
2003, ApJ, 586, 296, doi: 10.1086/367557
Grosso, N., Alves, J., Wood, K., et al. 2003, ApJ, 586, 296, doi: 10.1086/367557
2003 doi
-
[13]
2015, ApJ, 798, 34, doi: 10.1088/0004-637X/798/1/34
Gundlach, B., & Blum, J. 2015, ApJ, 798, 34, doi: 10.1088/0004-637X/798/1/34
2015 doi
-
[14]
Habing, H. J. 1968, BAN, 19, 421
1968
-
[15]
Haworth, T. J. 2021, MNRAS, 503, 4172, doi: 10.1093/mnras/stab728
2021 doi
-
[16]
J., Bergin, E
Hollenbach, D., Kaufman, M. J., Bergin, E. A., & Melnick, G. J. 2009, ApJ, 690, 1497, doi: 10.1088/0004-637X/690/2/1497
2009 doi
-
[17]
J., & Lada, E
Lada, C. J., & Lada, E. A. 2003, ARA&A, 41, 57, doi: 10.1146/annurev.astro.41.011802.094844 Ma ´ ız Apell´ aniz, J., Barb´ a, R. H., Fern´ andez Aranda, R., et al. 2022, A&A, 657, A131, doi: 10.1051/0004-6361/202142364
2003 arXiv
-
[18]
K., Di Francesco, J., Johnstone, D., et al
Mann, R. K., Di Francesco, J., Johnstone, D., et al. 2014, ApJ, 784, 82, doi: 10.1088/0004-637X/784/1/82
2014 doi
-
[19]
P., & Dalle Ore, C
Cruikshank, D. P., & Dalle Ore, C. M. 2009, ApJ, 701, 1347, doi: 10.1088/0004-637X/701/2/1347
2009 doi
-
[20]
2011, ApJ, 727, 90, doi: 10.1088/0004-637X/727/2/90
McCabe, C., Duchˆ ene, G., Pinte, C., et al. 2011, ApJ, 727, 90, doi: 10.1088/0004-637X/727/2/90
2011 doi
-
[21]
J., & O’Dell, C
McCaughrean, M. J., & O’Dell, C. R. 1996, AJ, 111, 1977, doi: 10.1086/117934
1996 doi
- [22]
-
[23]
J., Chen, H., Bally, J., et al
McCaughrean, M. J., Chen, H., Bally, J., et al. 1998, ApJL, 492, L157, doi: 10.1086/311110
1998 doi
-
[24]
K., Rocha, W
McClure, M. K., Rocha, W. R. M., Pontoppidan, K. M., et al. 2023, Nature Astronomy, 7, 431, doi: 10.1038/s41550-022-01875-w
2023 doi
-
[25]
Miotello, A., Robberto, M., Potenza, M. A. C., & Ricci, L. 2012, ApJ, 757, 78, doi: 10.1088/0004-637X/757/1/78
2012 doi
-
[26]
A., Fraser, H
Noble, J. A., Fraser, H. J., Aikawa, Y., Pontoppidan, K. M., & Sakon, I. 2013, ApJ, 775, 85, doi: 10.1088/0004-637X/775/2/85 ¨Oberg, K. I., Linnartz, H., Visser, R., & van Dishoeck, E. F. 2009, ApJ, 693, 1209, doi: 10.1088/0004-637X/693/2/1209 O’Dell, C. R., & Wen, Z. 1994, Ap...
2013 doi
-
[27]
M., Dullemond, C
Pontoppidan, K. M., Dullemond, C. P., van Dishoeck, E. F., et al. 2005, ApJ, 622, 463, doi: 10.1086/427688
2005 doi
-
[28]
Rocha, W. R. M., Rachid, M. G., McClure, M. K., He, J., & Linnartz, H. 2024, A&A, 681, A9, doi: 10.1051/0004-6361/202347437
2024 doi
-
[29]
Y., Bally, J., Morris, M., & Throop, H
Shuping, R. Y., Bally, J., Morris, M., & Throop, H. 2003, ApJL, 587, L109, doi: 10.1086/375334
2003 doi
-
[30]
G., Sellgren, K., & Tokunaga, A
Smith, R. G., Sellgren, K., & Tokunaga, A. T. 1989, ApJ, 344, 413, doi: 10.1086/167809
1989 doi
-
[31]
A., McClure, M
Sturm, J. A., McClure, M. K., Beck, T. L., et al. 2023a, A&A, 679, A138, doi: 10.1051/0004-6361/202347512
-
[32]
A., McClure, M
Sturm, J. A., McClure, M. K., Bergner, J. B., et al. 2023b, A&A, 677, A18, doi: 10.1051/0004-6361/202346053
-
[33]
A., McClure, M
Sturm, J. A., McClure, M. K., Harsono, D., et al. 2024, A&A, 689, A92, doi: 10.1051/0004-6361/202450865
2024 doi
-
[34]
Terada, H., & Tokunaga, A. T. 2012, ApJ, 753, 19, doi: 10.1088/0004-637X/753/1/19
2012 doi
-
[35]
T., Kobayashi, N., et al
Terada, H., Tokunaga, A. T., Kobayashi, N., et al. 2007, ApJ, 667, 303, doi: 10.1086/520951
2007 doi
-
[36]
T., Pyo, T.-S., et al
Terada, H., Tokunaga, A. T., Pyo, T.-S., et al. 2012, AJ, 144, 175, doi: 10.1088/0004-6256/144/6/175
2012 doi
- [37]
-
[38]
B., Bally, J., Esposito, L
Throop, H. B., Bally, J., Esposito, L. W., & McCaughrean, M. J. 2001, Science, 292, 1686, doi: 10.1126/science.1059093
2001 doi
-
[39]
R., Duchˆ ene, G., et al
Villenave, M., Stapelfeldt, K. R., Duchˆ ene, G., et al. 2024, ApJ, 961, 95, doi: 10.3847/1538-4357/ad0c4b
2024 doi
-
[40]
G., & Brandt, R
Warren, S. G., & Brandt, R. E. 2008, Journal of Geophysical Research (Atmospheres), 113, D14220, doi: 10.1029/2007JD009744
2008 doi
- [41]
-
[42]
Weidenschilling, S. J. 1977, MNRAS, 180, 57, doi: 10.1093/mnras/180.2.57
1977 doi
-
[43]
J., & Haworth, T
Winter, A. J., & Haworth, T. J. 2022, European Physical Journal Plus, 137, 1132, doi: 10.1140/epjp/s13360-022-03314-1
2022 doi
-
[44]
2016, A&A, 586, A103, doi: 10.1051/0004-6361/201526538
Woitke, P., Min, M., Pinte, C., et al. 2016, A&A, 586, A103, doi: 10.1051/0004-6361/201526538
2016 doi
-
[45]
D., Pontoppidan, K
Yang, Y.-L., Green, J. D., Pontoppidan, K. M., et al. 2022, ApJL, 941, L13, doi: 10.3847/2041-8213/aca28910.48550/ arXiv.2208.10673
2022
-
[46]
G., Mennella, V., Colangeli, L., & Bussoletti, E
Zubko, V. G., Mennella, V., Colangeli, L., & Bussoletti, E. 1996, MNRAS, 282, 1321, doi: 10.1093/mnras/282.4.1321 16 Ballering et al. F115W 20 40 60 80 100 F140M 20 40 60 80 100 120 F162M 50 100 150 F182M 50 100 150 F187N 100 200 300 400 500 F212N 50 100 150 F277W 50 100 150 F...
1996 doi
Reviewed August 11, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.