As part of the ALMA Large Program “Early Planet Formation in Embedded Disks,” ^12 CO (2–1) was observed towards 19 nearby low-mass protostars. Of these objects, 15 sources are found to show molecular outflow emission. Based on their morphological and kinematical structures, the CO outflows are classified into three types: a wind-driven shell, where ambient material is swept up by a wide-angle wind from the star, a bow shock, and a slow disk wind, which is a conical or parabolic flow with onion-like velocity structure. We categorize 11 outflows as a slow disk wind, 7 as a wind-driven shell, and 1 as a bow shock. Four of these outflows were found to show signs of both slow disk wind and wind-driven shell characteristics. Five objects show misalignment between the red- and blueshifted outflows. Seven objects show significant misalignment between the outflow axis (either or both of the red- and blueshifted outflows) and the minor axis of the dust continuum emission around the protostar. For the objects showing wind-driven shell emission, we compare simple parametrized models with the observations to derive physical properties of the observed shells, such as their dynamical ages. This shows evidence of a time variability in the outflows, such as changes in their direction. In some objects, large differences are seen between the properties of the red- and blueshifted outflows, possibly indicating differences in the properties of the ambient medium with which the outflow interacts.

We present ∼8–40 μ m SOFIA-FORCAST images of seven regions of “clustered” star formation as part of the SOFIA Massive Star Formation Survey. We identify a total of 34 protostar candidates and build their spectral energy distributions (SEDs). We fit these SEDs with a grid of radiative transfer models based on the turbulent core accretion (TCA) theory to derive key protostellar properties, including initial core mass, M c , clump environment mass surface density, Σ cl , and current protostellar mass, m * . We also carry out empirical graybody (GB) estimation of Σ cl , which allows a case of restricted SED fitting within the TCA model grid. We also release version 2.0 of the open-source Python package sedcreator , which is designed to automate the aperture photometry and SED building and fitting process for sources in clustered environments, where flux contamination from close neighbors typically complicates the process. Using these updated methods, SED fitting yields values of M c ∼ 30–200 M ⊙ , Σ cl,SED ∼ 0.1–3 g cm −2 , and m * ∼ 4–50 M ⊙ . The GB fitting yields smaller values of Σ cl,GB ≲ 1 g cm −2 . From these results, we do not find evidence for a critical Σ cl needed to form massive (≳8 M ⊙ ) stars. However, we do find tentative evidence for a dearth of the most massive ( m * ≳ 30 M ⊙ ) protostars in the clustered regions, suggesting a potential impact of environment on the stellar initial mass function.

With the exponential growth of astronomical data over time, finding the needles in the haystack is becoming increasingly difficult. The next frontier for science archives is to enable searches not only on observational metadata, but also on the content of the observations themselves. As a step in this direction, we have implemented morphological image similarity search into the ALMA Science Archive (ASA). To achieve this we use self-supervised contrastive affine-transformation-independent representation learning of source morphologies with a deep neural network. For a given image on the ASA web interface, astronomers are presented with a summary view of the morphologically most similar images. Each time an astronomer selects an additional image from that view, the display is instantly updated to show the images most similar to the combination of the selected images. Each selection thus refines the similarity display according to the scientific needs of the astronomer. This is the first time image similarity search has been offered in an astronomical science archive.

We analyze the dense gas kinematics in two class 0/I protostellar cores, Per 30 and NGC 1333 IRAS 7, in the Perseus Molecular Cloud to determine whether their velocity structures are indicative of rotation. We examine the hyperfine structure of the N _2 H ^+ J = 1–0 transition by combining 3″ (900 au) Atacama Large Millimeter/submillimeter Array measurements with 9″ (2700 au) measurements from the Green Bank Telescope. We use the CASA Feather method to combine these data in order to maximize our sensitivity across spatial scales. We fit the N _2 H ^+ spectra to constrain the centroid velocity of the gas at each pixel and use these values to calculate the linear velocity gradient and specific angular momentum within apertures centered on each protostar with radii ranging from 5″ to 60″. Our results indicate that the velocity structure probed by the N _2 H ^+ emission is likely not a result of core rotation. These findings are consistent with other studies in the literature that indicate rotation is often not evident on scales ≲1000 au. We instead suggest that the velocity structure we see is a result of torques caused by irregular density distributions in these protostellar systems.

