Search

This work aims to utilize the increased sensitivity and resolution of the JWST to quantify the HDO/H$_{2}$O ratio in ices toward young stellar objects (YSOs) and to determine if the HDO/H$_{2}$O ratios measured in the gas phase toward massive YSOs (MYSOs) are representative of the ratios in their ice envelopes. Two protostars observed in the Investigating Protostellar Accretion (IPA) program using JWST NIRSpec were analyzed: HOPS 370, an intermediate-mass YSO (IMYSO), and IRAS 20126+4104, a MYSO. The HDO ice toward these sources was detected above the 3$\sigma$ level and quantified via its 4.1 $\mu$m band. The contributions from the CH$_{3}$OH combination modes to the observed optical depth in this spectral region were constrained via the CH$_{3}$OH 3.53 $\mu$m band to ensure that the integrated optical depth of the HDO feature was not overestimated. H$_{2}$O ice was quantified via its 3 $\mu$m band. From these fits, ice HDO/H$_{2}$O abundance ratios of 4.6$\pm$1.8$\times$10$^{-3}$ and 2.6$\pm$1.2$\times$10$^{-3}$ are obtained for HOPS 370 and IRAS 20126+4104, respectively. The simultaneous detections of both crystalline HDO and crystalline H$_{2}$O corroborate the assignment of the observed feature at 4.1 $\mu$m to HDO ice. The ice HDO/H$_{2}$O ratios are similar to the highest reported gas HDO/H$_{2}$O ratios measured toward MYSOs as well as the hot inner regions of isolated low-mass protostars, suggesting that at least some of the gas HDO/H$_{2}$O ratios measured toward massive hot cores are representative of the HDO/H$_{2}$O ratios in ices. The need for an H$_{2}$O-rich CH$_{3}$OH component in the CH$_{3}$OH ice analysis supports recent experimental and observational results that indicate that some CH$_{3}$OH ice may form prior to the CO freeze-out stage in H$_{2}$O-rich ice layers., Comment: Accepted for publication in A&A. 23 pages, 17 figures, 10 tables

Using the MIRI/MRS spectrometer on JWST, we have detected pure rotational, suprathermal OH emissions from the vicinity of the intermediate-mass protostar HOPS 370 (OMC2/FIR3). These emissions are observed from shocked knots in a jet/outflow, and originate in states of rotational quantum number as high as 46 that possess excitation energies as large as $E_U/k = 4.65 \times 10^4$ K. The relative strengths of the observed OH lines provide a powerful diagnostic of the ultraviolet radiation field in a heavily-extinguished region ($A_V \sim 10 - 20$) where direct UV observations are impossible. To high precision, the OH line strengths are consistent with a picture in which the suprathermal OH states are populated following the photodissociation of water in its $\tilde B - X$ band by ultraviolet radiation produced by fast ($\sim 80\,\rm km\,s^{-1}$) shocks along the jet. The observed dominance of emission from symmetric ($A^\prime$) OH states over that from antisymmetric ($A^{\prime\prime}$) states provides a distinctive signature of this particular population mechanism. Moreover, the variation of intensity with rotational quantum number suggests specifically that Ly$\alpha$ radiation is responsible for the photodissociation of water, an alternative model with photodissociation by a 10$^4$ K blackbody being disfavored at a high level of significance. Using measurements of the Br$\alpha$ flux to estimate the Ly$\alpha$ production rate, we find that $\sim 4\%$ of the Ly$\alpha$ photons are absorbed by water. Combined with direct measurements of water emissions in the $\nu_2 = 1 -0$ band, the OH observations promise to provide key constraints on future models for the diffusion of Ly$\alpha$ photons in the vicinity of a shock front., Comment: 30 pages, 7 figures. Accepted for publication in ApJ Letters

The structure and composition of simple ices can be modified during stellar evolution by protostellar heating. Key to understanding the involved processes are thermal and chemical tracers that can diagnose the history and environment of the ice. The 15.2 $\mu$m bending mode of $^{12}$CO$_2$ has proven to be a valuable tracer of ice heating events but suffers from grain shape and size effects. A viable alternative tracer is the weaker $^{13}$CO$_2$ isotopologue band at 4.39 $\mu$m which has now become accessible at high S/N with the $\textit{James Webb}$ Space Telescope (JWST). We present JWST NIRSpec observations of $^{13}$CO$_2$ ice in five deeply embedded Class 0 sources spanning a wide range in luminosities (0.2 - 10$^4$ L$_{\odot}$ ) taken as part of the Investigating Protostellar Accretion Across the Mass Spectrum (IPA) program. The band profiles vary significantly, with the most luminous sources showing a distinct narrow peak at 4.38 $\mu$m. We first apply a phenomenological approach and show that a minimum of 3-4 Gaussian profiles are needed to fit the $^{13}$CO$_2$ absorption feature. We then combine these findings with laboratory data and show that a 15.2 $\mu$m $^{12}$CO$_2$ band inspired five-component decomposition can be applied for the isotopologue band where each component is representative of CO$_2$ ice in a specific molecular environment. The final solution consists of cold mixtures of CO$_2$ with CH$_3$OH, H$_2$O and CO as well as segregated heated pure CO$_2$ ice. Our results are in agreement with previous studies of the $^{12}$CO$_2$ ice band, further confirming that $^{13}$CO$_{2}$ is a useful alternative tracer of protostellar heating events. We also propose an alternative solution consisting only of heated CO$_2$:CH$_3$OH and CO$_2$:H$_2$O ices and warm pure CO$_2$ ice for decomposing the ice profiles of the two most luminous sources in our sample.

