The volatile composition of a planet is determined by the inventory of gas and ice in the parent disk. The volatile chemistry in the disk is expected to evolve over time, though this evolution is poorly constrained observationally. We present ALMA observations of C18O, C2H, and the isotopologues H13CN, HC15N, and DCN towards five Class 0/I disk candidates. Combined with a sample of fourteen Class II disks presented in Bergner et al. (2019b), this data set offers a view of volatile chemical evolution over the disk lifetime. Our estimates of C18O abundances are consistent with a rapid depletion of CO in the first ~0.5-1 Myr of the disk lifetime. We do not see evidence that C2H and HCN formation are enhanced by CO depletion, possibly because the gas is already quite under-abundant in CO. Further CO depletion may actually hinder their production by limiting the gas-phase carbon supply. The embedded sources show several chemical differences compared to the Class II stage, which seem to arise from shielding of radiation by the envelope (impacting C2H formation and HC15N fractionation) and sublimation of ices from infalling material (impacting HCN and C18O abundances). Such chemical differences between Class 0/I and Class II sources may affect the volatile composition of planet-forming material at different stages in the disk lifetime.

Star-forming regions show a rich and varied chemistry, including the presence of complex organic molecules - both in the cold gas distributed on large scales, and in the hot regions close to young stars where protoplanetary disks arise. Recent advances in observational techniques have opened new possibilities for studying this chemistry. In particular, the Atacama Large Millimeter/submillimeter Array (ALMA) has made it possible to study astrochemistry down to Solar System size scales, while also revealing molecules of increasing variety and complexity. In this review, we discuss recent observations of the chemistry of star-forming environments, with a particular focus on complex organic molecules, taking context from the laboratory experiments and chemical models that they have stimulated. The key takeaway points are: The physical evolution of individual sources plays a crucial role in their inferred chemical signatures, and remains an important area for observations and models to elucidate. Comparisons of the abundances measured toward different star-forming environments (high-mass versus low-mass, Galactic center versus Galactic disk) reveal a remarkable similarity, an indication that the underlying chemistry is relatively independent of variations in their physical conditions. Studies of molecular isotopologs in star-forming regions provide a link with measurements in our own Solar System, and thus may shed light on the chemical similarities and differences expected in other planetary systems.

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.

The collapse of the protostellar envelope results in the growth of the protostar and the development of a protoplanetary disk, playing a critical role during the early stages of star formation. Characterizing the gas infall in the envelope constrains the dynamical models of star formation. We present unambiguous signatures of infall, probed by optically thick molecular lines, toward an isolated embedded protostar, BHR 71 IRS1. The three dimensional radiative transfer calculations indicate that a slowly rotating infalling envelope model following the "inside-out" collapse reproduces the observations of both HCO$^{+}$ $J=4\rightarrow3$ and CS $J=7\rightarrow6$ lines, and the low velocity emission of the HCN $J=4\rightarrow3$ line. The envelope has a model-derived age of 12000$\pm$3000 years after the initial collapse. The envelope model underestimates the high velocity emission at the HCN $J=4\rightarrow3$ and H$^{13}$CN $J=4\rightarrow3$ lines, where outflows or a Keplerian disk may contribute. The ALMA observations serendipitously discover the emission of complex organic molecules (COMs) concentrated within a radius of 100 au, indicating that BHR 71 IRS1 harbors a hot corino. Eight species of COMs are identified, including CH$_{3}$OH and CH$_{3}$OCHO, along with H$_{2}$CS, SO$_{2}$ and HCN $v_{2}=1$. The emission of methyl formate and $^{13}$C-methanol shows a clear velocity gradient within a radius of 50 au, hinting at an unresolved Keplerian rotating disk.

We present ALMA observations of organic molecules towards five low-mass Class 0/I protostellar disk candidates in the Serpens cluster. Three sources (Ser-emb 1, Ser-emb 8, and Ser-emb 17) present emission of CH3OH as well as CH3OCH3, CH3OCHO, and CH2CO, while NH2CHO is detected in just Ser-emb 8 and Ser-emb 17. Detecting hot corino-type chemistry in three of five sources represents a high occurrence rate given the relative sparsity of these sources in the literature, and this suggests a possible link between protostellar disk formation and hot corino formation. For sources with CH3OH detections, we derive column densities of 10^{17}-10^{18} cm^{-2} and rotational temperatures of ~200-250 K. The CH3OH-normalized column density ratios of large, oxygen-bearing COMs in the Serpens sources and other hot corinos span two orders of magnitude, demonstrating a high degree of chemical diversity at the hot corino stage. Resolved observations of a larger sample of objects are needed to understand the origins of chemical diversity in hot corinos, and the relationship between different protostellar structural elements on disk-forming scales.

We have observed 29 transitions corresponding to 12 distinct species and 7 additional isotopologs toward the deeply embedded Class 0 young stellar object Ser-emb 1 with the Atacama Large Millimeter Array at 1 mm. The detected species include CH3OH and two complex organic molecules, CH3OCH3 and CH3OCHO. The emission of CH3OH and the two COMs is compact, and the CH3OH rotational temperature is 261 (46) K, implying that Ser-emb 1 hosts a hot corino. The derived CH3OH, CH3OCH3 and CH3OCHO column densities are at least 1.2 (0.4)E17 cm-2, 9.2 (3.8)E16 cm-2, and 9.1 (3.6)E16 cm-2, respectively, comparable to the values found for other Class 0 hot corinos. In addition, we observe evidence of rotation at compact scales: two of the more strongly detected lines, corresponding to C18O and H2CO present spatially resolved red- and blue-shifted compact emission orthogonal to the direction of a jet and outflow traced by CO, SiO, and several other molecules. The spatial coincidence of the hot corino emission and a possible disk in a compact region around the central protostar suggests that these structures may be physically and/or chemically related.

Modern astronomy has finally been able to observe protoplanetary disks in reasonable resolution and detail, unveiling the processes happening during planet formation. These observed processes are understood under the framework of disk-planet interaction, a process studied analytically and modeled numerically for over 40 years. Long a theoreticians' game, the wealth of observational data has been allowing for increasingly stringent tests of the theoretical models. Modeling efforts are crucial to support the interpretation of direct imaging analyses, not just for potential detections but also to put meaningful upper limits on mass accretion rates and other physical quantities in current and future large-scale surveys. This white paper addresses the questions of what efforts on the computational side are required in the next decade to advance our theoretical understanding, explain the observational data, and guide new observations. We identified the nature of accretion, ab initio planet formation, early evolution, and circumplanetary disks as major fields of interest in computational planet formation. We recommend that modelers relax the approximations of alpha-viscosity and isothermal equations of state, on the grounds that these models use flawed assumptions, even if they give good visual qualitative agreement with observations. We similarly recommend that population synthesis move away from 1D hydrodynamics. The computational resources to reach these goals should be developed during the next decade, through improvements in algorithms and the hardware for hybrid CPU/GPU clusters. Coupled with high angular resolution and great line sensitivity in ground based interferometers, ELTs and JWST, these advances in computational efforts should allow for large strides in the field in the next decade.

We present full spectral scans from 200-670$\mu$m of 26 Class 0+I protostellar sources, obtained with $Herschel$-SPIRE, as part of the "COPS-SPIRE" Open Time program, complementary to the DIGIT and WISH Key programs. Based on our nearly continuous, line-free spectra from 200-670 $\mu$m, the calculated bolometric luminosities ($L_{\rm bol}$) increase by 50% on average, and the bolometric temperatures ($T_{\rm bol}$) decrease by 10% on average, in comparison with the measurements without Herschel. Fifteen protostars have the same Class using $T_{\rm bol}$ and $L_{\rm bol}$/$L_{\rm submm}$. We identify rotational transitions of CO lines from J=4-3 to J=13-12, along with emission lines of $^{13}$CO, HCO$^+$, H$_{2}$O, and [CI]. The ratios of $^{12}$CO to $^{13}$CO indicate that $^{12}$CO emission remains optically thick for $J_{\rm up}$ < 13. We fit up to four components of temperature from the rotational diagram with flexible break points to separate the components. The distribution of rotational temperatures shows a primary population around 100 K with a secondary population at $\sim$350 K. We quantify the correlations of each line pair found in our dataset, and find the strength of correlation of CO lines decreases as the difference between $J$-level between two CO lines increases. The multiple origins of CO emission previously revealed by velocity-resolved profiles are consistent with this smooth distribution if each physical component contributes to a wide range of CO lines with significant overlap in the CO ladder. We investigate the spatial extent of CO emission and find that the morphology is more centrally peaked and less bipolar at high-$J$ lines. We find the CO emission observed with SPIRE related to outflows, which consists two components, the entrained gas and shocked gas, as revealed by our rotational diagram analysis as well as the studies with velocity-resolved CO emission.

CASTAway is a mission concept to explore our Solar System's main asteroid belt. Asteroids and comets provide a window into the formation and evolution of our Solar System and the composition of these objects can be inferred from space-based remote sensing using spectroscopic techniques. Variations in composition across the asteroid populations provide a tracer for the dynamical evolution of the Solar System. The mission combines a long-range (point source) telescopic survey of over 10,000 objects, targeted close encounters with 10 to 20 asteroids and serendipitous searches to constrain the distribution of smaller (e.g. 10 m) size objects into a single concept. With a carefully targeted trajectory that loops through the asteroid belt, CASTAway would provide a comprehensive survey of the main belt at multiple scales. The scientific payload comprises a 50 cm diameter telescope that includes an integrated low-resolution (R = 30 to 100) spectrometer and visible context imager, a thermal (e.g. 6 to 16 microns) imager for use during the flybys, and modified star tracker cameras to detect small (approx. 10 m) asteroids. The CASTAway spacecraft and payload have high levels of technology readiness and are designed to fit within the programmatic and cost caps for a European Space Agency medium class mission, whilst delivering a significant increase in knowledge of our Solar System.

We investigate the alignment between outflow axes in nine of the youngest binary/multiple systems in the Perseus Molecular Cloud. These systems have typical member spacing larger than 1000 AU. For outflow identification, we use 12CO(2-1) and 12CO(3-2) data from a large survey with the Sub- millimeter Array: Mass Assembly of Stellar Systems and their Evolution with the SMA (MASSES). The distribution of outflow orientations in the binary pairs is consistent with random or preferentially anti-aligned distributions, demonstrating that these outflows are misaligned. This result suggests that these systems are possibly formed in environments where the distribution of angular momentum is complex and disordered, and these systems do not come from the same co-rotating structures or from an initial cloud with aligned vectors of angular momentum.

