Search

Context. Historically, due to instrumental limitations and a lack of disk detections, the structure of the transition from the envelope to the rotationally supported disk has been poorly studied. This is now possible with ALMA through observations of CO isotopologues and tracers of freezeout. Class 0 sources are ideal for such studies given their almost intact envelope and young disk. Aims. The structure of the disk-envelope interface of the prototypical Class 0 source, VLA1623A, which has a confirmed Keplerian disk, is constrained through modeling and analysis of ALMA observations of DCO+ (3-2) and C18O (2-1) rotational lines. Methods. The physical structure of VLA1623 is obtained from the large-scale spectral energy distribution (SED) and continuum radiative transfer. An analytic model using a simple network coupled with radial density and temperature profiles is used as input for a 2D line radiative transfer calculation for comparison with the ALMA Cycle 0 12-m array and Cycle 2 ACA observations of VLA1623. Results. The DCO+ emission shows a clumpy structure bordering VLA1623A's Keplerian disk. This suggests a cold ring-like structure at the disk-envelope interface. The radial position of the observed DCO+ peak is reproduced in our model only if the region's temperature is between 11 K and 16 K, lower than expected from models constrained by continuum data and source SED. Altering the density profile has little e ect on the DCO+ peak position, but increased density is needed to reproduce the observed C18O tracing the disk. Conclusions. The observed DCO+ (3-2) emission around VLA1623A is the product of shadowing of the envelope by the disk observed in C18O. Disk-shadowing causes a drop in the gas temperature outside of the disk on >200 AU scales, encouraging the production of deuterated molecules. This indicates that the physical structure of the disk-envelope interface differs from the rest of the envelope, highlighting the drastic impact that the disk has on the envelope and temperature structure. The results presented here show that DCO+ is an excellent cold temperature tracer. [ABSTRACT FROM AUTHOR]

Context. Transitional disks with large dust cavities are important laboratories in which to study planet formation and disk evolution. Cold gas may still be present inside these cavities, but quantying this gas is challenging. The gas content is important for constraining the origin of the dust cavity. Aims. We use Atacama Large Millimeter/submillimeter Array (ALMA) observations of 12CO 6-5 and 690 GHz (Band 9) continuum of five well-studied transitional disks. In addition, we analyze previously published Band 7 observations of a disk in the 12CO 3-2 line and 345 GHz continuum. The observations are used to set constraints on the gas and dust surface density profiles, in particular, the drop δgas of the gas density inside the dust cavity. Methods. The physical-chemical modeling code DALI was used to simultaneously analyze the gas and dust images. We modeled SR21, HD135344B, LkCa15, SR24S, and RXJ1615-3255 (Band 9) and J1604-2130 (Band 7). The spectral energy distribution and continuum visibility curve constrain the dust surface density. Then we used the same model to calculate the 12CO emission, which we compared with the observations through spectra and intensity cuts. The amount of gas inside the cavity was quantified by varying the δgas parameter. Results. Model fits to the dust and gas indicate that gas is still present inside the dust cavity for all disks, but at a reduced level. The gas surface density drops inside the cavity by at least a factor 10, while the dust density drops by at least a factor 1000. Disk masses are comparable with previous estimates from the literature, cavity radii are found to be smaller than in the data obtained with the 345 GHz SubMillimeter Array. Conclusions. The derived gas surface density profiles suggest that the cavity was cleared by one or more companions in all cases, which trapped the millimeter-sized dust at the edge of the cavity. [ABSTRACT FROM AUTHOR]

