Protostellar System Research Articles

PRODIGE – envelope to disk with NOEMA

Context. The formation of stars has been subject to extensive studies in the past decades on scales from molecular clouds to proto-planetary disks. It is still not fully understood how the surrounding material in a protostellar system, which often shows asymmetric structures with complex kinematic properties, feeds the central protostar(s) and their disk(s). Aims. We study the spatial morphology and kinematic properties of the molecular gas surrounding the IRS3A and IRS3B protostellar systems in the L1448N region located in the Perseus molecular cloud. Methods. We present 1 mm Northern Extended Millimeter Array (NOEMA) observations of the large program PROtostars & DIsks: Global Evolution (PRODIGE). We analyzed the kinematic properties of molecular lines. Because the spectral profiles are complex, the lines were fit with up to three Gaussian velocity components. The clustering algorithm called density-based spatial clustering of applications with noise (DBSCAN) was used to disentangle the velocity components in the underlying physical structure. Results. We discover an extended gas bridge (≈3000 au) surrounding both the IRS3A and IRS3B systems in six molecular line tracers (C18O, SO, DCN, H2CO, HC3N, and CH3OH). This gas bridge is oriented in the northeast-southwest direction and shows clear velocity gradients on the order of 100 km s−1 pc−1 toward the IRS3A system. We find that the observed velocity profile is consistent with analytical streamline models of gravitational infall toward IRS3A. The high-velocity C18O (2-1) emission toward IRS3A indicates a protostellar mass of ≈1.2 M⊙. Conclusions. While high angular resolution continuum data often show IRS3A and IRS3B in isolation, molecular gas observations reveal that these systems are still embedded within a large-scale mass reservoir, whose spatial morphology and velocity profiles are complex. The kinematic properties of the extended gas bridge are consistent with gravitational infall toward the protostar IRS3A.

Multiple Components of the Outflow in the Protostellar System HH 212: Outer Outflow Shell, Rotating Wind, Shocked Wind, and Jet

We present the Atacama Large Millimeter/submillimeter Array Band 7 observations of the CO (J = 3−2) line emission of the protostellar system HH 212 at ∼24 au spatial resolution and compare them to those of the SiO (J = 8−7) and SO (J = 8−7) line emission reported in the literature. We find that the CO line traces four distinct regions: (1) an outer outflow shell, (2) a rotating wind region between the SiO and CO shells, (3) the shocked and wide-angle inner X-wind inside a SiO shell, and (4) the jet. The origin of the CO outer outflow shell could be associated with the entrained material of the envelope, or an extended disk wind. The rotating wind, which is shocked, is launched from a radius of 9–15 au, slightly exterior to that of the previously detected SO shell, which marks the boundary where the wide-angle X-wind is interacting with and shocking the disk wind. Additionally, the SO is found to be mixed with the CO emission within the thick and extended rotating wind region. The large-scale CO-shocked wind coexists with the SO emission near the upper portion of the inner shocked region converged on top of the inner SiO knots. The CO jet is traced by a chain of knots with roughly equal intervals, exhibiting quasiperiodicity, as reported in other jets in the literature.

JWST Observations of Young protoStars (JOYS). Overview of gaseous molecular emission and absorption in low-mass protostars

