To understand the formation of planetary systems, one needs to understand the initial conditions of planet formation, i.e. the young gas-rich planet forming disks. Spatially resolved high-contrast observations are of particular interest, since substructures in disks, linked to planet formation, can be detected and close companions or even planets in formation embedded in the disk can be revealed. In this study we present the first result of the DESTINYS survey (Disk Evolution Study Through Imaging of Nearby Young Stars). DESTINYS is an ESO/SPHERE large program that aims at studying disk evolution in scattered light, mainly focusing on a sample of low-mass stars (<1$M_\odot$) in nearby (~200 pc) star-forming regions. In this particular study we present the observations of the ET Cha (RECX 15) system, a nearby 'old' classical T Tauri star (5-8 Myr, ~100 pc), which is still strongly accreting. We use SPHERE/IRDIS in H-band polarimetric imaging mode to obtain high contrast images of the ET Cha system to search for scattered light from the circumstellar disk as well as thermal emission from close companions. We additionally employ VLT/NACO total intensity archival data taken in 2003. We report here the discovery of a low-mass (sub)stellar companion with SPHERE/IRDIS to ET Cha. We are estimating the mass of this new companion based on photometry. Depending on the system age it is a 5 Myr, 50 $M_{Jup}$ brown dwarf or an 8 Myr, 0.10 $M_\odot$ M-type pre-main-sequence star. We explore possible orbital solutions and discuss the recent dynamic history of the system. Independent of the precise companion mass we find that the presence of the companion likely explains the small size of the disk around ET Cha. The small separation of the binary pair indicates that the disk around the primary component is likely clearing from the outside in, explaining the high accretion rate of the system.

The evolution of protoplanetary disks is dominated by the conservation of angular momentum, where the accretion of material onto the central star is driven by viscous expansion of the outer disk or by disk winds extracting angular momentum without changing the disk size. Studying the time evolution of disk sizes allows us therefore to distinguish between viscous stresses or disk winds as the main mechanism of disk evolution. Observationally, estimates of the disk gaseous outer radius are based on the extent of the CO rotational emission, which, during the evolution, is also affected by the changing physical and chemical conditions in the disk. We use physical-chemical DALI models to study how the extent of the CO emission changes with time in a viscously expanding disk and investigate to what degree this observable gas outer radius is a suitable tracer of viscous spreading and whether current observations are consistent with viscous evolution. We find that the gas outer radius (R_co) measured from our models matches the expectations of a viscously spreading disk: R_co increases with time and for a given time R_co is larger for a disk with a higher viscosity alpha_visc. However, in the extreme case where the disk mass is low (less than 10^-4 Msun) and alpha_visc is high (larger than 10^-2), R_co will instead decrease with time as a result of CO photodissociation in the outer disk. For most disk ages R_co is up to 12x larger than the characteristic size R_c of the disk, and R_co/R_c is largest for the most massive disk. As a result of this difference, a simple conversion of R_co to alpha_visc will overestimate the true alpha_visc of the disk by up to an order of magnitude. We find that most observed gas outer radii in Lupus can be explained using a viscously evolving disk that starts out small (R_c = 10 AU) and has a low viscosity (alpha_visc = 10^-4 - 10^-3).

Recent ALMA surveys of protoplanetary disks have shown that for most disks the extent of the gas emission is greater than the extent of the thermal emission of the millimeter-sized dust. Both line optical depth and the combined effect of radially dependent grain growth and radial drift may contribute to this observed effect. For a sample of 10 disks from the Lupus survey we investigate how well dust-based models without radial dust evolution reproduce the observed 12CO outer radius, and determine whether radial dust evolution is required to match the observed gas-dust size difference. We used the thermochemical code DALI to obtain 12CO synthetic emission maps and measure gas and dust outer radii (Rco, Rmm) using the same methods as applied to the observations, which were compared to observations on a source-by-source basis. For 5 disks we find that the observed gas-dust size difference is larger than the gas-dust size difference due to optical depth, indicating that we need both dust evolution and optical depth effects to explain the observed gas-dust size difference. For the other 5 disks the observed gas-dust size difference can be explained using only line optical depth effects. We also identify 6 disks not included in our initial sample but part of a survey of the same star-forming region that show significant 12CO emission beyond 4 x Rmm. These disks, for which no Rco is available, likely have gas-dust size differences greater than 4 and are difficult to explain without substantial dust evolution. Our results suggest that radial drift and grain growth are common features among both bright and fain disks. The effects of radial drift and grain growth can be observed in disks where the dust and gas radii are significantly different, while more detailed models and deeper observations are needed to see this effect in disks with smaller differences.

Water is a key volatile that provides insights into the initial stages of planet formation. The low water abundances inferred from water observations toward low-mass protostellar objects may point to a rapid locking of water as ice by large dust grains during star and planet formation. However, little is known about the water vapor abundance in newly formed planet-forming disks. We aim to determine the water abundance in embedded Keplerian disks through spatially-resolved observations of H$_2^{18}$O lines to understand the evolution of water during star and planet formation. We present H$_2^{18}$O line observations with ALMA and NOEMA millimeter interferometers toward five young stellar objects. NOEMA observed the 3$_{1,3}$ - $2_{2,0}$ line (E$_{\rm up}$ = 203.7 K) while ALMA targeted the $4_{1,4}$ - $3_{2,1}$ line (E$_{\rm up}$ = 322.0 K). Water column densities are derived considering optically thin and thermalized emission. Our observations are sensitive to the emission from the known Keplerian disks around three out of the five Class I objects in the sample. No H$_2^{18}$O emission is detected toward any of our five Class I disks. We report upper limits to the integrated line intensities. The inferred water column densities in Class I disks are N < 10$^{15}$ cm$^{-2}$ on 100 au scales which include both disk and envelope. The upper limits imply a disk-averaged water abundance of $\lesssim 10^{-6}$ with respect to H$_2$ for Class I objects. After taking into account the physical structure of the disk, the upper limit to the water abundance averaged over the inner warm disk with $T>$ 100 K is between 10$^{-7}$ up to 10$^{-5}$. Water vapor is not abundant in warm protostellar envelopes around Class I protostars. Upper limits to the water vapor column densities in Class I disks are at least two orders magnitude lower than values found in Class 0 disk-like structures.

Hydrogen deuteride (HD) rotational line emission can provide reliable protoplanetary disk gas mass measurements, but it is difficult to observe and detections have been limited to three T-Tauri disks. No new data have been available since the \emph{Herschel} Space Observatory mission ended in 2013. We set out to obtain new disk gas mass constraints by analysing upper limits on HD 1-0 emission in \emph{Herschel}/PACS archival data from the DIGIT key programme. With a focus on the Herbig Ae/Be disks, whose stars are more luminous than T Tauris, we determine upper limits for HD in data previosly analysed for its line detections. Their significance is studied with a grid of models run with the DALI physical-chemical code, customised to include deuterium chemistry. Nearly all the disks are constrained to $M_{\rm gas}\leq0.1\,$M$_{\odot}$, ruling out global gravitational instability. A strong constraint is obtained for the HD 163296 disk mass, $M_{\rm gas} \leq 0.067\,$M$_{\odot}$, implying $\Delta_{\rm g/d}\leq100$. This HD-based mass limit is towards the low end of CO-based mass estimates for the disk, highlighting the large uncertainty in using only CO and suggesting that gas-phase CO depletion in HD 163296 is at most a factor of a few. The $M_{\rm gas}$ limits for HD 163296 and HD 100546, both bright disks with massive candidate protoplanetary systems, suggest disk-to-planet mass conversion efficiencies of $M_{\rm p}/(M_{\rm gas} + M_{\rm p})\approx10$ to $40\,$% for present-day values. Near-future observations with SOFIA/HIRMES will be able to detect HD in the brightest Herbig~Ae/Be disks within $150\,$pc with $\approx10\,$h integration time.

We present ALMA CO(1-0) observations toward the dust lane of the nearest elliptical and radio galaxy, NGC 5128 (Centaurus A), with high angular resolution ($\sim$ 1 arcsec, or 18 pc), including information from large to small spatial scales and total flux. We find a total molecular gas mass of 1.6$\times$10$^9$ $M_\odot$ and we reveal the presence of filamentary components more extended than previously seen, up to a radius of 4 kpc. We find that the global star formation rate is $\sim$1 \Msol yr$^{-1}$, which yields a star formation efficiency (SFE) of 0.6 Gyr$^{-1}$ (depletion time $\tau =$1.5 Gyr), similar to those in disk galaxies. We show the most detailed view to date (40\,pc resolution) of the relation between molecular gas and star formation within the stellar component of an elliptical galaxy, from several kpc scale to the circumnuclear region close to the powerful radio jet. Although on average the SFEs are similar to those of spiral galaxies, the circumnuclear disk (CND) presents SFEs of 0.3 Gyr$^{-1}$, lower by a factor of 4 than the outer disk. The low SFE in the CND is in contrast to the high SFEs found in the literature for the circumnuclear regions of some nearby disk galaxies with nuclear activity, probably as a result of larger shear motions and longer AGN feedback. The higher SFEs in the outer disk suggests that only central molecular gas or filaments with sufficient density and strong shear motions will remain in $\sim$1 Gyr, which will later result in the compact molecular distributions and low SFEs usually seen in other giant ellipticals with cold gas.

We present ALMA observations of 12CO, 13CO, and C18O J=2--1 lines and the 230 GHz continuum for the FU Ori-type object (FUor) V900 Mon (d~1.5 kpc), for which the accretion burst was triggered between 1953 and 2009. We identified CO emission associated with a molecular bipolar outflow extending up to a ~10^4 au scale and a rotating molecular envelope extending over >10^4 au. The interaction with the hot energetic FUor wind, which was observed using optical spectroscopy, appears limited to a region within ~400 au of the star. The envelope mass and the collimation of the extended CO outflow suggest that the progenitor of this FUor is a low-mass Class I young stellar object (YSO). These parameters for V900 Mon, another FUor, and a few FUor-like stars are consistent with the idea that FUor outbursts are associated with normal YSOs. The continuum emission is marginally resolved in our observations with a 0."2x0."15 (~300x225 au) beam, and a Gaussian model provides a deconvolved FWHM of ~90 au. The emission is presumably associated with a dusty circumstellar disk, plus a possible contribution from a wind or a wind cavity close to the star. The warm compact nature of the disk continuum emission could be explained with viscous heating of the disk, while gravitational fragmentation in the outer disk and/or a combination of grain growth and their inward drift may also contribute to its compact nature.

Observations of young stellar objects (YSOs) in centimeter bands can probe the continuum emission from growing dust grains, ionized winds, and magnetospheric activity, which are intimately connected to the evolution of protoplanetary disks and the formation of planets. We have carried out sensitive continuum observations toward the Ophiuchus A star-forming region using the Karl G. Jansky Very Large Array (VLA) at 10 GHz over a field-of-view of 6$'$ with a spatial resolution of $\theta_{maj}$ $\times$ $\theta_{min}$ $\sim$ 0.4$''$ $\times$ 0.2$''$. We achieved a 5 $\mu$Jy beam$^{-1}$ root-mean-square noise level at the center of our mosaic field of view. Among the eighteen sources we detected, sixteen are YSOs (three Class 0, five Class I, six Class II, and two Class III) and two are extragalactic candidates. We find that thermal dust emission generally contributes less that 30% of the emission at 10 GHz. The radio emission is dominated by other types of emission such as gyro-synchrotron radiation from active magnetospheres, free-free emission from thermal jets, free-free emission from the outflowing photo-evaporated disk material, and/or synchrotron emission from accelerated cosmic-rays in jet or protostellar surface shocks. These different types of emission could not be clearly disentangled. Our non-detections towards Class II/III disks suggest that extreme UV-driven photoevaporation is insufficient to explain the disk dispersal, assuming that the contribution of UV photoevaporating stellar winds to radio flux does not evolve with time. The sensitivity of our data cannot exclude photoevaporation due to X-ray photons as an efficient mechanism for disk dispersal. Deeper surveys with the Square Kilometre Array will be able to provide strong constraints on disk photoevaporation.

We present the highest angular resolution (20x15mas - 44x33au) Atacama Large Millimeter/sub-millimeter Array (ALMA) observations currently possible of the proto-O-star G17.64+0.16 in Band 6. The Cycle 5 observations with baselines out to 16km probes scales <50au and reveal the rotating disc around G17.64+0.16, a massive forming O-type star. The disc has a ring-like enhancement in the dust emission, especially visible as arc structures to the north and south. The Keplerian kinematics are most prominently seen in the vibrationally excited water line, H2O (Eu=3461.9K). The mass of the central source found by modelling the Keplerian rotation is consistent with 45+/-10Mo. The H30alpha (231.9GHz) radio-recombination line and the SiO (5-4) molecular line were detected at up to the 10 sigma$ level. The estimated disc mass is 0.6-2.6Mo under the optically thin assumption. Analysis of the Toomre Q parameter, in the optically thin regime, indicates that the disc stability is highly dependent on temperature. The disc currently appears stable for temperatures >150K, this does not preclude that the substructures formed earlier through disc fragmentation.

The chemical composition of high-mass protostars reflects the physical evolution associated with different stages of star formation. In this study, the molecular inventory of the forming high-mass star AFGL 4176 is studied in detail at high angular resolution (~0.35 arcsec) using ALMA. This high resolution makes it possible to separate the emission associated with the inner hot envelope and disk around the forming star from that of its cool outer envelope. In addition, the high sensitivity of ALMA makes it possible to identify weak and optically thin lines and allows for many isotopologues to be detected, providing a more complete and accurate inventory of the source. A total of 23 different molecular species and their isotopologues are detected in the spectrum towards AFGL 4176. The most abundant species is methanol (CH3OH), remaining species are present at levels between 0.003 % and 15 % with respect to CH3OH. Hints that N-bearing species peak slightly closer to the location of the peak continuum emission than the O-bearing species are seen. AFGL 4176 comprises a rich chemical inventory including many complex species present on disk-scales. On average, the derived column density ratios with respect to methanol of O-bearing species are higher than those derived for N-bearing species by a factor of three. This may indicate that AFGL 4176 is a relatively young source since nitrogen chemistry generally takes longer to evolve in the gas-phase. Taking methanol as a reference, the composition of AFGL 4176 more closely resembles that of the low-mass protostar IRAS 16293-2422B than that of high-mass star-forming regions located near the Galactic centre. This similarity hints that the chemical composition of complex species is already set in the cold cloud stage and implies that AFGL 4176 is a young source whose chemical composition has not yet been strongly processed by the central protostar.

In recent years, the disk populations in a number of young star-forming regions have been surveyed with ALMA. Understanding the disk properties and their correlation with those of the central star is critical to understand planet formation. In particular, a decrease of the average measured disk dust mass with the age of the region has been observed. We conducted high-sensitivity continuum ALMA observations of 43 Class II young stellar objects in CrA at 1.3 mm (230 GHz). The typical spatial resolution is 0.3". The continuum fluxes are used to estimate the dust masses of the disks, and a survival analysis is performed to estimate the average dust mass. We also obtained new VLT/X-Shooter spectra for 12 of the objects in our sample. 24 disks are detected, and stringent limits have been put on the average dust mass of the non-detections. Accounting for the upper limits, the average disk mass in CrA is $6\pm3\,\rm M_\oplus$, significantly lower than that of disks in other young (1-3 Myr) star forming regions (e.g. Lupus) and appears consistent with the 5-10 Myr old Upper Sco. The position of the stars in our sample on the HR diagram, however, seems to confirm that that CrA has age similar to Lupus. Neither external photoevaporation nor a lower than usual stellar mass distribution can explain the low disk masses. On the other hand, a low-mass disk population could be explained if the disks are small, which could happen if the parent cloud has a low temperature or intrinsic angular momentum, or if the the angular momentum of the cloud is removed by some physical mechanism such as magnetic braking. In order to fully explain and understand the dust mass distribution of protoplanetary disks and their evolution, it may also be necessary to take into consideration the initial conditions of star and disk formation process, which may vary from region to region, and affect planet formation.

Understanding the origin of the astonishing diversity of exoplanets is a key question for the coming decades. ALMA has revolutionized our view of the dust emission from protoplanetary disks, demonstrating the prevalence of ring and spiral structures that are likely sculpted by young planets in formation. To detect kinematic signatures of these protoplanets and to probe the chemistry of their gas accretion reservoir will require the imaging of molecular spectral line emission at high angular and spectral resolution. However, the current sensitivity of ALMA limits these important spectral studies to only the nearest protoplanetary disks. Although some promising results are emerging, including the identification of the snowlines of a few key molecules and the first attempt at detecting a protoplanet's spiral wake, it is not yet possible to search for these important signatures in a population of disks in diverse environments and ages. Harnessing the tremendous power of (sub)mm observations to pinpoint and characterize the chemistry of planets in formation will require a major increase of ALMA's spectral sensitivity (5-10x), increase in instantaneous bandwidth (2x) at high spectral resolution, and improved angular resolution (2x) in the 2030 era.

The extent of the gas in protoplanetary disks is observed to be universally larger than the extent of the dust. This is often attributed to radial drift and grain growth of the mm grains, but line optical depth produces a similar observational signature. We investigate in what parts of the disk structure parameter space dust evolution and line optical depth are the dominant drivers of the observed gas and dust size difference. Using the thermochemical model DALI with dust evolution included we ran a grid of models aimed at reproducing the observed gas and dust size dichotomy. The relation between R_dust and dust evolution is non-monotonic and depends on the disk structure. R_gas is directly related to the radius where the CO column density drops below 10^15 cm^-2 and CO becomes photodissociated. R_gas is not affected by dust evolution but scales with the total CO content of the disk. R_gas/R_dust > 4 is a clear sign for dust evolution and radial drift in disks, but these cases are rare in current observations. For disks with a smaller R_gas/R_dust, identifying dust evolution from R_gas/R_dust requires modelling the disk structure including the total CO content. To minimize the uncertainties due to observational factors requires FWHM_beam < 1x the characteristic radius and a peak SNR > 10 on the 12CO emission moment zero map. For the dust outer radius to enclose most of the disk mass, it should be defined using a high fraction (90-95%) of the total flux. For the gas, any radius enclosing > 60% of the 12CO flux will contain most of the disk mass. To distinguish radial drift and grain growth from line optical depth effects based on size ratios requires disks to be observed at high enough angular resolution and the disk structure should to be modelled to account for the total CO content of the disk.

In the search for the building blocks of life, nitrogen-bearing molecules are of particular interest since nitrogen-containing bonds are essential for the linking of amino acids and ultimately the formation of larger biological structures. The elusive molecule methylamine (CH$_3$NH$_2$) is thought to be a key pre-biotic species but has so far only been securely detected in the giant molecular cloud Sgr B2. We identify CH$_3$NH$_2$ and other simple nitrogen-bearing species towards three hot cores in NGC 6334I. Column density ratios are derived in order to investigate the relevance of the individual species as precursors of biotic molecules. Observations obtained with ALMA were used to study transitions of CH$_3$NH$_2$, CH$_2$NH, NH$_2$CHO, and the $^{13}$C- and $^{15}$N-methyl cyanide (CH$_3$CN) isotopologues. Column densities are derived for each species assuming LTE and excitation temperatures in the range 220-340 K for CH$_3$NH$_2$, 70-110 K for the CH$_3$CN isotopologues, and 120-215 K for NH$_2$CHO and CH$_2$NH. We report the first detections of CH$_3$NH$_2$ towards NGC 6334I with column density ratios with respect to CH$_3$OH of 5.9$\times$10$^{-3}$, 1.5$\times$10$^{-3}$, and 5.4$\times$10$^{-4}$ for the three hot cores MM1, MM2, and MM3, respectively. These values are slightly lower than the values derived for Sgr B2 but higher by more than order of magnitude as compared with the values derived for the low-mass protostar IRAS 16293-2422B. The detections of CH$_3$NH$_2$ in the hot cores of NGC 6334I hint that CH$_3$NH$_2$ is generally common in the interstellar medium, albeit high-sensitivity observations are essential for its detection. The good agreement between model predictions of CH$_3$NH$_2$ ratios and the observations towards NGC 6334I indicate a main formation pathway via radical recombination on grain surfaces.

Methanol (CH$_3$OH) is at the root of organic ice chemistry in protoplanetary disks. However, its weak emission has made detections difficult. To date, gas-phase CH$_3$OH has been detected in only one Class II disk, TW Hya. We use the Atacama Large Millimeter/submillimeter Array (ALMA) to search for a total of four CH$_3$OH emission lines in bands 6 and 7 toward the disk around the young Herbig Ae star HD 163296. The disk-averaged column density of methanol and its related species formaldehyde (H$_2$CO) are estimated assuming optically thin emission in local thermodynamic equilibrium. We compare these results to the gas-phase column densities of the TW Hya disk. No targeted methanol lines were detected individually nor after line stacking. The 3$\sigma$ disk-integrated intensity upper limits are $< 51$ mJy km s$^{-1}$ for the band 6 lines and $< 26$ mJy km s$^{-1}$ for the band 7 lines. The band 7 lines provide the strictest 3$\sigma$ upper limit on disk-averaged column density with $N_{\mathrm{avg}} < 5.0 \times 10^{11}$ cm$^{-2}$. The methanol-to-formaldehyde ratio is CH$_3$OH/H$_2$CO $< 0.24$ in the HD 163296 disk compared to a ratio of $1.27$ in the TW Hya disk. Differences in the stellar irradiation of Herbig disks compared to T Tauri disks likely influence the gaseous methanol and formaldehyde content. Possible reasons for the lower HD 163296 methanol-to-formaldehyde ratio include: a higher than expected gas-phase formation of H$_2$CO in the HD 163296 disk, uncertainties in the grain surface formation efficiency of CH$_3$OH and H$_2$CO, and differences in the disk structure and/or CH$_3$OH and H$_2$CO desorption processes that release the molecules from ice mantles back into the gas phase. These results provide observational evidence that the gas-phase chemical complexity found in disks may be strongly influenced by the spectral type of the host star.

The cyanide radical CN is abundant in protoplanetary disks, with line fluxes often comparable to those of $^{13}$CO. It is known to be sensitive to UV irradiation of the upper disk atmosphere, with models predicting ring-shaped emission. We seek to characterize the CN emission from 94 Class II disks in the Lupus star forming region, compare it to observations in other regions, and interpret our observations with a grid of models. The CN emission morphology is discussed for two primordial disks, Sz 71 and Sz 98, and modeled in more detail. ALMA observed CN N=3-2 in Lupus disks down to sensitivities better than previous surveys. Models constructed with the physico-chemical code DALI are used to study the integrated fluxes of the disks and resolved emission of CN in disks without (dust) substructures. We find that CN N=3-2 is bright, and detected in 38% of sources, but its disk-integrated flux is not strongly correlated to either $^{13}$CO or continuum flux. Compared to pre-ALMA single-dish surveys, no significant difference in the CN flux distributions in Lupus and Taurus-Auriga is found, although $\rho$ Ophiuchus disks may be fainter on average. We find ring-shaped CN emission with peak radii of ~50 AU in two resolved disks. We conclude that a large fraction of sources is faint in CN; only exponential gas surface density cutoffs at $R_{c} \leq 15$ AU can reconcile observations with models. This is the first observational evidence of such a compact gas disk population in Lupus. Absolute intensities and the emission morphology of CN are reproduced by DALI models without the need for any continuum substructure; they are unrelated to the CO snowline location. These observations and the successful modeling of these rings provide a new probe of the structure and conditions in disks, and particularly their incident UV radiation field, if disk size is determined by the data.

We present high angular resolution 0.2 arcsec continuum and molecular emission line Atacama Large Millimeter/sub-millimeter Array (ALMA) observations of G17.64+0.16 in Band 6 (220GHz) taken as part of a campaign in search of circumstellar discs around (proto)-O-stars. At a resolution of 400au the main continuum core is essentially unresolved and isolated from other strong and compact emission peaks. At a resolution of 400au the main continuum core is essentially unresolved and isolated from other strong and compact emission peaks. We detect SiO (5-4) emission that is marginally resolved and elongated in a direction perpendicular to the large-scale outflow seen in the 13CO (2-1) line using the main ALMA array in conjunction with the Atacama Compact Array (ACA). Morphologically, the SiO appears to represent a disc-like structure. Using parametric models we show that the position-velocity profile of the SiO is consistent with the Keplerian rotation of a disc around an object between 10-30Mo in mass, only if there is also radial expansion from a separate structure. The radial motion component can be interpreted as a disc wind from the disc surface. Models with a central stellar object mass between 20 and 30Mo are the most consistent with the stellar luminosity (100000 Lo) and indicative of an O-type star. The H30a millimetre recombination line (231.9GHz) is also detected, but spatially unresolved, and is indicative of a very compact, hot, ionised region co-spatial with the dust continuum core. Accounting for all observables, we suggest that G17.64 is consistent with a O-type young stellar object in the final stages of protostellar assembly, driving a wind, but that has not yet developed into a compact HII region. The existance and detection of the disc in G17.64 is likely related to its isolated and possibly more evolved nature, traits which may underpin discs in similar sources.

VLA 1623$-$2417 is a triple protostellar system deeply embedded in Ophiuchus A. Sources A and B have a separation of 1.1", making their study difficult beyond the submillimeter regime. Lack of circumstellar gas emission suggested that VLA 1623$-$2417 B has a very cold envelope and is much younger than source A, generally considered the prototypical Class 0 source. We explore the consequences of new ALMA Band 9 data on the spectral energy distribution (SED) of VLA 1623$-$2417 and their inferred nature. Using dust continuum observations spanning from centimeter to near-infrared wavelengths, the SED of each component in VLA 1623$-$2417 is constructed and analysed. The ALMA Band 9 data presented here show that the SED of VLA 1623$-$2417 B does not peak at 850 $\mu$m as previously expected, but instead presents the same shape as VLA 1623$-$2417 A at wavelengths shorter than 450 $\mu$m. The results presented here indicate that the previous assumption that the flux in $Herschel$ and Spitzer observations is solely dominated by VLA 1623$-$2417 A is not valid, and instead, VLA 1623$-$2417 B most likely contributes a significant fraction of the flux at $\lambda~<$ 450 $\mu$m. These results, however, do not explain the lack of circumstellar gas emission and puzzling nature of VLA 1623$-$2417 B.

We aim to reproduce the DCO$^+$ emission in the disk around HD163296 using a simple 2D chemical model for the formation of DCO$^+$ through the cold deuteration channel and a parametric treatment of the warm deuteration channel. We use data from ALMA in band 6 to obtain a resolved spectral imaging data cube of the DCO$^+$ $J$=3--2 line in HD163296 with a synthesized beam of 0."53$\times$ 0."42. We adopt a physical structure of the disk from the literature that reproduces the spectral energy distribution. We then apply a simplified chemical network for the formation of DCO$^+$ that uses the physical structure of the disk as parameters along with a CO abundance profile, a constant HD abundance and a constant ionization rate. Finally, from the resulting DCO$^+$ abundances, we calculate the non-LTE emission using the 3D radiative transfer code LIME. The observed DCO$^+$ emission is reproduced by a model with cold deuteration producing abundances up to $1.6\times 10^{-11}$. Warm deuteration, at a constant abundance of $3.2\times 10^{-12}$, becomes fully effective below 32 K and tapers off at higher temperatures, reproducing the lack of DCO$^+$ inside 90 AU. Throughout the DCO$^+$ emitting zone a CO abundance of $2\times 10^{-7}$ is found, with $\sim$99\% of it frozen out below 19 K. At radii where both cold and warm deuteration are active, warm deuteration contributes up to 20\% of DCO$^+$, consistent with detailed chemical models. The decrease of DCO$^+$ at large radii is attributed to a temperature inversion at 250 AU, which raises temperatures above values where cold deuteration operates. Increased photodesorption may also limit the radial extent of DCO$^+$. The corresponding return of the DCO$^+$ layer to the midplane, together with a radially increasing ionization fraction, reproduces the local DCO$^+$ emission maximum at $\sim$260 AU.

A wide variety of ring-like dust structures has been detected in protoplanetary disks, but their origin and frequency are still unclear. We characterize the structure of an extended, multi-ringed disk discovered serendipitously in the ALMA Lupus disk survey and put it in the context of the Lupus disk population. ALMA observations in Band 6 at 234 GHz and Band 7 at 328 GHz at 0.3" resolution toward the K6 star V1094 Sco in Lupus III are presented, and its disk structure is analyzed. The spectral index $\alpha_{mm}$ is determined in the inner 150 AU of the disk. The ALMA continuum data show a very extended disk with two gap/ring pairs. The gaps are located at 100 AU and 170 AU, the bright rings at 130 AU and 220 AU. Continuum emission is detected out to a 300 AU distance, similar to IM Lup but a factor of 5 larger than typically found for Lupus disks at this sensitivity and resolution. The bright central region of the disk (within 35 AU) is possibly optically thick at 1 mm wavelengths, and has a brightness temperature of only 13 K. The spectral index increases between the inner disk and the first ring, at the location of the first gap. Due to the low temperature of the disk midplane, snow lines can be excluded as the drivers behind the ring and gap formation in this disk. Disks the size of V1094 Sco are rare, and only 2.1+-1.5% of disks in Lupus show continuum emission beyond 200 AU. Possible connections between the large primordial disk population, transition disks, and exoplanets are discussed.

