Search

This work focuses on the ice features toward the binary protostellar system Ced 110 IRS 4A and 4B, and observed with JWST as part of the Early Release Science Ice Age collaboration. We aim to explore the JWST observations of the binary protostellar system Ced~110~IRS4A and IRS4B to unveil and quantify the ice inventories toward these sources. We compare the ice abundances with those found for the same molecular cloud. The analysis is performed by fitting or comparing laboratory infrared spectra of ices to the observations. Spectral fits are carried out with the ENIIGMA fitting tool that searches for the best fit. For Ced~110~IRS4B, we detected the major ice species H$_2$O, CO, CO$_2$ and NH$_3$. All species are found in a mixture except for CO and CO$_2$, which have both mixed and pure ice components. In the case of Ced~110~IRS4A, we detected the same major species as in Ced~110~IRS4B, as well as the following minor species CH$_4$, SO$_2$, CH$_3$OH, OCN$^-$, NH$_4^+$ and HCOOH. Tentative detection of N$_2$O ice (7.75~$\mu$m), forsterite dust (11.2~$\mu$m) and CH$_3^+$ gas emission (7.18~$\mu$m) in the primary source are also presented. Compared with the two lines of sight toward background stars in the Chameleon I molecular cloud, the protostar has similar ice abundances, except in the case of the ions that are higher in IRS4A. The clearest differences are the absence of the 7.2 and 7.4~$\mu$m absorption features due to HCOO$^-$ and icy complex organic molecules in IRS4A and evidence of thermal processing in both IRS4A and IRS4B as probed by the CO$_2$ ice features. We conclude that the binary protostellar system Ced~110~IRS4A and IRS4B has a large inventory of icy species. The similar ice abundances in comparison to the starless regions in the same molecular cloud suggest that the chemical conditions of the protostar were set at earlier stages in the molecular cloud., Comment: 33 pages, 19 Figures. Accepted for publication in Astronomy & Astrophysics

The MIRI-MRS instrument onboard JWST allows for probing the molecular gas composition at mid-IR wavelengths at unprecedented resolution and sensitivity. It is important to study these features in low-mass embedded protostellar systems since the formation of planets is thought to start in this phase. We present JWST/MIRI-MRS data of 18 low-mass protostellar systems in the JOYS program, focusing on gas-phase molecular lines in spectra extracted from the central protostellar positions. Besides H2, the most commonly detected molecules are H2O, CO2, CO, and OH. Other molecules such as 13CO2, C2H2, 13CCH, HCN, C4H2, CH4, and SO2 are detected only toward at most three of the sources. The JOYS data also yield the surprising detection of SiO gas toward two sources (BHR71-IRS1, L1448-mm) and for the first time CS and NH3 at mid-IR wavelengths toward a low-mass protostar (B1-c). The temperatures derived for the majority of the molecules are 100-300 K, much lower than what is typically derived toward more evolved Class II sources (>500 K). Toward three sources (e.g., TMC1-W), hot (~1000 K) H2O is detected, indicative of the presence of hot molecular gas in the embedded disks, but such warm emission from other molecules is absent. The agreement in abundance ratios with respect to H2O between ice and gas point toward ice sublimation in a hot core for a few sources (e.g., B1-c) whereas their disagreement and velocity offsets hint at high-temperature (shocked) conditions toward other sources (e.g., L1448-mm, BHR71-IRS1). The typical temperatures of the gas-phase molecules of 100-300 K are consistent with both ice sublimation in hot cores as well as high-temperature gas phase chemistry. Molecular features originating from the inner embedded disks are not commonly detected, likely because they are too extincted even at mid-IR wavelengths by small not-settled dust grains in upper layers of the disk., Comment: 25 pages, 11 figures, accepted for publication in A&A

The carbon isotope ratio is a powerful tool for studying the evolution of stellar systems. Recent detections of CO isotopologues in disks and exoplanet atmospheres pointed towards significant fractionation in these systems. In order to understand the evolution of this quantity, it is crucial to trace the isotope abundance from stellar nurseries to planetary systems. During the protostellar stage the multiple vibrational modes of CO$_2$ and CO ice provide a unique opportunity to examine the carbon isotope ratio in the solid state. Now with the sensitivity of the \textit{James Webb Space Telescope}, these absorption features have become accessible at high S/N in Solar-mass systems. We quantify the $^{12}$CO$_2$/$^{13}$CO$_2$ and the $^{12}$CO/$^{13}$CO isotope ratios in 17 class 0/I low mass protostars from the $^{12}$CO$_2$ combination modes (2.70 ${\mu}$m and 2.77 ${\mu}$m), the $^{12}$CO$_2$ stretching mode (4.27 ${\mu}$m), the $^{13}$CO$_2$ stretching mode (4.39 ${\mu}$m), the $^{12}$CO$_2$ bending mode (15.2 ${\mu}$m), the $^{12}$CO stretching mode (4.67 ${\mu}$m) and the $^{13}$CO stretching mode (4.78 ${\mu}$m) using JWST observations. We also report a detection of the $^{12}$CO overtone mode at 2.35 ${\mu}$m. The $^{12}$CO$_2$/$^{13}$CO$_2$ ratios are in agreement and we find mean ratios of 85 $\pm$ 23, 76 $\pm$ 12 and 97 $\pm$ 17 for the 2.70 ${\mu}$m, 4.27 ${\mu}$m and the 15.2 ${\mu}$m bands, respectively. The main source of uncertainty stem from the error on the band strengths. The $^{12}$CO/$^{13}$CO ratios derived from the 4.67 ${\mu}$m bands are consistent, albeit elevated with respect to the $^{12}$CO$_2$/$^{13}$CO$_2$ ratios and we find a mean ratio of 165 $\pm$ 52. These findings indicate that ices leave the pre-stellar stage with elevated carbon isotope ratios relative to the interstellar medium and that fractionation becomes significant during the later stages.

