Complex organic molecules with three carbon atoms are found in the earliest stages of star formation. In particular, propenal (C$_2$H$_3$CHO) is a species of interest due to its implication in the formation of more complex species and even biotic molecules. This study aims to search for the presence of C$_2$H$_3$CHO and other three-carbon species such as propylene (C$_3$H$_6$) in the hot corino region of the low-mass protostellar binary IRAS 16293--2422 to understand their formation pathways. We use ALMA observations in Band 6 and 7 from various surveys to search for the presence of C$_3$H$_6$ and C$_2$H$_3$CHO towards the protostar IRAS 16293--2422 B (IRAS 16293B). We report the detection of both C$_3$H$_6$ and C$_2$H$_3$CHO towards IRAS 16293B, however, no unblended lines were found towards the other component of the binary system, IRAS 16293A. We derive column density upper limits for C$_3$H$_8$, HCCCHO, n-C$_3$H$_7$OH, i-C$_3$H$_7$OH, C$_3$O, and cis-HC(O)CHO towards IRAS 16293B. We then use a three-phase chemical model to simulate the formation of these species in a typical prestellar environment followed by its hydrodynamical collapse until the birth of the central protostar. Different formation paths, such as successive hydrogenation and radical-radical additions on grain surfaces, are tested and compared to the observational results. The simulations reproduce the abundances within one order of magnitude from those observed towards IRAS 16293B, with the best agreement found for a rate of $10^{-12}$ cm$^3$ s$^{-1}$ for the gas-phase reaction C$_3$ + O $\rightarrow$ C$_2$ + CO. Successive hydrogenations of C$_3$, HC(O)CHO, and CH$_3$OCHO on grain surfaces are a major and crucial formation route of complex organics molecules, whereas both successive hydrogenation pathways and radical-radical addition reactions contribute to the formation of C$_2$H$_5$CHO.

The Square Kilometre Array (SKA) will answer fundamental questions about the origin, evolution, properties, and influence of magnetic fields throughout the Universe. Magnetic fields can illuminate and influence phenomena as diverse as star formation, galactic dynamics, fast radio bursts, active galactic nuclei, large-scale structure, and Dark Matter annihilation. Preparations for the SKA are swiftly continuing worldwide, and the community is making tremendous observational progress in the field of cosmic magnetism using data from a powerful international suite of SKA pathfinder and precursor telescopes. In this contribution, we revisit community plans for magnetism research using the SKA, in the light of these recent rapid developments. We focus in particular on the impact that new radio telescope instrumentation is generating, thus advancing our understanding of key SKA magnetism science areas, as well as the new techniques that are required for processing and interpreting the data. We discuss these recent developments in the context of the ultimate scientific goals for the SKA era.

This paper presents models to estimate the structure of density and magnetic field strength in spheroidal condensations, from maps of their column density and their polarization of magnetically aligned dust grains. The density model is obtained by fitting a column density map with an embedded p = 2 Plummer spheroid of any aspect ratio and inclination. The magnetic properties are based on the density model, on the Davis-Chandrasekhar-Fermi (DCF) model of Alfv\'enic fluctuations, and on the Spheroid Flux Freezing (SFF) model of mass and flux conservation in Paper I. The field strength model has the resolution of the column density map, which is finer than the resolution of the DCF estimate of field strength. The models are applied to ALMA observations of the envelope of the protostar BHR71 IRS1. Column density fits give the density model, from (2.0 +- 0.4) x 10^5 cm^-3 to (7 +- 1) x 10^7 cm^-3 . The density model predicts the field directions map, which fits the polarization map best within 1100 au, with standard deviation of angle differences 17{\deg}. In this region the DCF mean field strength is 0.7 +- 0.2 mG and the envelope mass is supercritical, with ratio of mass to magnetic critical mass 1.5 +- 0.4. The SFF field strength profile scales with the DCF field strength, from 60 x 10{\mu}G to 4+-1 mG. The spatial resolution of the SFF field strength estimate is finer than the DCF resolution by a factor ~7, and the peak SFF field strength exceeds the DCF field strength by a factor ~5.

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.

The Square Kilometre Array (SKA), currently under design, will be a transformational facility for studying the Universe at centimetre and metre wavelengths in the next decade and beyond. This paper provides the current best estimate of the anticipated performance of SKA Phase 1 (SKA1), using detailed design work, before actual on-sky measurements have been made. It will be updated as new information becomes available. The information contained in this paper takes precedent over any previous documents.

We present and release the full dataset for the Mass Assembly of Stellar Systems and their Evolution with the SMA (MASSES) survey. This survey used the Submillimeter Array (SMA) to image the 74 known protostars within the Perseus molecular cloud. The SMA was used in two array configurations to capture outflows for scales $>$30$^{\prime\prime}$ ($>$9000 au) and to probe scales down to $\sim$1$^{\prime\prime}$ ($\sim$300 au). The protostars were observed with the 1.3 mm and 850 $\mu$m receivers simultaneously to detect continuum at both wavelengths and molecular line emission from CO(2-1), $^{13}$CO(2-1), C$^{18}$O(2-1), N$_2$D$^+$(3-2), CO(3-2), HCO$^+$(4-3), and H$^{13}$CO$^+$(4-3). Some of the observations also used the SMA's recently upgraded correlator, SWARM, whose broader bandwidth allowed for several more spectral lines to be observed (e.g., SO, H$_2$CO, DCO$^+$, DCN, CS, CN). Of the main continuum and spectral tracers observed, 84% of the images and cubes had emission detected. The median C$^{18}$O(2-1) linewidth is $\sim$1.0 km s$^{-1}$, which is slightly higher than those measured with single-dish telescopes at scales of 3000-20000 au. Of the 74 targets, six are suggested to be first hydrostatic core candidates, and we suggest that L1451-mm is the best candidate. We question a previous continuum detection toward L1448 IRS2E. In the SVS13 system, SVS13A certainly appears to be the most evolved source, while SVS13C appears to be hotter and more evolved than SVS13B. The MASSES survey is the largest publicly available interferometric continuum and spectral line protostellar survey to date, and is largely unbiased as it only targets protostars in Perseus. All visibility ($uv$) data and imaged data are publicly available at https://dataverse.harvard.edu/dataverse/full_MASSES/.

We present 1.3 mm ALMA observations of polarized dust emission toward the wide-binary protostellar system BHR 71 IRS1 and IRS2. IRS1 features what appears to be a natal, hourglass-shaped magnetic field. In contrast, IRS2 exhibits a magnetic field that has been affected by its bipolar outflow. Toward IRS2, the polarization is confined mainly to the outflow cavity walls. Along the northern edge of the redshifted outflow cavity of IRS2, the polarized emission is sandwiched between the outflow and a filament of cold, dense gas traced by N$_2$D$^+$, toward which no dust polarization is detected. This suggests that the origin of the enhanced polarization in IRS2 is the irradiation of the outflow cavity walls, which enables the alignment of dust grains with respect to the magnetic field -- but only to a depth of ~300 au, beyond which the dust is cold and unpolarized. However, in order to align grains deep enough in the cavity walls, and to produce the high polarization fraction seen in IRS2, the aligning photons are likely to be in the mid- to far-infrared range, which suggests a degree of grain growth beyond what is typically expected in very young, Class 0 sources. Finally, toward IRS1 we see a narrow, linear feature with a high (10-20%) polarization fraction and a well ordered magnetic field that is not associated with the bipolar outflow cavity. We speculate that this feature may be a magnetized accretion streamer; however, this has yet to be confirmed by kinematic observations of dense-gas tracers.

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

[Abridged] The majority of stars form in binary or higher order systems. The Class 0 protostellar system IRAS16293-2422 contains two protostars, 'A' and 'B', separated by ~600 au and embedded in a single, 10^4 au scale envelope. Their relative evolutionary stages have been debated. We aim to study the relation and interplay between the two protostars A and B at spatial scales of 60 to ~1000 au. We selected molecular gas line transitions of CO, H2CO, HCN, CS, SiO, and CCH from the ALMA-PILS spectral imaging survey (329-363 GHz) and used them as tracers of kinematics, density, and temperature in the IRAS16293-2422 system. The angular resolution of the PILS data set allows us to study these quantities at a resolution of 0.5 arcsec (60 au [..]). Line-of-sight velocity maps of both optically thick and optically thin molecular lines reveal: (i) new manifestations of previously known outflows emanating from protostar A; (ii) a kinematically quiescent bridge of dust and gas spanning between the two protostars, with an inferred density between 4 10^4 and 3 10^7 cm^-3; and (iii) a separate, straight filament seemingly connected to protostar B seen only in CCH, with a flat kinematic signature. Signs of various outflows, all emanating from source A, are evidence of high-density and warmer gas; none of them coincide spatially and kinematically with the bridge. We hypothesize that the bridge arc is a remnant of filamentary substructure in the protostellar envelope material from which protostellar sources A and B have formed. One particular morphological structure appears to be due to outflowing gas impacting the quiescent bridge material. The continuing lack of clear outflow signatures unambiguously associated to protostar B and the vertically extended shape derived for its disk-like structure lead us to conclude that source B may be in an earlier evolutionary stage than source A.

In the standard picture for low-mass star formation, a dense molecular cloud undergoes gravitational collapse to form a protostellar system consisting of a new central star, a circumstellar disk, and a surrounding envelope of remaining material. The mass distribution of the system evolves as matter accretes from the large-scale envelope through the disk and onto the protostar. While this general picture is supported by simulations and indirect observational measurements, the specific timescales related to disk growth and envelope dissipation remain poorly constrained. In this paper we conduct a rigorous test of a method introduced by J{\o}rgensen et al. (2009) to obtain mass measurements of disks and envelopes around embedded protostars with observations that do not resolve the disk (resolution of $\sim$1000\,AU). Using unresolved data from the recent Mass Assembly of Stellar Systems and their Evolution with the SMA (MASSES) survey, we derive disk and envelope mass estimates for $59$ protostellar systems in the Perseus molecular cloud. We compare our results to independent disk mass measurements from the VLA Nascent Disk and Multiplicity (VANDAM) survey and find a strong linear correlation, suggesting that accurate disk masses can be measured from unresolved observations. Then, leveraging the size of the MASSES sample, we find no significant trend in protostellar mass distribution as a function of age, as approximated from bolometric temperatures. These results may indicate that the disk mass of a protostar is set near the onset of the Class 0 protostellar stage and remains roughly constant throughout the Class I protostellar stage.

Stars like our Sun form in self-gravitating dense and cold structures within interstellar clouds, called pre-stellar cores. Although much is known about the physical structure of dense clouds just before and soon after the switch-on of a protostar, the central few thousand astronomical units (au) of pre-stellar cores are unexplored. It is within these central regions that stellar systems assemble and fragmentation may take place, with the consequent formation of binaries and multiple systems. We present ALMA Band 6 observations (ACA and 12m array) of the dust continuum emission of the 8 Msun pre-stellar core L1544, with angular resolution of 2'' x 1.6'' (linear resolution 270 au x 216 au). Within the primary beam, a compact region of 0.1 Msun, which we call a "kernel", has been unveiled. The kernel is elongated, with a central flat zone with radius Rker ~ 10'' (~ 1400 au). The average number density within Rker is ~1 x 10^6 cm^{-3}, with possible local density enhancements. The region within Rker appears to have fragmented, but detailed analysis shows that similar substructure can be reproduced by synthetic interferometric observations of a smooth centrally concentrated dense core with a similar central flat zone. The presence of a smooth kernel within a dense core is in agreement with non-ideal magneto-hydro-dynamical simulations of a contracting cloud core with a peak number density of 1 x 10^7 cm^{-3}. Dense cores with lower central densities are completely filtered out when simulated 12m-array observations are carried out. These observations demonstrate that the kernel of dynamically evolved dense cores can be investigated at high angular resolution with ALMA.

A good constraint of when the growth of dust grains from sub-micrometer to millimeter sizes occurs, is crucial for planet formation models. This provides the first step towards the production of pebbles and planetesimals in protoplanetary disks. Currently, it is well established that Class II objects have large dust grains. However, it is not clear when in the star formation process this grain growth occurs. We use multi-wavelength millimeter observations of a Class I protostar to obtain the spectral index of the observed flux densities $\alpha_\mathrm{mm}$ of the unresolved disk and the surrounding envelope. Our goal is to compare our observational results with visibility modeling at both wavelengths simultaneously. We present data from NOEMA at 2.7 mm and SMA at 1.3 mm of the Class I protostar, Per-emb-50. We model the dust emission with a variety of parametric and radiative transfer models to deduce the grain size from the observed emission spectral index. We find a spectral index in the envelope of Per-emb-50 of $\alpha_{\rm env}$=$3.3\pm0.3$, similar to the typical ISM values. The radiative transfer modeling of the source confirms this value of $\alpha_{\rm env}$ with the presence of dust with a $a_\mathrm{max}$$\leq$100 $\mu$m. Additionally, we explore the backwarming effect, where we find that the envelope structure affects the millimeter emission of the disk. Our results reveal grains with a maximum size no larger than $100$ $\mu$m in the inner envelope of the Class I protostar Per-emb-50, providing an interesting case to test the universality of millimeter grain growth expected in these sources.

Studies of deuterated isotopologues of complex organic molecules can provide important constraints on their origin in regions of star formation. In particular, the abundances of deuterated species are very sensitive to the physical conditions in the environment where they form. Due to the low temperatures in regions of star formation, these isotopologues are enhanced to significant levels, making detections of multiply-deuterated species possible. However, for complex organic species, only the multiply-deuterated variants of methanol and methyl cyanide have been reported so far. The aim of this paper is to initiate the characterisation of multiply-deuterated variants of complex organic species with the first detection of doubly-deuterated methyl formate, CHD2OCHO. We use ALMA observations from the Protostellar Interferometric Line Survey (PILS) of the protostellar binary IRAS 16293-2422, in the spectral range of 329.1 GHz to 362.9 GHz. We report the first detection of doubly-deuterated methyl formate CHD2OCHO in the ISM. The D/H ratio of CHD2OCHO is found to be 2-3 times higher than the D/H ratio of CH2DOCHO for both sources, similar to the results for formaldehyde from the same dataset. The observations are compared to a gas-grain chemical network coupled to a dynamical physical model, tracing the evolution of a molecular cloud until the end of the Class 0 protostellar stage. The overall D/H ratio enhancements found in the observations are of the same order of magnitude as the predictions from the model for the early stages of Class 0 protostars. However, the higher D/H ratio of CHD2OCHO compared to the D/H ratio of CH2DOCHO is still not predicted by the model. This suggests that a mechanism is enhancing the D/H ratio of singly- and doubly-deuterated methyl formate that is not in the model, e.g. mechanisms for H-D substitutions.

We present Atacama Large Millimeter/submillimeter Array (ALMA) observations of 1.1 mm dust continuum and CO 2-1 emission toward six dense cores within the Ophiuchus molecular cloud. We detect compact, sub-arcsecond continuum structures toward three targets, two of which (Oph A N6 and SM1) are located in the Ophiuchus A ridge. Two targets, SM1 and GSS 30, contain two compact sources within the ALMA primary beam. We argue that several of the compact structures are small ($R \lesssim 80$ au) accretion disks around young protostars, due to their resolved, elongated structures, coincident radio and x-ray detections, or bipolar outflow detections. While CO line wings extend to $\pm 10-20$ km s$^{-1}$ for the more evolved sources GSS 30 IRS3 and IRS1, CO emission toward other sources, where detected, only extends a few km s$^{-1}$ from the cloud $v_\mathrm{LSR}$. The dust spectral index toward the compact objects suggests that the disks are either optically thick at 1.1 mm, or that significant grain growth has already occurred. We identify, for the first time, a single compact continuum source ($R \sim 100$ au) toward N6 embedded within a larger continuum structure. SM1N is extended in the continuum but is highly centrally concentrated, with a density profile that follows a $r^{-1.3}$ power law within 200 au, and additional structure suggested by the uv-data. Both N6 and SM1N show no clear bipolar outflows with velocities greater than a few km s$^{-1}$ from the cloud velocity. These sources are candidates to be the youngest protostars or first hydrostatic cores in the Ophiuchus molecular cloud.

We present a characterization of the binary protostar system that is forming within a dense core in the isolated dark cloud BHR71. The pair of protostars, IRS1 and IRS2, are both in the Class 0 phase, determined from observations that resolve the sources from 1 um out to 250 um and from 1.3 mm to 1.3cm. The resolved observations enable the luminosities of IRS1 and IRS2 to be independently measured (14.7 and 1.7L_sun, respectively), in addition to the bolometric temperatures 68~K, and 38~K, respectively. The surrounding core was mapped in NH3 (1,1) with the Parkes radio telescope, and followed with higher-resolution observations from ATCA in NH3 (1,1) and 1.3cm continuum. The protostars were then further characterized with ALMA observations in the 1.3~mm continuum along with N2D+ (J=3-2), 12CO, 13CO, and C18O (J=2-1) molecular lines. The Parkes observations find evidence for a velocity gradient across the core surrounding the two protostars, while ATCA reveals more complex velocity structure toward the protostars within the large-scale gradient. The ALMA observations then reveal that the two protostars are at the same velocity in C18O, and N2H+ exhibits a similar velocity structure as NH3. However, the C18O kinematics reveal that the rotation on scales $<$1000~AU around IRS1 and IRS2 are in opposite directions. Taken with the lack of a systematic velocity difference between the pair, it is unlikely that their formation resulted from rotational fragmentation. We instead conclude that the binary system most likely formed via turbulent fragmentation of the core.

Continuum observations of molecular clouds have revealed a surprising amount of substructure in the form of filaments of a few pc length and cores of ~0.1 pc diameter. Understanding the evolution of these substructures towards star formation requires the kinematic and dynamical insights provided uniquely by sensitive line observations at high angular and spectral resolution. In this short paper, we describe how an ngVLA can probe effectively the dynamics of filaments and cores in nearby star-forming molecular clouds using the NH3 rotation-inversion transitions at 24 GHz. Such emission has been proven to trace well the high column density environments of star-forming cores and filaments but higher-resolution observations are needed to reveal important details of how dense gas is flowing within and onto these substructures. In particular, we describe how 150 x 18-m antennas with a maximum baseline of 1 km can be used to map sensitively NH3 emission across high column density locations in clouds in roughly an order of magnitude less time than with the current Jansky VLA.

Hydroxylamine (NH$_{2}$OH) and methylamine (CH$_{3}$NH$_{2}$) have both been suggested as precursors to the formation of amino acids and are therefore of interest to prebiotic chemistry. Their presence in interstellar space and formation mechanisms, however, are not well established. We aim to detect both amines and their potential precursor molecules NO, N$_{2}$O and CH$_{2}$NH towards the low-mass protostellar binary IRAS 16293--2422, in order to investigate their presence and constrain their interstellar formation mechanisms around a young Sun-like protostar. ALMA observations from the unbiased, high angular resolution and sensitivity Protostellar Interferometric Line Survey (PILS) are used. Spectral transitions of the molecules under investigation are searched for with the CASSIS line analysis software. CH$_2$NH and N$_{2}$O are detected for the first time towards a low-mass source, the latter molecule through confirmation with the single-dish TIMASSS survey. NO is also detected. CH$_{3}$NH$_{2}$ and NH$_{2}$OH are not detected and stringent upper limit column densities are determined. The non-detection of CH$_{3}$NH$_{2}$ and NH$_{2}$OH limits the importance of formation routes to amino acids involving these species. The detection of CH$_{2}$NH makes amino acid formation routes starting from this molecule plausible. The low abundances of CH$_2$NH and CH$_{3}$NH$_{2}$ compared to Sgr B2 indicate that different physical conditions influence their formation in low- and high-mass sources.

