Search

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.

Observational constraints on dust properties in protoplanetary disks are key to better understanding disks' evolution. We continue our exploration of the protoplanetary disk around AB Aur by characterizing its dust properties. We present ALMA observations at 2.2 mm and VLA observations at 6.8 mm. Together with previous ALMA and NOEMA observations at 0.87 and 1.1 mm, these observations are used to compute global spectral index profiles as well as spectral index maps. On the interpretation side, we present the results of a simple isothermal slab model to help constrain dust properties along the ring of continuum emission. We also present results of dust radiative transfer calculations based on a disc-planet hydrodynamical simulation to explain how the azimuthal contrast ratio of the ring emission varies with millimeter wavelength. The spectral energy distribution and the radial profiles of the spectral index indicate that the radiation from the compact source towards the center is not dominated by dust thermal emission, but most likely by free-free emission originating in the radio jet: it constitutes 93% of the emission at 6.8 mm, and 37% at 0.87 mm. The protoplanetary disk has a typical spectral index of 2.3. We estimate a dust disk mass of 8$\rm \times 10^{-5}$ M$\rm _{\odot}$ which, assuming a mean gas-to-dust ratio of 40, gives a total disk mass of 3.2$\rm \times 10^{-3}$ M$\rm _{\odot}$. The azimuthal contrast ratio of the ring outside the millimeter cavity is smaller at 2.2 mm than at 1.1 mm, in agreement with previous findings. The VLA image shows several knots of $5\sigma$ emission all along the ring, which are consistent with the ring emission being nearly axisymmetric at that wavelength. The decrease in the azimuthal contrast ratio of the ring emission from 0.87 to 6.8 mm can be explained by a dust-losing decaying vortex at the outer edge of a planet gap., Comment: 13 pages, 11 figures

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 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