Planet formation is strongly influenced by the composition and distribution of volatiles within protoplanetary disks. With JWST, it is now possible to obtain direct observational constraints on disk ices, as recently demonstrated by the detection of ice absorption features towards the edge-on HH 48 NE disk as part of the Ice Age Early Release Science program. Here, we introduce a new radiative transfer modeling framework designed to retrieve the composition and mixing status of disk ices using their band profiles, and apply it to interpret the H2O, CO2, and CO ice bands observed towards the HH 48 NE disk. We show that the ices are largely present as mixtures, with strong evidence for CO trapping in both H2O and CO2 ice. The HH 48 NE disk ice composition (pure vs. polar vs. apolar fractions) is markedly different from earlier protostellar stages, implying thermal and/or chemical reprocessing during the formation or evolution of the disk. We infer low ice-phase C/O ratios around 0.1 throughout the disk, and also demonstrate that the mixing and entrapment of disk ices can dramatically affect the radial dependence of the C/O ratio. It is therefore imperative that realistic disk ice compositions are considered when comparing planetary compositions with potential formation scenarios, which will fortunately be possible for an increasing number of disks with JWST., Comment: Accepted to ApJ. 24 pages, 15 figures

A rich inventory of complex organic molecules (COMs) has been observed in high abundances in the gas phase toward Class 0 protostars. These molecules are suggested to be formed in ices and sublimate in the warm inner envelope close to the protostar. However, only the most abundant COM, methanol (CH3OH), has been firmly detected in ices before the era of James Webb Space Telescope (JWST). Now it is possible to detect the interstellar ices of other COMs and constrain their ice column densities quantitatively. We aim to determine the column densities of several oxygen-bearing COMs (O-COMs) in both gas and ice for two low-mass protostellar sources, NGC 1333 IRAS 2A and B1-c, as case studies in our JWST Observations of Young protoStars (JOYS+) program. By comparing the column density ratios w.r.t. CH3OH between both phases measured in the same sources, we can probe into the evolution of COMs from ice to gas in the early stages of star formation. We are able to fit the fingerprints range of COM ices between 6.8 and 8.8 um in the JWST/MIRI-MRS spectra of B1-c using similar components as recently used for IRAS 2A. We claim detection of CH4, OCN-, HCOO-, HCOOH, CH3CHO, C2H5OH, CH3OCH3, CH3OCHO, and CH3COCH3 in B1-c, and upper limits are estimated for SO2, CH3COOH, and CH3CN. The comparison of O-COM ratios w.r.t CH3OH between ice and gas shows two different cases. 1) the column density ratios of CH3OCHO and CH3OCH3 match well between the two phases, which may be attributed to a direct inheritance from ice to gas or strong chemical links with CH3OH. 2) the ice ratios of CH3CHO and C2H5OH w.r.t. CH3OH are higher than the gas ratios by 1-2 orders of magnitudes. This difference can be explained by the gas-phase reprocessing following sublimation, or different spatial distributions of COMs in the envelope., Comment: 42 pages (22 main text, 20 appendix); 27 figures (12 in main text, 15 in appendix); 5 tables (2 in main text, 3 in appendix) Accepted for publication in A&A

Ice-coated dust grains provide the main reservoir of volatiles that play an important role in planet formation processes and may become incorporated into planetary atmospheres. However, due to observational challenges, the ice abundance distribution in protoplanetary disks is not well constrained. We present JWST/MIRI observations of the edge-on disk HH 48 NE carried out as part of the IRS program Ice Age. We detect CO$_2$, NH$_3$, H$_2$O and tentatively CH$_4$ and NH$_4^+$. Radiative transfer models suggest that ice absorption features are produced predominantly in the 50-100 au region of the disk. The CO$_2$ feature at 15 micron probes a region closer to the midplane (z/r = 0.1-0.15) than the corresponding feature at 4.3 micron (z/r = 0.2-0.6), but all observations trace regions significantly above the midplane reservoirs where we expect the bulk of the ice mass to be located. Ices must reach a high scale height (z/r ~ 0.6; corresponding to modeled dust extinction Av ~ 0.1), in order to be consistent with the observed vertical distribution of the peak ice optical depths. The weakness of the CO$_2$ feature at 15 micron relative to the 4.3 micron feature and the red emission wing of the 4.3 micron CO$_2$ feature are both consistent with ices being located at high elevation in the disk. The retrieved NH$_3$ abundance and the upper limit on the CH$_3$OH abundance relative to H$_2$O are significantly lower than those in the interstellar medium (ISM), but consistent with cometary observations. Full wavelength coverage is required to properly study the abundance distribution of ices in disks. To explain the presence of ices at high disk altitudes, we propose two possible scenarios: a disk wind that entrains sufficient amounts of dust, thus blocking part of the stellar UV radiation, or vertical mixing that cycles enough ices into the upper disk layers to balance ice photodesorption., Comment: 16 pages, 11 figures, accepted in A&A

