Variation of the disk thickness across ice bands: A method to determine ice abundances in highly inclined protoplanetary disks
The James Webb Space Telescope provides unprecedented information to study ices in protoplanetary disks. However, the saturation of ice bands in highly inclined disks hinders the measurement of ice abundances using classical spectroscopy. This is unfortunate as the presence and more importantly abundance of ices plays a key role in, e.g., the evolution of dust (because it modifies the sticking properties) and the composition of planetesimals and exoplanetary atmospheres. To overcome this issue and quantify the ice abundance within disks, we introduce a new method based on measuring the changes in the apparent disk thickness as a function of wavelength, which is directly and quantitatively related to the grain opacity. Specifically, we expect i) that the increased opacity within ice bands should result in a thicker disk than in the adjacent continuum, and ii) the thickness variations to be proportional to the abundance of ice. We extracted the disk thickness in model images of edge-on disks containing different abundances of water ice, as well as in James Webb Space Telescope spectral imaging of four edge-on disks. For both models and observations, the disk thickness decreases toward longer wavelengths except across the positions of ice absorption features where the thickness is enhanced across the band. In the model images, we demonstrate that this effect increases with ice abundance without any hint of saturation. This definitely demonstrates the presence of the ice species within each disk and confirms our expectation that this method can be applied to estimate ice abundances. Thanks to this method, it will thus be possible to constrain the ice abundance in highly inclined disks with disks model fitting. Unlike spectroscopic analysis, this method is not subject to saturation and should therefore be more robust and applicable to all disks for which the two surfaces can be resolved.