We present continuum and molecular line observations at 230 GHz and 345 GHz from the Sub-millimeter Array (SMA) toward three protostars in the Perseus L1448N region. The data are from the large project "Mass Assembly of Stellar Systems and their Evolution with the SMA" (MASSES). Three dust continuum sources, Source B, Source NW, and Source A, are detected at both frequencies. These sources have corresponding emission peaks in C18O (J=2-1), 13CO (J=2-1), and HCO+ (J=4-3), and have offsets with N2D+ (J=3-2) peaks. High angular resolution data from a complimentary continuum survey with the Karl G. Jansky Very Large Array show that Source B is associated with three 8 mm continuum objects, Source NW with two, and Source A remains single. These results suggest that multiplicity in L1448N exists at different spatial scales from a few thousand AU to < 100 AU. Velocity gradients in each source obtained from two-dimensional fits to the SMA C18O emission are found to be perpendicular to within 20 degrees of the outflow directions as revealed by 12CO (J=2-1). We have observed that Sources B and NW with multiplicity have higher densities than Source A without multiplicity. This suggests that thermal Jeans fragmentation can be relevant in the fragmentation process. However, we have not observed a difference in the ratio between rotational and gravitational energy between sources with and without multiplicity. We also have not observed a trend between non-thermal velocity dispersions and the level of fragmentation. Our study has provided the first direct and comprehensive comparison between multiplicity and core properties in low-mass protostars, although based on small number statistics.

The accretion histories of embedded protostars are an integral part of descriptions of their physical and chemical evolution. In particular, are the accretion rates smoothly declining from the earlier toward later stages or in fact characterized by variations such as intermittent bursts? We aim to characterize the impact of possible accretion variations in a sample of embedded protostars by measuring the size of the inner regions of their envelopes where CO is sublimated and relate those to their temperature profiles dictated by their current luminosities. Using observations from the Submillimeter Array we measure the extents of the emission from the C18O isotopologue toward 16 deeply embedded protostars. We compare these measurements to the predicted extent of the emission given the current luminosities of the sources through dust and line radiative transfer calculations. Eight out of sixteen sources show more extended C18O emission than predicted by the models. The modeling shows that the likely culprit for these signatures is sublimation due to increases in luminosities of the sources by about a factor five or more during the recent 10,000 years - the time it takes for CO to freeze-out again on dust grains. For four of those sources the increase would have had to have been a factor 10 or more. The compact emission seen toward the other half of the sample suggests that C18O only sublimates when the temperature exceeds 30 K - as one would expect if CO is mixed with H2O in the grain ice-mantles. The small-number statistics from this survey suggest that protostars undergo significant bursts about once every 20,000 years. This also illustrates the importance of taking the physical evolutionary histories into account for descriptions of the chemical structures of embedded protostars.

Aims: Accretion rates in low-mass protostars can be highly variable in time. Each accretion burst is accompanied by a temporary increase in luminosity, heating up the circumstellar envelope and altering the chemical composition of the gas and dust. This paper aims to study such chemical effects and discusses the feasibility of using molecular spectroscopy as a tracer of episodic accretion rates and timescales. Methods: We simulate a strong accretion burst in a diverse sample of 25 spherical envelope models by increasing the luminosity to 100 times the observed value. Using a comprehensive gas-grain network, we follow the chemical evolution during the burst and for up to 10^5 yr after the system returns to quiescence. The resulting abundance profiles are fed into a line radiative transfer code to simulate rotational spectra of C18O, HCO+, H13CO+, and N2H+ at a series of time steps. We compare these spectra to observations taken from the literature and to previously unpublished data of HCO+ and N2H+ 6-5 from the Herschel Space Observatory. Results: The bursts are strong enough to evaporate CO throughout the envelope, which in turn enhances the abundance of HCO+ and reduces that of N2H+. After the burst, it takes 10^3-10^4 yr for CO to refreeze and for HCO+ and N2H+ to return to normal. The chemical effects of the burst remain visible in the rotational spectra for as long as 10^5 yr after the burst has ended, highlighting the importance of considering luminosity variations when analyzing molecular line observations in protostars. The spherical models are currently not accurate enough to derive robust timescales from single-dish observations. As follow-up work, we suggest that the models be calibrated against spatially resolved observations in order to identify the best tracers to be used for statistically significant source samples.

Abridged: Recent simulations have explored different ways to form accretion disks around low-mass stars. We aim to present observables to differentiate a rotationally supported disk from an infalling rotating envelope toward deeply embedded young stellar objects and infer their masses and sizes. Two 3D magnetohydrodynamics (MHD) formation simulations and 2D semi-analytical model are studied. The dust temperature structure is determined through continuum radiative transfer RADMC3D modelling. A simple temperature dependent CO abundance structure is adopted and synthetic spectrally resolved submm rotational molecular lines up to $J_{\rm u} = 10$ are simulated. All models predict similar compact components in continuum if observed at the spatial resolutions of 0.5-1$"$ (70-140 AU) typical of the observations to date. A spatial resolution of $\sim$14 AU and high dynamic range ($> 1000$) are required to differentiate between RSD and pseudo-disk in the continuum. The peak-position velocity diagrams indicate that the pseudo-disk shows a flatter velocity profile with radius than an RSD. On larger-scales, the CO isotopolog single-dish line profiles are similar and are narrower than the observed line widths of low-$J$ lines, indicating significant turbulence in the large-scale envelopes. However a forming RSD can provide the observed line widths of high-$J$ lines. Thus, either RSDs are common or a higher level of turbulence ($b \sim 0.8 \ {\rm km \ s^{-1}}$ ) is required in the inner envelope compared with the outer part. Multiple spatially and spectrally resolved molecular line observations are needed. The continuum data give a better estimate on disk masses whereas the disk sizes can be estimated from the spatially resolved molecular lines observations. The general observable trends are similar between the 2D semi-analytical models and 3D MHD RSD simulations.

We are conducting a Jansky VLA Ka-band (8 mm and 1 cm) and C-band (4 cm and 6.4 cm) survey of all known protostars in the Perseus Molecular Cloud, providing resolution down to $\sim$0.06'' and $\sim$0.35" in Ka-band and C-band, respectively. Here we present first results from this survey that enable us to examine the source NGC 1333 IRAS2A in unprecedented detail and resolve it into a proto-binary system separated by 0.621"$\pm$0.006" ($\sim$143 AU) at 8 mm, 1 cm, and 4 cm. These 2 sources (IRAS2A VLA1 and VLA2) are likely driving the two orthogonal outflows known to originate from IRAS2A. The brighter source IRAS2A VLA1 is extended perpendicular to its outflow in the VLA data, with a deconvolved size of 0.055" ($\sim$13 AU), possibly tracing a protostellar disk. The recently reported candidate companions (IRAS2A MM2 and MM3) are not detected in either our VLA data, CARMA 1.3 mm data, or SMA 850 $\mu$m data. SMA CO ($J=3\rightarrow2$), CARMA CO ($J=2\rightarrow1$), and lower resolution CARMA CO ($J=1\rightarrow0$) observations are used to examine the outflow origins and the nature of the candidate companions to IRAS2A VLA1. The CO ($J=3\rightarrow2$) and ($J=2\rightarrow1$) data show that IRAS2A MM2 is coincident with a bright CO emission spot in the east-west outflow, and IRAS2A MM3 is within the north-south outflow. In contrast, IRAS2A VLA2 lies at the east-west outflow symmetry point. We propose that IRAS2A VLA2 is the driving source of the East-West outflow and a true companion to IRAS2A VLA1, whereas IRAS2A MM2 and MM3 may not be protostellar.

We have surveyed a sample of massive star-forming regions located over a range of distances from the Galactic centre for methyl formate, HCOOCH$_{3}$, and its isotopologues H$^{13}$COOCH$_{3}$ and HCOO$^{13}$CH$_{3}$. The observations were carried out with the APEX telescope in the frequency range 283.4--287.4~GHz. Based on the APEX observations, we report tentative detections of the $^{13}$C-methyl formate isotopologue HCOO$^{13}$CH$_{3}$ towards the following four massive star-forming regions: Sgr~B2(N-LMH), NGC~6334~IRS~1, W51 e2 and G19.61-0.23. In addition, we have used the 1~mm ALMA science verification observations of Orion-KL and confirm the detection of the $^{13}$C-methyl formate species in Orion-KL and image its spatial distribution. Our analysis shows that the $^{12}$C/$^{13}$C isotope ratio in methyl formate toward Orion-KL Compact Ridge and Hot Core-SW components (68.4$\pm$10.1 and 71.4$\pm$7.8, respectively) are, for both the $^{13}$C-methyl formate isotopologues, commensurate with the average $^{12}$C/$^{13}$C ratio of CO derived toward Orion-KL. Likewise, regarding the other sources, our results are consistent with the $^{12}$C/$^{13}$C in CO. We also report the spectroscopic characterization, which includes a complete partition function, of the complex H$^{13}$COOCH$_{3}$ and HCOO$^{13}$CH$_{3}$ species. New spectroscopic data for both isotopomers H$^{13}$COOCH$_{3}$ and HCOO$^{13}$CH$_{3}$, presented in this study, has made it possible to measure this fundamentally important isotope ratio in a large organic molecule for the first time.