Nitrogen-bearing complex organic molecules have been commonly detected in the gas phase but not yet in interstellar ices. This has led to the long-standing question of whether these molecules form in the gas phase or in ices. $\textit{James Webb}$ Space Telescope ($\textit{JWST}$) offers the sensitivity, spectral resolution, and wavelength coverage needed to detect them in ices and investigate whether their abundance ratios are similar in gas and ice. We report the first tentative detection of CH$_3$CN, C$_2$H$_5$CN, and the simple molecule, N$_2$O, based on the CN-stretch band in interstellar ices toward three (HOPS 153, HOPS 370, and IRAS 20126+4104) out of the five protostellar systems observed as part of the Investigating Protostellar Accretion (IPA) GO program with $\textit{JWST}$-NIRSpec. We also provide upper limits for the two other sources with smaller luminosities in the sample. We detect OCN$^-$ in the ices of all sources with typical CH$_3$CN/OCN$^-$ ratios of around 1. Ice and gas column density ratios of the nitrogen-bearing species with respect to each other are better matched than those with respect to methanol, which are a factor of ${\sim}5$ larger in the ices than the gas. We attribute the elevated ice column densities with respect to methanol to the difference in snowline locations of nitrogen-bearing molecules and of methanol, biasing the gas-phase observations toward fewer nitrogen-bearing molecules. Moreover, we find tentative evidence for enhancement of OCN$^-$, CH$_3$CN, and C$_2$H$_5$CN in warmer ices, although formation of these molecules likely starts along with methanol in the cold prestellar phase. Future surveys combining NIRSpec and MIRI, and additional laboratory spectroscopic measurements of C$_2$H$_5$CN ice, are necessary for robust detection and conclusions on the formation history of complex cyanides., Comment: Accepted for publication in A&A

Complex organic molecules (COMs) detected in the gas phase are thought to be mostly formed on icy grains, but no unambiguous detection of icy COMs larger than CH3OH has been reported so far. Exploring this matter in more detail has become possible with the JWST the critical 5-10 $\mu$m range. In the JOYS+ program, more than 30 protostars are being observed with the MIRI/MRS. This study explores the COMs ice signatures in the low and high-mass protostar, IRAS 2A and IRAS 23385, respectively. We fit continuum and silicate subtracted observational data with IR laboratory ice spectra. We use the ENIIGMA fitting tool to find the best fit between the lab data and the observations and to performs statistical analysis of the solutions. We report the best fits for the spectral ranges between 6.8 and 8.6 $\mu$m in IRAS 2A and IRAS 23385, originating from simple molecules, COMs, and negative ions. The strongest feature in this range (7.7 $\mu$m) is dominated by CH4 and has contributions of SO2 and OCN-. Our results indicate that the 7.2 and 7.4 $\mu$m bands are mostly dominated by HCOO-. We find statistically robust detections of COMs based on multiple bands, most notably CH3CHO, CH3CH2OH, and CH3OCHO. The likely detection of CH3COOH is also reported. The ice column density ratios between CH3CH2OH and CH3CHO of IRAS 2A and IRAS 23385, suggests that these COMs are formed on icy grains. Finally, the derived ice abundances for IRAS 2A correlate well with those in comet 67P/GC within a factor of 5. Based on the MIRI/MRS data, we conclude that COMs are present in interstellar ices, thus providing additional proof for a solid-state origin of these species in star-forming regions. The good correlation between the ice abundances in comet 67P and IRAS 2A is in line with the idea that cometary COMs can be inherited from the early protostellar phases., Comment: Accepted for publication in A&A

(abridged) We study and compare the warm (>100 K) and cold (<100 K) material toward the high-mass star-forming region IRAS 23385+6053 (IRAS 23385 hereafter) combining high angular resolution observations in the mid-infrared (MIR) with the JWST Observations of Young protoStars (JOYS) project and with the NOEMA at mm wavelengths at angular resolutions of 0.2"-1". The spatial morphology of atomic and molecular species is investigated by line integrated intensity maps. The temperature and column density of different gas components is estimated using H2 transitions (warm and hot component) and a series of CH3CN transitions as well as 3 mm continuum emission (cold component). Toward the central dense core in IRAS 23385 the material consists of relatively cold gas and dust (~50 K), while multiple outflows create heated and/or shocked H2 and show enhanced temperatures (~400 K) along the outflow structures. An energetic outflow with enhanced emission knots of [Fe II] and [Ni II] hints at J-type shocks, while two other outflows have enhanced emission of only H2 and [S I] caused by C-type shocks. The latter two outflows are also more prominent in molecular line emission at mm wavelengths (e.g., SiO, SO, H2CO, and CH3OH). Even higher angular resolution data are needed to unambiguously identify the outflow driving sources given the clustered nature of IRAS 23385. While most of the forbidden fine structure transitions are blueshifted, [Ne II] and [Ne III] peak at the source velocity toward the MIR source A/mmA2 suggesting that the emission is originating from closer to the protostar., Comment: 15 pages, 7 figures, accepted for publication in A&A