We present ALMA Band 6 observations of a complete sample of protoplanetary disks in the young (1-3 Myr) Lupus star-forming region, covering the 1.33 mm continuum and the 12CO, 13CO, and C18O J=2-1 lines. The spatial resolution is 0.25 arcsec with a medium 3-sigma continuum sensitivity of 0.30 mJy, corresponding to M_dust ~ 0.2 M_earth. We apply "Keplerian masking" to enhance the signal-to-noise ratios of our 12CO zero-moment maps, enabling measurements of gas disk radii for 22 Lupus disks; we find that gas disks are universally larger than mm dust disks by a factor of two on average, likely due to a combination of the optically thick gas emission as well as the growth and inward drift of the dust. Using the gas disk radii, we calculate the dimensionless viscosity parameter, alpha_visc, finding a broad distribution and no correlations with other disk or stellar parameters, suggesting that viscous processes have not yet established quasi-steady states in Lupus disks. By combining our 1.33 mm continuum fluxes with our previous 890 micron continuum observations, we also calculate the mm spectral index, alpha_mm, for 70 Lupus disks; we find an anti-correlation between alpha_mm and mm flux for low-mass disks (M_dust < 5), followed by a flattening as disks approach alpha_mm = 2, which could indicate faster grain growth in higher-mass disks, but may also reflect their larger optically thick components. In sum, this work demonstrates the continuous stream of new insights into disk evolution and planet formation that can be gleaned from unbiased ALMA disk surveys.

The abundance of deuterated molecules in a star-forming region is sensitive to the environment in which they are formed. Deuteration fractions therefore provide a powerful tool for studying the physical and chemical evolution of a star-forming system. While local low-mass star-forming regions show very high deuteration ratios, much lower fractions are observed towards Orion and the Galactic Centre. We derive methanol deuteration fractions at a number of locations towards the high-mass star-forming region NGC 6334I, located at a mean distance of 1.3 kpc, and discuss how these can shed light on the conditions prevailing during its formation. We use high sensitivity, high spatial and spectral resolution observations obtained with ALMA to study transitions of the less abundant, optically thin, methanol-isotopologues: (13)CH3OH, CH3(18)OH, CH2DOH and CH3OD, detected towards NGC 6334I. Assuming LTE and excitation temperatures of 120-330 K, we derive column densities for each of the species and use these to infer CH2DOH/CH3OH and CH3OD/CH3OH fractions. Interestingly, the column densities of CH3OD are consistently higher than those of CH2DOH throughout the region. All regions studied in this work show CH2DOH/CH3OH as well as CH2DOH/CH3OD ratios that are considerably lower than those derived towards low-mass star-forming regions and slightly lower than those derived for the high-mass star-forming regions in Orion and the Galactic Centre. The low ratios indicate a grain surface temperature during formation ~30 K, for which the efficiency of the formation of deuterated species is significantly reduced.

This work aims to understand which midplane conditions are probed by the DCO$^+$ emission in the disk around the Herbig Ae star HD 169142. We explore the sensitivity of the DCO$^+$ formation pathways to the gas temperature and the CO abundance. The DCO$^+$ $J$=3-2 transition was observed with ALMA at a spatial resolution of 0.3". The HD 169142 DCO$^+$ radial intensity profile reveals a warm, inner component at radii <30 AU and a broad, ring-like structure from ~50-230 AU with a peak at 100 AU just beyond the millimeter grain edge. We modeled DCO$^+$ emission in HD 169142 with a physical disk structure adapted from the literature, and employed a simple deuterium chemical network to investigate the formation of DCO$^+$ through the cold deuterium fractionation pathway via H$_2$D$^+$. Contributions from the warm deuterium fractionation pathway via CH$_2$D$^+$ are approximated using a constant abundance in the intermediate disk layers. Parameterized models show that alterations to the midplane gas temperature and CO abundance of the literature model are both needed to recover the observed DCO$^+$ radial intensity profile. The best-fit model contains a shadowed, cold midplane in the region z/r < 0.1 with an 8 K decrease in gas temperature and a factor of five CO depletion just beyond the millimeter grain edge, and a 2 K decrease in gas temperature for r > 120 AU. The warm deuterium fractionation pathway is implemented as a constant DCO$^+$ abundance of 2.0$\times$10$^{-12}$ between 30-70 K. The DCO$^+$ emission probes a reservoir of cold material in the HD 169142 outer disk that is not revealed by the millimeter continuum, the SED, nor the emission from the 12CO, 13CO, or C18O $J$=2-1 lines.

We present 850 and 450 micron observations of the dense regions within the Auriga-California molecular cloud using SCUBA-2 as part of the JCMT Gould Belt Legacy Survey to identify candidate protostellar objects, measure the masses of their circumstellar material (disk and envelope), and compare the star formation to that in the Orion A molecular cloud. We identify 59 candidate protostars based on the presence of compact submillimeter emission, complementing these observations with existing Herschel/SPIRE maps. Of our candidate protostars, 24 are associated with young stellar objects (YSOs) in the Spitzer and Herschel/PACS catalogs of 166 and 60 YSOs, respectively (177 unique), confirming their protostellar nature. The remaining 35 candidate protostars are in regions, particularly around LkHalpha 101, where the background cloud emission is too bright to verify or rule out the presence of the compact 70 micron emission that is expected for a protostellar source. We keep these candidate protostars in our sample but note that they may indeed be prestellar in nature. Our observations are sensitive to the high end of the mass distribution in Auriga-Cal. We find that the disparity between the richness of infrared star forming objects in Orion A and the sparsity in Auriga-Cal extends to the submillimeter, suggesting that the relative star formation rates have not varied over the Class II lifetime and that Auriga-Cal will maintain a lower star formation efficiency.

Transition disks with large dust cavities around young stars are promising targets for studying planet formation. Previous studies have revealed the presence of gas cavities inside the dust cavities hinting at recently formed, giant planets. However, many of these studies are biased towards the brightest disks in the nearby star forming regions, and it is not possible to derive reliable statistics that can be compared with exoplanet populations. We present the analysis of 11 transition disks with large cavities (>20 AU radius) from a complete disk survey of the Lupus star forming region, using ALMA Band 7 observations at 0.3" (22-30 AU radius) resolution of the 345 GHz continuum, 13CO and C18O 3-2 observations and the Spectral Energy Distribution of each source. Gas and dust surface density profiles are derived using the physical-chemical modeling code DALI. This is the first study of transition disks of large cavities within a complete disk survey within a star forming region. The dust cavity sizes range from 20-90 AU radius and in three cases, a gas cavity is resolved as well. The deep drops in gas density and large dust cavity sizes are consistent with clearing by giant planets. The fraction of transition disks with large cavities in Lupus is ~11%, which is inconsistent with exoplanet population studies of giant planets at wide orbits. Furthermore, we present a hypothesis of an evolutionary path for large massive disks evolving into transition disks with large cavities.

FU Orionis-type objects are young, low-mass stars with large outbursts in visible light that last for several years or decades. They are thought to represent an evolutionary phase during the life of every young star when accretion from the circumstellar disk is enhanced during recurring time periods. These outbursts are able to rapidly build up the star while affecting the circumstellar disk and thus the ongoing or future planet formation. In many models infall from a circumstellar envelope seems to be necessary to trigger the outbursts. We observed the J=1$-$0 rotational transition of $^{13}$CO and C$^{18}$O towards eight northern FU Orionis-type stars (V1057 Cyg, V1515 Cyg, V2492 Cyg, V2493 Cyg, V1735 Cyg, V733 Cep, RNO 1B and RNO 1C) and derive temperatures and envelope masses and discuss the morphology and kinematics of the circumstellar material. We detected extended CO emission associated with all our targets. Smaller scale CO clumps were found to be associated with five objects with radii of 2000$-$5000 AU and masses of 0.02$-$0.5 $M_{\odot}$; these are clearly heated by the central stars. Three of these envelopes are also strongly detected in the 2.7 mm continuum. No central CO clumps were detected around V733 Cep and V710 Cas but there are many other clumps in their environments. Traces of outflow activity were observed towards V1735 Cyg, V733 Cep and V710 Cas. The diversity of the observed envelopes enables us to set up an evolutionary sequence between the objects. We find their evolutionary state to range from early, embedded Class I stage to late, Class II-type objects with very low-mass circumstellar material. The results reinforce the idea of FU Orionis-type stars as representatives of a transitory stage between embedded Class I young stellar objects and classical T-Tauri stars.

A long-standing open issue of the paradigm of low-mass star formation is the luminosity problem: most protostars are less luminous than theoretically predicted. One possible solution is that the accretion process is episodic. FU Ori-type stars (FUors) are thought to be the visible examples for objects in the high accretion state. FUors are often surrounded by massive envelopes, which replenish the disk material and enable the disk to produce accretion outbursts. However, we have insufficient information on the envelope dynamics in FUors, about where and how mass transfer from the envelope to the disk happens. Here we present ALMA observations of the FUor-type star V346 Nor at 1.3 mm continuum and in different CO rotational lines. We mapped the density and velocity structure of its envelope and analyze the results using channel maps, position-velocity diagrams, and spectro-astrometric methods. We found that V346 Nor is surrounded by gaseous material on 10000 au scale in which a prominent outflow cavity is carved. Within the central $\sim$700 au, the circumstellar matter forms a flattened pseudo-disk where material is infalling with conserved angular momentum. Within $\sim$350 au, the velocity profile is more consistent with a disk in Keplerian rotation around a central star of 0.1 $M_{\odot}$. We determined an infall rate from the envelope onto the disk of 6$\times$10$^{-6}\,M_{\odot}$yr$^{-1}$, a factor of few higher than the quiescent accretion rate from the disk onto the star, hinting for a mismatch between the infall and accretion rates as the cause of the eruption.

The formation pathways of deuterated species trace different regions of protoplanetary disks and may shed light into their physical structure. We aim to constrain the radial extent of main deuterated species; we are particularly interested in spatially characterizing the high and low temperature pathways for enhancing deuteration of these species. We observed the disk surrounding the Herbig Ae star HD 163296 using ALMA in Band 6 and obtained resolved spectral imaging data of DCO$^+$ ($J$=3-2), DCN ($J$=3-2) and N$_2$D$^+$ ($J$=3-2). We model the radial emission profiles of DCO$^+$, DCN and N$_2$D$^+$, assuming their emission is optically thin, using a parametric model of their abundances and radial excitation temperature estimates. DCO$^+$ can be described by a three-region model, with constant-abundance rings centered at 70 AU, 150 AU and 260 AU. The DCN radial profile peaks at about ~60 AU and N$_2$D$^+$ is seen in a ring at ~160 AU. Simple models of both molecules using constant abundances reproduce the data. Assuming reasonable average excitation temperatures for the whole disk, their disk-averaged column densities (and deuterium fractionation ratios) are 1.6-2.6$\times 10^{12}$ cm$^{-2}$ (0.04-0.07), 2.9-5.2$\times 10^{12}$ cm$^{-2}$ ($\sim$0.02) and 1.6-2.5 $\times 10^{11}$ cm$^{-2}$ (0.34-0.45) for DCO$^+$, DCN and N$_2$D$^+$, respectively. Our simple best-fit models show a correlation between the radial location of the first two rings in DCO$^+$ and the DCN and N$_2$D$^+$ abundance distributions that can be interpreted as the high and low temperature deuteration pathways regimes. The origin of the third DCO$^+$ ring at 260 AU is unknown but may be due to a local decrease of ultraviolet opacity allowing the photodesorption of CO or due to thermal desorption of CO as a consequence of radial drift and settlement of dust grains.

The Atacama Large millimetre/submillimetre Array (ALMA) makes use of water vapour radiometers (WVR), which monitor the atmospheric water vapour line at 183 GHz along the line of sight above each antenna to correct for phase delays introduced by the wet component of the troposphere. The application of WVR derived phase corrections improve the image quality and facilitate successful observations in weather conditions that were classically marginal or poor. We present work to indicate that a scaling factor applied to the WVR solutions can act to further improve the phase stability and image quality of ALMA data. We find reduced phase noise statistics for 62 out of 75 datasets from the long-baseline science verification campaign after a WVR scaling factor is applied. The improvement of phase noise translates to an expected coherence improvement in 39 datasets. When imaging the bandpass source, we find 33 of the 39 datasets show an improvement in the signal-to-noise ratio (S/N) between a few to ~30 percent. There are 23 datasets where the S/N of the science image is improved: 6 by <1%, 11 between 1 and 5%, and 6 above 5%. The higher frequencies studied (band 6 and band 7) are those most improved, specifically datasets with low precipitable water vapour (PWV), <1mm, where the dominance of the wet component is reduced. Although these improvements are not profound, phase stability improvements via the WVR scaling factor come into play for the higher frequency (>450 GHz) and long-baseline (>5km) observations. These inherently have poorer phase stability and are taken in low PWV (<1mm) conditions for which we find the scaling to be most effective. A promising explanation for the scaling factor is the mixing of dry and wet air components, although other origins are discussed. We have produced a python code to allow ALMA users to undertake WVR scaling tests and make improvements to their data.

The formation of planets strongly depends on the total amount as well as on the spatial distribution of solids in protoplanetary disks. Thanks to the improvements in resolution and sensitivity provided by ALMA, measurements of the surface density of mm-sized grains are now possible on large samples of disks. Such measurements provide statistical constraints that can be used to inform our understanding of the initial conditions of planet formation. We analyze spatially resolved observations of 36 protoplanetary disks in the Lupus star forming complex from our ALMA survey at 890 micron, aiming to determine physical properties such as the dust surface density, the disk mass and size and to provide a constraint on the temperature profile. We fit the observations directly in the uv-plane using a two-layer disk model that computes the 890 micron emission by solving the energy balance at each disk radius. For 22 out of 36 protoplanetary disks we derive robust estimates of their physical properties. The sample covers stellar masses between ~0.1 and ~2 Solar masses, and we find no trend between the average disk temperatures and the stellar parameters. We find, instead, a correlation between the integrated sub-mm flux (a proxy for the disk mass) and the exponential cut-off radii (a proxy of the disk size) of the Lupus disks. Comparing these results with observations at similar angular resolution of Taurus-Auriga/Ophiuchus disks found in literature and scaling them to the same distance, we observe that the Lupus disks are generally fainter and larger at a high level of statistical significance. Considering the 1-2 Myr age difference between these regions, it is possible to tentatively explain the offset in the disk mass/disk size relation with viscous spreading, however with the current measurements other mechanisms cannot be ruled out.

Three formaldehyde lines were observed (H$_2$CO 3$_{03}$--2$_{02}$, H$_2$CO 3$_{22}$--2$_{21}$, and H$_2$CO 3$_{21}$--2$_{20}$) in the protoplanetary disk around the Herbig Ae star HD 163296 with ALMA at 0.5 arcsecond (60 AU) spatial resolution. H$_2$CO 3$_{03}$--2$_{02}$ was readily detected via imaging, while the weaker H$_2$CO 3$_{22}$--2$_{21}$ and H$_2$CO 3$_{21}$--2$_{20}$ lines required matched filter analysis to detect. H$_2$CO is present throughout most of the gaseous disk, extending out to 550 AU. An apparent 50 AU inner radius of the H$_2$CO emission is likely caused by an optically thick dust continuum. The H$_2$CO radial intensity profile shows a peak at 100 AU and a secondary bump at around 300 AU, suggesting increased production in the outer disk. Different parameterizations of the H$_2$CO abundance were compared to the observed visibilities with $\chi^2$ minimization, using either a characteristic temperature, a characteristic radius or a radial power law index to describe the H$_2$CO chemistry. Similar models were applied to ALMA Science Verification data of C$^{18}$O. In all modeling scenarios, fits to the H$_2$CO data show an increased abundance in the outer disk. The overall best-fit H$_2$CO model shows a factor of two enhancement beyond a radius of 270$\pm$20 AU, with an inner abundance of $2\!-\!5 \times 10^{-12}$. The H$_2$CO emitting region has a lower limit on the kinetic temperature of $T > 20$ K. The C$^{18}$O modeling suggests an order of magnitude depletion in the outer disk and an abundance of $4\!-\!12 \times 10^{-8}$ in the inner disk. The increase in H$_2$CO outer disk emission could be a result of hydrogenation of CO ices on dust grains that are then sublimated via thermal desorption or UV photodesorption, or more efficient gas-phase production beyond about 300 AU if CO is photodisocciated in this region.

We present the first high-resolution sub-mm survey of both dust and gas for a large population of protoplanetary disks. Characterizing fundamental properties of protoplanetary disks on a statistical level is critical to understanding how disks evolve into the diverse exoplanet population. We use ALMA to survey 89 protoplanetary disks around stars with $M_{\ast}>0.1~M_{\odot}$ in the young (1--3~Myr), nearby (150--200~pc) Lupus complex. Our observations cover the 890~$\mu$m continuum and the $^{13}$CO and C$^{18}$O 3--2 lines. We use the sub-mm continuum to constrain $M_{\rm dust}$ to a few Martian masses (0.2--0.4~$M_{\oplus}$) and the CO isotopologue lines to constrain $M_{\rm gas}$ to roughly a Jupiter mass (assuming ISM-like $\rm {[CO]/[H_2]}$ abundance). Of 89 sources, we detect 62 in continuum, 36 in $^{13}$CO, and 11 in C$^{18}$O at $>3\sigma$ significance. Stacking individually undetected sources limits their average dust mass to $\lesssim6$ Lunar masses (0.03~$M_{\oplus}$), indicating rapid evolution once disk clearing begins. We find a positive correlation between $M_{\rm dust}$ and $M_{\ast}$, and present the first evidence for a positive correlation between $M_{\rm gas}$ and $M_{\ast}$, which may explain the dependence of giant planet frequency on host star mass. The mean dust mass in Lupus is 3$\times$ higher than in Upper Sco, while the dust mass distributions in Lupus and Taurus are statistically indistinguishable. Most detected disks have $M_{\rm gas}\lesssim1~M_{\rm Jup}$ and gas-to-dust ratios $<100$, assuming ISM-like $\rm {[CO]/[H_2]}$ abundance; unless CO is very depleted, the inferred gas depletion indicates that planet formation is well underway by a few Myr and may explain the unexpected prevalence of super-Earths in the exoplanet population.

Comets formed in the outer and cold parts of the disk which eventually evolved into our Solar System. Assuming that the comets have undergone no major processing, studying their composition provides insight in the pristine composition of the Solar Nebula. We derive production rates for a number of volatile coma species and explore how molecular line ratios can help constrain the uncertainties of these rates. We analyse observations obtained with the Atacama Large Millimetre/Submillimetre Array of the volatile composition of the comae of comets C/2012 F6 (Lemmon) and C/2012 S1 (ISON) at heliocentric distances of ~1.45 AU and ~0.56 AU, respectively. Assuming a Haser profile with constant outflow velocity, we model the line intensity of each transition using a 3D radiative transfer code and derive molecular production rates and parent scale lengths. We report the first detection of CS in comet ISON obtained with the ALMA array and derive a parent scale length for CS of ~200 km. Due to the high spatial resolution of ALMA, resulting in a synthesised beam with a size slightly smaller than the derived parent scale length, we are able to tentatively identify CS as a daughter species, i.e., a species produced in the coma and/or sublimated from icy grains, rather than a parent species. In addition we report the detection of several CH3OH transitions and confirm the previously reported detections of HCN, HNC and H2CO as well as dust in the coma of each comet, and report 3sigma upper limits for HCO+. We derive molecular production rates relative to water of 0.2% for CS, 0.06-0.1% for HCN, 0.003-0.05% for HNC, 0.1-0.2% for H2CO and 0.5-1.0% for CH3OH, and show that the modelling uncertainties due to unknown collision rates and kinematic temperatures are modest and can be mitigated by available observations of different transitions of HCN.

We performed very deep searches for 2 ground-state water transitions in 13 protoplanetary disks with the HIFI instrument on-board the Herschel Space Observatory, with integration times up to 12 hours per line. Two other water transitions that sample warmer gas were also searched for with shallower integrations. The detection rate is low, and the upper limits provided by the observations are generally much lower than predictions of thermo-chemical models with canonical inputs. One ground-state transition is newly detected in the stacked spectrum of AA Tau, DM Tau, LkCa 15, and MWC 480. We run a grid of models to show that the abundance of gas-phase oxygen needs to be reduced by a factor of at least ~100 to be consistent with the observational upper limits (and positive detections) if a dust-to-gas mass ratio of 0.01 were to be assumed. As a continuation of previous ideas, we propose that the underlying reason for the depletion of oxygen (hence the low detection rate) is the freeze-out of volatiles such as water and CO onto dust grains followed by grain growth and settling/migration, which permanently removes these gas-phase molecules from the emissive upper layers of the outer disk. Such depletion of volatiles is likely ubiquitous among different disks, though not necessarily to the same degree. The volatiles might be returned back to the gas phase in the inner disk (within about 15 AU), which is consistent with current constraints. Comparison with studies on disk dispersal due to photoevaporation indicates that the timescale for volatile depletion is shorter than that of photoevaporation.

H2CO is one of the most readily detected organic molecules in protoplanetary disks. Yet its distribution and dominant formation pathway(s) remain largely unconstrained. To address these issues, we present ALMA observations of two H2CO lines (3_{12}-2_{11} and 5_{15}-4_{14}) at 0".5 (~30 au) spatial resolution toward the disk around the nearby T Tauri star TW Hya. Emission from both lines is spatially resolved, showing a central depression, a peak at 0".4 radius, and a radial decline at larger radii with a bump at ~1", near the millimeter continuum edge. We adopt a physical model for the disk and use toy models to explore the radial and vertical H2CO abundance structure. We find that the observed emission implies the presence of at least two distinct H2CO gas reservoirs: (1) a warm and unresolved inner component (<10 au), and (2) an outer component that extends from ~15 au to beyond the millimeter continuum edge. The outer component is further constrained by the line ratio to arise in a more elevated disk layer at larger radii. The inferred H2CO abundance structure agrees well with disk chemistry models, which predict efficient H2CO gas-phase formation close to the star, and cold H2CO grain surface formation, through H additions to condensed CO, followed by non-thermal desorption in the outer disk. The implied presence of active grain surface chemistry in the TW Hya disk is consistent with the recent detection of CH3OH emission, and suggests that more complex organic molecules are formed in disks, as well.

The protoplanetary system HD 169142 is one of the few cases where a potential candidate protoplanet has been recently detected via direct imaging. To study the interaction between the protoplanet and the disk itself observations of the gas and dust surface density structure are needed. This paper reports new ALMA observations of the dust continuum at 1.3\,mm, $^{12}$CO, $^{13}$CO and C$^{18}$O $J=2-1$ emission from the system HD 169142 at angular resolution of $\sim 0".18 - 0".28$ ($\sim 20\,$au$ - 33\,$au). The dust continuum emission reveals a double-ring structure with an inner ring between $0".17-0".28$ ($\sim 20 - 35\,$au) and an outer ring between $0".48-0".64$ ($\sim 56 - 83\,$au). The size and position of the inner ring is in good agreement with previous polarimetric observations in the near-infrared and is consistent with dust trapping by a massive planet. No dust emission is detected inside the inner dust cavity ($R \lesssim 20\,$au) or within the dust gap ($\sim 35 - 56\,$au). In contrast, the channel maps of the $J=2-1$ line of the three CO isotopologues reveal the presence of gas inside the dust cavity and dust gap. The gaseous disk is also much larger than the compact dust emission extending to $\sim 1'.5$ ($\sim 180\,$au) in radius. This difference and the sharp drop of the continuum emission at large radii point to radial drift of large dust grains ($>$ \micron-size). Using the thermo-chemical disk code \textsc{dali}, the continuum and the CO isotopologues emission are modelled to quantitatively measure the gas and dust surface densities. The resulting gas surface density is reduced by a factor of $\sim 30-40$ inward of the dust gap. The gas and dust distribution hint at the presence of multiple planets shaping the disk structure via dynamical clearing (dust cavity and gap) and dust trapping (double ring dust distribution).

We present 450 and 850 micron submillimetre continuum observations of the IC5146 star-forming region taken as part of the JCMT Gould Belt Survey. We investigate the location of bright submillimetre (clumped) emission with the larger-scale molecular cloud through comparison with extinction maps, and find that these denser structures correlate with higher cloud column density. Ninety-six individual submillimetre clumps are identified using FellWalker and their physical properties are examined. These clumps are found to be relatively massive, ranging from 0.5to 116 MSun with a mean mass of 8 MSun and a median mass of 3.7 MSun. A stability analysis for the clumps suggest that the majority are (thermally) Jeans stable, with M/M_J < 1. We further compare the locations of known protostars with the observed submillimetre emission, finding that younger protostars, i.e., Class 0 and I sources, are strongly correlated with submillimetre peaks and that the clumps with protostars are among the most Jeans unstable. Finally, we contrast the evolutionary conditions in the two major star-forming regions within IC5146: the young cluster associated with the Cocoon Nebula and the more distributed star formation associated with the Northern Streamer filaments. The Cocoon Nebula appears to have converted a higher fraction of its mass into dense clumps and protostars, the clumps are more likely to be Jeans unstable, and a larger fraction of these remaining clumps contain embedded protostars. The Northern Streamer, however, has a larger number of clumps in total and a larger fraction of the known protostars are still embedded within these clumps.

An era has started in which gas and dust can be observed independently in protoplanetary disks, thanks to the recent surveys with ALMA. The first near-complete high-resolution disk survey in both dust and gas in a single star-forming region has been carried out in Lupus, finding surprisingly low gas/dust ratios. The goal of this work is to fully exploit CO isotopologues observations in Lupus, comparing them with physical-chemical model results, in order to obtain gas masses for a large number of disks. We have employed physical-chemical models to analyze continuum and CO isotopologues observations of Lupus disks, including isotope-selective processes and freeze-out. Employing also the ALMA 13CO-only detections, disk gas masses have been calculated for a total of 34 sources, expanding the sample of 10 disks studied by Ansdell et al. (2016), where also C18O was detected. We confirm that overall gas-masses are very low, often smaller than 1 $M_{\rm J}$, if volatile carbon is not depleted. Accordingly, global gas/dust ratios predominantly between 1 and 10. Low CO-based gas masses and gas/dust ratios may indicate rapid loss of gas, or alternatively chemical evolution, e.g. via sequestering of carbon from CO to more complex molecules, or carbon locked up in larger bodies. Current ALMA observations cannot distinguish between these two hypotheses. We have simulated both scenarios, but chemical model results do not allow us to rule out one of the two. Assuming that all Lupus disks have evolved mainly due to viscous processes over the past few Myr, the observed correlation between the current mass accretion rate and dust mass found by Manara et al. (2016) implies a constant gas-to-dust ratio, which is close to 100 based on the observed $M_{\rm disk}/\dot{M}_{\rm acc}$ ratio. This in turn points to a scenario in which carbon depletion is responsible for the low CO isotopologue line luminosities.