Protostars interact with their surroundings through jets and winds impacting on the envelope and creating shocks, but the nature of these shocks is still poorly understood. Our aim is to survey far-infrared molecular line emission from a uniform and significant sample of deeply-embedded low-mass young stellar objects in order to characterize shocks and the possible role of ultraviolet radiation in the immediate protostellar environment. Herschel/PACS spectral maps of 22 objects in the Perseus molecular cloud were obtained as part of the `William Herschel Line Legacy' survey. Line emission from H$_\mathrm{2}$O, CO, and OH is tested against shock models from the literature. Observed line ratios are remarkably similar and do not show variations with source physical parameters. Observations show good agreement with the shock models when line ratios of the same species are compared. Ratios of various H$_\mathrm{2}$O lines provide a particularly good diagnostic of pre-shock gas densities, $n_\mathrm{H}\sim10^{5}$ cm$^{-3}$, in agreement with typical densities obtained from observations of the post-shock gas. The corresponding shock velocities, obtained from comparison with CO line ratios, are above 20 km\,s$^{-1}$. However, the observations consistently show one-to-two orders of magnitude lower H$_\mathrm{2}$O-to-CO and H$_\mathrm{2}$O-to-OH line ratios than predicted by the existing shock models. The overestimated model H$_\mathrm{2}$O fluxes are most likely caused by an overabundance of H$_\mathrm{2}$O in the models since the excitation is well-reproduced. Illumination of the shocked material by ultraviolet photons produced either in the star-disk system or, more locally, in the shock, would decrease the H$_\mathrm{2}$O abundances and reconcile the models with observations. Detections of hot H$_\mathrm{2}$O and strong OH lines support this scenario.

(Abridged) The terrestrial planets, comets, and meteorites are significantly enriched in 15N compared to the Sun and Jupiter. While the solar and jovian nitrogen isotope ratio is believed to represent the composition of the protosolar nebula, a still unidentified process has caused 15N-enrichment in the solids. Several mechanisms have been proposed to explain the variations, including chemical fractionation. However, observational results that constrain the fractionation models are scarce. While there is evidence of 15N-enrichment in prestellar cores, it is unclear how the signature evolves into the protostellar phases. Our aim is to measure the 14N/15N ratio around three nearby, embedded low-to-intermediate-mass protostars. Isotopologues of HCN and HNC were used to probe the 14N/15N ratio. A selection of H13CN, HC15N, HN13C, and H15NC transitions was observed with the APEX telescope. The 14N/15N ratios were derived from the integrated intensities assuming a standard 12C/13C ratio. The assumption of optically thin emission was verified using radiative transfer modeling and hyperfine structure fitting. Two sources, IRAS 16293A and R CrA IRS7B, show 15N-enrichment by a factor of around 1.5-2.5 in both HCN and HNC with respect to the solar composition. Solar composition cannot be excluded for the third source, OMC-3 MMS6. Furthermore, there are indications of a trend toward increasing 14N/15N ratios with increasing outer envelope temperature. The enhanced 15N abundances in HCN and HNC found in two Class~0 sources (14N/15N of 160-290) and the tentative trend toward a temperature-dependent 14N/15N ratio are consistent with the chemical fractionation scenario, but 14N/15N ratios from additional tracers are indispensable for testing the models. Spatially resolved observations are needed to distinguish between chemical fractionation and isotope-selective photochemistry.

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

Low-mass protostars have been suggested to show highly variable accretion rates through-out their evolution. Such changes in accretion, and related heating of their ambient envelopes, may trigger significant chemical variations on different spatial scales and from source-to-source. We present images of emission from C17O, H13CO+, CH3OH, C34S and C2H toward the low-mass protostar IRAS 15398-3359 on 0.5" (75 AU diameter) scales with the Atacama Large Millimeter/submillimeter Array (ALMA) at 340 GHz. The resolved images show that the emission from H13CO+ is only present in a ring-like structure with a radius of about 1-1.5" (150-200 AU) whereas the CO and other high dipole moment molecules are centrally condensed toward the location of the central protostar. We propose that HCO+ is destroyed by water vapor present on small scales. The origin of this water vapor is likely an accretion burst during the last 100-1000 years increasing the luminosity of IRAS 15398-3359 by a factor of 100 above its current luminosity. Such a burst in luminosity can also explain the centrally condensed CH3OH and extended warm carbon-chain chemistry observed in this source and furthermore be reflected in the relative faintness of its compact continuum emission compared to other protostars.

PRISMA is a demonstration mission for formation-flying and on-orbit-servicing critical technologies that involves two spacecraft launched in low Earth orbit in June 2010 and still in operation. Funded by the Swedish National Space Board, PRISMA mission has been developed by OHB Sweden with important contributions from the German Aerospace Centre (DLR/GSOC), the French Space Agency (CNES), and the Technical University of Denmark (DTU). The paper focuses on the last CNES experiment achieved in September 2012 that was devoted to the preparation of future astrometry missions illustrated by the NEAT and microNEAT mission concepts. The experiment consisted in performing the type of formation maneuvers required to point the two-satellite axis to a celestial target and maintain it fixed during the observation period. Achieving inertial pointing for a LEO formation represented a new challenge given the numerous constraints from propellant usage to star tracker blinding. The paper presents the experiment objectives in relation with the NEAT/microNEAT mission concept, describes its main design features along with the guidance and control algorithms evolutions and discusses the results in terms of performances achieved during the two rehearsals

We present high angular resolution 1.3 mm and 850 um dust continuum data obtained with the Submillimeter Array toward 33 Class 0 protostars in nearby clouds (distance < 500 pc), which represents so far the largest survey toward protostellar binary/multiple systems. The median angular resolution in the survey is 2.5 arcsec, while the median linear resolution is approximately 600 AU. Compact dust continuum emission is observed from all sources in the sample. Twenty-one sources in the sample show signatures of binarity/multiplicity, with separations ranging from 50 to 5000 AU. The numbers of singles, binaries, triples, and quadruples in the sample are 12, 14, 5, and 2, respectively. The derived multiplicity frequency (MF) and companion star fraction (CSF) for Class 0 protostars are 0.64+/-0.08 and 0.91+/-0.05, respectively, with no correction for completeness. The derived MF and CSF in this survey are approximately two times higher than the values found in the binary surveys toward Class I YSOs, and approximately three (for MF) and four (for CSF) times larger than the values found among MS stars, with a similar range of separations. Furthermore, the observed fraction of high order multiple systems to binary systems in Class 0 protostars (0.50+/-0.09) is also larger than the fractions found in Class I YSOs (0.31+/-0.07) and MS stars (< 0.2). These results suggest that binary properties evolve as protostars evolve, as predicted by numerical simulations. The distribution of separations for Class 0 protostellar binary systems shows a general trend in which companion star fraction increases with decreasing companion separation. We find that 67%+/-8% of the protobinary systems have circumstellar mass ratios below 0.5, implying that unequal-mass systems are preferred in the process of binary star formation. We suggest an empirical sequential fragmentation picture for binary star formation.

Evaporation of water ice above 100 K in the inner few 100 AU of low-mass embedded protostars (the so-called hot core) should produce quiescent water vapor abundances of ~10^-4 relative to H2. Observational evidence so far points at abundances of only a few 10^-6. However, these values are based on spherical models, which are known from interferometric studies to be inaccurate on the relevant spatial scales. Are hot cores really that much drier than expected, or are the low abundances an artifact of the inaccurate physical models? We present deep velocity-resolved Herschel-HIFI spectra of the 3(12)-3(03) lines of H2-16O and H2-18O (1097 GHz, Eup/k = 249 K) in the low-mass Class 0 protostar NGC1333 IRAS2A. A spherical radiative transfer model with a power-law density profile is unable to reproduce both the HIFI data and existing interferometric data on the H2-18O 3(13)-2(20) line (203 GHz, Eup/k = 204 K). Instead, the HIFI spectra likely show optically thick emission from a hot core with a radius of about 100 AU. The mass of the hot core is estimated from the C18O J=9-8 and 10-9 lines. We derive a lower limit to the hot water abundance of 2x10^-5, consistent with the theoretical predictions of ~10^-4. The revised HDO/H2O abundance ratio is 1x10^-3, an order of magnitude lower than previously estimated.

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

Glycolaldehyde (HCOCH2OH) is the simplest sugar and an important intermediate in the path toward forming more complex biologically relevant molecules. In this paper we present the first detection of 13 transitions of glycolaldehyde around a solar-type young star, through Atacama Large Millimeter Array (ALMA) observations of the Class 0 protostellar binary IRAS 16293-2422 at 220 GHz (6 transitions) and 690 GHz (7 transitions). The glycolaldehyde lines have their origin in warm (200-300 K) gas close to the individual components of the binary. Glycolaldehyde co-exists with its isomer, methyl formate (HCOOCH3), which is a factor 10-15 more abundant toward the two sources. The data also show a tentative detection of ethylene glycol, the reduced alcohol of glycolaldehyde. In the 690 GHz data, the seven transitions predicted to have the highest optical depths based on modeling of the 220 GHz lines all show red-shifted absorption profiles toward one of the components in the binary (IRAS16293B) indicative of infall and emission at the systemic velocity offset from this by about 0.2" (25 AU). We discuss the constraints on the chemical formation of glycolaldehyde and other organic species - in particular, in the context of laboratory experiments of photochemistry of methanol-containing ices. The relative abundances appear to be consistent with UV photochemistry of a CH3OH-CO mixed ice that has undergone mild heating. The order of magnitude increase in line density in these early ALMA data illustrate its huge potential to reveal the full chemical complexity associated with the formation of solar system analogs.

We present Spitzer IRAC (2.1 sq. deg.) and MIPS (6.5 sq. deg.) observations of star formation in the Ophiuchus North molecular clouds. This fragmentary cloud complex lies on the edge of the Sco-Cen OB association, several degrees to the north of the well-known rho Oph star-forming region, at an approximate distance of 130 pc. The Ophiuchus North clouds were mapped as part of the Spitzer Gould Belt project under the working name `Scorpius'. In the regions mapped, selected to encompass all the cloud with visual extinction AV>3, eleven Young Stellar Object (YSO) candidates are identified, eight from IRAC/MIPS colour-based selection and three from 2MASS K/MIPS colours. Adding to one source previously identified in L43 (Chen et al. 2009), this increases the number of YSOcs identified in Oph N to twelve. During the selection process, four colour-based YSO candidates were rejected as probable AGB stars and one as a known galaxy. The sources span the full range of YSOc classifications from Class 0/1 to Class III, and starless cores are also present. Twelve high-extinction (AV>10) cores are identified with a total mass of approx. 100 solar masses. These results confirm that there is little ongoing star formation in this region (instantaneous star formation efficiency <0.34%) and that the bottleneck lies in the formation of dense cores. The influence of the nearby Upper Sco OB association, including the 09V star zeta Oph, is seen in dynamical interactions and raised dust temperatures but has not enhanced levels of star formation in Ophiuchus North.

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