Ices are the main carriers of volatiles in protoplanetary disks and are crucial to our understanding of the chemistry that ultimately sets the organic composition of planets. The ERS program Ice Age on the JWST follows the ice evolution through all stages of star and planet formation. JWST/NIRSpec observations of the edge-on Class II protoplanetary disk HH~48~NE reveal spatially resolved absorption features of the major ice components H$_2$O, CO$_2$, CO, and multiple weaker signatures from less abundant ices NH$_3$, OCN$^-$, and OCS. Isotopologue $^{13}$CO$_2$ ice has been detected for the first time in a protoplanetary disk. Since multiple complex light paths contribute to the observed flux, the ice absorption features are filled in by ice-free scattered light. The $^{12}$CO$_2$/$^{13}$CO$_2$ ratio of 14 implies that the $^{12}$CO$_2$ feature is saturated, without the flux approaching 0, indicative of a very high CO$_2$ column density on the line of sight, and a corresponding abundance with respect to hydrogen that is higher than ISM values by a factor of at least a few. Observations of rare isotopologues are crucial, as we show that the $^{13}$CO$_2$ observation allows us to determine the column density of CO$_2$ to be at an order of magnitude higher than the lower limit directly inferred from the observed optical depth. Radial variations in ice abundance, e.g., snowlines, are significantly modified since all observed photons have passed through the full radial extent of the disk. CO ice is observed at perplexing heights in the disk, extending to the top of the CO-emitting gas layer. We argue that the most likely interpretation is that we observe some CO ice at high temperatures, trapped in less volatile ices like H$_2$O and CO$_2$. Future radiative transfer models will be required to constrain the implications on our current understanding of disk physics and chemistry., Comment: 16 pages, 8 figures, accepted for publication in A&A

During the journey from the cloud to the disc, the chemical composition of the protostellar envelope material can be either preserved or processed to varying degrees depending on the surrounding physical environment. This works aims to constrain the interplay of solid (ice) and gaseous methanol (CH$_3$OH) in the outer regions of protostellar envelopes located in the Coronet cluster in Corona Australis (CrA), and assess the importance of irradiation by the Herbig Ae/Be star R CrA. CH$_3$OH is a prime test-case as it predominantly forms as a consequence of the solid-gas interplay (hydrogenation of condensed CO molecules onto the grain surfaces) and it plays an important role in future complex molecular processing. We present 1.3 mm Submillimeter Array (SMA) and Atacama Pathfinder Experiment (APEX) observations towards the envelopes of four low-mass protostars in the Coronet. Eighteen molecular transitions of seven species are identified. We calculate CH$_3$OH gas-to-ice ratios in this strongly irradiated cluster and compare them with ratios determined towards protostars located in less irradiated regions such as the Serpens SVS 4 cluster in Serpens Main and the Barnard 35A cloud in the $\lambda$ Orionis region. The CH$_3$OH gas-to-ice ratios in the Coronet vary by one order of magnitude (from 1.2$\times$10$^{-4}$ to 3.1$\times$10$^{-3}$) which is similar to less irradiated regions as found in previous studies. We find that the CH$_3$OH gas-to-ice ratios estimated in these three regions are remarkably similar despite the different UV radiation field intensities and formation histories. This result suggests that the overall CH$_3$OH chemistry in the outer regions of low-mass envelopes is relatively independent of variations in the physical conditions and hence that it is set during the prestellar stage., Comment: 20 pages, 17 figures, accepted for publication in Astronomy & Astrophysics

Interstellar ices are largely composed of frozen water. It is important to derive fundamental parameters for H$_2$O ice such as absorption and scattering opacities for which accurate complex refractive indexes are needed. The primary goal of this work is to derive ice-grain opacities based on accurate H$_2$O ice complex refractive indexes and to assess their impact on the derivation of ice column densities and porosity in space. We use the \texttt{optool} code to derive ice-grain opacities values based on new mid-IR complex refractive index measurements of H$_2$O ice. Next, we use those opacities in the \texttt{RADMC-3D} code to run a radiative transfer simulation of a protostellar envelope containing H$_2$O ice. This is used to calculate water ice column densities. We find that the real refractive index in the mid-IR of H$_2$O ice at 30~K is $\sim$14\% lower than previously reported in the literature. This has a direct impact on the ice column densities derived from the simulations of embedded protostars. We find that ice porosity plays a significant role in the opacity of icy grains and that the H$_2$O libration mode can be used as a diagnostic tool to constrain the porosity level. Finally, the refractive indexes presented here allow us to estimate a grain size detection limit of 18~$\mu$m based on the 3~$\mu$m band whereas the 6~$\mu$m band allows tracing grain sizes larger than 20~$\mu$m. Based on radiative transfer simulations using new mid-IR refractive indexes, we conclude that H$_2$O ice leads to more absorption of infrared light than previously estimated. This implies that the 3 and 6~$\mu$m bands remain detectable in icy grains with sizes larger than 10~$\mu$m. Finally, we propose that also the H$_2$O ice libration band can be a diagnostic tool to constrain the porosity level of the interstellar ice, in addition to the OH dangling bond, which is routinely used for this purpose., Comment: Accepted for publication in A&A, 12 pages, 15 figures