Context. Recent simulations have explored different ways to form accretion disks around low-mass stars. However, it has been difficult to differentiate between the proposed mechanisms because of a lack of observable predictions from these numerical studies. Aims. We aim to present observables that can differentiate a rotationally supported disk from an infalling rotating envelope toward deeply embedded young stellar objects (Menv > Mdisk) and infer their masses and sizes. Methods. Two 3D magnetohydrodynamics (MHD) formation simulations are studied with a rotationally supported disk (RSD) forming in one but not the other (where a pseudo-disk is formed instead), together with the 2D semi-analytical model. We determine the dust temperature structure through continuum radiative transfer RADMC3D modeling. A simple temperature-dependent CO abundance structure is adopted and synthetic spectrally resolved submm rotational molecular lines up to Ju = 10 are compared with existing data to provide predictions for future ALMA observations. Results. The 3D MHD simulations and 2D semi-analytical model predict similar compact components in continuum if observed at the spatial resolutions of 0.5-1'' (70-140 AU) typical of the observations to date. A spatial resolution of ~14 AU and high dynamic range (>1000) are required in order to differentiate between RSD and pseudo-disk formation scenarios in the continuum. The first moment maps of the molecular lines show a blue- to red-shifted velocity gradient along the major axis of the flattened structure in the case of RSD formation, as expected, whereas it is along the minor axis in the case of a pseudo-disk. The peak position-velocity diagrams indicate that the pseudo-disk shows a flatter velocity profile with radius than does an RSD. On larger scales, the CO isotopolog line profiles within large (>9'') beams are similar and are narrower than the observed line widths of low-J (2-1 and 3-2) lines, indicating significant turbulence in the large-scale envelopes. However a forming RSD can provide the observed line widths of high-J (6-5, 9-8, and 10-9) lines. Thus, either RSDs are common or a higher level of turbulence (b ~ 0:8 km s-1) is required in the inner envelope compared with the outer part (0.4 km s-1). Conclusions. Multiple spatially and spectrally resolved molecular line observations can differentiate between the pseudo-disk and the RSD much better than continuum data. The continuum data give a better estimate of disk masses, whereas the disk sizes can be estimated from the spatially resolved molecular lines observations. The general observable trends are similar between the 2D semi-analytical models and 3D MHD RSD simulations. [ABSTRACT FROM AUTHOR]

Context. Simple molecules such as H2CO and CH3OH in protoplanetary disks are the starting point for the production of more complex organic molecules. So far, the observed chemical complexity in disks has been limited because of freeze-out of molecules onto grains in the bulk of the cold outer disk. Aims. Complex molecules can be studied more directly in transitional disks with large inner holes because these have a higher potential of detection through the UV heating of the outer disk and the directly exposed midplane at the wall. Methods. We used Atacama Large Millimeter/submillimeter Array (ALMA) Band 9 (~680 GHz) line data of the transitional disk Oph IRS 48, which was previously shown to have a large dust trap, to search for complex molecules in regions where planetesimals are forming. Results. We report the detection of the H2CO 9(1, 8)-8(1, 7) line at 674 GHz, which is spatially resolved as a semi-ring at ~60 AU radius centered south from the star. The inferred H2CO abundance is ~10-8, derived by combining a physical disk model of the source with a non-LTE excitation calculation. Upper limits for CH3OH lines in the same disk give an abundance ratio H2CO/CH3OH >0.3, which indicates that both ice formation and gas-phase routes play a role in the H2CO production. Upper limits on the abundances of H13CO+, CN and several other molecules in the disk were also derived and found to be consistent with full chemical models. Conclusions. The detection of the H2CO line demonstrates the start of complex organic molecules in a planet-forming disk. Future ALMA observations are expected to reduce the abundance detection limits of other molecules by 1-2 orders of magnitude and test chemical models of organic molecules in (transitional) disks. [ABSTRACT FROM AUTHOR]

Context. For stars to form, material must fall inwards from core scales through the envelope towards the central protostar. While theories of how this takes place have been around for some time, the velocity profile around protostars is poorly constrained. The combination of observations in multiple transitions of a tracer which is sensitive to kinematics and radiative transfer modelling of those lines has the potential to break this deadlock. Aims. Seven protostars observed with the Heterodyne Instrument for the Far-Infrared (HIFI) on board the Herschel Space Observatory as part of the "Water in star-forming regions with Herschel" (WISH) survey show infall signatures in water line observations. We aim to constrain the infall velocity and the radii over which infall is taking place within the protostellar envelopes of these sources. We will also use these data to constrain the chemistry of cold water. Methods. We use 1-D non-LTE ratran radiative transfer models of the observed water lines to constrain the infall velocity and chemistry in the protostellar envelopes of six Class 0 protostars and one Class I source. We assume a free-fall velocity profile and, having found the best fit, vary the radii over which infall takes place. Results. In the well-studied Class 0 protostar NGC1333-IRAS4A we find that our observations probe infall over the whole envelope to which our observations are sensitive (r ≳ 1000 AU). For L1527, L1157, BHR71 and IRAS 15398 infall takes place on core to envelope scales (i.e. ~10 000-3000 AU). In Serpens-SMM4 and GSS30 the inverse P-Cygni profiles seen in the ground-state lines are more likely due to larger-scale motions or foreground clouds. Models including a simple consideration of the chemistry are consistent with the observations, while using step abundance profiles are not. The non-detection of excited water in the inner envelope in six out of seven protostars is further evidence that water must be heavily depleted from the gas-phase at these radii. Conclusions. Infall in four of the sources is supersonic and in all sources must take place at the outer edge of the envelope, which may be evidence that collapse is global or outside-in rather than inside-out. The mass infall rate in NGC1333-IRAS4A is large (≳10-4 M☉ yr-1), higher than the mass outflow rate and expected mass accretion rates onto the star. This suggests that any flattened disk-like structure on small scales will be gravitationally unstable, potentially leading to rotational fragmentation and/or episodic accretion. [ABSTRACT FROM AUTHOR]