Methyl isocyanide (CH$_3$NC) is the isocyanide with the largest number of atoms confirmed in the interstellar medium (ISM), but it is not an abundant molecule, having only been detected towards a handful of objects. Conversely, its isomer, methyl cyanide (CH$_3$CN), is one of the most abundant complex organic molecules detected in the ISM, with detections in a variety of low- and high-mass sources. We use ALMA observations from the Protostellar Interferometric Line Survey (PILS) to search for methyl isocyanide and compare its abundance with that of its isomer methyl cyanide. We use a new line catalogue from the Cologne Database for Molecular Spectroscopy (CDMS) to identify methyl isocyanide lines. We also model the chemistry with an updated version of the three-phase chemical kinetics model {\em MAGICKAL}, presenting the first chemical modelling of methyl isocyanide to date. We detect methyl isocyanide for the first time in a solar-type protostar, IRAS 16293$-$2422 B, and present upper limits for its companion protostar, IRAS 16293$-$2422 A. Methyl isocyanide is found to be at least 20 times more abundant in source B compared to source A, with a CH$_3$CN/CH$_3$NC abundance ratio of 200 in IRAS 16293--2422 B and >5517 in IRAS 16293$-$2422 A. We also present the results of a chemical model of methyl isocyanide chemistry in both sources, and discuss the implications on methyl isocyanide formation mechanisms and the relative evolutionary stages of both sources.

We present the Mass Assembly of Stellar Systems and their Evolution with the SMA (MASSES) survey, which uses the Submillimeter Array (SMA) interferometer to map the continuum and molecular lines for all 74 known Class 0/I protostellar systems in the Perseus molecular cloud. The primary goal of the survey is to observe an unbiased sample of young protostars in a single molecular cloud so that we can characterize the evolution of protostars. This paper releases the MASSES 1.3 mm data from the subcompact configuration ($\sim$4$^{\prime\prime}$ or $\sim$1000 au resolution), which is the SMA's most compact array configuration. We release both $uv$ visibility data and imaged data for the spectral lines CO(2-1),$^{13}$CO(2-1), C$^{18}$O(2-1), and N$_2$D$^+$(3-2), as well as for the 1.3 mm continuum. We identify the tracers that are detected toward each source. We also show example images of continuum and CO(2-1) outflows, analyze C$^{18}$O(2-1) spectra, and present data from the SVS 13 star-forming region. The calculated envelope masses from the continuum show a decreasing trend with bolometric temperature (a proxy for age). Typical C$^{18}$O(2-1) linewidths are 1.45 km s$^{-1}$, which is higher than the C$^{18}$O linewidths detected toward Perseus filaments and cores. We find that N$_2$D$^+$(3-2) is significantly more likely to be detected toward younger protostars. We show that the protostars in SVS 13 are contained within filamentary structures as traced by C$^{18}$O(2-1) and N$_2$D$^+$(3-2). We also present the locations of SVS 13A's high velocity (absolute line-of-sight velocities $>$150 km s$^{-1}$) red and blue outflow components. Data can be downloaded from https://dataverse.harvard.edu/dataverse/MASSES .

The magnetic field plays an important role in every stage of the star-formation process from the collapse of the initial protostellar core to the star's arrival on the main sequence. Consequently, the goal of this science case is to explore a wide range of magnetic phenomena that can be investigated using the polarization capabilities of the Next Generation Very Large Array (ngVLA). These include (1) magnetic fields in protostellar cores via polarized emission from aligned dust grains, including in regions optically thick at wavelengths observable by the Atacama Large Millimeter/submillimeter Array (ALMA); (2) magnetic fields in both protostellar cores and molecular outflows via spectral-line polarization from the Zeeman and Goldreich-Kylafis effects; (3) magnetic fields in protostellar jets via polarized synchrotron emission; and (4) gyrosynchrotron emission from magnetospheres around low-mass stars.

Complex organic molecules are readily detected in the inner regions of the gaseous envelopes of forming protostars. In particular, molecules that contain nitrogen are interesting due to the role nitrogen plays in the development of life and the compact scales such molecules have been found to trace around forming protostars. The goal of this work is to determine the inventory of one family of nitrogen-bearing organic molecules, complex nitriles (molecules with a $-$CN functional group) towards two hot corino sources in the low-mass protostellar binary IRAS 16293$-$2422. This work explores the abundance differences between the two sources, the isotopic ratios, and the spatial extent derived from molecules containing the nitrile functional group. Using data from the Protostellar Interferometric Line Survey (PILS) obtained with ALMA we determine abundances and excitation temperatures for the detected nitriles. We also present a new method for determining the spatial structure of sources with high line density and large velocity gradients $-$ Velocity-corrected INtegrated emission (VINE) maps. We detect methyl cyanide (CH$_3$CN) as well as 5 of its isotopologues, including the detection of CHD$_2$CN which is the first detection in the ISM. We also detect ethyl cyanide (C$_2$H$_5$CN), vinyl cyanide (C$_2$H$_3$CN), and cyanoacetylene (HC$_3$N). We find that abundances are similar between IRAS 16293A and IRAS 16293B on small scales except for vinyl cyanide which is only detected towards the latter source. This suggests an important difference between the sources either in their evolutionary stage or warm-up timescales. We also detect a spatially double-peaked emission for the first time in molecular emission in the A source, suggesting that this source is showing structure related to a rotating toroid of material.

We present a Chandra study of the deeply embedded Serpens South star-forming region, examining cluster structure and disk properties at the earliest stages. In total, 152 X-ray sources are detected. Combined with Spitzer and 2MASS photometry, 66 X-ray sources are reliably matched to an IR counterpart. We identify 21 class I, 6 flat spectrum, 16 class II, and 18 class III young stars; 5 were unclassified. Eighteen sources were variable in X-rays, 8 exhibiting flare-like emission, and one periodic source. The cluster X-ray luminosity distance was estimated, the best match was to the nearer distance of 260pc for the front of the Aquila Rift complex. The $N_{H}$ vs. $A_{K}$ ratio is found to be $\sim$0.68x10$^{22}$, similar to that measured in other young low mass regions, but lower than that measured in the ISM and high mass clusters ($\sim$1.6-2x10$^{22}$). We find the spatial distribution closely follows that of the dense filament from which the stars have formed, with the class II population still strongly associated with the filament. There are four sub-clusters in the field, with three forming knots in the filament, and a fourth to the west, which may not be associated but may be contributing to the distributed class III population. A high percentage of diskless class IIIs (upper limit 30% of classified X-ray sources) in such a young cluster could indicate that processing of disks is influenced by the cluster environment and is not solely time-scale dependent.

We present a characterization of the protostar embedded within the BHR7 dark cloud, based on both photometric measurements from the near-infrared to millimeter and interferometric continuum and molecular line observations at millimeter wavelengths. We find that this protostar is a Class 0 system, the youngest class of protostars, measuring its bolometric temperature to be 50.5~K, with a bolometric luminosity of 9.3~L$_{\odot}$. The near-infrared and \textit{Spitzer} imaging show a prominent dark lane from dust extinction separating clear bipolar outflow cavities. Observations of $^{13}$CO ($J=2\rightarrow1$), C$^{18}$O ($J=2\rightarrow1$), and other molecular lines with the Submillimeter Array (SMA) exhibit a clear rotation signature on scales $<$1300~AU. The rotation can be traced to an inner radius of $\sim$170~AU and the rotation curve is consistent with an R$^{-1}$ profile, implying that angular momentum is being conserved. Observations of the 1.3~mm dust continuum with the SMA reveal a resolved continuum source, extended in the direction of the dark lane, orthogonal to the outflow. The deconvolved size of the continuum indicates a radius of $\sim$100~AU for the continuum source at the assumed distance of 400~pc. The visibility amplitude profile of the continuum emission cannot be reproduced by an envelope alone and needs a compact component. Thus, we posit that the resolved continuum source could be tracing a Keplerian disk in this very young system. If we assume that the continuum radius traces a Keplerian disk (R$\sim$120~AU) the observed rotation profile is consistent with a protostar mass of 1.0~$M_{\odot}$.

Searches for the prebiotically-relevant cyanamide (NH$_2$CN) towards solar-type protostars have not been reported in the literature. We here present the first detection of this species in the warm gas surrounding two solar-type protostars, using data from the Atacama Large Millimeter/Submillimeter Array Protostellar Interferometric Line Survey (PILS) of IRAS 16293-2422 B and observations from the IRAM Plateau de Bure Interferometer of NGC1333 IRAS2A. We furthermore detect the deuterated and $^{13}$C isotopologues of NH$_2$CN towards IRAS 16293-2422 B. This is the first detection of NHDCN in the interstellar medium. Based on a local thermodynamic equilibrium analysis, we find that the deuteration of cyanamide ($\sim$ 1.7%) is similar to that of formamide (NH$_2$CHO), which may suggest that these two molecules share NH$_2$ as a common precursor. The NH$_2$CN/NH$_2$CHO abundance ratio is about 0.2 for IRAS 16293-2422 B and 0.02 for IRAS2A, which is comparable to the range of values found for Sgr B2. We explored the possible formation of NH$_2$CN on grains through the NH$_2$ + CN reaction using the chemical model MAGICKAL. Grain-surface chemistry appears capable of reproducing the gas-phase abundance of NH$_2$CN with the correct choice of physical parameters.

Context. The Class 0 protostellar binary IRAS 16293-2422 is an interesting target for (sub)millimeter observations due to, both, the rich chemistry toward the two main components of the binary and its complex morphology. Its proximity to Earth allows the study of its physical and chemical structure on solar system scales using high angular resolution observations. Such data reveal a complex morphology that cannot be accounted for in traditional, spherical 1D models of the envelope. Aims. The purpose of this paper is to study the environment of the two components of the binary through 3D radiative transfer modeling and to compare with data from the Atacama Large Millimeter/submillimeter Array. Such comparisons can be used to constrain the protoplanetary disk structures, the luminosities of the two components of the binary and the chemistry of simple species. Methods. We present 13CO, C17O and C18O J=3-2 observations from the ALMA Protostellar Interferometric Line Survey (PILS), together with a qualitative study of the dust and gas density distribution of IRAS 16293-2422. A 3D dust and gas model including disks and a dust filament between the two protostars is constructed which qualitatively reproduces the dust continuum and gas line emission. Results and conclusions. Radiative transfer modeling of source A and B, with the density solution of an infalling, rotating collapse or a protoplanetary disk model, can match the constraints for the disk-like emission around source A and B from the observed dust continuum and CO isotopologue gas emission. If a protoplanetary disk model is used around source B, it has to have an unusually high scale-height in order to reach the dust continuum peak emission value, while fulfilling the other observational constraints. Our 3D model requires source A to be much more luminous than source B; LA ~ 18 $L_\odot$ and LB ~ 3 $L_\odot$.

We present a study of hierarchical structure in the Perseus molecular cloud, from the scale of the entire cloud ($\gtrsim$10 pc) to smaller clumps ($\sim$1 pc), cores ($\sim$0.05-0.1 pc), envelopes ($\sim$300-3000 AU) and protostellar objects ($\sim$15 AU). We use new observations from the Submillimeter Array (SMA) large project "Mass Assembly of Stellar Systems and their Evolution with the SMA (MASSES)" to probe the envelopes, and recent single-dish and interferometric observations from the literature for the remaining scales. This is the first study to analyze hierarchical structure over five scales in the same cloud complex. We compare the number of fragments with the number of Jeans masses in each scale to calculate the Jeans efficiency, or the ratio of observed to expected number of fragments. The velocity dispersion is assumed to arise either from purely thermal motions, or from combined thermal and non-thermal motions inferred from observed spectral line widths. For each scale, thermal Jeans fragmentation predicts more fragments than observed, corresponding to inefficient thermal Jeans fragmentation. For the smallest scale, thermal plus non-thermal Jeans fragmentation also predicts too many protostellar objects. However at each of the larger scales thermal plus non-thermal Jeans fragmentation predicts fewer than one fragment, corresponding to no fragmentation into envelopes, cores, and clumps. Over all scales, the results are inconsistent with complete Jeans fragmentation based on either thermal or thermal plus non-thermal motions. They are more nearly consistent with inefficient thermal Jeans fragmentation, where the thermal Jeans efficiency increases from the largest to the smallest scale.

We present intensity-corrected Herschel maps at 100 um, 160 um, 250 um, 350 um, and 500 um for 56 isolated low-mass clouds. We determine the zero-point corrections for Herschel PACS and SPIRE maps from the Herschel Science Archive (HSA) using Planck data. Since these HSA maps are small, we cannot correct them using typical methods. Here, we introduce a technique to measure the zero-point corrections for small Herschel maps. We use radial profiles to identify offsets between the observed HSA intensities and the expected intensities from Planck. Most clouds have reliable offset measurements with this technique. In addition, we find that roughly half of the clouds have underestimated HSA-SPIRE intensities in their outer envelopes relative to Planck, even though the HSA-SPIRE maps were previously zero-point corrected. Using our technique, we produce corrected Herschel intensity maps for all 56 clouds and determine their line-of-sight average dust temperatures and optical depths from modified black body fits. The clouds have typical temperatures of ~ 14-20 K and optical depths of ~ 1e-5 - 1e-3. Across the whole sample, we find an anti-correlation between temperature and optical depth. We also find lower temperatures than what was measured in previous Herschel studies, which subtracted out a background level from their intensity maps to circumvent the zero-point correction. Accurate Herschel observations of clouds are key to obtain accurate density and temperature profiles. To make such future analyses possible, intensity-corrected maps for all 56 clouds are publicly available in the electronic version.

[abridged] The enhanced degrees of deuterium fractionation observed in envelopes around protostars demonstrate the importance of chemistry at low temperatures, relevant in pre- and protostellar cores. formaldehyde is an important species in the formation of methanol and more complex molecules. Here, we present the first study of formaldehyde deuteration on small scales around the prototypical low-mass protostar IRAS 16293-2422 using high spatial and spectral resolution Atacama Large Millimeter/submillimeter Array (ALMA) observations. Numerous isotopologues of formaldehyde are detected, among them H$_2$C$^{17}$O, and D$_2^{13}$CO for the first time in the ISM. The large range of upper energy levels covered by the HDCO lines help constrain the excitation temperature to 106$\pm$13 K. Using the derived column densities, formaldehyde shows a deuterium fractionation of HDCO/H$_2$CO=6.5$\pm$1%, D$_2$CO/HDCO=12.8$^{+3.3}_{-4.1}$%, and D$_2$CO/H$_2$CO=0.6(4)$\pm$0.1%. The isotopic ratios derived are $^{16}$O/$^{18}$O=805, $^{18}$O/$^{17}$O=3.2 and $^{12}$C/$^{13}$C=56. The HDCO/H$_2$CO ratio is lower than found in previous studies, highlighting the uncertainties involved in interpreting single dish observations of the inner warm regions. The D$_2$CO/HDCO ratio is only slightly larger than the HDCO/H$_2$CO ratio. This is consistent with formaldehyde forming in the ice as soon as CO has frozen onto the grains, with most of the deuteration happening towards the end of the prestellar core phase. A comparison with available time-dependent chemical models indicates that the source is in the early Class 0 stage.

Young stellar objects (YSOs) may undergo periods of active accretion (outbursts), during which the protostellar accretion rate is temporarily enhanced by a few orders of magnitude. Whether or not these accretion outburst YSOs possess similar dust/gas reservoirs to each other, and whether or not their dust/gas reservoirs are similar as quiescent YSOs, are issues not yet clarified. The aim of this work is to characterize the millimeter thermal dust emission properties of a statistically significant sample of long and short duration accretion outburst YSOs (i.e., FUors and EXors) and the spectroscopically identified candidates of accretion outbursting YSOs (i.e., FUor-like objects). We have carried out extensive Submillimeter Array (SMA) observations mostly at $\sim$225 GHz (1.33 mm) and $\sim$272 GHz (1.10 mm), from 2008 to 2017. We covered accretion outburst YSOs located at $<$1 kpc distances from the solar system. We analyze all the existing SMA data of such objects, both published and unpublished, in a coherent way to present a millimeter interferometric database of 29 objects. We obtained 21 detections at $>$3-$\sigma$ significance. Detected sources except for the two cases of V883 Ori and NGC 2071 MM3 were observed with $\sim$1$"$ angular resolution. Overall our observed targets show a systematically higher millimeter luminosity distribution than those of the $M_{*}>$0.3 $M_{\odot}$ Class II YSOs in the nearby ($\lesssim$400 pc) low-mass star-forming molecular clouds (e.g., Taurus, Lupus, Upp Scorpio, and Chameleon I). In addition, at 1 mm our observed confirmed binaries or triple-system sources are systematically fainter than the rest of the sources even though their 1 mm fluxes are broadly distributed. We may have detected $\sim$30-60\% millimeter flux variability from V2494 Cyg and V2495 Cyg, from the observations separated by $\sim$1 year.

We present new Submillimeter Array (SMA) observations of CO(2-1) outflows toward young, embedded protostars in the Perseus molecular cloud as part of the Mass Assembly of Stellar Systems and their Evolution with the SMA (MASSES) survey. For 57 Perseus protostars, we characterize the orientation of the outflow angles and compare them with the orientation of the local filaments as derived from $Herschel$ observations. We find that the relative angles between outflows and filaments are inconsistent with purely parallel or purely perpendicular distributions. Instead, the observed distribution of outflow-filament angles are more consistent with either randomly aligned angles or a mix of projected parallel and perpendicular angles. A mix of parallel and perpendicular angles requires perpendicular alignment to be more common by a factor of $\sim$3. Our results show that the observed distributions probably hold regardless of the protostar's multiplicity, age, or the host core's opacity. These observations indicate that the angular momentum axis of a protostar may be independent of the large-scale structure. We discuss the significance of independent protostellar rotation axes in the general picture of filament-based star formation.

Context: Understanding how accretion proceeds is a key question of star formation, with important implications for both the physical and chemical evolution of young stellar objects. In particular, very little is known about the accretion variability in the earliest stages of star formation. Aims: To characterise protostellar accretion histories towards individual sources by utilising sublimation and freeze-out chemistry of CO. Methods: A sample of 24 embedded protostars are observed with the Submillimeter Array (SMA) in context of the large program "Mass Assembly of Stellar Systems and their Evolution with the SMA" (MASSES). The size of the C$^{18}$O emitting region, where CO has sublimated into the gas-phase, is measured towards each source and compared to the expected size of the region given the current luminosity. The SMA observations also include 1.3 mm continuum data, which are used to investigate whether a link can be established between accretion bursts and massive circumstellar disks. Results: Depending on the adopted sublimation temperature of the CO ice, between 20% and 50% of the sources in the sample show extended C$^{18}$O emission indicating that the gas was warm enough in the past that CO sublimated and is currently in the process of refreezing; something which we attribute to a recent accretion burst. Given the fraction of sources with extended C$^{18}$O emission, we estimate an average interval between bursts of 20000-50000 yr, which is consistent with previous estimates. No clear link can be established between the presence of circumstellar disks and accretion bursts, however the three closest known binaries in the sample (projected separations <20 AU) all show evidence of a past accretion burst, indicating that close binary interactions may also play a role in inducing accretion variability.