The abundance and distribution of ice in protoplanetary disks (PPD) is critical to understand the linkage between the composition of circumstellar matter and the composition of exoplanets. Edge-on PPDs are a useful tool to constrain such ice composition and its location in the disk, as ice spectral signatures can be observed in absorption against the continuum emission arising from the warmer central disk regions. The aim of this work is to model ice absorption features in PPDs and determine how well the abundance of the main ice species across the disk can be determined within the uncertainty of the physical parameter space. The edge-on PPD around HH 48 NE, a target of the JWST ERS program IceAge, is used as a reference system. We use RADMC-3D to raytrace the mid-infrared continuum. Using a constant parameterized ice abundance, ice opacities are added to the dust opacity in regions wherever the disk is cold enough for the main carbon, oxygen and nitrogen carriers to freeze out. The global abundance of the main ice carriers in HH 48 NE can be determined within a factor of 3, when taking the uncertainty of the physical parameters into account. Ice features in PPDs can be saturated at an optical depth <1, due to local saturation. Spatially observed ice optical depths cannot be directly related to column densities due to radiative transfer effects. Vertical snowlines will not be a clear transition due to the radially increasing height of the snowsurface, but their location may be constrained from observations using radiative transfer modeling. Radial snowlines are not really accesible. Not only the ice abundance, but also inclination, settling, grain size distribution and disk mass have strong impact on the observed ice absorption features in disks. Relative changes in ice abundance can be inferred from observations only if the source structure is well constrained, Comment: 18 pages, 12 figures, accepted for publication in Astronomy & Astrophysics

Context. Observations of edge-on disks are an important tool for constraining general protoplanetary disk properties that cannot be determined in any other way. However, most radiative transfer models cannot simultaneously reproduce the spectral energy distributions (SEDs) and resolved scattered light and submillimeter observations of these systems, due to the differences in geometry and dust properties at different wavelengths. Aims. We simultaneously constrain the geometry of the edge-on protoplanetary disk HH 48 NE and the characteristics of the host star. HH 48 NE is part of the JWST early release science program Ice Age. This work serves as a stepping stone towards a better understanding of the disk physical structure and icy chemistry in this particular source. This kind of modeling lays the groundwork for studying other edge-on sources to be observed with the JWST. Methods. We fit a parameterized dust model to HH 48 NE by coupling the radiative transfer code RADMC-3D and an MCMC framework. The dust structure was fitted independently to a compiled SED, a scattered light image at 0.8 ${\mu}$m and an ALMA dust continuum observation at 890 ${\mu}$m. Results. We find that 90% of the dust mass in HH 48 NE is settled to the disk midplane, less than in average disks, and that the atmospheric layers of the disk contain exclusively large grains (0.3-10 ${\mu}$m). The exclusion of small grains in the upper atmosphere likely has important consequences for the chemistry due to the deep penetration of high-energy photons. The addition of a relatively large cavity (ca. 50 au in radius) is necessary to explain the strong mid-infrared emission, and to fit the scattered light and continuum observations simultaneously., Comment: 16 pages, 8 figures, accepted for publication in Astronomy & Astrophysics

Aims: The JWST program JOYS (JWST Observations of Young protoStars) aims at characterizing the physical and chemical properties of young high- and low-mass star-forming regions, in particular the unique mid-infrared diagnostics of the warmer gas and solid-state components. We present early results from the high-mass star formation region IRAS23385+6053. Methods: The JOYS program uses the MIRI MRS with its IFU to investigate a sample of high- and low-mass star-forming protostellar systems. Results: The 5 to 28mum MIRI spectrum of IRAS23385+6053 shows a plethora of features. While the general spectrum is typical for an embedded protostar, we see many atomic and molecular gas lines boosted by the higher spectral resolution and sensitivity compared to previous space missions. Furthermore, ice and dust absorption features are also present. Here, we focus on the continuum emission, outflow tracers like the H2, [FeII] and [NeII] lines as well as the potential accretion tracer Humphreys alpha HI(7--6). The short-wavelength MIRI data resolve two continuum sources A and B, where mid-infrared source A is associated with the main mm continuum peak. The combination of mid-infrared and mm data reveals a young cluster in its making. Combining the mid-infrared outflow tracer H2, [FeII] and [NeII] with mm SiO data shows a complex interplay of at least three molecular outflows driven by protostars in the forming cluster. Furthermore, the Humphreys alpha line is detected at a 3-4sigma level towards the mid-infrared sources A and B. Following Rigliaco et al. (2015), one can roughly estimate accretion luminosities and corresponding accretion rates between ~2.6x10^-6 and ~0.9x10^-4 M_sun/yr. This is discussed in the context of the observed outflow rates. Conclusions: The analysis of the MIRI MRS observations for this young high-mass star-forming region reveals connected outflow and accretion signatures., Comment: 12 pages, 9 figures, accepted for Astronomy & Astrophysics, the paper is also available at https://www2.mpia-hd.mpg.de/homes/beuther/papers.html