The Mid-InfraRed Instrument (MIRI) on board the James Webb Space Telescope (JWST) allows one to probe the molecular gas composition at mid-infrared (mid-IR) wavelengths with 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. Previous studies were sensitive primarily to high-mass protostars. The aim of this paper is to derive the physical conditions of all gas-phase molecules detected toward a sample of 18 low-mass protostars as part of the JWST Observations of Young protoStars (JOYS) program and to determine the origin of the molecular emission and absorption features. This includes molecules such as CO$_2$, C$_2$H$_2$, and CH$_4$ that cannot be studied at millimeter wavelengths. We present JWST/MIRI data taken with the Medium Resolution Spectrometer (MRS) of 18 low-mass protostellar systems, focusing on gas-phase molecular lines in spectra extracted from the central protostellar positions. The column densities and excitation temperatures were derived for each molecule using local thermodynamic equilibrium (LTE) slab models. Ratios of the column densities (absorption) or total number of molecules (emission) were taken with respect to H$_2$O in order to compare these to ratios derived in interstellar ices. Continuum emission is detected across the full MIRI-MRS wavelength toward 16/18 sources; the other two sources (NGC 1333 IRAS 4B and Ser-S68N-S) are too embedded to be detected. The MIRI-MRS spectra show a remarkable richness in molecular features across the full wavelength range, in particular toward B1-c (absorption) and L1448-mm (emission). Besides H$_2$, which is not considered here, water is the most commonly detected molecule (12/16) toward the central continuum positions followed by CO$_2$ (11/16), CO (8/16), and OH (7/16). Other molecules such as 13CO$_2$, C$_2$H$_2$, 13CCH$_2$, HCN, C$_4$H$_2$, CH$_4$, and SO$_2$ are detected only toward at most three of the sources, particularly toward B1-c and L1448-mm. 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 NH$_3$ 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 ($ K). Toward three sources (e.g., TMC1-W), hot ($ K) H$_2$O 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 H$_2$O between ice and gas points 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). Molecular emission and absorption features trace various warm components in young protostellar systems, from the hot core regions to shocks in the outflows and disk winds. 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, unsettled dust grains in upper layers of the disk.

Accurate proper motions of the protostellar binary system L 1551 IRS 5

Abstract We present an extensive astrometric study of the protostellar binary system L 1551 IRS 5, utilising nearly four decades of interferometric observations obtained between 1983 and 2022 with the Karl G. Jansky Very Large Array (VLA) and the Atacama Large Millimeter/submillimeter Array (ALMA). We focus on observations with sufficient angular resolution to separate the two protostars (L 1551 IRS 5 N and S) in the system and derive accurate absolute proper motions for the two sources, as well as the relative proper motion between them. The absolute proper motion is dominated by the solar motion with only a modest contribution from L 1551 IRS 5’s peculiar velocity, as expected for a young stellar object. The relative proper motions enable us to constrain the orbit and derive a total mass of 0.96 ± 0.17 M⊙ for the system. While the emission of both sources at wavelengths shorter than about 1.3 cm is compact, the emission at longer wavelengths (λ ≳ 2 cm) is often affected by a free-free contribution from nearby shock features. The results presented here demonstrate that, when appropriate care is taken to combine the observations, interferometric data collected with different facilities, at different frequencies, and with different gain calibrators can be combined to obtain accurate astrometry. Observations of L 1551 IRS 5 over the next several decades with the VLA, ALMA and eventually the ngVLA and SKA ought to improve its dynamical mass measurement down to an accuracy of a few percent. Similar observations of other young multiple systems have the unique potential to provide dynamical mass estimates for the youngest known stellar objects.

JOYS+: The link between the ice and gas of complex organic molecules

Context. A rich inventory of complex organic molecules (COMs) has been observed in high abundances in the gas phase toward Class 0 protostars. It has been suggested that these molecules are formed in ices and sublimate in the warm inner envelope close to the protostar. However, only the most abundant COM, methanol (CH3OH), had been firmly detected in ices before the era of the James Webb Space Telescope (JWST). Now, it is possible to detect the interstellar ices of other COMs and constrain their ice column densities quantitatively. Aims. 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 (hereafter IRAS 2A) and B1-c, as case studies in our JWST Observations of Young proto-Stars (JOYS+) program. By comparing the column density ratios with respect to CH3OH between both phases measured in the same sources, we can probe the evolution of COMs from ice to gas in the early stages of star formation. Methods. The column densities of COMs in gas and ice were derived by fitting the spectra observed by the Atacama Large Millimeter/submillimeter Array (ALMA) and the JWST/Mid-InfraRed Instrument-Medium Resolution Spectroscopy (MIRI-MRS), respectively. The gas-phase emission lines were fit using local thermal equilibrium models, and the ice absorption bands were fit by matching the infrared spectra measured in laboratories. The column density ratios of four O-COMs (CH3CHO, C2H5OH, CH3OCH3, and CH3OCHO) with respect to CH3OH were compared between ice and gas in IRAS 2A and B1-c. Results. We were able to fit the fingerprint range of COM ices between 6.8 and 8.8 μm in the JWST/MIRI-MRS spectra of B1-c using similar components to the ones recently used for NGC 1333 IRAS 2A. We claim detection of CH4, OCN−, HCOO−, HCOOH, CH3CHO, C2H5OH, CH3OCH3, CH3OCHO, and CH3COCH3 in B1-c, and upper limits have been estimated for SO2, CH3COOH, and CH3CN. The total abundance of O-COM ices is constrained to be 15% with respect to H2O ice, 80% of which is dominated by CH3OH. The comparison of O-COM ratios with respect to CH3OH between ice and gas shows two different cases. On the one hand, 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. On the other hand, the ice ratios of CH3CHO and C2H5OH with respect to CH3OH are higher than the gas ratios by 1–2 orders of magnitude. This difference can be explained by gas-phase reprocessing following sublimation, or different spatial distributions of COMs in the envelope, which is an observational effect resulting from ALMA and JWST tracing different components in a protostellar system. Conclusions. The firm detection of COM ices other than CH3OH is reported in another well-studied low-mass protostar, B1-c, following the recent detection in NGC 1333 IRAS 2A. The column density ratios of four O-COMs with respect to CH3OH show both similarities and differences between gas and ice. Although the straightforward explanations would be the direct inheritance from ice to gas and the gas-phase reprocessing, respectively, other possibilities such as different spatial distributions of molecules cannot be excluded.