We present observations and analysis of the low-mass star-forming region, Taurus Molecular Cloud-1 (TMC-1). CS ($J$=2-1)/N$_2$H$^+$ ($J$=1-0) and C$^{17}$O ($J$=2-1)/C$^{18}$O ($J$=2-1) were observed with FCRAO (Five College Radio Astronomy Observatory) and SRAO (Seoul Radio Astronomy Observatory), respectively. In addition, Spitzer infrared data and 1.2 mm continuum data observed with MAMBO (Max-Planck Millimetre Bolometer) are used. We also perform chemical modeling to investigate the relative molecular distributions of the TMC-1 filament. Based on Spitzer observations, there is no young stellar object along the TMC-1 filament, while five Class II and one Class I young stellar objects are identified outside the filament. The comparison between column densities calculated from dust continuum and C$^{17}$O 2-1 line emission shows that CO is depleted much more significantly in the ammonia peak than in the cyanopolyyne peak, while the column densities calculated from the dust continuum are similar at the two peaks. N$_2$H$^+$ is not depleted much in either peak. According to our chemical calculation, the differential chemical distribution in the two peaks can be explained by different timescales required to reach the same density, i.e., by different dynamical processes.

Compact substructure is expected to arise in a starless core as mass becomes concentrated in the central region likely to form a protostar. Additionally, multiple peaks may form if fragmentation occurs. We present ALMA Cycle 2 observations of 60 starless and protostellar cores in the Ophiuchus molecular cloud. We detect eight compact substructures which are >15 arcsec from the nearest Spitzer YSO. Only one of these has strong evidence for being truly starless after considering ancillary data, e.g., from Herschel and X-ray telescopes. An additional extended emission structure has tentative evidence for starlessness. The number of our detections is consistent with estimates from a combination of synthetic observations of numerical simulations and analytical arguments. This result suggests that a similar ALMA study in the Chamaeleon I cloud, which detected no compact substructure in starless cores, may be due to the peculiar evolutionary state of cores in that cloud.

One of the open questions in astrochemistry is how complex organic and prebiotic molecules are formed. Aims. Our aim is to start the process of compiling an inventory of oxygen-bearing complex organic molecules toward the solar-type Class 0 protostellar binary IRAS16293-2422 from an unbiased spectral survey with ALMA (PILS). Here we focus on the new detections of ethylene oxide (c-C$_2$H$_4$O), acetone (CH$_3$COCH$_3$), and propanal (C$_2$H$_5$CHO). Methods. With ALMA, we surveyed the spectral range from 329 to 363 GHz at 0.5$"$ (60 AU diameter) resolution. Using a simple model for the molecular emission in LTE, the excitation temperatures and column densities of each species were constrained. Results. We successfully detect propanal (44 lines), ethylene oxide (20 lines) and acetone (186 lines) toward one component of the protostellar binary, IRAS16293B. The high resolution maps demonstrate that the emission for all investigated species originates from the compact central region close to the protostar. This, along with a derived common excitation temperature of $\sim$ 125 K, is consistent with a coexistence of these molecules in the same gas. Conclusions. The observations mark the first detections of acetone, propanal and ethylene oxide toward a low-mass protostar. The relative abundance ratios of the two sets of isomers (CH$_3$COCH$_3$/C$_2$H$_5$CHO $\sim$ 8 and CH$_3$CHO/c-C$_2$H$_4$O $\sim$ 12) are comparable to previous observations toward high-mass protostars. The majority of observed abundance ratios from these results as well as those measured toward high-mass protostars are up to an order of magnitude above the predictions from chemical models. This may reflect either missing reactions or uncertain rates in the chemical networks. The physical conditions, such as temperatures or densities, used in the models, may not be applicable to solar-type protostars either.

Filamentary structures are ubiquitous in molecular clouds, and have been recently argued to play an important role in regulating the size and mass of embedded clumps through fragmentation and mass accretion. Here, we reveal the dynamical state and fragmentation of filamentary molecular gas associated with the Serpens South protocluster through analysis of wide (~4 x 4 pc) observations of NH3 (1,1) and (2,2) inversion transitions with the Green Bank Telescope. Detailed modeling of the NH3 lines reveals that the kinematics of the cluster and surrounding filaments are complex. We identify hierarchical structure using a dendrogram analysis of the NH3 emission. The distance between neighbour structures that are embedded within the same parent structure is generally greater than expected from a spherical Jeans analysis, and is in better agreement with cylindrical fragmentation models. The NH3 line width-size relation is flat, and average gas motions are sub- or trans-sonic over all physical scales observed. Subsonic regions extend far beyond the typical 0.1 pc scale previously identified in star-forming cores. As a result, we find a strong trend of decreasing virial parameter with increasing structure mass in Serpens South. Extremely low virial parameters on the largest scales probed by our data suggest that the previously observed, ordered magnetic field is insufficient to support the region against collapse, in agreement with large radial infall motions previously measured toward some of the filaments. A more complex magnetic field configuration in the dense gas, however, may be able to support the filaments.

Although surveys of infall motions in dense cores have been carried out for years, few surveys have focused on mapping infall across cores using multiple spectral line observations. To fill this gap, we present IRAM 30-m Telescope maps of N2H+(1-0), DCO+(2-1), DCO+(3-2), and HCO+(3-2) emission towards two prestellar cores (L492 and L694-2) and one protostellar core (L1521F). We find that the measured infall velocity varies with position across each core and choice of molecular line, likely as a result of radial variations in core chemistry and dynamics. Line-of-sight infall speeds estimated from DCO+(2-1) line profiles can decrease by 40-50 m/s when observing at a radial offset >= 0.04 pc from the core's dust continuum emission peak. Median infall speeds calculated from all observed positions across a core can also vary by as much as 65 m/s depending on the transition. These results show that while single-pointing, single-transition surveys of core infall velocities may be good indicators of whether a core is either contracting or expanding, the magnitude of the velocities they measure are significantly impacted by the choice of molecular line, proximity to the core center, and core evolutionary state.

The inner regions of the envelopes surrounding young protostars are characterised by a complex chemistry, with prebiotic molecules present on the scales where protoplanetary disks eventually may form. This paper introduces a systematic survey, "Protostellar Interferometric Line Survey (PILS)" of the Class 0 protostellar binary IRAS 16293-2422 using the Atacama Large Millimeter/submillimeter Array (ALMA). The survey covers the full frequency range from 329 to 363 GHz (0.8 mm) with additional targeted observations at 3.0 and 1.3 mm. More than 10,000 features are detected toward one component in the protostellar binary. Glycolaldehyde, its isomers, methyl formate and acetic acid, and its reduced alcohol, ethylene glycol, are clearly detected. For ethylene glycol both lowest state conformers, aGg' and gGg', are detected, the latter for the first time in the ISM. The abundance of glycolaldehyde is comparable to or slightly larger than that of ethylene glycol. In comparison to the Galactic Center, these two species are over-abundant relative to methanol, possibly an indication of formation at low temperatures in CO-rich ices. Both 13C and deuterated isotopologues of glycolaldehyde are detected, also for the first time ever in the ISM. For the deuterated species, a D/H ratio of approximately 5% is found with no differences between the deuteration in the different functional groups of glycolaldehyde. Measurements of the 13C-species lead to a 12C:13C ratio of approximately 30, lower than the typical ISM value. This low ratio may reflect an enhancement of 13CO in the ice due to either ion-molecule reactions in the gas before freeze-out or differences in the temperatures where 12CO and 13CO ices sublimate. The results reinforce the importance of low-temperature grain surface chemistry for the formation of prebiotic molecules seen here in the gas after sublimation of the entire ice mantle.

Formamide (NH$_2$CHO) has previously been detected in several star-forming regions and is thought to be a precursor for different prebiotic molecules. Its formation mechanism is still debated, however. Observations of formamide, related species and their isopotologues may provide useful clues to the chemical pathways leading to their formation. The Protostellar Interferometric Line Survey (PILS) represents an unbiased high angular resolution and sensitivity spectral survey of the low-mass protostellar binary IRAS 16293-2422 with the Atacama Large Millimeter/submillimeter Array (ALMA). We detect for the first time the three singly deuterated forms of NH$_2$CHO (NH$_2$CDO, cis- and trans-NHDCHO) as well as DNCO towards the component B of this binary source. The images reveal that the different isotopologues are all present in the same region. Based on observations of the $^{13}$C isotopologues of formamide and a standard $^{12}$C/$^{13}$C ratio, the deuterium fractionation is found to be similar for the three different forms with a value of about 2%. The DNCO/HNCO ratio is also comparable to the D/H ratio of formamide ($\sim$1%). These results are in agreement with the hypothesis that NH$_2$CHO and HNCO are chemically related through grain surface formation.

We present a catalog of low-mass dense cores observed with the SHARC-II instrument at 350 microns. Our observations have an effective angular resolution of 10", approximately 2.5 times higher than observations at the same wavelength obtained with the Herschel Space Observatory, albeit with lower sensitivity, especially to extended emission. The catalog includes 81 maps covering a total of 164 detected sources. For each detected source, we tabulate basic source properties including position, peak intensity, flux density in fixed apertures, and radius. We examine the uncertainties in the pointing model applied to all SHARC-II data and conservatively find that the model corrections are good to within ~3", approximately 1/3 of the SHARC-II beam. We examine the differences between two array scan modes and find that the instrument calibration, beam size, and beam shape are similar between the two modes. We also show that the same flux densities are measured when sources are observed in the two different modes, indicating that there are no systematic effects introduced into our catalog by utilizing two different scan patterns during the course of taking observations. We find a detection rate of 95% for protostellar cores but only 45% for starless cores, and demonstrate the existence of a SHARC-II detection bias against all but the most massive and compact starless cores. Finally, we discuss the improvements in protostellar classification enabled by these 350 micron observations.

We present an Atacama Large Millimeter/submillimeter Array (ALMA) 106 GHz (Band 3) continuum survey of the complete population of dense cores in the Chamaeleon I molecular cloud. We detect a total of 24 continuum sources in 19 different target fields. All previously known Class 0 and Class I protostars in Chamaeleon I are detected, whereas all of the 56 starless cores in our sample are undetected. We show that the Spitzer+Herschel census of protostars in Chamaeleon I is complete, with the rate at which protostellar cores have been misclassified as starless cores calculated as <1/56, or < 2%. We use synthetic observations to show that starless cores collapsing following the turbulent fragmentation scenario are detectable by our ALMA observations when their central densities exceed ~10^8 cm^-3, with the exact density dependent on the viewing geometry. Bonnor-Ebert spheres, on the other hand, remain undetected to central densities at least as high as 10^10 cm^-3. Our starless core non-detections are used to infer that either the star formation rate is declining in Chamaeleon I and most of the starless cores are not collapsing, matching the findings of previous studies, or that the evolution of starless cores are more accurately described by models that develop less substructure than predicted by the turbulent fragmentation scenario, such as Bonnor-Ebert spheres. We outline future work necessary to distinguish between these two possibilities.

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.

This paper discusses compelling science cases for a future long-baseline interferometer operating at millimeter and centimeter wavelengths, like the proposed Next Generation Vary Large Array (ngVLA). We report on the activities of the Cradle of Life science working group, which focused on the formation of low- and high-mass stars, the formation of planets and evolution of protoplanetary disks, the physical and compositional study of Solar System bodies, and the possible detection of radio signals from extraterrestrial civilizations. We propose 19 scientific projects based on the current specification of the ngVLA. Five of them are highlighted as possible Key Science Projects: (1) Resolving the density structure and dynamics of the youngest HII regions and high-mass protostellar jets, (2) Unveiling binary/multiple protostars at higher resolution, (3) Mapping planet formation regions in nearby disks on scales down to 1 AU, (4) Studying the formation of complex molecules, and (5) Deep atmospheric mapping of giant planets in the Solar System. For each of these projects, we discuss the scientific importance and feasibility. The results presented here should be considered as the beginning of a more in-depth analysis of the science enabled by such a facility, and are by no means complete or exhaustive.

There is much evidence that planet formation is occurring in the disk around the Herbig Be star HD100546. To learn more about the processes occurring in this disk we conducted high resolution imaging at 43/45 GHz with the Australia Telescope Compact Array (ATCA). Multiple array configurations were used, providing a best spatial resolution of $\sim$ 0.15 arcsec, or 15 AU at HD100546's distance of $\sim$ 100 pc. Significant structure is revealed, but its precise form is dependent on the $u-v$ plane sampling used for the image reconstruction. At a resolution of $\leq$ 30 AU we detected an inner gap in the disk with a radius of $\sim$ 25 AU and a position angle approximately along the known disk major axis. With different weighting, and an achieved resolution of $\sim$ 15 AU, emission appears at the centre and the disk takes on the shape of an incomplete ring, much like a horseshoe, again with a gap radius of $\sim$ 25 AU. The position angle of the disk major axis and its inclination from face-on are determined to be $140^{\circ}\pm5^{\circ}$ and $40^{\circ}\pm5^{\circ}$ respectively. The $\sim$ 25 AU gap radius is confirmed by a null in the real part of the binned visibilities at 320$\pm$10 k$\lambda$, whilst the non-axisymmetric nature is also confirmed through significant structure in the imaginary component. The emission mechanism at the central peak is most likely to be free-free emission from a stellar or disk wind. Overall our data support the picture of at least one, but probably several, giant planets orbiting HD100546 within 25 AU.

We present a CARMA 1.3 mm continuum survey toward 9 Class 0 protostars in the Perseus molecular cloud at $\sim$0.3$^{\prime\prime}$ (70 AU) resolution. This study approximately doubles the number of Class 0 protostars observed with spatial resolutions $<$ 100 AU at millimeter wavelengths, enabling the presence of protostellar disks and proto-binary systems to be probed. We detect flattened structures with radii $>$ 100 AU around 2 sources (L1448 IRS2 and Per-emb-14) and these sources may be strong disk candidates. Marginally-resolved structures with position angles within 30$^{\circ}$ of perpendicular to the outflow are found toward 3 protostars (L1448 IRS3C, IRAS 03282+3035, and L1448C) and are considered disk candidates. Two others (L1448 IRS3B and IRAS 03292+3039) have resolved structure, possibly indicative of massive inner envelopes or disks; L1448 IRS3B also has a companion separated by 0.9$^{\prime\prime}$ ($\sim$210 AU). IC348-MMS does not have well-resolved structure and the candidate first hydrostatic core L1451-MMS is marginally resolved on 1$^{\prime\prime}$ scales. The strong disk candidate sources were followed-up with C$^{18}$O ($J=2\rightarrow1$) observations, detecting velocity gradients consistent with rotation, but it is unclear if the rotation is Keplerian. We compare the observed visibility amplitudes to radiative transfer models, finding that visibility amplitude ratios suggest a compact component (possibly a disk) is necessary for 5 of 9 Class 0 sources; envelopes alone may explain the other 4 systems. We conclude that there is evidence for the formation of large disks in the Class 0 phase with a range of radii and masses dependent upon their initial formation conditions.

We present a multi-wavelength observational study of a low-mass star-forming region, L1251-C, with observational results at wavelengths from the near-infrared to the millimeter. Spitzer Space Telescope observations confirmed that IRAS 22343+7501 is a small group of protostellar objects. The extended emission to east-west direction with its intensity peak at the center of L1251A has been detected at 350 and 850 micron with the CSO and JCMT telescopes, tracing dense envelope materials around L1251A. The single-dish data from the KVN and TRAO telescopes show inconsistencies between the intensity peaks of several molecular line emission and that of the continuum emission, suggesting complex distributions of molecular abundances around L1251A. The SMA interferometer data, however, show intensity peaks of CO 2-1 and 13CO 2-1 located at the position of IRS 1, which is both the brightest source in IRAC image and the weakest source in the 1.3 mm dust continuum map. IRS 1 is the strongest candidate for the driving source of the newly detected compact CO 2-1 outflow. Over the whole region (14' by 14') of L125l-C, 3 Class I and 16 Class II sources have been detected, including three YSOs in L1251A. A comparison with the average projected distance among 19 YSOs in L1251-C and that among 3 YSOs in L1251A suggests L1251-C is an example of low-mass cluster formation, where protostellar objects are forming in a small group.

The measurement of Zeeman splitting in spectral lines---both in emission and absorption---can provide direct estimates of the magnetic field strength and direction in atomic and molecular clouds, both in our own Milky Way and in external galaxies. This method will probe the magnetic field in the warm and cold neutral components of the interstellar medium, providing a complement to the extensive SKA Faraday studies planning to probe the field in the ionized components.

Very Long Baseline Interferometry (VLBI) at radio wavelengths can provide astrometry accurate to 10 micro-arcseconds or better (i.e. better than the target GAIA accuracy) without being limited by dust obscuration. This means that unlike GAIA, VLBI can be applied to star-forming regions independently of their internal and line-of-sight extinction. Low-mass young stellar objects (particularly T Tauri stars) are often non-thermal compact radio emitters, ideal for astrometric VLBI radio continuum experiments. Existing observations for nearby regions (e.g. Taurus, Ophiuchus, or Orion) demonstrate that VLBI astrometry of such active T Tauri stars enables the reconstruction of both the regions' 3D structure (through parallax measurements) and their internal kinematics (through proper motions, combined with radial velocities). The extraordinary sensitivity of the SKA telescope will enable similar "tomographic mappings" to be extended to regions located several kpc from Earth, in particular to nearby spiral arm segments. This will have important implications for Galactic science, galactic dynamics and spiral structure theories.

We present Atacama Large Millimeter/submillimeter Array (ALMA) observations of the submillimeter dust continuum and H2D+ 1_{10}-1_{11} emission toward two evolved, potentially protostellar cores within the Ophiuchus molecular cloud, Oph A SM1 and SM1N. The data reveal small-scale condensations within both cores, with mass upper limits of M <~ 0.02M_Sun (~ 20 M_Jup). The SM1 condensation is consistent with a nearly-symmetric Gaussian source with a width of only 37 AU. The SM1N condensation is elongated, and extends 500 AU along its major axis. No evidence for substructure is seen in either source. A Jeans analysis indicates these sources are unlikely to fragment, suggesting that both will form single stars. H2D+ is only detected toward SM1N, offset from the continuum peak by ~150-200 AU. This offset may be due to either heating from an undetected, young, low luminosity protostellar source or first hydrostatic core, or HD (and consequently H2D+) depletion in the cold centre of the condensation. We propose that SM1 is protostellar, and that the condensation detected by ALMA is a warm (T ~ 30-50 K) accretion disk. The less concentrated emission of the SM1N condensation suggests that it is still starless, but we cannot rule out the presence of a low-luminosity source, perhaps surrounded by a pseudodisk. These data reveal observationally the earliest stages of the formation of circumstellar accretion regions, and agree with theoretical predictions that disk formation can occur very early in the star formation process, coeval with or just after the formation of a first hydrostatic core or protostar.

Context: The initial conditions for the gravitational collapse of molecular cloud cores and the subsequent birth of stars are still not well constrained. The characteristic cold temperatures (about 10 K) in such regions require observations at sub-millimetre and longer wavelengths. The Herschel Space Observatory and complementary ground-based observations presented in this paper have the unprecedented potential to reveal the structure and kinematics of a prototypical core region at the onset of stellar birth. Aims: This paper aims to determine the density, temperature, and velocity structure of the star-forming Bok globule CB 17. This isolated region is known to host (at least) two sources at different evolutionary stages: a dense core, SMM1, and a Class I protostar, IRS. Methods: We modeled the cold dust emission maps from 100 micron to 1.2 mm with both a modified blackbody technique to determine the optical depth-weighted line-of-sight temperature and column density and a ray-tracing technique to determine the core temperature and volume density structure. Furthermore, we analysed the kinematics of CB17 using the high-density gas tracer N2H+. Results: From the ray-tracing analysis, we find a temperature in the centre of SMM1 of 10.6 K, a flat density profile with radius 9500 au, and a central volume density of n(H) = 2.3x10^5 cm-3. The velocity structure of the N2H+ observations reveal global rotation with a velocity gradient of 4.3 km/s/pc. Superposed on this rotation signature we find a more complex velocity field, which may be indicative of differential motions within the dense core. Conclusions: SMM is a core in an early evolutionary stage at the verge of being bound, but the question of whether it is a starless or a protostellar core remains unanswered.

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