Methanol is the most complex molecule securely identified in interstellar ices and is a key chemical species for understanding chemical complexity in astrophysical environments. Important aspects of the methanol ice photochemistry are still unclear such as the branching ratios and photo-dissociation cross-sections at different temperatures and irradiation fluxes. This work aims at a quantitative agreement between laboratory experiments and astrochemical modelling of the CH3OH ice UV photolysis. This work also allows us to better understand which processes govern the methanol ice photochemistry present in laboratory experiments. We use ProDiMo to simulate the conditions of laboratory measurements. The simulations start with simple chemistry consisting only of methanol ice and helium to mimic the residual gas in the experimental chamber. A surface chemical network enlarged by photo-dissociation reactions is used to study the chemical reactions within the ice. Additionally, different surface chemistry parameters (surface competition, tunnelling, thermal diffusion and reactive desorption) are adopted to check those that reproduce the experimental results. The chemical models with ProDiMo can reproduce the methanol ice destruction via UV photodissociation at temperatures of 20, 30, 50 and 70 K as observed in the experiments. We note that the results are sensitive to different branching ratios after photolysis and to the mechanisms of reactive desorption. In the simulations of a molecular cloud at 20 K, we observed an increase in the methanol gas abundance of one order of magnitude, with a similar decrease in the solid-phase abundance. Comprehensive astrochemical models provide new insights into laboratory experiments as the quantitative understanding of the processes that govern the reactions within the ice. Ultimately, those insights can help to better interpret astronomical observations., Comment: Accepted for publication in A&A on 08-Feb-2023

Icy grain mantles are the main reservoir of the volatile elements that link chemical processes in dark, interstellar clouds with the formation of planets and composition of their atmospheres. The initial ice composition is set in the cold, dense parts of molecular clouds, prior to the onset of star formation. With the exquisite sensitivity of JWST, this critical stage of ice evolution is now accessible for detailed study. Here we show the first results of the Early Release Science program "Ice Age" that reveal the rich composition of these dense cloud ices. Weak ices, including, $^{13}$CO$_2$, OCN$^-$, $^{13}$CO, OCS, and COMs functional groups are now detected along two pre-stellar lines of sight. The $^{12}$CO$_2$ ice profile indicates modest growth of the icy grains. Column densities of the major and minor ice species indicate that ices contribute between 2 and 19% of the bulk budgets of the key C, O, N, and S elements. Our results suggest that the formation of simple and complex molecules could begin early in a water-ice rich environment., Comment: To appear in Nature Astronomy on January 23rd, 2023. 33 pages, 16 figures, 3 tables; includes extended and supplemental data sections. Part of the JWST Ice Age Early Release Science program's science enabling products. Enhanced spectra downloadable on Zenodo at the following DOI: 10.5281/zenodo.7501239

We present an analysis of multi-wavelength observations of the dark globule DC\,314.8--5.1, using data from the Gaia optical, 2MASS near-infrared, and WISE mid-infrared surveys, dedicated imaging with the Spitzer Space Telescope, and X-ray data obtained with the Swift-XRT Telescope (XRT). The main goal was to identify possible pre-main sequence stars (PMSs) and young stellar objects (YSOs) associated with the globule. For this, we studied the infrared colors of all point sources within the boundaries of the cloud. After removing sources with non-stellar spectra, we investigated the Gaia parallaxes for the YSO candidates, and found that none are physically related to DC\,314.8--5.1. In addition, we searched for X-ray emission from pre-main sequence stars with Swift-XRT, and found no 0.5--10\,keV emission down to a luminosity level $\lesssim 10^{31}$erg\,s$^{-1}$, typical of a PMS with mass\,$\ge 2 M_\odot$. Our detailed inspection therefore supports a very young, ``pre-stellar core'' evolutionary stage for the cloud. Based on archival Planck and IRAS data, we moreover identify the presence of hot dust, with temperatures $\gtrsim 100$\,K, in addition to the dominant dust component at 14\,K, originating with the associated reflection nebula., Comment: Accepted to AJ

High quality vibrational spectra of solid-phase molecules in ice mixtures and for temperatures of astrophysical relevance are needed to interpret infrared observations toward protostars and background stars. Over the last 25 years, the Laboratory for Astrophysics at Leiden Observatory has provided more than 1100 spectra of diverse ice samples. Timely with the recent launch of the James Webb Space Telescope, we have fully upgraded the Leiden Ice Database for Astrochemistry (LIDA) adding recently measured spectra. The goal of this manuscript is to describe what options exist to get access to and work with a large collection of IR spectra, and the UV/vis to mid-infrared refractive index of H2O ice and astronomy-oriented online tools to support the interpretation of IR ice observations. LIDA uses Flask and Bokeh for generating the web pages and graph visualization, respectively, SQL for searching ice analogues within the database and Jmol for 3D molecule visualization. The infrared data in the database are recorded via transmission spectroscopy of ice films condensed on cryogenic substrates. The real UV/vis refractive indices of H2O ice are derived from interference fringes created from the simultaneous use of a monochromatic HeNe laser beam and a broadband Xe-arc lamp, whereas the real and imaginary mid-IR values are theoretically calculated. LIDA also offers online tools. The first tool, SPECFY, used to create a synthetic spectrum of ices towards protostars. The second tool aims at the calculation of mid-infrared refractive index values. LIDA allows to search, download and visualize experimental data of astrophysically relevant molecules in the solid phase, as well as to provide the means to support astronomical observations. As an example, we analyse the spectrum of the protostar AFGL 989 using the resources available in LIDA and derive the column densities of H2O, CO and CO2 ices., Comment: Accepted for publication in A&A. 29 pages, 9 figures. Online database at https://icedb.strw.leidenuniv.nl/