Early Planet Formation in Embedded Disks. XI. A High-resolution View Toward the BHR 71 Class 0 Protostellar Wide Binary

We present Atacama Large Millimeter/submillimeter Array (ALMA) observations of the binary Class 0 protostellar system BHR 71 IRS1 and IRS2 as part of the Early Planet Formation in Embedded Disks (eDisk) ALMA Large Program. We describe the 12CO (J = 2–1), 13CO (J = 2–1), C18O (J = 2–1), H2CO (J = 32,1–22,0), and SiO (J = 5–4) molecular lines along with the 1.3 mm continuum at high spatial resolution (∼0.″08 or ∼5 au). Dust continuum emission is detected toward BHR 71 IRS1 and IRS2, with a central compact component and extended continuum emission. The compact components are smooth and show no sign of substructures such as spirals, rings, or gaps. However, there is a brightness asymmetry along the minor axis of the presumed disk in IRS1, possibly indicative of an inclined geometrically and optically thick disk-like component. Using a position–velocity diagram analysis of the C18O line, clear Keplerian motions were not detected toward either source. If Keplerian rotationally supported disks are present, they are likely deeply embedded in their envelope. However, we can set upper limits of the central protostellar mass of 0.46 M ⊙ and 0.26 M ⊙ for BHR 71 IRS1 and BHR 71 IRS2, respectively. Outflows traced by 12CO and SiO are detected in both sources. The outflows can be divided into two components, a wide-angle outflow and a jet. In IRS1, the jet exhibits a double helical structure, reflecting the removal of angular momentum from the system. In IRS2, the jet is very collimated and shows a chain of knots, suggesting episodic accretion events.

Simulated analogues I: apparent and physical evolution of young binary protostellar systems

ABSTRACT Protostellar binaries harbour complex environment morphologies. Observations represent a snapshot in time, and projection and optical depth effects impair our ability to interpret them. Careful comparison with high-resolution models that include the larger star-forming region can help isolate the driving physical processes and give context in the time domain to the observations. We carry out four zoom-in simulations with au scale resolution that result in three binaries and a single star. For the first time ever, we follow the detailed evolution of a protobinary in a full molecular cloud context until a circumbinary disc forms. We investigate the gas dynamics around the young stars and extract disc sizes. Using radiative transfer, we obtain the evolutionary tracer Tbol of the binary systems. We find that the centrifugal radius in prestellar cores is a poor estimator of the resulting disc size due to angular momentum transport at all scales. For binaries, the disc sizes are regulated periodically by the binary orbit, having larger radii close to the apastron. The bolometric temperature differs systematically between edge-on and face-on views and shows a high-frequency time dependence correlated with the binary orbit and a low-frequency time dependence with larger episodic accretion events. These oscillations can cause the appearance of the system to change rapidly from class 0 to class I and, for short periods, even bring it to class II. The highly complex structure in early stages, as well as the binary orbit itself, affects the classical interpretation of protostellar classes, and the direct translation to evolutionary stages has to be done with caution and include other evolutionary indicators such as the extent of envelope material.

