The edge-on protoplanetary disk HH 48 NE. II. Modeling ices and silicates

Context. The abundance and distribution of ice in protoplanetary disks is critical for an understanding of the link between the composition of circumstellar matter and the composition of exoplanets. Edge-on protoplanetary disks are a useful tool for constraining this ice composition and its location in the disk because the spectral signatures of the ice can be observed in absorption against the continuum emission that arises from the warmer regions in the central disk.
Aims: The aim of this work is to model ice absorption features in protoplanetary disks and to determine how well the abundance of the main ice species throughout the disk can be determined within the uncertainty of the physical parameter space. The edge-on proto-planetary disk around HH 48 NE, a target of the James Webb Space Telescope Early Release program Ice Age, is used as a reference system.
Methods: We used the full anisotropic scattering capabilities of the radiative transfer code RADMC-3D to ray-trace the mid-infrared continuum. Using a constant parameterized ice abundance, we added ice opacities to the dust opacity in regions in which the disk was cold enough for the main carbon, oxygen, and nitrogen carriers to freeze out.
Results: The global abundance relative to the dust content of the main ice carriers in HH 48 NE can be determined within a factor of 3 when the uncertainty of the physical parameters is taken into account. Ice features in protoplanetary disks can be saturated at an optical depth of ≲1 due to local saturation. Ices are observed at various heights in the disk model, but in this model, spatial information is lost for features at wavelengths >7 µm when observing with James Webb Space Telescope because the angular resolution decreases towards longer wavelengths. Spatially observed ice optical depths cannot be directly related to column densities, as would be the case for direct absorption against a bright continuum source, because of radiative transfer effects. Vertical snowlines will not be a clear transition because the height of the snow surface increases radially, but their location may be constrained from observations using radiative transfer modeling. Radial snowlines are not really accessible. Not only the ice abundance, but also the inclination, the settling, the grain size distribution, and the disk mass have a strong impact on the observed ice absorption features in disks. Relative changes in the ice abundance can only be inferred from observations if the source structure is well constrained.