We investigate the bright CO fundamental emission in the central regions of five protostars in their primary mass assembly phase using new observations from JWST's Near-Infrared Spectrograph (NIRSpec) and Mid-Infrared Instrument (MIRI). CO line emission images and fluxes are extracted for a forest of $\sim$150 ro-vibrational transitions from two vibrational bands, $v=1-0$ and $v=2-1$. However, ${}^{13}$CO is undetected, indicating that ${}^{12}$CO emission is optically thin. We use H$_2$ emission lines to correct fluxes for extinction and then construct rotation diagrams for the CO lines with the highest spectral resolution and sensitivity to estimate rotational temperatures and numbers of CO molecules. Two distinct rotational temperature components are required for $v=1$ ($\sim600$ to 1000 K and 2000 to $\sim 10^4$ K), while one hotter component is required for $v=2$ ($\gtrsim 3500$ K). ${}^{13}$CO is depleted compared to the abundances found in the ISM, indicating selective UV photodissociation of ${}^{13}$CO; therefore, UV radiative pumping may explain the higher rotational temperatures in $v=2$. The average vibrational temperature is $\sim 1000$ K for our sources and is similar to the lowest rotational temperature components. Using the measured rotational and vibrational temperatures to infer a total number of CO molecules, we find that the total gas masses range from lower limits of $\sim10^{22}$ g for the lowest mass protostars to $\sim 10^{26}$ g for the highest mass protostars. Our gas mass lower limits are compatible with those in more evolved systems, which suggest the lowest rotational temperature component comes from the inner disk, scattered into our line of sight, but we also cannot exclude the contribution to the CO emission from disk winds for higher mass targets., Comment: 31 pages, 7 figures, 4 tables, received to ApJ December 10 2023, accepted to ApJ August 4 2024

Investigating Protostellar Accretion (IPA) is a JWST Cycle~1 GO program that uses NIRSpec IFU and MIRI MRS to obtain 2.9--28~$\mu$m spectral cubes of young, deeply embedded protostars with luminosities of 0.2 to 10,000~L$_{\odot}$ and central masses of 0.15 to 12~M$_{\odot}$. In this Letter, we report the discovery of a highly collimated atomic jet from the Class~0 protostar IRAS~16253$-$2429, the lowest luminosity source ($L_\mathrm{bol}$ = 0.2 $L_\odot$) in the IPA program. The collimated jet is detected in multiple [Fe~II] lines, [Ne~II], [Ni~II], and H~I lines, but not in molecular emission. The atomic jet has a velocity of about 169~$\pm$~15~km\,s$^{-1}$, after correcting for inclination. The width of the jet increases with distance from the central protostar from 23 to~60 au, corresponding to an opening angle of 2.6~$\pm$~0.5\arcdeg. By comparing the measured flux ratios of various fine structure lines to those predicted by simple shock models, we derive a shock {speed} of 54~km\,s$^{-1}$ and a preshock density of 2.0$\times10^{3}$~cm$^{-3}$ at the base of the jet. {From these quantities and using a suite of jet models and extinction laws we compute a mass loss rate between $0.4 -1.1\times10^{-10}~M_{\odot}$~yr~$^{-1}$.} The low mass loss rate is consistent with simultaneous measurements of low mass accretion rate ($2.4~\pm~0.8~\times~10^{-9}~M_{\odot}$~yr$^{-1}$) for IRAS~16253$-$2429 from JWST observations (Watson et al. in prep), indicating that the protostar is in a quiescent accretion phase. Our results demonstrate that very low-mass protostars can drive highly collimated, atomic jets, even during the quiescent phase., Comment: Accepted to ApJL. Comments and feedback welcome

Investigating Protostellar Accretion is a Cycle 1 JWST program using the NIRSpec+MIRI integral field units to obtain 2.9--28 $\mu$m spectral cubes of five young protostars with luminosities of 0.2-10,000 L$_{\odot}$ in their primary accretion phase. This paper introduces the NIRSpec 2.9--5.3 $\mu$m data of the inner 840-9000 au with spatial resolutions from 28-300 au. The spectra show rising continuum emission; deep ice absorption; emission from H$_{2}$, H~I, and [Fe~II]; and the CO fundamental series in emission and absorption. Maps of the continuum emission show scattered light cavities for all five protostars. In the cavities, collimated jets are detected in [Fe~II] for the four $< 320$~L$_{\odot}$ protostars, two of which are additionally traced in Br-$\alpha$. Knots of [Fe~II] emission are detected toward the most luminous protostar, and knots of [FeII] emission with dynamical times of $< 30$~yrs are found in the jets of the others. While only one jet is traced in H$_2$, knots of H$_2$ and CO are detected in the jets of four protostars. H$_2$ is seen extending through the cavities, showing that they are filled by warm molecular gas. Bright H$_2$ emission is seen along the walls of a single cavity, while in three cavities narrow shells of H$_2$ emission are found, one of which has an [Fe~II] knot at its apex. These data show cavities containing collimated jets traced in atomic/ionic gas surrounded by warm molecular gas in a wide-angle wind and/or gas accelerated by bow shocks in the jets., Comment: 26 pages, 12 figures