Simulated analogues II: a new methodology for non-parametric matching of models to observations

ABSTRACT Star formation is a multiscale problem, and only global simulations that account for the connection from the molecular cloud-scale gas flow to the accreting protostar can reflect the observed complexity of protostellar systems. Star-forming regions are characterized by supersonic turbulence, and as a result, it is not possible to simultaneously design models that account for the larger environment and in detail reproduce observed stellar systems. Instead, the stellar inventories can be matched statistically, and the best matches found that approximate specific observations. Observationally, a combination of single-dish telescopes and interferometers are now able to resolve the nearest protostellar objects on all scales from the protostellar core to the inner $10\, \mathrm{au}$. We present a new non-parametric methodology which uses high-resolution simulations and post-processing methods to match simulations and observations using deep learning. Our goal is to perform a down-selection from large data sets of synthetic images to a ranked list of best-matching candidates with respect to the observation. This is particularly useful for binary and multiple stellar systems that form in turbulent environments. The objective is to accelerate the rate at which we can do such comparisons, remove biases from hand-picking matches, and contribute to identifying the underlying physical processes that drive the creation and evolution of observed protostellar systems.

Chemical Abundance Gradients of Organic Molecules within a Protostellar Disk

Observations of low-mass protostellar systems show evidence of rich, complex organic chemistry. Their low luminosity, however, makes determining abundance distributions of complex organic molecules within the water snowline challenging. However, the excitation conditions sampled by differing molecular distributions may produce substantive changes in the resulting emission. Thus, molecular excitation may recover spatial information from spatially unresolved data. By analyzing spatially unresolved NOrthern Extended Millimeter Array observations of CH3OH and CH3CN, we aim to determine if CH3OH and CH3CN are distributed differently in the protostellar disk around HOPS-370, a highly luminous intermediate-mass protostar. Rotational diagram analysis of CH3OH and CH3CN yields rotational temperatures of 198 ± 1.2 K and 448 ± 19 K, respectively, suggesting the two molecules have different spatial distributions. Source-specific 3D LTE radiative transfer models are used to constrain the spatial distribution of CH3OH and CH3CN within the disk. A uniform distribution with an abundance of 4 × 10−8 reproduces the CH3OH observations. In contrast, the spatial distribution of CH3CN needs to be either more compact (within ∼120 au versus ∼240 au for CH3OH) or exhibit a factor of ≳15 increase in abundance in the inner ∼55 au. A possible explanation for the difference in spatial abundance distributions of CH3OH and CH3CN is carbon-grain sublimation.

The factors that influence protostellar multiplicity

Context. Protostellar multiplicity is common at all stages and mass ranges. However, the factors that determine the multiplicity of protostellar systems have not been systematically characterized through their molecular gas. Aims. We characterize the physical properties of the Perseus molecular cloud at ≥5000 AU scales by mapping the diagnostic molecular lines. Methods. We used Nobeyama 45m Radio Observatory (NRO) on-the-fly maps of HCN, HNC, HCO+, and N2H+ (J=1–0) toward five subregions in Perseus, complemented with single-pointing Atacama Pathfinder Experiment (APEX) observations of HNC (J = 4–3), to derive the physical parameters of the dense gas. The spatial resolutions of both observations were ~18″, which is equivalent to ~5000 AU scales at the distance of Perseus. The kinetic gas temperature was derived from the I(HCN)/I(HNC) J ratio, and the H2 density was obtained from the HNC J=4–3/J=1–0 ratio. These parameters were used to obtain the N2H+ (cold) and HCO+ (warm) gas masses. The inferred and derived parameters were then compared to source the parameters, including protostellar multiplicity, bolometric luminosity, and dust envelope mass. Results. The inferred mean kinetic gas temperature (I(HCN)/I(HNC) J=1–0 ratio; ranging between 15 and 26 K), and H2 volumetric density (HNC J=4–3/J=1–0; 105−106 cm−3) are not correlated with multiplicity in Perseus. The derived gas and dust masses, 1.3 to 16 × 10−9 M⊙ for the cold-gas mass (N2H+), 0.1 to 25 M⊙ for the envelope dust masses (850 μm), and 0.8 to 10 × 10−10 M⊙ for the warm-gas mass (HCO+), are correlated to multiplicity and to the number of protostellar components. The warm-gas masses are lower by a factor of 16 than the cold-gas masses. Conclusions. The gas and dust mass is correlated to multiplicity at ~5000 AU scales in Perseus. Higher-order multiples tend to have higher gas and dust masses in general, while close binaries (separations ≤7″) and single protostars have similar gas and dust mass distributions. On the other hand, the H2 density and kinetic gas temperature are not correlated with multiplicity.