Recent high angular resolution ($\sim$0.2") ALMA observations of the 1.1 mm continuum and of HCO+ J=3-2 and HCN J=3-2 gas towards the binary protostar IRS 43 reveal multiple Keplerian disks which are significantly misaligned ($\gt$ 60$^\circ$), both in inclination and position angle and also with respect to the binary orbital plane. Each stellar component has an associated circumstellar disk while the binary is surrounded by a circumbinary disk. Together with archival VLA measurements of the stellar positions over 25 years, and assuming a circular orbit, we use our continuum measurements to determine the binary separation, a = 74 $\pm$ 4 AU, and its inclination, i $\lt$ 30$^\circ$. The misalignment in this system suggests that turbulence has likely played a major role in the formation of IRS 43.

We present observations of the Cepheus Flare obtained as part of the James Clerk Maxwell Telescope (JCMT) Gould Belt Legacy Survey (GBLS) with the SCUBA-2 instrument. We produce a catalogue of sources found by SCUBA-2, and separate these into starless cores and protostars. We determine masses and densities for each of our sources, using source temperatures determined by the Herschel Gould Belt Survey. We compare the properties of starless cores in four different molecular clouds: L1147/58, L1172/74, L1251 and L1228. We find that the core mass functions for each region typically show shallower-than-Salpeter behaviour. We find that L1147/58 and L1228 have a high ratio of starless cores to Class II protostars, while L1251 and L1174 have a low ratio, consistent with the latter regions being more active sites of current star formation, while the former are forming stars less actively. We determine that, if modelled as thermally-supported Bonnor-Ebert spheres, most of our cores have stable configurations accessible to them. We estimate the external pressures on our cores using archival $^{13}$CO velocity dispersion measurements and find that our cores are typically pressure-confined, rather than gravitationally bound. We perform a virial analysis on our cores, and find that they typically cannot be supported against collapse by internal thermal energy alone, due primarily to the measured external pressures. This suggests that the dominant mode of internal support in starless cores in the Cepheus Flare is either non-thermal motions or internal magnetic fields.

We present an overview of the scientific potential of MATISSE, the Multi Aperture mid-Infrared SpectroScopic Experiment for the Very Large Telescope Interferometer. For this purpose we outline selected case studies from various areas, such as star and planet formation, active galactic nuclei, evolved stars, extrasolar planets, and solar system minor bodies and discuss strategies for the planning and analysis of future MATISSE observations. Moreover, the importance of MATISSE observations in the context of complementary high-angular resolution observations at near-infrared and submillimeter/millimeter wavelengths is highlighted.

The Orion A molecular cloud is one of the most well-studied nearby star-forming regions, and includes regions of both highly clustered and more dispersed star formation across its full extent. Here, we analyze dense, star-forming cores identified in the 850 {\mu}m and 450 {\mu}m SCUBA-2 maps from the JCMT Gould Belt Legacy Survey. We identify dense cores in a uniform manner across the Orion A cloud and analyze their clustering properties. Using two independent lines of analysis, we find evidence that clusters of dense cores tend to be mass segregated, suggesting that stellar clusters may have some amount of primordial mass segregation already imprinted in them at an early stage. We also demonstrate that the dense core clusters have a tendency to be elongated, perhaps indicating a formation mechanism linked to the filamentary structure within molecular clouds.

MATISSE is the second-generation mid-infrared spectrograph and imager for the Very Large Telescope Interferometer (VLTI) at Paranal. This new interferometric instrument will allow significant advances by opening new avenues in various fundamental research fields: studying the planet-forming region of disks around young stellar objects, understanding the surface structures and mass loss phenomena affecting evolved stars, and probing the environments of black holes in active galactic nuclei. As a first breakthrough, MATISSE will enlarge the spectral domain of current optical interferometers by offering the L and M bands in addition to the N band. This will open a wide wavelength domain, ranging from 2.8 to 13 um, exploring angular scales as small as 3 mas (L band) / 10 mas (N band). As a second breakthrough, MATISSE will allow mid-infrared imaging - closure-phase aperture-synthesis imaging - with up to four Unit Telescopes (UT) or Auxiliary Telescopes (AT) of the VLTI. Moreover, MATISSE will offer a spectral resolution range from R ~ 30 to R ~ 5000. Here, we present one of the main science objectives, the study of protoplanetary disks, that has driven the instrument design and motivated several VLTI upgrades (GRA4MAT and NAOMI). We introduce the physical concept of MATISSE including a description of the signal on the detectors and an evaluation of the expected performances. We also discuss the current status of the MATISSE instrument, which is entering its testing phase, and the foreseen schedule for the next two years that will lead to the first light at Paranal.

HBC 722 (V2493 Cyg) is a young eruptive star in outburst since 2010. It is an FU Orionis-type object with an atypically low outburst luminosity. Because it was well characterized in the pre-outburst phase, HBC 722 is one of the few FUors where we can learn about the physical changes and processes associated with the eruption. We monitored the source in the BVRIJHKs bands from the ground, and at 3.6 and 4.5 $\mu$m from space with the Spitzer Space Telescope. We analyzed the light curves and the spectral energy distribution by fitting a series of steady accretion disk models at many epochs. We also analyzed the spectral properties of the source based on new optical and infrared spectra. We also mapped HBC 722 and its surroundings at millimeter wavelengths. From the light curve analysis we concluded that the first peak of the outburst in 2010 September was due to an abrupt increase of the accretion rate in the innermost part of the system. This was followed by a long term process, when the brightening was mainly due to a gradual increase of the accretion rate and the emitting area. Our new observations show that the source is currently in a constant plateau phase. We found that around the peak the continuum was bluer and the H$\alpha$ profile changed significantly between 2012 and 2013. The source was not detected in the millimeter continuum, but we discovered a flattened molecular gas structure with a diameter of 1700 au and mass of 0.3 M$_{\odot}$ centered on HBC 722. While the first brightness peak could be interpreted as a rapid fall of piled-up material from the inner disk onto the star, the later monotonic flux rise suggests the outward expansion of a hot component according to the theory of Bell & Lin (1994). Our study of HBC 722 demonstrated that accretion-related outbursts can occur in young stellar objects even with very low mass disks, in the late Class II phase.

We present the JCMT Gould Belt Survey's first look results of the southern extent of the Orion A Molecular Cloud ($\delta \leq -5\mathrm{:}31\mathrm{:}27.5$). Employing a two-step structure identification process, we construct individual catalogues for large-scale regions of significant emission labelled as islands and smaller-scale subregions called fragments using the 850 $\mu$m continuum maps obtained using SCUBA-2. We calculate object masses, sizes, column densities, and concentrations. We discuss fragmentation in terms of a Jeans instability analysis and highlight interesting structures as candidates for follow-up studies. Furthermore, we associate the detected emission with young stellar objects (YSOs) identified by Spitzer and Herschel. We find that although the population of active star-forming regions contains a wide variety of sizes and morphologies, there is a strong positive correlation between the concentration of an emission region and its calculated Jeans instability. There are, however, a number of highly unstable subregions in dense areas of the map that show no evidence of star formation. We find that only $\sim$72\% of the YSOs defined as Class 0+I and flat-spectrum protostars coincide with dense 850 $\mu$m emission structures (column densities $>3.7\times10^{21}\mathrm{\:cm}^{-2}$). The remaining 28\% of these objects, which are expected to be embedded in dust and gas, may be misclassified. Finally, we suggest that there is an evolution in the velocity dispersion of young stellar objects such that sources which are more evolved are associated with higher velocities.

The dust emissivity spectral index, $\beta$, is a critical parameter for deriving the mass and temperature of star-forming structures, and consequently their gravitational stability. The $\beta$ value is dependent on various dust grain properties, such as size, porosity, and surface composition, and is expected to vary as dust grains evolve. Here we present $\beta$, dust temperature, and optical depth maps of the star-forming clumps in the Perseus Molecular Cloud determined from fitting SEDs to combined Herschel and JCMT observations in the 160 $\mu$m, 250 $\mu$m, 350 $\mu$m, 500 $\mu$m, and 850 $\mu$m bands. Most of the derived $\beta$, and dust temperature values fall within the ranges of 1.0 - 2.7 and 8 - 20 K, respectively. In Perseus, we find the $\beta$ distribution differs significantly from clump to clump, indicative of grain growth. Furthermore, we also see significant, localized $\beta$ variations within individual clumps and find low $\beta$ regions correlate with local temperature peaks, hinting at the possible origins of low $\beta$ grains. Throughout Perseus, we also see indications of heating from B stars and embedded protostars, as well evidence of outflows shaping the local landscape.

The composition of planetary solids and gases is largely rooted in the processing of volatile elements in protoplanetary disks. To shed light on the key processes, we carry out a comparative analysis of the gas-phase carbon abundance in two systems with a similar age and disk mass, but different central stars: HD 100546 and TW Hya. We combine our recent detections of C$^{0}$ in these disks with observations of other carbon reservoirs (CO, C$^{+}$, C$_{2}$H) and gas mass and warm gas tracers (HD, O$^{0}$), as well as spatially resolved ALMA observations and the spectral energy distribution. The disks are modelled with the DALI 2D physical-chemical code. Stellar abundances for HD 100546 are derived from archival spectra. Upper limits on HD emission from HD 100546 place an upper limit on the total disk mass of $\leq0.1\,M_{\odot}$. The gas-phase carbon abundance in the atmosphere of this warm Herbig disk is at most moderately depleted compared to the interstellar medium, with [C]/[H]$_{\rm gas}=(0.1-1.5)\times 10^{-4}$. HD 100546 itself is a $\lambda\,$Bo\"{o}tis star, with solar abundances of C and O but a strong depletion of rock-forming elements. In the gas of the T Tauri disk TW Hya, both C and O are strongly underabundant, with [C]/[H]$_{\rm gas}=(0.2-5.0)\times 10^{-6}$ and C/O${>}1$. We discuss evidence that the gas-phase C and O abundances are high in the warm inner regions of both disks. Our analytical model, including vertical mixing and a grain size distribution, reproduces the observed [C]/[H]$_{\rm gas}$ in the outer disk of TW Hya and allows to make predictions for other systems.

We present SCUBA-2 450{\mu}m and 850{\mu}m observations of the W40 complex in the Serpens-Aquila region as part of the James Clerk Maxwell Telescope (JCMT) Gould Belt Survey (GBS) of nearby star-forming regions. We investigate radiative heating by constructing temperature maps from the ratio of SCUBA-2 fluxes using a fixed dust opacity spectral index, {\beta} = 1.8, and a beam convolution kernel to achieve a common 14.8" resolution. We identify 82 clumps ranging between 10 and 36K with a mean temperature of 20{\pm}3K. Clump temperature is strongly correlated with proximity to the external OB association and there is no evidence that the embedded protostars significantly heat the dust. We identify 31 clumps that have cores with densities greater than 105cm{^{-3}}. Thirteen of these cores contain embedded Class 0/I protostars. Many cores are associated with bright-rimmed clouds seen in Herschel 70 {\mu}m images. From JCMT HARP observations of the 12CO 3-2 line, we find contamination of the 850{\mu}m band of up to 20 per cent. We investigate the free-free contribution to SCUBA-2 bands from large-scale and ultracompact H ii regions using archival VLA data and find the contribution is limited to individual stars, accounting for 9 per cent of flux per beam at 450 {\mu}m or 12 per cent at 850 {\mu}m in these cases. We conclude that radiative heating has potentially influenced the formation of stars in the Dust Arc sub-region, favouring Jeans stable clouds in the warm east and fragmentation in the cool west.

The young star beta Pictoris is well known for its dusty debris disk, produced through the grinding down by collisions of planetesimals, kilometre-sized bodies in orbit around the star. In addition to dust, small amounts of gas are also known to orbit the star, likely the result from vaporisation of violently colliding dust grains. The disk is seen edge on and from previous absorption spectroscopy we know that the gas is very rich in carbon relative to other elements. The oxygen content has been more difficult to assess, however, with early estimates finding very little oxygen in the gas at a C/O ratio 20x higher than the cosmic value. A C/O ratio that high is difficult to explain and would have far-reaching consequences for planet formation. Here we report on observations by the far-infrared space telescope Herschel, using PACS, of emission lines from ionised carbon and neutral oxygen. The detected emission from C+ is consistent with that previously reported being observed by the HIFI instrument on Herschel, while the emission from O is hard to explain without assuming a higher-density region in the disk, perhaps in the shape of a clump or a dense torus, required to sufficiently excite the O atoms. A possible scenario is that the C/O gas is produced by the same process responsible for the CO clump recently observed by ALMA in the disk, and that the re-distribution of the gas takes longer than previously assumed. A more detailed estimate of the C/O ratio and the mass of O will have to await better constraints on the C/O gas spatial distribution.

The gas temperature structure of protoplanetary disks is a key ingredient for interpreting various disk observations and for quantifying the subsequent evolution of these systems. The comparison of low- and mid-$J$ CO rotational lines is a powerful tool to assess the temperature gradient in the warm molecular layer of disks. Spectrally resolved high-$J$ ($J_{\rm u} > 14$) CO lines probe intermediate distances and heights from the star that are not sampled by (sub-)millimeter CO spectroscopy. This paper presents new {\it Herschel}/HIFI and archival PACS observations of $^{12}$CO, $^{13}$CO and \cii \ emission in 4 Herbig AeBe (HD 100546, HD 97048, IRS 48, HD 163296) and 3 T Tauri (AS 205, S CrA, TW Hya) disks. In the case of the T Tauri systems AS 205 and S CrA, the CO emission has a single-peaked profile, likely due to a slow wind. For all other systems, the {\it Herschel} CO spectra are consistent with pure disk emission and the spectrally-resolved lines (HIFI) and the CO rotational ladder (PACS) are analyzed simultaneously assuming power-law temperature and column density profiles, using the velocity profile to locate the emission in the disk. The temperature profile varies substantially from disk to disk. In particular, $T_{\rm gas}$ in the disk surface layers can differ by up to an order of magnitude among the 4 Herbig AeBe systems with HD 100546 being the hottest and HD 163296 the coldest disk of the sample. Clear evidence of a warm disk layer where $T_{\rm gas} > T_{\rm dust}$ is found in all the Herbig Ae disks. The observed CO fluxes and line profiles are compared to predictions of physical-chemical models. The primary parameters affecting the disk temperature structure are the flaring angle, the gas-to-dust mass ratio the scale height and the dust settling.

Nitrogen chemistry in protoplanetary disks and the freeze-out on dust particles is key to understand the formation of nitrogen bearing species in early solar system analogs. So far, ammonia has not been detected beyond the snowline in protoplanetary disks. We aim to find gas-phase ammonia in a protoplanetary disk and characterize its abundance with respect to water vapor. Using HIFI on the Herschel Space Observatory we detect, for the first time, the ground-state rotational emission of ortho-NH$_3$ in a protoplanetary disk, around TW Hya. We use detailed models of the disk's physical structure and the chemistry of ammonia and water to infer the amounts of gas-phase molecules of these species. We explore two radial distributions ( confined to $<$60 au like the millimeter-sized grains) and two vertical distributions (near the midplane where water is expected to photodesorb off icy grains) to describe the (unknown) location of the molecules. These distributions capture the effects of radial drift and vertical settling of ice-covered grains. We use physical-chemical models to reproduce the fluxes with assuming that water and ammonia are co-spatial. We infer ammonia gas-phase masses of 0.7-11.0 $\times$10$^{21}$ g. For water, we infer gas-phase masses of 0.2-16.0 $\times$10$^{22}$ g. This corresponds to NH$_3$/H$_2$O abundance ratios of 7\%-84\%, assuming that water and ammonia are co-located. Only in the most compact and settled adopted configuration is the inferred NH$_3$/H$_2$O consistent with interstellar ices and solar system bodies of $\sim$ 5\%-10\%. Volatile release in the midplane may occur via collisions between icy bodies if the available surface for subsequent freeze-out is significantly reduced, e.g., through growth of small grains into pebbles or larger.

Hydrides are simple compounds containing one or a few hydrogen atoms bonded to a heavier atom. They are fundamental precursor molecules in cosmic chemistry and many hydride ions have become observable in high quality for the first time thanks to the Herschel Space Observatory. Ionized hydrides, such as CH+ and OH+, and also HCO+ that affect the chemistry of molecules such as water, provide complementary information on irradiation by far UV (FUV) or X-rays and gas temperature. The targeted lines of CH+, OH+, H2O+, C+ and CH are detected mostly in blue-shifted absorption. H3O+ and SH+ are detected in emission and only toward some high-mass objects. The observed line parameters and correlations suggest two different origins, related to gas entrained by the outflows and to the circumstellar envelope. The column density ratios of CH+/OH+ are estimated from chemical slab models, assuming that the H2 density is given by the specific density model of each object at the beam radius. For the low-mass YSOs the observed ratio can be reproduced for an FUV flux of 2-400 times the ISRF at the location of the molecules. In two high-mass objects, the UV flux is 20-200 times the ISRF derived from absorption lines, and 300-600 ISRF using emission lines. If the FUV flux required for low-mass objects originates at the central protostar, a substantial FUV luminosity, up to 1.5 L_sun, is required. There is no molecular evidence for X-ray induced chemistry in the low-mass objects on the observed scales of a few 1000 AU. For high-mass regions, the FUV flux required to produce the observed molecular ratios is smaller than the unattenuated flux expected from the central object(s) at the Herschel beam radius. This is consistent with an FUV flux reduced by circumstellar extinction or by bloating of the protostar.

The JCMT Gould Belt Legacy Survey obtained SCUBA-2 observations of dense cores within three sub-regions of Orion B: LDN 1622, NGC 2023/2024, and NGC 2068/2071, all of which contain clusters of cores. We present an analysis of the clustering properties of these cores, including the two-point correlation function and Cartwright's Q parameter. We identify individual clusters of dense cores across all three regions using a minimal spanning tree technique, and find that in each cluster, the most massive cores tend to be centrally located. We also apply the independent M-Sigma technique and find a strong correlation between core mass and the local surface density of cores. These two lines of evidence jointly suggest that some amount of mass segregation in clusters has happened already at the dense core stage.

An unsolved problem in step-wise core-accretion planet formation is that rapid radial drift in gas-rich protoplanetary disks should drive mm-/meter-sized particles inward to the central star before large bodies can form. One promising solution is to confine solids within small scale structures. Here we investigate dust structures in the (sub)mm continuum emission of four disks (TW Hya, HL Tau, HD 163296 and DM Tau), a sample of disks with the highest spatial resolution ALMA observations to date. We retrieve the surface brightness distributions using synthesized images and fitting visibilities with analytical functions. We find that the continuum emission of the four disks is ~axi-symmetric but rich in 10-30AU-sized radial structures, possibly due to physical gaps, surface density enhancements or localized dust opacity variations within the disks. These results suggest that small scale axi-symmetric dust structures are likely to be common, as a result of ubiquitous processes in disk evolution and planet formation. Compared with recent spatially resolved observations of CO snowlines in these same disks, all four systems show enhanced continuum emission from regions just beyond the CO condensation fronts, potentially suggesting a causal relationship between dust growth/trapping and snowlines.

In recent years there has been a paradigm shift from centralised to geographically distributed resources. Individual entities are no longer able to host or afford the necessary expertise in-house, and, as a consequence, society increasingly relies on widespread collaborations. Although such collaborations are now the norm for scientific projects, more technical structures providing support to a distributed scientific community without direct financial or other material benefits are scarce. The network of European ALMA Regional Centre (ARC) nodes is an example of such an internationally distributed user support network. It is an organised effort to provide the European ALMA user community with uniform expert support to enable optimal usage and scientific output of the ALMA facility. The network model for the European ARC nodes is described in terms of its organisation, communication strategies and user support.

The gas-solid budget of carbon in protoplanetary disks is related to the composition of the cores and atmospheres of the planets forming in them. The key gas-phase carbon carriers CO, C$^{0}$ and C$^{+}$ can now be observed in disks. The gas-phase carbon abundance in disks has not yet been well characterized, we aim to obtain new constraints on the [C]/[H] ratio in a sample of disks, and to get an overview of the strength of [CI] and warm CO emission. We carried out a survey of the CO$\,6$--$5$ and [CI]$\,1$--$0$ and $2$--$1$ lines towards $37$ disks with APEX, and supplemented it with [CII] data from the literature. The data are interpreted using a grid of models produced with the DALI code. We also investigate how well the gas-phase carbon abundance can be determined in light of parameter uncertainties. The CO$\,6$--$5$ line is detected in $13$ out of $33$ sources, the [CI]$\,1$--$0$ in $6$ out of $12$, and the [CI]$\,2$--$1$ in $1$ out of $33$. With deep integrations, the first unambiguous detections of [CI]~$1$--$0$ in disks are obtained, in TW~Hya and HD~100546. Gas-phase carbon abundance reductions of a factor $5$--$10$ or more can be identified robustly based on CO and [CI] detections. The atomic carbon detection in TW~Hya confirms a factor $100$ reduction of [C]/[H]$_{\rm gas}$ in that disk, while the data are consistent with an ISM-like carbon abundance for HD~100546. In addition, BP~Tau, T~Cha, HD~139614, HD~141569, and HD~100453 are either carbon-depleted or gas-poor disks. The low [CI]~$2$--$1$ detection rates in the survey mostly reflect insufficient sensitivity to detect T~Tauri disks. The Herbig~Ae/Be disks with CO and [CII] upper limits below the models are debris disk like systems. A roughly order of magnitude increase in sensitivity compared to our survey is required to obtain useful constraints on the gas-phase [C]/[H] ratio in most of the targeted systems.

Grain growth in planet-forming disks is the first step toward the formation of planets. The growth of grains and their inward drift leaves a distinct imprint on the dust surface-density distribution and the resulting surface-brightness profile of the thermal continuum emission. We determine the surface-brightness profile of the continuum emission using resolved observations at millimeter wavelengths of the disk around TW Hya, and infer the signature of dust evolution on the surface density and dust opacity. Archival ALMA observations at 820 micron on baselines up to 410 kilolambda are compared to parametrized disk models to determine the surface-brightness profile. Under the assumption of a constant dust opacity, a broken radial power law best describes the dust surface density, with a slope of -0.53 +/- 0.01 from the 4.1 au radius of the (already known) inner hole to a turn-over radius of 47.1 +/- 0.2 au, steepening to -8.0 +/- 0.1 at larger radii. The emission drops below the detection limit beyond ~60 au. The shape of the dust surface density is consistent with theoretical expectations for grain growth, fragmentation, and drift, but its total dust content and its turn-over radius are too large for TW Hya's age of 8-10 Myr even when taking into account a radially varying dust opacity. Higher resolution imaging with ALMA of TW Hya and other disks is required to establish if unseen gaps associated with, e.g., embedded planets trap grains at large radii or if locally enhanced grain growth associated with the CO snow line explains the extent of the millimeter-continuum surface brightness profile. In the latter case, population studies should reveal a correlation between the location of the CO snow line and the extent of the millimeter continuum. In the former case, and if CO freeze out promotes planet formation, this correlation should extend to the location of gaps as well.

We present a first look at the SCUBA-2 observations of three sub-regions of the Orion B molecular cloud: LDN 1622, NGC 2023/2024, and NGC 2068/2071, from the JCMT Gould Belt Legacy Survey. We identify 29, 564, and 322 dense cores in L1622, NGC 2023/2024, and NGC 2068/2071 respectively, using the SCUBA-2 850 micron map, and present their basic properties, including their peak fluxes, total fluxes, and sizes, and an estimate of the corresponding 450 micron peak fluxes and total fluxes, using the FellWalker source extraction algorithm. Assuming a constant temperature of 20 K, the starless dense cores have a mass function similar to that found in previous dense core analyses, with a Salpeter-like slope at the high-mass end. The majority of cores appear stable to gravitational collapse when considering only thermal pressure; indeed, most of the cores which have masses above the thermal Jeans mass are already associated with at least one protostar. At higher cloud column densities, above 1-2 x 10^23 cm^-2, most of the mass is found within dense cores, while at lower cloud column densities, below 1 x 10^23 cm^-2, this fraction drops to 10% or lower. Overall, the fraction of dense cores associated with a protostar is quite small (<8%), but becomes larger for the densest and most centrally concentrated cores. NGC 2023 / 2024 and NGC 2068/2071 appear to be on the path to forming a significant number of stars in the future, while L1622 has little additional mass in dense cores to form many new stars.

Performing ground-based submillimetre observations is a difficult task as the measurements are subject to absorption and emission from water vapour in the Earth's atmosphere and time variation in weather and instrument stability. Removing these features and other artifacts from the data is a vital process which affects the characteristics of the recovered astronomical structure we seek to study. In this paper, we explore two data reduction methods for data taken with the Submillimetre Common-User Bolometer Array-2 (SCUBA-2) at the James Clerk Maxwell Telescope (JCMT). The JCMT Legacy Reduction 1 (JCMT LR1) and The Gould Belt Legacy Survey Legacy Release 1 (GBS LR1) reduction both use the same software, Starlink, but differ in their choice of data reduction parameters. We find that the JCMT LR1 reduction is suitable for determining whether or not compact emission is present in a given region and the GBS LR1 reduction is tuned in a robust way to uncover more extended emission, which better serves more in-depth physical analyses of star-forming regions. Using the GBS LR1 method, we find that compact sources are recovered well, even at a peak brightness of only 3 times the noise, whereas the reconstruction of larger objects requires much care when drawing boundaries around the expected astronomical signal in the data reduction process. Incorrect boundaries can lead to false structure identification or it can cause structure to be missed. In the JCMT LR1 reduction, the extent of the true structure of objects larger than a point source is never fully recovered.

We present far-infrared and submillimeter maps from the Herschel Space Observatory and the James Clerk Maxwell Telescope of the debris disk host star AU Microscopii. Disk emission is detected at 70, 160, 250, 350, 450, 500 and 850 micron. The disk is resolved at 70, 160 and 450 micron. In addition to the planetesimal belt, we detect thermal emission from AU Mic's halo for the first time. In contrast to the scattered light images, no asymmetries are evident in the disk. The fractional luminosity of the disk is $3.9 \times 10^{-4}$ and its mm-grain dust mass is 0.01 MEarth (+/- 20%). We create a simple spatial model that reconciles the disk SED as a blackbody of 53 +/- 2 K (a composite of 39 and 50 K components) and the presence of small (non-blackbody) grains which populate the extended halo. The best fit model is consistent with the "birth ring" model explored in earlier works, i.e., an edge-on dust belt extending from 8.8-40 AU, but with an additional halo component with an $r^{-1.5}$ surface density profile extending to the limits of sensitivity (140 AU). We confirm that AU Mic does not exert enough radiation force to blow out grains. For stellar mass loss rates of 10-100x solar, compact (zero porosity) grains can only be removed if they are very small, consistently with previous work, if the porosity is 0.9, then grains approaching 0.1 micron can be removed via corpuscular forces (i.e., the stellar wind).

Intermediate mass protostarsprovide a bridge between theories of low- and high-mass star formation. Emerging molecular outflows can be used to determine the influence of fragmentation and multiplicity on protostellar evolution through the correlation of outflow forces of intermediate mass protostars with the luminosity. The aim of this paper is to derive outflow forces from outflows of six intermediate mass protostellar regions and validate the apparent correlation between total luminosity and outflow force seen in earlier work, as well as remove uncertainties caused by different methodology. By comparing CO 6--5 observations obtained with APEX with non-LTE radiative transfer model predictions, optical depths, temperatures, densities of the gas of the molecular outflows are derived. Outflow forces, dynamical timescales and kinetic luminosities are subsequently calculated. Outflow parameters, including the forces, were derived for all sources. Temperatures in excess of 50 K were found for all flows, in line with recent low-mass results. However, comparison with other studies could not corroborate conclusions from earlier work on intermediate mass protostars which hypothesized that fragmentation enhances outflow forces in clustered intermediate mass star formation. Any enhancement in comparison with the classical relation between outflow force and luminosity can be attributed the use of a higher excitation line and improvement in methods; They are in line with results from low-mass protostars using similar techniques. The role of fragmentation on outflows is an important ingredient to understand clustered star formation and the link between low and high-mass star formation. However, detailed information on spatial scales of a few 100 AU, covering all individual members is needed to make the necessary progress.

An abundance decrease in carbon- and oxygen-bearing species relative to dust has been frequently found in planet-forming disks, which can be attributed to an overall reduction of gas mass. However, in the case of TW Hya, the only disk with gas mass measured directly with HD rotational lines, the inferred gas mass ($\lesssim$0.005 solar mass) is significantly below the directly measured value ($\gtrsim$0.05 solar mass). We show that this apparent conflict can be resolved if the elemental abundances of carbon and oxygen are reduced in the upper layers of the outer disk but are normal elsewhere (except for a possible enhancement of their abundances in the inner disk). The implication is that in the outer disk, the main reservoir of the volatiles (CO, water, ...) resides close to the midplane, locked up inside solid bodies that are too heavy to be transported back to the atmosphere by turbulence. An enhancement in the carbon and oxygen abundances in the inner disk can be caused by inward migration of these solid bodies. This is consistent with estimates based on previous models of dust grain dynamics. Indirect measurements of the disk gas mass and disk structure from species such as CO will thus be intertwined with the evolution of dust grains, and possibly also with the formation of planetesimals.

A major goal of the Atacama Large Millimeter/submillimeter Array (ALMA) is to make accurate images with resolutions of tens of milliarcseconds, which at submillimeter (submm) wavelengths requires baselines up to ~15 km. To develop and test this capability, a Long Baseline Campaign (LBC) was carried out from September to late November 2014, culminating in end-to-end observations, calibrations, and imaging of selected Science Verification (SV) targets. This paper presents an overview of the campaign and its main results, including an investigation of the short-term coherence properties and systematic phase errors over the long baselines at the ALMA site, a summary of the SV targets and observations, and recommendations for science observing strategies at long baselines. Deep ALMA images of the quasar 3C138 at 97 and 241 GHz are also compared to VLA 43 GHz results, demonstrating an agreement at a level of a few percent. As a result of the extensive program of LBC testing, the highly successful SV imaging at long baselines achieved angular resolutions as fine as 19 mas at ~350 GHz. Observing with ALMA on baselines of up to 15 km is now possible, and opens up new parameter space for submm astronomy.

We describe the design, construction, and characterization of the Band 9 heterodyne receivers (600-720 GHz) for the Atacama Large Millimeter / submillimeter Array (ALMA). The ALMA Band 9 receiver units ("cartridges"), which are installed in the telescope's front end, have been designed to detect and down-convert two orthogonal linear polarization components of the light collected by the ALMA antennas. The light entering the front end is refocused with a compact arrangement of mirrors, which is fully contained within the cartridge. The arrangement contains a grid to separate the polarizations and two beam splitters to combine each resulting beam with a local oscillator signal. The combined beams are fed into independent double-sideband mixers, each with a corrugated feedhorn coupling the radiation by way of a waveguide with backshort cavity into an impedance-tuned SIS junction that performs the heterodyne down-conversion. Finally, the generated intermediate frequency signals are amplified by cryogenic and room-temperature HEMT amplifiers and exported to the telescope's back end for further processing and, finally, correlation. The receivers have been constructed and tested in the laboratory and they show excellent performance, complying with ALMA requirements. Performance statistics on all 73 Band 9 receivers are reported. On-sky characterization and tests of the performance of the Band 9 cartridges are presented using commissioning data.

We present 850$\mu$m and 450$\mu$m data from the JCMT Gould Belt Survey obtained with SCUBA-2 and characterise the dust attributes of Class I, Class II and Class III disk sources in L1495. We detect 23% of the sample at both wavelengths, with the detection rate decreasing through the Classes from I--III. The median disk mask is 1.6$\times 10^{-3}$M$_{\odot}$, and only 7% of Class II sources have disk masses larger than 20 Jupiter masses. We detect a higher proportion of disks towards sources with stellar hosts of spectral type K than spectral type M. Class II disks with single stellar hosts of spectral type K have higher masses than those of spectral type M, supporting the hypothesis that higher mass stars have more massive disks. Variations in disk masses calculated at the two wavelengths suggests there may be differences in dust opacity and/or dust temperature between disks with hosts of spectral types K to those with spectral type M.

In this paper we present the first observations of the Ophiuchus molecular cloud performed as part of the James Clerk Maxwell Telescope (JCMT) Gould Belt Survey (GBS) with the SCUBA-2 instrument. We demonstrate methods for combining these data with previous HARP CO, Herschel, and IRAM N$_{2}$H$^{+}$ observations in order to accurately quantify the properties of the SCUBA-2 sources in Ophiuchus. We produce a catalogue of all of the sources found by SCUBA-2. We separate these into protostars and starless cores. We list all of the starless cores and perform a full virial analysis, including external pressure. This is the first time that external pressure has been included in this level of detail. We find that the majority of our cores are either bound or virialised. Gravitational energy and external pressure are on average of a similar order of magnitude, but with some variation from region to region. We find that cores in the Oph A region are gravitationally bound prestellar cores, while cores in the Oph C and E regions are pressure-confined. We determine that N$_{2}$H$^{+}$ is a good tracer of the bound material of prestellar cores, although we find some evidence for N$_{2}$H$^{+}$ freeze-out at the very highest core densities. We find that non-thermal linewidths decrease substantially between the gas traced by C$^{18}$O and that traced by N$_{2}$H$^{+}$, indicating the dissipation of turbulence at higher densities. We find that the critical Bonnor-Ebert stability criterion is not a good indicator of the boundedness of our cores. We detect the pre-brown dwarf candidate Oph B-11 and find a flux density and mass consistent with previous work. We discuss regional variations in the nature of the cores and find further support for our previous hypothesis of a global evolutionary gradient across the cloud from southwest to northeast, indicating sequential star formation across the region.

During the embedded stage of star formation, bipolar molecular outflows and UV radiation from the protostar are important feedback processes. Our aim is to quantify the feedback, mechanical and radiative, for a large sample of low-mass sources. The outflow activity is compared to radiative feedback in the form of UV heating by the accreting protostar to search for correlations and evolutionary trends. Large-scale maps of 26 young stellar objects, which are part of the Herschel WISH key program are obtained using the CHAMP+ instrument on the APEX (12CO and 13CO 6-5), and the HARP-B instrument on the JCMT (12CO and 13CO 3-2). Maps are used to determine outflow parameters and envelope models are used to quantify the amount of UV-heated gas and its temperature from 13CO 6-5 observations. All sources in our sample show outflow activity and the outflow force, F_CO, is larger for Class 0 sources than for Class I sources, even if their luminosities are comparable. The outflowing gas typically extends to much greater distances than the power-law envelope and therefore influences the surrounding cloud material directly. Comparison of the CO 6-5 results with Herschel-HIFI H2O and PACS high-J CO lines, both tracing currently shocked gas, shows that the two components are linked, even though the transitions do not probe the same gas. The link does not extend down to CO 3-2. The conclusion is that CO 6-5 depends on the shock characteristics (density and velocity), whereas CO 3-2 is more sensitive to conditions in the surrounding environment (density). The radiative feedback is responsible for increasing the gas temperature by a factor of two, up to 30-50 K, on scales of a few thousand AU, particularly along the direction of the outflow. The mass of the UV heated gas exceeds the mass contained in the entrained outflow in the inner ~3000 AU and is therefore at least as important on small scales.

Context: Outflows are an important part of the star formation process as both the result of ongoing active accretion and one of the main sources of mechanical feedback on small scales. Water is the ideal tracer of these effects because it is present in high abundance in various parts of the protostar. Method: We present \textit{Herschel} HIFI spectra of multiple water-transitions towards 29 nearby Class 0/I protostars as part of the WISH Survey. These are decomposed into different Gaussian components, with each related to one of three parts of the protostellar system; quiescent envelope, cavity shock and spot shocks in the jet and at the base of the outflow. We then constrain the excitation conditions present in the two outflow-related components. Results: Water emission is optically thick but effectively thin, with line ratios that do not vary with velocity, in contrast to CO. The physical conditions of the cavity and spot shocks are similar, with post-shock H$_{2}$ densities of order 10$^{5}-$10$^{8}$\,cm$^{-3}$ and H$_{2}$O column densities of order 10$^{16}-$10$^{18}$\,cm$^{-2}$. H$_{2}$O emission originates in compact emitting regions: for the spot shocks these correspond to point sources with radii of order 10-200\,AU, while for the cavity shocks these come from a thin layer along the outflow cavity wall with thickness of order 1-30\,AU. Conclusions: Water emission at the source position traces two distinct kinematic components in the outflow; J shocks at the base of the outflow or in the jet, and C shocks in a thin layer in the cavity wall. Class I sources have similar excitation conditions to Class 0 sources, but generally smaller line-widths and emitting region sizes. We suggest that it is the velocity of the wind driving the outflow, rather than the decrease in envelope density or mass, that is the cause of the decrease in H$_{2}$O intensity between Class 0 and I.

We present SCUBA-2 450micron and 850micron observations of the Serpens MWC 297 region, part of the JCMT Gould Belt Survey of nearby star-forming regions. Simulations suggest that radiative feedback influences the star-formation process and we investigate observational evidence for this by constructing temperature maps. Maps are derived from the ratio of SCUBA-2 fluxes and a two component model of the JCMT beam for a fixed dust opacity spectral index of beta = 1.8. Within 40 of the B1.5Ve Herbig star MWC 297, the submillimetre fluxes are contaminated by free-free emission with a spectral index of 1.03+-0.02, consistent with an ultra-compact HII region and polar winds/jets. Contamination accounts for 73+-5 per cent and 82+-4 per cent of peak flux at 450micron and 850micron respectively. The residual thermal disk of the star is almost undetectable at these wavelengths. Young Stellar Objects are confirmed where SCUBA-2 850micron clumps identified by the fellwalker algorithm coincide with Spitzer Gould Belt Survey detections. We identify 23 objects and use Tbol to classify nine YSOs with masses 0.09 to 5.1 Msun. We find two Class 0, one Class 0/I, three Class I and three Class II sources. The mean temperature is 15+-2K for the nine YSOs and 32+-4K for the 14 starless clumps. We observe a starless clump with an abnormally high mean temperature of 46+-2K and conclude that it is radiatively heated by the star MWC 297. Jeans stability provides evidence that radiative heating by the star MWC 297 may be suppressing clump collapse.

We present high angular resolution Submillimeter Array observations ofthe outbursting Jupiter family comet 17P/Holmes on 2007 October 26-29, achieving a spatial resolution of 2.5", or ~3000 km at the comet distance. The observations resulted in detections of the rotational lines CO 3-2, HCN 4-3, H$^{13}$CN 4-3, CS 7-6, H$_2$CO 3$_{1,2}$-2$_{1,1}$, H$_2$S 2$_{2,0}$-2$_{1,1}$, and multiple CH$_3$OH lines, along with the associated dust continuum at 221 and 349 GHz. The continuum has a spectral index of 2.7$\pm$0.3, slightly steeper than blackbody emission from large dust particles. From the imaging data, we identify two components in the molecular emission. One component is characterized by a relatively broad line width (~1 km s$^{-1}$ FWHM) exhibiting a symmetric outgassing pattern with respect to the nucleus position. The second component has a narrower line width (<0.5 km s$^{-1}$ FWHM) with the line center red-shifted by 0.1-0.2 km s$^{-1}$ (cometocentric frame), and shows a velocity shift across the nucleus position with the position angle gradually changing from 66 to 30 degrees within the four days of observations. We determine distinctly different CO/HCN ratios for each of the components. For the broad-line component we find CO/HCN <7, while in the narrow-line component, CO/HCN = 40$\pm$5. We hypothesize that the narrow-line component originates from the ice grain halo found in near-nucleus photometry, believed to be created by sublimating recently released ice grains around the nucleus during the outburst. In this interpretation, the high CO/HCN ratio of this component reflects the more pristine volatile composition of nucleus material released in the outburst.

CO, $^{13}$CO and C$^{18}$O ${\it J}$ = 3--2 observations are presented of the Ophiuchus molecular cloud. The $^{13}$CO and C$^{18}$O emission is dominated by the Oph A clump, and the Oph B1, B2, C, E, F and J regions. The optically thin(ner) C$^{18}$O line is used as a column density tracer, from which the gravitational binding energy is estimated to be $4.5 \times 10^{39}$ J (2282 $M_\odot$ km$^2$ s$^{-2}$). The turbulent kinetic energy is $6.3 \times 10^{38}$ J (320 $M_\odot$ km$^2$ s$^{-2}$), or 7 times less than this, and therefore the Oph cloud as a whole is gravitationally bound. Thirty protostars were searched for high velocity gas, with eight showing outflows, and twenty more having evidence of high velocity gas along their lines-of-sight. The total outflow kinetic energy is $1.3 \times 10^{38}$ J (67 $M_\odot$ km$^2$ s$^{-2}$), corresponding to 21$\%$ of the cloud's turbulent kinetic energy. Although turbulent injection by outflows is significant, but does ${\it not}$ appear to be the dominant source of turbulence in the cloud. 105 dense molecular clumplets were identified, which had radii $\sim$ 0.01--0.05 pc, virial masses $\sim$ 0.1--12 $M_\odot$, luminosities $\sim$ 0.001--0.1 K~km s$^{-1}$ pc$^{-2}$, and excitation temperatures $\sim$ 10--50K. These are consistent with the standard GMC based size-line width relationships, showing that the scaling laws extend down to size scales of hundredths of a parsec, and to sub solar-mass condensations. There is however no compelling evidence that the majority of clumplets are undergoing free-fall collapse, nor that they are pressure confined.

HD 100546 is a well-studied Herbig Be star-disk system that likely hosts a close-in companion with compelling observational evidence for an embedded protoplanet at 68 AU. We present ALMA observations of the HD 100546 disk which resolve the gas and dust structure at (sub)mm wavelengths. The CO emission (at 345.795 GHz) originates from an extensive molecular disk (390+/-20 AU in radius) whereas the continuum emission is more compact (230+/-20 AU in radius) suggesting radial drift of the mm-sized grains. The CO emission is similar in extent to scattered light images indicating well-mixed gas and um-sized grains in the disk atmosphere. Assuming azimuthal symmetry, a single-component power-law model cannot reproduce the continuum visibilities. The visibilities and images are better reproduced by a double-component model: a compact ring with a width of 21 AU centered at 26 AU and an outer ring with a width of 75+/-3 AU centered at 190+/-3 AU. The influence of a companion and protoplanet on the dust evolution is investigated. The companion at 10 AU facilitates the accumulation of mm-sized grains within a compact ring, ~ 20 - 30 AU, by ~ 10 Myr. The injection of a protoplanet at 1 Myr hastens the ring formation (~ 1.2 Myr) and also triggers the development of an outer ring (~ 100 - 200 AU). These observations provide additional evidence for the presence of a close-in companion and hint at dynamical clearing by a protoplanet in the outer disk.

A fundamental question about the early evolution of low-mass protostars is when circumstellar disks may form. High angular resolution observations of molecular transitions in the (sub)millimeter wavelength windows make it possible to investigate the kinematics of the gas around newly-formed stars, for example to identify the presence of rotation and infall. IRAS 16293-2422 was observed with the extended Submillimeter Array (eSMA) resulting in subarcsecond resolution (0.46" x 0.29", i.e. $\sim$ 55 $\times$ 35~AU) images of compact emission from the C$^{17}$O (3-2) and C$^{34}$S (7-6) transitions at 337~GHz (0.89~mm). To recover the more extended emission we have combined the eSMA data with SMA observations of the same molecules. The emission of C$^{17}$O (3-2) and C$^{34}$S (7-6) both show a velocity gradient oriented along a northeast-southwest direction with respect to the continuum marking the location of one of the components of the binary, IRAS16293A. Our combined eSMA and SMA observations show that the velocity field on the 50--400~AU scales is consistent with a rotating structure. It cannot be explained by simple Keplerian rotation around a single point mass but rather needs to take into account the enclosed envelope mass at the radii where the observed lines are excited. We suggest that IRAS 16293-2422 could be among the best candidates to observe a pseudo-disk with future high angular resolution observations.

Protoplanetary disks composed of dust and gas are ubiquitous around young stars and are commonly recognized as nurseries of planetary systems. Their lifetime, appearance, and structure are determined by an interplay between stellar radiation, gravity, thermal pressure, magnetic field, gas viscosity, turbulence, and rotation. Molecules and dust serve as major heating and cooling agents in disks. Dust grains dominate the disk opacities, reprocess most of the stellar radiation, and shield molecules from ionizing UV/X-ray photons. Disks also dynamically evolve by building up planetary systems which drastically change their gas and dust density structures. Over the past decade significant progress has been achieved in our understanding of disk chemical composition thanks to the upgrade or advent of new millimeter/Infrared facilities (SMA, PdBI, CARMA, Herschel, e-VLA, ALMA). Some major breakthroughs in our comprehension of the disk physics and chemistry have been done since PPV. This review will present and discuss the impact of such improvements on our understanding of the disk physical structure and chemical composition.

(Abridged) Disks are observed around pre-main sequence stars, but how and when they form is still heavily debated. While disks around young stellar objects have been identified through thermal dust emission, spatially and spectrally resolved molecular line observations are needed to determine their nature. We present subarcsecond observations of dust and gas toward four Class I low-mass young stellar objects in Taurus. The 13CO and C18O J=2-1 transitions at 220 GHz were observed with the Plateau de Bure Interferometer at a spatial resolution of ~0.8'' and analyzed using uv-space position velocity diagrams to determine the nature of their observed velocity radient. Rotationally supported disks (RSDs) are detected around 3 of the 4 Class I sources studied. The derived masses identify them as Stage I objects; i.e., their stellar mass is higher than their envelope and disk masses. The outer radii of the Keplerian disks toward our sample of Class I sources are <~ 100 AU. The lack of on-source C18O emission for TMR1 puts an upper limit of 50 AU on its size. Flattened structures at radii > 100 AU around these sources are dominated by infalling motion (v propto r^-1). A large-scale envelope model is required to estimate the basic parameters of the flattened structure from spatially resolved continuum data. Similarities and differences between the gas and dust disk are discussed. Combined with literature data, the sizes of the RSDs around Class I objects are best described with evolutionary models with an initial rotation of 10^-14 Hz and slow sound speeds. Based on the comparison of gas and dust disk masses, little CO is frozen out within 100 AU in these disks. RSDs with radii up to 100 AU are present around Class I embedded objects. Larger surveys of both Class 0 and I objects are needed to determine whether most disks form late or early in the embedded phase.

Abridged abstract: For stars to form, material must fall inwards from core scales through the envelope towards the central protostar. The velocity profile around protostars is poorly constrained. 6 Class 0 protostars and one Class I protostars observed with HIFI on board Herschel as part of the "Water in star-forming regions with Herschel" (WISH) survey show infall signatures in water line observations. 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 these sources. We assume a free-fall velocity profile and, having found the best fit, vary the radii over which infall takes place. In the well-studied Class 0 protostar NGC1333-IRAS4A we find that infall takes place over the whole envelope to which our observations are sensitive (r>~1000 AU). For 4 sources infall takes place on core to envelope scales (i.e. ~10000-3000 AU). In 2 sources 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 6/7 protostars is further evidence that water must be heavily depleted from the gas-phase at these radii. Infall in four of the sources is supersonic and infall 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 IRAS4A is large (>~10^-4 msol\yr), higher than the mass outflow rate and expected mass accretion rates onto the star, suggesting that any flattened disk-like structure on small scales will be gravitationally unstable, potentially leading to rotational fragmentation and/or episodic accretion.

Herschel/HIFI spectroscopic observations of CO J=10-9, CO J=16-15 and [CII] towards HD 100546 are presented. The objective is to resolve the velocity profile of the lines to address the emitting region of the transitions and directly probe the distribution of warm gas in the disk. The spectra reveal double-peaked CO line profiles centered on the systemic velocity, consistent with a disk origin. The J=16-15 line profile is broader than that of the J=10-9 line, which in turn is broader than those of lower J transitions (6-5, 3-2, observed with APEX), thus showing a clear temperature gradient of the gas with radius. A power-law flat disk model is used to fit the CO line profiles and the CO rotational ladder simultaneously, yielding a temperature of T_0=1100 \pm 350 K (at r_0 = 13 AU) and an index of q=0.85 \pm 0.1 for the temperature radial gradient. This indicates that the gas has a steeper radial temperature gradient than the dust (mean q_{dust} ~ 0.5), providing further proof of the thermal decoupling of gas and dust at the disk heights where the CO lines form. The [CII] line profile shows a strong single-peaked profile red-shifted by 0.5 km s-1 compared to the systemic velocity. We conclude that the bulk of the [CII] emission has a non-disk origin (e.g., remnant envelope or diffuse cloud).

[Context] Two competing models describe the formation of massive stars in objects like the Orion Trapezium. In the turbulent core accretion model, the resulting stellar masses are directly related to the mass distribution of the cloud condensations. In the competitive accretion model, the gravitational potential of the protocluster captures gas from the surrounding cloud for which the individual cluster members compete. [Aims] With high resolution submillimeter observations of the structure, kinematics, and chemistry of the proto-Trapezium cluster W3 IRS5, we aim to determine which mode of star formation dominates. [Methods] We present 354 GHz Submillimeter Array observations at resolutions of 1"-3" (1800-5400 AU) of W3 IRS5. ...... [Results] The observations show five emission peaks (SMM1-5). SMM1 and SMM2 contain massive embedded stars (~20 Msun); SMM3-5 are starless or contain low-mass stars (<8 Msun). The inferred densities are high, >= 10^7 cm^-3, but the core masses are small, 0.2-0.6 Msun. The detected molecular emission reveals four different chemical zones. ...... [Conclusions] The proto-Trapezium cluster W3 IRS5 is an ideal test case to discriminate between models of massive star formation. Either the massive stars accrete locally from their local cores; in this case the small core masses imply that W3 IRS5 is at the very end stages (1000 yr) of infall and accretion, or the stars are accreting from the global collapse of a massive, cluster forming core. We find that the observed masses, densities and line widths observed toward W3 IRS 5 and the surrounding cluster forming core are consistent with the competitive accretion of gas at rates of Macc~10^-4 Msun yr^-1 by the massive young forming stars. ......

Planets form in the disks around young stars. Their formation efficiency and composition are intimately linked to the protoplanetary disk locations of "snow lines" of abundant volatiles. We present chemical imaging of the CO snow line in the disk around TW Hya, an analog of the solar nebula, using high spatial and spectral resolution Atacama Large Millimeter/Submillimeter Array (ALMA) observations of N2H+, a reactive ion present in large abundance only where CO is frozen out. The N2H+ emission is distributed in a large ring, with an inner radius that matches CO snow line model predictions. The extracted CO snow line radius of ~ 30 AU helps to assess models of the formation dynamics of the Solar System, when combined with measurements of the bulk composition of planets and comets.

The high spatial and line sensitivity of ALMA opens the possibility of resolving emission from molecules in circumstellar disks. With an understanding of physical conditions under which molecules have high abundance, they can be used as direct tracers of distinct physical regions. In particular, DCO+ is expected to have an enhanced abundance within a few Kelvin of the CO freezeout temperature of 19 K, making it a useful probe of the cold disk midplane. We compare ALMA line observations of HD 163296 to a grid of models. We vary the upper- and lower-limit temperatures of the region in which DCO+ is present as well as the abundance of DCO+ in order to fit channel maps of the DCO+ J=5-4 line. To determine the abundance enhancement compared to the general interstellar medium, we carry out similar fitting to HCO+ J=4-3 and H13CO+ J=4-3 observations. ALMA images show centrally peaked extended emission from HCO+ and H13CO+. DCO+ emission lies in a resolved ring from ~110 to 160 AU. The outer radius approximately corresponds to the size of the CO snowline as measured by previous lower resolution observations of CO lines in this disk. The ALMA DCO+ data now resolve and image the CO snowline directly. In the best fitting models, HCO+ exists in a region extending from the 19 K isotherm to the photodissociation layer with an abundance of 3x10^-10 relative to H2. DCO+ exists within the 19-21 K region of the disk with an abundance ratio [DCO+] / [HCO+] = 0.3. This represents a factor of 10^4 enhancement of the DCO+ abundance within this narrow region of the HD 163296 disk. Such a high enhancement has only previously been seen in prestellar cores. The inferred abundances provide a lower limit to the ionization fraction in the midplane of the cold outer disk (approximately greater than 4x10^-10), and suggest the utility of DCO+ as a tracer of its parent molecule H2D+. Abridged

The evolution of deeply embedded young stellar objects (YSOs) has traditionally been traced through dust continuum spectral energy distributions (SEDs), but the use of CO excitation as an evolutionary probe has not yet been explored due to lack of high-J CO observations. The aim is to constrain the physical characteristics (excitation, kinematics, column density) of the warm gas toward 26 low-mass Class 0 and I YSOs using spectrally-resolved Herschel Space Observatory data of high-J lines of CO. Data are complemented by ground-based observations from APEX and the JCMT to compare those with the colder gas traced by lower-J CO lines. This is the first large spectrally resolved high-J CO survey conducted for these types of sources. The median excitation temperatures for 12CO, 13CO and C18O derived from single-temperature fits to the J_u=2-10 integrated intensities are ~70K, 48K and 37K, respectively, with no significant difference between Class 0 and I sources and no trend with M_env or L_bol. Thus, in contrast with the continuum SEDs, the spectral line energy distributions (SLEDs) do not show an evolution during the embedded stage. In contrast, the integrated line intensities of all CO isotopologs show a clear decrease with evolutionary stage as the envelope is dispersed. The H2O 110-101/CO 10-9 intensity ratio does not change significantly with velocity, in contrast with the H2O/CO 3-2 ratio, indicating that CO 10-9 is the lowest transition for which the line wings probe the same warm shocked gas as H2O. Modeling of the full suite of C18O lines indicates an abundance profile for Class 0 sources that is consistent with a freeze-out zone below 25 K and evaporation at higher temperatures, with but some fraction of the CO transformed into other species in the cold phase. In contrast, observations for two Class I sources in Ophiuchus are consistent with a constant high CO abundance profile.

Disk winds have been postulated as a mechanism for angular momentum release in protostellar systems for decades. HD 163296 is a Herbig Ae star surrounded by a disk and has been shown to host a series of HH knots (HH 409) with bow shocks associated with the farthest knots. Here we present ALMA Science Verification data of CO J=2-1 and J=3-2 emission which are spatially coincident with the blue shifted jet of HH knots, and offset from the disk by -18.6 km/s. The emission has a double corkscrew morphology and extends more than 10" from the disk with embedded emission clumps coincident with jet knots. We interpret this double corkscrew as emission from material in a molecular disk wind, and that the compact emission near the jet knots is being heated by the jet which is moving at much higher velocities. We show that the J=3-2 emission is likely heavily filtered by the interferometer, but the J=2-1 emission suffers less due to the larger beam and measurable angular scales. Excitation analysis suggests temperatures exceeding 900 K in these compact features, with the wind mass, momentum and energy being of order 10^{-5} M_sun, 10^{-4} M_sun km/s and 10^{40} erg respectively. The high mass loss rate suggests that this star is dispersing the disk faster than it is funneling mass onto the star.

(Abridged) Far-infrared Herschel-PACS spectra of 18 low-mass protostars of various luminosities and evolutionary stages are studied. We quantify their far-infrared line emission and the contribution of different atomic and molecular species to the gas cooling budget during protostellar evolution. We also determine the spatial extent of the emission and investigate the underlying excitation conditions. Most of the protostars in our sample show strong atomic and molecular far-infrared emission. Water is detected in 17 objects, including 5 Class I sources. The high-excitation H2O line at 63.3 micron is detected in 7 sources. CO transitions from J=14-13 up to 49-48 are found and show two distinct temperature components on Boltzmann diagrams with rotational temperatures of ~350 K and ~700 K. H2O has typical excitation temperatures of ~150 K. Emission from both Class 0 and I sources is usually spatially extended along the outflow direction but with a pattern depending on the species and the transition. The H2O line fluxes correlate strongly with those of the high-J CO lines, as well as with the bolometric luminosity and envelope mass. They correlate less strongly with OH and not with [OI] fluxes. The PACS data probe at least two physical components. The H2O and CO emission likely arises in non-dissociative (irradiated) shocks along the outflow walls with a range of pre-shock densities. Some OH is also associated with this component, likely resulting from H2O photodissociation. UV-heated gas contributes only a minor fraction to the CO emission observed by PACS, based on the strong correlation between the shock-dominated CO 24-23 line and the CO 14-13 line. [OI] and some of the OH emission probe dissociative shocks in the inner envelope. The total far-infrared cooling is dominated by H2O and CO, with [OI] increasing for Class I sources.

Observations of higher-excited transitions of abundant molecules such as CO are important for determining where energy in the form of shocks is fed back into the parental envelope of forming stars. The nearby prototypical and protobinary low-mass hot core, IRAS16293-2422 (I16293) is ideal for such a study. The source was targeted with ALMA for science verification purposes in band 9, which includes CO J=6-5 (E_up/k_B ~ 116 K), at an unprecedented spatial resolution (~0.2", 25 AU). I16293 itself is composed of two sources, A and B, with a projected distance of 5". CO J=6-5 emission is detected throughout the region, particularly in small, arcsecond-sized hotspots, where the outflow interacts with the envelope. The observations only recover a fraction of the emission in the line wings when compared to data from single-dish telescopes, with a higher fraction of emission recovered at higher velocities. The very high angular resolution of these new data reveal that a bow shock from source A coincides, in the plane of the sky, with the position of source B. Source B, on the other hand, does not show current outflow activity. In this region, outflow entrainment takes place over large spatial scales, >~ 100 AU, and in small discrete knots. This unique dataset shows that the combination of a high-temperature tracer (e.g., CO J=6-5) and very high angular resolution observations is crucial for interpreting the structure of the warm inner environment of low-mass protostars.

Here we report that the activity of comet Hale--Bopp ceased between late 2007 and March, 2009, at about 28 AU distance from the Sun. At that time the comet resided at a distance from the Sun that exceeded the freeze-out distance of regular comets by an order of magnitude. A Herschel Space Observatory PACS scan was taken in mid-2010, in the already inactive state of the nucleus. The albedo has been found to be surprisingly large (8.1$\pm$0.9%{}), which exceeds the value known for any other comets. With re-reduction of archive HST images from 1995 and 1996, we confirm that the pre-perihelion albedo resembled that of an ordinary comet, and was smaller by a factor of two than the post-activity albedo. Our further observations with the Very Large Telescope (VLT) also confirmed that the albedo increased significantly by the end of the activity. We explain these observations by proposing gravitational redeposition of icy grains towards the end of the activity. This is plausible for such a massive body in a cold environment, where gas velocity is lowered to the range of the escape velocity. These observations also show that giant comets are not just the upscaled versions of the comets we know but can be affected by processes that are yet to be fully identified.

We present a new publicly available tool (DustPol) aimed to model the polarised thermal dust emission. The module DustPol, which is publicly available, is part of the ARTIST (Adaptable Radiative Transfer Innovations for Submillimetre Telescopes) package, which also offers tools for modelling the polarisation of line emission together with a model library and a Python-based user interface. DustPol can easily manage analytical as well as pre-gridded models to generate synthetic maps of the Stokes I, Q, and U parameters. These maps are stored in FITS format which is straightforwardly read by the data reduction software used, e.g., by the Atacama Large Millimeter Array (ALMA). This turns DustPol into a powerful engine for the prediction of the expected polarisation features of a source observed with ALMA or the Planck satellite as well as for the interpretation of existing submillimetre observations obtained with other telescopes. DustPol allows the parameterisation of the maximum degree of polarisation and we find that, in a prestellar core, if there is depolarisation, this effect should happen at densities of 10^6 cm-3 or larger. We compare a model generated by DustPol with the observational polarisation data of the low-mass Class 0 object NGC 1333 IRAS 4A, finding that the total and the polarised emission are consistent.

Fomalhaut is a young, nearby star that is suspected to harbor an infant planetary system, interspersed with one or more belts of dusty debris. We present far-infrared images obtained with the Herschel Space Observatory with an angular resolution between 5.7 and 36.7 arcsec at wavelengths between 70 and 500 micrometer. The images show the main debris belt in great detail. Even at high spatial resolution, the belt appears smooth. The region in between the belt and the central star is not devoid of material; thermal emission is observed here as well. Also at the location of the star, excess emission is detected. We use a dynamical model together with radiative-transfer tools to derive the parameters of the debris disk. We include detailed models of the interaction of the dust grains with radiation, for both the radiation pressure and the temperature determination. Comparing these models to the spatially resolved temperature information contained in the images allows us to place strong constraints on the presence of grains that will be blown out of the system by radiation pressure. We use this to derive the dynamical parameters of the system. The appearance of the belt points towards a remarkably active system in which dust grains are produced at a very high rate by a collisional cascade in a narrow region filled with dynamically excited planetesimals. Dust particles with sizes below the blow-out size are abundantly present. The equivalent of 2000 one-km-sized comets are destroyed every day, out of a cometary reservoir amounting to 110 Earth masses. From comparison of their scattering and thermal properties, we find evidence that the dust grains are fluffy aggregates, which indicates a cometary origin. The excess emission at the location of the star may be produced by hot dust with a range of temperatures, but may also be due to gaseous free-free emission from a stellar wind.

[Context] Recent detections of disks around young high-mass stars support the idea of massive star formation through accretion rather than coalescence, but the detailed kinematics in the equatorial region of the disk candidates is not well known, which limits our understanding of the accretion process. [Aims] This paper explores the kinematics of the gas around a young massive star with millimeter-wave interferometry to improve our understanding of the formation of massive stars though accretion. [Methods] We use Plateau de Bure interferometric images to probe the environment of the nearby (~1 kpc) and luminous (~20000 Lsun) high-mass (10-16 Msun) young star AFGL 2591-VLA3 in continuum and in lines of HDO, H218O and SO2 in the 115 and 230 GHz bands. Radiative transfer calculations are employed to investigate the kinematics of the source. [Results] At ~0.5" (500 AU) resolution, the line images clearly resolve the velocity field of the central compact source (diameter of ~ 800 AU) and show linear velocity gradients in the northeast-southwest direction. Judging from the disk-outflow geometry, the observed velocity gradient results from rotation and radial expansion in the equatorial region of VLA3. Radiative transfer calculations suggest that the velocity field is consistent with sub-Keplerian rotation plus Hubble-law like expansion. The line profiles of the observed molecules suggest a layered structure, with HDO emission arising from the disk mid-plane, H218O from the warm mid-layer, and SO2 from the upper disk. [Conclusions] We propose AFGL 2591-VLA3 as a new massive disk candidate, with peculiar kinematics. The rotation of this disk is sub-Keplerian, probably due to magnetic braking, while the stellar wind may be responsible for the expansion of the disk. The expansion motion [...]

(Abridged) Water is a key tracer of dynamics and chemistry in low-mass protostars, but spectrally resolved observations have so far been limited in sensitivity and angular resolution. In this first systematic survey of spectrally resolved water emission in low-mass protostellar objects, H2O was observed in the ground-state transition at 557 GHz with HIFI on Herschel in 29 embedded Class 0 and I protostars. Complementary far-IR and sub-mm continuum data (including PACS data from our program) are used to constrain the spectral energy distribution of each source. H2O intensities are compared to inferred envelope and outflow properties and CO 3-2 emission. H2O emission is detected in all objects except one. The line profiles are complex and consist of several kinematic components. The profiles are typically dominated by a broad Gaussian emission feature, indicating that the bulk of the water emission arises in outflows, not the quiescent envelope. Several sources show multiple shock components in either emission or absorption, thus constraining the internal geometry of the system. Furthermore, the components include inverse P-Cygni profiles in 7 sources (6 Class 0, 1 Class I) indicative of infalling envelopes, and regular P-Cygni profiles in 4 sources (3 Class I, 1 Class 0) indicative of expanding envelopes. "Bullets" moving at >50 km/s are seen in 4 Class 0 sources; 3 of these are new detections. In the outflow, the H2O/CO abundance ratio as a function of velocity is nearly the same for all sources, increasing from 10^-3 at <5 km/s to >10^-1 at >10 km/s. The H2O abundance in the outer envelope is low, ~10^-10. The different H2O profile components show a clear evolutionary trend: in the Class 0 sources, emission is dominated by outflow components originating inside an infalling envelope. When the infall diminishes during the Class I phase, the outflow weakens and H2O emission disappears.

The Gould Belt Legacy Survey will map star-forming regions within 500 pc, using HARP (Heterodyne Array Receiver Programme), SCUBA-2 (Submillimetre Common-User Bolometer Array 2) and POL-2 (Polarimeter 2) on the James Clerk Maxwell Telescope (JCMT). This paper describes HARP observations of the J = 3-2 transitions of 13CO and C18O towards Orion A. The 1500-resolution observations cover 5 pc of the Orion filament, including OMC1 (inc. BN-KL and Orion Bar), OMC 2/3 and OMC 4, and allow a comparative study of the molecular gas properties throughout the star-forming cloud. The filament shows a velocity gradient of ~1 km/s /pc between OMC 1, 2 and 3, and high velocity emission is detected in both isotopologues. The Orion Nebula and Bar have the largest masses and line widths, and dominate the mass and energetics of the high velocity material. Compact, spatially resolved emission from CH3CN, 13CH3OH, SO, HCOOCH3, C2H5OH, CH3CHO and CH3OCHO is detected towards the Orion Hot Core. The cloud is warm, with a median excitation temperature of ~24 K; the Orion Bar has the highest excitation temperature gas, at >80 K. The C18O excitation temperature correlates well with the dust temperature (to within 40%). The C18O emission is optically thin, and the 13CO emission is marginally optically thick; despite its high mass, OMC 1 shows the lowest opacities. A virial analysis indicates that Orion A is too massive for thermal or turbulent support, but is consistent with a model of a filamentary cloud that is threaded by helical magnetic fields. The variation of physical conditions across the cloud is reflected in the physical characteristics of the dust cores....continued

NGC 1333 IRAS 4A and IRAS 4B sources are among the best studied Stage 0 low-mass protostars which are driving prominent bipolar outflows. Most studies have so far concentrated on the colder parts (T<30K) of these regions. The aim is to characterize the warmer parts of the protostellar envelope in order to quantify the feedback of the protostars on their surroundings in terms of shocks, UV heating, photodissociation and outflow dispersal. Fully sampled large scale maps of the region were obtained; APEX-CHAMP+ was used for 12CO 6-5, 13CO 6-5 and [CI] 2-1, and JCMT-HARP-B for 12CO 3-2 emissions. Complementary Herschel-HIFI and ground-based lines of CO and its isotopologs, from 1-0 upto 10-9 (Eu/k 300K), are collected at the source positions. Radiative-transfer models of the dust and lines are used to determine temperatures and masses of the outflowing and UV-heated gas and infer the CO abundance structure. Broad CO emission line profiles trace entrained shocked gas along the outflow walls, with typical temperatures of ~100K. At other positions surrounding the outflow and the protostar, the 6-5 line profiles are narrow indicating UV excitation. The narrow 13CO 6-5 data directly reveal the UV heated gas distribution for the first time. The amount of UV-heated and outflowing gas are found to be comparable from the 12CO and 13CO 6-5 maps, implying that UV photons can affect the gas as much as the outflows. Weak [CI] emission throughout the region indicates a lack of CO dissociating photons. Modeling of the C18O lines indicates the necessity of a "drop" abundance profile throughout the envelope where the CO freezes out and is reloaded back into the gas phase, thus providing quantitative evidence for the CO ice evaporation zone around the protostars. The inner abundances are less than the canonical value of CO/H_2=2.7x10^-4, indicating some processing of CO into other species on the grains.

The distributions of deuterated molecules in protoplanetary disks are expected to depend on the molecular formation pathways. We use observations of spatially resolved DCN emission from the disk around TW Hya, acquired during ALMA Science verification with a ~3" synthesized beam, together with comparable DCO+ observations from the Submillimeter Array, to investigate differences in the radial distributions of these species and hence differences in their formation chemistry. In contrast to DCO+, which shows an increasing column density with radius, DCN is better fit by a model that is centrally peaked. We infer that DCN forms at a smaller radii and thus at higher temperatures than DCO+. This is consistent with chemical network model predictions of DCO+ formation from H2D+ at T<30 K and DCN formation from additional pathways involving CH2D+ at higher temperatures. We estimate a DCN/HCN abundance ratio of ~0.017, similar to the DCO+/HCO+ abundance ratio. Deuterium fractionation appears to be efficient at a range of temperatures in this protoplanetary disk. These results suggest caution in interpreting the range of deuterium fractions observed in Solar System bodies, as multiple formation pathways should be taken into account.

Icy bodies may have delivered the oceans to the early Earth, yet little is known about water in the ice-dominated regions of extra-solar planet-forming disks. The Heterodyne Instrument for the Far-Infrared on-board the Herschel Space Observatory has detected emission from both spin isomers of cold water vapor from the disk around the young star TW Hydrae. This water vapor likely originates from ice-coated solids near the disk surface hinting at a water ice reservoir equivalent to several thousand Earth Oceans in mass. The water's ortho-to-para ratio falls well below that of Solar System comets, suggesting that comets contain heterogeneous ice mixtures collected across the entire solar nebula during the early stages of planetary birth.

(Abridged) As the stellar X-ray and UV light penetration of a protoplanetary disk depends sensitively on the dust properties, trace molecular species like HCO+, HCN, and CN are expected to show marked differences from photoprocessing effects as the dust content in the disk evolves. We investigate the evolution of the UV irradiation of the molecular gas in a sample of classical T Tauri stars in Taurus that exhibit a wide range in grain growth and dust settling properties. We obtained HCO+ (J=3-2), HCN (J=3-2), and CN (J=2-1) observations of 13 sources with the JCMT. Our sample has 1.3mm fluxes in excess of 75mJy, indicating the presence of significant dust reservoirs; a range of dust settling as traced through their spectral slopes between 6, 13, and 25 microns; and varying degrees of grain growth as extrapolated from the strength of their 10-micron silicate emission features. We compare the emission line strengths with the sources' continuum flux and infrared features, and use detailed modeling based on two different model prescriptions to compare typical disk abundances for HCO+, HCN, and CN with the gas-line observations for our sample. We detected HCO+ (3-2) toward 6 disks, HCN (3-2) from 0 disks, and CN (2-1) toward 4 disks. For the complete sample, there is no correlation between the gas-line strengths or their ratios and either the sources' dust continuum flux or infrared slopes.

Herschel-HIFI observations of high-J lines (up to J_u=10) of 12CO, 13CO and C18O are presented toward three deeply embedded low-mass protostars in NGC1333. The observations show several energetic components including shocked and quiescent gas. Radiative transfer models are used to quantify the C18O envelope abundance which require a jump in the abundance at an evaporation temperature, T_ev ~25 K, providing new direct evidence of a CO ice evaporation zone around protostars. The abundance in the outermost part of the envelope, X_0, is within the canonical value of 2x10^-4; however the inner abundance, X_in, is found around a factor of 3-5 lower than X_0.

The young high-eccentricity binary DQ Tau exhibits powerful recurring millimeter-band (mm) flaring attributed to collisions between the two stellar magnetospheres near periastron, when the stars are separated by only ~8Rstar. These magnetospheric interactions are expected to have scales and magnetic field strengths comparable to those of large X-ray flares from single pre-main-sequence (PMS) stars observed in the Chandra Orion Ultradeep Project (COUP). To search for X-rays arising from processes associated with colliding magnetospheres, we performed simultaneous X-ray and mm observations of DQ Tau near periastron phase. We report here several results. 1) As anticipated, DQ Tau was caught in a flare state in both mm and X-rays. A single long X-ray flare spanned the entire 16.5 hour Chandra exposure. 2) The inferred morphology, duration, and plasma temperature of the X-ray flare are typical of those of large flares from COUP stars. 3) However, our study provides three lines of evidence that this X-ray flare likely arises from colliding magnetospheres: the chance of capturing a large COUP-like flare within the span of our observation is small; the relative timing of the X-ray and mm flares indicates the Neupert effect and is consistent with a common coronal structure; the size of the emitting coronal structure (4-5Rstar) inferred from our analysis (which is admittedly model-dependent and should be considered with caution) is comparable to half the binary separation. 4) The peak flare X-ray luminosity is in agreement with an estimate of the power dissipated by magnetic reconnection within the framework of a simple model of interacting magnetospheres.