It remains a key challenge to establish the molecular content of different components of low-mass protostars, like their envelopes and disks, and how this depends on the evolutionary stage and/or environment of the young stars. Observations at submillimeter wavelengths provide a direct possibility to study the chemical composition of low-mass protostars through transitions probing temperatures up to a few hundred K in the gas surrounding these sources. This paper presents a large molecular line survey of the deeply embedded protostellar binary IRAS 16293-2422 from the Submillimeter Array (SMA) - including images of individual lines down to approximately 1.5-3" (190-380 AU) resolution. More than 500 individual transitions are identified related to 54 molecular species (including isotopologues) probing temperatures up to about 550 K. Strong chemical differences are found between the two components in the protostellar system with a separation between, in particular, the sulfur- and nitrogen-bearing species and oxygen-bearing complex organics. The action of protostellar outflow on the ambient envelope material is seen in images of CO and SiO and appear to influence a number of other species, including (deuterated) water, HDO. The effects of cold gas-phase chemistry is directly imaged through maps of CO, N2D+ and DCO+, showing enhancements of first DCO+ and subsequently N2D+ in the outer envelope where CO freezes-out on dust grains.

We present Spitzer Space Telescope IRAC and MIPS observations of a 0.85 deg^2 field including the Corona Australis (CrA) star-forming region. At a distance of 130 pc, CrA is one of the closest regions known to be actively forming stars, particularly within its embedded association, the Coronet. Using the Spitzer data, we identify 51 young stellar objects (YSOs) in CrA which include sources in the well-studied Coronet cluster as well as distributed throughout the molecular cloud. Twelve of the YSOs discussed are new candidates, one of which is located in the Coronet. Known YSOs retrieved from the literature are also added to the list, and a total of 116 candidate YSOs in CrA are compiled. Based on these YSO candidates, the star formation rate is computed to be 12 M_o Myr^-1, similar to that of the Lupus clouds. A clustering analysis was also performed, finding that the main cluster core, consisting of 68 members, is elongated (having an aspect ratio of 2.36), with a circular radius of 0.59 pc and mean surface density of 150 pc^-2. In addition, we analyze outflows and jets in CrA by means of new CO and H_2 data. We present 1.3 mm interferometric continuum observations made with the Submillimeter Array (SMA) covering R CrA, IRS 5, IRS 7, and IRAS 18595-3712 (IRAS 32). We also present multi-epoch H_2 maps and detect jets and outflows, study their proper motions, and identify exciting sources. The Spitzer and ISAAC/VLT observations of IRAS 32 show a bipolar precessing jet, which drives a CO (2-1) outflow detected in the SMA observations. There is also clear evidence for a parsec-scale precessing outflow, E-W oriented, and originating in the SMA 2 region, likely driven by SMA 2 or IRS 7A.

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 HDO/H2O abundance ratio is thought to be a key diagnostic on the evolution of water during the star- and planet-formation process and thus on its origin on Earth. We here present millimeter-wavelength high angular resolution observations of the deeply embedded protostar NGC1333-IRAS4B from the Submillimeter Array targeting the 3(1,2)-2(2,1) transition of HDO at 225.6 GHz (Eu = 170 K). We do not (or only very tentatively) detect the HDO line toward the central protostar, contrasting the previous prominent detection of a line from another water isotopologue, H2-18-O, with similar excitation properties using the IRAM Plateau de Bure Interferometer. The non-detection of the HDO line provides a direct, model independent, upper limit to the HDO/H2O abundance ratio of 6e(-4) (3sigma) in the warm gas associated with the central protostar. This upper limit suggests that the HDO/H2O abundance ratio is not significantly enhanced in the inner 50 AU around the protostar relative to what is seen in comets and Earth's oceans and does not support previous suggestions of a generally enhanced HDO/H2O ratio in these systems.

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.

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.

Gas-phase complex organic molecules have been detected toward a range of high- and low-mass star-forming regions at abundances which cannot be explained by any known gas-phase chemistry. Recent laboratory experiments show that UV irradiation of CH3OH-rich ices may be an important mechanism for producing complex molecules and releasing them into the gas-phase. To test this ice formation scenario we mapped the B1-b dust core and nearby protostar in CH3OH gas using the IRAM 30m telescope to identify locations of efficient non-thermal ice desorption. We find three CH3OH abundance peaks tracing two outflows and a quiescent region on the side of the core facing the protostar. The CH3OH gas has a rotational temperature of ~10 K at all locations. The quiescent CH3OH abundance peak and one outflow position were searched for complex molecules. Narrow, 0.6-0.8 km s-1 wide, HCOOCH3 and CH3CHO lines originating in cold gas are clearly detected, CH3OCH3 is tentatively detected and C2H5OH and HOCH2CHO are undetected toward the quiescent core, while no complex molecular lines were found toward the outflow. The core abundances with respect to CH3OH are ~2.3% and 1.1% for HCOOCH3 and CH3CHO, respectively, and the upper limits are 0.7-1.1%, which is similar to most other low-mass sources. The observed complex molecule characteristics toward B1-b and the pre-dominance of HCO-bearing species suggest a cold ice (below 25 K, the sublimation temperature of CO) formation pathway followed by non-thermal desorption through e.g. UV photons traveling through outflow cavities. The observed complex gas composition together with the lack of any evidence of warm gas-phase chemistry provide clear evidence of efficient complex molecule formation in cold interstellar ices.

We present Spitzer IRS mid-infrared (5-35 micron) spectra of a complete flux-limited sample (> 3 mJy at 8 micron) of young stellar object (YSO) candidates selected on the basis of their infrared colors in the Serpens Molecular Cloud. Spectra of 147 sources are presented and classified. Background stars (with slope consistent with a reddened stellar spectrum and silicate features in absorption), galaxies (with redshifted PAH features) and a planetary nebula (with high ionization lines) amount to 22% of contamination in this sample, leaving 115 true YSOs. Sources with rising spectra and ice absorption features, classified as embedded Stage I protostars, amount to 18% of the sample. The remaining 82% (94) of the disk sources are analyzed in terms of spectral energy distribution shapes, PAHs and silicate features. The presence, strength and shape of these silicate features are used to infer disk properties for these systems. About 8% of the disks have 30/13 micron flux ratios consistent with cold disks with inner holes or gaps, and 3% of the disks show PAH emission. Comparison with models indicates that dust grains in the surface of these disks have sizes of at least a few \mu\m. The 20 micron silicate feature is sometimes seen in absence of the 10 micron feature, which may be indicative of very small holes in these disks. No significant difference is found in the distribution of silicate feature shapes and strengths between sources in clusters and in the field. Moreover, the results in Serpens are compared with other well-studied samples: the c2d IRS sample distributed over 5 clouds and a large sample of disks in the Taurus star-forming region. The remarkably similar distributions of silicate feature characteristics in samples with different environment and median ages - if significant - imply that the dust population in the disk surface results from an equilibrium between dust growth and destructive collision processes that are maintained over a few million years for any YSO population irrespective of environment.

Water is one of the key molecules in the physical and chemical evolution of star- and planet-forming regions. We here report the first spatially resolved observation of thermal emission of (an isotopologue of) water with the Plateau de Bure Interferometer toward the deeply embedded Class 0 protostar NGC 1333-IRAS4B. The observations of the H2-18-O 3_13-2_20 transition at 203.4 GHz resolve the emission of water toward this source with an extent of about 0.2" corresponding to the inner 25 AU (radius). The H2-18-O emission reveals a tentative velocity gradient perpendicular to the extent of the protostellar outflow/jet probed by observations of CO rotational transitions and water masers. The line is narrow, about 1 km/s (FWHM), significantly less than what would be expected for emission from an infalling envelope or accretion shock, but consistent with emission from a disk seen at a low inclination angle. The water column density inferred from these data suggests that the water emitting gas is a thin warm layer containing about 25 M_Earth of material, 0.03% of the total disk mass traced by continuum observations.

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.

Observations of Lynds Dark Nebula 1221 from the Spitzer Space Telescope are presented. These data show three candidate protostars towards L1221, only two of which were previously known. The infrared observations also show signatures of outflowing material, an interpretation which is also supported by radio observations with the Very Large Array. In addition, molecular line maps from the Five College Radio Astronomy Observatory are shown. One-dimensional dust continuum modelling of two of these protostars, IRS1 and IRS3, is described. These models show two distinctly different protostars forming in very similar environments. IRS1 shows a higher luminosity and larger inner radius of the envelope than IRS3. The disparity could be caused by a difference in age or mass, orientation of outflow cavities, or the impact of a binary in the IRS1 core.

Aims: To investigate the chemical relations between complex organics based on their spatial distributions and excitation conditions in the low-mass young stellar objects IRAS 16293-2422 A and B. Methods: Interferometric observations with the Submillimeter Array have been performed at 5''x3'' resolution revealing emission lines of HNCO, CH3CN, CH2CO, CH3CHO and C2H5OH. Rotational temperatures are determined from rotational diagrams when a sufficient number of lines are detected. Results: Compact emission is detected for all species studied here. For HNCO and CH3CN it mostly arises from source A, CH2CO and C2H5OH have comparable strength for both sources and CH3CHO arises exclusively from source B. HNCO, CH3CN and CH3CHO have rotational temperatures >200 K. The (u,v)-visibility data reveal that HNCO also has extended cold emission. Conclusions: The abundances of the molecules studied here are very similar within factors of a few to those found in high-mass YSOs. Thus the chemistry between high- and low-mass objects appears to be independent of luminosity and cloud mass. Bigger abundance differences are seen between the A and B source. The HNCO abundance relative to CH3OH is ~4 times higher toward A, which may be due to a higher initial OCN- ice abundances in source A compared to B. Furthermore, not all oxygen-bearing species are co-existent. The different spatial behavior of CH2CO and C2H5OH compared with CH3CHO suggests that hydrogenation reactions on grain-surfaces are not sufficient to explain the observed gas phase abundances. Selective destruction of CH3CHO may result in the anti-coincidence of these species in source A. These results illustrate the power of interferometric compared with single dish data in terms of testing chemical models.

We present a census of the population of deeply embedded young stellar objects (YSOs) in the Ophiuchus molecular cloud complex based on a combination of Spitzer Space Telescope mid-infrared data from the "Cores to Disks" (c2d) legacy team and JCMT/SCUBA submillimeter maps from the COMPLETE team. We have applied a method developed for identifying embedded protostars in Perseus to these datasets and in this way construct a relatively unbiased sample of 27 candidate embedded protostars with envelopes more massive than our sensitivity limit (about 0.1 M_sun). Embedded YSOs are found in 35% of the SCUBA cores - less than in Perseus (58%). On the other hand the mid-infrared sources in Ophiuchus have less red mid-infrared colors, possibly indicating that they are less embedded. We apply a nearest neighbor surface density algorithm to define the substructure in each of the clouds and calculate characteristic numbers for each subregion - including masses, star formation efficiencies, fraction of embedded sources etc. Generally the main clusters in Ophiuchus and Perseus (L1688, NGC1333 and IC348) are found to have higher star formation efficiencies than small groups such as B1, L1455 and L1448, which on the other hand are completely dominated by deeply embedded protostars. We discuss possible explanations for the differences between the regions in Perseus and Ophiuchus, such as different evolutionary timescales for the YSOs or differences, e.g., in the accretion in the two clouds.