We present observations of the Auriga-California Molecular Cloud (AMC) at 3.6, 4.5, 5.8, 8.0, 24, 70 and 160 micron observed with the IRAC and MIPS detectors as part of the Spitzer Gould Belt Legacy Survey. The total mapped areas are 2.5 sq-deg with IRAC and 10.47 sq-deg with MIPS. This giant molecular cloud is one of two in the nearby Gould Belt of star-forming regions, the other being the Orion A Molecular Cloud (OMC). We compare source counts, colors and magnitudes in our observed region to a subset of the SWIRE data that was processed through our pipeline. Using color-magnitude and color-color diagrams, we find evidence for a substantial population of 166 young stellar objects (YSOs) in the cloud, many of which were previously unknown. Most of this population is concentrated around the LkHalpha 101 cluster and the filament extending from it. We present a quantitative description of the degree of clustering and discuss the fraction of YSOs in the region with disks relative to an estimate of the diskless YSO population. Although the AMC is similar in mass, size and distance to the OMC, it is forming about 15 - 20 times fewer stars.

We present observations of three Class 0/I protostars (L1157-mm, CB230 IRS1, and L1165-SMM1) using the Karl G. Jansky Very Large Array (VLA) and observations of two (L1165-SMM1 and CB230 IRS1) with the Combined Array for Research in Millimeter-wave Astronomy (CARMA). The VLA observations were taken at wavelengths of $\lambda = 7.3$ mm, 1.4 cm, 3.3 cm, 4.0 cm, and 6.5 cm with a best resolution of $\sim$0\farcs06 (18 AU) at 7.3 mm. The L1165-SMM1 CARMA observations were taken at $\lambda = 1.3$ mm with a best resolution of $\sim0\farcs3$ (100 AU), and the CB230 IRS1 observations were taken at $\lambda = 3.4$ mm with a best resolution of $\sim$3\arcsec\ (900 AU). We find that L1165-SMM1 and CB230 IRS1 have probable binary companions at separations of $\sim$0\farcs3 (100 AU) from detections of secondary peaks at multiple wavelengths. The position angles of these companions are nearly orthogonal to the direction of the observed bipolar outflows, consistent with the expected protostellar disk orientations. We suggest that these companions may have formed from disk fragmentation; turbulent fragmentation would not preferentially arrange the binary companions to be orthogonal to the outflow direction. For L1165-SMM1, both the 7.3 mm and 1.3 mm emission show evidence of a large (R $>$ 100 AU) disk. For the L1165-SMM1 primary protostar and the CB230 IRS1 secondary protostar, the 7.3 mm emission is resolved into structures consistent with $\sim20$ AU radius disks. For the other protostars, including L1157-mm, the emission is unresolved, suggesting disks with radii $< 20$ AU.