Young protostars that undergo episodic accretion can provide insight into the impact on their circumstellar environments while matter is accreted from the disk onto the protostar. IRAS 22343+7501 is a four component protostar system with one of those being a fading outbursting protostar referred to as L1251 VLA 6. Given the rarity of YSOs undergoing this type of accretion, L1251 VLA 6 can elucidate the fading phase of the post-outburst process. Here we examine structure in the disk around L1251 VLA 6 at frequencies of 10 GHz and 33 GHz with the Karl G. Jansky Very Large Array (VLA). We model the disk structure using Markov chain Monte Carlo (MCMC). This method is then combined with a parametric ray-tracing code to generate synthetic model images of an axisymmetric disk, allowing us to characterize the radial distribution of dust in the system. The results of our MCMC fit show that the most probable values for the mass and radius are consistent with values typical of Class I objects. We find that the total mass of the disk is $0.070^{+0.031}_{-0.2} \rm ~ M_{\sun}$ and investigate the conditions that could cause the accretion outburst. We conclude that the eruption is not caused by gravitational instability and consider alternative explanations and trigger mechanisms.

We present the results of the observations made within the ALMA Large Program called Early Planet Formation in Embedded disks of the Class 0 protostar GSS30 IRS3. Our observations included the 1.3 mm continuum with a resolution of 0 ''.05 (7.8 au) and several molecular species, including (CO)-C-12, (CO)-C-13, (CO)-O-18, H2CO, and c-C3H2. The dust continuum analysis unveiled a disk-shaped structure with a major axis of similar to 200 au. We observed an asymmetry in the minor axis of the continuum emission suggesting that the emission is optically thick and the disk is flared. On the other hand, we identified two prominent bumps along the major axis located at distances of 26 and 50 au from the central protostar. The origin of the bumps remains uncertain and might be an embedded substructure within the disk or the temperature distribution and not the surface density because the continuum emission is optically thick. The (CO)-C-12 emission reveals a molecular outflow consisting of three distinct components: a collimated component, an intermediate-velocity component exhibiting an hourglass shape, and a wider angle low-velocity component. We associate these components with the coexistence of a jet and a disk wind. The (CO)-O-18 emission traces both a circumstellar disk in Keplerian rotation and the infall of the rotating envelope. We measured a stellar dynamical mass of 0.35 +/- 0.09 M-circle dot.

Molecular abundances in protoplanetary disks are highly sensitive to the local physical conditions, including gas temperature, gas density, radiation field, and dust properties. Often multiple factors are intertwined, impacting the abundances of both simple and complex species. We present a new approach to understanding these chemical and physical interdependencies using machine learning. Specifically, we explore the case of CO modeled under the conditions of a generic disk and build an explanatory regression model to study the dependence of CO spatial density on the gas density, gas temperature, cosmic-ray ionization rate, X-ray ionization rate, and UV flux. Our findings indicate that combinations of parameters play a surprisingly powerful role in regulating CO abundance compared to any singular physical parameter. Moreover, in general we find the conditions in the disk are destructive toward CO. CO depletion is further enhanced in an increased cosmic-ray environment and in disks with higher initial C/O ratios. These dependencies uncovered by our new approach are consistent with previous studies, which are more modeling intensive and computationally expensive. Our work thus shows that machine learning can be a powerful tool not only for creating efficient predictive models, but also for enabling a deeper understanding of complex chemical processes.