We present c2d Spitzer/IRAC observations of the Lupus I, III and IV dark clouds and discuss them in combination with optical and near-infrared and c2d MIPS data. With the Spitzer data, the new sample contains 159 stars, 4 times larger than the previous one. It is dominated by low- and very-low mass stars and it is complete down to M $\approx$ 0.1M$_\odot$. We find 30-40 % binaries with separations between 100 to 2000 AU with no apparent effect in the disk properties of the members. A large majority of the objects are Class II or Class III objects, with only 20 (12%) of Class I or Flat spectrum sources. The disk sample is complete down to ``debris''-like systems in stars as small as M $\approx$ 0.2 M$_\odot$ and includes sub-stellar objects with larger IR excesses. The disk fraction in Lupus is 70 -- 80%, consistent with an age of 1 -- 2 Myr. However, the young population contains 20% optically thick accretion disks and 40% relatively less flared disks. A growing variety of inner disk structures is found for larger inner disk clearings for equal disk masses. Lupus III is the most centrally populated and rich, followed by Lupus I with a filamentary structure and by Lupus IV, where a very high density core with little star-formation activity has been found. We estimate star formation rates in Lupus of 2 -- 10 M$_\odot$ Myr$^{-1}$ and star formation efficiencies of a few percent, apparently correlated with the associated cloud masses.

(abbreviated) We aim to determine the masses of the envelopes, disks, and central stars of young stellar objects (YSOs) in the Class I stage. We observed the embedded Class I objects IRS 63 and Elias 29 in the rho Ophiuchi star-forming region with the Submillimeter Array (SMA) at 1.1 mm. IRS 63 and Elias 29 are both clearly detected in the continuum, with peak fluxes of 459 resp. 47 mJy/beam. The continuum emission toward Elias 29 is clearly resolved, whereas IRS 63 is consistent with a point source down to a scale of 3 arcsec (400 AU). The SMA data are combined with single-dish data, and disk masses of 0.055 and >= 0.007 MSun and envelope masses of 0.058 and >= 0.058 MSun are empirically determined for IRS 63 and Elias 29, respectively. The disk+envelope systems are modelled with the axisymmetric radiative-transfer code RADMC, yielding disk and envelope masses that differ from the empirical results by factors of a few. HCO+ J = 3-2 is detected toward both sources, HCN J = 3-2 is not. The HCO+ position-velocity diagrams are indicative of Keplerian rotation. For a fiducial inclination of 30 degrees, we find stellar masses of 0.37 +/- 0.13 and 2.5 +/- 0.6 MSun for IRS 63 and Elias 29, respectively. We conclude that the sensitivity and spatial resolution of the SMA at 1.1 mm allow a good separation of the disks around Class I YSOs from their circumstellar envelopes and environments, and the spectral resolution makes it possible to resolve their dynamical structure and estimate the masses of the central stars. The ratios of the envelope and disk masses are found to be 0.2 and 6 for IRS 63 and Elias 29, respectively. This is lower than the values for Class 0 sources, which have Menv/Mdisk >= 10, suggesting that this ratio is a tracer of the evolutionary stage of a YSO.

With the goal to study the physical and chemical evolution of ices in solar-mass systems, a spectral survey is conducted of a sample of 41 low luminosity YSOs using 3-38 um Spitzer and ground-based spectra. The long-known 6.0 and 6.85 um bands are detected toward all sources, with the Class 0-type YSOs showing the deepest bands ever observed. In almost all sources the 6.0 um band is deeper than expected from the bending mode of pure solid H2O. The depth and shape variations of the remaining 5-7 um absorption indicate that it consists of 5 independent components, which, by comparison to laboratory studies, must be from at least 8 different carriers. Simple species are responsible for much of the absorption in the 5-7 um region, at abundances of 1-30% for CH3OH, 3-8% for NH3, 1-5% for HCOOH, ~6% for H2CO, and ~0.3% for HCOO- with respect to solid H2O. The 6.85 um band likely consists of one or two carriers, of which one is less volatile than H2O because its abundance relative to H2O is enhanced at lower H2O/tau_9.7 ratios. It does not survive in the diffuse interstellar medium (ISM), however. The similarity of the 6.85 um bands for YSOs and background stars indicates that its carrier(s) must be formed early in the molecular cloud evolution. If an NH4+ salt is the carrier its abundance with respect to solid H2O is typically 7%, and low temperature acid-base chemistry or cosmic ray induced reactions must have been involved in its formation. Possible origins are discussed for the carrier of an enigmatic, very broad absorption between 5 and 8 um. Finally, all the phenomena observed for ices toward massive YSOs are also observed toward low mass YSOs, indicating that processing of the ices by internal ultraviolet radiation fields is a minor factor in the early chemical evolution of the ices. [abridged]

We present maps of 14.4 deg^2 of the Ophiuchus dark clouds observed by the Spitzer Space Telescope Multiband Imaging Photometer for Spitzer (MIPS). These high quality maps depict both numerous point sources as well as extended dust emission within the star-forming and non-star-forming portions of these clouds. Using PSF-fitting photometry, we detect 5779 sources at 24 um and 81 sources at 70 um at the 10 sigma level of significance. Three hundred twenty-three candidate young stellar objects (YSOs) were identified according to their positions on the MIPS/2MASS K versus K$-$[24] color-magnitude diagrams as compared to 24 um detections in the SWIRE extragalactic survey. We find that more than half of the YSO candidates, and almost all the ones with protostellar Class I spectral energy distributions, are confined to the known cluster and aggregates.

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.

Submillimeter observations with ALMA will be the essential next step in our understanding of how stars and planets form. Key projects range from detailed imaging of the collapse of pre-stellar cores and measuring the accretion rate of matter onto deeply embedded protostars, to unravelling the chemistry and dynamics of high-mass star-forming clusters and high-spatial resolution studies of protoplanetary disks down to the 1 AU scale.