`Water In Star-forming regions with Herschel' (WISH) is a key program on the Herschel Space Observatory designed to probe the physical and chemical structure of young stellar objects using water and related molecules and to follow the water abundance from collapsing clouds to planet-forming disks. About 80 sources are targeted covering a wide range of luminosities and evolutionary stages, from cold pre-stellar cores to warm protostellar envelopes and outflows to disks around young stars. Both the HIFI and PACS instruments are used to observe a variety of lines of H2O, H218O and chemically related species. An overview of the scientific motivation and observational strategy of the program is given together with the modeling approach and analysis tools that have been developed. Initial science results are presented. These include a lack of water in cold gas at abundances that are lower than most predictions, strong water emission from shocks in protostellar environments, the importance of UV radiation in heating the gas along outflow walls across the full range of luminosities, and surprisingly widespread detection of the chemically related hydrides OH+ and H2O+ in outflows and foreground gas. Quantitative estimates of the energy budget indicate that H2O is generally not the dominant coolant in the warm dense gas associated with protostars. Very deep limits on the cold gaseous water reservoir in the outer regions of protoplanetary disks are obtained which have profound implications for our understanding of grain growth and mixing in disks.

Recent observations of the low-mass pre-main sequence, eccentric spectroscopic binaries DQ Tau and V773 Tau A reveal that their millimeter spectrum is occasionally dominated by flares from non-thermal emission processes. The transient activity is believed to be synchrotron in nature, resulting from powerful magnetic reconnection events when the separate magnetic structures of the binary components are capable of interacting and forced to reorganize, typically near periastron. We conducted the first systematic study of the millimeter variability toward a sample of 12 PMS spectroscopic binaries with the aim to characterize the proliferation of flares amongst sources likely to experience similar interbinary reconnection events. The source sample consists of short-period, close-separation binaries that possess either a high orbital eccentricity or a circular orbit. Using the MAMBO2 array on the IRAM 30m telescope, we carried out continuous monitoring at 1.25 mm over a 4-night period during which all of the high-eccentricity binaries approached periastron. We also obtained simultaneous optical VRI measurements, since a strong link is often observed between stellar reconnection events and optical brightenings. UZ Tau E is the only source to be detected at millimeter wavelengths: it exhibited significant variation; it is also the only source to undergo strong simultaneous optical variability. The binary possesses the largest orbital eccentricity in the current sample, a predicted factor in star-star magnetic interaction events. With orbital parameters and variable accretion activity similar to DQ Tau, the millimeter behavior of UZ Tau E draws many parallels to the DQ Tau model for colliding magnetospheres. However, on the basis of our observations alone, we cannot determine whether the variability is repetitive, or if it could also be due to variable free-free emission in an ionized wind.

We present arcsecond-scale Submillimeter Array observations of the CO(3-2) line emission from the disks around the young stars HD 163296 and TW Hya at a spectral resolution of 44 m/s. These observations probe below the ~100 m/s turbulent linewidth inferred from lower-resolution observations, and allow us to place constraints on the turbulent linewidth in the disk atmospheres. We reproduce the observed CO(3-2) emission using two physical models of disk structure: (1) a power-law temperature distribution with a tapered density distribution following a simple functional form for an evolving accretion disk, and (2) the radiative transfer models developed by D'Alessio et al. that can reproduce the dust emission probed by the spectral energy distribution. Both types of models yield a low upper limit on the turbulent linewidth (Doppler b-parameter) in the TW Hya system (<40 m/s), and a tentative (3-sigma) detection of a ~300 m/s turbulent linewidth in the upper layers of the HD 163296 disk. These correspond to roughly <10% and 40% of the sound speed at size scales commensurate with the resolution of the data. The derived linewidths imply a turbulent viscosity coefficient, alpha, of order 0.01 and provide observational support for theoretical predictions of subsonic turbulence in protoplanetary accretion disks.

We present Herschel observations of the water molecule in the massive star-forming region W3 IRS5. The o-H17O 110-101, p-H18O 111-000, p-H2O 22 202-111, p-H2O 111-000, o-H2O 221-212, and o-H2O 212-101 lines, covering a frequency range from 552 up to 1669 GHz, have been detected at high spectral resolution with HIFI. The water lines in W3 IRS5 show well-defined high-velocity wings that indicate a clear contribution by outflows. Moreover, the systematically blue-shifted absorption in the H2O lines suggests expansion, presumably driven by the outflow. No infall signatures are detected. The p-H2O 111-000 and o-H2O 212-101 lines show absorption from the cold material (T ~ 10 K) in which the high-mass protostellar envelope is embedded. One-dimensional radiative transfer models are used to estimate water abundances and to further study the kinematics of the region. We show that the emission in the rare isotopologues comes directly from the inner parts of the envelope (T > 100 K) where water ices in the dust mantles evaporate and the gas-phase abundance increases. The resulting jump in the water abundance (with a constant inner abundance of 10^{-4}) is needed to reproduce the o-H17O 110-101 and p-H18O 111-000 spectra in our models. We estimate water abundances of 10^{-8} to 10^{-9} in the outer parts of the envelope (T < 100 K). The possibility of two protostellar objects contributing to the emission is discussed.

Herschel-HIFI observations of high-J lines (up to J_u=10) of 12CO, 13CO and C18O are presented toward three deeply embedded low-mass protostars, NGC 1333 IRAS 2A, IRAS 4A, and IRAS 4B, obtained as part of the Water In Star-forming regions with Herschel (WISH) key program. The spectrally-resolved HIFI data are complemented by ground-based observations of lower-J CO and isotopologue lines. The 12CO 10-9 profiles are dominated by broad (FWHM 25-30 km s^-1) emission. Radiative transfer models are used to constrain the temperature of this shocked gas to 100-200 K. Several CO and 13CO line profiles also reveal a medium-broad component (FWHM 5-10 km s^-1), seen prominently in H2O lines. Column densities for both components are presented, providing a reference for determining abundances of other molecules in the same gas. The narrow C18O 9-8 lines probe the warmer part of the quiescent envelope. Their intensities require a jump in the CO abundance at an evaporation temperature around 25 K, thus providing new direct evidence for a CO ice evaporation zone around low-mass protostars.

We present a new code for solving the molecular and atomic excitation and radiation transfer problem in a molecular gas and predicting emergent spectra. This code works in arbitrary three dimensional geometry using unstructured Delaunay latices for the transport of photons. Various physical models can be used as input, ranging from analytical descriptions over tabulated models to SPH simulations. To generate the Delaunay grid we sample the input model randomly, but weigh the sample probability with the molecular density and other parameters, and thereby we obtain an average grid point separation that scales with the local opacity. Our code does photon very efficiently so that the slow convergence of opaque models becomes traceable. When convergence between the level populations, the radiation field, and the point separation has been obtained, the grid is ray-traced to produced images that can readily be compared to observations. Because of the high dynamic range in scales that can be resolved using this type of grid, our code is particularly well suited for modeling of ALMA data. Our code can furthermore deal with overlapping lines of multiple molecular and atomic species.

HERSCHEL-HIFI observations of water from the intermediate mass protostar NGC7129 FIRS 2 provide a powerful diagnostic of the physical conditions in this star formation environment. Six spectral settings, covering four H216O and two H218O lines, were observed and all but one H218O line were detected. The four H2 16 O lines discussed here share a similar morphology: a narrower, \approx 6 km/s, component centered slightly redward of the systemic velocity of NGC7129 FIRS 2 and a much broader, \approx 25 km/s component centered blueward and likely associated with powerful outflows. The narrower components are consistent with emission from water arising in the envelope around the intermediate mass protostar, and the abundance of H2O is constrained to \approx 10-7 for the outer envelope. Additionally, the presence of a narrow self-absorption component for the lowest energy lines is likely due to self-absorption from colder water in the outer envelope. The broader component, where the H2O/CO relative abundance is found to be \approx 0.2, appears to be tracing the same energetic region that produces strong CO emission at high J.

Observations of the T Tauri spectroscopic binary DQ Tau in April 2008 captured an unusual flare at 3 mm, which peaked at an observed max flux of 0.5 Jy (about 27x the quiescent value). Here we present follow-up mm observations that demonstrate a periodicity to the phenomenon. While monitoring 3 new periastron encounters, we detect flares within 17.5 hrs (or 4.6%) of the orbital phase of the first reported flare, and we constrain the main emitting region to a stellar height of 3.7-6.8 Rstar. The recorded activity is consistent with the proposed picture for synchrotron emission initiated by a magnetic reconnection event when the two stellar magnetospheres of the highly eccentric (e=0.556) binary are believed to collide near periastron as the stars approach a minimum separation of 8 Rstar (~13 Rsolar). The similar light curve decay profiles allow us to estimate an average flare duration of 30 hrs. Assuming one mm flare per orbit, DQ Tau could spend approximately 8% of its 15.8-d orbital period in an elevated flux state. Our analysis of the mm emission provides an upper limit of 5% on the linear polarization. We discuss the extent to which a severely entangled magnetic field structure and Faraday rotation effects are likely to reduce the observed polarization fraction. We also predict that, for the current picture, the stellar magnetospheres must be misaligned at a significant angle or, alternatively, that the topologies of the outer magnetospheres are poorly described by a well-ordered dipole inside a radius of 7 Rstar. Finally, to investigate whether reorganization of the magnetic field during the interaction affects mass accretion, we also present simultaneous optical (VRI) monitoring, as an established tracer of accretion activity in this system. We find that an accretion event can occur coincident in both time and duration with the synchrotron fallout of a magnetic reconnection event.

Early results from the Herschel Space Observatory revealed the water cation H2O+ to be an abundant ingredient of the interstellar medium. Here we present new observations of the H2O and H2O+ lines at 1113.3 and 1115.2 GHz using the Herschel Space Observatory toward a sample of high-mass star-forming regions to observationally study the relation between H2O and H2O+ . Nine out of ten sources show absorption from H2O+ in a range of environments: the molecular clumps surrounding the forming and newly formed massive stars, bright high-velocity outflows associated with the massive protostars, and unrelated low-density clouds along the line of sight. Column densities per velocity component of H2 O+ are found in the range of 10^12 to a few 10^13 cm-2 . The highest N(H2O+) column densities are found in the outflows of the sources. The ratios of H2O+/H2O are determined in a range from 0.01 to a few and are found to differ strongly between the observed environments with much lower ratios in the massive (proto)cluster envelopes (0.01-0.1) than in outflows and diffuse clouds. Remarkably, even for source components detected in H2O in emission, H2O+ is still seen in absorption.

The Heterodyne Instrument for the Far Infrared (HIFI) onboard the Herschel Space Observatory allows the first observations of light diatomic molecules at high spectral resolution and in multiple transitions. Here, we report deep integrations using HIFI in different lines of hydrides towards the high-mass star forming region AFGL 2591. Detected are CH, CH+, NH, OH+, H2O+, while NH+ and SH+ have not been detected. All molecules except for CH and CH+ are seen in absorption with low excitation temperatures and at velocities different from the systemic velocity of the protostellar envelope. Surprisingly, the CH(JF,P = 3/2_2,- - 1/2_1,+) and CH+(J = 1 - 0, J = 2 - 1) lines are detected in emission at the systemic velocity. We can assign the absorption features to a foreground cloud and an outflow lobe, while the CH and CH+ emission stems from the envelope. The observed abundance and excitation of CH and CH+ can be explained in the scenario of FUV irradiated outflow walls, where a cavity etched out by the outflow allows protostellar FUV photons to irradiate and heat the envelope at larger distances driving the chemical reactions that produce these molecules.

The late stages of evolution of the primordial circumstellar disks surrounding young stars are poorly understood, yet vital to constrain theories of planet formation. We consider basic structural models for the disks around two ~10 Myr-old members of the nearby RCrA association, RX J1842.9-3532 and RX J1852.3-3700. We present new arcsecond-resolution maps of their 230 GHz continuum emission from the Submillimeter Array and unresolved CO(3-2) spectra from the Atacama Submillimeter Telescope Experiment. By combining these data with broadband fluxes from the literature and infrared fluxes and spectra from the catalog of the Formation and Evolution of Planetary Systems (FEPS) Legacy program on the Spitzer Space Telescope, we assemble a multiwavelength data set probing the gas and dust disks. Using the Monte Carlo radiative transfer code RADMC to model simultaneously the SED and millimeter continuum visibilities, we derive basic dust disk properties and identify an inner cavity of radius 16 AU in the disk around RX J1852.3-3700. We also identify an optically thin 5 AU cavity in the disk around RX J1842.9-3532, with a small amount of optically thick material close to the star. The molecular line observations suggest an intermediate disk inclination in RX J1842.9-3532, consistent with the continuum emission. In combination with the dust models, the molecular data allow us to derive a lower CO content than expected, suggesting that the process of gas clearing is likely underway in both systems, perhaps simultaneously with planet formation.

'Water In Star-forming regions with Herschel' (WISH) is a key programme dedicated to studying the role of water and related species during the star-formation process and constraining the physical and chemical properties of young stellar objects. The Heterodyne Instrument for the Far-Infrared (HIFI) on the Herschel Space Observatory observed three deeply embedded protostars in the low-mass star-forming region NGC1333 in several H2-16O, H2-18O, and CO transitions. Line profiles are resolved for five H16O transitions in each source, revealing them to be surprisingly complex. The line profiles are decomposed into broad (>20 km/s), medium-broad (~5-10 km/s), and narrow (<5 km/s) components. The H2-18O emission is only detected in broad 1_10-1_01 lines (>20 km/s), indicating that its physical origin is the same as for the broad H2-16O component. In one of the sources, IRAS4A, an inverse P Cygni profile is observed, a clear sign of infall in the envelope. From the line profiles alone, it is clear that the bulk of emission arises from shocks, both on small (<1000 AU) and large scales along the outflow cavity walls (~10 000 AU). The H2O line profiles are compared to CO line profiles to constrain the H2O abundance as a function of velocity within these shocked regions. The H2O/CO abundance ratios are measured to be in the range of ~0.1-1, corresponding to H2O abundances of ~10-5-10-4 with respect to H2. Approximately 5-10% of the gas is hot enough for all oxygen to be driven into water in warm post-shock gas, mostly at high velocities.

We performed a sensitive search for the ground-state emission lines of ortho- and para-water vapor in the DM Tau protoplanetary disk using the Herschel/HIFI instrument. No strong lines are detected down to 3sigma levels in 0.5 km/s channels of 4.2 mK for the 1_{10}--1_{01} line and 12.6 mK for the 1_{11}--0_{00} line. We report a very tentative detection, however, of the 1_{10}--1_{01} line in the Wide Band Spectrometer, with a strength of T_{mb}=2.7 mK, a width of 5.6 km/s and an integrated intensity of 16.0 mK km/s. The latter constitutes a 6sigma detection. Regardless of the reality of this tentative detection, model calculations indicate that our sensitive limits on the line strengths preclude efficient desorption of water in the UV illuminated regions of the disk. We hypothesize that more than 95-99% of the water ice is locked up in coagulated grains that have settled to the midplane.

SWAS and Odin provided stringent upper limits on the gas phase water abundance of dark clouds (x(H2O) < 7x10^-9). We investigate the chemistry of water vapor in starless cores beyond the previous upper limits using the highly improved angular resolution and sensitivity of Herschel and measure the abundance of water vapor during evolutionary stages just preceding star formation. High spectral resolution observations of the fundamental ortho water (o-H2O) transition (557 GHz) were carried out with Herschel HIFI toward two starless cores: B68, a Bok globule, and L1544, a prestellar core embedded in the Taurus molecular cloud complex. The rms in the brightness temperature measured for the B68 and L1544 spectra is 2.0 and 2.2 mK, respectively, in a velocity bin of 0.59 km s^-1. The continuum level is 3.5+/-0.2 mK in B68 and 11.4+/-0.4 mK in L1544. No significant feature is detected in B68 and the 3 sigma upper limit is consistent with a column density of o-H2O N(o-H2O) < 2.5x10^13 cm^-2, or a fractional abundance x(o-H2O) < 1.3x10^-9, more than an order of magnitude lower than the SWAS upper limit on this source. The L1544 spectrum shows an absorption feature at a 5 sigma level from which we obtain the first value of the o-H2O column density ever measured in dark clouds: N(o-H2O) = (8+/-4)x10^12 cm^-2. The corresponding fractional abundance is x(o-H2O) ~ 5x10^-9 at radii > 7000 AU and ~2x10^-10 toward the center. The radiative transfer analysis shows that this is consistent with a x(o-H2O) profile peaking at ~10^-8, 0.1 pc away from the core center, where both freeze-out and photodissociation are negligible. Herschel has provided the first measurement of water vapor in dark regions. Prestellar cores such as L1544 (with their high central densities, strong continuum, and large envelopes) are very promising tools to finally shed light on the solid/vapor balance of water in molecular clouds.

The Gould Belt Legacy Survey (GBS) on the JCMT has observed a region of 260 square arcminutes in 12CO J=3--2 emission, and a 190 square arcminute subset of this in 13CO and C18O towards the Serpens molecular cloud. We examine the global velocity structure of the non-outflowing gas, and calculate excitation temperatures and opacities. The large scale mass and energetics of the region are evaluated, with special consideration for high velocity gas. We find the cloud to have a mass of 203 solar masses, and to be gravitationally bound, and that the kinetic energy of the outflowing gas is approximately seventy percent of the turbulent kinetic energy of the cloud. We identify compact outflows towards some of the submillimetre Class 0/I sources in the region

Aims: We present preliminary results of the first Herschel spectroscopic observations of NGC7129 FIRS2, an intermediate mass star-forming region. We attempt to interpret the observations in the framework of an in-falling spherical envelope. Methods: The PACS instrument was used in line spectroscopy mode (R=1000-5000) with 15 spectral bands between 63 and 185 microns. This provided good detections of 26 spectral lines seen in emission, including lines of H2O, CO, OH, O I, and C II. Results: Most of the detected lines, particularly those of H2O and CO, are substantially stronger than predicted by the spherical envelope models, typically by several orders of magnitude. In this paper we focus on what can be learned from the detected CO emission lines. Conclusions: It is unlikely that the much stronger than expected line emission arises in the (spherical) envelope of the YSO. The region hot enough to produce such high excitation lines within such an envelope is too small to produce the amount of emission observed. Virtually all of this high excitation emission must arise in structures such as as along the walls of the outflow cavity with the emission produced by a combination of UV photon heating and/or non-dissociative shocks.

In the framework of the Water in Star-forming regions with Herschel (WISH) key program, maps in water lines of several outflows from young stars are being obtained, to study the water production in shocks and its role in the outflow cooling. This paper reports the first results of this program, presenting a PACS map of the o-H2O 179 um transition obtained toward the young outflow L1157. The 179 um map is compared with those of other important shock tracers, and with previous single-pointing ISO, SWAS, and Odin water observations of the same source that allow us to constrain the water abundance and total cooling. Strong H2O peaks are localized on both shocked emission knots and the central source position. The H2O 179 um emission is spatially correlated with emission from H2 rotational lines, excited in shocks leading to a significant enhancement of the water abundance. Water emission peaks along the outflow also correlate with peaks of other shock-produced molecular species, such as SiO and NH3. A strong H2O peak is also observed at the location of the proto-star, where none of the other molecules have significant emission. The absolute 179 um intensity and its intensity ratio to the H2O 557 GHz line previously observed with Odin/SWAS indicate that the water emission originates in warm compact clumps, spatially unresolved by PACS, having a H2O abundance of the order of 10^-4. This testifies that the clumps have been heated for a time long enough to allow the conversion of almost all the available gas-phase oxygen into water. The total water cooling is ~10^-1 Lo, about 40% of the cooling due to H2 and 23% of the total energy released in shocks along the L1157 outflow.

We obtained Herschel PACS and SPIRE images of the thermal emission of the debris disk around the A5V star {\beta} Pic. The disk is well resolved in the PACS filters at 70, 100, and 160 {\mu}m. The surface brightness profiles between 70 and 160 {\mu}m show no significant asymmetries along the disk, and are compatible with 90% of the emission between 70 and 160 {\mu}m originating in a region closer than 200 AU to the star. Although only marginally resolving the debris disk, the maps obtained in the SPIRE 250 - 500 {\mu}m filters provide full-disk photometry, completing the SED over a few octaves in wavelength that had been previously inaccessible. The small far-infrared spectral index ({\beta} = 0.34) indicates that the grain size distribution in the inner disk (<200AU) is inconsistent with a local collisional equilibrium. The size distribution is either modified by non-equilibrium effects, or exhibits a wavy pattern, caused by an under-abundance of impactors which have been removed by radiation pressure.

We present five band imaging of the Vega debris disc obtained using the Herschel Space Observatory. These data span a wavelength range of 70-500 um with full-width half-maximum angular resolutions of 5.6-36.9". The disc is well resolved in all bands, with the ring structure visible at 70 and 160 um. Radial profiles of the disc surface brightness are produced, and a disc radius of 11" (~ 85 AU) is determined. The disc is seen to have a smooth structure thoughout the entire wavelength range, suggesting that the disc is in a steady state, rather than being an ephemeral structure caused by the recent collision of two large planetesimals.

We present far-infrared spectroscopic observations, taken with the Photodetector Array Camera and Spectrometer (PACS) on the Herschel Space Observatory, of the protoplanetary disk around the pre-main-sequence star HD 100546. These observations are the first within the DIGIT Herschel key program, which aims to follow the evolution of dust, ice, and gas from young stellar objects still embedded in their parental molecular cloud core, through the final pre-main-sequence phases when the circumstellar disks are dissipated. Our aim is to improve the constraints on temperature and chemical composition of the crystalline olivines in the disk of HD 100546 and to give an inventory of the gas lines present in its far-infrared spectrum. The 69 \mu\m feature is analyzed in terms of position and shape to derive the dust temperature and composition. Furthermore, we detected 32 emission lines from five gaseous species and measured their line fluxes. The 69 \mu\m emission comes either from dust grains with ~70 K at radii larger than 50 AU, as suggested by blackbody fitting, or it arises from ~200 K dust at ~13 AU, close to the midplane, as supported by radiative transfer models. We also conclude that the forsterite crystals have few defects and contain at most a few percent iron by mass. Forbidden line emission from [CII] at 157 \mu\m and [OI] at 63 and 145 \mu\m, most likely due to photodissociation by stellar photons, is detected. Furthermore, five H2O and several OH lines are detected. We also found high-J rotational transition lines of CO, with rotational temperatures of ~300 K for the transitions up to J=22-21 and T~800 K for higher transitions.

DK Cha is an intermediate-mass star in transition from an embedded configuration to a star plus disk stage. We aim to study the composition and energetics of the circumstellar material during this pivotal stage. Using the Range Scan mode of PACS on the Herschel Space Observatory, we obtained a spectrum of DK Cha from 55 to 210 micron as part of the DIGIT Key Program. Almost 50 molecular and atomic lines were detected, many more than the 7 lines detected in ISO-LWS. Nearly the entire ladder of CO from J=14-13 to 38-37 (E_u/k = 4080 K), water from levels as excited as E_u/k = 843 K, and OH lines up to E_u/k = 290 K were detected. The continuum emission in our PACS SED scan matches the flux expected from a model consisting of a star, a surrounding disk of 0.03 Solar mass, and an envelope of a similar mass, supporting the suggestion that the object is emerging from its main accretion stage. Molecular, atomic, and ionic emission lines in the far-infrared reveal the outflow's influence on the envelope. The inferred hot gas can be photon-heated, but some emission could be due to C-shocks in the walls of the outflow cavity.

'Water in Star-forming regions with Herschel' (WISH) is a Herschel Key Programme aimed at understanding the physical and chemical structure of young stellar objects (YSOs) with a focus on water and related species. The low-mass protostar HH 46 was observed with the Photodetector Array Camera and Spectrometer (PACS) on the Herschel Space Observatory to measure emission in H2O, CO, OH, [OI], and [CII] lines located between 63 and 186 um. The excitation and spatial distribution of emission can disentangle the different heating mechanisms of YSOs, with better spatial resolution and sensitivity than previously possible. Far-IR line emission is detected at the position of the protostar and along the outflow axis. The OH emission is concentrated at the central position, CO emission is bright at the central position and along the outflow, and H2O emission is concentrated in the outflow. In addition, [OI] emission is seen in low-velocity gas, assumed to be related to the envelope, and is also seen shifted up to 170 km/s in both the red- and blue-shifted jets. Envelope models are constructed based on previous observational constraints. They indicate that passive heating of a spherical envelope by the protostellar luminosity cannot explain the high-excitation molecular gas detected with PACS, including CO lines with upper levels at >2500 K above the ground state. Instead, warm CO and H2O emission is probably produced in the walls of an outflow-carved cavity in the envelope, which are heated by UV photons and non-dissociative C-type shocks. The bright OH and [OI] emission is attributed to J-type shocks in dense gas close to the protostar. In the scenario described here, the combined cooling by far-IR lines within the central spatial pixel is estimated to be 2 \times 10-2 L_sun, with 60-80% attributed to J- and C-type shocks produced by interactions between the jet and the envelope.

The disc around the Herbig Ae/Be star HD 100546 is one of the most extensively studied discs in the southern sky. Although there is a wealth of information about its dust content and composition, not much is known about its gas and large scale kinematics. We detect and study the molecular gas in the disc at spatial resolution from 7.7" to 18.9" using the APEX telescope. The lines 12CO J=7-6, J=6-5, J=3-2, 13CO J=3-2 and [C I] 3P2-3P1 are observed, diagnostic of disc temperature, size, chemistry, and kinematics. We use parametric disc models that reproduce the low-J 12CO emission from Herbig~Ae stars and vary the basic disc parameters - temperature, mass and size. Using the molecular excitation and radiative transfer code RATRAN we fit the observed spectral line profiles. Our observations are consistent with more than 0.001 Msun of molecular gas in a disc of approximately 400 AU radius in Keplerian rotation around a 2.5 Msun star, seen at an inclination of 50 degrees. The detected 12CO lines are dominated by gas at 30-70~K. The non-detection of the [C I] line indicates excess ultraviolet emission above that of a B9 type model stellar atmosphere. Asymmetry in the 12CO line emission suggests that one side of the outer disc is colder by 10-20~K than the other, possibly due to a shadow by a warped geometry of the inner disc. Pointing offsets, foreground cloud absorption and asymmetry in the disc extent are excluded scenarios. Efficient heating of the outer disc ensures that low- and high-J 12CO lines are dominated by the outermost disc regions, indicating a 400 AU radius. The 12CO J=6--5 line arises from a disc layer higher above disc midplane, and warmer by 15-20~K than the layer emitting the J=3--2 line. The existing models of discs around Herbig Ae stars, assuming a B9.5 type model stellar atmosphere overproduce the [CI] 3P2--3P1 line intensity from HD 100546 by an order of magnitude.

Observations of Serpens have been performed at the JCMT using Harp-B. Maps over a 4.5'x5.4' region were made in a frequency window around 338 GHz, covering the 7-6 transitions of methanol. Emission is extended over each source, following the column density of H2 but showing up also particularly strongly around outflows. The rotational temperature is low, 15-20 K, and does not vary with position within each source. The abundance is typically 10^-9 - 10^-8 with respect to H2 in the outer envelope, whereas "jumps" by factors of up to 10^2 -10^3 inside the region where the dust temperature exceeds 100 K are not excluded. A factor of up to ~ 10^3 enhancement is seen in outflow gas. In one object, SMM4, the ice abundance has been measured to be ~ 3x10^-5 with respect to H2 in the outer envelope, i.e., a factor of 10^3 larger than the gas-phase abundance. Comparison with C18O J=3-2 emission shows that strong CO depletion leads to a high gas-phase abundance of CH3OH not just for the Serpens sources, but for a larger sample of protostars. The observations illustrate the large-scale, low-level desorption of CH3OH from dust grains, extending out to and beyond 7500 AU from each source, a scenario which is consistent with non-thermal (photo-)desorption from the ice. The observations also illustrate the usefulness of CH3OH as a tracer of energetic input in the form of outflows, where methanol is sputtered from the grain surfaces. Finally, the observations provide further evidence of CH3OH formation through CO hydrogenation proceeding on grain surfaces in low-mass envelopes.

As part of a JCMT Legacy Survey of star formation in the Gould Belt, we present early science results for Taurus. CO J=3-2 maps have been secured along the north-west ridge and bowl, collectively known as L 1495, along with deep 13CO and C18O J=3-2 maps in two sub-regions. With these data we search for molecular outflows, and use the distribution of flows, HH objects and shocked H2 line emission features, together with the population of young stars, protostellar cores and starless condensations to map star formation across this extensive region. In total 21 outflows are identified. It is clear that the bowl is more evolved than the ridge, harbouring a greater population of T Tauri stars and a more diffuse, more turbulent ambient medium. By comparison, the ridge contains a much younger, less widely distributed population of protostars which, in turn, is associated with a greater number of molecular outflows. We estimate the ratio of the numbers of prestellar to protostellar cores in L 1495 to be ~ 1.3-2.3, and of gravitationally unbound starless cores to (gravitationally bound) prestellar cores to be ~ 1. If we take previous estimates of the protostellar lifetime of ~ 5 x 10^5 yrs, this indicates a prestellar lifetime of 9(+/-3) x 10^5 yrs. From the number of outflows we also crudely estimate the star formation efficiency in L 1495, finding it to be compatible with a canonical value of 10-15 %. We note that molecular outflow-driving sources have redder near-IR colours than their HH jet-driving counterparts. We also find that the smaller, denser cores are associated with the more massive outflows, as one might expect if mass build-up in the flow increases with the collapse and contraction of the protostellar envelope.

Methods. We use the recently developed disk code ProDiMo to calculate the physico-chemical structure of protoplanetary disks and apply the Monte-Carlo line radiative transfer code RATRAN to predict observable line profiles and fluxes. We consider a series of Herbig Ae type disk models ranging from 10^-6 M_Sun to 2.2 10^-2 M_Sun (between 0.5 and 700 AU) to discuss the dependency of the line fluxes and ratios on disk mass for otherwise fixed disk parameters. Results. We find the [CII] 157.7 mum line to originate in LTE from the surface layers of the disk, where Tg > Td . The total emission is dominated by surface area and hence depends strongly on disk outer radius. The [OI] lines can be very bright (> 10^-16 W/m^2) and form in slightly deeper and closer regions under non-LTE conditions. The high-excitation [OI] 145.5 mum line, which has a larger critical density, decreases more rapidly with disk mass than the 63.2 mum line. Therefore, the [OI] 63.2 mum/145.5 mum ratio is a promising disk mass indicator, especially as it is independent of disk outer radius for Rout > 200 AU. CO is abundant only in deeper layers A_V >~ 0.05. For too low disk masses (M_disk <~10^-4 M_Sun) the dust starts to become transparent, and CO is almost completely photo-dissociated. For masses larger than that the lines are an excellent independent tracer of disk outer radius and can break the outer radius degeneracy in the [OI] 63.2 mum/[CII]157.7 mum line ratio. Conclusions. The far-IR fine-structure lines of [CII] and [OI] observable with Herschel provide a promising tool to measure the disk gas mass, although they are mainly generated in the atomic surface layers. In spatially unresolved observations, none of these lines carry much information about the inner, possibly hot regions < 30 AU.

Since the main cooling lines of the gas phase are important tracers of the interstellar medium in Galactic and extragalactic sources, proper and detailed understanding of their emission, and the ambient conditions of the emitting gas, is necessary for a robust interpretation of the observations. With high resolution (7"-9") maps (~3x3 pc^2) of mid-J molecular lines we aim to probe the physical conditions and spatial distribution of the warm (50 to few hundred K) and dense gas (n(H_2)>10^5 cm^-3) across the interface region of M17 SW nebula. We have used the dual color multiple pixel receiver CHAMP+ on APEX telescope to obtain a 5'.3x4'.7 map of the J=6-5 and J=7-6 transitions of 12CO, the 13CO J=6-5 line, and the {^3P_2}-{^3P_1} 370 um fine-structure transition of [C I] in M17 SW. LTE and non-LTE radiative transfer models are used to constrain the ambient conditions. The warm gas extends up to a distance of ~2.2 pc from the M17 SW ridge. The 13CO J=6-5 and [C I] 370 um lines have a narrower spatial extent of about 1.3 pc along a strip line at P.A=63 deg. The structure and distribution of the [C I] {^3P_2}-{^3P_1} 370 um map indicate that its emission arises from the interclump medium with densities of the order of 10^3 cm^-3. The warmest gas is located along the ridge of the cloud, close to the ionization front. An LTE approximation indicates that the excitation temperature of the embedded clumps goes up to ~120 K. The non-LTE model suggests that the kinetic temperature at four selected positions cannot exceed 230 K in clumps of density n(H_2)~5x10^5 cm^-3, and that the warm T_k>100 K and dense (n(H_2)>10^4 cm^-3) gas traced by the mid-J 12CO lines represent just about 2% of the bulk of the molecular gas. The clump volume filling factor ranges between 0.04 and 0.11 at these positions.

The Gould Belt Legacy Survey will survey nearby star-forming regions (within 500 pc), using HARP (Heterodyne Array Receiver Programme), SCUBA-2 (Submillimetre Common- User Bolometer Array 2) and POL-2 (Polarimeter 2) on the James Clerk Maxwell Telescope (JCMT). This paper describes the initial data obtained using HARP to observe 12CO, 13CO and C18O J = 3 - 2 towards two regions in Orion B, NGC 2024 and NGC 2071. We describe the physical characteristics of the two clouds, calculating temperatures and opacities utilizing all three isotopologues. We find good agreement between temperatures calculated from CO and from dust emission in the dense, energetic regions. We determine the mass and energetics of the clouds, and of the high-velocity material seen in 12CO emission, and compare the relative energetics of the high- and low-velocity material in the two clouds. We present a CLUMPFIND analysis of the 13CO condensations. The slope of the condensation mass functions, at the high-mass ends, is similar to the slope of the initial mass function.

Observations of dense molecular gas lie at the basis of our understanding of the density and temperature structure of protostellar envelopes and molecular outflows. We aim to characterize the properties of the protostellar envelope, molecular outflow and surrounding cloud, through observations of high excitation molecular lines within a sample of 16 southern sources presumed to be embedded YSOs. Observations of submillimeter lines of CO, HCO+ and their isotopologues, both single spectra and small maps were taken with the FLASH and APEX-2a instruments mounted on APEX to trace the gas around the sources. The HARP-B instrument on the JCMT was used to map IRAS 15398-3359 in these lines. HCO+ mapping probes the presence of dense centrally condensed gas, a characteristic of protostellar envelopes. The rare isotopologues C18O and H13CO+ are also included to determine the optical depth, column density, and source velocity. The combination of multiple CO transitions, such as 3-2, 4-3 and 7-6, allows to constrain outflow properties, in particular the temperature. Archival submillimeter continuum data are used to determine envelope masses. Eleven of the sixteen sources have associated warm and/or dense quiescent as characteristic of protostellar envelopes, or an associated outflow. Using the strength and degree of concentration of the HCO+ 4-3 and CO 4-3 lines as a diagnostic, five sources classified as Class I based on their spectral energy distributions are found not to be embedded YSOs. The C18O 3-2 lines show that for none of the sources, foreground cloud layers are present. Strong molecular outflows are found around six sources, .. (continued in paper)

The key question about early protostellar evolution is how matter is accreted from the large-scale molecular cloud, through the circumstellar disk onto the central star. A sample of 20 Class 0 and I protostars has been observed in continuum at (sub)millimeter wavelengths at high angular resolution with the Submillimeter Array. Using detailed dust radiative transfer models, we have developed a framework for disentangling the continuum emission from the envelopes and disks, and from that estimated their masses. For the Class I sources in the sample, HCO+ 3-2 line emission has furthermore been observed with the Submillimeter Array. Four of these sources show signs of Keplerian rotation, constraining the masses of the central stars. Both Class 0 and I protostars are surrounded by disks with typical masses of about 0.05 M_sun. No evidence is found for a correlation between the disk mass and evolutionary stage of the young stellar objects. This contrasts the envelope mass, which decreases sharply from 1 M_sun in the Class 0 stage to <0.1 M_sun in the Class I stage. Typically, the disks have masses that are 1-10% of the corresponding envelope masses in the Class 0 stage and 20-60% in the Class I stage. For the Class I sources for which Keplerian rotation is seen, the central stars contain 70-98% of the total mass in the star-disk-envelope system, confirming that these objects are late in their evolution through the embedded protostellar stages. Theoretical models tend to overestimate the disk masses relative to the stellar masses in the late Class I stage. The results argue in favor of a picture in which circumstellar disks are formed early during the protostellar evolution and rapidly process material accreted from the larger scale envelope onto the central star.

The origin and heating mechanisms of warm (50<T<200 K) molecular gas in low-mass young stellar objects (YSOs) are strongly debated. Both passive heating of the inner collapsing envelope by the protostellar luminosity as well as active heating by shocks and by UV associated with the outflows or accretion have been proposed. We aim to characterize the warm gas within protosteller objects, and disentangle contributions from the (inner) envelope, bipolar outflows and the quiescent cloud. High-J CO maps (12CO J=6--5 and 7--6) of the immediate surroundings (up to 10,000 AU) of eight low-mass YSOs are obtained with the CHAMP+ 650/850 GHz array receiver mounted on the APEX telescope. In addition, isotopologue observations of the 13CO J=6--5 transition and [C I] 3P_2-3P_1 line were taken. Strong quiescent narrow-line 12CO 6--5 and 7--6 emission is seen toward all protostars. In the case of HH~46 and Ced 110 IRS 4, the on-source emission originates in material heated by UV photons scattered in the outflow cavity and not just by passive heating in the inner envelope. Warm quiescent gas is also present along the outflows, heated by UV photons from shocks. Shock-heated warm gas is only detected for Class 0 flows and the more massive Class I sources such as HH~46. Outflow temperatures, estimated from the CO 6--5 and 3--2 line wings, are ~100 K, close to model predictions, with the exception of the L~1551 IRS 5 and IRAS 12496-7650, for which temperatures <50 K are found. APEX-CHAMP+ is uniquely suited to directly probe a protostar's feedback on its accreting envelope gas in terms of heating, photodissociation, and outflow dispersal by mapping 1'x1' regions in high-J CO and [C I] lines.

The spatial origin and detectability of rotational H2O emission lines from Herbig Ae type protoplanetary disks beyond 70 micron is discussed. We use the recently developed disk code ProDiMo to calculate the thermo-chemical structure of a Herbig Ae type disk and apply the non-LTE line radiative transfer code Ratran to predict water line profiles and intensity maps. The model shows three spatially distinct regions in the disk where water concentrations are high, related to different chemical pathways to form the water: (1) a big water reservoir in the deep midplane behind the inner rim, (2) a belt of cold water around the distant icy midplane beyond the snow-line r>20AU, and (3) a layer of irradiated hot water at high altitudes z/r=0.1...0.3, extending from about 1AU to 30AU, where the kinetic gas temperature ranges from 200K to 1500K. Although region 3 contains only little amounts of water vapour (~3x10^-5 M_Earth), we find this warm layer to be almost entirely responsible for the rotational water emission lines, execpt for the 3 lowest excitation lines. Thus, Herschel will probe first and foremost the conditions and radial extension of region 3, where water is predominantly formed via neutral-neutral reactions and the gas is thermally decoupled from the dust T_gas>T_dust. The observations do not allow for a determination of the snow-line, because the snow-line truncates the radial extension of region 1, whereas the lines originate from region 3. Different line transfer approximations (LTE, escape probability, Monte Carlo) are discussed. A non-LTE treatment is required in most cases, and the results obtained with the escape probability method are found to underestimate the Monte Carlo results by 2%...45%.

Low-mass star formation is described by gravitational collapse of dense cores of gas and dust. At some point during the collapse, a disk is formed around the protostar and the disk will spin up and grow in size as the core contracts because of angular momentum conservation. The question is how early the disk formation process occurs. In this paper we aim to characterize the kinematical state of a deeply embedded, Class 0 young stellar object, NGC1333-IRAS2A, based on high angular resolution (< 1$''$ $\approx$ 200 AU) interferometric observations of HCN and H$^{13}$CN J = 4-3 from the Submillimeter Array, and test whether a circumstellar disk can be detected based on gas kinematic features. We adopt a physical model which has been shown to describe the object well and obtain a fit of a parameterized model of the velocity field, using a two-dimensional axis-symmetric radiation transfer code. The parameterization and fit to the high angular resolution data characterize the central dynamical mass and the ratio of infall velocity to rotation velocity. We find a large amount of infall and very little rotation on all scales. The central object has a relatively low mass of 0.25 M$_\odot$ . As an object with a low stellar mass compared to the envelope mass, we conclude that NGC1333-IRAS2A is consistent with the suggestion that, as a Class 0 object, it represents the earliest stages of star formation. The large amount of infall relative to rotation also suggests that this is a young object. We do however find the need of a central compact component on scales of a few hundred AU based on the continuum data, which suggests that disk formation happens shortly after the initial gravitational collapse. The data do not reveal a distinct velocity field for this 0.1 M$_\odot$ component.

AIMS Our aim is to characterize the size, mass, density and temperature profiles of the protostellar envelope of HH~46 IRS 1 and its surrounding cloud material as well as the effect the outflow has on its environment.METHODS The CHAMP+ and LABOCA arrays on the APEX telescope, combined with lower frequency line receivers, are used to obtain a large continuum map and smaller heterodyne maps in various isotopologues of CO and HCO+. The high-J lines of CO (6--5 and 7--6) and its isotopologues together with [C I] 2--1, observed with CHAMP+, are used to probe the warm molecular gas in the inner few hundred AU and in the outflowing gas. The data are interpreted with continuum and line radiative transfer models. RESULTS Broad outflow wings are seen in CO low- and high-J lines at several positions, constraining the gas temperatures to a constant value of ~100 K along the red outflow axis and to ~60 K for the blue outflow. The derived outflow mass is of order 0.4--0.8 M_sol, significantly higher than previously found. The bulk of the strong high-J CO line emission has a surprisingly narrow width, however, even at outflow positions. These lines cannot be fit by a passively heated model of the HH 46 IRS envelope. We propose that it originates from photon heating of the outflow cavity walls by ultraviolet photons originating in outflow shocks and the accretion disk boundary layers. At the position of the bow shock itself, the UV photons are energetic enough to dissociate CO. The envelope mass of ~5 M_sol is strongly concentrated towards HH 46 IRS with a density power law of -1.8.

We study the distribution and physical properties of molecular gas in the disc around the T Tauri star IM Lup on scales close to 200 AU. We investigate how well the gas and dust distributions compare and work towards a unified disc model that can explain both gas and dust emission. 12CO, 13CO, and C18O J=2-1 line emission, as well as the dust continuum at 1.3 mm, is observed at 1.8" resolution towards IM Lup using the Submillimeter Array. A detailed disc model based on the dust emission is tested against these observations with the aid of a molecular excitation and radiative transfer code. Apparent discrepancies between the gas and dust distribution are investigated by adopting simple modifications to the existing model. The disc is seen at an inclination of 54+/-3 degrees and is in Keplerian rotation around a 0.8-1.6 Msun star. The outer disc radius traced by molecular gas emission is 900 AU, while the dust continuum emission and scattered light images limit the amount of dust present beyond 400 AU and are consistent with the existing model that assumes a 400 AU radius. Our observations require a drastic density decrease close to 400 AU with the vertical gas column density at 900 AU in the range of 5.d20 - 1.d22 cm-2. We derive a gas-to-dust mass ratio of 100 or higher in disc regions beyond 400 AU. Within 400 AU from the star our observations are consistent with a gas-to-dust ratio of 100 but other values are not ruled out.

We present 0.3 arcsec resolution observations of the disk around GM Aurigae with the Submillimeter Array (SMA) at a wavelength of 860 um and with the Plateau de Bure Interferometer at a wavelength of 1.3 mm. These observations probe the distribution of disk material on spatial scales commensurate with the size of the inner hole predicted by models of the spectral energy distribution. The data clearly indicate a sharp decrease in millimeter optical depth at the disk center, consistent with a deficit of material at distances less than ~20 AU from the star. We refine the accretion disk model of Calvet et al. (2005) based on the unresolved spectral energy distribution (SED) and demonstrate that it reproduces well the spatially resolved millimeter continuum data at both available wavelengths. We also present complementary SMA observations of CO J=3-2 and J=2-1 emission from the disk at 2" resolution. The observed CO morphology is consistent with the continuum model prediction, with two significant deviations: (1) the emission displays a larger CO J=3-2/J=2-1 line ratio than predicted, which may indicate additional heating of gas in the upper disk layers; and (2) the position angle of the kinematic rotation pattern differs by 11 +/- 2 degrees from that measured at smaller scales from the dust continuum, which may indicate the presence of a warp. We note that photoevaporation, grain growth, and binarity are unlikely mechanisms for inducing the observed sharp decrease in opacity or surface density at the disk center. The inner hole plausibly results from the dynamical influence of a planet on the disk material. Warping induced by a planet could also potentially explain the difference in position angle between the continuum and CO data sets.

The Ophiuchus clouds, in particular L~1688, are an excellent region to study the embedded phases of star formation, due to the relatively large number of protostars. However, the standard method of finding and characterizing embedded young stellar objects (YSOs) through just their infrared spectral slope does not yield a reliable sample. This may affect the age determinations, often derived from the statistics on the total number of embedded YSOs and pre-main sequence stars within a cloud.Our aim is to characterize the structure of protostellar envelopes on an individual basis and to correctly identify the embedded YSO population of L1688. Spectral maps of the HCO+ J=4--3 and C18O J=3--2 lines using the HARP-B array on the James Clerk Maxwell Telescope and SCUBA 850 micron dust maps are obtained of all sources in the L1688 region with infrared spectral slopes consistent with, or close to, that of embedded YSOs. Selected 350 micron maps obtained with the Caltech Submillimeter Observatory are presented as well. The properties, extent and variation of dense gas, column density and dust on scalesup to 1' are probed at 15" resolution. Using the spatial variation of the gas and dust, together with the intensity of the HCO+ J=4--3 line, we are able to accurately identify the truly embedded YSOs and determine their properties. RESULTS The protostellar envelopes range from 0.05 to 0.5 Msun in mass. The concentration of HCO+ emission (~0.5 to 0.9) is generally higher than that of the dust concentration. Combined with absolute intensities, HCO+ proves to be a better tracer of protostellar envelopes than dust, which can contain disk and cloud contributions. Our total sample of 45 sources, including all previously classified Class I sources, several flat-spectrum sources and some known disks, was re-classified using the ....

The physical conditions in a collapsing cloud can be traced by observations of molecular lines. To correctly interpret these observations the abundance distributions of the observed species need to be derived. The chemistry in a collapsing molecular cloud is not in a steady state as the density and temperature evolve. We therefore need to follow chemical reactions, both in the gas phase and on dust grains, as well as gas-grain interactions, to predict the abundance distributions. Our aim is to model the abundances of molecules, in the gas phase and on grain mantles in the form of ice, from prestellar core collapse to disk formation. We use a 2-dimensional hydrodynamical simulation as a physical model from which we take the density, temperature, and the flow of the gas. Trace particles, moving along with the gas, are used to follow the chemistry during prestellar core collapse and disk formation. The evolution of the abundances and the composition of ices on grain mantles were compared to observations. We also investigated the initial abundances to be adopted in more detailed modeling of protoplanetary disks by following the chemical evolution of trace particles accreting onto the disk. Fractional abundances of HCO+, N2H+, H2CO, HC3N, and CH3OH from our model with grain surface reactions provide a good match to observations, while abundances of CO, CS, SO, HCN, and HNC show better agreement without grain surface reactions. The observed mantle composition of dust grains is best reproduced when we include surface reactions. The initial chemical abundances to be used for detailed modeling of a protoplanetary disk are found to be different from those in dark interstellar clouds. Ices with a binding energy lower than about 1200 K sublimate before accreting onto the disk, while those with a higher binding energy do not.

We report the first science observations and results obtained with the "extended" SMA (eSMA), which is composed of the SMA (Submillimeter Array), JCMT (James Clerk Maxwell Telescope) and CSO (Caltech Submillimeter Observatory). Redshifted absorptions at z=0.886 of CI (^3P_1 - ^3P_0) were observed with the eSMA with an angular resolution of 0.55"x0.22" at 1.1 mm toward the southwestern image of the remarkable lensed quasar PKS 1830-211, but not toward the northeastern component at a separation of ~1". Additionally, SMA observations of CO, 13CO and C18O (all J=4-3) were obtained toward this object: CO was also detected toward the SW component, but none of the isotopologues were. This is the first time [CI] is detected in this object, allowing the first direct determination of relative abundances of neutral atomic carbon to CO in the molecular clouds of a spiral galaxy at z>0.1. The [CI] and CO profiles can be decomposed into two and three velocity components respectively. We derive C/CO column density ratios ranging from <0.5 (representative of dense cores) to ~2.5 (close to translucent clouds values). This could indicate that we are seeing environments with different physical conditions or that we are witnessing chemical evolution of regions where C has not completely been converted into CO.

For several hours on 2008 April 19 the T Tauri spectroscopic binary DQ Tau was observed to brighten, reaching a maximum detected flux of 468 mJy and likely making it (briefly) the brightest object at 3 mm in the Taurus star-forming region. We present the light curve of a rarely before observed millimeter flare originating in the region around a pre-main-sequence star, and the first from a classical T Tauri star. We discuss the properties and nature of the flaring behavior in the context of pulsed accretion flows (the current picture based largely on studies of this object's optically variable spectrum), as well as magnetospheric re-connection models (a separate theory that predicts millimeter flares for close binaries of high orbital eccentricity). We believe that the flare mechanism is linked to the binary orbit, and therefore periodic. DQ Tau makes a strong case for multi-wavelength follow-up studies, performed in parallel, of future flares to help determine whether magnetospheric and dynamical interactions in a proto-binary system are independent.

The eSMA ("extended SMA") combines the SMA, JCMT and CSO into a single facility, providing enhanced sensitivity and spatial resolution owing to the increased collecting area at the longest baselines. Until ALMA early science observing (2011), the eSMA will be the facility capable of the highest angular resolution observations at 345 GHz. The gain in sensitivity and resolution will bring new insights in a variety of fields, such as protoplanetary/transition disks, high-mass star formation, solar system bodies, nearby and high-z galaxies. Therefore the eSMA is an important facility to prepare the grounds for ALMA and train scientists in the techniques. Over the last two years, and especially since November 2006, there has been substantial progress toward making the eSMA into a working interferometer. In particular, (i) new 345-GHz receivers, that match the capabilities of the SMA system, were installed at the JCMT and CSO; (ii) numerous tests have been performed for receiver, correlator and baseline calibrations in order to determine and take into account the effects arising from the differences between the three types of antennas; (iii) first fringes at 345 GHz were obtained on August 30 2007, and the array has entered the science-verification stage. We report on the characteristics of the eSMA and its measured performance at 230 GHz and that expected at 345 GHz. We also present the results of the commissioning and some initial science-verification observations, including the first absorption measurement of the C/CO ratio in a galaxy at z=0.89, located along the line of sight to the lensed quasar PKS1830-211, and on the imaging of the vibrationally excited HCN line towards IRC+10216.

Understanding the gas abundance distribution is essential when tracing star formation using molecular line observations. Changing density and temperature conditions cause gas to freeze-out onto dust grains, and this needs to be taken into account when modeling a collapsing molecular cloud. This study aims to provide a realistic estimate of the CO abundance distribution throughout the collapse of a molecular cloud. We provide abundance profiles and synthetic spectral lines which can be compared to observations. We use a 2D hydrodynamical simulation of a collapsing cloud and subsequent formation of a protoplanetary disk as input for the chemical calculations. From the resulting abundances, synthetic spectra are calculated using a molecular excitation and radiation transfer code. We compare three different methods to calculate the abundance of CO. Our models also consider cosmic ray desorption and the effects of an increased CO binding energy. The resulting abundance profiles are compared to observations from the literature and are found to agree well. The resulting abundance profiles agree well with analytic approximations, and the corresponding line fluxes match observational data. Our developed method to calculate abundances in hydrodynamical simulations should greatly aid in comparing these to observations, and can easily be generalized to include gas-phase reaction networks.

When low-mass stars form, the collapsing cloud of gas and dust goes through several stages which are usually characterized by the shape of their spectral energy distributions. Such classification is based on the cloud morphology only and does not address the dynamical state of the object. In this paper we investigate the initial cloud collapse and subsequent disk formation through the dynamical behavior as reflected in the sub-millimeter spectral emission line profiles. If a young stellar object is to be characterized by its dynamical structure it is important to know how accurately information about the velocity field can be extracted and which observables provide the best description of the kinematics. Of particular interest is the transition from infalling envelope to rotating disk, because this provides the initial conditions for the protoplanetary disk, such as mass and size. We use a hydrodynamical model, describing the collapse of a core and formation of a disk, to produce synthetic observables which we compare to calculated line profiles of a simple parameterized model. Because we know the velocity field from the hydrodynamical simulation we can determine in a quantitative way how well our best-fit parameterized velocity field reproduces the original. We use a molecular line excitation and radiation transfer code to produce spectra of both our hydro dynamical simulation as well as our parameterized model. We find that information about the velocity field can reasonably well be derived by fitting a simple model to either single-dish lines or interferometric data, but preferentially by using a combination of the two. Our result shows that it is possible to establish relative ages of a sample of young stellar objects using this method, independently of the details of the hydrodynamical model.

(Abridged) Classification schemes for YSOs are based on evaluating the degree of dissipation of the surrounding envelope, whose main effects are the extinction of the optical radiation from the central YSO and re-emission in the far-infrared. Since extinction is a property of column density along the line of sight, the presence of a protoplanetary disk may lead to a misclassification when the system is viewed edge-on. We performed radiative transfer calculations, using the axysimmetric 3D radiative transfer codes RADMC and RADICAL, to show the effects of different geometries on the main indicators of YSO evolutionary stage, like the slope of the flux between 2 and 24mum, the bolometric temperature and the optical depth of silicates and ices. We show that for systems viewed at intermediate angles the 'classical' indicators of evolution accurately trace the envelope column density, and they all agree with each other. On the other hand, edge-on system are misclassified for inclinations larger than ~65deg. In particular, silicate emission, typical of pre-main sequence stars with disks, turns into absorption when the disk column density reaches 1e22cm-2, corresponding e.g. to a 5e-3 Msun flaring disk viewed at 64deg. A similar effect is noticed in all the other classification indicators studied alpha, Tbol, and the H2O and CO2 ices absorption strengths. This misclassification has a big impact on the nature of the flat-spectrum sources (alpha ~0), whose number can be explained by simple geometrical arguments without invoking evolution. A reliable classification scheme using a minimal number of observations is constituted by observations of the mm-flux with both a single dish and an interferometer.

We investigate the physical structure of the gas component of the disk around the pre-main-sequence star HD169142. The 13CO and C18O J=2-1 line emission is observed from the disk with 1.4'' resolution using the Submillimeter Array. We adopt the disk physical structure derived from a model which fits the spectral energy distribution of HD169142. We obtain the full three-dimensional information on the CO emission with the aid of a molecular excitation and radiative transfer code. This information is used for the analysis of our observations and previous 12CO J=2-1 and 1.3 mm continuum data. The disk is in Keplerian rotation and seen at an inclination close to 13 deg from face-on. We conclude that the regions traced by different CO isotopologues are distinct in terms of their vertical location within the disk, their temperature and their column densities. With the given disk structure, we find that freeze-out is not efficient enough to remove a significant amount of CO from gas phase. Both observed lines match the model prediction both in flux and in the spatial structure of the emission. Therefore we use our data to derive the 13CO and C18O mass and consequently the 12CO mass using standard isotopic ratios. We constrain the total disk gas mass to (0.6-3.0)x10(-2) Msun. Adopting a maximum dust opacity of 2 cm2 per gram of dust we derive a minimum dust mass of 2.16x10(-4) Msun from the fit to the 1.3 mm data. Comparison of the derived gas and dust mass shows that the gas to dust mass ratio of 100 is only possible under the assumption of a dust opacity of 2 cm2/g and 12CO abundance of 10(-4) with respect to H2. However, our data are also compatible with a gas to dust ratio of 25, with a dust opacity of 1 cm2/g and 12CO abundance of 2x10(-4).

Within low-mass star formation, water vapor plays a key role in the chemistry and energy balance of the circumstellar material. The Herschel Space Observatory will open up the possibility to observe water lines originating from a wide range of excitation energies.Our aim is to simulate the emission of rotational water lines from envelopes characteristic of embedded low-mass protostars. A large number of parameters that influence the water line emission are explored: luminosity, density,density slope and water abundances.Both dust and water emission are modelled using full radiative transfer in spherical symmetry. The temperature profile is calculated for a given density profile. The H2O level populations and emission profiles are in turn computed with a non-LTE line code. The results are analyzed to determine the diagnostic value of different lines, and are compared with existing observations. Lines can be categorized in: (i) optically thick lines, including ground-state lines, mostly sensitive to the cold outer part; (ii) highly excited (E_u>200-250 K) optically thin lines sensitive to the abundance in the hot inner part; and (iii) lines which vary from optically thick to thin depending on the abundances. Dust influences the emission of water significantly by becoming optically thick at the higher frequencies, and by pumping optically thin lines. A good physical model of a source, including a correct treatment of dust, is a prerequisite to infer the water abundance structure and possible jumps at the evaporation temperature from observations. The inner warm (T>100 K) envelope can be probed byhighly-excited lines, while a combination of excited and spectrally resolved ground state lines probes the outer envelope. Observations of H218O lines, although weak, provide even stronger constraints on abundances.

The A star 49 Ceti, at a distance of 61 pc, is unusual in retaining a substantial quantity of molecular gas while exhibiting dust properties similar to those of a debris disk. We present resolved observations of the disk around 49 Ceti from the Submillimeter Array in the J=2-1 rotational transition of CO with a resolution of 1.0x1.2 arcsec. The observed emission reveals an extended rotating structure viewed approximately edge-on and clear of detectable CO emission out to a distance of ~90 AU from the star. No 1.3 millimeter continuum emission is detected at a 3-sigma sensitivity of 2.1 mJy/beam. Models of disk structure and chemistry indicate that the inner disk is devoid of molecular gas, while the outer gas disk between 40 and 200 AU from the star is dominated by photochemistry from stellar and interstellar radiation. We determine parameters for a model that reproduces the basic features of the spatially resolved CO J=2-1 emission, the spectral energy distribution, and the unresolved CO J=3-2 spectrum. We investigate variations in disk chemistry and observable properties for a range of structural parameters. 49 Ceti appears to be a rare example of a system in a late stage of transition between a gas-rich protoplanetary disk and a tenuous, virtually gas-free debris disk.

We present Submillimeter Array (SMA) observations of several deuterated species in the disk around the classical T Tauri star TW Hydrae at arcsecond scales, including detections of the DCN J=3-2 and DCO+ J=3-2 lines, and upper limits to the HDO 3(1,2)-2(2,1), ortho-H2D+ 1(1,0)-1(1,1) and para-D2H+ 1(1,0)-1(0,1) transitions. We also present observations of the HCN J=3-2, HCO+ J=3-2 and H13CO+ J=4-3 lines for comparison with their deuterated isotopologues. We constrain the radial and vertical distributions of various species in the disk by fitting the data using a model where the molecular emission from an irradiated accretion disk is sampled with a 2D Monte Carlo radiative transfer code. We find that the distribution of DCO+ differs markedly from that of HCO+. The D/H ratios inferred change by at least one order of magnitude (0.01 to 0.1) for radii <30 AU to >70 AU and there is a rapid falloff of the abundance of DCO+ at radii larger than 90 AU. Using a simple analytical chemical model, we constrain the degree of ionization, x(e-)=n(e-)/n(H2), to be ~10^-7 in the disk layer(s) where these molecules are present. Provided the distribution of DCN follows that of HCN, the ratio of DCN to HCN is determined to be 1.7\pm0.5 \times 10^-2; however, this ratio is very sensitive to the poorly constrained vertical distribution of HCN. The resolved radial distribution of DCO+ indicates that {\it in situ} deuterium fractionation remains active within the TW Hydrae disk and must be considered in the molecular evolution of circumstellar accretion disks.

We investigate the apparent discrepancy between gas and dust outer radii derived from millimeter observations of protoplanetary disks. Using 230 and 345 GHz continuum and CO J=3-2 data from the Submillimeter Array for four nearby disk systems (HD 163296, TW Hydrae, GM Aurigae, and MWC 480), we examine models of circumstellar disk structure and the effects of their treatment of the outer disk edge. We show that for these disks, models described by power laws in surface density and temperature that are truncated at an outer radius are incapable of reproducing both the gas and dust emission simultaneously: the outer radius derived from the dust continuum emission is always significantly smaller than the extent of the molecular gas disk traced by CO emission. However, a simple model motivated by similarity solutions of the time evolution of accretion disks that includes a tapered exponential edge in the surface density distribution (and the same number of free parameters) does much better at reproducing both the gas and dust emission. While this analysis does not rule out the disparate radii implied by the truncated power-law models, a realistic alternative disk model, grounded in the physics of accretion, provides a consistent picture for the extent of both the gas and dust.

Circumstellar disks are expected to form early in the process that leads to the formation of a young star, during the collapse of the dense molecular cloud core. It is currently not well understood at what stage of the collapse the disk is formed or how it subsequently evolves. We aim to identify whether an embedded Keplerian protoplanetary disk resides in the L1489 IRS system. Given the amount of envelope material still present, such a disk would respresent a very young example of a protoplanetary disk. Using the Submillimeter Array (SMA) we have observed the HCO$^+$ $J=$ 3--2 line with a resolution of about 1$''$. At this resolution a protoplanetary disk with a radius of a few hundred AUs should be detectable, if present. Radiative transfer tools are used to model the emission from both continuum and line data. We find that these data are consistent with theoretical models of a collapsing envelope and Keplerian circumstellar disk. Models reproducing both the SED and the interferometric continuum observations reveal that the disk is inclined by 40$^\circ$ which is significantly different to the surrounding envelope (74$^\circ$). This misalignment of the angular momentum axes may be caused by a gradient within the angular momentum in the parental cloud or if L1489 IRS is a binary system rather than just a single star. In the latter case, future observations looking for variability at sub-arcsecond scales may be able to constrain these dynamical variations directly. However, if stars form from turbulent cores, the accreting material will not have a constant angular momentum axis (although the average is well defined and conserved) in which case it is more likely to have a misalignment of the angular momentum axes of the disk and the envelope.

This paper describes a James Clerk Maxwell Telescope (JCMT) legacy survey that has been awarded roughly 500 hrs of observing time to be carried out from 2007 to 2009. In this survey we will map with SCUBA-2 (Submillimetre Common User Bolometer Array 2) almost all of the well-known low-mass and intermediate-mass star-forming regions within 0.5 kpc that are accessible from the JCMT. Most of these locations are associated with the Gould Belt. From these observations we will produce a flux-limited snapshot of star formation near the Sun, providing a legacy of images, as well as point-source and extended-source catalogues, over almost 700 square degrees of sky. The resulting images will yield the first catalogue of prestellar and protostellar sources selected by submillimetre continuum emission, and should increase the number of known sources by more than an order of magnitude. We will also obtain CO maps with the array receiver HARP (Heterodyne Array Receiver Programme), in three CO isotopologues, of a large typical sample of prestellar and protostellar sources. We will then map the brightest hundred sources with the SCUBA-2 polarimeter (POL-2), producing the first statistically significant set of polarization maps in the submillimetre. The images and source catalogues will be a powerful reference set for astronomers, providing a detailed legacy archive for future telescopes, including ALMA, Herschel and JWST.

We present Very Large Array observations at 7 millimeters wavelength that resolve the dust emission structure in the disk around the young star TW Hydrae at the scale of the ~4 AU (~0.16") radius inner hole inferred from spectral energy distribution modeling. These high resolution data confirm directly the presence of an inner hole in the dust disk and reveal a high brightness ring that we associate with the directly illuminated inner edge of the disk. The clearing of the inner disk plausibly results from the dynamical effects of a giant planet in formation. In an appendix, we develop an analytical framework for the interpretation of visibility curves from power-law disk models with inner holes.

This paper presents an overview of several modeling tools for analyzing molecular line observations at submillimeter wavelengths. These tools are already proving to be very valuable for the interpretation of data from current telescopes, and will be indispensable for data obtained with ALMA. The tools are: (1) the Leiden Atomic and Molecular DAtabase (LAMDA), a collection of spectroscopic data and collisional excitation rates; (2) RADEX, an on-line and off-line program to calculate non-LTE excitation and emission from a homogeneous medium, based on the escape probability approximation; (3) RATRAN, an accelerated Monte Carlo program to solve molecular excitation and radiative transfer in spherical and cylindrical symmetry. The paper presents examples of how to use these tools in conjunction with existing data reduction packages to quantitatively interpret submillimeter single-dish and interferometric observations. The described tools are publically available at http://www.strw.leidenuniv.nl/~moldata . The paper concludes with a discussion of future needs in the fields of molecular data and radiative transfer.

Stars form in the densest, coldest, most quiescent regions of molecular clouds. Molecules provide the only probes which can reveal the dynamics, physics, chemistry and evolution of these regions, but our understanding of the molecular inventory of sources and how this is related to their physical state and evolution is rudimentary and incomplete. The Spectral Legacy Survey (SLS) is one of seven surveys recently approved by the JCMT Board. Starting in 2007, the SLS will produce a spectral imaging survey of the content and distribution of all the molecules detected in the 345 GHz atmospheric window (between 332 GHz and 373 GHz) towards a sample of 5 sources. Our intended targets are: a low mass core (NGC1333 IRAS4), 3 high mass cores spanning a range of star forming environments and evolutionary states (W49, AFGL2591, and IRAS20126), and a PDR (the Orion Bar). The SLS will use the unique spectral imaging capabilities of HARP-B/ACSIS to study the molecular inventory and the physical structure of these objects, which span different evolutionary stages and physical environments, to probe their evolution during the star formation process. As its name suggests, the SLS will provide a lasting data legacy from the JCMT that is intended to benefit the entire astronomical community. As such, the entire data set (including calibrated spectral datacubes, maps of molecular emission, line identifications, and calculations of the gas temperature and column density) will be publicly available.

Aims: To calculate chemistry and gas temperature of evolving protoplanetary disks with decreasing mass or dust settling, and to explore the sensitivity of gas-phase tracers. Methods: The density and dust temperature profiles for a range of models of flaring and self-shadowed disks around a typical Herbig Ae star are used together with 2-dimensional ultraviolet (UV) radiative transfer to calculate the chemistry and gas temperature. In each model the line profiles and intensities for the fine structure lines of [O I], [C II] and [C I] and the pure rotational lines of CO, CN, HCN and HCO+ are determined. Results: The chemistry shows a strong correlation with disk mass. Molecules that are easily dissociated, like HCN, require high densities and large extinctions before they can become abundant. The products of photodissociation, like CN and C2H, become abundant in models with lower masses. Dust settling mainly affects the gas temperature, and thus high temperature tracers like the O and C+ fine structure lines. The carbon chemistry is found to be very sensitive to the adopted PAH abundance. The line ratios CO/13CO, CO/HCO+ and [O I] 63 um/146 um can be used to distinguish between disks where dust growth and settling takes place, and disks that undergo overall mass loss.

We investigate the velocity transition in the low-mass protostar L1489 IRS, which is known to be embedded in a flattened, disc-like structure that shows both infall and rotation. We construct a model for L1489 IRS consisting of an flattened envelope and a velocity field that can vary from pure infall to pure rotation. We obtain best-fit parameters by comparison to 24 molecular transitions from the literature, and using a molecular excitation code and a Voronoi optimisation algorithm. We test the model against existing millimeter interferometric observations, near-infrared scattered light imaging, and 12CO ro-vibrational lines.We find that L1489 IRS is well described by a central stellar mass of 1.3M$_\odot$ surrounded by a 0.10M$_\odot$ flattened envelope with approximate scale height h\approx 0.57 R, inclined at 74^\circ. The velocity field is strongly dominated by rotation, with the velocity vector making an angle of 15^\circ with the azimuthal direction. Reproducing low-excitation transitions requires that the emission and absorption by the starless core 1' (8400 AU) east of L1489 IRS is included properly, implying that L1489 IRS is located partially behind this core. We speculate that L1489 IRS was originally formed closer to the center of this core, but has migrated to its current position over the past few times 10^5 yr, consistent with their radial velocity difference of 0.4 kms-1. This suggests that L1489 IRS' unusual appearance may be result of its migration, and that it would appear as a `normal' embedded protostar if it were still surrounded by an extended cloud core. Conversely, we hypothesize that the inner envelopes of embedded protostars resemble the rotating structure seen around L1489 IRS.