We have detected bright HC7N J = 21-20 emission toward multiple locations in the Serpens South cluster-forming region using the K-Band Focal Plane Array at the Robert C. Byrd Green Bank Telescope. HC7N is seen primarily toward cold filamentary structures that have yet to form stars, largely avoiding the dense gas associated with small protostellar groups and the main central cluster of Serpens South. Where detected, the HC7N abundances are similar to those found in other nearby star forming regions. Toward some HC7N `clumps', we find consistent variations in the line centroids relative to NH3 (1,1) emission, as well as systematic increases in the HC7N non-thermal line widths, which we argue reveal infall motions onto dense filaments within Serpens South with minimum mass accretion rates of M ~ 2-5 M_sun Myr^-1. The relative abundance of NH3 to HC7N suggests that the HC7N is tracing gas that has been at densities n ~ 10^4 cm^-3, for timescales t < 1-2 x 10^5 yr. Since HC7N emission peaks are rarely co-located with those of either NH3 or continuum, it is likely that Serpens South is not particularly remarkable in its abundance of HC7N, but instead the serendipitous mapping of HC7N simultaneously with NH3 has allowed us to detect HC7N at low abundances in regions where it otherwise may not have been looked for. This result extends the known star-forming regions containing significant HC7N emission from typically quiescent regions, like the Taurus molecular cloud, to more complex, active environments.

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

We report interferometric observations of the VeLLO L1521F-IRS, in 12CO (2-1) line emission and the 1.3 mm continuum emission, using the Submillimeter Array (SMA). With the 12CO high-resolution observations, we have spatially resolved a compact but poorly collimated molecular outflow associated with L1521F-IRS for the first time. The blueshifted and redshifted lobes are aligned along the east and west side of L1521F-IRS with a lobe size of ~1000 AU. The estimated outflow mass, maximum outflow velocity, and outflow force are (9.0-80)x10^{-4} M_{sol}, 7.2 km s^{-1}, and (7.4-66)x10^{-7} M_{sol} km s^{-1} yr^{-1}, respectively. The estimated outflow parameters such as size, mass, and momentum rate are similar to values derived for other VeLLOs, and are located at the lower end of values compared to previously studied outflows associated with low- to high-mass star forming regions. Low-velocity less collimated (1.5 km s^{-1}/1200 AU) and higher-velocity compact (4.0 km s^{-1}/920 AU) outflow components are suggested by the data. These velocity structures are not consistent with those expected in the jet driven or wind driven outflow models, perhaps suggesting a remnant outflow from the FHSC as well as an undeveloped outflow from the protostar. Detection of an infrared source and compact millimeter continuum emission suggest the presence of the protostar, while its low bolometric luminosity (0.034-0.07 L$_{\odot}$), and small outflow, suggests that L1521F is in the earliest protostellar stage ($<10^{4}$ yr) and contains a substellar mass object. The bolometric (or internal) luminosity of L1521F-IRS suggests that the current mass accretion rate is an order-of-magnitude lower than expected in the standard mass accretion model (~10^{-6} M_{sol} yr^{-1}), which may imply that L1521F-IRS is currently in a low activity phase.

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.

We present observations of the Zeeman effect in OH thermal absorption main lines at 1665 and 1667 MHz taken with the Very Large Array (VLA) toward the star forming region S88B. The OH absorption profiles toward this source are complicated, and contain several blended components toward a number of positions. Almost all of the OH absorbing gas is located in the eastern parts of S88B, toward the compact continuum source S88B-2 and the eastern parts of the extended continuum source S88B-1. The ratio of 1665/1667 MHz OH line intensities indicates the gas is likely highly clumped, in agreement with other molecular emission line observations in the literature. S88-B appears to present a similar geometry to the well-known star forming region M17, in that there is an edge-on eastward progression from ionized to molecular gas. The detected magnetic fields appear to mirror this eastward transition; we detected line-of-sight magnetic fields ranging from 90-400 \mu G, with the lowest values of the field to the southwest of the S88B-1 continuum peak, and the highest values to its northeast. We used the detected fields to assess the importance of the magnetic field in S88B by a number of methods; we calculated the ratio of thermal to magnetic pressures, we calculated the critical field necessary to completely support the cloud against self-gravity and compared it to the observed field, and we calculated the ratio of mass to magnetic flux in terms of the critical value of this parameter. All these methods indicated that the magnetic field in S88B is dynamically significant, and should provide an important source of support against gravity. Moreover, the magnetic energy density is in approximate equipartition with the turbulent energy density, again pointing to the importance of the magnetic field in this region.

Observations of the circumstellar disk in the Bok globule CB 26 at 110, 230, and 270 GHz are presented together with the results of the simulations and estimates of the disk parameters. These observations were obtained using the SMA, IRAM Plateau de Bure, and OVRO interferometers. The maps have relatively high angular resolutions (0.4-1"), making it possible to study the spatial structure of the gas-dust disk. The disk parameters are reconstructed via a quantitative comparison of observational and theoretical intensity maps. The disk model used to construct the theoretical maps is based on the assumption of hydrostatic and radiative equilibrium in the vertical direction, while the radial surface density profile is described phenomenologically. The system of equations for the transfer of the infrared and ultraviolet radiation is solved in the vertical direction, in order to compute the thermal structure of the disk. The disk best-fit parameters are derived for each map and all the maps simultaneously, using a conjugate gradient method. The degrees of degeneracy of the parameters describing the thermal structure and density distribution of the disk are analyzed in detail. All three maps indicate the presence of an inner dust-free region with a radius of approximately 35 AU, in agreement with the conclusions of other studies. The inclination of the disk is 78 deg, which is smaller than the value adopted in our earlier study of rotating molecular outflows from CB 26. The model does not provide any evidence for the growth of dust particles above a_max=0.02 cm.

One puzzle in understanding how stars form in clusters is the source of mass -- is all of the mass in place before the first stars are born, or is there an extended period when the cluster accretes material which can continuously fuel the star formation process? We use a multi-line spectral survey of the southern filament associated with the Serpens South embedded cluster-forming region in order to determine if mass is accreting from the filament onto the cluster, and whether the accretion rate is significant. Our analysis suggests that material is flowing along the filament's long axis at a rate of ~30Msol/Myr (inferred from the N2H+ velocity gradient along the filament), and radially contracting onto the filament at ~130Msol/Myr (inferred from HNC self-absorption). These accretion rates are sufficient to supply mass to the central cluster at a similar rate to the current star formation rate in the cluster. Filamentary accretion flows may therefore be very important in the ongoing evolution of this cluster.

We have mapped the Auriga/California molecular cloud with the Herschel PACS and SPIRE cameras and the Bolocam 1.1 mm camera on the Caltech Submillimeter Observatory (CSO) with the eventual goal of quantifying the star formation and cloud structure in this Giant Molecular Cloud (GMC) that is comparable in size and mass to the Orion GMC, but which appears to be forming far fewer stars. We have tabulated 60 compact 70/160um sources that are likely pre-main-sequence objects and correlated those with Spitzer and WISE mid-IR sources. At 1.1 mm we find 18 cold, compact sources and discuss their properties. The most important result from this part of our study is that we find a modest number of additional compact young objects beyond those identified at shorter wavelengths with Spitzer. We also describe the dust column density and temperature structure derived from our photometric maps. The column density peaks at a few x 10^22 cm^-2 (N_H2) and is distributed in a clear filamentary structure along which nearly all the pre-main-sequence objects are found. We compare the YSO surface density to the gas column density and find a strong non-linear correlation between them. The dust temperature in the densest parts of the filaments drops to ~10K from values ~ 14--15K in the low density parts of the cloud. We also derive the cumulative mass fraction and probability density function of material in the cloud which we compare with similar data on other star-forming clouds.

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.

IRAS 20050+2720 is young star forming region at a distance of 700 pc without apparent high mass stars. We present results of our multiwavelength study of IRAS 20050+2720 which includes observations by Chandra and Spitzer, and 2MASS and UBVRI photometry. In total, about 300 YSOs in different evolutionary stages are found. We characterize the distribution of young stellar objects (YSOs) in this region using a minimum spanning tree (MST) analysis. We newly identify a second cluster core, which consists mostly of class II objects, about 10 arcmin from the center of the cloud. YSOs of earlier evolutionary stages are more clustered than more evolved objects. The X-ray luminosity function (XLF) of IRAS 20050+2720 is roughly lognormal, but steeper than the XLF of the more massive Orion nebula complex. IRAS 20050+2720 shows a lower N_H/A_K ratio compared with the diffuse ISM.

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.

We present 230 GHz Submillimeter Array continuum and molecular line observations of the newly discovered FUor candidate HBC722. We report the detection of seven 1.3 mm continuum sources in the vicinity of HBC722, none of which correspond to HBC722 itself. We compile infrared and submillimeter continuum photometry of each source from previous studies and conclude that three are Class 0 embedded protostars, one is a Class I embedded protostar, one is a Class I/II transition object, and two are either starless cores or very young, very low luminosity protostars or first hydrostatic cores. We detect a northwest-southeast outflow, consistent with the previous detection of such an outflow in low-resolution, single-dish observations, and note that its axis may be precessing. We show that this outflow is centered on and driven by one of the nearby Class 0 sources rather than HBC722, and find no conclusive evidence that HBC722 itself is driving an outflow. The non-detection of HBC722 in the 1.3 mm continuum observations suggests an upper limit of 0.02 solar masses for the mass of the circumstellar disk. This limit is consistent with typical T Tauri disks and with a disk that provides sufficient mass to power the burst.

We present high angular resolution SMA and Spitzer observations toward the Bok globule CB17. SMA 1.3mm dust continuum images reveal within CB17 two sources with an angular separation of about 21" (about 5250 AU at a distance of 250 pc). The northwestern continuum source, referred to as CB17 IRS, dominates the infrared emission in the Spitzer images, drives a bipolar outflow extending in the northwest-southeast direction, and is classified as a low luminosity Class0/I transition object (L_bol ~ 0.5 L_sun). The southeastern continuum source, referred to as CB17 MMS, has faint dust continuum emission in the SMA 1.3mm observations (about 6 sigma detection; ~3.8 mJy), but is not detected in the deep Spitzer infrared images at wavelengths from 3.6 to 70 micron. Its bolometric luminosity and temperature, estimated from its spectral energy distribution, are less than 0.04 L_sun and 16 K, respectively. The SMA CO(2-1) observations suggest that CB17 MMS may drive a low-velocity molecular outflow (about 2.5 km/s), extending in the east-west direction. Comparisons with prestellar cores and Class0 protostars suggest that CB17 MMS is more evolved than prestellar cores but less evolved than Class0 protostars. The observed characteristics of CB17 MMS are consistent with the theoretical predictions from radiative/magneto hydrodynamical simulations of a first hydrostatic core, but there is also the possibility that CB17 MMS is an extremely low luminosity protostar deeply embedded in an edge-on circumstellar disk. Further observations are needed to study the properties of CB17 MMS and to address more precisely its evolutionary stage.

We present high spatial (<300 AU) and spectral (0.07 km/s) resolution Submillimeter Array observations of the dense starless cluster core Oph A-N6, in the 1 mm dust continuum and the 3-2 line of N2H+ and N2D+. The dust continuum observations reveal a compact source not seen in single-dish observations, of size ~1000 AU and mass 0.005-0.01 M\odot. The combined line and single-dish observations reveal a core of size 3000 \times 1400 AU elongated in a NW-SE direction, with almost no variation in either line width or line center velocity across the map, and very small non-thermal motions. The deuterium fraction has a peak value of ~0.15 and is >0.05 over much of the core. The N2H+ column density profile across the major axis of Oph A-N6 is well represented by an isothermal cylinder, with temperature 20 K, peak density 7.1 \times 10^6 cm^{-3}, and N2H+ abundance 2.7 \times 10^{-10}. The mass of Oph A-N6 is estimated to be 0.29 M\odot, compared to a value of 0.18 M\odot from the isothermal cylinder analysis, and 0.63 M\odot for the critical mass for fragmentation of an isothermal cylinder. Compared to isolated low-mass cores, Oph A-N6 shows similar narrow line widths and small velocity variation, with a deuterium fraction similar to "evolved" dense cores. It is significantly smaller than isolated cores, with larger peak column and volume density. The available evidence suggests Oph A-N6 has formed through the fragmentation of the Oph A filament and is the precursor to a low-mass star. The dust continuum emission suggests it may already have begun to form a star.

We report Spitzer observations of five newly identified bow shocks in the massive star-forming region RCW 38. Four are visible at IRAC wavelengths, the fifth is visible only at 24 microns. Chandra X-ray emission indicates that winds from the central O5.5 binary, IRS~2, have caused an outflow to the NE and SW of the central subcluster. The southern lobe of hot ionised gas is detected in X-rays; shocked gas and heated dust from the shock-front are detected with Spitzer at 4.5 and 24 microns. The northern outflow may have initiated the present generation of star formation, based on the filamentary distribution of the protostars in the central subcluster. Further, the bow-shock driving star, YSO 129, is photo-evaporating a pillar of gas and dust. No point sources are identified within this pillar at near- to mid-IR wavelengths. We also report on IRAC 3.6 & 5.8 micron observations of the cluster DBS2003-124, NE of RCW 38, where 33 candidate YSOs are identified. One star associated with the cluster drives a parsec-scale jet. Two candidate HH objects associated with the jet are visible at IRAC and MIPS wavelengths. The jet extends over a distance of ~3 pc. Assuming a velocity of 100 km/s for the jet material gives an age of about 30,000 years, indicating that the star (and cluster) are likely to be very young, with a similar or possibly younger age than RCW 38, and that star formation is ongoing in the extended RCW 38 region.

We present a study of the structure of the high mass star-forming region RCW~38 and the spatial distribution of its young stellar population. Spitzer IRAC photometry 3-8um are combined with 2MASS near-IR data to identify young stellar objects by IR-excess emission from their circumstellar material. Chandra X-ray data are used to identify class III pre-main sequence stars lacking circumstellar material. We identify 624 YSOs: 23 class 0/I and 90 flat spectrum protostars, 437 Class II stars, and 74 Class III stars. We also identify 29 (27 new) O star candidates over the IRAC field. Seventy-two stars exhibit IR-variability, including seven class 0/I and 12 flat spectrum YSOs. A further 177 tentative candidates are identified by their location in the IRAC [3.6] vs. [3.6]-[5.8] cmd. We find strong evidence of subclustering in the region. Three subclusters were identified surrounding the central cluster, with massive and variable stars in each subcluster. The central region shows evidence of distinct spatial distributions of the protostars and pre-main sequence stars. A previously detected IR cluster, DB2001_Obj36, has been established as a subcluster of RCW 38. This suggests that star formation in RCW 38 occurs over a more extended area than previously thought. The gas to dust ratio is examined using the X-ray derived hydrogen column density, N_H and the K-band extinction, and found to be consistent with the diffuse ISM, in contrast with Serpens & NGC1333. We posit that the high photoionising flux of massive stars in RCW 38 affects the agglomeration of the dust grains.

We present the detection of a dust continuum source at 3-mm (CARMA) and 1.3-mm (SMA), and 12CO(2-1) emission (SMA) towards the L1451-mm dense core. These detections suggest a compact object and an outflow where no point source at mid-infrared wavelengths is detected using Spitzer. An upper limit for the dense core bolometric luminosity of 0.05 Lsun is obtained. By modeling the broadband SED and the continuum interferometric visibilities simultaneously, we confirm that a central source of heating is needed to explain the observations. This modeling also shows that the data can be well fitted by a dense core with a YSO and disk, or by a dense core with a central First Hydrostatic Core (FHSC). Unfortunately, we are not able to decide between these two models, which produce similar fits. We also detect 12CO(2-1) emission with red- and blue-shifted emission suggesting the presence of a slow and poorly collimated outflow, in opposition to what is usually found towards young stellar objects but in agreement with prediction from simulations of a FHSC. This presents the best candidate, so far, for a FHSC, an object that has been identified in simulations of collapsing dense cores. Whatever the true nature of the central object in L1451-mm, this core presents an excellent laboratory to study the earliest phases of low-mass star formation.

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 new 230 GHz Submillimeter Array observations of the candidate first hydrostatic core Per-Bolo 58. We report the detection of a 1.3 mm continuum source and a bipolar molecular outflow, both centered on the position of the candidate first hydrostatic core. The continuum detection has a total flux density of 26.6 +/- 4.0 mJy, from which we calculate a total (gas and dust) mass of 0.11 +/- 0.05 Msun and a mean number density of 2.0 +/- 1.6 X 10^7 cm-3. There is some evidence for the existence of an unresolved component in the continuum detection, but longer-baseline observations are required in order to confirm the presence of this component and determine whether its origin lies in a circumstellar disk or in the dense inner envelope. The bipolar molecular outflow is observed along a nearly due east-west axis. The outflow is slow (characteristic velocity of 2.9 km/s), shows a jet-like morphology (opening semi-angles ~8 degrees for both lobes), and extends to the edges of the primary beam. We calculate the kinematic and dynamic properties of the outflow in the standard manner and compare them to several other protostars and candidate first hydrostatic cores with similarly low luminosities. We discuss the evidence both in support of and against the possibility that Per-Bolo 58 is a first hydrostatic core, and we outline future work needed to further evaluate the evolutionary status of this object.

We report the detection of a compact (of order 5 arcsec; about 1800 AU projected size) CO outflow from L1148-IRS. This confirms that this Spitzer source is physically associated with the nearby (about 325 pc) L1148 dense core. Radiative transfer modeling suggests an internal luminosity of 0.08 to 0.13 L_sun. This validates L1148-IRS as a Very Low Luminosity Object (VeLLO; L < 0.1 L_sun). The L1148 dense core has unusually low densities and column densities for a star-forming core. It is difficult to understand how L1148-IRS might have formed under these conditions. Independent of the exact final mass of this VeLLO (which is likely < 0.24 M_sun), L1148-IRS and similar VeLLOs might hold some clues about the isolated formation of brown dwarfs.

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.

We present Gould's Belt (GB) Spitzer IRAC and MIPS observations of the Lupus V and VI clouds and discuss them in combination with near-infrared (2MASS) data. Our observations complement those obtained for other Lupus clouds within the frame of the Spitzer "Core to Disk" (c2d) Legacy Survey. We found 43 Young Stellar Object (YSO) candidates in Lupus V and 45 in Lupus VI, including 2 transition disks, using the standard c2d/GB selection method. None of these sources was classified as a pre-main sequence star from previous optical, near-IR and X-ray surveys. A large majority of these YSO candidates appear to be surrounded by thin disks (Class III; ~79% in Lupus V and ~87% in Lupus VI). These Class III abundances differ significantly from those observed for the other Lupus clouds and c2d/GB surveyed star-forming regions, where objects with optically thick disks (Class II) dominate the young population. We investigate various scenarios that can explain this discrepancy. In particular, we show that disk photo-evaporation due to nearby OB stars is not responsible for the high fraction of Class III objects. The gas surface densities measured for Lupus V and VI lies below the star-formation threshold (AV {\eqsim}8.6 mag), while this is not the case for other Lupus clouds. Thus, few Myrs older age for the YSOs in Lupus V and VI with respect to other Lupus clouds is the most likely explanation of the high fraction of Class III objects in these clouds, while a higher characteristic stellar mass might be a contributing factor. Better constraints on the age and binary fraction of

We present new observations of the CB130 region, composed of three separate cores. Using the \textit{Spitzer Space Telescope} we detected a Class 0 and a Class II object in one of these, CB130-1. The observed photometric data from \textit{Spitzer} and ground-based telescopes are used to establish the physical parameters of the Class 0 object. SED fitting with a radiative transfer model shows that the luminosity of the Class 0 object is 0.14 $-$ 0.16 L$_{\odot}$, which is a low luminosity for a protostellar object. In order to constrain the chemical characteristics of the core having the low luminosity object, we compare our molecular line observations to models of lines including abundance variations. We tested both ad hoc step function abundance models and a series of self-consistent chemical evolution models. In the chemical evolution models, we consider a continuous accretion model and an episodic accretion model to explore how variable luminosity affects the chemistry. The step function abundance models can match observed lines reasonably well. The best fitting chemical evolution model requires episodic accretion and the formation of CO$_2$ ice from CO ice during the low luminosity periods. This process removes C from the gas phase, providing a much improved fit to the observed gas-phase molecular lines and the CO$_2$ ice absorption feature. Based on the chemical model result, the low luminosity of CB130-1 is explained better as a quiescent stage between episodic accretion bursts rather than being at the first hydrostatic core stage.

Class I methanol masers are believed to be produced in the shock-excited environment around star-forming regions. Many authors have argued that the appearance of various subsets of class I masers may be indicative of specific evolutionary stages of star formation or excitation conditions. Until recently, however, no major interferometer was capable of imaging the important 36 GHz transition. We report on Expanded Very Large Array observations of the 36 GHz methanol masers and Submillimeter Array observations of the 229 GHz methanol masers in DR21(OH), DR21N, and DR21W. The distribution of 36 GHz masers in the outflow of DR21(OH) is similar to that of the other class I methanol transitions, with numerous multitransition spatial overlaps. At the site of the main continuum source in DR21(OH), class I masers at 36 and 229 GHz are found in virtual overlap with class II 6.7 GHz masers. To the south of the outflow, the 36 GHz masers are scattered over a large region but usually do not appear coincident with 44 GHz masers. In DR21W we detect an "S-curve" signature in Stokes V that implies a large value of the magnetic field strength if interpreted as due to Zeeman splitting, suggesting either that class I masers may exist at higher densities than previously believed or that the direct Zeeman interpretation of S-curve Stokes V profiles in class I masers may be incorrect. We find a diverse variety of different maser phenomena in these sources, suggestive of differing physical conditions among them.

We present early results of our multiwavelength study of the star forming region IRAS 20050+2720. While we use X-rays and IR to classify young stellar objects, the optical data can be used to exclude foreground objects. The dataset set contains 57 class I sources, 183 class II sources and 183 X-ray sources, of which 140 are class III candidates. Within IRAS 20050+2720 four subclusters are found. Subcluster A and B are the densest regions, which contain stars of all evolutionary stages. Subcluster C is much younger than the other two. It has not formed any class III objects yet. We newly identify a fourth subcluster, which consists mostly of class II objects and is located about 10 arcmin from the center of the cloud.

We present a global study of low mass, young stellar object (YSO) surface densities in nearby (< 500 pc) star forming regions based on a comprehensive collection of Spitzer Space Telescope surveys. We show that the distribution of YSO surface densities in the solar neighbourhood is a smooth distribution, being adequately described by a lognormal function from a few to 10^3 YSOs per pc^2, with a peak at 22 stars/pc^2 and a dispersion of 0.85. We do not find evidence for multiple discrete modes of star-formation (e.g. clustered and distributed). Comparing the observed surface density distribution to previously reported surface density threshold definitions of clusters, we find that the fraction of stars in clusters is crucially dependent on the adopted definitions, ranging from 40 to 90%. However, we find that only a low fraction (< 26%) of stars are formed in dense environments where their formation/evolution (along with their circumstellar disks and/or planets) may be affected by the close proximity of their low-mass neighbours.

We present new infrared, submillimeter, and millimeter observations of the dense core L673-7 and report the discovery of a low-luminosity, embedded Class 0 protostar driving a molecular outflow. L673-7 is seen in absorption against the mid-infrared background in 5.8, 8, and 24 micron Spitzer images, allowing for a derivation of the column density profile and total enclosed mass of L673-7, independent of dust temperature assumptions. Estimates of the core mass from these absorption profiles range from 0.2-4.5 solar masses. Millimeter continuum emission indicates a mass of about 2 solar masses, both from a direct calculation assuming isothermal dust and from dust radiative transfer models constrained by the millimeter observations. We use dust radiative transfer models to constrain the internal luminosity of L673-7, defined to be the luminosity of the central source and excluding the luminosity from external heating, to be 0.01-0.045 solar luminosities, with 0.04 solar luminosities the most likely value. L673-7 is thus classified as a very low luminosity object (VeLLO), and is among the lowest luminosity VeLLOs yet studied. We calculate the kinematic and dynamic properties of the molecular outflow in the standard manner, and we show that the expected accretion luminosity based on these outflow properties is greater than or equal to 0.36 solar luminosities. The discrepancy between this expected accretion luminosity and the internal luminosity derived from dust radiative transfer models indicates that the current accretion rate is much lower than the average rate over the lifetime of the outflow. Although the protostar embedded within L673-7 is consistent with currently being substellar, it is unlikely to remain as such given the substantial mass reservoir remaining in the core.

We present N2D+ 3-2 (IRAM) and H2D+ 1_11 - 1_10 and N2H+ 4-3 (JCMT) maps of the small cluster-forming Ophiuchus B2 core in the nearby Ophiuchus molecular cloud. In conjunction with previously published N2H+ 1-0 observations, the N2D+ data reveal the deuterium fractionation in the high density gas across Oph B2. The average deuterium fractionation R_D = N(N2D+)/N(N2H+) ~ 0.03 over Oph B2, with several small scale R_D peaks and a maximum R_D = 0.1. The mean R_D is consistent with previous results in isolated starless and protostellar cores. The column density distributions of both H2D+ and N2D+ show no correlation with total H2 column density. We find, however, an anticorrelation in deuterium fractionation with proximity to the embedded protostars in Oph B2 to distances >= 0.04 pc. Destruction mechanisms for deuterated molecules require gas temperatures greater than those previously determined through NH3 observations of Oph B2 to proceed. We present temperatures calculated for the dense core gas through the equating of non-thermal line widths for molecules (i.e., N2D+ and H2D+) expected to trace the same core regions, but the observed complex line structures in B2 preclude finding a reasonable result in many locations. This method may, however, work well in isolated cores with less complicated velocity structures. Finally, we use R_D and the H2D+ column density across Oph B2 to set a lower limit on the ionization fraction across the core, finding a mean x_e, lim >= few x 10^{-8}. Our results show that care must be taken when using deuterated species as a probe of the physical conditions of dense gas in star-forming regions.

The protostellar jet driven by L1448C was observed in the SiO J=8-7 and CO J=3-2 lines and 350 GHz dust continuum at ~1" resolution with the Submillimeter Array (SMA). A narrow jet from the northern source L1448C(N) was observed in the SiO and the high-velocity CO. The jet consists of a chain of emission knots with an inter-knot spacing of ~2" (500 AU) and a semi-periodic velocity variation. The innermost pair of knots, which are significant in the SiO map but barely seen in the CO, are located at ~1" (250 AU) from the central source, L1448C(N). Since the dynamical time scale for the innermost pair is only ~10 yr, SiO may have been formed in the protostellar wind through the gas-phase reaction, or been formed on the dust grain and directly released into the gas phase by means of shocks. It is found that the jet is extremely active with a mechanical luminosity of ~7 L_sun, which is comparable to the bolometric luminosity of the central source (7.5 L_sun). The mass accretion rate onto the protostar derived from the mass-loss rate is ~10^{-5} M_sun/yr. Such a high mass accretion rate suggests that the mass and the age of the central star are 0.03-0.09 M_sun and (4-12)x10^3 yr, respectively, implying that the central star is in the very early stage of protostellar evolution. The low-velocity CO emission delineates two V-shaped shells with a common apex at L1448C(N). The kinematics of these shells are reproduced by the model of a wide opening angle wind. The co-existence of the highly-collimated jets and the wide-opening angle shells can be explained by the unified X-wind model" in which highly-collimated jet components correspond to the on-axis density enhancement of the wide-opening angle wind. The CO $J$=3--2 map also revealed the second outflow driven by the southern source L1448C(S) located at ~8.3" (2000 AU) from L1448C(N).

Intermediate between the prestellar and Class 0 protostellar phases, the first core is a quasi-equilibrium hydrostatic object with a short lifetime and an extremely low luminosity. Recent MHD simulations suggest that the first core can even drive a molecular outflow before the formation of the second core (i.e., protostar). Using the Submillimeter Array and the Spitzer Space Telescope, we present high angular resolution observations towards the embedded dense core IRS2E in L1448. We find that source L1448 IRS2E is not visible in the sensitive Spitzer infrared images (at wavelengths from 3.6 to 70 um), and has weak (sub-)millimeter dust continuum emission. Consequently, this source has an extremely low bolometric luminosity (< 0.1 L_sun). Infrared and (sub-)millimeter observations clearly show an outflow emanating from this source; L1448 IRS2E represents thus far the lowest luminosity source known to be driving a molecular outflow. Comparisons with prestellar cores and Class 0 protostars suggest that L1448 IRS2E is more evolved than prestellar cores but less evolved than Class 0 protostars, i.e., at a stage intermediate between prestellar cores and Class 0 protostars. All these results are consistent with the theoretical predictions of the radiative/magneto hydrodynamical simulations, making L1448 IRS2E the most promising candidate of the first hydrostatic core revealed so far.

We present the results of a comprehensive infrared, submillimetre, and millimetre continuum emission study of isolated low-mass star-forming cores in 32 Bok globules, with the aim to investigate the process of star formation in these regions. The submillimetre and millimetre dust continuum emission maps together with the spectral energy distributions are used to model and derive the physical properties of the star-forming cores, such as luminosities, sizes, masses, densities, etc. Comparisons with ground-based near-infrared and space-based mid and far-infrared images from Spitzer are used to reveal the stellar content of the Bok globules, association of embedded young stellar objects with the submm dust cores, and the evolutionary stages of the individual sources. Submm dust continuum emission was detected in 26 out of the 32 globule cores observed. For 18 globules with detected (sub)mm cores we derive evolutionary stages and physical parameters of the embedded sources. We identify nine starless cores, most of which are presumably prestellar, nine Class 0 protostars, and twelve Class I YSOs. Specific source properties like bolometric temperature, core size, and central densities are discussed as function of evolutionary stage. We find that at least two thirds (16 out of 24) of the star-forming globules studied here show evidence of forming multiple stars on scales between 1,000 and 50,000 AU. However, we also find that most of these small prototstar and star groups are comprised of sources with different evolutionary stages, suggesting a picture of slow and sequential star formation in isolated globules

Young stars are formed within dusty discs. The grains in the disc are originally of the same size as interstellar dust. Models predict that these grains will grow in size through coagulation. Observations of the silicate features at micron wavelengths are consistent with growth to micron sizes whereas the slope of the SED at longer wavelengths traces growth up to mm sizes. We here look for a correlation between these two grain growth indicators. A large sample of T-Tauri and Herbig-Ae/Be stars was observed with the Spitzer Space Telescope at 5-13 micron; a subsample was observed at mm wavelengths. We complement this subsample with data from the literature to maximise the overlap between micron and mm observations and search for correlations. Synthetic spectra are produced to determine which processes may produce the dust evolution. Dust disc masses in the range <1 to 7 x 10^-4 MSun are obtained. Most sources have a mm spectral slope consistent with grain growth. There is a tentative correlation between the 10-micron silicate feature and the mm slope of the SED. The observed sources seem to be grouped per star-forming region in the micron-vs-mm diagram. The modelling results show that the 10-micron feature becomes flatter and subsequently the mm slope becomes shallower. Grain size distributions shallower than that of the ISM and/or bright central stars are required to explain specific features. Settling of larger grains towards the disc midplane affects the 10-micron feature, but hardly the mm slope. The tentative correlation between the strength of the 10-micron feature and the mm slope suggests that the inner and outer disc evolve simultaneously. Dust with a mass dominated by mm-sized grains is required to explain the shallowest mm slopes. Other processes besides grain growth may also be responsible for the removal of small grains.

A long infrared jet has been discovered by the Spitzer c2d legacy program in core A of L1251. It is associated with a very embedded Class 0 object with an accretion luminosity of about 0.9 Lsun derived by radiative transfer model fitting to the observed SED. Comparing the observed IRAC colors along the infrared jet with those calculated from a model of an admixture of gas with a power-law temperature distribution indicates that the jet is possibly created by a paraboloidal bow shock propagating into the ambient medium of n(H_2)=10^5 cm^{-3}. In addition, the variation of the power-law index along the jet suggests that the portion of hot gas decreases with distance from the jet engine. The molecular outflow in this region has been mapped for the first time using CO data. From the calculated outflow momentum flux, a very strong lower limit to the average accretion luminosity is 3.6 sin i/cos^3 i Lsun, indicative of a decrease in the accretion rate with time.

Fully sampled degree-scale maps of the 13CO 2-1 and CO 4-3 transitions toward three members of the Lupus Molecular Cloud Complex - Lupus I, III, and IV - trace the column density and temperature of the molecular gas. Comparison with IR extinction maps from the c2d project requires most of the gas to have a temperature of 8-10 K. Estimates of the cloud mass from 13CO emission are roughly consistent with most previous estimates, while the line widths are higher, around 2 km/s. CO 4-3 emission is found throughout Lupus I, indicating widespread dense gas, and toward Lupus III and IV. Enhanced line widths at the NW end and along the edge of the B228 ridge in Lupus I, and a coherent velocity gradient across the ridge, are consistent with interaction between the molecular cloud and an expanding HI shell from the Upper-Scorpius subgroup of the Sco-Cen OB Association. Lupus III is dominated by the effects of two HAe/Be stars, and shows no sign of external influence. Slightly warmer gas around the core of Lupus IV and a low line width suggest heating by the Upper-Centaurus-Lupus subgroup of Sco-Cen, without the effects of an HI shell.

We present results from the Spitzer Space Telescope and molecular line observations of 9 species toward the dark cloud L43. The Spitzer images and molecular line maps suggest it has a starless core and a Class I protostar evolving in the same environment. CO depletion is seen in both sources, and DCO+ lines are stronger toward the starless core. With a goal of testing the chemical characteristics from pre- to protostellar stages, we adopt an evolutionary chemical model to calculate the molecular abundances and compare with our observations. Among the different model parameters we tested, the best-fit model suggests a longer total timescale at the pre-protostellar stage, but with faster evolution at the later steps with higher densities.

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.

Circumstellar discs are expected to be the nursery of planets. Grain growth within such discs is the first step in the planet formation process. The Bok globule CB 26 harbours such a young disc. We present a detailed model of the edge-on circumstellar disc and its envelope in the Bok globule CB 26. The model is based on HST near-infrared maps in the I, J, H, and K bands, OVRO and SMA radio maps at 1.1mm, 1.3mm and 2.7mm, and the spectral energy distribution (SED) from 0.9 microns to 3mm. New photometric and spectroscopic data from the Spitzer Space Telescope and the Caltech Submilimeter Observatory have been obtained and are part of our analysis. Using the self-consistent radiative transfer code MC3D, the model we construct is able to discriminate parameter sets and dust properties of both its parts, namely envelope and disc. We find that the disc has an inner hole with a radius of 45 +/- 5 AU. Based on a dust model including silicate and graphite the maximum grain size needed to reproduce the spectral millimetre index is 2.5 microns. Features seen in the near-infrared images, dominated by scattered light, can be described as a result of a rotating envelope. Successful employment of ISM dust in both the disc and envelope hint that grain growth may not yet play a significant role for the appearance of this system. A larger inner hole gives rise to the assumption that CB 26 is a circumbinary disc.

Observations of the high-mass star forming region AFGL 2591 reveal a large abundance of CO+, a molecule known to be enhanced by far UV (FUV) and X-ray irradiation. In chemical models assuming a spherically symmetric envelope, the volume of gas irradiated by protostellar FUV radiation is very small due to the high extinction by dust. The abundance of CO+ is thus underpredicted by orders of magnitude. In a more realistic model, FUV photons can escape through an outflow region and irradiate gas at the border to the envelope. Thus, we introduce the first 2D axi-symmetric chemical model of the envelope of a high-mass star forming region to explain the CO+ observations as a prototypical FUV tracer. The model assumes an axi-symmetric power-law density structure with a cavity due to the outflow. The local FUV flux is calculated by a Monte Carlo radiative transfer code taking scattering on dust into account. A grid of precalculated chemical abundances, introduced in the first part of this series of papers, is used to quickly interpolate chemical abundances. This approach allows to calculate the temperature structure of the FUV heated outflow walls self-consistently with the chemistry. Synthetic maps of the line flux are calculated using a raytracer code. Single-dish and interferometric observations are simulated and the model results are compared to published and new JCMT and SMA observations. The two-dimensional model of AFGL 2591 is able to reproduce the JCMT single-dish observations and also explains the non-detection by the SMA. We conclude that the observed CO+ line flux and its narrow width can be interpreted by emission from the warm and dense outflow walls irradiated by protostellar FUV radiation.

HH 211 is a nearby young protostellar system with a highly collimated jet. We have mapped it in 352 GHz continuum, SiO (J=8-7), and HCO+ (J=4-3) emission at up to ~ 0.2" resolution with the Submillimeter Array (SMA). The continuum source is now resolved into two sources, SMM1 and SMM2, with a separation of ~ 84 AU. SMM1 is seen at the center of the jet, probably tracing a (inner) dusty disk around the protostar driving the jet. SMM2 is seen to the southwest of SMM1 and may trace an envelope-disk around a small binary companion. A flattened envelope-disk is seen in HCO+ around SMM1 with a radius of ~ 80 AU perpendicular to the jet axis. Its velocity structure is consistent with a rotation motion and can be fitted with a Keplerian law that yields a mass of ~ 50+-15 Jupiter mass (a mass of a brown dwarf) for the protostar. Thus, the protostar could be the lowest mass source known to have a collimated jet and a rotating flattened envelope-disk. A small-scale (~ 200 AU) low-speed (~ 2 km/s) outflow is seen in HCOP+ around the jet axis extending from the envelope-disk. It seems to rotate in the same direction as the envelope-disk and may carry away part of the angular momentum from the envelope-disk. The jet is seen in SiO close to ~ 100 AU from SMM1. It is seen with a "C-shaped" bending. It has a transverse width of <~ 40 AU and a velocity of ~ 170+-60 km/s. A possible velocity gradient is seen consistently across its innermost pair of knots, with ~ 0.5 km/s at ~ 10 AU, consistent with the sense of rotation of the envelope-disk. If this gradient is an upper limit of the true rotational gradient of the jet, then the jet carries away a very small amount of angular momentum of ~ 5 AU km/s and thus must be launched from the very inner edge of the disk near the corotation radius.

We present adaptive optics (AO) near-infrared (JHKs) observations of the deeply embedded massive cluster RCW 38 using NACO on the VLT. Narrowband AO observations centered at wavelengths of 1.28, 2.12, and 2.17 micron were also obtained. The area covered by these observations is about 0.5 pc square, centered on the O star RCW 38 IRS2. We use the JHKs colors to identify young stars with infrared excess. Through a detailed comparison to a nearby control field, we find that most of the 337 stars detected in all three infrared bands are cluster members (~317), with essentially no contamination due to background or foreground sources. Five sources have colors suggestive of deeply embedded protostars, while 53 sources are detected at Ks only; their spatial distribution with respect to the extinction suggests they are highly reddened cluster members. Detectable Ks-band excess is found toward 29 +/- 3 % of the stars. For comparison to a similar area of Orion observed in the near-infrared, mass and extinction cuts are applied, and the excess fractions redetermined. The resulting excesses are then 25 +/- 5 % for RCW 38, and 42 +/- 8 % for Orion. RCW 38 IRS2 is shown to be a massive star binary with a projected separation of ~500 AU. Two regions of molecular hydrogen emission are revealed through the 2.12 micron imaging. One shows a morphology suggestive of a protostellar jet, and is clearly associated with a star only detected at H and Ks, previously identified as a highly obscured X-ray source. Three spatially extended cometary-like objects, suggestive of photoevaporating disks, are identified, but only one is clearly directly influenced by RCW 38 IRS2. A King profile provides a reasonable fit to the cluster radial density profile and a nearest neighbor distance analysis shows essentially no sub-clustering.

We present Spitzer IRAC (~2 deg^2) and MIPS (~8 deg^2) observations of the Cepheus Flare which is associated with the Gould Belt, at an approximate distance of ~300 pc. Around 6500 sources are detected in all four IRAC bands, of which ~900 have MIPS 24 micron detections. We identify 133 YSO candidates using color-magnitude diagram techniques, a large number of the YSO candidates are associated with the NGC 7023 reflection nebula. Cross identifications were made with the Guide Star Catalog II and the IRAS Faint Source Catalog, and spectral energy distributions (SED) were constructed. SED modeling was conducted to estimate the degree of infrared excess. It was found that a large majority of disks were optically thick accreting disks, suggesting that there has been little disk evolution in these sources. Nearest-neighbor clustering analysis identified four small protostellar groups (L1228, L1228N, L1251A, and L1251B) with 5-15 members each and the larger NGC 7023 association with 32 YSO members. The star formation efficiency for cores with clusters of protostars and for those without clusters was found to be ~8% and ~1% respectively. The cores L1155, L1241, and L1247 are confirmed to be starless down to our luminosity limit of L_bol=0.06 L_sol.

Most stars are born in rich young stellar clusters (YSCs) embedded in giant molecular clouds. The most massive stars live out their short lives there, profoundly influencing their natal environments by ionizing HII regions, inflating wind-blown bubbles, and soon exploding as supernovae. Thousands of lower-mass pre-main sequence stars accompany the massive stars, and the expanding HII regions paradoxically trigger new star formation as they destroy their natal clouds. While this schematic picture is established, our understanding of the complex astrophysical processes involved in clustered star formation have only just begun to be elucidated. The technologies are challenging, requiring both high spatial resolution and wide fields at wavelengths that penetrate obscuring molecular material and remove contaminating Galactic field stars. We outline several important projects for the coming decade: the IMFs and structures of YSCs; triggered star formation around YSC; the fate of OB winds; the stellar populations of Infrared Dark Clouds; the most massive star clusters in the Galaxy; tracing star formation throughout the Galactic Disk; the Galactic Center region and YSCs in the Magellanic Clouds. Programmatic recommendations include: developing a 30m-class adaptive optics infrared telescope; support for high-resolution and wide field X-ray telescopes; large-aperture sub-millimeter and far-infrared telescopes; multi-object infrared spectrographs; and both numerical and analytical theory.

To better understand the initial conditions of the high-mass star formation process, it is crucial to study at high-angular resolution the morphology, the kinematics, and eventually the interactions of the coldest condensations associated with intermediate-/high-mass star forming regions. The paper studies the cold condensations in the intermediate-/high-mass proto-cluster IRAS 05345+3157, focusing the attention on the interaction with the other objects in the cluster. We have performed millimeter high-angular resolution observations, both in the continuum and several molecular lines, with the PdBI and the SMA. In a recent paper, we have already published part of these data. The main finding of that work was the detection of two cold and dense gaseous condensations, called N and S (masses ~2 and ~9 M_sun), characterised by high values of the deuterium fractionation (~0.1 in both cores). In this paper, we present a full report of the observations, and a complete analysis of the data obtained. The millimeter maps reveal the presence of 3 cores inside the interferometers primary beam, called C1-a, C1-b and C2. None of them are associated with cores N and S. C1-b is very likely associated with a newly formed early-B ZAMS star embedded inside a hot-core, while C1-a is more likely associated with a class 0 intermediate-mass protostar. The nature of C2 is unclear. Both C1-a and C1-b are good candidates as driving sources of a powerful CO outflow, which strongly interacts with N and S, as demonstrated by the velocity gradient across both condensations. Our major conclusion is that the chemical properties of these pre-stellar cores are similar to those observed in low-mass isolated ones, while the kinematics is dominated by the turbulence triggered by the CO outflow and can influece their evolution.

Grains in disks around young stars grow from interstellar submicron sizes to planetesimals over the course of several Myr. Thermal emission of large grains or pebbles can be best observed at cm wavelengths. However, other emission mechanisms can contribute. We aim to determine the mechanisms of cm emission for 3 T Tauri stars. WW Cha and RU Lup were recently found to have grain growth at least up to mm sizes in their circumstellar disks. CS Cha has similar indications for grain growth in its circumbinary disk. The T Tauri stars WW Cha and RU Lup were monitored over several years at mm and cm wavelengths, using ATCA. The new ATCA 7 mm system was also used to observe CS Cha. WW Cha was detected on several occasions at 7 and 16 mm. We obtained one detection of WW Cha at 3 cm and upper limits only at 6 cm. The emission at 16 mm was stable over days, months and years, but the emission at 3 cm is found to be variable. RU Lup was detected at 7 mm. It was observed at 16 mm 3 times and at 3 and 6 cm 4 times and found to be variable in all 3 wavebands. CS Cha was detected at 7 mm, but the S/N was too low to resolve the gap in the circumbinary disk. The emission at 3, 7 and 16 mm for WW Cha is well explained by thermal emission from mm and cm-sized pebbles. The cm spectral index is consistent with the emission from an optically-thick ionised wind, but the high variability of the cm emission points to a non-thermal contribution. The SEDs of RU Lup and CS Cha from 1 to 7 mm are consistent with thermal emission from mm-sized grains. The variability of the longer-wavelength emission for RU Lup and the negative spectral index suggest non-thermal emission.

This paper reports the discovery of a Very Low Luminosity Object (VeLLO) in the "starless" dense core L328, using the Spitzer Space Telescope and ground based observations from near-infrared to millimeter wavelengths. The Spitzer 8 micron image indicates that L328 consists of three subcores of which the smallest one may harbor a source, L328-IRS while two other subcores remain starless. L328-IRS is a Class 0 protostar according to its bolometric temperature (44 K) and the high fraction ~72 % of its luminosity emitted at sub-millimeter wavelengths. Its inferred "internal luminosity" (0.04 - 0.06 Lsun) using a radiative transfer model under the most plausible assumption of its distance as 200 pc is much fainter than for a typical protostar, and even fainter than other VeLLOs studied previously. Note, however, that its inferred luminosity may be uncertain by a factor of 2-3 if we consider two extreme values of the distance of L328-IRS (125 or 310 pc). Low angular resolution observations of CO do not show any clear evidence of a molecular outflow activity. But broad line widths toward L328, and Spitzer and near-infrared images showing nebulosity possibly tracing an outflow cavity, strongly suggest the existence of outflow activity. Provided that an envelope of at most ~0.1 Msunis the only mass accretion reservoir for L328-IRS, and the star formation efficiency is close to the canonical value ~30%, L328-IRS has not yet accreted more than 0.05 Msun. At the assumed distance of 200 pc, L328-IRS is destined to be a brown dwarf.

We present Spitzer infrared observations of the starless core L429. The IR images of this core show an absorption feature, caused by the dense core material, at wavelengths <= 70 micron. The core has a steep density profile, and reaches A_V > 35 mag near the center. We show that L429 is either collapsing or in a near-collapse state.

We present interferometric observations in the 12CO (2-1) line and at 1.3 mm dust continuum of the low-mass protostellar binary system in the cometary globule CG30, using the Submillimeter Array. The dust continuum images resolve two compact sources (CG30N and CG30S), with a linear separation of ~8700 AU and total gas masses of ~1.4 and ~0.6 M_sun, respectively. With the CO images, we discover two high-velocity bipolar molecular outflows, driven by the two sources. The two outflows are nearly perpendicular to each other, showing a quadrupolar morphology. The northern bipolar outflow extends along the southeast (redshifted, with a velocity up to ~23 km/s) and northwest (blueshifted, velocity up to ~30 km/s) directions, while the southern pair has an orientation from southwest (blueshifted, velocity up to 13 km/s) to northeast (redshifted, velocity up to ~41 km/s). The outflow mass of the northern pair, driven by the higher mass source CG30N, is ~9 times larger than that of the southern pair. The discovery of the quadrupolar molecular outflow in the CG30 protobinary system, as well as the presence of other quadrupolar outflows associated with binary systems, demonstrate that the disks in (wide) binary systems are not necessarily co-aligned after fragmentation.

RCW~38 is a uniquely young ($<$1 Myr), embedded ($A_V \sim 10$) stellar cluster surrounding a pair of early O stars ($\sim$O5.5) and is one of the few regions within 2 kpc other than Orion to contain over 1000 members. X-ray and deep near-infrared observations reveal a dense cluster with over 200 X-ray sources and 400 infrared sources embedded in a diffuse hot plasma within a 1 pc diameter. The central O star has evacuated its immediate surroundings of dust, creating a wind bubble $\sim$0.1 pc in radius that is confined by the surrounding molecular cloud, as traced by millimeter continuum and molecular line emission. The interface between the bubble and cloud is a region of warm dust and ionized gas, which shows evidence for ongoing star formation. Extended warm dust is found throughout a 2--3 pc region and coincides with extended X-ray plasma. This is evidence that the influence of the massive stars reaches beyond the confines of the O star bubble. RCW~38 appears similar in structure to RCW~49 and M~20 but is at an earlier evolutionary phase. RCW~38 appears to be a blister compact H{\small II} region lying just inside the edge of a giant molecular cloud.

The Ara OB1a association is one of the closest sites where triggered star formation is visible for multiple generations of massive stars. At about 1.3 kpc distance, it contains complex environments including cleared young clusters, embedded infrared clusters, CO clouds with no evidence of star formation, and clouds with evidence of ongoing star formation. In this review we discuss the research on this region spanning the last half-century. It has been proposed that the current configuration is the result of an expanding wave of neutral gas set in motion between 10--40 million years ago in combination with photoionization from the current epoch.

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 search for all embedded protostars with internal luminosities < 1 solar luminosity in the sample of nearby, low-mass star-forming regions surveyed by the Spitzer Space Telescope c2d Legacy Project. The internal luminosity (Lint) of a source is the luminosity of the central source and excludes luminosity arising from external heating. On average, the c2d data are sensitive to embedded protostars with Lint > 4E-3 (d/140 pc)^2 solar luminosities, a factor of 25 better than the sensitivity of IRAS to such objects. We present selection criteria used to identify candidates from the Spitzer data and examine complementary data to decide whether each candidate is truly an embedded protostar. We find a tight correlation between the 70 micron flux and internal luminosity of a protostar, an empirical result based on observations and two-dimensional radiative transfer models of protostars. We identify 50 embedded protostars with Lint < 1 solar luminosities; 15 have Lint < 0.1 solar luminosities. The intrinsic distribution of source luminosities increases to lower luminosities. While we find sources down to the above sensitivity limit, indicating that the distribution may extend to luminosities lower than probed by these observations, we are able to rule out a continued rise in the distribution below 0.1 solar luminosities. Between 75-85% of cores classified as starless prior to being observed by Spitzer remain starless to our luminosity sensitivity; the remaining 15-25% harbor low-luminosity, embedded protostars. We compile complete Spectral Energy Distributions for all 50 objects and calculate standard evolutionary signatures, and argue that these objects are inconsistent with the simplest picture of star formation wherein mass accretes from the core onto the protostar at a constant rate.

AIMS: To study the structure of nearby (< 500 pc) dense starless and star-forming cores with the particular goal to identify and understand evolutionary trends in core properties, and to explore the nature of Very Low Luminosity Objects (< 0.1 L_sun; VeLLOs). METHODS: Using the MAMBO bolometer array, we create maps unusually sensitive to faint (few mJy per beam) extended (approx. 5 arcmin) thermal dust continuum emission at 1.2 mm wavelength. Complementary information on embedded stars is obtained from Spitzer, IRAS, and 2MASS. RESULTS: Our maps are very rich in structure, and we characterize extended emission features (``subcores'') and compact intensity peaks in our data separately to pay attention to this complexity. We derive, e.g., sizes, masses, and aspect ratios for the subcores, as well as column densities and related properties for the peaks. Combination with archival infrared data then enables the derivation of bolometric luminosities and temperatures, as well as envelope masses, for the young embedded stars. CONCLUSIONS: (abridged) Starless and star-forming cores occupy the same parameter space in many core properties; a picture of dense core evolution in which any dense core begins to actively form stars once it exceeds some fixed limit in, e.g., mass, density, or both, is inconsistent with our data. Comparison of various evolutionary indicators for young stellar objects in our sample (e.g., bolometric temperatures) reveals inconsistencies between some of them, possibly suggesting a revision of some of these indicators.

We present interferometric observations with resolution of ~3 arcsecs of the isolated, low-mass protostellar double cores CG30 and BHR71 in the N2$H+(1-0) line and at 3mm dust continuum, using the Australian Telescope Compact Array (ATCA). The results are complemented by infrared data from the Spitzer Space Telescope. In CG30, the 3mm dust continuum images resolve two compact sources with a separation of ~21.7 arcsecs (~8700 AU). In BHR71, strong dust continuum emission is detected at the position of the mid-infrared source IRS1, while only weak emission is detected from the secondary mid-infrared source IRS2. Assuming optically thin 3mm dust continuum emission, we derive hydrogen gas masses of 0.05--2.1 $M_\odot$ for the four sub-cores. N2H+(1-0) line emission is detected in both CG30 and BHR71, and is spatially associated with the thermal dust continuum emission. We derive the velocity fields and find symmetric velocity gradients in both sources. Assuming that these gradients are due to core rotation, we estimate the specific angular momenta and ratios of rotational energy to gravitational energy for all cores. We also find that the N2H+ emission is strongly affected by the outflows, both in terms of entrainment and molecule destruction. $Spitzer$ images show the mid-infrared emission from all four sub-cores. All four sources appear to drive their own outflows. Based on the ATCA and $Spitzer$ observations, we construct spectral energy distributions (SEDs) and derive temperatures and luminosities for all cores. Based on the morphology and velocity structure, we suggest that the sub-cores in CG30 were formed by initial fragmentation of a filamentary prestellar core, while those in BHR71 could originate from rotational fragmentation of a single collapsing protostellar core.

We report the discovery of a nearby, embedded cluster of young stellar objects, associated filamentary infrared dark cloud, and 4.5 micron shock emission knots from outflows detected in Spitzer/IRAC mid-infrared imaging of the Serpens-Aquila Rift obtained as part of the Spitzer Gould Belt Legacy Survey. We also present radial velocity measurements of the region from molecular line observations obtained with the Submillimeter Array (SMA) that suggest the cluster is co-moving with the Serpens Main embedded cluster 3 degrees to the north. We therefore assign it the same distance, 260 pc. The core of the new cluster, which we call Serpens South, is composed of an unusually large fraction of protostars (77%) at high mean surface density (>430 pc^-2) and short median nearest neighbor spacing (3700 AU). We perform basic cluster structure characterization using nearest neighbor surface density mapping of the YSOs and compare our findings to other known clusters with equivalent analyses available in the literature.

We report on the results of an approximately 90 ks Chandra observation of a complex region that hosts multiple sites of recent and active star formation in ARA OB1a. The field is centered on the embedded cluster RCW 108-IR and includes and a large portion of the open cluster NGC 6193. We detect over 420 X-ray sources in the field and combined these data with deep near-IR, Spitzer/IRAC and MSX mid-IR data. We find about 360 of the X-ray sources have near--IR counterparts. We divide the region into 5 parts based on the X-ray point source characteristics and extended 8 micron emission. The most clearly defined regions are the central region - identified by embedded sources with high luminosities in the both the near-IR and X-ray as well as high X-ray temperatures (about 3 keV) and the eastern region - identified by low extinction and 1 keV X-ray temperatures. Other regions, identified by their directional relationship to RCW 108-IR are less uniform - representing combinations of the first two regions, independent star formation epochs, or both. Over 18% percent of the cluster members with over 100 counts exhibit flares. Overall about 50% of the stars appear to have optically thick disks when IRAC data are employed. The largest fraction of X-ray sources are best described as possessing some disk material via a more detailed extinction fitting. We estimate that the total number of pre--main sequence stars in the field is about 1600. Approximately 800 are confined to (1.1 pc) central region.

We have observed the deuterated gas in the high-mass star formation region IRAS 05345+3157 at high-angular resolution, in order to determine the morphology and the nature of such gas. We have mapped the N2H+ (1-0) line with the Plateau de Bure Interferometer, and the N2D+ (3-2) and N2H+ (3-2) lines with the Submillimeter Array. The N2D+ (3-2) integrated emission is concentrated in two condensations, with masses of 2-3 and 9 M_sun and diameters of 0.05 and 0.09 pc, respectively. The high deuterium fractionation (0.1) and the line parameters in the N2D+ condensations indicate that they are likely low- to intermediate-mass pre-stellar cores, even though other scenarios are possible.

We have mapped the proto-binary source IRAS 16293-2422 in CO 2-1, 13CO 2-1, and CO 3-2 with the Submillimeter Array (SMA). The maps with resolution of 1".5-5" reveal a single small scale (~3000 AU) bipolar molecular outflow along the east-west direction. We found that the blueshifted emission of this small scale outflow mainly extends to the east and the redshifted emission to the west from the position of IRAS 16293A. A comparison with the morphology of the large scale outflows previously observed by single-dish telescopes at millimeter wavelengths suggests that the small scale outflow may be the inner part of the large scale (~15000 AU) E-W outflow. On the other hand, there is no clear counterpart of the large scale NE-SW outflow in our SMA maps. Comparing analytical models to the data suggests that the morphology and kinematics of the small scale outflow can be explained by a wide-angle wind with an inclination angle of ~30-40 degrees with respect to the plane of the sky. The high resolution CO maps show that there are two compact, bright spots in the blueshifted velocity range. An LVG analysis shows that the one located 1" to the east of source A is extremely dense, n(H_2)~10^7 cm^-3, and warm, T_kin >55 K. The other one located 1" southeast of source B has a higher temperature of T_kin >65 K but slightly lower density of n(H_2)~10^6 cm^-3. It is likely that these bright spots are associated with the hot core-like emission observed toward IRAS 16293. Since both two bright spots are blueshifted from the systemic velocity and are offset from the protostellar positions, they are likely formed by shocks.

CONTEXT: The chemistry in the inner few thousand AU of accreting envelopes around young stellar objects is predicted to vary greatly with far-UV and X-ray irradiation by the central star. Aim We search for molecular tracers of high-energy irradiation by the protostar in the hot inner envelope. METHODS: The Submillimeter Array (SMA) has observed the high-mass star forming region AFGL 2591 in lines of CS, SO, HCN, HCN(v2=1), and HC15N with 0.6" resolution at 350 GHz probing radial scales of 600-3500 AU for an assumed distance of 1 kpc. The SMA observations are compared with the predictions of a chemical model fitted to previous single-dish observations. RESULTS: The CS and SO main peaks are extended in space at the FWHM level, as predicted in the model assuming protostellar X-rays. However, the main peak sizes are found smaller than modeled by nearly a factor of 2. On the other hand, the lines of CS, HCN, and HC15N, but not SO and HCN(v2=1), show pedestal emissions at radii of about 3500 AU that are not predicted. All lines except SO show a secondary peak within the approaching outflow cone. A dip or null in the visibilities caused by a sharp decrease in abundance with increasing radius is not observed in CS and only tentatively in SO. CONCLUSIONS: The emission of protostellar X-rays is supported by the good fit of the modeled SO and CS amplitude visibilities including an extended main peak in CS. The broad pedestals can be interpreted by far-UV irradiation in a spherically non-symmetric geometry, possibly comprising outflow walls on scales of 3500 -- 7000 AU. The extended CS and SO main peaks suggest sulfur evaporation near the 100 K temperature radius.

We present molecular line observations of L1251B, a small group of pre- and protostellar objects, and its immediate environment in the dense C18O core L1251E. These data are complementary to near-infrared, submillimeter and millimeter continuum observations reported by Lee et al. (2006, ApJ, 648, 491; Paper I). The single-dish data of L1251B described here show very complex kinematics including infall, rotation and outflow motions, and the interferometer data reveal these in greater detail. Interferometer data of N2H+ 1-0 suggest a very rapidly rotating flattened envelope between two young stellar objects, IRS1 and IRS2. Also, interferometer data of CO 2-1 resolve the outflow associated with L1251B seen in single-dish maps into a few narrow and compact components. Furthermore, the high resolution data support recent theoretical studies of molecular depletions and enhancements that accompany the formation of protostars within dense cores. Beyond L1251B, single-dish data are also presented of a dense core located ~150" to the east that, in Paper I, was detected at 850 micron but has no associated point sources at near- and mid-infrared wavelengths. The relative brightness between molecules, which have different chemical timescales, suggests it is less chemically evolved than L1251B. This core may be a site for future star formation, however, since line profiles of HCO+, CS, and HCN show asymmetry with a stronger blue peak, which is interpreted as an infall signature.

We have mapped the protostellar jet HH 211 in 342 GHz continuum, SiO ($J=8-7$), and CO ($J=3-2$) emission at $\sim$ \arcs{1} resolution with the Submillimeter Array (SMA). Thermal dust emission is seen in continuum at the center of the jet, tracing an envelope and a possible optically thick compact disk (with a size $<$ 130 AU) around the protostar. A knotty jet is seen in CO and SiO as in \H2{}, but extending closer to the protostar. It consists of a chain of knots on each side of the protostar, with an interknot spacing of $\sim$ \arcs{2}$-$\arcs{3} or 600$-$900 AU and the innermost pair of knots at only $\sim$ \arcsa{1}{7} or 535 AU from the protostar. These knots likely trace unresolved internal (bow) shocks (i.e., working surfaces) in the jet, with a velocity range up to $\sim$ 25 \vkm{}. The two-sided mass-loss rate of the jet is estimated to be $\sim (0.7-2.8)\times 10^{-6}$ \solarmass{} yr$^{-1}$. The jet is episodic, precessing, and bending. A velocity gradient is seen consistently across two bright SiO knots (BK3 and RK2) perpendicular to the jet axis, with $\sim$ 1.5$\pm$0.8 \vkm{} at $\sim$ 30$\pm$15 AU, suggesting a presence of a jet rotation. The launching radius of the jet, derived from the potential jet rotation, is $\sim$ 0.15$-$0.06 AU in the inner disk.

We present maps of 7.78 square degrees of the Lupus molecular cloud complex at 24, 70, and $160\:\mu$m. They were made with the Spitzer Space Telescope's Multiband Imaging Photometer for Spitzer (MIPS) instrument as part of the Spitzer Legacy Program, ``From Molecular Cores to Planet-Forming Disks'' (c2d). The maps cover three separate regions in Lupus, denoted I, III, and IV. We discuss the c2d pipeline and how our data processing differs from it. We compare source counts in the three regions with two other data sets and predicted star counts from the Wainscoat model. This comparison shows the contribution from background galaxies in Lupus I. We also create two color magnitude diagrams using the 2MASS and MIPS data. From these results, we can identify background galaxies and distinguish them from probable young stellar objects. The sources in our catalogs are classified based on their spectral energy distribution (SED) from 2MASS and Spitzer wavelengths to create a sample of young stellar object candidates. From 2MASS data, we create extinction maps for each region and note a strong corresponence between the extinction and the $160\:\mu$m emission. The masses we derived in each Lupus cloud from our extinction maps are compared to masses estimated from $^{13}$CO and C$^{18}$O and found to be similar to our extinction masses in some regions, but significantly different in others. Finally, based on our color-magnitude diagrams, we selected 12 of our reddest candidate young stellar objects for individual discussion. Five of the 12 appear to be newly-discovered YSOs.

Observations by the Cores to Disk Legacy Team with the Spitzer Space Telescope have identified a low luminosity, mid-infrared source within the dense core, Lynds 1014, which was previously thought to harbor no internal source. Followup near-infrared and submillimeter interferometric observations have confirmed the protostellar nature of this source by detecting scattered light from an outflow cavity and a weak molecular outflow. In this paper, we report the detection of cm continuum emission with the VLA. The emission is characterized by a quiescent, unresolved 90 uJy 6 cm source within 0.2" of the Spitzer source. The spectral index of the quiescent component is $\alpha = 0.37\pm 0.34$ between 6 cm and 3.6 cm. A factor of two increase in 6 cm emission was detected during one epoch and circular polarization was marginally detected at the $5\sigma$ level with Stokes {V/I} $= 48 \pm 16$% . We have searched for 22 GHz H2O maser emission toward L1014-IRS, but no masers were detected during 7 epochs of observations between June 2004 and December 2006. L1014-IRS appears to be a low-mass, accreting protostar which exhibits cm emission from a thermal jet or a wind, with a variable non-thermal emission component. The quiescent cm radio emission is noticeably above the correlation of 3.6 cm and 6 cm luminosity versus bolometric luminosity, indicating more radio emission than expected. We characterize the cm continuum emission in terms of observations of other low-mass protostars, including updated correlations of centimeter continuum emission with bolometric luminosity and outflow force, and discuss the implications of recent larger distance estimates on the physical attributes of the protostar and dense molecular core.

Forty molecular cloud cores in the southern hemisphere from the initial Spitzer Space Telescope Cores-to-Disks (c2d) Legacy program source list have been surveyed in 13CO(2-1), 12CO(4-3), and 12CO(7-6) with the Antarctic Submillimeter Telescope and Remote Observatory (AST/RO). The cores, ten of which contain embedded sources, are located mostly in the Vela, Ophiuchus, Lupus, Chamaeleon, Musca, and Scorpius complexes. 12CO(7-6) emission was undetected in all 40 clouds. We present data of 40 sources in 13CO(2-1) and 12CO(4-3), significant upper limits of 12CO(7-6), as well as a statistical analysis of the observed properties of the clouds. We find the typical 13CO(2-1) linewidth to be 2.0 km/s for cores with embedded stars, and 1.8 km/s for all others. The typical 12CO(4-3) linewidth is 2.6 to 3.7 km/s for cores with known embedded sources, and 1.6 to 2.3 km/s for all others. The average 13CO column density derived from the line intensities was found to be 1.9 x 10^15 cm^(-2) for cores with embedded stars, and 1.5 x 10^15 cm^(-2) for all others. The average kinetic temperature in the molecular cores, determined through a Large Velocity Gradient analysis of a set of nine cores, has an average lower limit of 16 K and an average upper limit of 26 K. The average molecular hydrogen density has an average lower limit of 10^2.9 cm^(-3) and an average upper limit of 10^3.3 cm^(-3) for all cores. For a different subset of nine cores, we have derived masses. They range from 4 to 255 M_sun. Overall, our c2d sample of southern molecular cores has a range of properties (linewidth, column density, size, mass, embedded stars) similar to those of past studies.

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

We present observations of 10.6 square degrees of the Perseus molecular cloud at 24, 70, and 160 microns with the Spitzer Space Telescope Multiband Imaging Photometer for Spitzer (MIPS). The image mosaics show prominent, complex extended emission dominated by illuminating B stars on the East side of the cloud, and by cold filaments of 160 micron emission on the West side. Of 3950 point sources identified at 24 microns, 1141 have 2MASS counterparts. A quarter of these populate regions of the Ks vs. Ks-[24] diagram that are distinct from stellar photospheres and background galaxies, and thus are likely to be cloud members with infrared excess. Nearly half (46%) of these 24 micron excess sources are distributed outside the IC 348 and NGC 1333 clusters. NGC 1333 shows the highest fraction of stars with flat or rising spectral energy distributions (28%), while Class II SEDs are most common in IC 348. These results are consistent with previous relative age determinations for the two clusters. The intercluster region contains several tightly clumped (r~0.1 pc) young stellar aggregates whose members exhibit a wide variety of infrared spectral energy distributions characteristic of different circumstellar environments. One possible explanation is a significant age spread among the aggregate members, such that some have had time to evolve more than others. Alternatively, if the aggregate members all formed at roughly the same time, then remarkably rapid circumstellar evolution would be required to account for the association of Class I and Class III sources at ages <~1 Myr. We highlight important results for several other objects as well (full abstract in the paper).

HH 212 is a nearby (460 pc) protostellar jet discovered in H$_2$ powered by a Class 0 source, IRAS 05413-0104, in the L1630 cloud of Orion. It is highly collimated and symmetric with matched pairs of bow shocks on either side of the source. We have mapped it in 850 $\mu$m continuum, SiO ($J=8-7$), CO ($J=3-2$), SO ($N_J=8_9-7_8$), HCO$^+$ ($J = 4-3$), and H$^{13}$CO$^+$ ($J = 4-3$) emission simultaneously at $\sim$ 1$''$ resolution with the Submillimeter Array (SMA). Thermal dust emission is seen in continuum around the source, mainly arising from an inner envelope (i.e., the inner part of a previously seen flattened envelope) and a possible disk. The inner envelope is also seen with rotation in CO, HCO$^+$, and probably SO. Like H$_2$ emission, CO and SiO emission are seen along the jet axis but extending closer to the source, tracing the bow shocks with a broad range of velocities and the continuous structures in between. SO emission is seen only around the source, forming a jetlike structure extending along the jet axis from the source, likely tracing the jet near the launching region. The jet is episodic and bending. It may also be slightly precessing as the jetlike SO structure shows a slight S-shaped symmetry about the source. A hint of jet rotation is also seen across the jet axis. Internal outflow shells are seen in CO and HCO$^+$, associated with the bow shocks in the inner part of the jet. The bases of the HCO$^+$ shells are seen with a hint of rotation similar to that seen in the inner envelope, probably consisted mainly of the material extended from the inner envelope and even the possible disk.

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 present the results of a submillimeter survey of 53 low-mass dense cores with the Submillimeter High Angular Resolution Camera II (SHARC-II). The survey is a follow-up project to the Spitzer Legacy Program ``From Molecular Cores to Planet-Forming Disks'', with the purpose being to create a complete data set of nearby low-mass dense cores from the infrared to the millimeter. We present maps of 52 cores at 350 microns and three cores at 450 microns, two of which were observed at both wavelengths. Of these 52 cores, 41 were detected by SHARC-II: 32 contained one submillimeter source while 9 contained multiple sources. For each submillimeter source detected, we report various source properties including source position, fluxes in various apertures, size, aspect ratio, and position angle. For the 12 cores that were not detected we present upper limits. The sources detected by SHARC-II have, on average, smaller sizes at the 2sigma contours than those derived from longer-wavelength bolometer observations. We conclude that this is not caused by a failure to integrate long enough to detect the full extent of the core; instead it arises primarily from the fact that the observations presented in this survey are insensitive to smoothly varying extended emission. We find that SHARC-II observations of low-mass cores are much better suited to distinguishing between starless and protostellar cores than observations at longer wavelengths. Very Low Luminosity Objects, a new class of objects being discovered by the Spitzer Space Telescope in cores previously classified as starless, look very similar at 350 microns to other cores with more luminous protostars.

We present IRAC (3.6, 4.5, 5.8, and 8.0 micron) observations of the Chamaeleon II molecular cloud. The observed area covers about 1 square degree defined by $A_V >2$. Analysis of the data in the 2005 c2d catalogs reveals a small number of sources (40) with properties similar to those of young stellaror substellar objects (YSOs). The surface density of these YSO candidates is low, and contamination by background galaxies appears to be substantial, especially for sources classified as Class I or flat SED. We discuss this problem in some detail and conclude that very few of the candidate YSOs in early evolutionary stages are actually in the Cha II cloud. Using a refined set of criteria, we define a smaller, but more reliable, set of 24 YSO candidates.

Infrared images of the dark cloud core B59 were obtained with the Spitzer Space Telescope as part of the "Cores to Disks" Legacy Science project. Photometry from 3.6-70 microns indicates at least 20 candidate low-mass young stars near the core, more than doubling the previously known population. Out of this group, 13 are located within about 0.1 pc in projection of the molecular gas peak, where a new embedded source is detected. Spectral energy distributions span the range from small excesses above photospheric levels to rising in the mid-infrared. One other embedded object, probably associated with the millimeter source B59-MMS1, with a bolometric luminosity L(bol) roughly 2 L(sun), has extended structure at 3.6 and 4.5 microns, possibly tracing the edges of an outflow cavity. The measured extinction through the central part of the core is A(V) greater than of order 45 mag. The B59 core is producing young stars with a high efficiency.

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 observed the clustered star forming complex NGC1333 with the BIMA and FCRAO telescopes in the transitions HCO+(1-0) and N2H+(1-0) over an 11'x11' area with resolution ~10" (0.015pc). The N2H+ emission follows very closely the submillimeter dust continuum emission, while HCO+ emission appears more spatially extended and also traces outflows. We have identified 93 N2H+ cores using the CLUMPFIND algorithm, and we derive N2H+ core masses between 0.05 and 2.5M_sun, with uncertainties of a factor of a few, dominated by the adopted N2H+ abundance. From a comparison with virial masses, we argue that most of these N2H+ cores are likely to be bound, even at the lowest masses, suggesting that the cores do not trace transient structures, and implies the entire mass distribution consists of objects that can potentially form stars. We find that the mass distribution of N2H+ cores resembles the field star IMF, which suggests that the IMF is locked in at the pre-stellar stage of evolution. We find that the N2H+ cores associated with stars identified from Spitzer infrared images have a flat mass distribution. This might be because lower mass cores lose a larger fraction of their mass when forming a star. Even in this clustered environment, we find no evidence for ballistic motions of the cores relative to their lower density surroundings traced by isotopic CO emission, though this conclusion must remain tentative until the surroundings are observed at the same high resolution as the N2H+.

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.

We present Spitzer Space Telescope observations of the "evolved starless core" L1521F which reveal the presence of a very low luminosity object (L < 0.07 Lsun). The object, L1521F-IRS, is directly detected at mid-infrared wavelengths (>5 micron) but only in scattered light at shorter infrared wavelengths, showing a bipolar nebula oriented east-west which is probably tracing an outflow cavity. The nebula strongly suggests that L1521F-IRS is embedded in the L1521F core. Thus L1521F-IRS is similar to the recently discovered L1014-IRS and the previously known IRAM 04191 in its substellar luminosity and dense core environment. However these objects differ significantly in their core density, core chemistry, and outflow properties, and some may be destined to be brown dwarfs rather than stars.

We report the first detections of the Class 0 protostellar source IRAM 04191+1522 at wavelengths shortward of 60 microns with the Spitzer Space Telescope. We see extended emission in the Spitzer images that suggests the presence of an outflow cavity in the circumstellar envelope. We combine the Spitzer observations with existing data to form a complete dataset ranging from 3.6 to 1300 microns and use these data to construct radiative transfer models of the source. We conclude that the internal luminosity of IRAM 04191+1522, defined to be the sum of the luminosity from the internal sources (a star and a disk), is L_int = 0.08 +/- 0.04 L_sun, placing it among the lowest luminosity protostars known. Though it was discovered before the launch of the Spitzer Space Telescope, IRAM 04191+1522 falls within a new class of Very Low Luminosity Objects being discovered by Spitzer. Unlike the two other well-studied objects in this class, which are associated either with weak, compact outflows or no outflows at all, IRAM 04191+1522 has a well-defined molecular outflow with properties consistent with those expected based on relations derived from higher luminosity (L_int > 1 L_sun) protostars. We discuss the difficulties in understanding IRAM 04191+1522 in the context of the standard model of star formation, and suggest a possible explanation for the very low luminosity of this source.

We present a comprehensive study of a low-mass star-forming region,L1251B, at wavelengths from the near-infrared to the millimeter. L1251B, where only one protostar, IRAS 22376+7455, was known previously, is confirmed to be a small group of protostars based on observations with the Spitzer Space Telescope. The most luminous source of L1251B is located 5" north of the IRAS position. A near-infrared bipolar nebula, which is not associated with the brightest object and is located at the southeast corner of L1251B, has been detected in the IRAC bands. OVRO and SMA interferometric observations indicate that the brightest source and the bipolar nebula source in the IRAC bands are deeply embedded disk sources.Submillimeter continuum observations with single-dish telescopes and the SMA interferometric observations suggest two possible prestellar objects with very high column densities. Outside of the small group, many young stellar object candidates have been detected over a larger region of 12' x 12'. Extended emission to the east of L1251B has been detected at 850 micron; this "east core" may be a site for future star formation since no point source has been detected with IRAC or MIPS. This region is therefore a possible example of low-mass cluster formation, where a small group of pre- and protostellar objects (L1251B) is currently forming, alongside a large starless core (the east core).

We report on results of a 96.7 ks Chandra observation of one of the youngest, most embedded and massive young stellar clusters studied in X-rays - RCW 38. We detect 460 sources in the field of which 360 are confirmed to be associated with the RCW~38 cluster. The cluster members range in luminosity from ~10^30 ergs s^-1 to 10^33.5 ergs s^-1. Over ten percent of the cluster members with over 100 counts exhibit flares while about 15% of cluster members with over 30 counts are variable. Of the sources identified as cluster members, 160 have near-infrared (NIR) counterparts either in the 2MASS database or detected via VLT observations. Of these about 20% appear to have optically thick disks. An additional 353 members are identified through NIR observations of which at least 50% possess optically thick disks. We fit over 100 X-ray sources as absorbed Raymond-Smith type plasmas and find the column to the cluster members varies from 10^21.5 to 10^23 per sq cm. We compare the gas to dust absorption signatures in these stars and find Nh= A_v x 2 x 10^21 cm^-2. We find that the cluster contains 31 candidate OB stars and is centered about 10" (0.1 pc) west of the primary source of the ionization, the O5 star IRS 2. The cluster has a peak central density of about 400 X-ray sources pc^-2. We estimate that the total cluster membership exceeds 2000 stars.

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 observed the HH 211 jet in the submillimeter continuum and the CO(3-2) and SiO(8-7) transitions with the Submillimeter Array. The continuum source detected at the center of the outflow shows an elongated morphology, perpendicular to the direction of the outflow axis. The high-velocity emission of both molecules shows a knotty and highly collimated structure. The SiO(8-7) emission at the base of the outflow, close to the driving source, spans a wide range of velocities, from -20 up to 40 km s^{-1}. This suggests that a wide-angle wind may be the driving mechanism of the HH 211 outflow. For distances greater than 5" (1500 AU) from the driving source, emission from both transitions follows a Hubble-law behavior, with SiO(8-7) reaching higher velocities than CO(3-2), and being located upstream of the CO(3-2) knots. This indicates that the SiO(8-7) emission is likely tracing entrained gas very close to the primary jet, while the CO(3-2) is tracing less dense entrained gas. From the SiO(5-4) data of Hirano et al. we find that the SiO(8-7)/SiO(5-4) brightness temperature ratio along the jet decreases for knots far from the driving source. This is consistent with the density decreasing along the jet, from (3-10)x10^6 cm^{-3} at 500 AU to (0.8-4)x10^6 cm^{-3} at 5000 AU from the driving source.

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

HH 212 is a highly collimated jet discovered in H2 powered by a young Class 0 source, IRAS 05413-0104, in the L1630 cloud of Orion. We have mapped around it in 1.33 mm continuum, 12CO ($J=2-1$), 13CO ($J=2-1$), C18O ($J=2-1$), and SO ($J_K = 6_5-5_4$) emission at $\sim$ \arcs{2.5} resolution with the Submillimeter Array. A dust core is seen in the continuum around the source. A flattened envelope is seen in C18O around the source in the equator perpendicular to the jet axis, with its inner part seen in 13CO. The structure and kinematics of the envelope can be roughly reproduced by a simple edge-on disk model with both infall and rotation. In this model, the density of the disk is assumed to have a power-law index of $p=-1.5$ or -2, as found in other low-mass envelopes. The envelope seems dynamically infalling toward the source with slow rotation because the kinematics is found to be roughly consistent with a free fall toward the source plus a rotation of a constant specific angular momentum. A 12CO outflow is seen surrounding the H2 jet, with a narrow waist around the source. Jetlike structures are also seen in 12CO near the source aligned with the H2 jet at high velocities. The morphological relationship between the H2 jet and the 12CO outflow, and the kinematics of the 12CO outflow along the jet axis are both consistent with those seen in a jet-driven bow shock model. SO emission is seen around the source and the H2 knotty shocks in the south, tracing shocked emission around them.

Recently, the Spitzer Space Telescope discovered L1014-IRS, a mid-infrared source with protostellar colors, toward the heretofore "starless" core L1014. We present deep near-infrared observations that show a scattered light nebula extending from L1014-IRS. This nebula resembles those typically associated with protostars and young stellar objects, tracing envelope cavities presumably evacuated by an outflow. The northern lobe of the nebula has an opening angle of ~100 degrees, while the southern lobe is barely detected. Its morphology suggests that the bipolar cavity and inferred protostellar disk is not inclined more than 30 degrees from an edge-on orientation. The nebula extends at least 8" from the source at Ks, strongly suggesting that L1014-IRS is embedded within L1014 at a distance of 200 pc rather than in a more distant cloud associated with the Perseus arm at 2.6 kpc. In this case, the apparently low luminosity of L1014-IRS, 0.090 Lsun, is consistent with it having a substellar mass. However, if L1014-IRS is obscured by a circumstellar disk, its luminosity and inferred mass may be greater. Using near-infrared colors of background stars, we investigate characteristics of the L1014 molecular cloud core. We determine a mass of 3.6 Msun for regions of the core with Av > 2 magnitudes. A comparison of the radial extinction profile of L1014 with other cores suggests that L1014 may be among the most centrally condensed cores known, perhaps indicative of the earliest stages of brown dwarf or star formation processes.

We report ``infall asymmetry'' in the HCO$^+$ (1--0) and (3--2) lines toward NGC 1333, extended over $\sim 0.39 {\rm pc}^2$, a larger extent than has been reported be fore, for any star-forming region. The infall asymmetry extends over a major portion of the star-forming complex, and is not limited to a single protostar, or to a single dense core, or to a single spectral line. It seems likely that the infall asymmetry represents inward motions, and that these motions are physically associated with the complex. Both blue-asymmetric and red-asymmetric lines are seen, but in both the (3--2) and (1--0) lines of HCO$^+$ the vast majority of the asymmetric lines are blue, indicating inward motions. The (3--2) line, tracing denser gas, has the spectra with the strongest asymmetry and these spectra are associated with the protostars IRAS 4A and 4B, which most likely indicates a warm central source is affecting the line profiles. The (3--2) and (1--0) lines usually have the same sense of asymmetry in common positions, but their profiles differ significantly, and the (1--0) line appears to trace motions on much larger spatial scales than does the (3--2) line. Line profile models fit the spectra well, but do not strongly constrain their parameters. The mass accretion rate of the inward motions is of order 10$^{-4}$ M$_\odot$/yr, similar to the ratio of stellar mass to cluster age.

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

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.

We present the results of a survey of a small sample (14) of low-mass protostars (L_IR < 10^3 Lsun) for 6.7 GHz methanol maser emission performed using the ATNF Parkes radio telescope. No new masers were discovered. We find that the lower luminosity limit for maser emission is near 10^3 Lsun, by comparison of the sources in our sample with previously detected methanol maser sources. We examine the IRAS properties of our sample and compare them with sources previously observed for methanol maser emission, almost all of which satisfy the Wood & Churchwell criterion for selecting candidate UCHII regions. We find that about half of our sample satisfy this criterion, and in addition almost all of this subgroup have integrated fluxes between 25 and 60 microns that are similar to sources with detectable methanol maser emission. By identifying a number of low-mass protostars in this work and from the literature that satisfy the Wood & Churchwell criterion for candidate UCHII regions, we show conclusively for the first time that the fainter flux end of their sample is contaminated by lower-mass non-ionizing sources, confirming the suggestion by van der Walt and Ramesh & Sridharan.

Using the recently-commissioned Submillimeter Array (SMA), we have detected several complex organic molecules, including (CH3)2O, C2H5OH, C2H5CN, and tentatively CH2CDCN, toward the protostellar hot cores of IRAS 16293-2422 at arcsecond-resolution (~< 400 AU in radius). Vibrationally excited transitions of SO, SO2 and HCN with energy levels up to 1800 K were also observed. In addition to the other organic molecules (HC3N, CH2CO, CH3OH, CH2CHCN and HCOOCH3) previously reported by us (Kuan et al. 2004), these results clearly indicate the existence of a rich organic chemistry in low-mass `hot corinos'. From the observation of optically thin HC15N emission, we conclude I16293A is a rotating circumstellar disk lying along the north-south direction ~10 deg to the east and with an inclination ~30 deg to the sky. We suggest that the observed vibrational SO and SO2 emission may originate from shock waves near or in the circumstellar disks. Between the two cores, we find a strong anticorrelation in emission from C2H5OH and C2H5CN. The relative contribution of gas phase and grain-surface chemistries to the production of the observed complex molecules is discussed. We point out the shortcomings underlying recent claims that all the O-bearing organics are formed on grains. The presence of so many well-known interstellar molecules in solar-type hot corinos strengthens the link between molecular cloud chemistry, the starting materials of protoplanetary disks such as the protosolar nebula, and the composition of comets. Establishing the fine details of this connection is crucial in answering fundamental questions concerning the importance of galactic astrochemistry for astrobiology.

