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.

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.

We observed 146 Galactic clumps in HCN (4-3) and CS (7-6) with the Atacama Submillimeter Telescope Experiment (ASTE) 10-m telescope. A tight linear relationship between star formation rate and gas mass traced by dust continuum emission was found for both Galactic clumps and the high redshift (z>1) star forming galaxies (SFGs), indicating a constant gas depletion time of ~100 Myr for molecular gas in both Galactic clumps and high z SFGs. However, low z galaxies do not follow this relation and seem to have a longer global gas depletion time. The correlations between total infrared luminosities (L_TIR) and molecular line luminosities (L'_mol) of HCN (4-3) and CS (7-6) are tight and sublinear extending down to clumps with LTIR~10^{3} L_sun. These correlations become linear when extended to external galaxies. A bimodal behavior in the LTIR--L'mol correlations was found for clumps with different dust temperature, luminosity-to-mass ratio, and sigmaline/sigmavir. Such bimodal behavior may be due to evolutionary effects. The slopes of LTIR--L'mol correlations become more shallow as clumps evolve. We compared our results with lower J transition lines in wu et al. (2010). The correlations between clump masses and line luminosities are close to linear for low effective excitation density tracers but become sublinear for high effective excitation density tracers for clumps with LTIR larger than LTIR~10^4.5 Lsun. High effective excitation density tracers cannot linearly trace the total clump masses, leading to a sublinear correlations for both Mclump-L'mol and LTIR-L'mol relations.

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 present new measurements of the dust emissivity index, beta, for the high-mass, star-forming OMC 2/3 filament. We combine 160-500 um data from Herschel with long-wavelength observations at 2 mm and fit the spectral energy distributions across a ~ 2 pc long, continuous section of OMC 2/3 at 15000 AU (0.08 pc) resolution. With these data, we measure beta and reconstruct simultaneously the filtered-out large-scale emission at 2 mm. We implement both variable and fixed values of beta, finding that beta = 1.7 - 1.8 provides the best fit across most of OMC 2/3. These beta values are consistent with a similar analysis carried out with filtered Herschel data. Thus, we show that beta values derived from spatial filtered emission maps agree well with those values from unfiltered data at the same resolution. Our results contradict the very low beta values (~ 0.9) previously measured in OMC 2/3 between 1.2 mm and 3.3 mm data, which we attribute to elevated fluxes in the 3.3 mm observations. Therefore, we find no evidence or rapid, extensive dust grain growth in OMC 2/3. Future studies with Herschel data and complementary ground-based long-wavelength data can apply our technique to obtain robust determinations of beta in nearby cold molecular clouds.

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

To test the theoretical understanding that finding bright CO emission depends primarily on dust shielding, we investigate the relationship between CO emission ($I_{\rm CO}$) and the amount of dust (estimated from IR emission and expressed as "$A_V$") across the Large Magellanic Cloud, the Small Magellanic Cloud, and the Milky Way. We show that at our common resolution of 10 pc scales, $I_{\rm CO}$ given a fixed line-of-sight $A_V$ is similar across all three systems despite the difference in metallicity. We find some evidence for a secondary dependence of $I_{\rm CO}$ on radiation field; in the LMC, $I_{\rm CO}$ at a given $A_V$ is smaller in regions of high $T_{\rm dust}$, perhaps because of an increased photodissociating radiation field. We suggest a simple but useful picture in which the CO-to-H$_2$ conversion factor (\xco) depends on two separable factors: (1) the distribution of gas column densities, which maps to an extinction distribution via a dust-to-gas ratio; and (2) the dependence of $I_{\rm CO}$ on $A_V$. Assuming that the probability distribution function (PDF) of local Milky Way clouds is universal, this approach predicts a dependence of \xco\ on $Z$ between $Z^{-1}$ and $Z^{-2}$ above about a third solar metallicity. Below this metallicity, CO emerges from only the high column density parts of the cloud and so depends very sensitively on the adopted PDF and the H$_2$/{\sc Hi} prescription. The PDF of low metallicity clouds is thus of considerable interest and the uncertainty associated with even an ideal prescription for \xco\ at very low metallicity will be large.

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.

The physical and chemical properties of prestellar cores, especially massive ones, are still far from being well understood due to the lack of a large sample. The low dust temperature ($<$14 K) of Planck cold clumps makes them promising candidates for prestellar objects or for sources at the very initial stages of protostellar collapse. We have been conducting a series of observations toward Planck cold clumps (PCCs) with ground-based radio telescopes. In general, when compared with other star forming samples (e.g. infrared dark clouds), PCCs are more quiescent, suggesting that most of them may be in the earliest phase of star formation. However, some PCCs are associated with protostars and molecular outflows, indicating that not all PCCs are in a prestellar phase. We have identified hundreds of starless dense clumps from the mapping survey with the Purple Mountain Observatory (PMO) 13.7-m telescope. Follow-up observations suggest that these dense clumps are ideal targets to search for prestellar objects.

Core accretion models of massive star formation require the existence of stable massive starless cores, but robust observational examples of such objects have proven elusive. We report subarcsecond-resolution SMA 1.3 mm, 1.1 mm, and 0.88 mm and VLA 1.3 cm observations of an excellent massive starless core candidate, G11.92-0.61-MM2, initially identified in the course of studies of GLIMPSE Extended Green Objects (EGOs). Separated by ~7.2" from the nearby MM1 protostellar hot core, MM2 is a strong, compact dust continuum source (submillimeter spectral index alpha=2.6+/-0.1), but is devoid of star formation indicators. In contrast to MM1, MM2 has no masers, no centimeter continuum, and no (sub)millimeter wavelength line emission in ~24 GHz of bandwidth observed with the SMA, including N2H+(3-2), HCO+(3-2), and HCN(3-2). Additionally, there is no evidence for an outflow driven by MM2. The (sub)millimeter spectral energy distribution (SED) of MM2 is best fit with a dust temperature of ~17-19 K and luminosity of ~5-7 L_sun. The combined physical properties of MM2, as inferred from its dust continuum emission, are extreme: M>30 M_sun within a radius<1000 AU, N(H2)>10^25 cm^-2 and n(H2)>10^9 cm^-3. Comparison of the molecular abundance limits derived from our SMA observations with gas-grain chemical models indicates that extremely dense (n(H)>>10^8 cm^-3), cold (<20 K) conditions are required to explain the lack of observed (sub)millimeter line emission, consistent with the dust continuum results. Our data suggest that G11.92-0.61-MM2 is the best candidate for a bonafide massive prestellar core found to date, and a promising target for future, higher-sensitivity observations.

We present a new 3.3 mm continuum map of the OMC-2/3 region. When paired with previously published maps of 1.2mm continuum and NH3-derived temperature, we derive the emissivity spectral index of dust emission in this region, tracking its changes across the filament and cores. We find that the median value of the emissivity spectral index is 0.9, much shallower than previous estimates in other nearby molecular clouds. We find no significant difference between the emissivity spectral index of dust in the OMC-2/3 filament and the starless or protostellar cores. Furthermore, the temperature and emissivity spectral index, beta, are anti-correlated at the 4 sigma level. The low values of the emissivity spectral index found in OMC-2/3 can be explained by the presence of millimeter-sized dust grains in the dense regions of the filaments to which these maps are most sensitive. Alternatively, a shallow dust emissivity spectral index may indicate non-powerlaw spectral energy distributions, significant free-free emission, or anomalous microwave emission. We discuss the possible implications of millimeter-sized dust grains compared to the alternatives.

The Atacama Large mm and sub-mm Array (ALMA) radio observatory is one of the world's largest astronomical projects. After the very successful conclusion of the first observation cycles Early Science Cycles 0 and 1, the ALMA project can report many successes and lessons learned. The science data taken interleaved with commissioning tests for the still continuing addition of new capabilities has already resulted in numerous publications in high-profile journals. The increasing data volume and complexity are challenging but under control. The radio-astronomical data analysis package "Common Astronomy Software Applications" (CASA) has played a crucial role in this effort. This article describes the implementation of the ALMA data quality assurance system, in particular the level 2 which is based on CASA, and the lessons learned.

Users of the Atacama Large Millimeter/submillimeter Array (ALMA) are provided with calibration and imaging products in addition to raw data. In Cycle 0 and Cycle 1, these products are produced by a team of data reduction experts spread across Chile, East Asia, Europe, and North America. This article discusses the lines of communication between the data reducers and ALMA users that enable this model of distributed data reduction. This article also discusses the calibration and imaging scripts that have been provided to ALMA users in Cycles 0 and 1, and what will be different in future Cycles.

Interpreting the nature of starless cores has been a prominent goal in star formation for many years. In order to characterise the evolutionary stages of these objects, we perform synthetic observations of a numerical simulation of a turbulent molecular cloud. We find that nearly all cores that we detect are associated with filaments and eventually form protostars. We conclude that observed starless cores which appear Jeans unstable are only marginally larger than their respective Jeans masses (within a factor of 3). We note single dish observations such as those performed with the JCMT appear to miss significant core structure on small scales due to beam averaging. Finally, we predict that interferometric observations with ALMA Cycle 1 will resolve the important small scale structure, which has so far been missed by mm-wavelength observations.

We present a survey of HCO+ (3-2) observations pointed towards dense cores with previous measurements of N(N2D+)/N(N2H+). Of the 26 cores in this survey, five show the spectroscopic signature of outward motion, nine exhibit neither inward nor outward motion, eleven appear to be infalling, and one is not detected. We compare the degree of deuterium fractionation with infall velocities calculated from the HCO+ spectra and find that those cores with [D]/[H] > 0.1 are more likely to have the signature of inward motions than cores with smaller [D]/[H] ratios. Infall motions are also much more common in cores with masses exceeding their thermal Jeans masses. The fastest infall velocity measured belongs to one of the two protostellar cores in our survey, L1521F, and the observed motions are typically on the order of the sound speed.

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.

We present observations of outflows in the star-forming region NGC 1333 using the Combined Array for Research in Millimeter-Wave Astronomy (CARMA). We combined the 12CO and 13CO (1-0) CARMA mosaics with data from the 14-m Five College Radio Astronomy Observatory (FCRAO) to probe the central, most dense and active region of this protostellar cluster at scales from 5'' to 7' (or 1000 AU to 0.5 pc at a distance of 235 pc). We map and identify 12CO outflows, and along with 13CO data we estimate their mass, momentum and energy. Within the 7'x7' map, the 5'' resolution allows for a detailed study of morphology and kinematics of outflows and outflow candidates, some of which were previously confused with other outflow emission in the region. In total, we identify 22 outflow lobes, as well as 9 dense circumstellar envelopes marked by continuum emission, of which 6 drive outflows. We calculate a total outflow mass, momentum and energy within the mapped region of 6 Msun, 19 Msun km/s, and 7x10^44 erg, respectively. Within this same region, we compare outflow kinematics with turbulence and gravitational energy, and we suggest that outflows are likely important agents for the maintenance of turbulence in this region. In the earliest stages of star formation, outflows do not yet contribute enough energy to totally disrupt the clustered region where most star formation is happening, but have the potential to do so as the protostellar sources evolve. Our results can be used to constrain outflow properties, such as outflow strength, in numerical simulations of outflow-driven turbulence in clusters.

We study the structure and kinematics of nine 0.1 pc-scale cores in Orion with the IRAM 30-m telescope and at higher resolution eight of the cores with CARMA, using CS(2-1) as the main tracer. The single-dish moment zero maps of the starless cores show single structures with central column densities ranging from 7 to 42 times 10^23 cm^-2 and LTE masses from 20 solar masses to 154 solar masses. However, at the higher CARMA resolution (5 arcsec), all of the cores except one fragment into 3 - 5 components. The number of fragments is small compared to that found in some turbulent fragmentation models, although inclusion of magnetic fields may reduce the predicted fragment number and improve the model agreement. This result demonstrates that fragmentation from parsec-scale molecular clouds to sub-parsec cores continues to take place inside the starless cores. The starless cores and their fragments are embedded in larger filamentary structures, which likely played a role in the core formation and fragmentation. Most cores show clear velocity gradients, with magnitudes ranging from 1.7 to 14.3 km/s/pc. We modeled one of them in detail, and found that its spectra are best explained by a converging flow along a filament toward the core center; the gradients in other cores may be modeled similarly. We infer a mass inflow rate of ~ 2 x 10^{-3} Msolar/yr, which is in principle high enough to overcome radiation pressure and allow for massive star formation. However, the core contains multiple fragments, and it is unclear whether the rapid inflow would feed the growth of primarily a single massive star or a cluster of lower mass objects. We conclude that fast, supersonic converging flow along filaments play an important role in massive star and cluster formation.

We present SMA and CARMA continuum and spectral line observations of five dense cores located in the Perseus and Ophiuchus molecular clouds whose masses exceed their thermal Jeans masses. Three of these cores have previously been identified as being starless and two have been classified as being possibly protostellar. We find that one core is certainly protostellar. The other four cores, however, are starless and undetected in both C18O (2-1) and 1.3mm continuum emission. These four starless cores have flat density profiles out to at least 0.006 pc, which is typical for starless cores in general. Density profiles predicted by some collapse models, especially in the early stages of infall, are consistent with our observations. Archival data reveal that these starless cores have significant non-thermal support against collapse, although they may still be unstable.

In this paper, we present the results of CARMA continuum and spectral line observations of the dense core Per-Bolo 45. Although this core has previously been classified as starless, we find evidence for an outflow and conclude that Per-Bolo 45 is actually an embedded, low-luminosity protostar. We discuss the impact of newly discovered, low-luminosity, embedded objects in the Perseus molecular cloud on starless core and protostar lifetimes. We estimate that the starless core lifetime has been overestimated by 4-18% and the Class 0/I protostellar lifetime has been underestimated by 5-20%. Given the relatively large systematic uncertainties involved in these calculations, variations on the order of 10% do not significantly change either core lifetimes or the expected protostellar luminosity function. Finally, we suggest that high resolution (sub)millimeter surveys of known cores lacking near-infrared and mid-infrared emission are necessary to make an accurate census of starless cores.

Determining the initial stellar multiplicity is a challenging problem since protostars are faint and deeply embedded at early times; once formed, multiple protostellar systems may significantly dynamically evolve before they are optically revealed. Interferometers such as CARMA and ALMA make it possible to probe the scales at which turbulent fragmentation occurs in dust continuum emission, potentially constraining early stellar multiplicity. In this Letter we present synthetic observations of starless and protostellar cores undergoing fragmentation on scales of a few thousand AU to produce wide binary systems. We show that interferometric observations of starless cores by CARMA should be predominantly featureless at early stages, although wide protostellar companions should be apparent. The enhanced capabilities of ALMA improve the detection of core morphology so that it may be possible to detect substructure at earlier times. In either case, spatial filtering from interferometry reduces the observed core substructure and often eradicates traces of existing filamentary morphology on scales down to 0.025 pc. However, some missing structure may be recaptured by combining data from the ALMA full science and Atacama compact arrays.

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.

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 results of a pilot multiwavelength survey in the radio continuum (X, Ka and Q bands, i.e., from 3.6 cm to 7 mm) carried out with the EVLA in order to confirm the presence of very large dust grains in dusty disks and torii around the central stars in a small sample of post-AGB objects, as inferred from millimeter and sub-millimeter observations. Supporting mm-wave observations were also obtained with CARMA towards three of our sources. Our EVLA survey has resulted in a robust detection of our most prominent submm emission source, the pre-planetary nebula IRAS 22036+5306, in all three bands, and the disk-prominent post-AGB object, RV Tau, in one band. The observed fluxes are consistent with optically-thin free-free emission, and since they are insignificant compared to their submm/mm fluxes, we conclude that the latter must come from substantial masses of cool, large (mm-sized) grains. We find that the power-law emissivity in the cm-to-submm range for the large grains in IRAS22036 is nu^{beta}, with beta=1-1.3. Furthermore, the value of beta in the 3 to 0.85 mm range for the three disk-prominent post-AGB sources (beta < or = 0.4) is significantly lower than that of IRAS22036, suggesting that the grains in post-AGB objects with circumbinary disks are likely larger than those in the dusty waists of pre-planetary nebulae.

The first hydrostatic core (FHSC) represents a very early phase in the low-mass star formation process, after collapse of the parent core has begun but before a true protostar has formed. This large (few AU), cool (100 K), pressure supported core of molecular hydrogen is expected from theory, but has yet to be observationally verified. Here we present observations of an excellent candidate for the FHSC phase: Per-Bolo 58, a dense core in Perseus that was previously believed to be starless. The 70 micron flux of 65 mJy, from new deep Spitzer MIPS observations, is consistent with that expected for the FHSC. A low signal-to-noise detection at 24 micron leaves open the possibility that Per-Bolo 58 could be a very low luminosity protostar, however. We utilize radiative transfer models to determine the best-fitting FHSC and protostar models to the spectral energy distribution and 2.9 mm visibilities of Per-Bolo 58. The source is consistent with a FHSC with some source of lower opacity through the envelope allowing 24 micron emission to escape; a small outflow cavity or a cavity in the envelope are both possible. While we are unable to rule out the presence of a protostar, if present it would be one of the lowest luminosity protostellar objects yet observed, with an internal luminosity of approximately 0.01 Lsun.

In this paper we present the results of a high resolution (5") CARMA and SZA survey of the 3mm continuum emission from 11 of the brightest (at 1.1mm) starless cores in the Perseus molecular cloud. We detect 2 of the 11 cores, both of which are composed of single structures, and the median 3 sigma upper limit for the non-detections is 0.2 M_sun in a 5" beam. These results are consisent with, and as stringent as, the low detection rate of compact 3mm continuum emission in dense cores in Perseus reported by Olmi et al. (2005). From the non-detection of multiple components in any of the eleven cores we conclude that starless core mass functions derived from bolometer maps at resolutions from 10"-30" (e.g. with MAMBO, SCUBA or Bolocam) are unlikely to be significantly biased by the blending of lower mass cores with small separations. These observations provide additional evidence that the majority of starless cores in Perseus have inner density profiles shallower than r^-2.

In this paper we present a dust emission map of the starless core TMC-1C taken at 2100 microns. Along with maps at 160, 450, 850 and 1200 microns, we study the dust emissivity spectral index from the (sub)millimeter spectral energy distribution, and find that it is close to the typically assumed value of beta = 2. We also map the dust temperature and column density in TMC-1C, and find that at the position of the dust peak (A_V ~ 50), the line-of-sight-averaged temperature is ~7 K. Employing simple Monte Carlo modeling, we show that the data are consistent with a constant value for the emissivity spectral index over the whole map of TMC-1C.

We present polarization observations of the dust continuum emission from IRAS 16293 which is a Class 0 protostar and is known to have at least two components, source A and B. These measurements were conducted by the Submillimeter Array (SMA) at a frequency of ~341 GHz and with high angular resolution (~2-3 arcseconds). We find that the large scale global direction of the field, which is perpendicular to the observed polarization, appears to be along the dust ridge where the emission peaks. On smaller scales we find that the field structure is significantly different for the two components of the binary. The first component, source A, shows a magnetic field structure which is "hourglass" shaped as predicted from theoretical models of low mass star formation in the presence of strong magnetic fields. However, the other component, source B, shows a relatively ordered magnetic field with no evidence of any deformation. There is an observed decrease in polarization towards the center. Our calculations show that in IRAS 16293 the magnetic energy is stronger than the turbulent energy but is approximately similar to the centrifugal energy. Our results provide additional evidence to show that the two protostars appear to be in different stages during their evolution.

We present 1" resolution CARMA observations of the 3mm continuum and 95 GHz methanol masers toward 14 candidate high mass protostellar objects (HMPOs). Dust continuum emission is detected toward seven HMPOs, and methanol masers toward 5 sources. The 3mm continuum sources have diameters < 2x10^4 AU, masses between 21 and 1200 M_sun, and volume densities > 10^8 cm^-3. Most of the 3mm continuum sources are spatially coincident with compact HII regions and/or water masers, and are presumed to be formation sites of massive stars. A strong correlation exists between the presence of 3mm continuum emission, 22 GHz water masers, and 95 GHz methanol masers. However, no 3mm continuum emission is detected toward ultracompact HII regions lacking maser emission. These results are consistent with the hypothesis that 22 GHz water masers and methanol masers are signposts of an early phase in the evolution of an HMPO before an expanding HII region destroys the accretion disk.

We investigate the effect of line of sight temperature variations and noise on two commonly used methods to determine dust properties from dust continuum observations of dense cores. One method employs a direct fit to a modified blackbody SED; the other involves a comparison of flux ratios to an analytical prediction. Fitting fluxes near the SED peak produces inaccurate temperature and dust spectral index estimates due to the line of sight temperature (and density) variations. Longer wavelength fluxes in the Rayleigh-Jeans part of the spectrum (>~ 600 micron for typical cores) may more accurately recover the spectral index, but both methods are very sensitive to noise. The temperature estimate approaches the density weighted temperature, or "column temperature," of the source as short wavelength fluxes are excluded. An inverse temperature - spectral index correlation naturally results from SED fitting, due to the inaccurate isothermal assumption, as well as noise uncertainties. We show that above some "threshold" temperature, the temperatures estimated through the flux ratio method can be highly inaccurate. In general, observations with widely separated wavelengths, and including shorter wavelengths, result in higher threshold temperatures; such observations thus allow for more accurate temperature estimates of sources with temperatures less than the threshold temperature. When only three fluxes are available, a constrained fit, where the spectral index is fixed, produces less scatter in the temperature estimate when compared to the estimate from the flux ratio method.

We investigate how uncertainties in flux measurements affect the results from modified blackbody SED fits. We show that an inverse correlation between the dust temperature T and spectral index (beta) naturally arises from least squares fits due to the uncertainties, even for sources with a single T and beta. Fitting SEDs to noisy fluxes solely in the Rayleigh-Jeans regime produces unreliable T and beta estimates. Thus, for long wavelength observations (lambda >~ 200 micron), or for warm sources (T >~ 60 K), it becomes difficult to distinguish sources with different temperatures. We assess the role of noise in recent observational results that indicate an inverse and continuously varying T - beta relation. Though an inverse and continuous T - beta correlation may be a physical property of dust in the ISM, we find that the observed inverse correlation may be primarily due to noise.

We use the COMPLETE Survey's observations of the Perseus star-forming region to assess and intercompare three methods for measuring column density in molecular clouds: extinction mapping (NIR); thermal emission mapping (FIR); and mapping the intensity of CO isotopologues. The structures shown by all three tracers are morphologically similar, but important differences exist. Dust-based measures give similar, log-normal, distributions for the full Perseus region, once careful calibration corrections are made. We also compare dust- and gas-based column density distributions for physically-meaningful sub-regions of Perseus, and we find significant variations in the distributions for those regions. Even though we have used 12CO data to estimate excitation temperatures, and we have corrected for opacity, the 13CO maps seem unable to give column distributions that consistently resemble those from dust measures. We have edited out the effects of the shell around the B-star HD 278942. In that shell's interior and in the parts where it overlaps the molecular cloud, there appears to be a dearth of 13CO, likely due either to 13CO not yet having had time to form in this young structure, and/or destruction of 13CO in the molecular cloud. We conclude that the use of either dust or gas measures of column density without extreme attention to calibration and artifacts is more perilous than even experts might normally admit. And, the use of 13CO to trace total column density in detail, even after proper calibration, is unavoidably limited in utility due to threshold, depletion, and opacity effects. If one's main aim is to map column density, then dust extinction seems the best probe. Linear fits amongst column density tracers are given, quantifying the inherent uncertainties in using one tracer (when compared with others). [abridged]

In this paper we study the determinants of starless core temperatures in the Perseus molecular cloud. We use NH3 (1,1) and (2,2) observations to derive core temperatures (T_kin) and data from the COMPLETE Survey of Star Forming Regions and the c2d Spitzer Legacy Survey for observations of the other core and molecular cloud properties. The kinetic temperature distribution probed by NH3 is in the fairly narrow range of 9 - 15 K. We find that cores within the clusters IC348 and NGC1333 are significantly warmer than "field" starless cores, and T_kin is higher within regions of larger extinction-derived column density. Starless cores in the field are warmer when they are closer to class O/I protostars, but this effect is not seen for those cores in clusters. For field starless cores, T_kin is higher in regions in which the 13CO linewidth and the 1.1mm flux from the core are larger, and T_kin is lower when the the peak column density within the core and average volume density of the core are larger. There is no correlation between T_kin and 13CO linewidth, 1.1mm flux, density or peak column density for those cores in clusters. The temperature of the cloud material along the line of sight to the core, as measured by CO or far-infrared emission from dust, is positively correlated with core temperature when considering the collection of cores in the field and in clusters, but this effect is not apparent when the two subsamples of cores are considered separately.

Using far-infrared emission maps taken by IRAS and Spitzer and a near-infrared extinction map derived from 2MASS data, we have made dust temperature and column density maps of the Perseus molecular cloud. We show that the emission from transiently heated very small grains and the big grain dust emissivity vary as a function of extinction and dust temperature, with higher dust emissivities for colder grains. This variable emissivity can not be explained by temperature gradients along the line of sight or by noise in the emission maps, but is consistent with grain growth in the higher density and lower temperature regions. By accounting for the variations in the dust emissivity and VSG emission, we are able to map the temperature and column density of a nearby molecular cloud with better accuracy than has previously been possible.

We have mapped the starless core TMC-1C in a variety of molecular lines with the IRAM 30m telescope. High density tracers show clear signs of self-absorption and sub-sonic infall asymmetries are present in N2H+ (1-0) and DCO+ (2-1) lines. The inward velocity profile in N2H+ (1-0) is extended over a region of about 7,000 AU in radius around the dust continuum peak, which is the most extended ``infalling'' region observed in a starless core with this tracer. The kinetic temperature (~12 K) measured from C17O and C18O suggests that their emission comes from a shell outside the colder interior traced by the mm continuum dust. The C18O (2-1) excitation temperature drops from 12 K to ~10 K away from the center. This is consistent with a volume density drop of the gas traced by the C18O lines, from ~4x10^4 cm^-3 towards the dust peak to ~6x10^3 cm^-3 at a projected distance from the dust peak of 80" (or 11,000 AU). The column density implied by the gas and dust show similar N2H+ and CO depletion factors (f_D < 6). This can be explained with a simple scenario in which: (i) the TMC-1C core is embedded in a relatively dense environment (H2 ~10^4 cm^-3), where CO is mostly in the gas phase and the N2H+ abundance had time to reach equilibrium values; (ii) the surrounding material (rich in CO and N2H+) is accreting onto the dense core nucleus; (iii) TMC-1C is older than 3x10^5 yr, to account for the observed abundance of N2H+ across the core (~10^-10 w.r.t. H2); and (iv) the core nucleus is either much younger (~10^4 yr) or ``undepleted'' material from the surrounding envelope has fallen towards it in the past 10,000 yr.

We have mapped the central 10'x10' of the dense core TMC-1C at 450, 850 and 1200 microns using SCUBA on the James Clerk Maxwell Telescope and MAMBO on the IRAM 30m telescope. We show that although one can, in principle, use images at these wavelengths to map the emissivity spectral index, temperature and column density independently, noise and calibration errors would have to be less than 2% to accurately derive these three quantities from a set of three emission maps. Because our data are not this free of errors, we use our emission maps to fit the dust temperature and column density assuming a constant value of the emissivity spectral index and explore the effects of noise on the derived physical parameters. We find that the derived extinction values for TMC-1C are large for a starless core (80 mag Av), and the derived temperatures are low (6 K) in the densest regions of the core, using our derived value of beta = 1.8.

We present an overview of data available for the Ophiuchus and Perseus molecular clouds from ``Phase I'' of the COMPLETE Survey of Star-Forming Regions. This survey provides a range of data complementary to the Spitzer Legacy Program ``From Molecular Cores to Planet Forming Disks.'' Phase I includes: Extinction maps derived from 2MASS near-infrared data using the NICER algorithm; extinction and temperature maps derived from IRAS 60 and 100um emission; HI maps of atomic gas; 12CO and 13CO maps of molecular gas; and submillimetre continuum images of emission from dust in dense cores. Not unexpectedly, the morphology of the regions appears quite different depending on the column-density tracer which is used, with IRAS tracing mainly warmer dust and CO being biased by chemical, excitation and optical depth effects. Histograms of column-density distribution are presented, showing that extinction as derived from 2MASS/NICER gives the closest match to a log-normal distribution as is predicted by numerical simulations. All the data presented in this paper, and links to more detailed publications on their implications are publically available at the COMPLETE website.

We have used a numerical simulation of a turbulent cloud to synthesize maps of the thermal emission from dust at a variety of far-IR and sub-mm wavelengths. The average column density and external radiation field in the simulation is well matched to clouds such as Perseus and Ophiuchus. We use pairs of single-wavelength emission maps to derive the dust color temperature and column density, and we compare the derived column densities with the true column density. We demonstrate that longer wavelength emission maps yield less biased estimates of column density than maps made towards the peak of the dust emission spectrum. We compare the scatter in the derived column density with the observed scatter in Perseus and Ophiuchus. We find that while in Perseus all of the observed scatter in the emission-derived versus the extinction-derived column density can be attributed to the flawed assumption of isothermal dust along each line of sight, in Ophiuchus there is additional scatter above what can be explained by the isothermal assumption. Our results imply that variations in dust emission properties within a molecular cloud are not necessarily a major source of uncertainty in column density measurements.

The Perseus molecular cloud complex is a ~30pc long chain of molecular clouds most well-known for the two star-forming clusters NGC1333 and IC348 and the well-studied outflow source in B5. However, when studied at mid- to far-infrared wavelengths the region is dominated by a ~10pc diameter shell of warm dust, likely generated by an HII region caused by the early B-star HD278942. Using a revised calibration technique the COMPLETE team has produced high-sensitivity temperature and column-density maps of the Perseus region from IRAS Sky Survey Atlas (ISSA) 60 and 100um data. In this paper, we combine the ISSA based dust-emission maps with other observations collected as part of the COMPLETE Survey, along with archival H-alpha and MSX observations. Molecular line observations from FCRAO and extinction maps constructed by applying the NICER method to the 2MASS catalog provide independent estimates of the ``true'' column-density of the shell. H-alpha emission in the region of the shell confirms that it is most likely an HII region located behind the cloud complex, and 8um data from MSX indicates that the shell may be interacting with the cloud. Finally, the two polarisation components previously seen towards background stars in the region can be explained by the association of the stronger component with the shell. If confirmed, this would be the first observation of a parsec-scale swept-up magnetic field.

We have created new dust temperature and column density maps of Perseus, Ophiuchus and Serpens using 60 and 100 micron data from the IRIS recalibration of IRAS data. We describe an optimized method for finding the dust temperature, emissivity spectral index, and optical depth using optical and NIR extinction maps. The creation of these temperature and extinction maps (covering tens of square degrees of molecular clouds) is one of the first results from the ongoing COordinated Molecular Probe Line Extinction Thermal Emission (COMPLETE) Survey of Star Forming regions. However, while the extinctions derived from the IRIS emission maps are globally accurate, we warn that FIR emission is not a good proxy for extinction on the scale of one pixel (5').

We have mapped the central 10'x10' of the dense core TMC-1C at 450 and 850 microns using SCUBA on the James Clerk Maxwell Telescope. The unusually high quality of the 450 micron map allows us to make a detailed analysis of the temperature and column density profiles of the core. We find that the dust temperature at the center of TMC-1C is 7 K, rising to 11 K at the edges. We discuss the possibility and effects of a variable emissivity spectral index on the derived mass profile. The low dust temperature of TMC-1C results in a high derived mass for the core, significantly larger than the virial mass estimated from the linewidth of the N2H+ (1-0) transition. This result is valid within a wide range of dust properties and ellipticities of the core. The N2H+ (1-0) spectra, taken with the IRAM 30m telescope, show signs of self-absorption, which provide evidence of sub-sonic infall motions. The derived density profile and infall velocity is compared to the predictions of several popular star formation models, and the Bonnor-Ebert model isthe best fit analytic model.