We investigate the structure of the core surrounding the recently identified deeply embedded young stellar object Barnard 1c which has an unusual polarization pattern as traced in submillimeter dust emission. Barnard 1c lies within the Perseus molecular cloud at a distance of 250 pc. It is a deeply embedded core of 2.4 solar masses (Kirk et al.) and a luminosity of 4 +/- 2 solar luminosities. Observations of CO, 13CO, C18O, HCO+ and N2H+ were obtained with the BIMA array, together with the continuum at 3.3 mm and 2.7 mm. Single-dish measurements of N2H+ and HCO+ with FCRAO reveal the larger scale emission in these lines, The CO and HCO+ emission traces the outflow, which coincides in detail with the S-shaped jet recently found in Spitzer IRAC imaging. The N2H+ emission, which anticorrelates spatially with the C18O emission, originates from a rotating envelope with effective radius ~ 2400 AU and mass 2.1 - 2.9 solar masses. N2H+ emission is absent from a 600 AU diameter region around the young star. The remaining N2H+ emission may lie in a coherent torus of dense material. With its outflow and rotating envelope, B1c closely resembles the previously studied object L483-mm, and we conclude that it is a protostar in an early stage of evolution. We hypothesize that heating by the outflow and star has desorbed CO from grains which has destroyed N2H+ in the inner region and surmise that the presence of grains without ice mantles in this warm inner region can explain the unusual polarization signature from B1c.