We present maps of over 1.5 square degrees in Chamaeleon (Cha) II at 24, 70, and 160 micron observed with the Spitzer Space Telescope Multiband Imaging Photometer for Spitzer (MIPS) and a 1.2 square degree millimeter map from SIMBA on the Swedish-ESO Submillimetre Telescope (SEST). The c2d Spitzer Legacy Team's data reduction pipeline is described in detail. Over 1500 24 micron sources and 41 70 micron sources were detected by MIPS with fluxes greater than 10-sigma. More than 40 potential YSOs are identified with a MIPS and 2MASS color-color diagram and by their spectral indices, including two previously unknown sources with 24 micron excesses. Our new SIMBA millimeter map of Cha II shows that only a small fraction of the gas is in compact structures with high column densities. The extended emission seen by MIPS is compared with previous CO observations. Some selected interesting sources, including two detected at 1 mm, associated with Cha II are discussed in detail and their SEDs presented. The classification of these sources using MIPS data is found to be consistent with previous studies.

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.

We have undertaken a survey of N2H+ and N2D+ towards 31 low-mass starless cores using the IRAM 30m telescope. Our main objective has been to determine the abundance ratio of N2D+ and N2H+ towards the nuclei of these cores and thus to obtain estimates of the degree of deuterium enrichment, a symptom of advanced chemical evolution according to current models. We find that the N(N2D+)/N(N2H+) ratio is larger in more "centrally concentrated cores" with larger peak H2 and N2H+ column density than the sample mean. The deuterium enrichment in starless cores is presently ascribed to depletion of CO in the high density (> 3*10^4 cm-3) core nucleus. To substantiate this picture, we compare our results with observations in dust emission at 1.2 mm and in two transitions of C18O. We find a good correlation between deuterium fractionation and N(C18O)/N(H2) for the nuclei of 14 starless cores. We, thus, identified a set of properties that characterize the most evolved, or "pre-stellar", starless cores. These are: higher N2H+ and N2D+ column densities, higher N(N2D+)/N(N2H+), more pronounced CO depletion, broader N2H+ lines with infall asymmetry, higher central H2 column densities and a more compact density profile than in the average core. We conclude that this combination of properties gives a reliable indication of the evolutionary state of the core. Seven cores in our sample (L1521F, OphD, L429, L694, L183, L1544 and TMC2) show the majority of these features and thus are believed to be closer to forming a protostar than are the other members of our sample. Finally, we note that the subsample of Taurus cores behaves more homogeneously than the total sample, an indication that the external environment could play an important role in the core evolution.