We present a two-epoch Hubble Space Telescope (HST) study of NGC 2071 IR highlighting HOPS 361-C, a protostar producing an arced 0.2 parsec-scale jet. Proper motions for the brightest knots decrease from 350 to 100 km/s with increasing distance from the source. The [Fe II] and Pa$\beta$ emission line intensity ratio gives a velocity jump through each knot of 40--50 km/s. A new [O I] 63 \mic\ spectrum, taken with the German REciever for Astronomy at Terahertz frequencies (GREAT) instrument aboard Stratospheric Observatory for Infrared Astronomy (SOFIA), shows a low line-of-sight velocity indicative of high jet inclination. Proper motions and jump velocities then estimate 3D flow speed for knots. Subsequently, we model knot positions and speeds with a precessing jet that decelerates. Measurements are matched with a precession period of 1,000--3,000 years and half opening angle of $15^\circ$. The [Fe II] 1.26-to-1.64 \mic\ line intensity ratio determines visual extinction to each knot from 5--30 mag. Relative to $\sim$14 mag of extinction through the cloud from $\rm{C}^{18}$O emission maps, the jet is embedded at a 1/5 to 4/5 fractional cloud depth. Our model suggests the jet is dissipated over a 0.2 pc arc. This short distance may result from the jet sweeping through a wide angle, allowing the cloud time to fill cavities opened by the jet. Precessing jets contrast with nearly unidirectional protostellar jets that puncture host clouds and can propagate significantly further., Comment: 25 pages, 10 figures, submitted to ApJ on 22-11-2022, revised until 8-3-2023, accepted 9-3-2023

We present an 870 $\mu$m continuum survey of 300 protostars from the Herschel Orion Protostar Survey using the Atacama Compact Array (ACA). These data measure protostellar flux densities on envelope scales $\leq$8000 au (20") and resolve the structure of envelopes with 1600 au (4") resolution, a factor of 3-5 improvement in angular resolution over existing single-dish 870 $\mu$m observations. We compare the ACA observations to Atacama Large Millimeter/submillimeter Array 12 m array observations at 870 $\mu$m with $\sim$0.1 (40 au) resolution. Using the 12 m data to measure the fluxes from disks and the ACA data within 2500 au to measure the combined disk plus envelope fluxes, we calculate the 12 m/ACA 870 $\mu$m flux ratios. Our sample shows a clear evolution in this ratio. Class 0 protostars are mostly envelope-dominated with ratios $<$0.5. In contrast, Flat Spectrum protostars are primarily disk-dominated with ratios near 1, although with a number of face-on protostars dominated by their envelopes. Class I protostars span the range from envelope to disk-dominated. The increase in ratio is accompanied by a decrease in the envelope fluxes and estimated mass infall rates. We estimate that 80$\%$ of the mass is accreted during the envelope-dominated phase. We find that the 12 m/ACA flux ratio is an evolutionary indicator that largely avoids the inclination and foreground extinction dependence of spectral energy distribution-based indicators.

We present a Spitzer/Herschel focused survey of the Aquila molecular clouds ($d \sim 436$~pc) as part of the eHOPS (extension of HOPS Out to 500 ParSecs) census of nearby protostars. For every source detected in the Herschel/PACS bands, the eHOPS-Aquila catalog contains 1-850~$\mu$m SEDs assembled from 2MASS, Spitzer, Herschel, WISE, and JCMT/SCUBA-2 data. Using a newly developed set of criteria, we classify objects by their SEDs as protostars, pre-ms sequence stars with disks, and galaxies. A total of 172 protostars are found in Aquila, tightly concentrated in the molecular filaments that thread the clouds. Of these, 71 (42\%) are Class 0 protostars, 54 (31\%) are Class I protostars, 43 (25\%) are flat-spectrum protostars, and 4 (2\%) are Class II sources. Ten of the Class 0 protostars are young PACS Bright Red Sources similar to those discovered in Orion. We compare the SEDs to a grid of radiative transfer models to constrain the luminosities, envelope densities, and envelope masses of the protostars. A comparison of the eHOPS-Aquila to the HOPS protostars in Orion finds that the protostellar luminosity functions in the two star-forming regions are statistically indistinguishable, the bolometric temperatures/envelope masses of eHOPS-Aquila protostars are shifted to cooler temperatures/higher masses, and the eHOPS-Aquila protostars do not show the decline in luminosity with evolution found in Orion. We briefly discuss whether these differences are due to biases between the samples, diverging star formation histories, or the influence of environment on protostellar evolution., Comment: Accepted in ApJS. NASA IPAC/IRSA Data Release DOI: 10.26131/IRSA553