We present submillimeter observations of dark clouds that are part of the Spitzer Legacy Program, From Molecular Cores to Planet-Forming Disks (c2d). We used the Submillimetre Common User's Bolometer Array to map the regions observed by Spitzer by the c2d program to create a census of dense molecular cores including data from the infrared to the submillimeter. In this paper, we present the basic data from these observations: maps, fluxes, and source attributes. We also show data for an object just outside the Perseus cloud that was serendipitously observed in our program. We propose that this object is a newly discovered, evolved protostar.

AIMS : To characterize the molecular environment of classical T Tauri stars in Lupus observed with the Spitzer Space Telescope and to search for gas-rich disks toward these sources. METHODS : Submillimeter observations of 12CO, 13CO and C18O in the J=3-2 and 2-1 lines were obtained with the James Clerk Maxwell Telescope toward a sample of 21 T Tauri stars with disks in the Lupus molecular clouds. Pointings at the sources and at selected off-positions are presented in order to disentangle material associated with disks from ambient cloud material. RESULTS : One source, IM Lup (Sz 82), has been found with double-peaked 12CO and 13CO profiles characteristic of a large rotating gas disk. The inclination of the disk is found to be ~20 degrees, with an outer radius of 400-700 AU. For most other sources, including GQ Lup with its substellar companion, the single-dish 12CO lines are dominated by extended cloud emission with a complex velocity structure. No evidence for molecular outflows is found. Compact C18O emission due to dense circumstellar material has been detected toward only two sources. Future searches for gas-rich disks in Lupus should either use interferometers or perform very deep single-dish integrations in dense gas tracers to separate the cloud and disk emission.