Arcsecond-resolution spectral observations toward the protobinary system IRAS 16293-2422 at 344 and 354 GHz were conducted using the Submillimeter Array. Complex organic molecules such as CH3OH and HCOOCH3 were detected. Together with the rich organic inventory revealed, it clearly indicates the existence of two, rather than one, compact hot molecular cores (smaller than or equal to 400 AU in radius) associated with each of the protobinary components identified by their dust continuum emission in the inner star-forming core.

We present observations of L1014, a dense core in the Cygnus region previously thought to be starless, but data from the Spitzer Space Telescope shows the presence of an embedded source. We propose a model for this source that includes a cold core, heated by the interstellar radiation field, and a low-luminosity internal source. The low luminosity of the internal source suggests a substellar object. If L1014 is representative, other "starless" cores may turn out to harbor central sources.

We have imaged the circumstellar envelope around the binary protostar L1551 IRS 5 in CS (J=7-6) and 343 GHz continuum emission at ~ 3 arcsec resolution using the Submillimeter Array. The continuum emission shows an elongated structure (~ 220 x 100 AU) around the binary perpendicular to the axis of the associated radio jet. The CS emission extends over ~ 400 AU, appears approximately circularly symmetric, and shows a velocity gradient from southeast (blueshifted) to northwest (redshifted). The direction of the velocity gradient is different from that observed in C18O (J=1-0). This may be because rotation is more dominant in the CS envelope than the C18O envelope, in which both infall and rotation exist. The CS emission may be divided into two velocity components: (1) a "high" velocity disk-like structure surrounding the protostar, +-1.0 - 1.5 km s-1 from the systemic velocity, and (2) a "low" velocity structure, located southwest of the protostar, < 1.0 km s -1 from the systemic velocity. The high-velocity component traces warm and dense gas with kinematics consistent with rotation around the protostar. The low-velocity component may arise from dense gas entrained in the outflow. Alternatively, this component may trace infalling and rotating gas in an envelope with a vertical structure.