We present ALMA band 6/7 (1.3 mm/0.87 mm) and VLA Ka band (9 mm) observations toward NGC 2071 IR, an intermediate-mass star forming region. We characterize the continuum and associated molecular line emission towards the most luminous protostars, i.e., IRS1 and IRS3, on ~100 au (0. 2") scales. IRS1 is partly resolved in millimeter and centimeter continuum, which shows a potential disk. IRS3 has a well resolved disk appearance in millimeter continuum and is further resolved into a close binary system separated by ~40 au at 9 mm. Both sources exhibit clear velocity gradients across their disk major axes in multiple spectral lines including C18O, H2CO, SO, SO2, and complex organic molecules like CH3OH, 13CH3OH and CH3OCHO. We use an analytic method to fit the Keplerian rotation of the disks, and give constraints on physical parameters with a MCMC routine. The IRS3 binary system is estimated to have a total mass of 1.4-1.5$M_\odot$. IRS1 has a central mass of 3-5$M_\odot$ based on both kinematic modeling and its spectral energy distribution, assuming that it is dominated by a single protostar. For both IRS1 and IRS3, the inferred ejection directions from different tracers, including radio jet, water maser, molecular outflow, and H2 emission, are not always consistent, and for IRS1, these can be misaligned by ~50$^{\circ}$. IRS3 is better explained by a single precessing jet. A similar mechanism may be present in IRS1 as well but an unresolved multiple system in IRS1 is also possible., Comment: 36 pages, 21 figures, accepted by ApJ

Observed changes in protostellar brightness can be complicated to interpret. In our JCMT~Transient monitoring survey, we discovered that a young binary protostar, HOPS 373, is undergoing a modest $30\%$ brightness increase at 850 $\mu$m, caused by a factor of 1.8$-$3.3 enhancement in the accretion rate. The initial burst occurred over a few months, with a sharp rise and then shallower decay. A second rise occurred soon after the decay, and the source is still bright one year later. The mid-IR emission, the small-scale CO outflow mapped with ALMA, and the location of variable maser emission indicate that the variability is associated with the SW component. The near-infrared and NEOWISE $W1$ and $W2$ emission is located along the blueshifted CO outflow, spatially offset by $\sim3$ to $4^{\prime\prime}$ from the SW component. The $K$-band emission imaged by UKIRT shows a compact H$_2$ emission source at the edge of the outflow, with a tail tracing the outflow back to the source. The $W1$ emission, likely dominated by scattered light, brightens by 0.7 mag, consistent with expectations based on the sub-mm lightcurve. The signal of continuum variability in $K$-band and $W2$ is masked by stable H$_2$ emission, as seen in our Gemini/GNIRS spectrum, and perhaps by CO emission. These differences in emission sources complicate infrared searches for variability of the youngest protostars., Comment: 29 pages, 14 figures, accepted for publication in ApJ

At least half of a protostar's mass is accreted in the Class 0 phase, when the central protostar is deeply embedded in a dense, infalling envelope. We present the first systematic search for outbursts from Class 0 protostars in the Orion clouds. Using photometry from Spitzer/IRAC spanning 2004 to 2017, we detect three outbursts from Class 0 protostars with $\ge 2$ mag changes at 3.6 or 4.5 $\mu$m. This is comparable to the magnitude change of a known protostellar FU Ori outburst. Two are newly detected bursts from the protostars HOPS 12 and 124. The number of detections implies that Class 0 protostars burst every 438 yr, with a 95% confidence interval of 161 to 1884 yr. Combining Spitzer and WISE/NEOWISE data spanning 2004-2019, we show that the bursts persist for more than nine years with significant variability during each burst. Finally, we use $19-100$ $\mu$m photometry from SOFIA, Spitzer and Herschel to measure the amplitudes of the bursts. Based on the burst interval, a duration of 15 yr, and the range of observed amplitudes, 3-100% of the mass accretion during the Class 0 phase occurs during bursts. In total, we show that bursts from Class 0 protostars are as frequent, or even more frequent, than those from more evolved protostars. This is consistent with bursts being driven by instabilities in disks triggered by rapid mass infall. Furthermore, we find that bursts may be a significant, if not dominant, mode of mass accretion during the Class 0 phase., Comment: Accepted to ApJL

We present VLA and ALMA observations of the close (0.3" = 90 au separation) protobinary system SVS 13. We detect two small circumstellar disks (radii $\sim$12 and $\sim$9 au in dust, and $\sim$30 au in gas) with masses of $\sim$0.004-0.009 $M_{sun}$ for VLA 4A (the western component) and $\sim$0.009-0.030 $M_{sun}$ for VLA 4B (the eastern component). A circumbinary disk with prominent spiral arms extending $\sim$500 au and a mass of $\sim$0.052 $M_{sun}$ appears to be in the earliest stages of formation. The dust emission is more compact and with a very high optical depth toward VLA 4B, while toward VLA 4A the dust column density is lower, allowing the detection of stronger molecular transitions. We infer rotational temperatures of $\sim$140 K, on scales of $\sim$30 au, across the whole source, and a rich chemistry. Molecular transitions typical of hot corinos are detected toward both protostars, being stronger toward VLA 4A, with several ethylene glycol transitions detected only toward this source. There are clear velocity gradients, that we interpret in terms of infall plus rotation of the circumbinary disk, and purely rotation of the circumstellar disk of VLA 4A. We measured orbital proper motions and determined a total stellar mass of 1 $M_{sun}$. From the molecular kinematics we infer the geometry and orientation of the system, and stellar masses of $\sim$0.26 $M_{sun}$ for VLA 4A and $\sim$0.60 $M_{sun}$ for VLA 4B., Comment: 59 pages, 28 Figures, 14 Tables, accepted by The Astrophysical Journal. Includes changes after the ApJ proof corrections

We characterize protostellar multiplicity in the Orion molecular clouds using ALMA 0.87~mm and VLA 9~mm continuum surveys toward 328 protostars. These observations are sensitive to projected spatial separations as small as $\sim$20~au, and we consider source separations up to 10$^4$~au as potential companions. The overall multiplicity fraction (MF) and companion fraction (CF) for the Orion protostars are 0.30$\pm$0.03 and 0.44$\pm$0.03, respectively, considering separations from 20 to 10$^4$~au. The MFs and CFs are corrected for potential contamination by unassociated young stars using a probabilistic scheme based on the surface density of young stars around each protostar. The companion separation distribution as a whole is double peaked and inconsistent with the separation distribution of solar-type field stars, while the separation distribution of Flat Spectrum protostars is consistent solar-type field stars. The multiplicity statistics and companion separation distributions of the Perseus star-forming region are consistent with those of Orion. Based on the observed peaks in the Class 0 separations at $\sim$100~au and $\sim$10$^3$~au, we argue that multiples with separations $<$500~au are likely produced by both disk fragmentation and turbulent fragmentation with migration, and those at $\ga$10$^3$~au result primarily from turbulent fragmentation. We also find that MFs/CFs may rise from Class 0 to Flat Spectrum protostars between 100 and 10$^3$~au in regions of high YSO density. This finding may be evidence for migration of companions from $>$10$^3$~au to $<$10$^3$~au, and that some companions between 10$^3$ and 10$^4$~au must be (or become) unbound., Comment: 76 pages, 20 Figures, 10 Tables, accepted by the Astrophysical Journal

We present results from the most comprehensive radio monitoring campaign towards the closest star to our Sun, Proxima Centauri. We report 1.1 to 3.1 GHz observations with the Australian Telescope Compact Array over 18 consecutive days in April 2017. We detect radio emission from Proxima Centauri for most of the observing sessions, which spanned $\sim$1.6 orbital periods of the planet Proxima b. The radio emission is stronger at the low-frequency band, centered around 1.6 GHz, and is consistent with the expected electron-cyclotron frequency for the known star's magnetic field intensity of about 600 Gauss. The 1.6 GHz light curve shows an emission pattern that is consistent with the orbital period of the planet Proxima b around the star Proxima, with its maxima of emission happening near the quadratures. We also observed two short-duration (a few minutes) flares and a long-duration (about three days) burst whose peaks happened close to the quadratures. We find that the frequency, large degree of circular polarization, change of the sign of circular polarization, and intensity of the observed radio emission are all consistent with expectations from electron cyclotron-maser emission arising from sub-Alfv\'enic star-planet interaction. We interpret our radio observations as signatures of interaction between the planet Proxima b and its host star Proxima. We advocate for monitoring other dwarf stars with planets to eventually reveal periodic radio emission due to star-planet interaction, thus opening a new avenue for exoplanet hunting and the study of a new field of exoplanet-star plasma interaction., Comment: 14 pages, 11 figures. Accepted on December 2, 2020 in A&A (main journal), several typos corrected

The degree to which the properties of protostars are affected by environment remains an open question. To investigate this, we look at the Orion A and B molecular clouds, home to most of the protostars within 500 pc. At ~400 pc, Orion is close enough to distinguish individual protostars across a range of environments in terms of both the stellar and gas projected densities. As part of the Herschel Orion Protostar Survey (HOPS), we used the Photodetector Array Camera and Spectrometer (PACS) to map 108 partially overlapping square fields with edge lengths of 5 arcmin or 8 arcmin and measure the 70 micron and 160 micron flux densities of 338 protostars within them. In this paper we examine how these flux densities and their ratio depend on evolutionary state and environment within the Orion complex. We show that Class 0 protostars occupy a region of the 70 micron flux density versus 160 micron to 70 micron flux density ratio diagram that is distinct from their more evolved counterparts. We then present evidence that the Integral-Shaped Filament (ISF) and Orion B contain protostars with more massive envelopes than those in the more sparsely populated LDN 1641 region. This can be interpreted as evidence for increasing star formation rates in the ISF and Orion B or as a tendency for more massive envelopes to be inherited from denser birth environments. We also provide technical details about the map-making and photometric procedures used in the HOPS program., Comment: Accepted by ApJ

We present ALMA (0.87~mm and 1.3~mm) and VLA (9~mm) observations toward the candidate intermediate-mass protostar OMC2-FIR3 (HOPS-370; L$_{bol}$~314~L$_{\odot}$) at $\sim$0.1" (40~au) resolution for the continuum emission and ~0.25" (100 au) resolution of nine molecular lines. The dust continuum observed with ALMA at 0.87~mm and 1.3~mm resolve a near edge-on disk toward HOPS-370 with an apparent radius of ~100 au. The VLA observations detect both the disk in dust continuum and free-free emission extended along the jet direction. The ALMA observations of molecular lines (H$_2$CO, SO, CH$_3$OH, $^{13}$CO, C$^{18}$O, NS, and H$^{13}$CN) reveal rotation of the apparent disk surrounding HOPS-370 orthogonal to the jet/outflow direction. We fit radiative transfer models to both the dust continuum structure of the disk and molecular line kinematics of the inner envelope and disk for the H$_2$CO, CH$_3$OH, NS, and SO lines. The central protostar mass is determined to be $\sim$2.5 M_sun with a disk radius of $\sim$94~au, when fit using combinations of the H$_2$CO, CH$_3$OH, NS, and SO lines, consistent with an intermediate-mass protostar. Modeling of the dust continuum and spectral energy distribution (SED) yields a disk mass of 0.035~M$_{\odot}$ (inferred dust+gas) and a dust disk radius of 62~au, thus the dust disk may have a smaller radius than the gas disk, similar to Class II disks. In order to explain the observed luminosity with the measured protostar mass, HOPS-370 must be accreting at a rate between 1.7 and 3.2$\times$10$^{-5}$~M$_{\odot}$~yr$^{-1}$., Comment: Accepted to ApJ; 51 pages, 12 Figures, 7 Tables

AB Aur is a Herbig Ae star that hosts a prototypical transition disk. The disk shows a plethora of features connected with planet formation mechanisms. Understanding the physical and chemical characteristics of these features is crucial to advancing our knowledge of planet formation. We aim to characterize the gaseous disk around the Herbig Ae star AB Aur. A complete spectroscopic study was performed using NOEMA to determine the physical and chemical conditions. We present new observations of the continuum and 12CO, 13CO, C18O, H2CO, and SO lines. We used the integrated intensity maps and stacked spectra to derive estimates of the disk temperature. By combining our 13CO and C18O observations, we computed the gas-to-dust ratio along the disk. We also derived column density maps for the different species and used them to compute abundance maps. The results of our observations were compared with Nautilus astrochemical models. We detected continuum emission in a ring that extends from 0.6 to 2.0 arcsec, peaking at 0.97 and with a strong azimuthal asymmetry. The molecules observed show different spatial distributions, and the peaks of the distributions are not correlated with the binding energy. Using H2CO and SO lines, we derived a mean disk temperature of 39 K. We derived a gas-to-dust ratio that ranges from 10 to 40. The comparison with Nautilus models favors a disk with a low gas-to-dust ratio (40) and prominent sulfur depletion. From a very complete spectroscopic study of the prototypical disk around AB Aur, we derived, for the first time, the gas temperature and the gas-to-dust ratio along the disk, providing information that is essential to constraining hydrodynamical simulations.Moreover, we explored the gas chemistry and, in particular, the sulfur depletion. The derived sulfur depletion is dependent on the assumed C/O ratio. Our data are better explained with C/O ~ 0.7 and S/H=8e-8., Comment: 13 figures, 6 tables

We present sensitive and high angular resolution ($\sim$0.2-0.3$''$) (sub)millimeter (230 and 345 GHz) continuum and CO(2$-$1)/CO(3$-$2) line archive observations of the disk star system in UX Tauri carried out with ALMA (The Atacama Large Millimeter/Submillimeter Array). These observations reveal the gas and dusty disk surrounding the young star UX Tauri A with a large signal-to-noise ratio ($>$400 in the continuum and $>$50 in the line), and for the first time is detected the molecular gas emission associated with the disk of UX Tauri C (with a size for the disk of $<$56 au). No (sub)millimeter continuum emission is detected at 5$\sigma$-level (0.2 mJy at 0.85 mm) associated with UX Tauri C. For the component UX Tauri C, we estimate a dust disk mass of $\leq$ 0.05 M$_\oplus$. Additionally, we report a strong tidal disk interaction between both disks UX Tauri A/C, separated 360 au in projected distance. The CO line observations reveal marked spiral arms in the disk of UX Tauri A and an extended redshifted stream of gas associated with the UX Tauri C disk. No spiral arms are observed in the dust continuum emission of UX Tauri A. Assuming a Keplerian rotation we estimate the enclosed masses (disk$+$star) from their radial velocities in 1.4 $\pm$ 0.6 M$_\odot$ for UX Tauri A, and 70 $\pm$ 30 / $\sin i$ Jupiter masses for UX Tauri C (the latter coincides with the mass upper limit value for a brown dwarf). The observational evidence presented here lead us to propose that UX Tauri C is having a close approach of a possible wide, evolving and eccentric orbit around the disk of UX Tauri A causing the formation of spiral arms and the stream of molecular gas falling towards UX Tauri C., Comment: Accepted by The Astrophysical Journal

We report ALMA and VLA continuum observations that potentially identify the four youngest protostars in the Orion Molecular Clouds taken as part of the Orion VANDAM program. These are distinguished by bright, extended, irregular emission at 0.87 mm and 8 mm and are optically thick at 0.87 mm. These structures are distinct from the disk or point-like morphologies seen toward the other Orion protostars. The 0.87 mm emission implies temperatures of 41-170 K, requiring internal heating. The bright 8 mm emission implies masses of 0.5 to 1.2 M_sun assuming standard dust opacity models. One source has a Class 0 companion, while another exhibits substructure indicating a companion-candidate. Three compact outflows are detected, two of which may be driven by companions, with dynamical times of ~300 to ~400 years. The slowest outflow may be driven by a first hydrostatic core. These protostars appear to trace an early phase when the centers of collapsing fragments become optically thick to their own radiation and compression raises the gas temperature. This phase is thought to accompany the formation of hydrostatic cores. A key question is whether these structures are evolving on free fall times of ~100 years, or whether they are evolving on Kelvin-Helmholtz times of several thousand years. The number of these sources imply a lifetime of ~6000 years, in closer agreement with the Kelvin-Helmholtz time. In this case, rotational and/or magnetic support could be slowing the collapse., Comment: accepted by ApJ, 22 pages, 11 figures, 5 tables

We have conducted a survey of 328 protostars in the Orion molecular clouds with ALMA at 0.87 mm at a resolution of $\sim$0.1" (40 au), including observations with the VLA at 9 mm toward 148 protostars at a resolution of $\sim$0.08" (32 au). This is the largest multi-wavelength survey of protostars at this resolution by an order of magnitude. We use the dust continuum emission at 0.87 mm and 9 mm to measure the dust disk radii and masses toward the Class 0, Class I, and Flat Spectrum protostars, characterizing the evolution of these disk properties in the protostellar phase. The mean dust disk radii for the Class 0, Class I, and Flat Spectrum protostars are 44.9$^{+5.8}_{-3.4}$, 37.0$^{+4.9}_{-3.0}$, and 28.5$^{+3.7}_{-2.3}$ au, respectively, and the mean protostellar dust disk masses are 25.9$^{+7.7}_{-4.0}$, 14.9$^{+3.8}_{-2.2}$, 11.6$^{+3.5}_{-1.9}$ Earth masses, respectively. The decrease in dust disk masses is expected from disk evolution and accretion, but the decrease in disk radii may point to the initial conditions of star formation not leading to the systematic growth of disk radii or that radial drift is keeping the dust disk sizes small. At least 146 protostellar disks (35% out of 379 detected 0.87 mm continuum sources plus 42 non-detections) have disk radii greater than 50 au in our sample. These properties are not found to vary significantly between different regions within Orion. The protostellar dust disk mass distributions are systematically larger than that of Class II disks by a factor of $>$4, providing evidence that the cores of giant planets may need to at least begin their formation during the protostellar phase., Comment: 107 pages, 21 Figures, 11 Tables, accepted to ApJ. Version with all source figures: https://www.cv.nrao.edu/~jtobin/Orion-disks-fullfigs.pdf Reduced data available from https://dataverse.harvard.edu/dataverse/VANDAMOrion

We present ALMA (0.87~mm) and VLA (9~mm) observations toward OMC2-FIR4 and OMC2-FIR3 within the Orion integral-shaped filament that are thought to be the nearest regions of intermediate mass star formation. We characterize the continuum sources within these regions on $\sim$40~AU (0\farcs1) scales and associated molecular line emission at a factor of $\sim$30 better resolution than previous observations at similar wavelengths. We identify six compact continuum sources within OMC2-FIR4, four in OMC2-FIR3, and one additional source just outside OMC2-FIR4. This continuum emission is tracing the inner envelope and/or disk emission on less than 100~AU scales. HOPS-108 is the only protostar in OMC2-FIR4 that exhibits emission from high-excitation transitions of complex organic molecules (e.g., methanol and other lines) coincident with the continuum emission. HOPS-370 in OMC2-FIR3 with L~$\sim$~360~\lsun, also exhibits emission from high-excitation methanol and other lines. The methanol emission toward these two protostars is indicative of temperatures high enough to thermally evaporate methanol from icy dust grains; overall these protostars have characteristics similar to hot corinos. We do not identify a clear outflow from HOPS-108 in \twco, but find evidence of interaction between the outflow/jet from HOPS-370 and the OMC2-FIR4 region. The multitude of observational constraints indicate that HOPS-108 is likely a low to intermediate-mass protostar in its main mass accretion phase and it is the most luminous protostar in OMC2-FIR4. The high resolution data presented here are essential for disentangling the embedded protostars from their surrounding dusty environments and characterizing them., Comment: 32 Pages, 10 Figures, 5 Tables, accepted to ApJ

Understanding planet formation requires to discern how dust grows in protoplanetary disks. An important parameter to measure in disks is the maximum dust grain size present. This is usually estimated through measurements of the dust opacity at different millimeter wavelengths assuming optically thin emission and dust opacity dominated by absorption. However, ALMA observations have shown that these assumptions might not be correct in the case of protoplanetary disks, leading to overestimation of particle sizes and to underestimation of the disk's mass. Here, we present an analysis of high quality ALMA and VLA images of the HL Tau protoplanetary disk, covering a wide range of wavelengths, from 0.8 mm to 1 cm, and with a physical resolution of $\sim$7.35 au. We describe a procedure to analyze a set of millimeter images without any assumption about the optical depth of the emission, and including the effects of absorption and scattering in the dust opacity. This procedure allows us to obtain the dust temperature, the dust surface density and the maximum particle size at each radius. In the HL Tau disk, we found that particles have already grown up to a few millimeters in size. We detect differences in the dust properties between dark and bright rings, with dark rings containing low dust density and small dust particles. Different features in the HL Tau disk seem to have different origins. Planet-disk interactions can explain substructure at the external half of the disk, but the internal rings seem to be associated to the presence of snow lines of several molecules.

We present a detailed multi-wavelength characterization of the multi-ring disk of HD 169142. We report new ALMA observations at 3 mm and analyze them together with archival 0.89 and 1.3 mm data. Our observations resolve three out of the four rings in the disk previously seen in high-resolution ALMA data. A simple parametric model is used to estimate the radial profile of the dust optical depth, temperature, density, and particle size distribution. We find that the multiple ring features of the disk are produced by annular accumulations of large particles, probably associated with gas pressure bumps. Our model indicates that the maximum dust grain size in the rings is $\sim1$ cm, with slightly flatter power-law size distributions than the ISM-like size distribution ($p\sim3.5$) found in the gaps. In particular, the inner ring ($\sim26$ au) is associated with a strong and narrow buildup of dust particles that could harbor the necessary conditions to trigger the streaming instability. According to our analysis, the snowlines of the most important volatiles do not coincide with the observed substructures. We explore different ring formation mechanisms and find that planet-disk interactions are the most likely scenario to explain the main features of HD 169142. Overall, our multi-wavelength analysis provides some of the first unambiguous evidence of the presence of radial dust traps in the rings of HD 169142. A similar analysis in a larger sample of disks could provide key insights on the impact that disk substructures have on the dust evolution and planet formation processes., Comment: 19 pages, 8 figures, accepted for publication in ApJ

We study dust concentration in axisymmetric gas rings in protoplanetary disks. Given the gas surface density, we derived an analytical total dust surface density by taking into account the differential concentration of all the grain sizes. This model allows us to predict the local dust-to-gas mass ratio and the slope of the particle size distribution, as a function of radius. We test this analytical model comparing it with a 3D magneto-hydrodynamical simulation of dust evolution in an accretion disk. The model is also applied to the disk around HD 169142. By fitting the disk continuum observations simultaneously at $\lambda = 0.87$, 1.3, 3.0 mm, we obtain a global dust-to-gas mass ratio $\epsilon_{\rm global} = 1.05 \times 10^{-2}$ and a viscosity coefficient $\alpha = 1.35 \times 10^{-2}$. This model can be easily implemented in numerical simulations of accretion disks.

Our understanding of protoplanetary disks is rapidly departing from the classical view of a smooth, axisymmetric disk. This is in part thanks to the high angular resolution that (sub)mm observations can provide. Here we present the combined results of ALMA (0.9 mm) and VLA (7 mm) dust continuum observations toward the protoplanetary disk around the solar analogue GM Aur. Both images clearly resolve the $\sim$35 au inner cavity. The ALMA observations also reveal a fainter disk that extends up to $\sim250$ au. We model our observations using two approaches: an analytical fit to the observed deprojected visibilities, and a physical disk model that fits the SED as well as the VLA and ALMA observations. Despite not being evident in the deconvolved images, the VLA and ALMA visibilities can only be fitted with two bright rings of radii $\sim$40 and $\sim$80 au. Our physical model indicates that this morphology is the result of an accumulation or trapping of large dust grains, probably due to the presence of two pressure bumps in the disk. Even though alternative mechanisms cannot be discarded, the multiple rings suggest that forming planets may have cleared at least two gaps in the disk. Finally, our analysis suggests that the inner cavity might display different sizes at 0.9 mm and 7 mm. This discrepancy could be caused by the presence of free-free emission close to the star at 7 mm, or by a more compact accumulation of the large dust grains at the edge of the cavity., Comment: 18 pages, 8 figures, accepted for publication in ApJ

Proxima Centauri, the star closest to our Sun, is known to host at least one terrestrial planet candidate in a temperate orbit. Here we report the ALMA detection of the star at 1.3 mm wavelength and the discovery of a belt of dust orbiting around it at distances ranging between 1 and 4 au, approximately. Given the low luminosity of the Proxima Centauri star, we estimate a characteristic temperature of about 40 K for this dust, which might constitute the dust component of a small-scale analog to our solar system Kuiper belt. The estimated total mass, including dust and bodies up to 50 km in size, is of the order of 0.01 Earth masses, which is similar to that of the solar Kuiper belt. Our data also show a hint of warmer dust closer to the star. We also find signs of two additional features that might be associated with the Proxima Centauri system, which, however, still require further observations to be confirmed: an outer extremely cold (about 10 K) belt around the star at about 30 au, whose orbital plane is tilted about 45 degrees with respect to the plane of the sky; and additionally, we marginally detect a compact 1.3 mm emission source at a projected distance of about 1.2 arcsec from the star, whose nature is still unknown., Comment: 9 pages, 4 figures, accepted for publication in The Astrophysical Journal Letters

The Herschel Orion Protostar Survey obtained well-sampled 1.2 - 870 micron spectral energy distributions (SEDs) of over 300 protostars in the Orion molecular clouds, home to most of the young stellar objects (YSOs) in the nearest 500 pc. We plot the bolometric luminosities and temperatures for 330 Orion YSOs, 315 of which have bolometric temperatures characteristic of protostars. The histogram of bolometric temperature is roughly flat; 29% of the protostars are in Class 0. The median luminosity decreases by a factor of four with increasing bolometric temperature; consequently, the Class 0 protostars are systematically brighter than the Class I protostars, with a median luminosity of 2.3 L_sun as opposed to 0.87 L_sun. At a given bolometric temperature, the scatter in luminosities is three orders of magnitude. Using fits to the SEDs, we analyze how the luminosities corrected for inclination and foreground reddening relate to the mass in the inner 2500 AU of the best-fit model envelopes. The histogram of envelope mass is roughly flat, while the median corrected luminosity peaks at 15 L_sun for young envelopes and falls to 1.7 L_sun for late-stage protostars with remnant envelopes. The spread in luminosity at each envelope mass is three orders of magnitude. Envelope masses that decline exponentially with time explain the flat mass histogram and the decrease in luminosity, while the formation of a range of stellar masses explains the dispersion in luminosity., Comment: Accepted by ApJ, 16 pages, 14 figures

We carried out multiwavelength (0.7-5 cm), multiepoch (1994-2015) Very Large Array (VLA) observations toward the region enclosing the bright far-IR sources FIR 3 (HOPS 370) and FIR 4 (HOPS 108) in OMC-2. We report the detection of 10 radio sources, seven of them identified as young stellar objects. We image a well-collimated radio jet with a thermal free-free core (VLA 11) associated with the Class I intermediate-mass protostar HOPS 370. The jet presents several knots (VLA 12N, 12C, 12S) of non-thermal radio emission (likely synchrotron from shock-accelerated relativistic electrons) at distances of ~7,500-12,500 au from the protostar, in a region where other shock tracers have been previously identified. These knots are moving away from the HOPS 370 protostar at ~ 100 km/s. The Class 0 protostar HOPS 108, which itself is detected as an independent, kinematically decoupled radio source, falls in the path of these non-thermal radio knots. These results favor the previously proposed scenario where the formation of HOPS 108 has been triggered by the impact of the HOPS 370 outflow with a dense clump. However, HOPS 108 presents a large proper motion velocity of ~ 30 km/s, similar to that of other runaway stars in Orion, whose origin would be puzzling within this scenario. Alternatively, an apparent proper motion could result because of changes in the position of the centroid of the source due to blending with nearby extended emission, variations in the source shape, and /or opacity effects., Comment: 16 pages, 4 figures, accepted for publication in The Astrophysical Journal

We report Very Large Array observations at 7 mm, 9 mm, and 3 cm toward the pre-transitional disk of the Herbig Ae star HD 169142. These observations have allowed us to study the mm emission of this disk with the highest angular resolution so far ($0\rlap."12\times0\rlap."09$, or 14 au$\times$11 au, at 7 mm). Our 7 and 9 mm images show a narrow ring of emission at a radius of $\sim25$ au tracing the outer edge of the inner gap. This ring presents an asymmetric morphology that could be produced by dynamical interactions between the disk and forming planets. Additionally, the azimuthally averaged radial intensity profiles of the 7 and 9 mm images confirm the presence of the previously reported gap at $\sim45$ au, and reveal a new gap at $\sim85$ au. We analyzed archival DCO$^+$(3-2) and C$^{18}$O(2-1) ALMA observations, showing that the CO snowline is located very close to this third outer gap. This suggests that growth and accumulation of large dust grains close to the CO snowline could be the mechanism responsible for this proposed outer gap. Finally, a compact source of emission is detected at 7 mm, 9 mm, and 3 cm toward the center of the disk. Its flux density and spectral index indicate that it is dominated by free-free emission from ionized gas, which could be associated with either the photoionization of the inner disk, an independent object, or an ionized jet., Comment: 12 pages, 6 figures, accepted for publication in ApJ

Photoevaporation is probably the main agent for gas dispersal during the last stages of protoplanetary disk evolution. However, the overall mass loss rate in the photoevaporative wind and its driving mechanism are still not well understood. Here we report multi-configuration Very Large Array observations at 0.7, 3, and 5 cm towards the transitional disk of GM Aur. Our radio continuum observations allow us to image and spatially resolve, for the first time, the three main components at work in this stage of the disk evolution: the disk of dust, the ionized radio jet perpendicular to it, and the photoevaporative wind arising from the disk. The mass loss rate inferred from the flux density of the radio jet is consistent with the ratio between ejection and accretion rates found in younger objects, suggesting that transitional disks can power collimated ejections of material apparently following the same physical mechanisms as much younger protostars. Our results indicate that extreme-UV (EUV) radiation is the main ionizing mechanism of the photoevaporative wind traced by the free-free emission. The required low EUV photon luminosity of $\sim6\times10^{40}$ s$^{-1}$ would produce a photoevaporation rate of only $\dot{M}_{\rm w,EUV}\simeq1.3\times10^{-10}~M_{\odot}$ yr$^{-1}$. Therefore, other mechanisms are required to disperse the disk in the timescale imposed by observations., Comment: 20 pages, 3 figures, accepted for publication in ApJ

We report the discovery of a dwarf protoplanetary disk around the star XZ Tau B that shows all the features of a classical transitional disk but on a much smaller scale. The disk has been imaged with the Atacama Large Millimeter/Submillimeter Array (ALMA), revealing that its dust emission has a quite small radius of ~ 3.4 au and presents a central cavity of ~ 1.3 au in radius that we attribute to clearing by a compact system of orbiting (proto)planets. Given the very small radii involved, evolution is expected to be much faster in this disk (observable changes in a few months) than in classical disks (observable changes requiring decades) and easy to monitor with observations in the near future. From our modeling we estimate that the mass of the disk is large enough to form a compact planetary system., Comment: To appear in ApJL, 15 pages

The first long-baseline ALMA campaign resolved the disk around the young star HL Tau into a number of axisymmetric bright and dark rings. Despite the very young age of HL Tau these structures have been interpreted as signatures for the presence of (proto)planets. The ALMA images triggered numerous theoretical studies based on disk-planet interactions, magnetically driven disk structures, and grain evolution. Of special interest are the inner parts of disks, where terrestrial planets are expected to form. However, the emission from these regions in HL Tau turned out to be optically thick at all ALMA wavelengths, preventing the derivation of surface density profiles and grain size distributions. Here, we present the most sensitive images of HL Tau obtained to date with the Karl G. Jansky Very Large Array at 7.0 mm wavelength with a spatial resolution comparable to the ALMA images. At this long wavelength the dust emission from HL Tau is optically thin, allowing a comprehensive study of the inner disk. We obtain a total disk dust mass of 0.001 - 0.003 Msun, depending on the assumed opacity and disk temperature. Our optically thin data also indicate fast grain growth, fragmentation, and formation of dense clumps in the inner densest parts of the disk. Our results suggest that the HL Tau disk may be actually in a very early stage of planetary formation, with planets not already formed in the gaps but in the process of future formation in the bright rings., Comment: Accepted by The Astrophysical Journal Letters

(abridged) The HH 80/81/80N jet extends from the HH 80 object to the recently discovered Source 34 and has a total projected jet size of 10.3 pc, constituting the largest collimated radio-jet system known so far. It is powered by IRAS 18162-2048 associated with a massive young stellar object. We report 6 cm JVLA observations that, compared with previous 6 cm VLA observations carried out in 1989, allow us to derive proper motions of the HH 80, HH 81 and HH 80N radio knots located about 2.5 pc away in projection from the powering source. For the first time, we measure proper motions of the optically obscured HH 80N object providing evidence that HH 81, 80 and 80N are associated with the same radio-jet. We derived tangential velocities of these HH objects between 260 and 350 km/s, significantly lower than those for the radio knots of the jet close to the powering source (600-1400 km/s) derived in a previous work, suggesting that the jet material is slowing down due to a strong interaction with the ambient medium. The HH 80 and HH 80N emission at 6 cm is, at least in part, probably synchrotron radiation produced by relativistic electrons in a magnetic field of 1 mG. If these electrons are accelerated in a reverse adiabatic shock, we estimate a jet total density of $\lesssim1000$ cm$^{-3}$. All these features are consistent with a jet emanating from a high mass protostar and make evident its capability of accelerating particles up to relativistic velocities., Comment: 15 pages, 3 figures, 1 table, accepted for publication in ApJ

We present L' and J-band high-contrast observations of HD169142, obtained with the VLT/NACO AGPM vector vortex coronagraph and the Gemini Planet Imager, respectively. A source located at 0".156+/-0".032 north of the host star (PA=7.4+/-11.3 degrees) appears in the final reduced L' image. At the distance of the star (~145 pc), this angular separation corresponds to a physical separation of 22.7+/-4.7 AU, locating the source within the recently resolved inner cavity of the transition disk. The source has a brightness of L'=12.2+/-0.5 mag, whereas it is not detected in the J band (J>13.8 mag). If its L' brightness arose solely from the photosphere of a companion and given the J-L' color constraints, it would correspond to a 28-32 MJupiter object at the age of the star, according to the COND models. Ongoing accretion activity of the star suggests, however, that gas is left in the inner disk cavity from which the companion could also be accreting. In this case the object could be lower in mass and its luminosity enhanced by the accretion process and by a circumplanetary disk. A lower mass object is more consistent with the observed cavity width. Finally, the observations enable us to place an upper limit on the L'-band flux of a second companion candidate orbiting in the disk annular gap at ~50 AU, as suggested by millimeter observations. If the second companion is also confirmed, HD169142 might be forming a planetary system, with at least two companions opening gaps and possibly interacting with each other., Comment: Accepted to ApJL, see also Biller et al. 2014

We present Very Large Array observations at 7 mm that trace the thermal emission of large dust grains in the HD 169142 protoplanetary disk. Our images show a ring of enhanced emission of radius ~25-30 AU, whose inner region is devoid of detectable 7 mm emission. We interpret this ring as tracing the rim of an inner cavity or gap, possibly created by a planet or a substellar companion. The ring appears asymmetric, with the western part significantly brighter than the eastern one. This azimuthal asymmetry is reminiscent of the lopsided structures that are expected to be produced as a consequence of trapping of large dust grains. Our observations also reveal an outer annular gap at radii from ~40 to ~70 AU. Unlike other sources, the radii of the inner cavity, the ring, and the outer gap observed in the 7 mm images, which trace preferentially the distribution of large (mm/cm sized) dust grains, coincide with those obtained from a previous near-infrared polarimetric image, which traces scattered light from small (micron- sized) dust grains. We model the broad-band spectral energy distribution and the 7 mm images to constrain the disk physical structure. From this modeling we infer the presence of a small (radius ~0.6 AU) residual disk inside the central cavity, indicating that the HD 169142 disk is a pre-transitional disk. The distribution of dust in three annuli with gaps in between them suggests that the disk in HD 169142 is being disrupted by at least two planets or substellar objects., Comment: Accepted by ApJ Letters, 16 pages, 3 figures, ApJ Letters 2014

Edge-on protostars are valuable for understanding the disk and envelope properties of embedded young stellar objects, since the disk, envelope, and envelope cavities are all distinctly visible in resolved images and well constrained in modeling. Comparing 2MASS, WISE, Spitzer, Herschel, APEX, and IRAM photometry and limits from 1.25 to 1200 microns, Spitzer spectroscopy from 5 to 40 microns, and high-resolution Hubble imaging at 1.60 and 2.05 microns to radiative transfer modeling, we determine envelope and disk properties for the Class I protostar HOPS 136, an edge-on source in Orion's Lynds 1641 region. The source has a bolometric luminosity of 0.8 L_sun, a bolometric temperature of 170 K, and a ratio of submillimeter to bolometric luminosity of 0.8%. Via modeling, we find a total luminosity of 4.7 L_sun (larger than the observed luminosity due to extinction by the disk), an envelope mass of 0.06 M_sun, and a disk radius and mass of 450 AU and 0.002 M_sun. The stellar mass is highly uncertain but is estimated to fall between 0.4 and 0.5 M_sun. To reproduce the flux and wavelength of the near-infrared scattered-light peak in the spectral energy distribution, we require 5.4 x 10^-5 M_sun of gas and dust in each cavity. The disk has a large radius and a mass typical of more evolved T Tauri stars in spite of the significant remaining envelope. HOPS 136 appears to be a key link between the protostellar and optically revealed stages of star formation., Comment: Accepted by ApJ, 12 pages, 6 figures

We present VLA NH3 and PdBI NH2D and H13NC observations of the star forming core ahead of HH 80N, the optically obscured northern counterpart of the Herbig-Haro objects HH 80/81. The main goal is to determine the kinematical information of the high density regions of the core ($n\lesssim 10^5$ cm$^{-3}$), missed in previous works due to the depletion of the species observed (e.g. CS). The obtained maps show different kinematical signatures between the eastern and western parts of the core, suggesting a possible dynamical interaction of the core with the HH 80/81/80N outflow. The analysis of the Position-Velocity (PV) plots of these species rules out a previous interpretation of having a molecular ring-like structure of $6\times 10^4$ AU of radius traced by CS infalling onto a central protostar found in the core (IRS1). High degree of deuteration, with respect to the central part of the core harboring IRS1, is derived in the eastern part, where a dust condensation (SE) is located. This deuteration trend of NH3 suggests that SE is in a prestellar evolutionary stage, earlier than that of the IRS1. Since SE is the closest condensation to the HH 80N/81/80N outflow, in case of having outflow-core dynamical interaction, it should be perturbed first and be the most evolved condensation in the core. Therefore, the derived evolutionary sequence for SE and IRS1 makes the outflow triggered star formation on IRS1 unlikely., Comment: 19 pages, 5 figures, 2 tables, accepted for publication in ApJ

We perform a census of the reddest, and potentially youngest, protostars in the Orion molecular clouds using data obtained with the PACS instrument onboard the Herschel Space Observatory and the LABOCA and SABOCA instruments on APEX as part of the Herschel Orion Protostar Survey (HOPS). A total of 55 new protostar candidates are detected at 70 um and 160 um that are either too faint (m24 > 7 mag) to be reliably classified as protostars or undetected in the Spitzer/MIPS 24 um band. We find that the 11 reddest protostar candidates with log (lambda F_lambda 70) / (lambda F_lambda 24) > 1.65 are free of contamination and can thus be reliably explained as protostars. The remaining 44 sources have less extreme 70/24 colors, fainter 70 um fluxes, and higher levels of contamination. Taking the previously known sample of Spitzer protostars and the new sample together, we find 18 sources that have log (lambda F_lambda 70) / (lambda F_lambda 24) > 1.65; we name these sources "PACS Bright Red sources", or PBRs. Our analysis reveals that the PBRs sample is composed of Class 0 like sources characterized by very red SEDs (T_bol < 45 K) and large values of sub-millimeter fluxes (L_smm/L_bol > 0.6%). Modified black-body fits to the SEDs provide lower limits to the envelope masses of 0.2 M_sun to 2 M_sun and luminosities of 0.7 L_sun to 10 L_sun. Based on these properties, and a comparison of the SEDs with radiative transfer models of protostars, we conclude that the PBRs are most likely extreme Class 0 objects distinguished by higher than typical envelope densities and hence, high mass infall rates., Comment: Accepted to Apj; 36 pages, 25 figures, 8 tables

Surveys with the Spitzer and Herschel space observatories are now enabling the discovery and characterization of large samples of protostars in nearby molecular clouds, providing the observational basis for a detailed understanding of star formation in diverse environments. We are pursuing this goal with the Herschel Orion Protostar Survey (HOPS), which targets 328 Spitzer-identified protostars in the Orion molecular clouds, the largest star-forming region in the nearest 500 pc. The sample encompasses all phases of protostellar evolution and a wide range of formation environments, from dense clusters to relative isolation. With a grid of radiative transfer models, we fit the 1-870 micron spectral energy distributions (SEDs) of the protostars to estimate their envelope densities, cavity opening angles, inclinations, and total luminosities. After correcting the bolometric luminosities and temperatures of the sources for foreground extinction and inclination, we find a spread of several orders of magnitude in luminosity at all evolutionary states, a constant median luminosity over the more evolved stages, and a possible deficit of high-inclination, rapidly infalling envelopes among the Spitzer-identified sample. We have detected over 100 new sources in the Herschel images; some of them may fill this deficit. We also report results from modeling the pre- and post-outburst 1-870 micron SEDs of V2775 Ori (HOPS 223), a known FU Orionis outburster in the sample. It is the least luminous FU Ori star with a protostellar envelope., Comment: Accepted to Astronomische Nachrichten (Cool Stars 17 proceedings)

Individual outbursting young stars are important laboratories for studying the physics of episodic accretion and the extent to which this phenomenon can explain the luminosity distribution of protostars. We present new and archival data for V2775 Ori (HOPS 223), a protostar in the L 1641 region of the Orion molecular clouds that was discovered by Caratti o Garatti et al. (2011) to have recently undergone an order-of-magnitude increase in luminosity. Our near-infrared spectra of the source have strong blueshifted He I 10830 absorption, strong H2O and CO absorption, and no H I emission, all typical of FU Orionis sources. With data from IRTF, 2MASS, HST, Spitzer, WISE, Herschel, and APEX that span from 1 to 70 microns pre-outburst and from 1 to 870 microns post-outburst, we estimate that the outburst began between 2005 April and 2007 March. We also model the pre- and post-outburst spectral energy distributions of the source, finding it to be in the late stages of accreting its envelope with a disk-to-star accretion rate that increased from about 2x10^-6 M_sun/yr to about 10^-5 M_sun/yr during the outburst. The post-outburst luminosity at the epoch of the FU Orionis-like near-IR spectra is 28 L_sun, making V2775 Ori the least luminous documented FU Orionis outburster with a protostellar envelope. The existence of low-luminosity outbursts supports the notion that a range of episiodic accretion phenomena can partially explain the observed spread in protostellar luminosities., Comment: Accepted to ApJ, 13 pages, 6 figures

We present new high-angular resolution observations toward the driving source of the HH 80-81 jet (IRAS 18162-2048). Continuum emission was observed with the Very Large Array at 7 mm and 1.3 cm, and with the Submillimeter Array at 860 microns, with angular resolutions of ~0"1 and ~0"8 respectively. Submillimeter observations of the sulfur oxide (SO) molecule are reported as well. At 1.3 cm the emission traces the well-known radio jet, while at 7 mm the continuum morphology is quadrupolar and seems to be produced by a combination of free-free and dust emission. An elongated structure perpendicular to the jet remains in the 7 mm image after subtraction of the free-free contribution. This structure is interpreted as a compact accretion disk of ~200 AU radius. Our interpretation is favored by the presence of rotation in our SO observations observed at larger scales. The observations presented here add to the small list of cases where the hundred-AU scale emission from a circumstellar disk around a massive protostar has been resolved., Comment: Accepted for publication in ApJ Letters

We examine eight young stellar objects in the OMC-2 star forming region based on observations from the SOFIA/FORCAST early science phase, the Spitzer Space Telescope, the Herschel Space Observatory, 2MASS, APEX, and other results in the literature. We show the spectral energy distributions of these objects from near-infrared to millimeter wavelengths, and compare the SEDs with those of sheet collapse models of protostars and circumstellar disks. Four of the objects can be modelled as protostars with infalling envelopes, two as young stars surrounded by disks, and the remaining two objects have double-peaked SEDs. We model the double-peaked sources as binaries containing a young star with a disk and a protostar. The six most luminous sources are found in a dense group within a 0.15 x 0.25 pc region; these sources have luminosities ranging from 300 L_sun to 20 L_sun. The most embedded source (OMC-2 FIR 4) can be fit by a class 0 protostar model having a luminosity of ~50 L_sun and mass infall rate of ~10^-4 solar masses per year., Comment: Accepted by ApJ Letters

We present observations of continuum emission in the MIR to mm wavelength range, complemented with ammonia observations, of the dense core ahead of the radio Herbig Haro object HH 80N, found in the GGD 27 region. The continuum emission in all the observed bands peaks at the same position, consistent with the presence of an embedded object, HH 80N-IRS1, within the core. The distribution of the VLA ammonia emission is well correlated with that of the dust, suggesting that photochemical effects caused by the nearby Herbig Haro object do not play an important role in shaping this particular molecular emission. In order to unveil the nature of HH 80N-IRS1 we analyzed the continuum data of this source, using self-consistent models of protostellar collapse. We find that a young protostar surrounded by a slowly rotating collapsing envelope of radius $\sim$0.08 pc and 20 \msun\ plus a circumstellar disk of radius $\sim$300 AU and 0.6 \msun\ provide a good fit to the observed spectral energy distribution and to the maps at 350 \micron, 1.2 mm and 3.5 mm of HH 80N-IRS1. Besides, the APEX and PdBI continuum maps at 350 \micron\ and 3.5 mm, respectively, reveal additional clumps in the continuum emission. Given the modeling results and the observed morphology of the emission, we propose a scenario consisting of a central embedded Class 0 object, HH 80N-IRS1, with the rest of material of the HH 80N core possibly undergoing fragmentation that may lead to the formation of several protostars, Comment: 33 pages, 8 figures, 5 tables, accepted for publication in ApJ

Synchrotron emission is commonly found in relativistic jets from active galactic nuclei (AGNs) and microquasars, but so far its presence in jets from young stellar objects (YSOs) has not been proved. Here, we present evidence of polarized synchrotron emission arising from the jet of a YSO. The apparent magnetic field, with strength of ~0.2 milligauss, is parallel to the jet axis, and the polarization degree increases towards the jet edges, as expected for a confining helical magnetic field configuration. These characteristics are similar to those found in AGN jets, hinting at a common origin of all astrophysical jets., Comment: First published in Science on Nov 26, 2010 (vol. 330, pp.1209-1212). Includes Supporting Online Material. For original publication, see http://www.sciencemag.org/content/330/6008/1209.abstract?sid=4c98088f-eee1-4754-bfff-b5b468351e4b

We present a model aimed to reproduce the observed spectral energy distribution (SED) as well as the ammonia line emission of the G31.41+0.31 hot core. The core is modeled as an infalling envelope onto a massive star that is undergoing an intense accretion phase. We assume an envelope with a density and velocity structure resulting from the dynamical collapse of a singular logatropic sphere. The stellar and envelope physical properties are determined by fitting the SED. From these physical conditions, the ammonia line emission is calculated and compared with subarcsecond resolution VLA data of the (4,4) transition. The only free parameter in this line fitting is the ammonia abundance. The observed properties of the NH3(4,4) lines and their spatial distribution can be well reproduced provided it is taken into account the steep increase of the abundance in the hotter (> 100 K), inner regions of the core produced by the sublimation of icy mantles where ammonia molecules are trapped. The model predictions for the (2,2), (4,4), and (5,5) transitions are also in reasonably agreement with the single-dish spectra available in the literature. The best fit is obtained for a model with a star of 25 Msun, a mass accretion rate of 0.003 Msun/yr, and a total luminosity of 200,000 Lsun. The gas-phase ammonia abundance ranges from 2 \times 10^{-8} in the outer region to 3 \times 10^{-6} in the inner region. To our knowledge, this is the first time that the dust and molecular line data of a hot molecular core, including subarcsecond resolution data that spatially resolve the structure of the core, have been simultaneously explained by a physically self-consistent model. This modeling shows that massive protostars are able to excite high excitation ammonia transitions up to the outer edge (30,000 AU) of the large scale envelope., Comment: 44 pages, 8 figures, To appear in The ApJ

We present the results of a detailed modeling aimed to reproduce the spectral energy distribution (SED) of dust and molecular line emission of massive protostars under the hypothesis that they form via an accretion process. We model the emission originated in the infalling envelopes at scales smaller than 0.1 pc from the central protostar. To do that, we compared our model results with observational data covering a wide range of wavelengths, paying special attention to the high angular resolution mid-infrared data obtained with the Gemini Observatory and the ammonia line emission observed with the VLA at centimeter wavelengths. We have explored two kind of model envelopes. In the first kind of models, spherical symmetry is assumed and the SED as well as the ammonia emission of the infalling envelope are calculated. In this way, the temperature, density, velocity, velocity dispersion, and ammonia abundance variations along the core can be obtained. The second approach takes into account deviations from the spherical symmetry, and parameters such as the rotation, degree of elongation of the core, or inclination of the system can be constrained through the SED fitting. Using these two approaches we have been able to model the formation of massive stars with a degree of detail similar to that reached for the low mass stars., Comment: 8 pages, 5 figures. Conference: "Massive Star Formation: Observations confront Theory" 2007, Heidelberg Germany. To appear in Astronomical Society of the Pacific

We present high angular resolution Very Large Array multi-epoch continuum observations at 3.6 cm and 7 mm towards the core of the L723 multipolar outflow revealing a multiple system of four radio sources suspected to be YSOs in a region of only ~4 arcsecs (1200 AU) in extent. The 3.6 cm observations show that the previously detected source VLA 2 contains a close (separation ~0.29 arcsecs or ~90 AU) radio binary, with components (A and B) along a position angle of ~150 degrees. The northern component (VLA 2A) of this binary system is also detected in the 7 mm observations, with a positive spectral index between 3.6 cm and 7 mm. In addition, the source VLA 2A is associated with extended emission along a position angle of ~115 degrees, that we interpret as outflowing shock-ionized gas that is exciting a system of HH objects with the same position angle. A third, weak 3.6 cm source, VLA 2C, that is detected also at 7 mm, is located ~0.7 arcsecs northeast of VLA 2A, and is possibly associated with the water maser emission in the region. The 7 mm observations reveal the presence of an additional source, VLA 2D, located ~3.5 arcsecs southeast of VLA 2A, and with a 1.35 mm counterpart. All these radio continuum sources have a positive spectral index, compatible with them being YSOs. We also propose that the high velocity CO emission observed in the region could be the superposition of multiple outflows (at least three independent bipolar outflows) excited by the YSOs located at the core, instead of the previous interpretations in terms of only one or two outflows., Comment: Accepted for publication in The Astrophysical Journal (2007 December 6)

This is a report on detailed modeling of young high-mass protostellar candidates during their most embedded and obscured phases. We performed narrowband mid-infrared imaging of three candidate high-mass protostellar objects in G11.94-0.62, G29.96-0.02, and G45.07+0.13 at Gemini Observatory using the Thermal-Region Camera and Spectrograph (T-ReCS). The sources were imaged through up to 11 narrowband filters, sampling their SEDs over the entire 2-25um infrared range. For the first time, we have fit the observed SEDs of massive protostars with models that take into account departures from spherical symmetry in the infalling envelopes. In this way, we have been able to back out of the models detailed physical parameters for these earliest stages of massive stellar life. Our detailed modeling suggests that massive star formation can proceed in a way very similar to the formation of low-mass stars., Comment: 23 pages, 7 figures, Accepted for publication by ApJ. For a version with high-resolution images go here http://www.ctio.noao.edu/~debuizer/

We present Very Large Array observations at 7 mm of the sources IRAS 2A, IRAS 2B, MMS2, MMS3 and SVS 13, in the NGC1333 region. SVS 13 is a young close binary system whose components are separated by 65 AU in projection. Our high angular resolution observations reveal that only one of the components of the SVS 13 system (VLA 4B) is associated with detectable circumstellar dust emission. This result is in contrast with the well known case of L1551 IRS5, a binary system of two protostars separated by 45 AU, where each component is associated with a disk of dust. Both in SVS 13 and in L1551 IRS5 the emission apparently arises from compact accretion disks, smaller than those observed around single stars, but still massive enough to form planetary systems like the solar one. These observational results confirm that the formation of planets can occur in close binary systems, either in one or in both components of the system, depending on the specific angular momentum of the infalling material., Comment: 12 pages, 3 figures, to be published in ApJL