The presence of H2D+ in dense cloud cores underlies ion-molecule reactions that strongly enhance the deuterium fractionation of many molecular species. We determine the H2D+ abundance in one starless core, Barnard 68, that has a particularly well established physical, chemical, and dynamical structure. We observed the ortho-H2D+ ground-state line 1_10-1_11, the N2H+ J=4-3 line, and the H13CO+ 4-3 line with the APEX telescope. We report the probable detection of the o-H2D+ line at an intensity Tmb=0.22 +- 0.08 K and exclusively thermal line width, and find only upper limits to the N2H+ 4-3 and H13CO+ 4-3 intensities. Within the uncertainties in the chemical reaction rates and the collisional excitation rates, chemical model calculations and excitation simulations reproduce the observed intensities and that of o-H2D+ in particular.

Aims:To obtain insight into the physical structure of the warm gas in the inner envelope of protostars and the interaction with the outflow. Methods: Sub-millimeter observations of 12CO, 13CO and/or C18O in J=3-2, J=4-3 and J=7-6 were obtained with the APEX Telescope towards IRAS 12496-7650, an intermediate mass young stellar object. The data are compared to ISO-LWS observations of CO J=14-13 up to J=19-18 lines to test the different proposed origins of the CO lines. Results: The outflow is prominently detected in the 3-2 and 4-3 lines, but not seen at similar velocities in the 7-6 line, constraining the temperature in the high-velocity (> 5 km/s from line center) gas to less than 50 K, much lower than inferred from the analysis of the ISO-LWS data. In addition, no isothermal gas model can reproduce the emission in both the 7-6 and the higher-J ISO-LWS lines.The 7-6 line probably originates in the inner (<250 AU) region of the envelope at ~150 K. Detailed radiative transfer calculations suggest that the ISO-LWS lines are excited by a different mechanism, possibly related to the larger-scale outflow. All possible mechanisms on scales smaller than 8 " are excluded. High-resolution continuum as well as high-J 12CO and isotopic line mapping are needed to better constrain the structure of the warm gas in the inner envelope and the interaction with the outflow.

We present an interferometric search for large molecules, including methanol, methyl cyanide, ethyl cyanide, ethanol, and methyl formate in comets LINEAR (C/2002 T7) and NEAT (C/2001 Q4) with the Berkeley-Illinois-Maryland Association (BIMA) array. In addition, we also searched for transitions of the simpler molecules CS, SiO, HNC, HN13C and 13CO . We detected transitions of methanol and CS around Comet LINEAR and one transition of methanol around Comet NEAT within a synthesized beam of ~20''. We calculated the total column density and production rate of each molecular species using the variable temperature and outflow velocity (VTOV) model described by Friedel et al.(2005).Considering the molecular production rate ratios with respect to water, Comet T7 LINEAR is more similar to Comet Hale-Bopp while Comet Q4 NEAT is more similar to Comet Hyakutake. It is unclear, however, due to such a small sample size, whether there is a clear distinction between a Hale-Bopp and Hyakutake class of comet or whether comets have a continuous range of molecular production rate ratios.

We present Submillimeter Array observations of the Herbig Ae star HD169142 in 1.3 millimeter continuum emission and 12CO J=2-1 line emission at 1.5 arcsecond resolution that reveal a circumstellar disk. The continuum emission is centered on the star position and resolved, and provides a mass estimate of about 0.02 solar masses for the disk. The CO images show patterns in position and velocity that are well matched by a disk in Keplerian rotation with low inclination to the line-of-sight. We use radiative transfer calculations based on a flared, passive disk model to constrain the disk parameters by comparison to the spectral line emission. The derived disk radius is 235 AU, and the inclination is 13 degrees. The model also necessitates modest depletion of the CO molecules, similar to that found in Keplerian disks around T Tauri stars.

We present the first images of the 691.473 GHz CO J=6-5 line in a protoplanetary disk, obtained along with the 690 GHz dust continuum, toward the classical T Tauri star TW Hya using the Submillimeter Array. Imaging in the CO J=6-5 line reveals a rotating disk, consistent with previous observations of CO J=3-2 and 2-1 lines. Using an irradiated accretion disk model and 2D Monte Carlo radiative transfer, we find that additional surface heating is needed to fit simultaneously the absolute and relative intensities of the CO J=6-5, 3-2 and 2-1 lines. In particular, the vertical gas temperature gradient in the disk must be steeper than that of the dust, mostly likely because the CO emission lines probe nearer to the surface of the disk. We have used an idealized X-ray heating model to fit the line profiles of CO J=2-1 and 3-2 with Chi-square analysis, and the prediction of this model yields CO J=6-5 emission consistent with the observations.

Spitzer Space Telescope spectra of the low mass young stellar object (YSO) IRS 46 (L_bol ~ 0.6 L_sun) in Ophiuchus reveal strong vibration-rotation absorption bands of gaseous C2H2, HCN, and CO2. This is the only source out of a sample of ~100 YSO's that shows these features and the first time they are seen in the spectrum of a solar-mass YSO. Analysis of the Spitzer data combined with Keck L- and M-band spectra gives excitation temperatures of > 350 K and abundances of 10(-6)-10(-5) with respect to H2, orders of magnitude higher than those found in cold clouds. In spite of this high abundance, the HCN J=4-3 line is barely detected with the James Clerk Maxwell Telescope, indicating a source diameter less than 13 AU. The (sub)millimeter continuum emission and the absence of scattered light in near-infrared images limits the mass and temperature of any remnant collapse envelope to less than 0.01 M_sun and 100 K, respectively. This excludes a hot-core type region as found in high-mass YSO's. The most plausible origin of this hot gas rich in organic molecules is in the inner (<6 AU radius) region of the disk around IRS 46, either the disk itself or a disk wind. A nearly edge-on 2-D disk model fits the spectral energy distribution (SED) and gives a column of dense warm gas along the line of sight that is consistent with the absorption data. These data illustrate the unique potential of high-resolution infrared spectroscopy to probe organic chemistry, gas temperatures and kinematics in the planet-forming zones close to a young star.

The structure and chemistry of protoplanetary disks depends strongly on the nature of the central star around which it has formed. The dust temperature is mainly set by the stellar luminosity, while the chemistry of the upper disk layers depends on the amount of intercepted UV and X-ray flux. We will study the differences in chemistry, thermal sturcture and line emission around Herbig Ae/Be, T Tauri stars and low mass M dwarfs. Predictions will be made for future observations with SOFIA and Herschel.