We observed the pre-stellar core L1521F in dust emission at 1.2mm and in two transitions each of N2H+, N2D+, C18O, and C17O in order to increase the sample of well studied centrally concentrated and chemically evolved starless cores, likely on the verge of star formation, and to determine the initial conditions for low--mass star formation in the Taurus Molecular Cloud. We derived in this object a molecular hydrogen number density n(H2) ~ 10^6 cm-3 and a CO depletion factor, integrated along the line of sight, fD ~ 15 in the central 20'', similar to the pre-stellar core L1544. However, the N(N2D+)/N(N2H+) column density ratio is ~0.1, a factor of about 2 lower than that found in L1544. The observed relation between the deuterium fractionation and the integrated CO depletion factor across the core can be reproduced by chemical models if N2H+ is slightly (factor of ~2 in fractional abundance) depleted in the central 3000 AU. The N2H+ and N2D+ linewidths in the core center are ~0.3 km/s, significantly larger than in other more quiescent Taurus starless cores but similar to those observed in the center of L1544. The kinematical behaviour of L1521F is more complex than seen in L1544, and a model of contraction due to ambipolar diffusion is only marginally consistent with the present data. Other velocity fields, perhaps produced by unresolved substructure, are present. Both chemical and kinematical analyses suggest that L1521F is less evolved than L1544, but, in analogy with L1544, it is approaching the ``critical'' state.

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

Magnetic fields are believed to play an important role in the evolution of molecular clouds, from their large scale structure to dense cores, protostellar envelopes, and protoplanetary disks. How important is unclear, and whether magnetic fields are the dominant force driving star formation at any scale is also unclear. In this review we examine the observational data which address these questions, with particular emphasis on high angular resolution observations. Unfortunately the data do not clarify the situation. It is clear that the fields are important, but to what degree we don't yet know. Observations to date have been limited by the sensitivity of available telescopes and instrumentation. In the future ALMA and the SKA in particular should provide great advances in observational studies of magnetic fields, and we discuss which observations are most desirable when they become available.

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.

BHR 71 is a well isolated Bok globule located at ~200 pc, which harbours a highly collimated bipolar outflow. The outflow is driven by a very young Class 0 protostar with a luminosity of ~9 L_sun. It is one of a very small number that show enhanced abundances of a number of molecular species, notably SiO and CH3OH, due to shock processing of the ambient medium. In this paper the properties of the globule and outflow are discussed.

New AAT near-infrared and SEST 12CO J=2-1 observations are combined with existing ISO mid-infrared and ATCA cm radio continuum observations to examine the protostellar content of the Bok globule BHR 71. Together with observations of Herbig-Haro objects, these data show: (1) Two protostellar sources, IRS1 and IRS2, with a separation of ~17 arcsec (3400 AU) are located within BHR 71. (2) Each protostar is driving its own molecular outflow. The outflow from IRS1 is much larger in extent, is more massive, and dominates the CO emission. (3) Both protostars are associated with Herbig-Haro objects and shock excited 2.122 micron H2 v=1-0S(1) emission, which coincide spatially with their CO outflows. (4) IRS1 is associated with cm continuum emission, with a flat or rising spectrum which is consistent with free-free emission, a signpost of protostellar origin.

We present the results of a new survey of 23 molecular clouds for the Zeeman effect in OH undertaken with the ATNF Parkes 64-m radio telescope and the NRAO Green Bank 43-m radio telescope. The Zeeman effect was clearly detected in the cloud associated with the HII region RCW 38, with a field strength of 38+/-3 micro-Gauss, and possibly detected in a cloud associated with the HII region RCW 57, with a field strength of -203+/-24 micro-Gauss. The remaining 21 measurements give formal upper limits to the magnetic field strength, with typical 1-sigma sensitivities <20 micro-Gauss. For 22 of the molecular clouds we are also able to determine thecolumn density of the gas in which we have made a sensitive search for the Zeeman effect. We combine these results with previous Zeeman studies of 29 molecular clouds, most of which were compiled by Crutcher (1999), for a comparsion of theoretical models with the data. This comparison implies that if the clouds can be modeled as initially spherical with uniform magnetic fields and densities that evolve to their final equilibrium state assuming flux-freezing then the typical cloud is magnetically supercritical, as was found by Crutcher (1999). If the clouds can be modeled as highly flattened sheets threaded by uniform perpendicular fields, then the typical cloud is approximately magnetically critical, in agreement with Shu et al. (1999), but only if the true values of the field for the non-detections are close to the 3-sigma upper limits. If instead these values are significantly lower (for example, similar to the 1-sigma limits), then the typical cloud is generally magnetically supercritical.

We present new observational data for IRAS 09371+1212, the Frosty Leo Nebula, in the form of infrared spectra from 2.2-2.5 microns which reveal photopheric bands of CO. The ^12CO/^13CO band ratio determined is similar to those exhibited by evolved K giant stars, and supports the proposal that the object is highly evolved. The equivalent width of the CO bands implies, however, that the spectral type lies in the range G5III to K0III, somewhat earlier than K7III, as derived from colours and optical spectra. The smaller equivalent widths of the CO bands may reflect a low metal abundance which could fit well with the considerable height of the source above the galactic plane (>0.9kpc).