Aims. Herschel-HIFI spectra of H2O towards low-mass protostars show a distinct velocity component not seen in observations from the ground of CO or other species. The aim is to characterise this component in terms of excitation conditions and physical origin. Methods. A velocity component with an offset of ~10 km s-1 detected in spectra of the H2O 110-101 557 GHz transition towards six low-mass protostars in the "Water in star-forming regions with Herschel" (WISH) programme is also seen in higher-excited H2O lines. The emission from this component is quantified and local excitation conditions are inferred using 1D slab models. Data are compared to observations of hydrides (high-J CO, OH+, CH+, C+, OH) where the same component is uniquely detected. Results. The velocity component is detected in all six targeted H2O transitions (Eup ~ 50-250 K), as well as in CO 16-15 towards one source, Ser SMM1. Inferred excitation conditions imply that the emission arises in dense (n ~ 5 × 106-108 cm-3) and hot (T ~ 750 K) gas. The H2O and CO column densities are ≳1016 and 1018 cm-2, respectively, implying a low H2O abundance of ~10-2 with respect to CO. The high column densities of ions such as OH+ and CH+ (both ≳1013 cm-2) indicate an origin close to the protostar where the UV field is strong enough that these species are abundant. The estimated radius of the emitting region is 100 AU. This component likely arises in dissociative shocks close to the protostar, an interpretation corroborated by a comparison with models of such shocks. Furthermore, one of the sources, IRAS 4A, shows temporal variability in the offset component over a period of two years which is expected from shocks in dense media. High-J CO gas detected with Herschel-PACS with Trot ~ 700 K is identified as arising in the same component and traces the part of the shock where H2 reforms. Thus, H2O reveals new dynamical components, even on small spatial scales in low-mass protostars. [ABSTRACT FROM AUTHOR]

Context. The OH radical is a key species in the water chemistry network of star-forming regions, because its presence is tightly related to the formation and destruction of water. Previous studies of the OH far-infrared emission from low- and intermediate-mass protostars suggest that the OH emission mainly originates from shocked gas and not from the quiescent protostellar envelopes. Aims. We aim to study the excitation of OH in embedded low- and intermediate-mass protostars, determine the influence of source parameters on the strength of the emission, investigate the spatial extent of the OH emission, and further constrain its origin. Methods. This paper presents OH observations from 23 low- and intermediate-mass young stellar objects obtained with the PACS integral field spectrometer on-board Herschel in the context of the "Water In Star-forming regions with Herschel" (WISH) key program. Radiative transfer codes are used to model the OH excitation. Results. Most low-mass sources have compact OH emission (≲5000 AU scale), whereas the OH lines in most intermediate-mass sources are extended over the whole 47."0 × 47."0 PACS detector field-of-view (≳20 000 AU). The strength of the OH emission is correlated with various source properties such as the bolometric luminosity and the envelope mass, but also with the [OI] and H2O emission. Rotational diagrams for sources with many OH lines show that the level populations of OH can be approximated by a Boltzmann distribution with an excitation temperature at around 70 K. Radiative transfer models of spherically symmetric envelopes cannot reproduce the OH emission fluxes nor their broad line widths, strongly suggesting an outflow origin. Slab excitation models indicate that the observed excitation temperature can either be reached if the OH molecules are exposed to a strong far-infrared continuum radiation field or if the gas temperature and density are sufficiently high. Using realistic source parameters and radiation fields, it is shown for the case of Ser SMM1 that radiative pumping plays an important role in transitions arising from upper level energies higher than 300 K. The compact emission in the low-mass sources and the required presence of a strong radiation field and/or a high density to excite the OH molecules points toward an origin in shocks in the inner envelope close to the protostar. [ABSTRACT FROM AUTHOR]