Correlations among complex organic molecules around protostars: Effects of physical structure

Context. Complex organic molecules have been observed toward many protostars. Their column density ratios are generally constant across protostellar systems, with some low-level scatter. However, the scatter in the column density ratio of formamide (NH2CHO) to methanol (CH3OH), NNH2CHO/NCH3OH, is one of the highest compared to other ratios. The larger scatter for NNH2CHO/NCH3OH (or weak correlation of these two molecules) is sometimes interpreted as evidence of gas-phase formation of NH2CHO. Aims. In this work, we propose an alternative interpretation in which this scatter is produced by differences in the snowline locations related to differences in binding energies of these species (formamide typically has a ≳2000 K larger binding energy than methanol) and the small-scale structure of the envelope and the disk system. Therefore, we do not include chemistry in our models in order to isolate the effect of physical factors. We also include CH3CN in our work as a control molecule, as it has a similar binding energy to CH3OH. Methods. We used radiative transfer models to calculate the emission from NH2CHO, CH3OH, and CH3CN in protostellar systems with and without disks. The abundances of these species were parameterized in our models, and we fit the calculated emission lines to find the column densities and excitation temperatures of these species, as done in real observations. Results. Given the difference in binding energies of NH2CHO and CH3OH, we find the gas-phase NNH2CHO/NCH3OH needs to be multiplied by a correction factor of approximately ten in order to give the true abundance ratio of these two species in the ices. This factor is much smaller (i.e., ~2) for NCH3CN/NCH3OH (the control molecule). We find that models with different disk sizes, luminosities, and envelope masses produce a scatter in this correction factor, and hence in NNH2CHO/NCH3OH comparable with that of observations. The scatter in NNH2CHO/NCH3OH is larger than that of NCH3CN/NCH3OH in models consistent with the observations. However, the scatter in the models for NCH3CN/NCH3OH is smaller than observations by a factor of around two, as expected from the similar binding energies of CH3OH and CH3CN pointing to the need for some chemical effects in the gas or ice to explain the observed ratios. We show that the scatter in NNH2CHO/NCH3OH will be lower than previously measured if we correct for the difference in sublimation temperatures of these two species in observations of ~40 protostellar systems with ALMA. Conclusions. The scatter in NNH2CHO/NCH3OH (or the ratio of any two molecules with a large binding energy difference) can be partially explained by the difference in their binding energies. Correction for this bias makes the scatter in this ratio similar to that in ratios of other complex organics in the observations, making NH2CHO a “normal” molecule. Therefore, we conclude that gas-phase chemistry routes for NH2CHO are not necessary to explain the larger scatter of NNH2CHO/NCH3OH compared with other ratios.

JWST Observations of Young protoStars (JOYS)