We present interferometric detections of HCN in comets LINEAR (C/2002 T7) and NEAT (C/2001 Q4) with the Berkeley-Illinois-Maryland Association (BIMA) Array. With a 25.4"x20.3" synthesized beam around Comet LINEAR and using a variable temperature and outflow velocity (VTOV) model, we found an HCN column density of N_T=6.4+-2.1x10^12 cm^-2, and a production rate of Q(HCN)=6.5+-2.2x10^26 s^-1, giving a production rate ratio of HCN relative to H_2O of ~3.3+-1.1x10^-3 and relative to CN of ~4.6+-1.5. With a 21.3"x17.5" synthesized beam around Comet NEAT and using a VTOV model, we found an HCN column density of N_T=8.5+-4.5x10^11 cm^-2, and a production rate of Q(HCN)=8.9+-4.7x10^25 s^-1, giving a production rate ratio of HCN relative to H_2O of ~7.4+-3.9x10^-4 and relative to CN of ~0.3+-0.2. For both comets, the production rates relative to H_2O are similar to those found in previous comet observations. For Comet LINEAR the production rate relative to CN is consistent with HCN being the primary parent species of CN, while for Comet NEAT it is too low for this to be the case.

Sensitive 5-38 um Spitzer Space Telescope (SST) and ground based 3-5 um spectra of the embedded low mass protostars B5 IRS1 and HH46 IRS show deep ice absorption bands superposed on steeply rising mid-infrared continua. The ices likely originate in the circumstellar envelopes. The CO2 bending mode at 15 um is a particularly powerful tracer of the ice composition and processing history. Toward these protostars, this band shows little evidence for thermal processing at temperatures above 50 K. Signatures of lower temperature processing are present in the CO and OCN- bands, however. The observed CO2 profile indicates an intimate mixture with H2O, but not necessarily with CH3OH, in contrast to some high mass protostars. This is consistent with the low CH3OH abundance derived from the ground based L band spectra. The CO2/H2O column density ratios are high in both B5 IRS1 and HH46 IRS (~35%). Clearly, the SST spectra are essential to study ice evolution in low mass protostellar environments, and to eventually determine the relation between interstellar and solar system ices.

We present ~2"-4" aperture synthesis observations of the circumstellar disk surrounding the nearby young star TW Hya in the CO J=2--1 and J=3--2 lines and associated dust continuum obtained with the partially completed Submillimeter Array. The extent and peak flux of the 230 and 345 GHz dust emission follow closely the predictions of the irradiated accretion disk model of Calvet et al. (2002). The resolved molecular line emission extends to a radius of at least 200 AU, the full extent of the disk visible in scattered light, and shows a clear pattern of Keplerian rotation. Comparison of the images with 2D Monte Carlo models constrains the disk inclination angle to 7+/-1 degrees. The CO emission is optically thick in both lines, and the kinetic temperature in the line formation region is ~20K. Substantial CO depletion, by an order of magnitude or more from canonical dark cloud values, is required to explain the characteristics of the line emission.

We present deep images of dust continuum emission at 450, 800, and 850 micron of the dark cloud LDN 1689N which harbors the low-mass young stellar objects (YSOs) IRAS 16293-2422A and B (I16293A and I16293B) and the cold prestellar object I16293E. Toward the positions of I16293A and E we also obtained spectra of CO-isotopomers and deep submillimeter observations of chemically related molecules with high critical densities. To I16293A we report the detection of the HDO 1_01 - 0_00 and H2O 1_10 - 1_01 ground-state transitions as broad self-reversed emission profiles with narrow absorption, and a tentative detection of H2D+ 1_10 - 1_11. To I16293E we detect weak emission of subthermally excited HDO 1_01 - 0_00. Based on this set of submillimeter continuum and line data we model the envelopes around I16293A and E. The density and velocity structure of I16293A is fit by an inside-out collapse model, yielding a sound speed of a=0.7 km/s, an age of t=(0.6--2.5)e4 yr, and a mass of 6.1 Msun. The density in the envelope of I16293E is fit by a radial power law with index -1.0+/-0.2, a mass of 4.4 Msun, and a constant temperature of 16K. These respective models are used to study the chemistry of the envelopes of these pre- and protostellar objects. The [HDO]/[H2O] abundance ratio in the warm inner envelope of I16293A of a few times 1e-4 is comparable to that measured in comets. This supports the idea that the [HDO]/[H2O] ratio is determined in the cold prestellar core phase and conserved throughout the formation process of low-mass stars and planets.

We present results from an observing campaign of the molecular content of the coma of comet C/1999 S4 (LINEAR) carried out jointly with the millimeter-arrays of the Berkeley-Illinois-Maryland Association (BIMA) and the Owens Valley Radio Observatory (OVRO). Using the BIMA array in autocorrelation (`single-dish') mode, we detected weak HCN J=1-0 emission from comet C/1999 S4 (LINEAR) at 14 +- 4 mK km/s averaged over the 143" beam. The three days over which emission was detected, 2000 July 21.9-24.2, immediately precede the reported full breakup of the nucleus of this comet. During this same period, we find an upper limit for HCN 1-0 of 144 mJy/beam km/s (203 mK km/s) in the 9"x12" synthesized beam of combined observations of BIMA and OVRO in cross-correlation (`imaging') mode. Together with reported values of HCN 1-0 emission in the 28" IRAM 30-meter beam, our data probe the spatial distribution of the HCN emission from radii of 1300 to 19,000 km. Using literature results of HCN excitation in cometary comae, we find that the relative line fluxes in the 12"x9", 28" and 143" beams are consistent with expectations for a nuclear source of HCN and expansion of the volatile gases and evaporating icy grains following a Haser model.

Molecular lines hold valuable information on the physical and chemical composition of disks around young stars, the likely progenitors of planetary systems. This invited contribution discusses techniques to calculate the molecular emission (and absorption) line spectrum based on models for the physical and chemical structure of protoplanetary disks. Four examples of recent research illutrate these techniques in practice: matching resolved molecular-line emission from the disk around LkCa15 with theoertical models for the chemistry; evaluating the two-dimensional transfer of ultraviolet radiation into the disk, and the effect on the HCN/CN ratio; far-infrared CO line emission from a superheated disk surface layer; and inward motions in the disk around L1489 IRS.

This paper presents a detailed study of the chemistry in the outflow associated with the low-mass protostar NGC1333-IRAS2A down to 3" (650 AU) scales. Millimeter-wavelength aperture-synthesis observations from the OVRO and BIMA interferometers and (sub)millimeter single-dish observations from the Onsala 20m telescope and CSO are presented. The interaction of the highly collimated protostellar outflow with a molecular condensation ~15000 AU from the central protostar is clearly traced by molecular species such as HCN, SiO, SO, CS, and CH3OH. Especially SiO traces a narrow high velocity component at the interface between the outflow and the molecular condensation. Multi-transition single-dish observations are used to distinguish the chemistry of the shock from that of the molecular condensation and to address the physical conditions therein. Statistical equilibrium calculations reveal temperatures of 20 and 70 K for the quiescent and shocked components, respectively, and densities near 10^6 cm^{-3}. Significant abundance enhancements of two to four orders of magnitude are found in the shocked region for molecules such as CH3OH, SiO and the sulfur-bearing molecules. HCO+ is seen only in the aftermath of the shock consistent with models where it is destroyed through release of H2O from grain mantles in the shock. N2H+ shows narrow lines, not affected by the outflow but rather probing the ambient cloud. Differences in abundances of HCN, H2CO and CS are seen between different outflow regions and are suggested to be related to differences in the atomic carbon abundance. Compared to the warm inner parts of protostellar envelopes, higher abundances of in particular CH3OH and SiO are found in the outflows, which may be related to density differences between the regions.

This paper investigates small-scale structures of dense gas and dust around the low-mass protostellar binary NGC1333-IRAS2 using millimeter-wavelength aperture-synthesis observations from the OVRO and BIMA interferometers. The detected 3 mm continuum emission is consistent with models of the envelope around IRAS2A, based on previously reported submillimeter continuum images, down to the 3", or 500 AU, resolution of the interferometer data. Our data constrain the contribution of an unresolved point source to 22 mJy. Within the accuracy of the parameters describing the envelope model, the point source flux has an uncertainty by up to 25%. We interpret this point source as a cold disk of mass \gtrsim 0.3 M_\odot. The same envelope model also reproduces aperture-synthesis line observations of the optically thin isotopic species C34S and H13CO+. The more optically thick main isotope lines show a variety of components in the protostellar environment: N2H+ is closely correlated with dust concentrations as seen at submillimeter wavelengths and is particularly strong toward the starless core IRAS2C. We hypothesize that N2H+ is destroyed through reactions with CO that is released from icy grains near the protostellar sources IRAS2A and B. CS, HCO+, and HCN have complex line shapes apparently affected by both outflow and infall. In addition to the east-west jet from IRAS2A, a north-south velocity gradient near this source indicates a second, perpendicular outflow. This suggests the presence of a binary companion within 0.3" (65 AU) from IRAS2A as driving source of this outflow. Alternative explanations of the velocity gradient such as rotation in a circumstellar envelope or a single, wide-angle outflow are less likely.

The origin of disks surrounding young stars has direct implications for our understanding of the formation of planetary systems. In the interstellar clouds from which star form, angular momentum is regulated by magnetic fields, preventing the spin up of contracting cores. When ~0.03 pc-sized dense cores decouple from the magnetic field and collapse dynamically, ~10^-3 km/s*pc of specific angular momentum is locked into the system. A viscous accretion disk is one of two possible mechanisms available for the necessary redistribution of angular momentum; the other one is the formation of a multiple stellar system. Recent observational results involving high-angular resolution observations are reviewed: the presence of disks deep inside collapsing envelopes; an accretion shock surrounding a disk; the velocity field in collapsing and slowly rotating envelopes; a possible transitional object, characterized as a large, contracting disk; and the velocity field in disks around T Tauri stars. Observational facilities becoming available over the next several years promise to offer significant progress in the study of the origin of protoplanetary disks.

We report detection of continuum emission at 850 and 450 micron from disks around four Classical T Tauri stars in the MBM 12 (L1457) young association. Using a simple model we infer masses of 0.0014-0.012 M_sun for the disk of LkHa 263 ABC, 0.005-0.021 M_sun for S18 ABab, 0.03-0.18 M_sun for LkHa 264 A, and 0.023-0.23 M_sun for LkHa 262. The disk mass found for LkHa 263 ABC is consistent with the 0.0018 M_sun inferred from the scattered light image of the edge-on disk around component C. Comparison to earlier 13CO line observations indicates CO depletion by up to a factor 300 with respect to dark-cloud values. The spectral energy distributions (SED) suggest grain growth, possibly to sizes of a few hundred micron, but our spatially unresolved data cannot rule out opacity as an explanation for the SED shape. Our observations show that these T Tauri stars are still surrounded by significant reservoirs of cold material at an age of 1-5 Myr. We conclude that the observed differences in disk mass are likely explained by binary separation affecting the initial value. With available accretion rate estimates we find that our data are consistent with theoretical expectations for viscously evolving disks having decreased their masses by ~30%.

Abundances for C18O, CS, NH3, H2CO, C3H2, and N2H+ and upper limits for the abundances of 13CO and HCO+ are derived for gas within the Bok globule B68 using our own BIMA array data, single-dish data from the literature, and Monte Carlo radiative transfer models. B68 has had its density structure well determined, removing a major uncertainty from abundance determinations. All abundances for B68 are lower than those derived for translucent and cold dense clouds, but perhaps only significantly for N2H+, NH3, and C3H2. CS depletion toward the extinction peak of B68 is hinted at by the large offset between the extinction peak and the position of maximum CS line brightness. C18O and N2H+ abundances derived here are consistent with recently determined values at positions observed in common.

This poster presents single-dish and aperture-synthesis observations of the J=1-0 (lambda~3 mm) transitions of HCO+, HCN, and N2H+ towards the Serpens star-forming region. Jets driven by young stars affect the structure and the chemistry of their surrounding cloud, and this work aims to assess the extent to which the emission of these three molecular lines is dominated by such processes. In Serpens I find that N2H+ 1-0 traces the total amount of material, except in two regions slightly ahead of shocks. In contrast, the HCO+ and, especially, HCN emission is dominated by regions impacted by outflows. One previously unknown, strongly shocked region is located ~0.1 pc northwest of the young stellar object SMM 4. There is a marked spatial offset between the peaks in the HCN and the N2H+ emission associated with shocked regions. I construct a simple, qualitative chemical model where the N2H+ emission increases in the magnetic precursor of a C-type shock, while N2H+ is destroyed deeper in the shock as the neutrals heat up and species like HCN and water are released from icy grain mantles. I conclude that N2H+ is a reliable tracer of cloud material, and that unresolved observations of HCO+ and HCN will be dominated by material impacted by outflows.

A new two-dimensional axi-symmetric ultraviolet radiative transfer code is presented, which is used to calculate photodissociation and ionization rates for use in chemistry models of flaring circumstellar disks. Scattering and absorption of photons from the central star and from the interstellar radiation field are taken into account. The molecules are effectively photodissociated in the surface layer of the disk, but can exist in the intermediate, moderately warm layers. A comparison has been made with an approximate 2D ray-tracing method and it was found that the latter underestimates the ultraviolet field and thus the molecular photodissociation rates below the disk surface. The full 2D results show significantly higher abundances of radicals such as CN and C2H than previous work, partly due to the fact that CO is dissociated to greater depths. Results for different stellar radiation fields are also presented. The CN/HCN ratio shows a strong dependence on the stellar spectrum, whereas other ratios such as HCO+/CO show only little variation.

We report results of a sensitive search for cold dust and molecular gas in the disks around 8 T Tauri stars in the high-latitude cloud MBM 12. Interferometric observations of 3 mm continuum emission in 5 fields containing 6 of the objects, and literature values for the remaining two, limit the disk masses to M_disk < 0.04 - 0.09 M_sun (gas+dust), for a gas:dust mass ratio of 100 and a distance of 275 pc. By coadding the 3 mm data of our five fields, we set an upper limit to the average disk mass of <M_disk>(N=5) < 0.03 M_sun. Simultaneous observation of the CS J=2-1 and the N2H+ 1-0 lines show no emission. Single-dish observations of the 13CO 2-1 line limit the disk mass to (5 - 10) * 10^-4 M_sun for a standard CO abundance of 2 * 10^-4. Depletion of CO by up to two orders of magnitude, through freezing out or photodissociation, can reconcile these limits. These mass limits lie within the range found in the Taurus-Auriga and rho Oph star-forming regions (0.001 - 0.3 M_sun), and preclude conclusions about possible decrease in disk mass over the 1 - 2 Myr age range spanned by the latter two regions and MBM 12. Our observations can exclude the presence in MBM 12 of T Tauri stars with relatively bright and massive disks such as T Tau, DG Tau, and GG Tau.

Comparison is made between a number of independent computer programs for radiative transfer in molecular rotational lines. The test models are spherically symmetric circumstellar envelopes with a given density and temperature profile. The first two test models have a simple power law density distribution, constant temperature and a fictive 2-level molecule, while the other two test models consist of an inside-out collapsing envelope observed in rotational transitions of HCO+. For the 2-level molecule test problems all codes agree well to within 0.2%, comparable to the accuracy of the individual codes, for low optical depth and up to 2% for high optical depths (tau=4800). The problem of the collapsing cloud in HCO+ has a larger spread in results, ranging up to 12% for the J=4 population. The spread is largest at the radius where the transition from collisional to radiative excitation occurs. The resulting line profiles for the HCO+ J=4-3 transition agree to within 10%, i.e., within the calibration accuracy of most current telescopes. The comparison project and the results described in this paper provide a benchmark for future code development, and give an indication of the typical accuracy of present day calculations of molecular line transfer.

Abundances for 5 molecules (C18O, CS, NH3, H2CO, and C3H2) and 1 molecular ion (N2H+) and upper limits for the abundances of 1 molecule (13CO) and 1 molecular ion (HCO+) are derived for gas within the Bok globule Barnard 68 (B68). The abundances were determined using our own BIMA millimeter interferometer data and single-dish data gathered from the literature, in conjunction with a Monte Carlo radiative transfer model. Since B68 is the only starless core to have its density structure strongly constrained via extinction mapping, a major uncertainty has been removed from these determinations. All abundances for B68 are lower than those derived for translucent and cold dense clouds, but perhaps only significantly for N2H+, NH3, and C3H2. Depletion of CS toward the extinction peak of B68 is hinted at by the large offset between the extinction peak and the position of maximum CS line brightness. Abundances derived here for C18O and N2H+ are consistent with other, recently determined values at positions observed in common.

A (sub-)millimeter line and continuum study of the class I protostar Elias 29 in the Rho Ophiuchi molecular cloud is presented, whose goals are to understand the nature of this source, and to locate the ices that are abundantly present along this line of sight. Within 15-60" beams, several different components contribute to the line emission. Two different foreground clouds are detected, an envelope/disk system and a dense ridge of HCO+ rich material. The latter two components are spatially separated in millimeter interferometer maps. We analyze the envelope/disk system by using inside-out collapse and flared disk models. The disk is in a relatively face-on orientation (<60 degrees), which explains many of the remarkable observational features of Elias 29, such as its flat SED, its brightness in the near infrared, the extended components found in speckle interferometry observations, and its high velocity molecular outflow. It cannot account for the ices seen along the line of sight, however. A small fraction of the ices is present in a (remnant) envelope of mass 0.12-0.33 Msun, but most of the ices (~70%) are present in cool (T<40 K) quiescent foreground clouds. This explains the observed absence of thermally processed ices (crystallized H2O) toward Elias 29. Nevertheless, the temperatures could be sufficiently high to account for the low abundance of apolar (CO, N2, O2) ices. This work shows that it is crucial to obtain spectrally and spatially resolved information from single-dish and interferometric molecular gas observations in order to determine the nature of protostars and to interpret infrared ISO satellite observations of ices and silicates along a pencil beam.

We explore the infrared M band (4.7 um) spectrum of the class I protostar L1489 IRS in the Taurus Molecular Cloud. This is the highest resolution wide coverage spectrum at this wavelength of a low mass protostar observed to date (R=25,000; Dv=12 km/s). Many narrow absorption lines of gas phase 12CO, 13CO, and C18O are detected, as well as a prominent band of solid 12CO. The gas phase 12CO lines have red shifted absorption wings (up to 100 km/s), likely originating from warm disk material falling toward the central object. The isotopes and the 12CO line wings are successfully fitted with a contracting disk model of this evolutionary transitional object (Hogerheijde 2001). This shows that the inward motions seen in millimeter wave emission lines continue to within ~0.1 AU from the star. The colder parts of the disk are traced by the prominent CO ice band. The band profile results from CO in 'polar' ices (CO mixed with H2O), and CO in 'apolar' ices. At the high spectral resolution, the 'apolar' component is, for the first time, resolved into two distinct components, likely due to pure CO and CO mixed with CO2, O2 and/or N2. The ices have probably experienced thermal processing in the upper disk layer traced by our pencil absorption beam: much of the volatile 'apolar' ices has evaporated and the depletion factor of CO onto grains is remarkably low (~7%). This study shows that high spectral resolution 4.7 um observations provide important and unique information on the dynamics and structure of protostellar disks and the evolution of ices in these disks.

Evidence for a transitional stage in the formation of a low-mass star is reported, intermediate between the fully embedded and the T Tauri phases. Millimeter aperture synthesis observations in the HCO+ J=1-0 and 3-2, HCN 1-0, 13CO 1-0, and C18O 1-0 transitions reveal distinctly different velocity fields around two embedded, low-mass young stellar objects. The 0.6 M(sun) of material around TMC 1 (IRAS 04381+2517) closely follows inside-out collapse in the presence of a small amount of rotation (~3 km/s/pc), while L1489 IRS (IRAS 04016+2610) is surrounded by a 2000 AU radius, flared disk containing 0.02 M(sun). This disk shows Keplerian rotation around a ~0.65 M(sun) star and infall at 1.3 (r/100 AU)^-0.5 km/s, or, equivalently, sub-Keplerian motions around a central object between 0.65 and 1.4 M(sun). Its density is characterized by a radial power law and an exponential vertical scale height. The different relative importance of infall and rotation around these two objects suggests that rotationally supported structures grow from collapsing envelopes over a few times 10^5 yr to sizes of a few thousand AU, and then decrease over a few times 10^4 yr to several hundred AU typical for T Tauri disks. In this scenario, L1489 IRS represents a transitional phase between embedded YSOs and T Tauri stars with disks. The expected duration of this phase of ~5% of the embedded stage is consistent with the current lack of other known objects like L1489 IRS. Alternative explanations cannot explain L1489 IRS's large disk, such as formation from a cloud core with an unusually large velocity gradient or a binary companion that prevents mass accretion onto small scales. It follows that the transfer and dissipation of angular momentum is key to understanding the formation of disks from infalling envelopes.

We present a numerical method and computer code to calculate the radiative transfer and excitation of molecular lines. Formulating the Monte Carlo method from the viewpoint of cells rather than photons allows us to separate local and external contributions to the radiation field. This separation is critical to accurate and fast performance at high optical depths (tau>100). The random nature of the Monte Carlo method serves to verify the independence of the solution to the angular, spatial, and frequency sampling of the radiation field. These features allow use of our method in a wide variety of astrophysical problems without specific adaptations: in any axially symmetric source model and for all atoms or molecules for which collisional rate coefficients are available. Continuum emission and absorption by dust is explicitly taken into account but scattering is neglected. We illustrate these features in calculations of (i) the HCO+ J=1-0 and 3-2 emission from a flattened protostellar envelope with infall and rotation, (ii) the CO, HCO+, CN and HCN emission from a protoplanetary disk and (iii) HCN emission from a high-mass young stellar object, where infrared pumping is important. The program can be used for optical depths up to 1000-10,000, depending on source model. We expect this program to be an important tool in analysing data from present and future infrared and (sub) millimetre telescopes.

This paper examines the density and velocity structure of envelopes around young stellar objects through submillimeter continuum imaging of four objects in Taurus and previously obtained molecular-line data. Observations carried out with the SCUBA on the JCMT at 850 and 450 micron of L1489 IRS, L1535 IRS, L1527 IRS, and TMC 1 reveal ~2000 AU elongated structures embedded in extended envelopes. The density distribution in these envelopes is equally well fit by a radial power-law of index p=1.0-2.0 or with a collapse model such as that of Shu (1997: ApJ, 214, 488). This inside-out collapse model predicts 13CO, C18O, HCO+, and H13CO+ line profiles which closely match observed spectra toward three of our four sources. This shows that the inside-out collapse model offers a good description of YSO envelopes, but also that reliable constraints on its parameters require independent measurements of the density and the velocity structure, e.g., through continuum and line observations. For the remaining source, L1489 IRS, we find that a model consisting of a 2000 AU radius, rotating, disk-like structure better describes the data. Possibly, this source is in transition between the embedded Class I and the optically revealed T Tauri phases. Two apparently starless cores are found at ~10,000 AU from L1489 IRS and L1535 IRS. They are cold, 10-15 K, contain 0.5-3.0 M_sol, and have flat density distributions characterized by a Gaussian of ~10,000 AU FWHM. The proximity of these cores shows that star formation in truly isolated cores is rare even in Taurus.

We present the first detection of the low-lying pure rotational emission lines of H2 from circumstellar disks around T~Tauri stars, using the Short Wavelength Spectrometer on the Infrared Space Observatory. These lines provide a direct measure of the total amount of warm molecular gas in disks. The J=2->0 S(0) line at 28.218 mum and the J=3->1 S(1) line at 17.035 mum have been observed toward the double binary system GG Tau. Together with limits on the J=5->3 S(3) and J=7->5 S(5) lines, the data suggest the presence of gas at T_kin=110+-10 K with a mass of (3.6+-2.0)x10^-3 M_sol (3sigma). This amounts to ~3% of the total gas + dust mass of the circumbinary disk as imaged by millimeter interferometry, but is larger than the estimated mass of the circumstellar disk(s). Possible origins for the warm gas seen in H2 are discussed in terms of photon and wind-shock heating mechanisms of the circumbinary material, and comparisons with model calculations are made.

The study of the chemical evolution of gas and dust from pre-stellar dense cores to circumstellar disks around young stars forms an essential part of understanding star- and planet formation. Throughout the collapse- and protostellar phases, simple and complex molecules are formed, many of which deplete onto cold grains and are eventually incorporated into the icy planetesimals of new solar systems. Tracing this chemical evolution provides a wealth of information, not only about the chemical processing in primitive solar nebulae, but also about physical processes which occur in the immediate surroundings of young stellar objects (YSOs). Here we review the chemical processes which occur in the protostellar environment, and models and observations of the chemical structure of the various stages of star formation. We briefly discuss the way in which molecular abundances are derived from observations, and conclude with two examples: the low- to intermediate mass YSOs in Serpens, and the massive YSOs in W3.

ISO [C II] 158 micron, [O I] 63 micron, and H_2 9 and 17 micron observations are presented of the reflection nebula Ced 201, which is a photon-dominated region illuminated by a B9.5 star with a color temperature of 10,000 K (a cool PDR). In combination with ground based [C I] 609 micron, CO, 13CO, CS and HCO+ data, the carbon budget and physical structure of the reflection nebula are constrained. The obtained data set is the first one to contain all important cooling lines of a cool PDR, and allows a comparison to be made with classical PDRs. To this effect one- and three-dimensional PDR models are presented which incorporate the physical characteristics of the source, and are aimed at understanding the dominant heating processes of the cloud. The contribution of very small grains to the photo-electric heating rate is estimated from these models and used to constrain the total abundance of PAHs and small grains. Observations of the pure rotational H_2 lines with ISO, in particular the S(3) line, indicate the presence of a small amount of very warm, approximately 330 K, molecular gas. This gas cannot be accommodated by the presented models.

Aperture synthesis and single-dish (sub) millimeter molecular lines and continuum observations reveal in great detail the envelope structure of deeply embedded young stellar objects (SMM1, SMM2, SMM3, SMM4) in the densely star-forming Serpens Molecular Cloud. Resolved millimeter continuum emission constrains the density structure to a radial power law with index -2.0 +/- 0.5, and envelope masses of 8.7, 3.0, and 5.3 M_sol for SMM1, SMM3, and SMM4. The core SMM2 does not seem to have a central condensation and may not have formed a star yet. The molecular line observations can be described by the same envelope model, if an additional, small amount of warm (100 K) material is included. This probably corresponds to the inner few hundred AU of the envelope were the temperature is high. In the interferometer beam, the molecular lines reveal the inner regions of the envelopes, as well as interaction of the outflow with the surrounding envelope. Bright HCO+ and HCN emission outlines the cavities, while SiO and SO trace the direct impact of the outflow on ambient gas. Taken together, these observations provide a first comprehensive view of the physical and chemical structure of the envelopes of deeply embedded young stellar objects in a clustered environment on scales between 1000 and 10,000 AU.

The young binary system T Tau was observed with the Owens Valley Millimeter Array at 267 GHz continuum and HCO+ J=3-2 emission at 0.8'' resolution, with the JCMT-CSO single-baseline interferometer in 357 GHz continuum, and with the W. M. Keck Telescope at 4 micron. The 267 GHz emission is unresolved with a flux of 397 +/- 35 mJy, located close to the position of the optical star T Tau N. An upper limit of 100 mJy is obtained toward the infrared companion T Tau S. The 357 GHz continuum emission is unresolved, with a flux of 1.35 +/- 0.68 Jy. HCO+ 3-2 was detected from a 2'' diameter core surrounding T Tau N and S. Both stars are detected ar 4 micron, but there is no evidence of the radio source T Tau R. We propose a model in which T Tau S is intrinsically similar to T Tau N, but is obscured by the outer parts of T Tau N's disk. A fit to the spectral energy distribution (SED) between 21 cm and 1.22 micron is constructed on this basis. Adopting a r^-1 surface density distribution and an exponentially truncated edge, disk masses of 0.04 +/- 0.01 Msol and 6e-5 - 3e-3 Msol are inferred for T Tau N and S, respectively. A 0.005-0.03 Msol circumbinary envelope is also required to fit the millimeter to mid-infrared SED.

Interferometer observations of millimeter-continuum (OVRO) and single-dish observations of HCO+ and H13CO+ J=1-0, 3-2, and 4-3 (JCMT, IRAM 30m) are presented of nine embedded low-mass young stellar objects (YSOs) in Taurus. All nine objects are detected at 3.4 and 2.7 mm, with fluxes of 4-200 mJy, and consist of unresolved (<3 arcsec) point sources, plus, toward about half of the objects, an extended envelope. The point sources likely are circumstellar disks, showing that these are established early in the embedded phase. Literature values of 1.1 mm continuum emission are used to trace the envelopes, carrying 0.001-0.26 M(sol). In HCO+, the 1-0 lines trace the surrounding clouds, while the 3-2 and 4-3 are concentrated toward the sources with intensities well correlated with the envelope flux. An HCO+/H2 abundance of 1.2e-8 is derived. The HCO+ line strengths and envelope fluxes can be fit simultaneously with the simple collapse model of Shu (1977), and related density power laws with slopes p=1-3. As an indicator of the relative evolutionary phase of a YSO, the ratio of HCO+ 3-2 line intensity over bolometric luminosity is proposed, which is roughly proportional to the current ratio of envelope over stellar mass. It is concluded that HCO+ 3-2 and 4-3 are excellent tracers of the early embedded phase of star formation.