Search

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.

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

We present a study connecting the physical properties of protostellar envelopes to the morphology of the envelope-scale magnetic field. We use the ALMA polarization observations of 55 young prtostars at 0.87 mm on $\sim400-3000$ au scales from the {\em B}-field Orion Protostellar Survey (BOPS) to infer the envelope-scale magnetic field and both dust and gas emission on comparable scales to measure the envelope properties. We find that the protostellar envelopes with compact polarized dust emission tend to have lower envelope masses, than the sources with more extended envelopes. We also find that protostars showing hourglass-field morphologies tend to have lower velocity dispersions in their envelopes, whereas systems with spiral-field morphologies have higher velocity dispersion. Combining with the disk properties taken from the Orion VLA/ALMA Nascent Disk and Multiplicity (VANDAM) survey, we connect envelope properties to fragmentation. Our results suggest that envelope mass may not correlate with fragmentation, whereas turbulence appears to promote fragmentation. On the other hand, we find that fragmentation is suppressed in systems with pinched magnetic fields, suggesting that the magnetic field play a role on providing additional support against gravitational collapse, and the formation of an hourglass-like field may coincide with enhanced magnetic braking that removes angular momentum and hinders the formation of embedded disks. Nevertheless, significant misalignment between magnetic field and outflow axes tends to reduce magnetic braking, leading to the formation of larger disks., Comment: 18 pages, 6 figures, submitted to ApJ

Iron inclusions embedded inside dust grains play a crucial role in both internal alignment via Barnett relaxation and external alignment via the MAgnetically Enhanced RAdiative Torque (MRAT) mechanism. Moreover, inelastic relaxation is predicted to dominate over Barnett relaxation in driving the internal alignment of micron-sized and very large grains above $10\mu m$ (VLGs). Yet, a detailed modeling of polarized thermal dust emission from Class 0/I Young Stellar Objects (YSOs) taking into account these effects and their observational constraints are still lacking. In this paper, we perform synthetic dust polarization modeling for MHD simulations of an intermediate-mass YSOs using our updated POLARIS code developed by Giang et al. (2022). Synthetic polarization results are then post-processed with CASA to confront ALMA observations. We found that to reproduce the high polarization degree of $p \sim 5-30\%$ observed in protostellar envelopes by ALMA, micron-sized and VLGs must contain iron inclusions with $N_{\rm cl} \sim 5 - 10^{4}$ iron atoms per cluster, assuming $30\%$ of iron abundance locked inside dust grains under the cluster form. Inside the inner $\sim 500$ au region, it requires inelastic relaxation, larger iron inclusions of $N_{\rm cl} \sim 10^{2}-10^{4}$, and grain growth beyond $\geq 10\mu m$ to reproduce $\sim 3-10\%$ of polarization observed by ALMA. But given such a combination, the internal alignment and MRAT efficiency acting on VLGs still decrease toward the center, which produces the decrease of $p(\%)$ with increasing gas density, reaching $p \sim 1\%$ inside the disk. Additionally, our results show that the high polarization by ALMA toward YSOs is not totally an artifact caused by the interferometric filtering, but it could be explained by the efficient magnetic alignment of SPM grains by MRAT mechanism., Comment: 26 pages, submitted to ApJ

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

We present new VLA observations, between 6.8mm and 66mm, of the edge-on Class~I disk IRAS04302+2247. Observations at 6.8mm and 9.2mm lead to the detection of thermal emission from the disk, while shallow observations at the other wavelengths are used to correct for emission from other processes. The disk radial brightness profile transitions from broadly extended in previous ALMA 0.9mm and 2.1mm observations to much more centrally brightened at 6.8mm and 9.2mm, which can be explained by optical depth effects. The radiative transfer modeling of the 0.9mm, 2.1mm, and 9.2mm data suggests that the grains are smaller than 1cm in the outer regions of the disk and allows us to obtain the first lower limit for the scale height of grains emitting at millimeter wavelengths in a protoplanetary disk. We find that the millimeter dust scale height is between 1au and 6au at a radius 100au from the central star, while the gas scale height is estimated to be about 7au, indicating a modest level of settling. The estimated dust height is intermediate between less evolved Class 0 sources, that are found to be vertically thick, and more evolved Class II sources, which show a significant level of settling. This suggests that we are witnessing an intermediate stage of dust settling., Comment: Accepted for publication in the Astrophysical Journal