Context. Accretion and ejection dictate the outcomes of star and planet formation processes. The mid-infrared (MIR) wavelength range offers key tracers of processes that have been difficult to detect and spatially resolve in protostars until now. Aims. We aim to characterize the interplay between accretion and ejection in the low-mass Class I protobinary system TMC1, comprising two young stellar objects: TMC1-W and TMC1-E at a 85 au separation. Methods. Using the James Webb Space Telescope (JWST) Mid-Infrared Instrument (MIRI) observations in 5–28 μm range, we measured the intensities of emission lines of H2, atoms, and ions, for instance, the [Fe II] and [Ne II], and HI recombination lines. We analyzed the spatial distribution of the different species using the MIRI Medium Resolution Spectrometer (MRS) capabilities to spatially resolve emission at 0″​​.2–0″​​.7 scales. we compared these results with the corresponding Atacama Large Millimeter/submillimeter Array (ALMA) maps tracing cold gas and dust. Results. We detected H2 outflow coming from TMC1-E, with no significant H2 emission from TMC1-W. The H2 emission from TMC1-E outflow appears narrow and extends to wider opening angles with decreasing Eup from S(8) to S(1) rotational transitions, indicating the disk wind as its origin. The outflow from TMC1-E protostar shows spatially extended emission lines of [Ne II], [Ne III], [Ar II], and [Ar III], with their line ratios consistent with UV radiation as a source of ionization. With ALMA, we detected an accretion streamer infalling from > 1000 au scales onto the TMC1-E component. The TMC1-W protostar powers a collimated jet, detected with [Fe II] and [Ni II], making it consistent with energetic flow. A much weaker ionized jet is observed from TMC1-E, and both jets appear strikingly parallel to each other, indicating that the disks are co-planar. TMC1-W is associated with strong emission from hydrogen recombination lines, tracing the accretion onto the young star. Conclusions. MIRI-MRS observations provide an unprecedented view of protostellar accretion and ejection processes on 20 au scales. Observations of a binary Class I protostellar system show that the two processes are clearly intertwined, with accretion from the envelope onto the disk influencing a wide-angle wind ejected on disk scales. Finally, the accretion from the protostellar disk onto the protostar is associated with the source launching a collimated high-velocity jet within the innermost regions of the disk.

The rotation-tunneling spectrum of 3-hydroxypropenal and confirmation of its detection toward IRAS 16293-2422 B

Context. 3-Hydroxypropenal (HOCHCHCHO) is the lower energy tautomer of malonaldehyde which displays a complex rotation-tunneling spectrum. It was detected somewhat tentatively toward the solar-type protostellar system IRAS 16293-2422 with ALMA in the framework of the Protostellar Interferometric Line Survey (PILS). Several transitions, however, had large residuals, preventing not only their detection, but also the excitation temperature of the species from being determined unambiguously. Aims. We want to extend the existing rotational line list of 3-hydroxypropenal to shed more light on the recent observational results and to facilitate additional radio astronomical searches for this molecule. Methods. We recorded and analyzed the rotation-tunneling spectrum of 3-hydroxypropenal in the frequency regions between 150 and 330 GHz and between 400 and 660 GHz. Transitions were searched for in the PILS observations of IRAS 16293-2422. Local thermodynamic equilibrium (LTE) models were carried out and compared to the observations to constrain the excitation temperature. Additional transitions were searched for in other ALMA archival data of the same source to confirm the presence of 3-hydroxypropenal. Results. More than 7500 different spectral lines, corresponding to more than 11 500 transitions, were assigned in the course of our investigation with quantum numbers 2 ≤ J ≤ 100, Ka ≤ 59, and Kc ≤ 97, resulting in a greatly improved set of spectroscopic parameters. The comparison between the LTE models and the observations yields an excitation temperature of 125 K with a column density N = 1.0 × 1015 cm−2 for this species. We identified seven additional lines of 3-hydroxypropenal that show a good agreement with the model in the ALMA archive data. Conclusions. The calculated rotation-tunneling spectrum of 3-hydroxypropenal has sufficient accuracy for radio astronomical searches. With the solution of the excitation temperature conundrum and the detection of seven more lines, we consider the detection of 3-hydroxypropenal toward IRAS 16293-2422 as secure.

Probing the eccentricity in protostellar discs: Modelling kinematics and morphologies