Context. Water is a key tracer of dynamics and chemistry in low-mass star-forming regions, but spectrally resolved observations have so far been limited in sensitivity and angular resolution, and only data from the brightest low-mass protostars have been published. Aims. The first systematic survey of spectrally resolved water emission in 29 low-mass (L <40 Lʘ) protostellar objects is presented. The sources cover a range of luminosities and evolutionary states. The aim is to characterise the line profiles to distinguish physical components in the beam and examine how water emission changes with protostellar evolution. Methods. H2O was observed in the ground-slate 110-101 transition at 557 GHz (Eup/kB -∼ 60K) as single-point observations with the Heterodyne Instrument for the Far-Infrared (HIFI) on Herschel in 29 deeply embedded Class 0 and I low-mass protostars. Complementary far-R and tub-mm continuum data (including PACS data from our programme) are used to constrain the spectral energy distribution (SED) of each source. H20 intensities are compared to inferred envelope properties, e.g., mass and density, outflow properties and CO 3-2 emission. Results. H2O emission is detected in all objects except one (TMCIA). The line profiles are complex and consist of several kinematic components tracing different physical regions in each system. In particular, the profiles are typically dominated by a broad Gaussian emission feature, indicating that the bulk of the water emission arises in outflows, not in the quiescent envelope. Several sources show multiple shock components appearing in either emission or absorption, thus constraining the internal geometry of the system. Furthermore, the components include inverse P-Cygni profiles in seven sources (nix Class 0, one Class 1) indicative of infalling envelopes, and regular P-Cygni profiles in four sources (three Class 1, one Class 0) indicative of expanding envelopes. Molecular "bullets" moving at ≳50 km s-1 with respect to the source are detected in four Class 0 sources; three of these sources were not known to harbour bullets previously. In the outflow, the H2O/CO abundance ratio an a function of velocity in nearly the name for all line wings, increasing from 10-3 at low velocities (<5 km s-1) to ≳ 10-1 at high velocities (>10 km s-1). The water abundance in the outer cold envelope in low, ≳10-10. The different H2O profile components show a clear evolutionary trend: in the younger Class 0 sources the emission in dominated by outflow components originating inside an infalling envelope. When large-scale infall diminishes during the Class I phase, the outflow weakens and H2O emission all but disappears. [ABSTRACT FROM AUTHOR]

Context. With the Herschel Space Observatory, lines of simple molecules (C+, O, and high-J lines of CO, Jup ≳ 14) have been observed in the atmosphere of protoplanetary disks. When combined with ground-based data on [C i], all principle forms of carbon can be studied. These data allow us to test model predictions for the main carbon-bearing species and verify the presence of a warm surface layer. The absence of neutral carbon [Ci], which is predicted by models to be strong, can then be interpreted together with ionized carbon [C ii] and carbon monoxide. Aims. We study the gas temperature, excitation, and chemical abundance of the simple carbon-bearing species C, C+, and CO, as well as O by the method of chemical-physical modeling. Using the models, we explore the sensitivity of the lines to the entering parameters and constrain the region from which the line radiation emerges. Methods. Numerical models of the radiative transfer in the lines and dust are used together with a chemical network simulation and a calculation of the gas energetics to obtain the gas temperature. We present our new model, which is based on our previous models but includes several improvements that we report in detail, together with the results of benchmark tests. Results. A model of the disk around the Herbig Be star HD 100546 is able to reproduce the CO ladder together with the atomic finestructure lines of [O i] and either [C i] or [Cii].We find that the high-J lines of CO can only be reproduced by a warm atmosphere with Tgas ⪢ Tdust. The low-J lines of CO, observable from the ground, are dominated by the outer disk with a radius of several 100 AU, while the high-J CO observable with Herschel-PACS are dominated from regions within some tens of AU. The spectral profiles of high-J lines of CO are predicted to be broader than those of the low-J lines. We study the effect of several parameters including the size of the disk, the gas mass of the disk, the PAH abundance and distribution, and the amount of carbon in the gas phase. Conclusions. The main conclusions of our work are (i) only a warm atmosphere with Tgas ⪢ Tdust can reproduce the CO ladder. (ii) The CO ladder together with [O i] and the upper limit to [Ci] can be reproduced by models with a high gas/dust ratio and a low abundance of volatile carbon. These models however produce too small amounts of [Cii]. Models with a low gas/dust ratio and more volatile carbon also reproduce CO and [O i], are in closer agreement with observations of [Cii], but overproduce [C i]. Owing to the uncertain origin of the [Cii] emission, we prefer the high gas/dust ratio models, indicating a low abundance of volatile carbon. [ABSTRACT FROM AUTHOR]