Polarized dust emission is a key tracer in the study of interstellar medium and of star formation. The observed polarization, however, is a product of magnetic field structure, dust grain properties and grain alignment efficiency, as well as their variations in the line of sight, making it difficult to interpret polarization unambiguously. The comparison of polarimetry at multiple wavelengths is a possible way of mitigating this problem. We use data from HAWC+/SOFIA and from SCUBA-2/POL-2 (from the BISTRO survey) to analyse the NGC 2071 molecular cloud at 154, 214 and 850 $\mu$m. The polarization angle changes significantly with wavelength over part of NGC 2071, suggesting a change in magnetic field morphology on the line of sight as each wavelength best traces different dust populations. Other possible explanations are the existence of more than one polarization mechanism in the cloud or scattering from very large grains. The observed change of polarization fraction with wavelength, and the 214-to-154 $\mu$m polarization ratio in particular, are difficult to reproduce with current dust models under the assumption of uniform alignment efficiency. We also show that the standard procedure of using monochromatic intensity as a proxy for column density may produce spurious results at HAWC+ wavelengths. Using both long-wavelength (POL-2, 850 $\mu$m) and short-wavelength (HAWC+, $\lesssim 200\, \mu$m) polarimetry is key in obtaining these results. This study clearly shows the importance of multi-wavelength polarimetry at submillimeter bands to understand the dust properties of molecular clouds and the relationship between magnetic field and star formation., Comment: Main article: 18 pages, 11 figures. Online supplemental material: 2 pages, 3 figures

Recent observational progress has challenged the dust grain-alignment theories used to explain the polarized dust emission routinely observed in star-forming cores. In an effort to improve our understanding of the dust grain alignment mechanism(s), we have gathered a dozen ALMA maps of (sub)millimeter-wavelength polarized dust emission from Class 0 protostars, and carried out a comprehensive statistical analysis of dust polarization quantities. We analyze the statistical properties of the polarization fraction P_frac and dispersion of polarization position angles S. More specifically, we investigate the relationship between S and P_frac as well as the evolution of the product S*P_frac as a function of the column density of the gas in the protostellar envelopes. We find a significant correlation in the polarized dust emission from protostellar envelopes seen with ALMA; the power-law index differs significantly from the one observed by Planck in star-forming clouds. The product S*P_frac, which is sensitive to the dust grain alignment efficiency, is approximately constant across three orders of magnitude in envelope column density. This suggests that the grain alignment mechanism producing the bulk of the polarized dust emission in star-forming cores may not depend systematically on the local conditions such as local gas density. Ultimately, our results suggest dust alignment mechanism(s) are efficient at producing dust polarized emission in the various local conditions typical of Class 0 protostars. The grain alignment efficiency found in these objects seems to be higher than the efficiency produced by the standard RAT alignment of paramagnetic grains. Further study will be needed to understand how more efficient grain alignment via, e.g., different irradiation conditions, dust grain characteristics, or additional grain alignment mechanisms can reproduce the observations.

It has been proposed that the magnetic field, pervasive in the ISM, plays an important role in the process of massive star formation. To better understand its impact at the pre and protostellar stages, high-angular resolution observations of polarized dust emission toward a large sample of massive dense cores are needed. To this end, we used the Atacama Large Millimeter Array in Band 6 (1.3 mm) in full polarization mode to map the polarized emission from dust grains at a physical scale of $\sim$2700 au in the massive protocluster W43-MM1. We used these data to measure the orientation of the magnetic field at the core scale. Then, we examined the relative orientations of the core-scale magnetic field, of the protostellar outflows determined from CO molecular line emission, and of the major axis of the dense cores determined from 2D Gaussian fit in the continuum emission. We found that the orientation of the dense cores is not random with respect to the magnetic field. Instead, the dense cores are compatible with being oriented 20-50$^\deg$ with respect to the magnetic field. The outflows could be oriented 50-70$^\deg$ with respect to the magnetic field, or randomly oriented with respect to the magnetic field, similar to current results in low-mass star-forming regions. In conclusion, the observed alignment of the position angle of the cores with respect to the magnetic field lines shows that the magnetic field is well coupled with the dense material; however, the 20-50$^\deg$ preferential orientation contradicts the predictions of the magnetically-controlled core-collapse models. The potential correlation of the outflow directions with respect to the magnetic field suggests that, in some cases, the magnetic field is strong enough to control the angular momentum distribution from the core scale down to the inner part of the circumstellar disks where outflows are triggered., Comment: Accepted by A&A, 8 pages, 4 figures