(Abridged) With the SMA we have made high angular-resolution (~1" = 160 AU) observations of the protobinary system IRAS 16293-2422 in the J = 4-3 lines of HCN and HC^15N, and in the continuum at 354.5 GHz. The HCN (4-3) line was also observed using the JCMT to supply missing short spacing information. Submillimeter continuum emission is detected from the individual binary components with a separation of ~5". The HC^15N (4-3) emission has revealed a compact (~500 AU) flattened structure (P.A. = -16 degree) associated with Source A, and shows a velocity gradient along the projected minor axis, which can be interpreted as an infalling gas motion. Our HCN image including the short-spacing information shows an extended (~3000 AU) circumbinary envelope as well as the compact structure associated with Source A. A toy model consisting of a flattened structure with radial infall towards a 1 Msun central star reproduces the HCN/HC^15N position-velocity diagram along the minor axis of the HC^15N emission. In the extended envelope there is also a North-East (Blue) to South-West (Red) velocity gradient across the binary alignment, which is likely to reflect gas motion in the swept-up dense gas associated with the molecular outflow from Source A. At Source B, there is only a weak compact structure with much narrower line widths (~2 km/s) seen in the optically-thin HC^15N emission than that at Source A (>10 km/s), and there is no clearly defined bipolar molecular outflow associated with Source B. These results imply the different evolutionary stages between Source A and B in the common circumbinary envelope.

A partial submillimeter line-survey was performed toward 7 high-mass young stellar objects (YSOs) aimed at detecting H2CO, CH3OH, CH2CO, CH3CHO, C2H5OH, HCOOH, HNCO and NH2CHO. In addition, other organic species such as CH3CN have been observed. The aim is to establish the chemical origin of a set of complex organic molecules thought to be produced by grain surface chemistry in high mass YSOs. Rotation temperatures and beam-averaged column densities are determined. Based on their rotation diagrams, molecules can be classified as either cold (<100 K) or hot (>100 K). This implies that complex organics are present in at least two distinct regions. Furthermore, the abundances of the hot oxygen-bearing species are correlated, as are those of HNCO and NH2CHO. This is suggestive of chemical relationships within, but not between, those two groups of molecules. The most likely explanation for the observed correlations of the various hot molecules is that they are ``first generation'' species that originate from solid-state chemistry. This includes H2CO, CH3OH, C2H5OH, HCOOCH3, CH3OCH3, HNCO, NH2CHO, and possibly CH3CN, and C2H5CN. The correlations between sources implies very similar conditions during their formation or very similar doses of energetic processing. Cold species such as CH2CO, CH3CHO, and HCOOH, some of which are seen as ices along the same lines of sight, are probably formed in the solid state as well, but appear to be destroyed at higher temperatures. A low level of non-thermal desorption by cosmic rays can explain their low rotation temperatures and relatively low abundances in the gas phase compared to the solid state. The CH3CCH abundances can be fully explained by low temperature gas phase chemistry. No cold N-containing molecules are found.

Submillimeter lines of CN, NO, CO+ and SO+, and upper limits on SH+ and N2O are observed with the James Clerk Maxwell Telescope in two high-mass and up to nine low-mass young stellar objects and compared with chemical models. Constant fractional abundances derived from radiative transfer modeling of the line strengths are x(CN) ~ a few x 10^{-11}-10^{-8}, x(NO) ~ 10^{-9}-10^{-8} and x(CO+) ~ 10^{-12}-10^{-10}. SO+ has abundances of a few x 10^{-11} in the high-mass objects and upper limits of ~ 10^{-12}-10^{-11} in the low-mass sources. All abundances are up to 1-2 orders of magnitude higher if the molecular emission is assumed to originate mainly from the inner region (< 1000 AU) of the envelope. For high-mass sources, the CN, SO+ and CO+ abundances and abundance ratios are best explained by an enhanced far-ultraviolet (FUV) field impacting gas at temperatures of a few hundred K. The observed column densities require that this region of enhanced FUV has scales comparable to the observing beam, such as in a geometry in which the enhanced FUV irradiates outflow walls. For low-mass sources, the required temperatures within the FUV models of T > 300 K are much higher than found in models, so that an X-ray enhanced region close to the protostar (r < 500 AU) is more plausible. The observed abundances imply X-ray fluxes for the Class 0 objects of L_X ~ 10^{29}-10^{31} erg s^{-1}, comparable to those observed from low-mass Class I protostars. Spatially resolved data are needed to clearly distinguish the effects of FUV and X-rays for individual species.

This paper presents a large spectral line and continuum survey of 8 deeply embedded, low-mass protostellar cores using the Submillimeter Array. Each source was observed in high excitation lines of some of the most common molecular species, CO, HCO+, CS, SO, H2CO, CH3OH and SiO. Line emission from 11 species originating from warm and dense gas have been imaged at high angular resolution (1-3"; typically 200-600 AU) together with continuum emission at 230 GHz (1.3 mm) and 345 GHz (0.8 mm). Compact continuum emission is observed for all sources which likely originates in marginally optically thick circumstellar disks, with typical lower limits to their masses of 0.1 M_sun (1-10% of the masses of their envelopes) and having a dust opacity law with beta approximately 1. Prominent outflows are present in CO 2-1 observations in all sources: the most diffuse outflows are found in the sources with the lowest ratios of disk-to-envelope mass, and it is suggested that these sources are in a phase where accretion of matter from the envelope has almost finished and the remainder of the envelope material is being dispersed by the outflows. Other characteristic dynamical signatures are found with inverse P Cygni profiles indicative of infalling motions seen in the 13CO 2-1 lines toward NGC1333-IRAS4A and -IRAS4B. Outflow-induced shocks are present on all scales in the protostellar environments and are most clearly traced by the emission of CH3OH in NGC1333-IRAS4A and -IRAS4B. These observations suggest that the emission of CH3OH and H2CO from these proposed "hot corinos" are related to the shocks caused by the protostellar outflows. Only one source, NGC1333-IRAS2A, has evidence for hot, compact CH3OH emission coincident with the embedded protostar.

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.

CONTEXT - Low-mass stars form with disks in which the coagulation of grains may eventually lead to the formation of planets. It is not known when and where grain growth occurs, as models that explain the observations are often degenerate. A way to break this degeneracy is to resolve the sources under study. AIMS - To find evidence for the existence of grains of millimetre sizes in disks around in T Tauri stars, implying grain growth. METHODS - The Australia Telescope Compact Array (ATCA) was used to observe 15 southern T Tauri stars, five in the constellation Lupus and ten in Chamaeleon, at 3.3 millimetre. The five Lupus sources were also observed with the Submillimeter Array (SMA) at 1.4 millimetre. Our new data are complemented with data from the literature to determine the slopes of the spectral energy distributions in the millimetre regime. RESULTS - Ten sources were detected at better than 3sigma with the ATCA, with sigma ~1-2 mJy, and all sources that were observed with the SMA were detected at better than 15sigma, with sigma ~4 mJy. Six of the sources in our sample are resolved to physical radii of ~100 AU. Assuming that the emission from such large disks is predominantly optically thin, the millimetre slope can be related directly to the opacity index. For the other sources, the opacity indices are lower limits. Four out of six resolved sources have opacity indices <~1, indicating grain growth to millimetre sizes and larger. The masses of the disks range from < 0.01 to 0.08 MSun, which is comparable to the minimum mass solar nebula. A tentative correlation is found between the millimetre slope and the strength and shape of the 10-micron silicate feature, indicating that grain growth occurs on similar (short) timescales in both the inner and outer disk.

We present a census of the population of deeply embedded young stellar objects (YSOs) in the Perseus molecular cloud complex based on a combination of Spitzer Space Telescope mid-IR data from the c2d legacy team and JCMT/SCUBA submillimeter maps from the COMPLETE team. The mid-IR sources detected at 24 micron and having [3.6]-[4.5] > 1 are located close to the center of the SCUBA cores, typically within 15" of their peaks. The narrowness of the spatial distribution of mid-IR sources around the peaks of the SCUBA cores suggests that no significant dispersal of the newly formed YSOs has occurred. This argues against the suggestion that motions of protostars regulate the time scales over which significant (Bondi-Hoyle) accretion can occur. The most deeply embedded YSOs are found in regions with high extinction, AV > 5, similar to the extinction threshold observed for the SCUBA cores. All the SCUBA cores with high concentrations have embedded YSOs, but not all cores with low concentrations are starless. An unbiased sample of 49 deeply embedded YSOs is constructed. Embedded YSOs are found in 40 of the 72 SCUBA cores with only three cores harboring multiple embedded YSOs within 15". The equal number of SCUBA cores with and without embedded YSOs suggests that the time scale for the evolution through the dense prestellar stages, where the cores are recognized in the submillimeter maps and have central densities of 5e4-1e5 cm^{-3}, is similar to the time scale for the embedded protostellar stages. The current star formation efficiency of cores is estimated to be approximately 10-15%. In contrast, the star formation efficiency averaged over the cloud life time and compared to the total cloud mass is only a few percent, reflecting also the efficiency in assembling cloud material into the dense cores forming stars.

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.

This paper presents new APEX submillimetre molecular line observations of three low-mass protostars, IRS7A, IRS7B, and IRAS32, in the R Coronae Australis molecular cloud complex. The molecular excitation analysis is performed using a statistical equilibrium radiative transfer code. The derived beam averaged fractional abundances vary by less than a factor of two among the three sources, except those of H2CO and CH3OH, which show differences of about an order of magnitude. The molecular abundances are similar to those typically found in other star-forming regions in the Galaxy, such as the $\rho$~Oph and Perseus molecular clouds. There is a marked difference in the kinetic temperatures derived for the protobinary source IRS7 from H2CO (40-60 K) and CH3OH (20 K), possibly indicating a difference in origin of the emission from these two molecules.

We present observations of 3.86 sq. deg. of the Perseus molecular cloud complex with the Spitzer Space Telescope Infrared Array Camera (IRAC). The maps show strong extended emission arising from shocked H2 in outflows in the region and from polycyclic aromatic hydrocarbon features. More than 120,000 sources are extracted toward the cloud. Based on their IRAC colors and comparison to off-cloud and extragalactic fields, we identify 400 candidate young stellar objects. About two thirds of these are associated with the young clusters IC348 and NGC1333, while the last third is distributed over the remaining cloud. We classify the young stellar objects according to the traditional scheme based on the slope of their spectral energy distributions. Significant differences are found for the numbers of embedded Class I objects relative to more evolved Class II objects in IC348, NGC1333 and the remaining cloud with the embedded Class I and "flat spectrum" YSOs constituting 14%, 36% and 47% of the total number of YSOs identified in each of these regions. The high number of Class I objects in the extended cloud (61% of the Class I objects in the entire cloud) suggests that a significant fraction of the current star formation is occuring outside the two main clusters. Finally we discuss a number of outflows and identify their driving sources, including the known deeply embedded Class 0 sources outside the two major clusters. The Class 0 objects are found to be detected by Spitzer and have very red [3.6]-[4.5] colors but do not show similarly red [5.8]-[8.0] colors. The Class 0 objects are easily identifiable in color-color diagrams plotting these two colors but are problematic to extract automatically due to the extended emission from shocked gas or scattered light in cavities related to the associated outflows.

We study the H2O chemistry in star-forming environments under the influence of a central X-ray source and a central far ultraviolet (FUV) radiation field. The gas-phase water chemistry is modeled as a function of time, hydrogen density and X-ray flux. To cover a wide range of physical environments, densities between n_H = 10^4-10^9 cm^-3 and temperatures between T = 10-1000 K are studied. Three different regimes are found: For T < 100 K, the water abundance is of order 10^-7-10^-6 and can be somewhat enhanced or reduced due to X-rays, depending on time and density. For 100 K < T < 250 K, H2O is reduced from initial x(H2O) ~ 10^-4 following ice evaporation to x(H2O) ~ 10^-6 for F_X > 10^-3 ergs s-1 cm^-2 (t = 10^4 yrs) and for F_X > 10^-4 ergs s^-1 cm^-2 (t = 10^5 yrs). At higher temperatures (T > 250 K) and hydrogen densities, water can persist with x(H2O) ~ 10^-4 even for high X-ray fluxes. The X-ray and FUV models are applied to envelopes around low-mass Class 0 and I young stellar objects (YSOs). Water is destroyed in both Class 0 and I envelopes on relatively short timescales (t ~ 5000 yrs) for realistic X-ray fluxes, although the effect is less prominent in Class 0 envelopes due to the higher X-ray absorbing densities there. FUV photons from the central source are not effective in destroying water. The average water abundance in Class I sources for L_X > 10^27 ergs s^-1 is predicted to be x(H2O) < 10^-6.

We discuss the effects of an enhanced interstellar radiation field (ISRF) on the observables of protostellar cores in the Orion cloud region. Dust radiative transfer is used to constrain the envelope physical structure by reproducing SCUBA 850 micron emission. Previously reported 13CO, C17O and H2CO line observations are reproduced through detailed Monte Carlo line radiative transfer models. It is found that the 13CO line emission is marginally optically thick and sensitive to the physical conditions in the outer envelope. An increased temperature in this region is needed in order to reproduce the 13CO line strengths and it is suggested to be caused by a strong heating from the exterior, corresponding to an ISRF in Orion 10^3 times stronger than the "standard" ISRF. The typical temperatures in the outer envelope are higher than the desorption temperature for CO. The C17O emission is less sensitive to this increased temperature but rather traces the bulk envelope material. The data are only fit by a model where CO is depleted, except in the inner and outermost regions where the temperature increases above 30-40 K. The fact that the temperatures do not drop below approximately 25 K in any of the envelopes whereas a significant fraction of CO is frozen-out suggest that the interstellar radiation field has changed through the evolution of the cores. The H2CO lines are successfully reproduced in the model of an increased ISRF with constant abundances of 3-5x10^{-10}.

Using a fully autonomous spacecraft - Bering - we propose to detect and study sub-km asteroids from an orbit within the asteroid Main Belt. The main purpose of the proposed Bering mission is to detect a statistically significant sample of an expected population of approximately 10^(10) main belt asteroids in the size range 1 m to 1 km. These asteroids are too faint to be observed using Earth-based telescopes. Sub-km asteroids can be detected from spacecraft at close range but due to the high relative velocities and the long communication times this requires a fully autonomous spacecraft. Using theoretical estimates of the distribution and abundance of sub-km asteroids we find that the Bering mission would detect approximately 6 new sub-km asteroids per day. With an expected lifetime for the mission of a few years we expect to detect and study several thousand sub-km asteroids. Results from the Bering mission would allow us to: 1) Provide further links between groups of meteorites and their parent asteroids. 2) Constrain the cratering rate at planetary surfaces and thus allow significantly improved cratering ages for terrains on Mars and other planets. 3) Constrain processes that transfer small asteroids from orbits in the main belt to the inner Solar System.