Context. Models of pure gas-phase chemistry in well-shielded regions of molecular clouds predict relatively high levels of molecular oxygen, O2, and water, H2O. These high abundances imply high cooling rates, leading to relatively short timescales for the evolution of gravitationally unstable dense cores, forming stars and planets. Contrary to expectations, the dedicated space missions SWAS and Odin typically found only very small amounts of water vapour and essentially no O2 in the dense star-forming interstellar medium. Aims. Only toward ρ OphA did Odin detect a very weak line of O2 at 119 GHz in a beam of size 10 arcmin. The line emission of related molecules changes on angular scales of the order of some tens of arcseconds, requiring a larger telescope aperture such as that of the Herschel Space Observatory to resolve the O2 emission and pinpoint its origin. Methods. We use the Heterodyne Instrument for the Far Infrared (HIFI) aboard Herschel to obtain high resolution O2 spectra toward selected positions in the ρOphAcore. These data are analysed using standard techniques for O2 excitation and compared to recent PDR-like chemical cloud models. Results. The NJ = 33-12 line at 487.2 GHz is clearly detected toward all three observed positions in the ρ Oph Acore. In addition, an oversampled map of the 54-34 transition at 773.8 GHz reveals the detection of the line in only half of the observed area. On the basis of their ratios, the temperature of the O2 emitting gas appears to vary quite substantially, with warm gas (≳50 K) being adjacent to a much colder region, of temperatures lower than 30 K. Conclusions. The exploited models predict that the O2 column densities are sensitive to the prevailing dust temperatures, but rather insensitive to the temperatures of the gas. In agreement with these models, the observationally determined O2 column densities do not seem to depend strongly on the derived gas temperatures, but fall into the range N(O2) = 3 to ≳6 ×1015 cm-2. Beam-averaged O2 abundances are about 5 ×10-8 relative to H2. Combining the HIFI data with earlier Odin observations yields a source size at 119 GHz in the range of 4 to 5 arcmin, encompassing the entire ρ OphAcore. We speculate that one of the reasons for the generally very low detection rate of O2 is the short period of time during which O2 molecules are reasonably abundant in molecular clouds. [ABSTRACT FROM AUTHOR]

Aims. Young stars interact vigorously with their surroundings, as evident from the highly rotationally excited CO (up to Eu/k = 4000 K) and H2O emission (up to 600 K) detected by the Herschel Space Observatory in embedded low-mass protostars. Our aim is to construct a model that reproduces the observations quantitatively, to investigate the origin of the emission, and to use the lines as probes of the various heating mechanisms. Methods. The model consists of a spherical envelope with a power-law density structure and a bipolar outflow cavity. Three heating mechanisms are considered: passive heating by the protostellar luminosity, ultraviolet irradiation of the outflow cavity walls, and small-scale C-type shocks along the cavity walls. Most of the model parameters are constrained from independent observations; the two remaining free parameters considered here are the protostellar UV luminosity and the shock velocity. Line fluxes are calculated for CO and H2O and compared to Herschel data and complementary ground-based data for the protostars NGC 1333 IRAS2A, HH 46 and DK Cha. The three sources are selected to span a range of evolutionary phases (early Stage 0 to late Stage I) and physical characteristics such as luminosity and envelope mass. Results. The bulk of the gas in the envelope, heated by the protostellar luminosity, accounts for 3-10% of the CO luminosity summed over all rotational lines up to J = 40-39; it is best probed by low-J CO isotopologue lines such as C18O 2-1 and 3-2. The UV-heated gas and the C-type shocks, probed by 12CO 10-9 and higher-J lines, contribute 20-80% each. The model fits show a tentative evolutionary trend: the CO emission is dominated by shocks in the youngest source and by UV-heated gas in the oldest one. This trend is mainly driven by the lower envelope density in more evolved sources. The total H2O line luminosity in all cases is dominated by shocks (>99%). The exact percentages for both species are uncertain by at least a factor of 2 due to uncertainties in the gas temperature as function of the incident UV flux. However, on a qualitative level and within the context of our model, both UV-heated gas and C-type shocks are needed to reproduce the emission in far-infrared rotational lines of CO and H2O. [ABSTRACT FROM AUTHOR]