Context. Protostellar discs are mostly modelled as circular structures of gas and dust orbiting a protostar. However, a number of physical mechanisms, for example, the presence of a (sub)stellar companion or initial axial asymmetry, can cause the gas and dust orbital motion to become eccentric. Theoretical studies have revealed that, when present, disc eccentricity is expected to occur with predictable profiles that can be long-lasting and potentially observable in protostellar systems. Aims. We construct an analytical model predicting the typical features of the kinematics and morphology of eccentric protostellar discs, with the final goal of characterising the observational appearance of eccentricity in discs. Methods. We validate the model using a numerical simulation of a circumbinary disc (where the binary makes the disc eccentric). We finally post-process the simulation with Monte Carlo radiative transfer to study how eccentric features would appear through the ‘eyes’ of ALMA. Results. Besides the motion of the material on eccentric Keplerian orbits in the disc orbital plane, the most characteristic eccentric feature emerging from the analytical model is strong vertical motion with a typical anti-symmetric pattern (with respect to the disc line of pericentres). A circumbinary disc with a ≈ 40 au eccentric cavity (ecav = 0.2), carved by an abin = 15 au binary, placed at a distance d = 130 pc, is expected to host in its upper emission surface vertical oscillations up to vz ~ 400 m s−1 close to the cavity edge, that is to say, well within ALMA spectral and spatial resolution capabilities. A residual spiral pattern in the vertical velocity Δvz ~ 150 m s−1 of the simulation cannot be captured by the theoretical model, we speculate it to be possibly linked to the presence of a companion in the system.

Feasibility of detecting shadows in disks induced by infall

Observations performed with high-resolution imaging techniques have revealed the existence of shadows in circumstellar disks that can be explained by the misalignment of an inner disk with respect to an outer one. The cause of misalignment, however, is still a matter of debate. In this study, we investigate the feasibility of observing shadows induced by one prominent scenario that may lead to misalignment, which involves the late infall of material onto a protostellar system. In particular, we used previously performed hydrodynamical simulations of such events and we generated flux maps in the visible, near-infrared, submillimeter, and millimeter wavelength ranges using Monte Carlo radiative transfer. Based on those results, we derived synthetic observations of these systems performed with the instruments SPHERE/VLT and ALMA, which we used as a basis for our subsequent analysis. We find that near-infrared observations with SPHERE are particularly well suited for detecting shadows via direct imaging alongside other features such as gaps, arcs, and streamers. On the contrary, performing a shadow detection based on reconstructed ALMA observations is very challenging due to the high sensitivity that is required for this task. Thus, in cases that allow for a detection, sophisticated analyses may be needed, for instance by the utilization of carefully constructed azimuthal profiles, aiding the search for potentially shallow shadows. Lastly, we conclude that a late infall-induced disk misalignment offers a plausible explanation for the emergence of shadows observed in various systems.

Magnetic Fields Observed along the East–West Outflow of IRAS 16293-2422

Magnetic fields likely play an important role in the formation of young protostars. Multiscale and multiwavelength dust polarization observations can reveal the inferred magnetic field from scales of the cloud to core to protostar. We present continuum polarization observations of the young protostellar triple system IRAS 16293-2422 at 89 μm using HAWC+ on SOFIA. The inferred magnetic field is very uniform with an average field angle of 89° ± 23° (E of N), which is different from the ∼170° field morphology seen at 850 μm at larger scales (≳2000 au) with JCMT POL-2 and at 1.3 mm on smaller scales (≲300 au) with Atacama Large Millimeter/submillimeter Array. The HAWC+ magnetic field direction is aligned with the known E-W outflow. This alignment difference suggests that the shorter wavelength HAWC+ data is tracing the magnetic field associated with warmer dust likely from the outflow cavity, whereas the longer wavelength data are tracing the bulk magnetic field from cooler dust. Also, we show in this source the dust emission peak is strongly affected by the observing wavelength. The dust continuum peaks closer to source B (northern source) at shorter wavelengths and progressively moves toward the southern A source with increasing wavelength (from 22 to 850 μm).

ALMA view of the L1448-mm protostellar system on disk scales: CH3OH and H13CN as new disk wind tracers