Infrared ~5--35 um spectra for 40 solar-mass T Tauri stars and 7 intermediate-mass Herbig Ae stars with circumstellar disks were obtained using the Spitzer Space Telescope as part of the c2d IRS survey. This work complements prior spectroscopic studies of silicate infrared emission from disks, which were focused on intermediate-mass stars, with observations of solar-mass stars limited primarily to the 10 um region. The observed 10 and 20 um silicate feature strengths/shapes are consistent with source-to-source variations in grain size. A large fraction of the features are weak and flat, consistent with um-sized grains indicating fast grain growth (from 0.1--1.0 um in radius). In addition, approximately half of the T Tauri star spectra show crystalline silicate features near 28 and 33 um indicating significant processing when compared to interstellar grains. A few sources show large 10-to-20 um ratios and require even larger grains emitting at 20 um than at 10 um. This size difference may arise from the difference in the depth into the disk probed by the two silicate emission bands in disks where dust settling has occurred. The 10 um feature strength vs. shape trend is not correlated with age or Halpha equivalent width, suggesting that some amount of turbulent mixing and regeneration of small grains is occurring. The strength vs. shape trend is related to spectral type, however, with M stars showing significantly flatter 10 um features (larger grain sizes) than A/B stars. The connection between spectral type and grain size is interpreted in terms of the variation in the silicate emission radius as a function of stellar luminosity, but could also be indicative of other spectral-type dependent factors (e.g, X-rays, UV radiation, stellar/disk winds, etc.).

Using the Submillimeter Array we report the discovery of a compact low mass bipolar molecular outflow from L1014-IRS and confirm its association with the L1014 dense core at 200 pc. Consequently, L1014-IRS is the lowest luminosity (L \~0.09 Lsun) and perhaps the lowest mass source known to be driving a bipolar molecular outflow, which is one of the smallest known in size (~500 AU), mass (< 10^{-4} Msun), and energetics (e.g., force < 10^{-7} Msun km/s/yr).

Mid-infrared (23-35 micron) emission from the deeply embedded "Class 0" protostar IRAS16293-2422 is detected with the Spitzer Space Telescope infrared spectrograph. A detailed radiative transfer model reproducing the full spectral energy distribution (SED) from 23 micron to 1.3 mm requires a large inner cavity of radius 600 AU in the envelope to avoid quenching the emission from the central sources. This is consistent with a previous suggestion based on high angular resolution millimeter interferometric data. An alternative interpretation using a 2D model of the envelope with an outflow cavity can reproduce the SED but not the interferometer visibilities. The cavity size is comparable to the centrifugal radius of the envelope and therefore appears to be a natural consequence of the rotation of the protostellar core, which has also caused the fragmentation leading to the central protostellar binary. With a large cavity such as required by the data, the average temperature at a given radius does not increase above 60-80 K and although hot spots with higher temperatures may be present close to each protostar, these constitute a small fraction of the material in the inner envelope. The proposed cavity will also have consequences for the interpretation of molecular line data, especially of complex species probing high temperatures in the inner regions of the envelope.

The Submillimeter Array has opened a new window to study the innermost warm and dense regions of the envelopes and disks around deeply embedded protostars. This paper presents high-angular resolution (< 2'') submillimeter observations of the class 0 young stellar object NGC1333-IRAS2A. Dust continuum emission and lines of complex organic molecules such as CH3OCH3 and CH3OCHO, high excitation CH3OH transitions, deuterated methanol CH3OD as well as lines of CO, HCN, H13CN, SO and SO2 are detected on < 200 AU scales. The observations are interpreted using detailed radiative transfer models of the physical and chemical structure, consistent with both single-dish and interferometer data. The continuum emission is explained by an extended envelope and a compact but resolved component, presumably a circumstellar disk with a diameter of 200-300 AU and a mass of a few times 0.01-0.1 M_sun. If related to the rotation of the envelope, then the size of this disk suggests a centrifugal barrier of 200-300 AU, which implies that the temperature in the envelope does not increase above 100 K. Its large size also suggests that the build-up of disks proceeds rapidly throughout the early protostellar stages. The smaller (< 100 AU) disks found around other deeply embedded protostars may be a result of tidal truncation. The high-resolution observations of SO can be explained with a simple constant abundance, ~1E-9, constrained through single-dish observations, whereas those of H13CN and the organic species require high abundances, increased by one to two orders of magnitude, or an additional compact source of emission at small scales. The compact molecular emission could originate in a hot core region of the inner envelope, but a more likely reservoir is the circumstellar disk.