We present an analysis of the mid-infrared spectra, obtained from the Spitzer Space Telescope, of the dark globule, DC 314.8--5.1, which is at the onset of low-mass star formation. The target has a serendipitous association with a B-type field star, which illuminates a reflection nebula in the cloud. We focus on the polycyclic aromatic hydrocarbon (PAH) emission features prevalent throughout the mid-infrared range. The analysis of the spectra with the PAHFIT software as well as pypahdb package, shows that (i) the intensities of PAH features decrease over distance from the ionizing star toward the cloud center, some however showing a saturation at larger distances; (ii) the relative intensities of the 6.2 and 8.6 features with respect to the 11.2 micron feature remain high throughout the globule, suggesting a larger cation-to-neutral PAH ratio of the order of unity; the breakdown from pypahdb confirms a high ionized fraction within the cloud; (iii) the pypahdb results display a decrease in large PAH fraction with increased distance from HD 130079, as well as a statistically significant correlation between the large size fraction and the ionized fraction across the globule; (iv) the 7.7 PAH feature displays a peak nearer to 7.8 microns, suggesting a chemically processed PAH population with a small fraction of UV-processed PAHs; (v) the H2 S(0) line is detected at larger distances from the ionizing star. All in all, our results suggest divergent physical conditions within the quiescent cloud DC 314.8--5.1 as compared to molecular clouds with ongoing starformation., Comment: Accepted for publication by the Astrophysical Journal

In this era of spatially resolved observations of planet forming disks with ALMA and large ground-based telescopes such as the VLT, Keck and Subaru, we still lack statistically relevant information on the quantity and composition of the material that is building the planets, such as the total disk gas mass, the ice content of dust, and the state of water in planetesimals. SPICA is an infrared space mission concept developed jointly by JAXA and ESA to address these questions. The key unique capabilities of SPICA that enable this research are (1) the wide spectral coverage 10-220 micron, (2) the high line detection sensitivity of (1-2) 10-19 W m-2 with R~2000-5000 in the far-IR (SAFARI) and 10-20 W m-2 with R~29000 in the mid-IR (SMI, spectrally resolving line profiles), (3) the high far-IR continuum sensitivity of 0.45 mJy (SAFARI), and (4) the observing efficiency for point source surveys. This paper details how mid- to far-IR infrared spectra will be unique in measuring the gas masses and water/ice content of disks and how these quantities evolve during the planet forming period. These observations will clarify the crucial transition when disks exhaust their primordial gas and further planet formation requires secondary gas produced from planetesimals. The high spectral resolution mid-IR is also unique for determining the location of the snowline dividing the rocky and icy mass reservoirs within the disk and how the divide evolves during the build-up of planetary systems. Infrared spectroscopy (mid- to far-IR) of key solid state bands is crucial for assessing whether extensive radial mixing, which is part of our Solar System history, is a general process occurring in most planetary systems and whether extrasolar planetesimals are similar to our Solar System comets/asteroids. ... (abbreviated), Comment: accepted for publication in PASA

Dust grains play an important role in the synthesis of molecules in the interstellar medium, from the simplest species to complex organic molecules. How some of these solid-state molecules are converted into gas-phase species is still a matter of debate. Our aim is to directly compare ice and gas abundances of methanol (CH$_3$OH) and CO, and to investigate the relationship between ice and gas in low-mass protostellar envelopes. We present Submillimeter Array and Atacama Pathfinder EXperiment observations of gas-phase CH$_3$OH and CO towards the multiple protostellar system IRAS05417+0907 located in the B35A cloud. We use archival AKARI ice data toward the same target to calculate CH$_3$OH and CO gas-to-ice ratios. The CO isotopologues emissions are extended, whereas the CH$_3$OH emission is compact and traces the giant outflow emanating from IRAS05417+0907. A discrepancy between submillimeter dust emission and H$_2$O ice column density is found for B35A$-$4 and B35A$-$5, similar to what has previously been reported. B35A$-$2 and B35A$-$3 are located where the submillimeter dust emission peaks and show H$_2$O column densities lower than for B35A$-$4. The difference between the submillimeter continuum emission and the infrared H$_2$O ice observations suggests that the distributions of dust and H$_2$O ice differ around the young stellar objects in this dense cloud. The reason for this may be that the sources are located in different environments resolved by the interferometric observations: B35A$-$2, B35A$-$3 and in particular B35A$-$5 are situated in a shocked region plausibly affected by sputtering and heating impacting the submillimeter dust emission pattern, while B35A$-$4 is situated in a more quiescent part of the cloud. Gas and ice maps are essential to connect small-scale variations in the ice composition with large-scale astrophysical phenomena probed by gas observations., Comment: 18 pages, 15 figures, accepted for publication in Astronomy & Astrophysics