(abridged) Protostellar outflows exhibit large variations in their structure depending on the observed gas emission. This study analyzes the atomic jet and molecular outflow in the Class I protostar, TMC1A to characterize morphology and identify previously undetected spatial features with JWST's NIRSpec IFU. In addition to identifying a large number of Fe II and H2 lines, we have detected the bipolar Fe jet by revealing, for the first time, the presence of a red-shifted atomic jet. Similarly, the red-shifted component of the H2 slower wide-angle outflow is observed. Both Fe II and H2 red-shifted emission exhibit significantly lower flux densities compared to their blue-shifted counterparts. Additionally, we report the detection of a collimated high-velocity (100 km s-1), blue-shifted H2 outflow, suggesting the presence of a molecular jet in addition to the well-known wider angle low-velocity structure. The Fe II and H2 jets show multiple intensity peaks along the jet axis, which may be associated with ongoing or recent outburst events. In addition to the variation in their intensities, the H2 wide-angle outflow exhibits a "ring"-like structure. The blue-shifted H2 outflow also shows a left-right brightness asymmetry likely due to interactions with the surrounding ambient medium and molecular outflows. Using the Fe II line ratios, the extinction along the atomic jet is estimated to be between Av = 10-30 on the blue-shifted side, with a trend of decreasing extinction with distance from the protostar. A similar Av is found for the red-shifted side, supporting the argument for an intrinsic red-blue outflow lobe asymmetry rather than environmental effects such as extinction. This intrinsic difference revealed by the unprecedented sensitivity of JWST, suggests that younger outflows already exhibit the red-blue side asymmetry more commonly observed towards jets associated with Class II disks.

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 (similar to 0.'' 08 or similar to 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 circle dot and 0.26 M circle dot 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.

Young protostellar binary systems, with expected ages less than $\sim$10$^5$ years, are little modified since birth, providing key clues to binary formation and evolution. We present a first look at the young, Class 0 binary protostellar system R CrA IRAS 32 from the Early Planet Formation in Embedded Disks (eDisk) ALMA large program, which observed the system in the 1.3 mm continuum emission, $^{12}$CO (2-1), $^{13}$CO (2-1), C$^{18}$O (2-1), SO (6$_5$-5$_4$), and nine other molecular lines that trace disk, envelope, shocks, and outflows. With a continuum resolution of $\sim$0.03$^{\prime\prime}$ ($\sim$5 au, at a distance of 150 pc), we characterize the newly discovered binary system with a separation of 207 au, their circumstellar disks, and a circumbinary disk-like structure. The circumstellar disk radii are 26.9$\pm$0.3 and 22.8$\pm$0.3 au for sources A and B, respectively, and their circumstellar disk dust masses are estimated as 22.5$\pm$1.1 and 12.4$\pm$0.6 M$_{\Earth}$. The circumstellar disks and the circumbinary structure have well aligned position angles and inclinations, indicating formation in a smooth, ordered process such as disk fragmentation. In addition, the circumstellar disks have a near/far-side asymmetry in the continuum emission suggesting that the dust has yet to settle into a thin layer near the midplane. Spectral analysis of CO isotopologues reveals outflows that originate from both of the sources and possibly from the circumbinary disk-like structure. Furthermore, we detect Keplerian rotation in the $^{13}$CO isotopologues toward both circumstellar disks and likely Keplerian rotation in the circumbinary structure; the latter suggests that it is probably a circumbinary disk.

Previous observations of B335 have presented evidence of ongoing infall in various molecular lines, e.g., HCO$^+$, HCN, CO. There have been no confirmed observations of a rotationally supported disk on scales greater than ~12~au. The presence of an outflow in B335 suggests that also a disk should be present or in formation. To constrain the earliest stages of protostellar evolution and disk formation, we aim to map the region where gas falls inwards and observationally constrain its kinematics. Furthermore, we aim to put strong limits on the size and orientation of any disk-like structure in B335. We use high angular resolution $^{13}$CO data from ALMA, and combine it with shorter-baseline archival data to produce a high-fidelity image of the infall in B335. We also revisit the imaging of high-angular resolution Band 6 continuum data to study the dust distribution in the immediate vicinity of B335. Continuum emission shows an elliptical structure (10 by 7 au) with a position angle 5 degrees east of north, consistent with the expectation for a forming disk in B335. A map of the infall velocity (as estimated from the $^{13}$CO emission), shows evidence of asymmetric infall, predominantly from the north and south. Close to the protostar, infall velocities appear to exceed free-fall velocities. 3D radiative transfer models, where the infall velocity is allowed to vary within the infall region, can explain the observed kinematics. The data suggests that a disk has started to form in B335 and that gas is falling towards that disk. However, kinematically-resolved line data towards the disk itself is needed to confirm the presence of a rotationally supported disk around this young protostar. The measured high infall velocities are not easily reconcilable with a magnetic braking scenario and suggest that there is a pressure gradient that allows the infall velocity to vary in the region.

We present Atacama Large Millimeter/submillimeter Array (ALMA) observations of the Class I source Oph IRS 63 in the context of the Early Planet Formation in Embedded Disks large program. Our ALMA observations of Oph IRS 63 show a myriad of protostellar features, such as a shell-like bipolar outflow (in 12 CO), an extended rotating envelope structure (in 13 CO), a streamer connecting the envelope to the disk (in C 18 O), and several small-scale spiral structures seen toward the edge of the dust continuum (in SO). By analyzing the velocity pattern of 13 CO and C 18 O, we measure a protostellar mass of M ⋆ = 0.5 ± 0.2 M ⊙ and confirm the presence of a disk rotating at almost Keplerian velocity that extends up to ∼260 au. These calculations also show that the gaseous disk is about four times larger than the dust disk, which could indicate dust evolution and radial drift. Furthermore, we model the C 18 O streamer and SO spiral structures as features originating from an infalling rotating structure that continuously feeds the young protostellar disk. We compute an envelope-to-disk mass infall rate of ∼10 −6 M ⊙ yr −1 and compare it to the disk-to-star mass accretion rate of ∼10 −8 M ⊙ yr −1 , from which we infer that the protostellar disk is in a mass buildup phase. At the current mass infall rate, we speculate that soon the disk will become too massive to be gravitationally stable.

We present an overview of the Large Program, “Early Planet Formation in Embedded Disks (eDisk),” conducted with the Atacama Large Millimeter/submillimeter Array (ALMA). The ubiquitous detections of substructures, particularly rings and gaps, in protoplanetary disks around T Tauri stars raise the possibility that at least some planet formation may have already started during the embedded stages of star formation. In order to address exactly how and when planet formation is initiated, the program focuses on searching for substructures in disks around 12 Class 0 and 7 Class I protostars in nearby (<200 pc) star-forming regions through 1.3 mm continuum observations at a resolution of ∼7 au (0.″04). The initial results show that the continuum emission, mostly arising from dust disks around the sample protostars, has relatively few distinctive substructures, such as rings and spirals, in marked contrast to Class II disks. The dramatic difference may suggest that substructures quickly develop in disks when the systems evolve from protostars to Class II sources, or alternatively that high optical depth of the continuum emission could obscure internal structures. Kinematic information obtained through CO isotopologue lines and other lines reveals the presence of Keplerian disks around protostars, providing us with crucial physical parameters, in particular, the dynamical mass of the central protostars. We describe the background of the eDisk program, the sample selection and their ALMA observations, and the data reduction, and we also highlight representative first-look results.

We present the first results from the Early Planet Formation in Embedded Disks Atacama Large Millimeter/submillimeter Array Large Program toward Oph IRS43, a binary system of solar mass protostars. The 1.3 mm dust continuum observations resolve a compact disk, ∼6 au radius, around the northern component and show that the disk around the southern component is even smaller, ≲3 au. CO, 13 CO, and C 18 O maps reveal a large cavity in a low-mass envelope that shows kinematic signatures of rotation and infall extending out to ∼2000 au. An expanding CO bubble centered on the extrapolated location of the source ∼130 yr ago suggests a recent outburst. Despite the small size of the disks, the overall picture is of a remarkably large and dynamically active region.

We present Atacama Large Millimeter/submillimeter Array observations of the $\sim$10 kAU environment surrounding 21 protostars in the Orion A molecular cloud tracing outflows. Our sample is composed of Class 0 to flat-spectrum protostars, spanning the full $\sim$1 Myr lifetime. We derive the angular distribution of outflow momentum and energy profiles and obtain the first two-dimensional instantaneous mass, momentum, and energy ejection rate maps using our new approach: the Pixel Flux-tracing Technique (PFT). Our results indicate that by the end of the protostellar phase, outflows will remove $\sim$2 to 4 M$_\odot$ from the surrounding $\sim$1 M$_\odot$ low-mass core. These high values indicate that outflows remove a significant amount of gas from their parent cores and continuous core accretion from larger scales is needed to replenish core material for star formation. This poses serious challenges to the concept of ``cores as well-defined mass reservoirs", and hence to the simplified core-to-star conversion prescriptions. Furthermore, we show that cavity opening angles, and momentum and energy distributions all increase with the protostar evolutionary stage. This is clear evidence that even garden-variety protostellar outflows: (a) effectively inject energy and momentum into their environments on $10$ kAU scales, and (b) significantly disrupt their natal cores, ejecting a large fraction of the mass that would have otherwise fed the nascent star. Our results support the conclusion that protostellar outflows have a direct impact on how stars get their mass, and that the natal sites of individual low-mass star formation are far more dynamic than commonly accepted theoretical paradigms.

We present observations of the Class 0 protostar IRAS 16544–1604 in CB 68 from the “Early Planet Formation in Embedded Disks (eDisk)” ALMA Large program. The ALMA observations target continuum and lines at 1.3 mm with an angular resolution of ∼5 au. The continuum image reveals a dusty protostellar disk with a radius of ∼30 au seen close to edge-on and asymmetric structures along both the major and minor axes. While the asymmetry along the minor axis can be interpreted as the effect of the dust flaring, the asymmetry along the major axis comes from a real nonaxisymmetric structure. The C 18 O image cubes clearly show the gas in the disk that follows a Keplerian rotation pattern around a ∼0.14 M ⊙ central protostar. Furthermore, there are ∼1500 au scale streamer-like features of gas connecting from northeast, north–northwest, and northwest to the disk, as well as the bending outflow as seen in the 12 CO (2–1) emission. At the apparent landing point of the NE streamer, there is SO (6 5 –5 4 ) and SiO (5–4) emission detected. The spatial and velocity structure of the NE streamer can be interpreted as a free-falling gas with a conserved specific angular momentum, and the detection of the SO and SiO emission at the tip of the streamer implies the presence of accretion shocks. Our eDisk observations have unveiled that the Class 0 protostar in CB 68 has a Keplerian-rotating disk with a flaring and nonaxisymmetric structure associated with accretion streamers and outflows.

Studying the physical and chemical conditions of young embedded disks is crucial to constrain the initial conditions for planet formation. Here we present Atacama Large Millimeter/submillimeter Array observations of dust continuum at ∼0.″06 (8 au) resolution and molecular line emission at ∼0.″17 (24 au) resolution toward the Class 0 protostar L1527 IRS from the Large Program eDisk (Early Planet Formation in Embedded Disks). The continuum emission is smooth without substructures but asymmetric along both the major and minor axes of the disk as previously observed. The detected lines of ^12 CO, ^13 CO, C ^18 O, H _2 CO, c-C _3 H _2 , SO, SiO, and DCN trace different components of the protostellar system, with a disk wind potentially visible in ^12 CO. The ^13 CO brightness temperature and the H _2 CO line ratio confirm that the disk is too warm for CO freezeout, with the snowline located at ∼350 au in the envelope. Both molecules show potential evidence of a temperature increase around the disk–envelope interface. SO seems to originate predominantly in UV-irradiated regions such as the disk surface and the outflow cavity walls rather than at the disk–envelope interface as previously suggested. Finally, the continuum asymmetry along the minor axis is consistent with the inclination derived from the large-scale (100″ or 14,000 au) outflow, but opposite to that based on the molecular jet and envelope emission, suggesting a misalignment in the system. Overall, these results highlight the importance of observing multiple molecular species in multiple transitions to characterize the physical and chemical environment of young disks.

In this Letter we report Very Long Baseline Array observations of 22 GHz water masers toward the protostar CARMA–6 , located at the center of the Serpens South young cluster. From the astrometric fits to maser spots, we derive a distance of 440.7±3.5 pc for the protostar (1% error). This represents the best direct distance determination obtained so far for an object this young and deeply embedded in this highly obscured region. Taking depth effects into account, we obtain a distance to the cluster of 440.7 ± 4.6 pc. Stars visible in the optical that have astrometric solutions in the Gaia Data Release 3 are, on the other hand, all located at the periphery of the cluster. Their mean distance of 437 −41 +51 pc is consistent within 1 σ with the value derived from maser astrometry. As the maser source is at the center of Serpens South, we finally solve the ambiguity of the distance to this region that has prevailed over the years.

Modern interferometers routinely provide radio-astronomical images down to subarcsecond resolution. However, interferometers filter out spatial scales larger than those sampled by the shortest baselines, which affects the measurement of both spatial and spectral features. Complementary single-dish data are vital for recovering the true flux distribution of spatially resolved astronomical sources with such extended emission. In this work, we provide an overview of the prominent available methods to combine single-dish and interferometric observations. We test each of these methods in the framework of the CASA data analysis software package on both synthetic continuum and observed spectral data sets. We develop a set of new assessment tools that are generally applicable to all radio-astronomical cases of data combination. Applying these new assessment diagnostics, we evaluate the methods' performance and demonstrate the significant improvement of the combined results in comparison to purely interferometric reductions. We provide combination and assessment scripts as add-on material. Our results highlight the advantage of using data combination to ensure high-quality science images of spatially resolved objects.

Protostellar disks are an ubiquitous part of the star formation process and the future sites of planet formation. As part of the Early Planet Formation in Embedded Disks large program, we present high angular resolution dust continuum (∼40 mas) and molecular line (∼150 mas) observations of the Class 0 protostar IRAS 15398–3359. The dust continuum is small, compact, and centrally peaked, while more extended dust structures are found in the outflow directions. We perform a 2D Gaussian fitting and find the deconvolved size and 2 σ radius of the dust disk to be 4.5 × 2.8 au and 3.8 au, respectively. We estimate the gas+dust disk mass assuming optically thin continuum emission to be 0.6 M J –1.8 M J , indicating a very low mass disk. The CO isotopologues trace components of the outflows and inner envelope, while SO traces a compact, rotating disk-like component. Using several rotation curve fittings on the position–velocity diagram of the SO emission, the lower limits of the protostellar mass and gas disk radius are 0.022 M ⊙ and 31.2 au, respectively, from our Modified 2 single power-law fitting. A conservative upper limit of the protostellar mass is inferred to be 0.1 M ⊙ . The protostellar mass accretion rate and the specific angular momentum at the protostellar disk edge are found to be in the range of (1.3–6.1) × 10 −6 M ⊙ yr −1 and (1.2–3.8) × 10 −4 km s −1 pc, respectively, with an age estimated between 0.4 × 10 4 yr and 7.5 × 10 4 yr. At this young age with no clear substructures in the disk, planet formation would likely not yet have started. This study highlights the importance of high-resolution observations and systematic fitting procedures when deriving dynamical properties of deeply embedded Class 0 protostars.