Protostellar disks are known to accrete; however, the exact mechanism that extracts the angular momentum and drives accretion in the low-ionization “dead” region of the disk is under debate. In recent years, magnetohydrodynamic (MHD) disk winds have become a popular solution. Even so, observations of these winds require both high spatial resolution (~10 s au) and high sensitivity, which has resulted in only a handful of MHD disk wind candidates to date. In this work we present high angular resolution (~30 au) ALMA observations of the emblematic L1448-mm protostellar system and find suggestive evidence for an MHD disk wind. The disk seen in dust continuum (~0.9 mm) has a radius of ~23 au. Rotating infall signatures in H13CO+ indicate a central mass of 0.4 ± 0.1 M⊙ and a centrifugal radius similar to the dust disk radius. Above the disk, we identify rotation signatures in the outflow traced by H13CN, CH3OH, and SO lines and find a kinematical structure consistent with theoretical predictions for MHD disk winds. This is the first detection of an MHD disk wind candidate in H13CN and CH3OH. The wind launching region estimated from cold MHD wind theory extends out to the disk edge. The magnetic lever arm parameter would be λϕ ≃ 1.7, in line with recent non-ideal MHD disk models. The estimated mass-loss rate is approximately four times the protostellar accretion rate (Ṁacc ≃ 2 × 10−6M⊙ yr−1) and suggests that the rotating wind could carry enough angular momentum to drive disk accretion.

The First Estimation of the Ambipolar Diffusivity Coefficient from Multi-scale Observations of the Class 0/I Protostar, HOPS-370

Protostars are born in magnetized environments. As a consequence, the formation of protostellar disks can be suppressed by the magnetic field, efficiently removing the angular momentum of the infalling material. Nonideal MHD effects are proposed as one way to allow protostellar disks to form. Thus, it is important to understand their contributions to observations of protostellar systems. We derive an analytical equation to estimate the ambipolar diffusivity coefficient at the edge of the protostellar disk in the Class 0/I protostar, HOPS-370, for the first time, under the assumption that the disk radius is set by ambipolar diffusion. Using previous results of the protostellar mass, disk mass, disk radius, density and temperature profiles, and magnetic field strength, we estimate the ambipolar diffusivity coefficient to be 1.7−1.4+1.5×1019cm2s−1 . We quantify the contribution of ambipolar diffusion by estimating its dimensionless Elsässer number to be ∼1.7−1.0+1.0 , indicating its dynamical importance in this region. We compare our results to those of the chemical calculations of the ambipolar diffusivity coefficient using the Non-Ideal Magnetohydrodynamics Coefficients and Ionization Library, which are consistent with our results. In addition, we compare our derived ambipolar diffusivity coefficient to the diffusivity coefficients for ohmic dissipation and the Hall effect, and find ambipolar diffusion is dominant in our density regime. These results demonstrate a new methodology to understand nonideal MHD effects in observations of protostellar disks. More detailed modeling of the magnetic field, envelope, and microphysics, along with a larger sample of protostellar systems, is needed to further understand the contributions of nonideal MHD.

Hunting for complex cyanides in protostellar ices with the JWST

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. The James Webb Space Telescope (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 CH3CN, C2H5CN and the simple molecule, N2O, 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 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 CH3CN/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 ~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 of the enhancement of OCN−, CH3CN, and C2H5CN in warmer ices; although, the 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 C2H5CN ice, are necessary for robust detection and conclusions on the formation history of complex cyanides.

HCN as a probe of the inner disc in a candidate proto-brown dwarf

ABSTRACT The detection of Keplerian rotation is rare among Class 0 protostellar systems. We have investigated the high-density tracer HCN as a probe of the inner disc in a Class 0 proto-brown dwarf candidate. Our ALMA high angular resolution observations show the peak in the HCN (3–2) line emission arises from a compact component near the proto-brown dwarf with a small bar-like structure and a deconvolved size of ∼50 au. Radiative transfer modelling indicates that this HCN feature is tracing the innermost, dense regions in the proto-brown dwarf where a small Keplerian disc is expected to be present. The limited velocity resolution of the observations, however, makes it difficult to confirm the rotational kinematics of this feature. A brightening in the HCN emission towards the core centre suggests that HCN can survive in the gas phase in the inner, dense regions of the proto-brown dwarf. In contrast, modelling of the HCO+ (3–2) line emission indicates that it originates from the outer pseudo-disc/envelope region and is centrally depleted. HCN line emission can reveal the small-scale structures and can be an efficient observational tool to study the inner disc properties in such faint compact objects where spatially resolving the disc is nearly impossible.