The interaction between dust, ice, and gas during the formation of stars produces complex organic molecules. While observations indicate that several species are formed on ice-covered dust grains and are released into the gas phase, the exact chemical interplay between solid and gas phases and their relative importance remain unclear. Our goal is to study the interplay in regions of low-mass star formation through ice- and gas-mapping and by directly measuring gas-to-ice ratios. This provides constraints on the routes that lead to the chemical complexity that is observed in both phases. We present observations of gas-phase methanol (CH$_3$OH) and carbon monoxide at 1.3 mm towards ten low-mass young protostars in the Serpens SVS4 cluster from the SubMillimeter Array and the Atacama Pathfinder EXperiment telescope. We used archival data from the Very Large Telescope to derive abundances of ice H$_2$O, CO, and CH$_3$OH towards the same region. Finally, we constructed gas-ice maps of SVS4 and directly measured CO and CH$_3$OH gas-to-ice ratios. The CH$_3$OH gas-to-ice ratio agrees with values that were previously reported for embedded Class 0/I low-mass protostars. The CO gas-maps trace an extended gaseous component that is not sensitive to the effect of freeze-out. We find that there is no straightforward correlation between CO and CH$_3$OH gas with their ice counterparts in the cluster. This is likely related to the complex morphology of SVS4: the Class 0 protostar SMM4 and its envelope lie in the vicinity, and the outflow associated with SMM4 intersects the cluster. This study serves as a pathfinder for future observations with ALMA and the James Webb Space Telescope that will provide high-sensitivity gas-ice maps of molecules more complex than methanol. Such comparative maps will be essential to constrain the chemical routes that regulate the chemical complexity in star-forming regions., Comment: 25 pages, 20 figures, accepted for publication by A&A

The physical evolution of Young Stellar Objects (YSOs) is accompanied by an enrichment of the molecular complexity, mainly triggered by the heating and energetic processing of the astrophysical ices. In this paper, a study of how the ice column density varies across the protostellar evolution has been performed. Tabulated data of H$_2$O, CO$_2$, CH$_3$OH, HCOOH observed by ground- and space-based telescopes toward 27 early-stage YSOs were taken from the literature. The observational data shows that ice column density and spectral index ($\alpha$), used to classify the evolutionary stage, are well correlated. A 2D continuum radiative transfer simulation containing bare and grains covered by ices at different levels of cosmic-ray processing were used to calculate the Spectral Energy Distributions (SEDs) in different angle inclinations between face-on and edge-on configuration. The H$_2$O:CO$_2$ ice mixture was used to address the H$_2$O and CO$_2$ column density variation whereas the CH$_3$OH and HCOOH are a byproduct of the virgin ice after the energetic processing. The simulated spectra were used to calculate the ice column densities of YSOs in an evolutionary sequence. As a result, the models show that the ice column density variation of HCOOH with $\alpha$ can be justified by the envelope dissipation and ice energetic processing. On the other hand, the ice column densities are mostly overestimated in the cases of H$_2$O, CO$_2$ and CH$_3$OH, even though the physical and cosmic-ray processing effects are taken into account., Comment: Accepted for publication in The Astrophysical Journal. 29 pages, 16 figures

In dense and cold regions of the interstellar medium (ISM), molecules may be adsorbed onto dust grains to form the ice mantles. Once formed, they can be processed by ionizing radiation coming from stellar or interstellar medium leading to formation of several new molecules in the ice. Among the different kind of ionizing radiation, cosmic rays play an important role in the solid-phase chemistry because of the large amount of energy deposited in the ices. The physicochemical changes induced by the energetic processing of astrophysical ices are recorded in a intrinsic parameter of the matter called complex refractive index (CRI). In this paper, we present for the first time a catalogue containing 39 complex refractive indices (n, k) in the infrared from 2.0 - 16.6 micrometer for 13 different water-containing ices processed in laboratory by cosmic ray analogs. The calculation was done by using the NKABS (acronym of determination of N and K from ABSorbance data) code, which employs the Lambert-Beer and Kramers-Kronig equations to calculate the values of n and k. The results are also available at the website: http://www1.univap.br/gaa/nkabs-database/data.htm. As test case, a H2O:NH3:CO2:CH4 ice was employed in a radiative transfer simulation of a prototoplanetary disk to show that these data are indispensable to reproduce the spectrum of YSOs containing ices.

We present the study of the effects of high energy cosmic rays (CRs) over the astrophysical ices, observed toward the embedded class I protostar Elias 29, by using computational modeling and laboratory data. Its spectrum was observed with {\it Infrared Space Observatory - ISO}, covering 2.3 - 190 $\mu$m. The modeling employed the three-dimensional Monte Carlo radiative transfer code RADMC-3D (Dullemond et al. 2012) and laboratory data of bombarded ice grains by CRs analogs, and unprocessed ices (not bombarded). We are assuming that Elias 29 has a self-irradiated disk with inclination $i =$ 60$^{\circ}$, surrounded by an envelope with bipolar cavity. The results show that absorption features toward Elias 29, are better reproduced by assuming a combination between unprocessed astrophysical ices at low temperature (H$_2$O, CO, CO$_2$) and bombarded ices (H$_2$O:CO$_2$) by high energy CRs. Evidences of the ice processing around Elias 29 can be observed by the good fitting around 5.5-8.0 $\mu$m, by polar and apolar ice segregation in 15.15-15.25 $\mu$m, and by presence of the CH$_4$ and HCOOH ices. Given that non-nitrogen compounds were employed in this work, we assume that absorption around 5.5-8.0 $\mu$m should not be associated with NH$_4^+$ ion (Shutte & Khanna2003), but more probably with aliphatic ethers (e.g. R1-OCH$_2$-R2), CH$_3$CHO and related species. The results obtained in this paper are important, because they show that the environment around protostars is better modeled considering processed samples and, consequently, demonstrates the chemical evolution of the astrophysical ices., Comment: 51 pages, 16 figures, 4 tables