This paper presents the third in a series of single-dish studies of molecular abundances in the envelopes around a large sample of 18 low-mass pre- and protostellar objects. It focuses on typical grain mantle products and organic molecules, including H2CO, CH3OH and CH3CN. With a few exceptions, all H2CO lines can be fit by constant abundances of 7e(-11)-8e(-9) throughout the envelopes if ortho- and para lines are considered independently. The current observational dataset does not require a large H2CO abundance enhancement in the inner warm regions, but this can also not be ruled out. Through comparison of the H2CO abundances of the entire sample, the H2CO ortho-para ratio is constrained to be 1.6+/-0.3 consistent with thermalization on grains at temperatures of 10-15 K. The H2CO abundances can be related to the empirical chemical network established on the basis of our previously reported survey of other species and is found to be closely correlated with that of the nitrogen-bearing molecules. These correlations reflect the freeze-out of molecules at low temperatures and high densities. An improved fit to the data is obtained with a `drop' abundance structure in which the abundance decrease when the temperature is lower than the evaporation temperature and the density high enough so that the timescale for depletion is less than the lifetime of the core. The derived abundance structures for CH3OH are consistent with a scenario where CH3OH probes the action of compact outflows on the envelopes, which is further supported by comparison to high frequency, high excitation CS J=10-9 and HDO line profiles. The extent to which the outflow dominates the abundance enhancements compared with the passively heated inner envelope depends on the filling factors of the two components in the observing beam.

Spitzer Space Telescope observations of a point-like source, L1014-IRS, close to the dust peak of the low-mass dense core L1014 have questioned its starless nature. The presence of an object with colors expected for an embedded protostar makes L1014-IRS the lowest luminosity isolated protostar known, and an ideal target with which to test star formation theories at the low mass end. In order to study its molecular content and to search for the presence of a molecular outflow, we mapped L1014 in at least one transition of 12CO, N2H+, HCO+, CS and of their isotopologues 13CO, C18O, C17O, N2D+ and H13CO+, using the FCRAO, the IRAM 30 meter and the CSO. The data show physical and chemical properties in L1014 typical of the less evolved starless cores: i.e. H2 central density of a few 10^5 molecules cm^-3, estimated mass of ~2M_sun, CO integrated depletion factor less than 10, N(N2H+)~6*10^12 cm^-2, N(N2D+)/N(N2H+) equal to 10% and relatively broad N2H+(1--0) lines (0.35 km/s). Infall signatures and significant velocity shifts between optically thick and optically thin tracers are not observed in the line profiles. No classical signatures of molecular outflow are found in the 12CO and 13CO observations. In particular, no high velocity wings are found, and no well-defined blue-red lobes of 12CO emission are seen in the channel maps. If sensitive, higher resolution observations confirm the absence of an outflow on a smaller scale than probed by our observations, L1014-IRS would be the only protostellar object known to be formed without driving an outflow.

Radiative transfer models of multi-transition observations are used to determine molecular abundances as functions of position in pre- and protostellar cores. The data require a "drop" abundance profile with radius, with high abundances in the outermost regions probed by low excitation 3 mm lines, and much lower abundances at intermediate zones probed by higher frequency lines. The results are illustrated by detailed analysis of CO and HCO+ lines for a subset of objects. We propose a scenario in which the molecules are frozen out in a region of the envelope where the temperature is low enough (<= 40 K) to prevent immediate desorption, but where the density is high enough (> 10^4-10^5 cm^{-3}) that the freeze-out timescales are shorter than the lifetime of the core. The size of the freeze-out zone is thereby a record of the thermal and dynamical evolution of the cores. Fits to CO data for a sample of 16 objects indicate that the size of the freeze-out zone decreases significantly between Class 0 and I objects, explaining the variations in, for example, CO abundances with envelope masses. However, the corresponding timescales are 10^{5+/-0.5} years, with no significant difference between Class 0 and I objects. These timescales suggest that the dense pre-stellar phase with heavy depletions lasts only a short time, of order 10^5 yr, in agreement with recent chemical-dynamical models.

We have studied the influence of far ultraviolet (UV) radiation from a massive young stellar object (YSO) on the chemistry of its own envelope by extending the models of Doty et al. (2002) to include a central source of UV radiation. The models are applied to the massive star-forming region AFGL 2591 for different inner UV field strengths. Depth-dependent abundance profiles for several molecules are presented and discussed. We predict enhanced column densities for more than 30 species, especially radicals and ions. Comparison between observations and models is improved with a moderate UV field incident on the inner envelope, corresponding to an enhancement factor G0~10-100 at 200 AU from the star with an optical depth tau~15-17. Subtle differences are found compared with traditional models of Photon Dominated Regions (PDRs) because of the higher temperatures and higher gas-phase H2O abundance caused by evaporation of ices in the inner region. In particular, the CN/HCN ratio is not a sensitive tracer of the inner UV field, in contrast with the situation for normal PDRs: for low UV fields, the extra CN reacts with H2 in the inner dense and warm region and produces more HCN. It is found that the CH+ abundance is strongly enhanced and grows steadily with increasing UV field. High-J lines of molecules like CN and HCN are most sensitive to the inner dense region where UV radiation plays a role. Thus, even though the total column density affected by UV photons is small, comparison of high-J and low-J lines can selectively trace and distinguish the inner UV field from the outer one. In addition, future Herschel-HIFI observations of hydrides can sensitively probe the inner UV field.

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.

This paper presents a millimeter wavelength aperture-synthesis study of the spatial variations of the chemistry in the envelope around the deeply embedded low-mass protostar L483-mm on 1000 AU (5") scales. Lines of 8 molecular species including CN, C18O, CS, C34S, HCN, H13CN, HCO+ and N2H+ have been observed using the Owens Valley Radio Observatory Millimeter Array. Continuum emission at 2.7-3.4 millimeter is well-fit by a pure envelope model down to the sensitivity of the interferometer without introducing a disk/compact source, in contrast to what is seen for other protostars. A velocity gradient in the densest material is traced by HCN, CS and N2H+, and is perpendicular to the large-scale CO outflow, with a pattern consistent with rotation around a 1 solar mass central object. Velocity gradients in the propagation direction of the outflow suggest a clear interaction between the outflowing material and "quiescent" core. The C18O observations are fit with a "drop" abundance profile where CO is frozen-out in a region of the envelope with temperatures lower than 40 K and densities higher than 1.5e5 cm^{-3}, which is also required to reproduce single-dish observations. The N2H+ emission strongly resembles that of NH3 and is found to be absent toward the central continuum source. This is a direct consequence of the high CO abundances in the inner region as illustrated by a chemical model for the L483 envelope. The observed CN emission forms a spatial borderline between the outflowing and quiescent material probed by, respectively, HCO+ and N2H+, and also shows intermediate velocities compared to these two species. A scenario is suggested in which CN is enhanced in the walls of an outflow cavity due to the impact of UV irradiation either from the central protostellar system or related to shocks caused by the outflow.

High angular resolution H2CO 218 GHz line observations have been carried out toward the low-mass protostars IRAS 16293-2422 and L1448-C using the Owens Valley Millimeter Array at ~2" resolution. Simultaneous 1.37 mm continuum data reveal extended emission which is compared with that predicted by model envelopes constrained from single-dish data. For L1448-C the model density structure works well down to the 400 AU scale to which the interferometer is sensitive. For IRAS 16293-2422, a known proto-binary object, the interferometer observations indicate that the binary has cleared much of the material in the inner part of the envelope, out to the binary separation of ~800 AU. The morphology and velocity structure of the H2CO array data have been used to investigate whether the abundance enhancements inferred from single-dish modelling are due to thermal evaporation of ices or due to liberation of the ice mantles by shocks in the inner envelope. A scenario in which the H2CO abundance drops in the cold dense part of the envelope where CO is frozen out but is undepleted in the outermost region provides good fits to the single-dish and interferometer data on short baselines for both sources. Emission on the longer baselines is best reproduced if the H2CO abundance is increased by about an order of magnitude from ~1E-10 to ~1E-9 in the inner parts of the envelope due to thermal evaporation when the temperature exceeds ~50K. Other scenarios, including weak outflow-envelope interactions and photon heating of the envelope, are discussed and predictions for future generation interferometers are presented, illustrating their potential in distinguishing these competing scenarios.

This paper presents the first substantial study of the chemistry of the envelopes around a sample of 18 low-mass pre- and protostellar objects for which physical properties have previously been derived from radiative transfer modeling of their dust continuum emission. Single-dish line observations of 24 transitions of 9 molecular species (not counting isotopes) including HCO+, N2H+, CS, SO, SO2, HCN, HNC, HC3N and CN are reported. The line intensities are used to constrain the molecular abundances by comparison to Monte Carlo radiative transfer modeling of the line strengths. An empirical chemical network is constructed on the basis of correlations between the abundances of various species. For example, it is seen that the HCO+ and CO abundances are linearly correlated, both increasing with decreasing envelope mass. Species such as CS, SO and HCN show no trend with envelope mass. In particular no trend is seen between ``evolutionary stage'' of the objects and the abundances of the main sulfur- or nitrogen-containing species. Among the nitrogen-bearing species abundances of CN, HNC and HC3N are found to be closely correlated, which can be understood from considerations of the chemical network. The CS/SO abundance ratio is found to correlate with the abundances of CN and HC3N, which may reflect a dependence on the atomic carbon abundance. An anti-correlation is found between the deuteration of HCO+ and HCN, reflecting different temperature dependences for gas-phase deuteration mechanisms. The abundances are compared to other protostellar environments. In particular it is found that the abundances in the cold outer envelope of the previously studied class 0 protostar IRAS16293-2422 are in good agreement with the average abundances for the presented sample of class 0 objects.

We present a survey of the formaldehyde emission in a sample of eight Class 0 protostars obtained with the IRAM and JCMT telescopes. The data have been analyzed with three different methods with increasing level of sophistication. We first analyze the observed emission in the LTE approximation, and derive rotational temperatures between 11 and 40 K, and column densities between 1 and 20 x 10^13 cm^-2. Second, we use a LVG code and derive larger kinetic temperatures, between 30 and 90 K, consistent with subthermally populated levels and densities from 1 to 6 x 10^5 cm^-3. The column densities from the LVG modeling are within a factor of 10 with respect to those derived in the LTE approximation. Finally, we analyze the observations based upon detailed models for the envelopes surrounding the protostars, using temperature and density profiles previously derived from continuum observations. We approximate the formaldehyde abundance across the envelope with a jump function, the jump occurring when the dust temperature reaches 100 K, the evaporation temperature of the grain mantles. The observed formaldehyde emission is well reproduced only if there is a jump, more than two orders of magnitude, in four sources. In the remaining four sources the data are consistent with a formaldehyde abundance jump, but the evidence is more marginal (~2 sigma). The inferred inner H2CO abundance varies between 1 x 10^-8 and 6 x 10^-6. We discuss the implications of these jumps for our understanding of the origin and evolution of ices in low mass star forming regions. Finally, we give predictions for the submillimeter H2CO lines, which are particularly sensitive to the abundance jumps.

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.

We present observations of the young stars TW Hya and HD 100546 made with the Australia Telescope Compact Array at 89 GHz with $\sim2''$ resolution and $\sim3$ mJy continuum sensitivity. Compact thermal dust continuum emission is detected from disks surrounding both stars. HD 100546 also shows hints of extended emission, presumably a residual protostellar envelope, which is also visible in scattered light at optical wavelengths. For TW Hya, HCO$^+$ J=1--0 line emission from the circumstellar disk is detected and spatially resolved. The observed size and intensity are in good agreement with model calculations based on an irradiated disk with substantial depletions derived previously from single dish observations of higher-J HCO$^+$ transitions.

High angular resolution H2CO line observations have been carried out for the low mass proto-binary star IRAS 16293-2422 using the OVRO millimeter array. Simultaneous continuum observations reveal that the source most probably contains two accretion disks, less than 250 AU in size, through which matter is fed onto the stars. The observations suggest that the binary has cleared most of the material in the inner part of the envelope, out to the binary separation. The H2CO interferometer observations indicate the presence of large scale rotation roughly perpendicular to the large scale CO outflow associated with this source. The H2CO emission is dominated by the compact dense and hot (T>100 K) gas close to the positions of the continuum peaks.

A study of the chemical structure of the envelopes around a sample protostars is introduced. Physical models for the envelopes derived using 1D radiative transfer modeling of their dust continuum emission are used as input for Monte Carlo modeling of single-dish line observations in order to establish the chemical inventories for the different sources. CO and HCO+ abundances are found to be correlated with envelope mass: the sources with the most massive envelopes show the lowest abundances, supporting the idea that the depletion of these molecules is most efficient in the colder, denser environments. Other molecules like CS and HCN do not show a similar trend. Deuterated species like DCO+ and DCN also show signs of correlations with envelope mass and temperature: the deuterium fractionation of HCO+ is most prominent for the coldest and most massive envelopes in the sample and at the same time anti-correlated with the fractionation of HCN. This puts constraints on the low temperature gas-phase deuterium chemistry.

A detailed radiative transfer analysis of the observed continuum and molecular line emission toward the deeply embedded young stellar object IRAS 16293-2422 is performed. The continuum modelling is used to constrain the temperature and density distributions in the envelope, enabling quantitative estimates of various molecular abundances. The molecular excitation analysis reveals that the emission from some molecular species is well reproduced assuming a constant fractional abundance throughout the envelope. The abundances and isotope ratios are generally close to typical values found in cold molecular clouds in these cases, and there is a high degree of deuterium fractionation. There are, however, a number of notable exceptions. Lines covering a wide range of excitation conditions indicate for some molecules, e.g., H2CO, CH3OH, SO, SO2 and OCS, a drastic increase in their abundances in the warm and dense inner region of the circumstellar envelope. The location at which this increase occurs is consistent with the radius at which ices are expected to thermally evaporate off the grains. In all, there is strong evidence for the presence of a `hot core' close to the protostar, whose physical properties are similar to those detected towards most high mass protostars except for a scaling factor. However, the small scale of the hot gas and the infalling nature of the envelope lead to very different chemical time scales between low mass and high mass hot cores, such that only very rapidly produced second-generation complex molecules can be formed in IRAS 16293-2422. Alternatively, the ices may be liberated due to grain-grain collisions in turbulent shear zones where the outflow interacts with the envelope.

We present 1D radiative transfer modelling of the envelopes of a sample of 18 low-mass protostars and pre-stellar cores with the aim of setting up realistic physical models, for use in a chemical description of the sources. The density and temperature profiles of the envelopes are constrained from their radial profiles obtained from SCUBA maps at 450 and 850 micron and from measurements of the source fluxes ranging from 60 micron to 1.3 mm. The densities of the envelopes within ~10000 AU can be described by single power-laws r^{-p} for the class 0 and I sources with p ranging from 1.3 to 1.9, with typical uncertainties of +/- 0.2. Four sources have flatter profiles, either due to asymmetries or to the presence of an outer constant density region. No significant difference is found between class 0 and I sources. The power-law fits fail for the pre-stellar cores, supporting recent results that such cores do not have a central source of heating. The derived physical models are used as input for Monte Carlo modelling of submillimeter C18O and C17O emission. It is found that class I objects typically show CO abundances close to those found in local molecular clouds, but that class 0 sources and pre-stellar cores show lower abundances by almost an order of magnitude implying that significant depletion occurs for the early phases of star formation. While the 2-1 and 3-2 isotopic lines can be fitted using a constant fractional CO abundance throughout the envelope, the 1-0 lines are significantly underestimated, possibly due to contribution of ambient molecular cloud material to the observed emission. The difference between the class 0 and I objects may be related to the properties of the CO ices.