We present ALMA maps of the starless molecular cloud core Ophiuchus/H-MM1 in the lines of deuterated ammonia (ortho-NH2D), methanol (CH3OH), and sulphur monoxide (SO). The dense core is seen in NH2D emission, whereas the CH3OH and SO distributions form a halo surrounding the core. Because methanol is formed on grain surfaces, its emission highlights regions where desorption from grains is particularly efficient. Methanol and sulphur monoxide are most abundant in a narrow zone that follows the eastern side of the core. This side is sheltered from the stronger external radiation field coming from the west. We show that photodissociation on the illuminated side can give rise to an asymmetric methanol distribution, but that the stark contrast observed in H-MM1 is hard to explain without assuming enhanced desorption on the shaded side. The region of the brightest emission has a wavy structure that rolls up at one end. This is the signature of Kelvin-Helmholtz instability occurring in sheared flows. We suggest that in this zone, methanol and sulphur are released as a result of grain-grain collisions induced by shear vorticity.

We report the results of a pilot program for a Green Bank Telescope (GBT) MUSTANG Galactic Plane survey at 3 mm (90 GHz), MGPS90. The survey achieves a typical $1\sigma$ depth of $1-2$ mJy beam$^{-1}$ with a 9" beam. We describe the survey parameters, quality assessment process, cataloging, and comparison with other data sets. We have identified 709 sources over seven observed fields selecting some of the most prominent millimeter-bright regions between $0\deg < \ell < 50\deg$ (total area $\approx 7.5 \deg^2$). The majority of these sources have counterparts at other wavelengths. By applying flux selection criteria to these sources, we successfully recovered several known hypercompact HII (HCHII) regions, but did not confirm any new ones. We identify 126 sources that have mm-wavelength counterparts but do not have cm-wavelength counterparts and are therefore candidate HCHII regions; of these, 10 are morphologically compact and are strong candidates for new HCHII regions. Given the limited number of candidates in the extended area in this survey compared to the relatively large numbers seen in protoclusters W51 and W49, it appears that most HCHII regions exist within dense protoclusters. Comparing the counts of HCHII to ultracompact HII (UCHII) regions, we infer the HCHII region lifetime is 16-46% that of the UCHII region lifetime. We additionally separated the 3 mm emission into dust and free-free emission by comparing with archival 870 $\mu$m and 20 cm data. In the selected pilot fields, most ($\gtrsim80$%) of the 3 mm emission comes from plasma, either through free-free or synchrotron emission.

The role played by magnetic field during star formation is an important topic in astrophysics. We investigate the correlation between the orientation of star-forming cores (as defined by the core major axes) and ambient magnetic field directions in 1) a 3D MHD simulation, 2) synthetic observations generated from the simulation at different viewing angles, and 3) observations of nearby molecular clouds. We find that the results on relative alignment between cores and background magnetic field in synthetic observations slightly disagree with those measured in fully 3D simulation data, which is partly because cores identified in projected 2D maps tend to coexist within filamentary structures, while 3D cores are generally more rounded. In addition, we examine the progression of magnetic field from pc- to core-scale in the simulation, which is consistent with the anisotropic core formation model that gas preferably flow along the magnetic field toward dense cores. When comparing the observed cores identified from the GBT Ammonia Survey (GAS) and Planck polarization-inferred magnetic field orientations, we find that the relative core-field alignment has a regional dependence among different clouds. More specifically, we find that dense cores in the Taurus molecular cloud tend to align perpendicular to the background magnetic field, while those in Perseus and Ophiuchus tend to have random (Perseus) or slightly parallel (Ophiuchus) orientations with respect to the field. We argue that this feature of relative core-field orientation could be used to probe the relative significance of the magnetic field within the cloud.

Recent observations of global velocity gradients across and along molecular filaments have been interpreted as signs of gas accreting onto and along these filaments, potentially feeding star-forming cores and proto-clusters. The behavior of velocity gradients in filaments, however, has not been studied in detail, particularly on small scales (< 0.1 pc). In this paper, we present MUFASA, an efficient, robust, and automatic method to fit ammonia lines with multiple velocity components, generalizable to other molecular species. We also present CRISPy, a Python package to identify filament spines in 3D images (e.g., position-position-velocity cubes), along with a complementary technique to sort fitted velocity components into velocity-coherent filaments. In NGC 1333, we find a wealth of velocity gradient structures on a beam-resolved scale of ~0.05 pc. Interestingly, these local velocity gradients are not randomly oriented with respect to filament spines and their perpendicular, i.e., radial, component decreases in magnitude towards the spine for many filaments. Together with remarkably constant velocity gradients on larger scales along many filaments, these results suggest a scenario in which gas falling onto filaments is progressively damped and redirected to flow along these filaments.

The detection of complex organic molecules (COMs) toward dense, collapsing prestellar cores has sparked interest in the fields of astrochemistry and astrobiology, yet the mechanisms for COM formation are still debated. It was originally believed that COMs initially form in ices which are then irradiated by UV radiation from the surrounding interstellar radiation field as well as forming protostars and subsequently photodesorbed into the gas-phase. However, starless and prestellar cores do not have internal protostars to heat-up and sublimate the ices. Alternative models using chemical energy have been developed to explain the desorption of COMs, yet in order to test these models robust measurements of COM abundances are needed toward representative samples of cores. We've conducted a large-sample survey of 31 starless and prestellar cores in the Taurus Molecular Cloud, detecting methanol (CH$_3$OH) in 100$\%$ of the cores targeted and acetaldehyde (CH$_3$CHO) in 70$\%$. At least two transition lines of each molecule were measured, allowing us to place tight constraints on excitation temperature, column density and abundance. Additional mapping of methanol revealed extended emission, detected down to A$_\mathrm{V}$ as low as $\sim$ 3 mag. We find complex organic molecules are detectable in the gas-phase and are being formed early, at least hundreds of thousands of years prior to star and planet formation. The precursor molecule, CH$_3$OH, may be chemically linked to the more complex CH$_3$CHO, however higher spatial resolution maps are needed to further test chemical models.

We present an analysis of the internal velocity structures of the newly identified sub-0.1 pc coherent structures, droplets, in L1688 and B18. By fitting 2D linear velocity fields to the observed maps of velocity centroids, we determine the magnitudes of linear velocity gradients and examine the potential rotational motions that could lead to the observed velocity gradients. The results show that the droplets follow the same power-law relation between the velocity gradient and size found for larger-scale dense cores. Assuming that rotational motion giving rise to the observed velocity gradient in each core is a solid-body rotation of a rotating body with a uniform density, we derive the "net rotational motions" of the droplets. We find a ratio between rotational and gravitational energies, $\beta$, of $\sim 0.046$ for the droplets, and when including both droplets and larger-scale dense cores, we find $\beta \sim 0.039$. We then examine the alignment between the velocity gradient and the major axis of each droplet, using methods adapted from the histogram of relative orientations (HRO) introduced by Soler et al. (2013). We find no definitive correlation between the directions of velocity gradients and the elongations of the cores. Lastly, we discuss physical processes other than rotation that may give rise to the observed velocity field.

(Abridged) The initial physical conditions of high-mass stars and protoclusters remain poorly characterized. To this end we present the first targeted ALMA 1.3mm continuum and spectral line survey towards high-mass starless clump candidates, selecting a sample of 12 of the most massive candidates ($400-4000\, M_\odot$) within 5 kpc. The joint 12+7m array maps have a high spatial resolution of $\sim 3000\, \mathrm{au}$ ($\sim 0.8^{\prime\prime}$) and have point source mass-completeness down to $\sim 0.3\, M_\odot$ at $6\sigma$ (or $1\sigma$ column density sensitivity of $1.1\times10^{22}\, \mathrm{cm^{-2}}$). We discover previously undetected signposts of low-luminosity star formation from CO (2-1) and SiO (5-4) bipolar outflows and other signatures towards 11 out of 12 clumps, showing that current MIR/FIR Galactic Plane surveys are incomplete to low- and intermediate-mass protostars ($\lesssim 50\, L_\odot$). We compare a subset of the observed cores with a suite of radiative transfer models of starless cores. We find a high-mass starless core candidate with a model-derived mass consistent with $29^{52}_{15}\, M_\odot$ when integrated over size scales of $2\times10^4\, \mathrm{au}$. Unresolved cores are poorly fit by starless core models, supporting the interpretation that they are protostellar even without detection of outflows. Substantial fragmentation is observed towards 10 out of 12 clumps. We extract sources from the maps using a dendrogram to study the characteristic fragmentation length scale. Nearest neighbor separations when corrected for projection are consistent with being equal to the clump average thermal Jeans length. Our findings support a hierarchical fragmentation process, where the highest density regions are not strongly supported against thermal gravitational fragmentation by turbulence or magnetic fields.

We present the observation and analysis of newly discovered coherent structures in the L1688 region of Ophiuchus and the B18 region of Taurus. Using data from the Green Bank Ammonia Survey (GAS), we identify regions of high density and near-constant, almost-thermal, velocity dispersion. Eighteen coherent structures are revealed, twelve in L1688 and six in B18, each of which shows a sharp "transition to coherence" in velocity dispersion around its periphery. The identification of these structures provides a chance to study the coherent structures in molecular clouds statistically. The identified coherent structures have a typical radius of 0.04 pc and a typical mass of 0.4 Msun, generally smaller than previously known coherent cores identified by Goodman et al. (1998), Caselli et al. (2002), and Pineda et al. (2010). We call these structures "droplets." We find that unlike previously known coherent cores, these structures are not virially bound by self-gravity and are instead predominantly confined by ambient pressure. The droplets have density profiles shallower than a critical Bonnor-Ebert sphere, and they have a velocity (VLSR) distribution consistent with the dense gas motions traced by NH3 emission. These results point to a potential formation mechanism through pressure compression and turbulent processes in the dense gas. We present a comparison with a magnetohydrodynamic simulation of a star-forming region, and we speculate on the relationship of droplets with larger, gravitationally bound coherent cores, as well as on the role that droplets and other coherent structures play in the star formation process.

We perform a virial analysis of starless dense cores in three nearby star-forming regions : L1688 in Ophiuchus, NGC 1333 in Perseus, and B18 in Taurus. Our analysis takes advantage of comprehensive kinematic information for the dense gas in all of these regions made publicly available through the Green Bank Ammonia Survey Data Release 1, which used to estimate internal support against collapse. We combine this information with ancillary data used to estimate other important properties of the cores, including continuum data from the James Clerk Maxwell Telescope Gould Belt Survey for core identification, core masses, and core sizes. Additionally, we used \textit{Planck} and \textit{Herschel}-based column density maps for external cloud weight pressure, and Five College Radio Astronomy Observatory $^{13}$CO observations for external turbulent pressure. Our self-consistent analysis suggests that many dense cores in all three star-forming regions are not bound by gravity alone, but rather require additional pressure confinement to remain bound. Unlike a recent, similar study in Orion~A, we find that turbulent pressure represents a significant portion of the external pressure budget. Our broad conclusion emphasizing the importance of pressure confinement in dense core evolution, however, agrees with earlier work.

We present a comparative study of physical properties derived from gas and dust emission in a sample of 1068 dense Galactic clumps. The sources are selected from the crossmatch of the Hi-GAL survey with 16 catalogues of NH$_3$ line emission in its lowest inversion (1,1) and (2,2) transitions. The sample covers a large range in masses and bolometric luminosities, with surface densities above $\Sigma=0.1$ g cm$^{-2}$ and with low virial parameters $\alpha<1$. The comparison between dust and gas properties shows an overall agreement between $T_{\textit{kin}}$ and $T_{\textit{dust}}$ at volumetric densities $n\gtrsim1.2\times10^{4}$ cm$^{-3}$, and a median fractional abundance $\chi$(NH$_3$)$=1.46\times10^{-8}$. While the protostellar clumps in the sample have small differences between $T_{\textit{kin}}$ and $T_{\textit{dust}}$, prestellar clumps have a median ratio $T_{\textit{kin}}/T_{\textit{dust}}=1.24$, suggesting that these sources are thermally decoupled. A correlation is found between the evolutionary tracer $L/M$ and the parameters $T_{\textit{kin}}/T_{\textit{dust}}$ and $\chi$(NH$_3$) in prestellar sources and protostellar clumps with $L/M<1$ L$_\odot$ M$_\odot^{-1}$. In addition, a weak correlation is found between non-thermal velocity dispersion and the $L/M$ parameter, possibly indicating an increase of turbulence with protostellar evolution in the interior of clumps. Finally, different processes are discussed to explain the differences between gas and dust temperatures in prestellar candidates, and the origin of non-thermal motions observed in the clumps.

We present deep CCS and HC$_7$N observations of the L1495-B218 filaments in the Taurus molecular cloud obtained using the K-band focal plane array on the 100m Green Bank Telescope. We observed the L1495-B218 filaments in CCS $J_N$ = 2$_1$$-$1$_0$ and HC$_7$N $J$ = 21$-$20 with a spectral resolution of 0.038 km s$^{-1}$ and an angular resolution of 31$''$. We observed strong CCS emission in both evolved and young regions and weak emission in two evolved regions. HC$_7$N emission is observed only in L1495A-N and L1521D. We find that CCS and HC$_7$N intensity peaks do not coincide with NH$_3$ or dust continuum intensity peaks. We also find that the fractional abundance of CCS does not show a clear correlation with the dynamical evolutionary stage of dense cores. Our findings and chemical modeling indicate that the fractional abundances of CCS and HC$_7$N are sensitive to the initial gas-phase C/O ratio, and they are good tracers of young condensed gas only when the initial C/O is close to solar value. Kinematic analysis using multiple lines including NH$_3$, HC$_7$N, CCS, CO, HCN, \& HCO$^+$ suggests that there may be three different star formation modes in the L1495-B218 filaments. At the hub of the filaments, L1495A/B7N has formed a stellar cluster with large-scale inward flows (fast mode), while L1521D, a core embedded in a filament, is slowly contracting due to its self-gravity (slow mode). There is also one isolated core that appears to be marginally stable and may undergo quasi-static evolution (isolated mode).

Post-starburst or "E+A" galaxies are rapidly transitioning from star-forming to quiescence. While the current star formation rate of post-starbursts is already at the level of early type galaxies, we recently discovered that many have large CO-traced molecular gas reservoirs consistent with normal star forming galaxies. These observations raise the question of why these galaxies have such low star formation rates. Here we present an ALMA search for the denser gas traced by HCN (1--0) and HCO+ (1--0) in two CO-luminous, quiescent post-starburst galaxies. Intriguingly, we fail to detect either molecule. The upper limits are consistent with the low star formation rates and with early-type galaxies. The HCN/CO luminosity ratio upper limits are low compared to star-forming and even many early type galaxies. This implied low dense gas mass fraction explains the low star formation rates relative to the CO-traced molecular gas and suggests the state of the gas in post-starburst galaxies is unusual, with some mechanism inhibiting its collapse to denser states. We conclude that post-starbursts galaxies are now quiescent because little dense gas is available, in contrast to the significant CO-traced lower density gas reservoirs that still remain.

Recent Galactic plane surveys of dust continuum emission at long wavelengths have identified a population of dense, massive clumps with no evidence for on-going star formation. These massive starless clump candidates are excellent sites to search for the initial phases of massive star formation before the feedback from massive star formation effects the clump. In this study, we search for the spectroscopic signature of inflowing gas toward starless clumps, some of which are massive enough to form a massive star. We observed 101 starless clump candidates identified in the Bolocam Galactic Plane Survey (BGPS) in HCO+ J = 1-0 using the 12m Arizona Radio Observatory telescope. We find a small blue excess of E = (Nblue - Nred)/Ntotal = 0.03 for the complete survey. We identified 6 clumps that are good candidates for inflow motion and used a radiative transfer model to calculate mass inflow rates that range from 500 - 2000 M /Myr. If the observed line profiles are indeed due to large-scale inflow motions, then these clumps will typically double their mass on a free fall time. Our survey finds that massive BGPS starless clump candidates with inflow signatures in HCO+ J = 1-0 are rare throughout our Galaxy.

We use Green Bank Ammonia Survey observations of NH$_3$ (1,1) and (2,2) emission with 32'' FWHM resolution from a ~ 10 pc$^{2}$ portion of the Cepheus-L1251 molecular cloud to identify hierarchical dense gas structures. Our dendrogram analysis of the NH$_3$ data results in 22 top-level structures, which reside within 13 lower-level, parent structures. The structures are compact (0.01 pc $\lesssim R_{eff} \lesssim$ 0.1 pc) and are spatially correlated with the highest H$_2$ column density portions of the cloud. We also compare the ammonia data to a catalog of dense cores identified by higher-resolution (18.2'' FWHM) Herschel Space Observatory observations of dust continuum emission from Cepheus-L1251. Maps of kinetic gas temperature, velocity dispersion, and NH$_3$ column density, derived from detailed modeling of the NH$_3$ data, are used to investigate the stability and chemistry of the ammonia-identified and Herschel-identified structures. We show that the dust and dense gas in the structures have similar temperatures, with median $T_{dust}$ and $T_K$ measurements of 11.7 $\pm$ 1.1 K and 10.3 $\pm$ 2.0 K, respectively. Based on a virial analysis, we find that the ammonia-identified structures are gravitationally dominated, yet may be in or near a state of virial equilibrium. Meanwhile, the majority of the Herschel-identified dense cores appear to be not bound by their own gravity and instead confined by external pressure. CCS $(2_0-1_0)$ and HC$_5$N $(9-8)$ emission from the region reveal broader line widths and centroid velocity offsets when compared to the NH$_3$ (1,1) emission in some cases, likely due to these carbon-based molecules tracing the turbulent outer layers of the dense cores.

Density profiles of isolated cores derived from thermal dust continuum emission rely on models of dust properties, such as mass opacity, which are poorly constrained. With complementary measures from near-infrared extinction maps, we can assess the reliability of commonly-used dust models. In this work, we compare Herschel-derived maps of the optical depth with equivalent maps derived from CFHT WIRCAM near-infrared observations for three isolated cores: CB68, L429, and L1552. We assess the dust opacities provided from four models: OH1a, OH5a, Orm1, and Orm4. Although the consistency of the models differs between the three sources, the results suggest that the optical properties of dust in the envelopes of the cores are best described by either silicate and bare graphite grains (e.g., Orm1) or carbonaceous grains with some coagulation and either thin or no ice mantles (e.g., OH5a). None of the models, however, individually produced the most consistent optical depth maps for every source. The results suggest that either the dust in the cores is not well described by any one dust property model, the application of the dust models cannot be extended beyond the very center of the cores, or more complex SED fitting functions are necessary.

We use gas temperature and velocity dispersion data from the Green Bank Ammonia Survey and core masses and sizes from the James Clerk Maxwell Telescope Gould Belt Survey to estimate the virial states of dense cores within the Orion A molecular cloud. Surprisingly, we find that almost none of the dense cores are sufficiently massive to be bound when considering only the balance between self-gravity and the thermal and non-thermal motions present in the dense gas. Including the additional pressure binding imposed by the weight of the ambient molecular cloud material and additional smaller pressure terms, however, suggests that most of the dense cores are pressure confined.

We present an overview of the first data release (DR1) and first-look science from the Green Bank Ammonia Survey (GAS). GAS is a Large Program at the Green Bank Telescope to map all Gould Belt star-forming regions with $A_V \gtrsim 7$ mag visible from the northern hemisphere in emission from NH$_3$ and other key molecular tracers. This first release includes the data for four regions in Gould Belt clouds: B18 in Taurus, NGC 1333 in Perseus, L1688 in Ophiuchus, and Orion A North in Orion. We compare the NH$_3$ emission to dust continuum emission from Herschel, and find that the two tracers correspond closely. NH$_3$ is present in over 60\% of lines-of-sight with $A_V \gtrsim 7$ mag in three of the four DR1 regions, in agreement with expectations from previous observations. The sole exception is B18, where NH$_3$ is detected toward ~ 40\% of lines-of-sight with $A_V \gtrsim 7$ mag. Moreover, we find that the NH$_3$ emission is generally extended beyond the typical 0.1 pc length scales of dense cores. We produce maps of the gas kinematics, temperature, and NH$_3$ column densities through forward modeling of the hyperfine structure of the NH$_3$ (1,1) and (2,2) lines. We show that the NH$_3$ velocity dispersion, ${\sigma}_v$, and gas kinetic temperature, $T_K$, vary systematically between the regions included in this release, with an increase in both the mean value and spread of ${\sigma}_v$ and $T_K$ with increasing star formation activity. The data presented in this paper are publicly available.

The calculation of the molecular column density from molecular spectral (rotational or ro-vibrational) transition measurements is one of the most basic quantities derived from molecular spectroscopy. Starting from first principles where we describe the basic physics behind the radiative and collisional excitation of molecules and the radiative transfer of their emission, we derive a general expression for the molecular column density. As the calculation of the molecular column density involves a knowledge of the molecular energy level degeneracies, rotational partition functions, dipole moment matrix elements, and line strengths, we include generalized derivations of these molecule-specific quantities. Given that approximations to the column density equation are often useful, we explore the optically thin, optically thick, and low-frequency limits to our derived general molecular column density relation. We also evaluate the limitations of the common assumption that the molecular excitation temperature is constant, and address the distinction between beam- and source-averaged column densities. We conclude our discussion of the molecular column density with worked examples for C$^{18}$O, C$^{17}$O, N$_2$H$^+$, NH$_3$, and H$_2$CO. Ancillary information on some subtleties involving line profile functions, conversion between integrated flux and brightness temperature, the calculation of the uncertainty associated with an integrated intensity measurement, the calculation of spectral line optical depth using hyperfine or isotopologue measurements, the calculation of the kinetic temperature from a symmetric molecule excitation temperature measurement, and relative hyperfine intensity calculations for NH$_3$ are presented in appendices. The intent of this document is to provide a reference for researchers studying astrophysical molecular spectroscopic measurements.

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.

Ammonia and its deuterated isotopologues probe physical conditions in dense molecular cloud cores. With the aim of testing the current understanding of the spin-state chemistry of these molecules, we observed spectral lines of NH3, NH2D, NHD2, ND3, and N2D+ towards a dense, starless core in Ophiuchus with the APEX, GBT, and IRAM 30-m telescopes. The observations were interpreted using a gas-grain chemistry model combined with radiative transfer calculations. The chemistry model distinguishes between the different nuclear spin states of light hydrogen molecules, ammonia, and their deuterated forms. High deuterium fractionation ratios with NH2D/NH3=0.4, NHD2/NH2D=0.2, and ND3/NHD2=0.06 were found in the core. The observed ortho/para ratios of NH2D and NHD2 are close to the corresponding nuclear spin statistical weights. The chemistry model can approximately reproduce the observed abundances, but predicts uniformly too low ortho/para-NH2D, and too large ortho/para-NHD2 ratios. The longevity of N2H+ and NH3 in dense gas, which is prerequisite to their strong deuteration, can be attributed to the chemical inertia of N2 on grain surfaces. The discrepancies between the chemistry model and the observations are likely to be caused by the fact that the model assumes complete scrambling in principal gas-phase deuteration reactions of ammonia, which means that all the nuclei are mixed in reactive collisions. If, instead, these reactions occur through proton hop/hydrogen abstraction processes, statistical spin ratios are to be expected. The present results suggest that while the deuteration of ammonia changes with physical conditions and time, the nuclear spin ratios of ammonia isotopologues do not probe the evolutionary stage of a cloud.

A well-instrumented Green Bank Telescope (GBT) operating at high frequency represents a unique scientific resource for the US community. As a filled-aperture, 100m-diameter telescope, the GBT is ideally suited to fast mapping of extended, low surface brightness emission with excellent instantaneous frequency coverage. This capability makes the GBT a key facility for a range of cutting edge science described in this document, only possible at these frequencies. We note that the ability to perform the necessary observations is unique and highly complementary to the capabilities offered by interferometers, and should be preserved. We argue that rather than divesting from this exceptional resource, it makes sense for the US community to invest moderately to maintain GBT operations and to instrument it in an optimal manner, enabling it to become an extraordinary complement to existing and future radio interferometers. Adequately instrumented, the GBT would be a pillar for 20-115 GHz science in the US and the world.

We present a 3 mm spectral line and continuum survey of L1451 in the Perseus Molecular Cloud. These observations are from the CARMA Large Area Star Formation Survey (CLASSy), which also imaged Barnard 1, NGC 1333, Serpens Main and Serpens South. L1451 is the survey region with the lowest level of star formation activity---it contains no confirmed protostars. HCO+, HCN, and N2H+ (J=1-0) are all detected throughout the region, with HCO+ the most spatially widespread, and molecular emission seen toward 90% of the area above N(H_2) column densities of 1.9x10^21 cm^-2. HCO+ has the broadest velocity dispersion, near 0.3 km/s on average, compared to ~0.15 km/s for the other molecules, thus representing a range from supersonic to subsonic gas motions. Our non-binary dendrogram analysis reveals that the dense gas traced by each molecule has similar hierarchical structure, and that gas surrounding the candidate first hydrostatic core (FHSC), L1451-mm, and other previously detected single-dish continuum clumps have similar hierarchical structure; this suggests that different sub-regions of L1451 are fragmenting on the pathway to forming young stars. We determined the three-dimensional morphology of the largest detectable dense gas structures to be relatively ellipsoidal compared to other CLASSy regions, which appeared more flattened at largest scales. A virial analysis shows the most centrally condensed dust structures are likely unstable against collapse. Additionally, we identify a new spherical, centrally condensed N2H+ feature that could be a new FHSC candidate. The overall results suggest L1451 is a young region starting to form its generation of stars within turbulent, hierarchical structures.

We sort $4683$ molecular clouds between $10^\circ< \ell <65^\circ$ from the Bolocam Galactic Plane Survey based on observational diagnostics of star formation activity: compact $70$ $\mu{\rm m}$ sources, mid-IR color-selected YSOs, ${\rm H_2O}$ and ${\rm CH_3OH}$ masers, and UCHII regions. We also present a combined ${\rm NH_3}$-derived gas kinetic temperature and ${\rm H_2O}$ maser catalog for $1788$ clumps from our own GBT 100m observations and from the literature. We identify a subsample of $2223$ ($47.5\%$) starless clump candidates, the largest and most robust sample identified from a blind survey to date. Distributions of flux density, flux concentration, solid angle, kinetic temperature, column density, radius, and mass show strong ($>1$ dex) progressions when sorted by star formation indicator. The median starless clump candidate is marginally sub-virial ($\alpha \sim 0.7$) with $>75\%$ of clumps with known distance being gravitationally bound ($\alpha < 2$). These samples show a statistically significant increase in the median clump mass of $\Delta M \sim 170-370$ M$_\odot$ from the starless candidates to clumps associated with protostars. This trend could be due to (i) mass growth of the clumps at $\dot{M}\sim200-440$ Msun Myr$^{-1}$ for an average free-fall $0.8$ Myr time-scale, (ii) a systematic factor of two increase in dust opacity from starless to protostellar phases, (iii) and/or a variation in the ratio of starless to protostellar clump lifetime that scales as $\sim M^{-0.4}$. By comparing to the observed number of ${\rm CH_3OH}$ maser containing clumps we estimate the phase-lifetime of massive ($M>10^3$ M$_\odot$) starless clumps to be $0.37 \pm 0.08 \ {\rm Myr} \ (M/10^3 \ {\rm M}_\odot)^{-1}$; the majority ($M<450$ M$_\odot$) have phase-lifetimes longer than their average free-fall time.

We present 13CO and C18O (1-0), (2-1), and (3-2) maps towards the core-forming Perseus B1-E clump using observations from the James Clerk Maxwell Telescope (JCMT), Submillimeter Telescope (SMT) of the Arizona Radio Observatory, and IRAM 30 m telescope. We find that the 13CO and C18O line emission both have very complex velocity structures, indicative of multiple velocity components within the ambient gas. The (1-0) transitions reveal a radial velocity gradient across B1-E of 1 km/s/pc that increases from north-west to south-east, whereas the majority of the Perseus cloud has a radial velocity gradient increasing from south-west to north-east. In contrast, we see no evidence of a velocity gradient associated with the denser Herschel-identified substructures in B1-E. Additionally, the denser substructures have much lower systemic motions than the ambient clump material, which indicates that they are likely decoupled from the large-scale gas. Nevertheless, these substructures themselves have broad line widths (0.4 km/s) similar to that of the C18O gas in the clump, which suggests they inherited their kinematic properties from the larger-scale, moderately dense gas. Finally, we find evidence of C18O depletion only toward one substructure, B1-E2, which is also the only object with narrow (transonic) line widths. We suggest that as prestellar cores form, their chemical and kinematic properties are linked in evolution, such that these objects must first dissipate their turbulence before they deplete in CO.

We use the distance probability density function (DPDF) formalism of Ellsworth-Bowers et al. (2013, 2015) to derive physical properties for the collection of 1,710 Bolocam Galactic Plane Survey (BGPS) version 2 sources with well-constrained distance estimates. To account for Malmquist bias, we estimate that the present sample of BGPS sources is 90% complete above 400 $M_\odot$ and 50% complete above 70 $M_\odot$. The mass distributions for the entire sample and astrophysically motivated subsets are generally fitted well by a lognormal function, with approximately power-law distributions at high mass. Power-law behavior emerges more clearly when the sample population is narrowed in heliocentric distance (power-law index $\alpha = 2.0\pm0.1$ for sources nearer than 6.5 kpc and $\alpha = 1.9\pm0.1$ for objects between 2 kpc and 10 kpc). The high-mass power-law indices are generally $1.85 \leq \alpha \leq 2.05$ for various subsamples of sources, intermediate between that of giant molecular clouds and the stellar initial mass function. The fit to the entire sample yields a high-mass power-law $\hat{\alpha} = 1.94_{-0.10}^{+0.34}$. Physical properties of BGPS sources are consistent with large molecular cloud clumps or small molecular clouds, but the fractal nature of the dense interstellar medium makes difficult the mapping of observational categories to the dominant physical processes driving the observed structure. The face-on map of the Galactic disk's mass surface density based on BGPS dense molecular cloud structures reveals the high-mass star-forming regions W43, W49, and W51 as prominent mass concentrations in the first quadrant. Furthermore, we present a 0.25-kpc resolution map of the dense gas mass fraction across the Galactic disk that peaks around 5%.

We present 107 maps of continuum emission at 350 microns from Galactic molecular clumps. Observed sources were mainly selected from the Bolocam Galactic Plane Survey (BGPS) catalog, with 3 additional maps covering star forming regions in the outer Galaxy. The higher resolution of the SHARC-II images (8.5" beam) compared with the 1.1 mm images from BGPS (33" beam) allowed us to identify a large population of smaller substructures within the clumps. A catalog is presented for the 1386 sources extracted from the 350 micron maps. The color temperature distribution of clumps based on the two wavelengths has a median of 13.3 K and mean of 16.3 +- 0.4 K, assuming an opacity law index of 1.7. For the structures with good determination of color temperatures, the mean ratio of gas temperature, determined from NH3 observations, to dust color temperature is 0.88 and the median ratio is 0.76. About half the clumps have more than two substructures and 22 clumps have more than 10. The fraction of the mass in dense substructures seen at 350 microns compared to the mass of their parental clump is ~0.19, and the surface densities of these substructures are, on average, 2.2 times those seen in the clumps identified at 1.1 mm. For a well-characterized sample, 88 structures (31%) exceed a surface density of 0.2 g cm^(-2), and 18 (6%) exceed 1.0 g cm^(-2), thresholds for massive star formation suggested by theorists.

The chemical evolution in high-mass star-forming regions is still poorly constrained. Studying the evolution of deuterated molecules allows to differentiate between subsequent stages of high-mass star formation regions due to the strong temperature dependence of deuterium isotopic fractionation. We observed a sample of 59 sources including 19 infrared dark clouds, 20 high-mass protostellar objects, 11 hot molecular cores and 9 ultra-compact HII regions in the (3-2) transitions of the four deuterated molecules, DCN, DNC, DCO+ and N2D+ as well as their non-deuterated counterpart. The overall detection fraction of DCN, DNC and DCO+ is high and exceeds 50% for most of the stages. N2D+ was only detected in a few infrared dark clouds and high-mass protostellar objects. It can be related to problems in the bandpass at the frequency of the transition and to low abundances in the more evolved, warmer stages. We find median D/H ratios of ~0.02 for DCN, ~0.005 for DNC, ~0.0025 for DCO+ and ~0.02 for N2D+. While the D/H ratios of DNC, DCO+ and N2D+ decrease with time, DCN/HCN peaks at the hot molecular core stage. We only found weak correlations of the D/H ratios for N2D+ with the luminosity of the central source and the FWHM of the line, and no correlation with the H2 column density. In combination with a previously observed set of 14 other molecules (Paper I) we fitted the calculated column densities with an elaborate 1D physico-chemical model with time-dependent D-chemistry including ortho- and para-H2 states. Good overall fits to the observed data have been obtained the model. It is one of the first times that observations and modeling have been combined to derive chemically based best-fit models for the evolution of high-mass star formation including deuteration.

We present deep NH$_3$ observations of the L1495-B218 filaments in the Taurus molecular cloud covering over a 3 degree angular range using the K-band focal plane array on the 100m Green Bank Telescope. The L1495-B218 filaments form an interconnected, nearby, large complex extending over 8 pc. We observed NH$_3$ (1,1) and (2,2) with a spectral resolution of 0.038 km/s and a spatial resolution of 31$"$. Most of the ammonia peaks coincide with intensity peaks in dust continuum maps at 350 $\mu$m and 500 $\mu$m. We deduced physical properties by fitting a model to the observed spectra. We find gas kinetic temperatures of 8 $-$ 15 K, velocity dispersions of 0.05 $-$ 0.25 km/s, and NH$_3$ column densities of 5$\times$10$^{12}$ $-$ 1$\times$10$^{14}$ cm$^{-2}$. The CSAR algorithm, which is a hybrid of seeded-watershed and binary dendrogram algorithms, identifies a total of 55 NH$_3$ structures including 39 leaves and 16 branches. The masses of the NH$_3$ sources range from 0.05 M$_\odot$ to 9.5 M$_\odot$. The masses of NH$_3$ leaves are mostly smaller than their corresponding virial mass estimated from their internal and gravitational energies, which suggests these leaves are gravitationally unbound structures. 9 out of 39 NH$_3$ leaves are gravitationally bound and 7 out of 9 gravitationally bound NH$_3$ leaves are associated with star formation. We also found that 12 out of 30 gravitationally unbound leaves are pressure-confined. Our data suggest that a dense core may form as a pressure-confined structure, evolve to a gravitationally bound core, and undergo collapse to form a protostar.

The optically thin critical densities and the effective excitation densities to produce a 1 K km/s (or 0.818 Jy km/s $(\frac{\nu_{jk}}{100 \rm{GHz}})^2 \, (\frac{\theta_{beam}}{10^{\prime\prime}})^2$) spectral line are tabulated for 12 commonly observed dense gas molecular tracers. The dependence of the critical density and effective excitation density on physical assumptions (i.e. gas kinetic temperature and molecular column density) is analyzed. Critical densities for commonly observed dense gas transitions in molecular clouds (i.e. HCN $1-0$, HCO$^+$ $1-0$, N$_2$H$^+$ $1-0$) are typically $1 - 2$ orders of magnitude larger than effective excitation densities because the standard definitions of critical density do not account for radiative trapping and 1 K km/s lines are typically produced when radiative rates out of the upper energy level of the transition are faster than collisional depopulation. The use of effective excitation density has a distinct advantage over the use of critical density in characterizing the differences in density traced by species such as NH$_3$, HCO$^+$, N$_2$H$^+$, and HCN as well as their isotpologues; but, the effective excitation density has the disadvantage that it is undefined for transitions when $E_u/k \gg T_k$, for low molecular column densities, and for heavy molecules with complex spectra (i.e. CH$_3$CHO).

Post-starburst (or "E+A") galaxies are characterized by low H$\alpha$ emission and strong Balmer absorption, suggesting a recent starburst, but little current star formation. Although many of these galaxies show evidence of recent mergers, the mechanism for ending the starburst is not yet understood. To study the fate of the molecular gas, we search for CO (1-0) and (2-1) emission with the IRAM 30m and SMT 10m telescopes in 32 nearby ($0.01<z<0.12$) post-starburst galaxies drawn from the Sloan Digital Sky Survey. We detect CO in 17 (53%). Using CO as a tracer for molecular hydrogen, and a Galactic conversion factor, we obtain molecular gas masses of $M(H_2)=10^{8.6}$-$10^{9.8} M_\odot$ and molecular gas mass to stellar mass fractions of $\sim10^{-2}$-$10^{-0.5}$, comparable to those of star-forming galaxies. The large amounts of molecular gas rule out complete gas consumption, expulsion, or starvation as the primary mechanism that ends the starburst in these galaxies. The upper limits on $M(H_2)$ for the 15 undetected galaxies range from $10^{7.7} M_\odot$ to $10^{9.7} M_\odot$, with the median more consistent with early-type galaxies than with star-forming galaxies. Upper limits on the post-starburst star formation rates (SFRs) are lower by $\sim10\times$ than for star-forming galaxies with the same $M(H_2)$. We also compare the molecular gas surface densities ($\Sigma_{\rm H_2}$) to upper limits on the SFR surface densities ($\Sigma_{\rm SFR}$), finding a significant offset, with lower $\Sigma_{\rm SFR}$ for a given $\Sigma_{\rm H_2}$ than is typical for star-forming galaxies. This offset from the Kennicutt-Schmidt relation suggests that post-starbursts have lower star formation efficiency, a low CO-to-H$_2$ conversion factor characteristic of ULIRGs, and/or a bottom-heavy initial mass function, although uncertainties in the rate and distribution of current star formation remain.

We present an expanded distance catalog for 1,710 molecular cloud structures identified in the Bolocam Galactic Plane Survey (BGPS) version 2, representing a nearly threefold increase over the previous BGPS distance catalog. We additionally present a new method for incorporating extant data sets into our Bayesian distance probability density function (DPDF) methodology. To augment the dense-gas tracers (e.g., HCO+(3-2), NH3(1,1)) used to derive line-of-sight velocities for kinematic distances, we utilize the Galactic Ring Survey 13CO(1-0) data to morphologically extract velocities for BGPS sources. The outline of a BGPS source is used to select a region of the GRS 13CO data, along with a reference region to subtract enveloping diffuse emission, to produce a line profile of 13CO matched to the BGPS source. For objects with a HCO+(3-2) velocity, \approx 95% of the new 13CO(1-0) velocities agree with that of the dense gas. A new prior DPDF for kinematic distance ambiguity (KDA) resolution, based on a validated formalism for associating molecular cloud structures with known objects from the literature, is presented. We demonstrate this prior using catalogs of masers with trigonometric parallaxes and HII regions with robust KDA resolutions. The distance catalog presented here contains well-constrained distance estimates for 20% of BGPS V2 sources, with typical distance uncertainties \lesssim 0.5 kpc. Approximately 75% of the well-constrained sources lie within 6 kpc of the Sun, concentrated in the Scutum-Centarus arm. Galactocentric positions of objects additionally trace out portions of the Sagittarius, Perseus, and Outer arms in the first and second Galactic quadrants, and we also find evidence for significant regions of interarm dense gas.

We present observations of N2H+(1-0), HCO+(1-0), and HCN(1-0) toward the Serpens Main molecular cloud from the CARMA Large Area Star Formation Survey (CLASSy). We mapped 150 square arcminutes of Serpens Main with an angular resolution of 7 arcsecs. The gas emission is concentrated in two subclusters (the NW and SE subclusters). The SE subcluster has more prominent filamentary structures and more complicated kinematics compared to the NW subcluster. The majority of gas in the two subclusters has subsonic to sonic velocity dispersions. We applied a dendrogram technique with N2H+(1-0) to study the gas structures; the SE subcluster has a higher degree of hierarchy than the NW subcluster. Combining the dendrogram and line fitting analyses reveals two distinct relations: a flat relation between nonthermal velocity dispersion and size, and a positive correlation between variation in velocity centroids and size. The two relations imply a characteristic depth of 0.15 pc for the cloud. Furthermore, we have identified six filaments in the SE subcluster. These filaments have lengths of 0.2 pc and widths of 0.03 pc, which is smaller than a characteristic width of 0.1 pc suggested by Herschel observations. The filaments can be classified into two types based on their properties. The first type, located in the northeast of the SE subcluster, has larger velocity gradients, smaller masses, and nearly critical mass-per-unit-length ratios. The other type, located in the southwest of the SE subcluster, has the opposite properties. Several YSOs are formed along two filaments which have supercritical mass per unit length ratios, while filaments with nearly critical mass-per-unit-length ratios are not associated with YSOs, suggesting that stars are formed on gravitationally unstable filaments.

We present observations of near-infrared 2.12 micro-meter molecular hydrogen outflows emerging from 1.1 mm dust continuum clumps in the North America and Pelican Nebula (NAP) complex selected from the Bolocam Galactic Plane Survey (BGPS). Hundreds of individual shocks powered by over 50 outflows from young stars are identified, indicating that the dusty molecular clumps surrounding the NGC 7000 / IC 5070 / W80 HII region are among the most active sites of on-going star formation in the Solar vicinity. A spectacular X-shaped outflow, MHO 3400, emerges from a young star system embedded in a dense clump more than a parsec from the ionization front associated with the Pelican Nebula (IC 5070). Suspected to be a binary, the source drives a pair of outflows with orientations differing by 80 degrees. Each flow exhibits S-shaped symmetry and multiple shocks indicating a pulsed and precessing jet. The `Gulf of Mexico' located south of the North America Nebula (NGC 7000), contains a dense cluster of molecular hydrogen objects (MHOs), Herbig-Haro (HH) objects, and over 300 YSOs, indicating a recent burst of star formation. The largest outflow detected thus far in the North America and Pelican Nebula complex, the 1.6 parsec long MHO 3417 flow, emerges from a 500 Solar mass BGPS clump and may be powered by a forming massive star. Several prominent outflows such as MHO 3427 appear to be powered by highly embedded YSOs only visible at a wavelength > 70 micro-meters. An `activity index' formed by dividing the number of shocks by the mass of the cloud containing their source stars is used to estimate the relative evolutionary states of Bolocam clumps. Outflows can be used as indicators of the evolutionary state of clumps detected in mm and sub-mm dust continuum surveys.

We present details of the CARMA Large Area Star Formation Survey (CLASSy), while focusing on observations of Barnard 1. CLASSy is a CARMA Key Project that spectrally imaged N2H+, HCO+, and HCN (J=1-0 transitions) across over 800 square arcminutes of the Perseus and Serpens Molecular Clouds. The observations have angular resolution near 7" and spectral resolution near 0.16 km/s. We imaged ~150 square arcminutes of Barnard 1, focusing on the main core, and the B1 Ridge and clumps to its southwest. N2H+ shows the strongest emission, with morphology similar to cool dust in the region, while HCO+ and HCN trace several molecular outflows from a collection of protostars in the main core. We identify a range of kinematic complexity, with N2H+ velocity dispersions ranging from ~0.05-0.50 km/s across the field. Simultaneous continuum mapping at 3 mm reveals six compact object detections, three of which are new detections. A new non-binary dendrogram algorithm is used to analyze dense gas structures in the N2H+ position-position-velocity (PPV) cube. The projected sizes of dendrogram-identified structures range from about 0.01-0.34 pc. Size-linewidth relations using those structures show that non-thermal line-of-sight velocity dispersion varies weakly with projected size, while rms variation in the centroid velocity rises steeply with projected size. Comparing these relations, we propose that all dense gas structures in Barnard 1 have comparable depths into the sky, around 0.1-0.2 pc; this suggests that over-dense, parsec-scale regions within molecular clouds are better described as flattened structures rather than spherical collections of gas. Science-ready PPV cubes for Barnard 1 molecular emission are available for download.

We present the N2H+(J=1-0) map of the Serpens South molecular cloud obtained as part of the CARMA Large Area Star Formation Survey (CLASSy). The observations cover 250 square arcminutes and fully sample structures from 3000 AU to 3 pc with a velocity resolution of 0.16 km/s, and they can be used to constrain the origin and evolution of molecular cloud filaments. The spatial distribution of the N2H+ emission is characterized by long filaments that resemble those observed in the dust continuum emission by Herschel. However, the gas filaments are typically narrower such that, in some cases, two or three quasi-parallel N2H+ filaments comprise a single observed dust continuum filament. The difference between the dust and gas filament widths casts doubt on Herschel ability to resolve the Serpens South filaments. Some molecular filaments show velocity gradients along their major axis, and two are characterized by a steep velocity gradient in the direction perpendicular to the filament axis. The observed velocity gradient along one of these filaments was previously postulated as evidence for mass infall toward the central cluster, but these kind of gradients can be interpreted as projection of large-scale turbulence.

Understanding the chemical evolution of young (high-mass) star-forming regions is a central topic in star formation research. Chemistry is employed as a unique tool 1) to investigate the underlying physical processes and 2) to characterize the evolution of the chemical composition. We observed a sample of 59 high-mass star-forming regions at different evolutionary stages varying from the early starless phase of infrared dark clouds to high-mass protostellar objects to hot molecular cores and, finally, ultra-compact HII regions at 1mm and 3mm with the IRAM 30m telescope. We determined their large-scale chemical abundances and found that the chemical composition evolves along with the evolutionary stages. On average, the molecular abundances increase with time. We modeled the chemical evolution, using a 1D physical model where density and temperature vary from stage to stage coupled with an advanced gas-grain chemical model and derived the best-fit chi^2 values of all relevant parameters. A satisfying overall agreement between observed and modeled column densities for most of the molecules was obtained. With the best-fit model we also derived a chemical age for each stage, which gives the timescales for the transformation between two consecutive stages. The best-fit chemical ages are ~10,000 years for the IRDC stage, ~60,000 years for the HMPO stage, ~40,000 years for the HMC stage, and ~10,000 years for the UCHII stage. The total chemical timescale for the entire evolutionary sequence of the high-mass star formation process is on the order of 10^5 years, which is consistent with theoretical estimates. Furthermore, based on the approach of a multiple-line survey of unresolved data, we were able to constrain an intuitive and reasonable physical and chemical model. The results of this study can be used as chemical templates for the different evolutionary stages in high-mass star formation.

The Bolocam Galactic Plane Survey (BGPS) is a 1.1 mm continuum survey of dense clumps of dust throughout the Galaxy covering 170 square degrees. We present spectroscopic observations using the Heinrich Hertz Submillimeter Telescope of the dense gas tracers, HCO+ and N2H+ 3-2, for all 6194 sources in the Bolocam Galactic Plane Survey v1.0.1 catalog between 7.5 <= l <= 194 degrees. This is the largest targeted spectroscopic survey of dense molecular gas in the Milky Way to date. We find unique velocities for 3126 (50.5%) of the BGPS v1.0.1 sources observed. Strong N2H+ 3-2 emission (T_{mb} > 0.5 K) without HCO+ 3-2 emission does not occur in this catalog. We characterize the properties of the dense molecular gas emission toward the entire sample. HCO+ is very sub-thermally populated and the 3-2 transitions are optically thick toward most BGPS clumps. The median observed line width is 3.3 km/s consistent with supersonic turbulence within BGPS clumps. We find strong correlations between dense molecular gas integrated intensities and 1.1 mm peak flux and the gas kinetic temperature derived from previously published NH3 observations. These intensity correlations are driven by the sensitivity of the 3-2 transitions to excitation conditions rather than by variations in molecular column density or abundance. We identify a subset of 113 sources with stronger N2H+ than HCO+ integrated intensity, but we find no correlations between the N2H+ / HCO+ ratio and 1.1 mm continuum flux density, gas kinetic temperature, or line width. Self-absorbed profiles are rare (1.3%).

We present a re-reduction and expansion of the Bolocam Galactic Plane Survey, first presented by Aguirre et al. (2011) and Rosolowsky et al. (2010). The BGPS is a 1.1 mm survey of dust emission in the Northern galactic plane, covering longitudes -10 < \ell < 90 and latitudes |b| < 0.5 with a typical 1-\sigma RMS sensitivity of 30-100 mJy in a 33" beam. Version 2 of the survey includes an additional 20 square degrees of coverage in the 3rd and 4th quadrants and 2 square degrees in the 1st quadrant. The new data release has improved angular recovery, with complete recovery out to 80" and partial recovery to 300", and reduced negative bowls around bright sources resulting from the atmospheric subtraction process. We resolve the factor of 1.5 flux calibration offset between the v1.0 data release and other data sets and determine that there is no offset between v2.0 and other data sets. The v2.0 pointing accuracy is tested against other surveys and demonstrated to be accurate and an improvement over v1.0. We present simulations and tests of the pipeline and its properties, including measurements of the pipeline's angular transfer function. The Bolocat cataloging tool was used to extract a new catalog, which includes 8594 sources, with 591 in the expanded regions. We have demonstrated that the Bolocat 40" and 80" apertures are accurate even in the presence of strong extended background emission. The number of sources is lower than in v1.0, but the amount of flux and area included in identified sources is larger.

We present a new distance estimation method for dust-continuum-identified molecular cloud clumps. Recent (sub-)millimeter Galactic plane surveys have cataloged tens of thousands of these objects, but detailed study of their physical properties requires robust distance determinations. We derive Bayesian distance probability density functions (DPDFs) for 770 objects from the Bolocam Galactic Plane Survey in the longitude range 7.5 < l < 65. The DPDF formalism is based on kinematic distances, and uses external data sets to place prior distance probabilities to resolve the kinematic distance ambiguity (KDA) for objects in the inner Galaxy. We present priors related to the mid-infrared absorption of dust in dense molecular regions and the Galactic distribution of molecular gas. By assuming a numerical model of Galactic mid-infrared emission and simple radiative transfer, we match the morphology of millimeter thermal dust emission with mid-infrared absorption to compute a prior DPDF for distance discrimination. Selecting objects first from (sub-)millimeter source catalogs avoids a bias towards the darkest infrared dark clouds (IRDCs) and extends the range of heliocentric distance and contrast range probed by mid-infrared extinction. We derive well-constrained KDA resolutions for 618 molecular cloud clumps, with approximately 15% placed at or beyond the tangent distance. Objects sufficiently dark to be cataloged as IRDCs are generally placed at the near kinematic distance. Distance comparisons with Galactic Ring Survey KDA resolutions yield a 92% agreement, and a Galactic face-on view reveals sections of the Sagittarius and Scutum-Centaurus Arms. This KDA-resolution method for large catalogs of sources through the combination of (sub-)millimeter and mid-infrared observations of molecular cloud clumps is generally applicable to other dust-continuum Galactic plane surveys.

In order to understand the collapse dynamics of observed low-mass starless cores, we revise the conventional stability condition of hydrostatic Bonnor-Ebert spheres to take internal motions into account. Because observed starless cores resemble Bonnor-Ebert density structures, the stability and dynamics of the starless cores are frequently analyzed by comparing to the conventional stability condition of a hydrostatic Bonnor-Ebert sphere. However, starless cores are not hydrostatic but have observed internal motions. In this study, we take gaseous spheres with a homologous internal velocity field and derive stability conditions of the spheres utilizing a virial analysis. We propose two limiting models of spontaneous gravitational collapse: the collapse of critical Bonnor-Ebert spheres and uniform density spheres. The collapse of these two limiting models are intended to provide the lower and the upper limits, respectively, of the infall speeds for a given density structure. The results of our study suggest that the stability condition sensitively depends on internal motions. A homologous inward motion with a transonic speed can reduce the critical size compared to the static Bonnor-Ebert sphere by more than a factor of two. As an application of the two limiting models of spontaneous gravitational collapse, we compare the density structures and infall speeds of the observed starless cores L63, L1544, L1689B, and L694-2 to the two limiting models. L1689B and L694-2 seem to have been perturbed to result in faster infall motions than for spontaneous gravitational collapse.

We have performed a pointed survey of N2D+ 2-1 and N2D+ 3-2 emission toward 64 N2H+-bright starless and protostellar cores in the Perseus molecular cloud using the Arizona Radio Observatory Submillimeter Telescope and Kitt Peak 12 m telescope. We find a mean deuterium fractionation in N2H+, R_D = N(N2D+)/N(N2H+), of 0.08, with a maximum R_D = 0.2. In detected sources, we find no significant difference in the deuterium fractionation between starless and protostellar cores, nor between cores in clustered or isolated environments. We compare the deuterium fraction in N2H+ with parameters linked to advanced core evolution. We only find significant correlations between the deuterium fraction and increased H_2 column density, as well as with increased central core density, for all cores. Towards protostellar sources, we additionally find a significant anti-correlation between R_D and bolometric temperature. We show that the Perseus cores are characterized by low CO depletion values relative to previous studies of star forming cores, similar to recent results in the Ophiuchus molecular cloud. We suggest that the low average CO depletion is the dominant mechanism that constrains the average deuterium fractionation in the Perseus cores to small values. While current equilibrium and dynamic chemical models are able to reproduce the range of deuterium fractionation values we find in Perseus, reproducing the scatter across the cores requires variation in parameters such as the ionization fraction or the ortho- to para-H_2 ratio across the cloud, or a range in core evolution timescales.

The temperature and density structure of molecular cloud cores are the most important physical quantities that determine the course of the protostellar collapse and the properties of the stars they form. Nevertheless, density profiles often rely either on the simplifying assumption of isothermality or on observationally poorly constrained model temperature profiles. With the aim of better constraining the initial physical conditions in molecular cloud cores at the onset of protostellar collapse, we initiated the Guaranteed Time Key Project (GTKP) "The Earliest Phases of Star Formation" (EPoS) with the Herschel satellite. This paper gives an overview of the low-mass sources in the EPoS project, including all observations, the analysis method, and the initial results of the survey. We study the thermal dust emission of 12 previously well-characterized, isolated, nearby globules using FIR and submm continuum maps at up to eight wavelengths between 100 micron and 1.2 mm. Our sample contains both globules with starless cores and embedded protostars at different early evolutionary stages. The dust emission maps are used to extract spatially resolved SEDs, which are then fit independently with modified blackbody curves to obtain line-of-sight-averaged dust temperature and column density maps. We find that the thermal structure of all globules is dominated by external heating from the interstellar radiation field and moderate shielding by thin extended halos. All globules have warm outer envelopes (14-20 K) and colder dense interiors (8-12 K). The protostars embedded in some of the globules raise the local temperature of the dense cores only within radii out to about 5000 AU, but do not significantly affect the overall thermal balance of the globules.

We present a systematic single-dish search for molecular outflows toward a sample of 9 candidate low-luminosity protostars and 30 candidate Very Low Luminosity Objects (VeLLOs; L_int < 0.1 L_sun). The sources are identified using data from the Spitzer Space Telescope catalogued by Dunham et al. toward nearby (D < 400 pc) star forming regions. Each object was observed in 12CO and 13CO J = 2-1 simultaneously using the sideband separating ALMA Band-6 prototype receiver on the Heinrich Hertz Telescope at 30 arcsecond resolution. Using 5-point grid maps we identify five new potential outflow candidates and make on-the-fly maps of the regions surrounding sources in the dense cores B59, L1148, L1228, and L1165. Of these new outflow candidates, only the map of B59 shows a candidate blue outflow lobe associated with a source in our survey. We also present larger and more sensitive maps of the previously detected L673-7 and the L1251-A IRS4 outflows and analyze their properties in comparison to other outflows from VeLLOs. The accretion luminosities derived from the outflow properties of the VeLLOs with detected CO outflows are higher than the observed internal luminosity of the protostars, indicating that these sources likely had higher accretion rates in the past. The known L1251-A IRS3 outflow is detected but not remapped. We do not detect clear, unconfused signatures of red and blue molecular wings toward the other 31 sources in the survey indicating that large-scale, distinct outflows are rare toward this sample of candidate protostars. Several potential outflows are confused with kinematic structure in the surrounding core and cloud. Interferometric imaging is needed to disentangle large-scale molecular cloud kinematics from these potentially weak protostellar outflows.

The millimeter-wave rotational transition line ($J_K = 1_0 - 0_0$) of deuterated methane CH$_3$D has tentatively been detected toward the low-mass Class 0 protostar IRAS 04368+2557 in L1527 with the Heinrich Hertz Submillimeter Telescope. This is the first detection of CH$_3$D in interstellar clouds, if confirmed. The column density and fractional abundance of CH$_3$D are determined to be $(9.1\pm3.4) \times 10^{15}$ cm$^{-2}$ and $(3.0\pm1.1) \times 10^{-7}$, respectively, where we assume the rotational temperature of 25 K. The column density and fractional abundance of the gaseous CH$_4$ are estimated to be $(1.3-4.6) \times 10^{17}$ cm$^{-2}$ and $(4.3-15.2) \times 10^{-6}$, respectively, by adopting the molecular D/H ratios of 2--7% reported for various molecules in L1527. The fractional abundance of CH$_4$ is higher than or comparable to that found in high-mass star-forming cores by infrared observations. It is sufficiently high for triggering the efficient production of various carbon-chain molecules in a lukewarm region near the protostar, which supports the picture of the warm carbon-chain chemistry.

We present an analysis of the relation between the star formation rate (SFR) and mass of dense gas in Galactic clumps and nearby galaxies. Using the bolometric luminosity as a measure of SFR and the molecular line luminosity of HCO+ (3-2) as a measure of dense gas mass, we find that the relation between SFR and M_{dense} is approximately linear. This is similar to published results derived using HCN (1-0) as a dense gas tracer. HCO+ (3-2) and HCN (1-0) have similar conditions for excitation. Our work includes 16 Galactic clumps that are in both the Bolocam Galactic Plane Survey and the Red MSX Survey, 27 water maser sources from the literature, and the aforementioned HCN (1-0) data. Our results agree qualitatively with predictions of recent theoretical models which state that the nature of the relation should depend on how the critical density of the tracer compares with the mean density of the gas.

We analyze the HCO+ 3-2 and H13CO+ 3-2 line profiles of 27 high-mass star-forming regions to identify asymmetries that are suggestive of mass inflow. Three quantitative measures of line asymmetry are used to indicate whether a line profile is blue, red or neither - the ratio of the temperature of the blue and red peaks, the line skew and the dimensionless parameter delta_v. We find nine HCO+ 3-2 line profiles with a significant blue asymmetry and four with significant red asymmetric profiles. Comparing our HCO+ 3-2 results to HCN 3-2 observations from Wu et al. (2003, 2010), we find that eight of the blue and three of red have profiles with the same asymmetry in HCN. The eight sources with blue asymmetries in both tracers are considered strong candidates for inflow. Quantitative measures of the asymmetry (e.g. delta_v) tend to be larger for HCN. This, combined with possible HCO+ abundance enhancements in outflows, suggests that HCN may be a better tracer of inflow. Understanding the behavior of common molecular tracers like HCO+ in clumps of different masses is important for properly analyzing the line profiles seen in a sample of sources representing a broad range of clump masses. Such studies will soon be possible with the large number of sources with possible self-absorption seen in spectroscopic follow-up observations of clumps identified in the Bolocam Galactic Plane Survey.

We present the results of observations of 1882 sources in the Bolocam Galactic Plane Survey (BGPS) at 1.1 mm with the 10m Heinrich Hertz Telescope simultaneously in HCO+ J=3-2 and N2H+ J=3-2. We detect 77% of these sources in HCO^+ and 51% in N2H+ at greater than 3$\sigma$. We find a strong correlation between the integrated intensity of both dense gas tracers and the 1.1 mm dust emission of BGPS sources. We determine kinematic distances for 529 sources (440 in the first quadrant breaking the distance ambiguity and 89 in the second quadrant) We derive the size, mass, and average density for this subset of clumps. The median size of BGPS clumps is 0.75 pc with a median mass of 330 M$_{\odot}$ (assuming T_{Dust}=20 K). The median HCO+ linewidth is 2.9 km s$^{-1}$ indicating that BGPS clumps are dominated by supersonic turbulence or unresolved kinematic motions. We find no evidence for a size-linewidth relationship for BGPS clumps. We analyze the effects of the assumed dust temperature on the derived clump properties with a Monte Carlo simulation and we find that changing the temperature distribution will change the median source properties (mass, volume-averaged number density, surface density) by factors of a few. The observed differential mass distribution has a power-law slope that is intermediate between that observed for diffuse CO clouds and the stellar IMF. BGPS clumps represent a wide range of objects (from dense cores to more diffuse clumps) and are typically characterized by larger sizes and lower densities than previously published surveys of high-mass star-forming regions. This collection of objects is a less-biased sample of star-forming regions in the Milky Way that likely span a wide range of evolutionary states.

We present observations of HCO+ and H^13CO+, N2H+, HCS+, HCN and HN^13C, SO and ^34SO, CCH, SO_2, and CH_3OH-E towards a sample of 27 high-mass clumps coincident with water maser emission. All transitions are observed with or convolved to nearly identical resolution (30"), allowing for inter-comparison of the clump properties derived from the mapped transitions. We find N2H+ emission is spatially differentiated compared to the dust and the other molecules towards a few very luminous cores (10 of 27) and the N2H+ integrated intensity does not correlate well with dust continuum flux. We calculate the effective excitation density, n_eff, the density required to excite a 1 K line in T_kin=20 K gas for each molecular tracer. The intensity of molecular tracers with larger effective excitation densities (n_eff > 10^5 cm^-3) appear to correlate more strongly with the submillimeter dust continuum intensity. The median sizes of the clumps are anti-correlated with the n_eff of the tracers (which span more than three orders of magnitude). Virial mass is not correlated with n_eff, especially where the lines are optically thick as the linewidths may be broadened significantly by non-virial motions. The median mass surface density and median volume density of the clumps is correlated with n_eff indicating the importance of understanding the excitation conditions of the molecular tracer when deriving the average properties of an ensemble of cores.

We present Very Large Array observations of ammonia (NH3) (1,1), (2,2), and CCS (2_1-1_0) emission toward the Infrared Dark Cloud (IRDC) G19.30+0.07 at ~22GHz. The NH3 emission closely follows the 8 micron extinction. The NH3 (1,1) and (2,2) lines provide diagnostics of the temperature and density structure within the IRDC, with typical rotation temperatures of ~10 to 20K and NH3 column densities of ~10^15 cm^-2. The estimated total mass of G19.30+0.07 is ~1130 Msun. The cloud comprises four compact NH3 clumps of mass ~30 to 160 Msun. Two coincide with 24 micron emission, indicating heating by protostars, and show evidence of outflow in the NH3 emission. We report a water maser associated with a third clump; the fourth clump is apparently starless. A non-detection of 8.4GHz emission suggests that the IRDC contains no bright HII regions, and places a limit on the spectral type of an embedded ZAMS star to early-B or later. From the NH3 emission we find G19.30+0.07 is composed of three distinct velocity components, or "subclouds." One velocity component contains the two 24 micron sources and the starless clump, another contains the clump with the water maser, while the third velocity component is diffuse, with no significant high-density peaks. The spatial distribution of NH3 and CCS emission from G19.30+0.07 is highly anti-correlated, with the NH3 predominantly in the high-density clumps, and the CCS tracing lower-density envelopes around those clumps. This spatial distribution is consistent with theories of evolution for chemically young low-mass cores, in which CCS has not yet been processed to other species and/or depleted in high-density regions.

We present here findings for C18O depletion in eight starless cores in Taurus: TMC-2, L1498, L1512, L1489, L1517B, L1521E, L1495A-S, and L1544. We compare observations of the C18O J=2-1 transition taken with the ALMA prototype receiver on the Heinrich Hertz Submillimeter Telescope to results of radiative transfer modeling using RATRAN. We use temperature and density profiles calculated from dust continuum radiative transfer models to model the C18O emission. We present modeling of three cores, TMC-2, L1489, and L1495A-S, which have not been modeled before and compare our results for the five cores with published models. We find that all of the cores but one, L1521E, are substantially depleted. We also find that varying the temperature profiles of these model cores has a discernable effect, and varying the central density has an even larger effect. We find no trends with depletion radius or depletion fraction with the density or temperature of these cores, suggesting that the physical structure alone is insufficient to fully constrain evolutionary state. We are able to place tighter constraints on the radius at which C18O is depleted than the absolute fraction of depletion. As the timeline of chemical depletion depends sensitively on the fraction of depletion, this difficulty in constraining depletion fraction makes comparison with other timescales, such as the free-fall timescale, very difficult. Keywords: Stars:Formation -- ISM:abundances-- ISM:clouds -- Individual Objects: TMC-2, L1498, L1512, L1489, L1517B, L1521E, L1495A-S, L1544

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.

Infrared extinction maps and submillimeter dust continuum maps are powerful probes of the density structure in the envelope of star-forming cores. We make a direct comparison between infrared and submillimeter dust continuum observations of the low-mass Class 0 core, B335, to constrain the ratio of submillimeter to infrared opacity (\kaprat) and the submillimeter opacity power-law index ($\kappa \propto \lambda^{-\beta}$). Using the average value of theoretical dust opacity models at 2.2 \micron, we constrain the dust opacity at 850 and 450 \micron . Using new dust continuum models based upon the broken power-law density structure derived from interferometric observations of B335 and the infall model derived from molecular line observations of B335, we find that the opacity ratios are $\frac{\kappa_{850}}{\kappa_{2.2}} = (3.21 - 4.80)^{+0.44}_{-0.30} \times 10^{-4}$ and $\frac{\kappa_{450}}{\kappa_{2.2}} = (12.8 - 24.8)^{+2.4}_{-1.3} \times 10^{-4}$ with a submillimeter opacity power-law index of $\beta_{smm} = (2.18 - 2.58)^{+0.30}_{-0.30}$. The range of quoted values are determined from the uncertainty in the physical model for B335. For an average 2.2 \micron\ opacity of $3800 \pm 700$ cm$^2$g$^{-1}$, we find a dust opacity at 850 and 450 \micron\ of $\kappa_{850} = (1.18 - 1.77)^{+0.36}_{-0.24}$ and $\kappa_{450} = (4.72 - 9.13)^{+1.9}_{-0.98}$ cm$^2$g$^{-1}$ of dust. These opacities are from $(0.65 - 0.97) \kappa^{\rm{OH}5}_{850}$ of the widely used theoretical opacities of Ossenkopf and Henning for coagulated ice grains with thin mantles at 850\micron.

We present observations of six Class 0 protostars at 3.3 mm (90 GHz) using the 64-pixel MUSTANG bolometer camera on the 100-m Green Bank Telescope. The 3.3 mm photometry is analyzed along with shorter wavelength observations to derive spectral indices (S_nu ~ nu^alpha) of the measured emission. We utilize previously published dust continuum radiative transfer models to estimate the characteristic dust temperature within the central beam of our observations. We present constraints on the millimeter dust opacity index, beta, between 0.862 mm, 1.25 mm, and 3.3 mm. Beta_mm typically ranges from 1.0 to 2.4 for Class 0 sources. The relative contributions from disk emission and envelope emission are estimated at 3.3 mm. L483 is found to have negligible disk emission at 3.3 mm while L1527 is dominated by disk emission within the central beam. The beta_mm^disk <= 0.8 - 1.4 for L1527 indicates that grain growth is likely occurring in the disk. The photometry presented in this paper may be combined with future interferometric observations of Class 0 envelopes and disks.

The Bolocam Galactic Plane Survey (BGPS) data for a six square degree region of the Galactic plane containing the Galactic center is analyzed and compared to infrared and radio continuum data. The BGPS 1.1 mm emission consists of clumps interconnected by a network of fainter filaments surrounding cavities, a few of which are filled with diffuse near-IR emission indicating the presence of warm dust or with radio continuum characteristic of HII regions or supernova remnants. New 350 {\mu}m images of the environments of the two brightest regions, Sgr A and B, are presented. Sgr B2 is the brightest mm-emitting clump in the Central Molecular Zone and may be forming the closest analog to a super star cluster in the Galaxy. The Central Molecular Zone (CMZ) contains the highest concentration of mm and sub-mm emitting dense clumps in the Galaxy. Most 1.1 mm features at positive longitudes are seen in silhouette against the 3.6 to 24 {\mu}m background observed by the Spitzer Space Telescope. However, only a few clumps at negative longitudes are seen in absorption, confirming the hypothesis that positive longitude clumps in the CMZ tend to be on the near-side of the Galactic center, consistent with the suspected orientation of the central bar in our Galaxy. Some 1.1 mm cloud surfaces are seen in emission at 8 {\mu}m, presumably due to polycyclic aromatic hydrocarbons (PAHs). A ~0.2\degree (~30 pc) diameter cavity and infrared bubble between l \approx 0.0\degree and 0.2\degree surrounds the Arches and Quintuplet clusters and Sgr A. The bubble contains several clumpy dust filaments that point toward Sgr A\ast; its potential role in their formation is explored. [abstract truncated]

We examine the physical properties and evolutionary stages of a sample of 17 clumps within 8 Infrared Dark Clouds (IRDCs) by combining existing infrared, millimeter, and radio data with new Bolocam Galactic Plane Survey (BGPS) 1.1 mm data, VLA radio continuum data, and HHT dense gas (HCO+ and N2H+) spectroscopic data. We combine literature studies of star formation tracers and dust temperatures within IRDCs with our search for ultra-compact (UC) HII regions to discuss a possible evolutionary sequence for IRDC clumps. In addition, we perform an analysis of mass tracers in IRDCs and find that 8 micron extinction masses and 1.1 mm Bolocam Galactic Plane Survey (BGPS) masses are complementary mass tracers in IRDCs except for the most active clumps (notably those containing UCHII regions), for which both mass tracers suffer biases. We find that the measured virial masses in IRDC clumps are uniformly higher than the measured dust continuum masses on the scale of ~1 pc. We use 13CO, HCO+, and N2H+ to study the molecular gas properties of IRDCs and do not see any evidence of chemical differentiation between hot and cold clumps on the scale of ~1 pc. However, both HCO+ and N2H+ are brighter in active clumps, due to an increase in temperature and/or density. We report the identification of four UCHII regions embedded within IRDC clumps and find that UCHII regions are associated with bright (>1 Jy) 24 micron point sources, and that the brightest UCHII regions are associated with "diffuse red clumps" (an extended enhancement at 8 micron). The broad stages of the discussed evolutionary sequence (from a quiescent clump to an embedded HII region) are supported by literature dust temperature estimates; however, no sequential nature can be inferred between the individual star formation tracers.

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.

We present the 1.1 millimeter Bolocam Galactic Plane Survey (BGPS) observations of the Gemini OB1 molecular cloud complex, and targeted ammonia observations of the BGPS sources. When paired with molecular spectroscopy of a dense gas tracer, millimeter observations yield physical properties such as masses, radii, mean densities, kinetic temperatures and line widths. We detect 34 distinct BGPS sources above 5-sigma=0.37 Jy/beam with corresponding 5-sigma detections in the ammonia (1,1) transition. Eight of the objects show water maser emission (20%). We find a mean millimeter source FWHM of 1.12 pc, and a mean kinetic temperature of 20 K for the sample of 34 BGPS sources. The observed ammonia line widths are dominated by non-thermal motions, typically found to be a few times the thermal sound speed expected for the derived kinetic temperature. We calculate the mass for each source from the millimeter flux assuming the sources are isothermal and find a mean isothermal mass within a 120" aperture of 230 +/- 180 solar masses. We find a total mass of 8,400 solar masses for all BGPS sources in the Gemini OB1 molecular cloud, representing 6.5% of the cloud mass. By comparing the millimeter isothermal mass to the virial mass within a radius equal to the mm source size calculated from the ammonia line widths, we find a mean virial parameter (M_vir/M_iso) of 1.0 +/- 0.9 for the sample. We find mean values for the distributions of column densities of 10^22 cm^-2 for H_2, and 3.0x10^14 cm^-2 for ammonia, giving a mean ammonia abundance of 3.0x10^-8 relative to H_2. We find volume-averaged densities on the order of 10^3-10^4 cm^-3. The sizes and densities suggest that in the Gem OB1 region the BGPS is detecting the clumps from which stellar clusters form, rather than smaller, higher density cores where single stars or small multiple systems form.

Aims: To explore the infrared and radio properties of one of the closest Galactic starburst regions. Methods: Images obtained with the Herschel Space Observatory at wavelengths of 70, 160, 250, 350, and 500 microns using the PACS and SPIRE arrays are analyzed and compared with radio continuum VLA data and 8 micron images from the Spitzer Space Telescope. The morphology of the far-infrared emission is combined with radial velocity measurements of millimeter and centimeter wavelength transitions to identify features likely to be associated with the W43 complex. Results: The W43 star-forming complex is resolved into a dense cluster of protostars, infrared dark clouds, and ridges of warm dust heated by massive stars. The 4 brightest compact sources with L > 1.5 x 10^4 Lsun embedded within the Z-shaped ridge of bright dust emission in W43 remain single at 4" (0.1 pc) resolution. These objects, likely to be massive protostars or compact clusters in early stages of evolution are embedded in clumps with masses of 10^3 to 10^4 Msun, but contribute only 2% to the 3.6 x 10^6 Lsun far-IR luminosity of W43 measured in a 16 by 16 pc box. The total mass of gas derived from the far-IR dust emission inside this region is ~10^6 Msun. Cometary dust clouds, compact 6 cm radio sources, and warm dust mark the locations of older populations of massive stars. Energy release has created a cavity blowing-out below the Galactic plane. Compression of molecular gas in the plane by the older HII region near G30.684-0.260 and the bipolar structure of the resulting younger W43 HII region may have triggered the current mini-star burst.

We have mapped over 50 massive, dense clumps with four dense gas tracers: HCN J=1-0 and 3-2; and CS J=2-1 and 7-6 transitions. Spectral lines of optically thin H^{13}CN 3-2 and C^{34}S 5-4 were also obtained towards the map centers. These maps usually demonstrate single well-peaked distributions at our resolution, even with higher J transitions. The size, virial mass, surface density, and mean volume density within a well-defined angular size (FWHM) were calculated from the contour maps for each transition. We found that transitions with higher effective density usually trace the more compact, inner part of the clumps but have larger linewidths, leading to an inverse linewidth-size relation using different tracers. The mean surface densities are 0.29, 0.33, 0.78, 1.09 g cm^{-2} within FWHM contours of CS 2-1, HCN 1-0, HCN 3-2 and CS 7-6, respectively. We find no correlation of L_{IR} with surface density and a possible inverse correlation with mean volume density, contrary to some theoretical expectations. We see no evidence in the data for the relation between L'_{mol} and mean density posited by modelers. The correlation between L'_{mol} and the virial mass is roughly linear for each dense gas tracer. A nearly linear correlation was found between the infrared luminosity and the line luminosity of all dense gas tracers for these massive, dense clumps, with a lower cutoff in luminosity at L_{IR}=10^{4.5} Lsun. The L_{IR}-L'_{HCN1-0} correlation agrees well with the one found in galaxies. These correlations indicate a constant star formation rate per unit mass from the scale of dense clumps to that of distant galaxies when the mass is measured for dense gas. These results support the suggestion that starburst galaxies may be understood as having a large fraction of gas in dense clumps.

We present a detailed analysis of the relation between infrared luminosity and molecular line luminosity, for a variety of molecular transitions, using a sample of 34 nearby galaxies spanning a broad range of infrared luminosities (10^{10} < L_{IR} < 10^{12.5} L_sun). We show that the power-law index of the relation is sensitive to the critical density of the molecular gas tracer used, and that the dominant driver in observed molecular line ratios in galaxies is the gas density. As most nearby ultraluminous infrared galaxies (ULIRGs) exhibit strong signatures of active galactic nuclei (AGN) in their center, we revisit previous claims questioning the reliability of HCN as a probe of the dense gas responsible for star formation in the presence of AGN. We find that the enhanced HCN(1-0)/CO(1-0) luminosity ratio observed in ULIRGs can be successfully reproduced using numerical models with fixed chemical abundances and without AGN-induced chemistry effects. We extend this analysis to a total of ten molecular line ratios by combining the following transitions: CO(1-0), HCO+(1-0), HCO+(3-2), HCN(1-0), and HCN(3-2). Our results suggest that AGNs reside in systems with higher dense gas fraction, and that chemistry or other effects associated with their hard radiation field may not dominate (NGC 1068 is one exception). Galaxy merger could be the underlying cause of increased dense gas fraction and the evolutionary stage of such mergers may be another determinant of the HCN/CO luminosity ratio.

We present Spitzer observations of a sample of 12 starless cores selected to have prominent 24 micron shadows. The Spitzer images show 8 and 24 micron shadows and in some cases 70 micron shadows; these spatially resolved absorption features trace the densest regions of the cores. We have carried out a 12CO (2-1) and 13CO (2-1) mapping survey of these cores with the Heinrich Hertz Telescope (HHT). We use the shadow features to derive optical depth maps. We derive molecular masses for the cores and the surrounding environment; we find that the 24 micron shadow masses are always greater than or equal to the molecular masses derived in the same region, a discrepancy likely caused by CO freeze--out onto dust grains. We combine this sample with two additional cores that we studied previously to bring the total sample to 14 cores. Using a simple Jeans mass criterion we find that ~ 2/3 of the cores selected to have prominent 24 micron shadows are collapsing or near collapse, a result that is supported by millimeter line observations. Of this subset at least half have indications of 70 micron shadows. All cores observed to produce absorption features at 70 micron are close to collapse. We conclude that 24 micron shadows, and even more so the 70 micron ones, are useful markers of cloud cores that are approaching collapse.

In an effort to resolve the discrepancy between two measurements of the fundamental constant mu, the proton to electron mass ratio, at early times in the universe we reanalyze the same data used in the earlier studies. Our analysis of the molecular hydrogen absorption lines in archival VLT/UVES spectra of the damped Lyman alpha systems in the QSOs Q0347-383 and Q0405-443 yields a combined measurement of a (Delta mu)/mu value of (-7 +/- 8) x 10^{-6}, consistent with no change in the value of mu over a time span of 11.5 gigayears. Here we define (Delta mu) as (mu_z - mu_0) where mu_z is the value of mu at a redshift of z and mu_0 is the present day value. Our null result is consistent with the recent measurements of King et al. 2009, (Delta mu)/u = (2.6 +/- 3.0) x 10^{-6}, and inconsistent with the positive detection of a change in mu by Reinhold et al. 2006. Both of the previous studies and this study are based on the same data but with differing analysis methods. Improvements in the wavelength calibration over the UVES pipeline calibration is a key element in both of the null results. This leads to the conclusion that the fundamental constant mu is unchanged to an accuracy of 10^{-5} over the last 80% of the age of the universe, well into the matter dominated epoch. This limit provides constraints on models of dark energy that invoke rolling scalar fields and also limits the parameter space of Super Symmetric or string theory models of physics. New instruments, both planned and under construction, will provide opportunities to greatly improve the accuracy of these measurements.

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.

We present combined interferometer and single dish telescope data of NH3 (J,K) = (1,1) and (2,2) emission towards the clustered star forming Ophiuchus B, C and F Cores at high spatial resolution (~1200 AU) using the Australia Telescope Compact Array, the Very Large Array, and the Green Bank Telescope. While the large scale features of the NH3 (1,1) integrated intensity appear similar to 850 micron continuum emission maps of the Cores, on 15" (1800 AU) scales we find significant discrepancies between the dense gas tracers in Oph B, but good correspondence in Oph C and F. Using the Clumpfind structure identifying algorithm, we identify 15 NH3 clumps in Oph B, and 3 each in Oph C and F. Only five of the Oph B NH3 clumps are coincident within 30" (3600 AU) of a submillimeter clump. We find v_LSR varies little across any of the Cores, and additionally varies by only ~1.5 km/s between them. The observed NH3 line widths within the Oph B and F Cores are generally large and often mildly supersonic, while Oph C is characterized by narrow line widths which decrease to nearly thermal values. We find several regions of localized narrow line emission (\Delta v < 0.4 km/s), some of which are associated with NH3 clumps. We derive the kinetic temperatures of the gas, and find they are remarkably constant across Oph B and F, with a warmer mean value (T_K = 15 K) than typically found in isolated regions and consistent with previous results in clustered regions. Oph C, however, has a mean T_K = 12 K, decreasing to a minimum T_K = 9.4 K towards the submillimeter continuum peak, similar to previous studies of isolated starless cores. There is no significant difference in temperature towards protostars embedded in the Cores. [Abridged]

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 infrared and millimeter observations of Barnard 335, the prototypical isolated Bok globule with an embedded protostar. Using Spitzer data we measure the source luminosity accurately; we also constrain the density profile of the innermost globule material near the protostar using the observation of an 8.0 um shadow. HHT observations of 12CO 2 --> 1 confirm the detection of a flattened molecular core with diameter ~10000 AU and the same orientation as the circumstellar disk (~100 to 200 AU in diameter). This structure is probably the same as that generating the 8.0 um shadow and is expected from theoretical simulations of collapsing embedded protostars. We estimate the mass of the protostar to be only ~5% of the mass of the parent globule.

We present observations made with the 10m Heinrich Hertz Sub-Millimeter Telescope of HCN (3-2) emission from a sample of 30 nearby galaxies ranging in infrared luminosity from 10^10 - 10^12.5 L_sun and HCN (3-2) luminosity from 10^6 - 10^9 K km s^-1 pc^2. We examine the correlation between the infrared luminosity and HCN (3-2) luminosity and find that the best fit linear regression has a slope (in log-log space) of 0.74+/-0.12. Including recently published data from Gracia-Carpio et al. tightens the constraints on the best-fit slope to 0.79+/-0.09. This slope below unity suggests that the HCN (3-2) molecular line luminosity is not linearly tracing the amount of dense gas. Our results are consistent with predictions from recent theoretical models that find slopes below unity when the line luminosity depends upon the average gas density with a power-law index greater than a Kennicutt-Schmidt index of 1.5.

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

We derive a physical model for the observed relations between star formation rate (SFR) and molecular line (CO and HCN) emission in galaxies, and show how these observed relations are reflective of the underlying star formation law. We do this by combining 3D non-LTE radiative transfer calculations with hydrodynamic simulations of isolated disk galaxies and galaxy mergers. We demonstrate that the observed SFR-molecular line relations are driven by the relationship between molecular line emission and gas density, and anchored by the index of the underlying Schmidt law controlling the SFR in the galaxy. Lines with low critical densities (e.g. CO J=1-0) are typically thermalized and trace the gas density faithfully. In these cases, the SFR will be related to line luminosity with an index similar to the Schmidt law index. Lines with high critical densities greater than the mean density of most of the emitting clouds in a galaxy (e.g. CO J=3-2, HCN J=1-0) will have only a small amount of thermalized gas, and consequently a superlinear relationship between molecular line luminosity and mean gas density. This results in a SFR-line luminosity index less than the Schmidt index for high critical density tracers. One observational consequence of this is a significant redistribution of light from the small pockets of dense, thermalized gas to diffuse gas along the line of sight, and prodigious emission from subthermally excited gas. At the highest star formation rates, the SFR-Lmol slope tends to the Schmidt index, regardless of the molecular transition. The fundamental relation is the Kennicutt-Schmidt law, rather than the relation between SFR and molecular line luminosity. We use these results to make imminently testable predictions for the SFR-molecular line relations of unobserved transitions.

I review the basic processes that may be used to develop a chemical evolutionary sequence for low-mass starless cores. I highlight observational results from the Arizona Radio Observatory-Green Bank Survey. Observations were performed with the SMT 10-m, ARO 12-m, and GBT 100-m toward a sample of 25 nearby (D < 400 pc) low-mass starless cores which have radiative transfer models of the 850 $\mu$m emission and observed SED (160 - 1300 um). The cores were observed in the lines of NH3 (1,1) and (2,2), o-NH2D 1_{11} - 1_{01}, C2S 1_2 - 2_1, C3S 4 - 3, HCN 1 - 0, HC5N 9 - 8, HC7N 21 - 20, C18O and C17O 2 - 1, and p-H2CO 1_{01} - 0_{00}.

We review the evidence for a constant star formation rate per unit mass in dense molecular gas in the Milky Way and the extragalactic correlations of L_IR with L' from observations of dense molecular gas. We discuss the connection between the constant SFR/M interpretation in dense gas and the global Schmidt-Kennicutt star formation law.

We present Spitzer observations of the dark globule CB190 (L771). We observe a roughly circular 24 micron shadow with a 70 arcsec radius. The extinction profile of this shadow matches the profile derived from 2MASS photometry at the outer edges of the globule and reaches a maximum of ~32 visual magnitudes at the center. The corresponding mass of CB190 is ~10 Msun. Our 12CO and 13CO J = 2-1 data over a 10 arcmin X 10 arcmin region centered on the shadow show a temperature ~10 K. The thermal continuum indicates a similar temperature for the dust. The molecular data also show evidence of freezeout onto dust grains. We estimate a distance to CB190 of 400 pc using the spectroscopic parallax of a star associated with the globule. Bonnor-Ebert fits to the density profile, in conjunction with this distance, yield xi_max = 7.2, indicating that CB190 may be unstable. The high temperature (56 K) of the best fit Bonnor-Ebert model is in contradiction with the CO and thermal continuum data, leading to the conclusion that the thermal pressure is not enough to prevent free-fall collapse. We also find that the turbulence in the cloud is inadequate to support it. However, the cloud may be supported by the magnetic field, if this field is at the average level for dark globules. Since the magnetic field will eventually leak out through ambipolar diffusion, it is likely that CB190 is collapsing or in a late pre-collapse stage.

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.

We present results from SMA observations of the star forming region Cepheus A-East at ~340 GHz (875 micron) with 0.7'' - 2'' resolution. At least four compact submm continuum sources have been detected, as well as a rich forest of hot core line emission. Two kinematically, chemically, and thermally distinct regions of molecular emission are present in the vicinity of the HW2 thermal jet, both spatially distinct from the submm counterpart to HW2. We propose that this emission is indicative of multiple protostars rather than a massive disk as reported by Patel et al. (2005).

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.

Low-mass dense cores represent the state of molecular gas associated with the earliest phases of low-mass star formation. Such cores are called "protostellar" or "starless," depending on whether they do or do not contain compact sources of luminosity. In this chapter, the first half of the review of low-mass dense cores, we describe the numerous inferences made about the nature of starless cores as a result of recent observations, since these reveal the initial conditions of star formation. We focus on the identification of isolated starless cores and their internal physical and chemical properties, including morphologies, densities, temperatures, kinematics, and molecular abundances. These objects display a wide range of properties since they are each at different points on evolutionary paths from ambient molecular cloud material to cold, contracting, and centrally concentrated configurations with significant molecular depletions and, in rare cases, enhancements.

Observations have revealed prodigious amounts of star formation in starburst galaxies as traced by dust and molecular emission, even at large redshifts. Recent work shows that for both nearby spiral galaxies and distant starbursts, the global star formation rate, as indicated by the infrared luminosity, has a tight and almost linear correlation with the amount of dense gas as traced by the luminosity of HCN. Our surveys of Galactic dense cores in HCN 1-0 emission show that this correlation continues to a much smaller scale, with nearly the same ratio of infrared luminosity to HCN luminosity found over 7-8 orders of magnitude in L_IR, with a lower cutoff around 10^{4.5} L_sun of infrared luminosity. The linear correlation suggests that we may understand distant star formation in terms of the known properties of local star-forming regions. Both the correlation and the luminosity cutoff can be explained if the basic unit of star formation in galaxies is a dense core, similar to those studied in our Galaxy.

We present observations of IRAS 16293-2422 using the Submillimeter Array (SMA) at 300 GHz, and the Very Large Array (VLA) at frequencies from 1.5 to 43 GHz. This source was already known to be a binary system with components A and B separated by ~5". The new SMA data now separate source A into two submm continuum components, denoted Aa and Ab. The strongest of these, Aa, peaks between the cm radio sources A1 and A2. Archival VLA data show proper motion of sources A and B of 17 mas/yr, associated with motion of the rho Ophiuchi cloud. We also find relative motion between A1 and A2 which excludes the possibility that they are gravitationally bound. A2 remains stationary relative to source B, and we identify it as the protostar driving the NE-SW CO outflow. A1 is shock-ionized gas where a precessing jet interacts with dense gas. This jet probably drives the E-W CO outflow. The origin of this jet must be located close to A2, and may be the submm source Aa. Thus source A is now shown to comprise three components within 1". Source B, however, is single, exhibits optically-thick dust emission even at 8 GHz, has a high luminosity, and yet shows no sign of outflow, indicating it is probably very young. The SMA spectrum of IRAS 16293-2422 reports the first astronomical identification of many lines at 300 and 310 GHz. The species detected are typical of hot cores, and are mainly associated with source A. We suggest the shock at A1 may explain the enhanced abundances measured, especially of sulphur-bearing molecules. We use the high brightness temperature of the submm dust emission from source B to demonstrate the unambiguous detection of infall by observing redshifted SO (7_7-6_6) absorption against the dust disk.

Lynds 1498 is a pre-protostellar core (PPC) and was one of the initial objects toward which molecular depletion and differentiation was detected. Despite the considerable scrutiny of L1498, there has not been a extensive study of the density and temperature structure as derived from radiative transfer modeling of dust continuum observations. We present deep SCUBA observations of L1498 at 850 and 450 micron, high resolution BEARS maps of the N2H+ 1-0 transition, CSO observations of the N2H+ 3-2 transition, and GBT observations of the C3S 4-3 transition. We also present a comparison of derived properties between L1498 and nearby PPCs that have been observed at far-infrared and submillimeter wavelengths. We present a more realistic treatment of PPC heating which varies the strength of the ISRF, Sisrf, and includes attenuation of the ISRF due to dust grains at the outer radius of the core, Av. The best-fitted model consists of a Bonner-Ebert sphere with a central density of 1 - 3 x 10^4 cm-3, R_o ~ 0.29 pc, 0.5 <= Sisrf <= 1, Av ~ 1 mag, and a nearly isothermal temperature profile of ~ 10.5 K for OH8 opacities. C3S emission shows a central depletion hole while N2H+ emission is centrally peaked. The observed depletions of C3S and H2CO, the modest N2H+ abundance, and a central density that is an order of magnitude lower than other modeled PPCs suggests that L1498 may be a forming PPC. Our derived temperature and density profile will improve modeling of molecular line observations that will explicate the core's kinematical and chemical state. (abridged)

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.

A series of self-consistent, three-dimensional continuum radiative transfer models are constructed of the pre-protostellar core L1544, with the results compared with existing SCUBA and ISO data. The source is well-fit by a prolate spheroid, having an ellipsoidal power-law density distribution of index m ~ 2 (1.75 < m < 2.25) in to at least r ~ 1600AU. For r<1600 AU, the data are consistent with either an extension of the power law to smaller radii, or a flattened (Bonner-Ebert like) density distribtion. We can further constrain the optical depth along the short axis at 1300um to be ~ 5e-3, the central luminosity to be L < 1e-3 solar luminosities, the long axis diameter D ~ 0.1 pc, the axis ratio to be q ~ 2, and the external ISRF to be similar to that defined by Mathis, Mezger, & Panagia (1983) to within 50 per cent. The outer diameter and axis ratio may each be somewhat larger due to potential on-source chopping in the observations, and the projection of the long axis onto the plane of the sky. While these results are similar to those inferred directly from observations or spherical modeling due to the source transparency at submillimeter wavelengths, we infer a smaller size, lower mass, and higher optical depth / column density, exposed to a stronger external radiation field than previously assumed. Finally, we find that both the spectral energy distribution (SED) and surface brightness distribution are necessary to constrain the source properties in this way.

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 observed the massive star forming region associated with the early B protostar G192.16-3.84 in NH3(1,1), 22.2 GHz water masers, 1.3 cm continuum emission, and at 850 microns. The dense gas associated with G192.16 is clumpy, optically thin, and has a mass of 0.9 Msun. The ammonia core is gravitationally unstable which may signal that the outflow phase of this system is coming to an end. Water masers trace an ionized jet 0.8" (1600 AU at a distance of 2 kpc) north of G192.16. Masers are also located within 500 AU of G192.16, their velocity distribution is consistent with but does not strongly support the interpretation that the maser emission arises in a 1000 AU rotating disk centered on G192.16. Roughly 30" south of G192.16 (0.3 pc) is a compact, optically thick (optical depth = 1.2) ammonia core (called G192 S3) with an estimated mass of 2.6 Msun. Based on the presence of 850 micron and 1.2 mm continuum emission, G192 S3 probably harbors a very young, low-mass protostar or proto-cluster. The dense gas in the G192 S3 core is likely to be gravitationally bound and may represent the next site of star formation in this region.

We have mapped 63 regions forming high-mass stars in CS J=5-4 using the CSO. The CS peak position was observed in C34S J=5-4 towards 57 cores and in 13CS J=5-4 towards the 9 brightest cores. The sample is a subset of a sample originally selected toward water masers; the selection on maser sources should favor sources in an early stage of evolution. The integrated intensity of CS J=5-4 correlates very well with the dust continuum emission at 350 microns. The distributions of size, virial mass, surface density, and luminosity are all peaked with a few cores skewed towards much larger values than the mean. We find a weak correlation between C34S linewidth and size, consistent with Dv ~ R^{0.3}. The linewidths are much higher than would be predicted by the usual relations between linewidth and size determined from regions of lower mass. These regions are very turbulent. The derived virial mass agrees within a factor of 2 to 3 with mass estimates from dust emission at 350 microns after corrections for the density structure are accounted for. The resulting cumulative mass spectrum of cores above 1000 solar masses can be approximated by a power law with a slope of about -0.9, steeper than that of clouds measured with tracers of lower density gas and close to that for the total masses of stars in OB associations. The median turbulent pressures are comparable to those in UCHII regions, and the pressures at small radii are similar to those in hypercompact-HII regions (P/k ~ 10^{10} K cm^{-3}). The filling factors for dense gas are substantial, and the median abundance of CS is about 10^{-9}. The ratio of bolometric luminosity to virial mass is much higher than the value found for molecular clouds as a whole, and the correlation of luminosity with mass is tighter. (Abridged).

We have compared the intensity distribution of molecular line emission with that of dust continuum emission, and modeled molecular line profiles in three different preprotostellar cores in order to test how dynamical evolution is related to chemical evolution, and whether we can use different chemical tracers to identify specific dynamical evolutionary stages. We used dust continuum emission to obtain the input density and temperature structures by calculating radiative transfer of dust emission. Our results show that chemical evolution is dependent on dynamical processes, which can give different evolutionary timescales, as well as the density structure of the core.

We have compared molecular line emission to dust continuum emission and modeled molecular lines using Monte Carlo simulations in order to study the depletion of molecules and the ionization fraction in three preprotostellar cores, L1512, L1544, and L1689B. L1512 is much less dense than L1544 and L1689B, which have similar density structures. L1689B has a different environment from those of L1512 and L1544. We used density and temperature profiles, calculated by modeling dust continuum emission in the submillimeter, for modeling molecular line profiles. In addition, we have used molecular line profiles and maps observed in several different molecules toward the three cores. We find a considerable diversity in chemical state among the three cores. The molecules include those sensitive to different timescales of chemical evolution such as CCS, the isotopes of CO and HCO+, DCO+, and N2H+. The CO molecule is significantly depleted in L1512 and L1544, but not in L1689B. CCS may be in the second enhancement of its abundance in L1512 and L1544 because of the significant depletion of CO molecules. N2H+ might already start to be depleted in L1512, but it traces very well the distribution of dust emission in L1544. On the other hand, L1689B may be so young that N2H+ has not reached its maximum yet. The ionization fraction has been calculated using H13CO+ and DCO+. This study suggests that chemical evolution depends on the absolute timescale during which a core stays in a given environment as well as its density structure.

We present results from the observations and modeling of seventeen Class I cores with the Submillimetre Common Users Bolometer Array (SCUBA) on the James Clerk Maxwell Telescope (JCMT). By modeling the transfer of radiation through the envelope for nine cores, we find, for a power law distribution n(r)=n_f(r/r_f)^-p, the average and standard deviation p=1.6 +/- 0.4 and a median of p=1.8. However, the inclusion of a disk or other point-like component can cause the derived p to be shallower by as much as 0.5. In addition, we test the Shu collapse model for our sources and discuss the application of simpler analyses that derive a density power law distribution directly from the slope of the intensity radial profile. The total mass of the envelope in our sample has a range from 0.04 to 5.0 M_sun, but these masses disagree with the virial masses derived from molecular line observations. Finally, we discuss the nature of these sources in light of various evolutionary indicators and find that T_bol and L_obs/L_smm are often inconsistent in distinguishing Class 0 from Class I cores.

Fifty-one dense cores associated with water masers were mapped at 350 micron. These cores are very luminous, 10^3 < Lbol/Lsun < 10^6, indicative of the formation of massive stars. Dust continuum contour maps and photometry are presented for these sources. The spectral energy distributions and normalized radial profiles of dust continuum emission were modeled for 31 sources using a one-dimensional dust radiative transfer code, assuming a power law density distribution in the envelope, n = n_f (r/r_f)^{-p}. The best fit density power law exponent, p, ranged from 0.75 to 2.5 with <p> = 1.8 +/- 0.4. The mean value of p is comparable to that found in regions forming only low mass stars. The mean p is incompatible with a logatropic sphere (p = 1), but other star formation models cannot be ruled out. Different mass estimates are compared and mean masses of gas and dust are reported within a half-power radius determined from the dust emission and within a radius where the total density exceeds 10^4 cm^3. Evolutionary indicators commonly used for low mass star formation may have some utility for regions forming massive stars. For comparison with extragalactic star formation studies, the luminosity to dust mass ratio is calculated for these sources with a method most parallel to that used in studies of distant galaxies and is found to be similar to that seen in high redshift starburst galaxies.

Studying the physical environments of low mass and high mass cores using dust continuum emission provides important observational constraints on theoretical models of star formation. The motivation and procedure for modeling dust continuum emission is reviewed and the results of recent surveys towards low mass and high mass star forming regions are compared.

Seven Class 0 sources mapped with SCUBA at 850 and 450 micron are modeled using a one dimensional radiative transfer code. The modeling takes into account heating from an internal protostar, heating from the ISRF, realistic beam effects, and chopping to model the normalized intensity profile and spectral energy distribution. Power law density models, n(r) ~ r^{-p}, fit all of the sources; best fit values are mostly p = 1.8 +/- 0.1, but two sources with aspherical emission contours have lower values (p ~ 1.1). Including all sources, <p> = 1.63 +/- 0.33. Based on studies of the sensitivity of the best-fit p to variations in other input parameters, uncertainties in p for an envelope model are \Delta p = +/- 0.2. If an unresolved source (e.g., a disk) contributes 70% of the flux at the peak, p is lowered in this extreme case and \Delta p = ^{+0.2}_{-0.6}. The models allow a determination of the internal luminosity (<L_{int}> = 4.0 \lsun) of the central protostar as well as a characteristic dust temperature for mass determination (<T_{iso}> = 13.8 +/- 2.4 K). We find that heating from the ISRF strongly affects the shape of the dust temperature profile and the normalized intensity profile, but does not contribute strongly to the overall bolometric luminosity of Class 0 sources. There is little evidence for variation in the dust opacity as a function of distance from the central source. The data are well-fitted by dust opacities for coagulated dust grains with ice mantles (Ossenkopf & Henning 1994). The density profile from an inside-out collapse model (Shu 1977) does not fit the data well, unless the infall radius is set so small as to make the density nearly a power-law.

We have observed 51 high mass star forming cores associated with water masers at 350 micron. The spectral energy distributions (SEDs) and dust continuum normalized radial intensity profiles were modeled for 28 sources using a one-dimensional dust radiative transfer code assuming a power law density distribution in the envelope n = n_0 (r/r_0)^{-p}. The best fit density power law exponent, p, ranged from 1.0 to 2.5 with a mean value of p = 1.72. We determined the dust masses for the modeled cores and found a mean mass of 209 M_{sun}.

We present the results of a multi-faceted study of the high-mass star-forming region S235. With the CalTech Submillimeter Observatory (CSO), we have observed and/or mapped over fifty molecular transitions and 350 micron dust emission in an attempt to learn more about the chemistry, kinematics, and mass distribution within this rich region. We have found two dust cores that are not resolved in our molecular line observations and suggest that one of these is a less evolved star-forming core.

Twenty-four cores have been mapped in CS J = 7 -> 6 at the CSO. From the spectra we determine core sizes and virial masses. Combining results from the CS and dust continuum studies for M8E, we use Monte Carlo simulations for the CS emission to determine radial profiles.

Our knowledge of the initial conditions and early stages of high mass star formation is very limited. We will review recent surveys of regions in the early stages of massive star formation using molecular tracers of high density and dust continuum emission and consider the status of evolutionary schemes. Comparison to the situation for low mass, relatively isolated star formation will be used to illustrate the outstanding issues in massive star formation. The problem of initial conditions is particularly acute because there is a lack of observational evidence for regions capable of forming massive stars BEFORE star formation actually begins. By analogy with the Pre-Protostellar Cores (PPCs) studied for low-mass star formation, one might call such regions Pre-Proto-cluster Cores (PPclCs). We will conclude with some speculation about what such cores might look like and possibilities for their detection.

In this survey we have systematically mapped 63 cloud cores in the CS J = 5 -> 4 line towards a sample of high mass star forming cores with water masers (Plume et al. 1992, Plume et al. 1997) using the CSO. From the CS spectra and maps we determine cloud core sizes, virial masses, a mass spectrum, and a size-linewidth relationship.

We have modeled the emission from dust in pre-protostellar cores, including a self-consistent calculation of the temperature distribution for each input density distribution. Model density distributions include Bonnor-Ebert spheres and power laws. The Bonnor-Ebert spheres fit the data well for all three cores we have modeled. The dust temperatures decline to very low values ($\Td \sim 7$ K) in the centers of these cores, strongly affecting the dust emission. Compared to earlier models that assume constant dust temperatures, our models indicate higher central densities and smaller regions of relatively constant density. Indeed, for L1544, a power-law density distribution, similar to that of a singular, isothermal sphere, cannot be ruled out. For the three sources modeled herein, there seems to be a sequence of increasing central condensation, from L1512 to L1689B to L1544. The two denser cores, L1689B and L1544, have spectroscopic evidence for contraction, suggesting an evolutionary sequence for pre-protostellar cores.

We have obtained 850 and 450 micron continuum maps of 21 low mass cores with SED's ranging from pre-protostellar to Class I (18K < Tbol < 370K), using SCUBA at the JCMT. In this paper we present the maps, radial intensity profiles, and photometry. Pre-protostellar cores do not have power-law intensity profiles, whereas the intensity profiles of Class 0 and Class I sources can be fitted with power laws over a large range of radii. A substantial number of sources have companion sources within a few arcminutes (2 out of 5 pre-protostellar cores, 9 out of 16 Class 0/I sources). The mean separation between sources is 10800 AU. The median separation is 18000 AU including sources without companions as a lower limit. The mean value of the spectral index between 450 and 850 microns is 2.8 +- 0.4, with PPCs having slightly lower spectral indices (2.5 +- 0.4). The mean mass of the sample, based on the dust emission in a 120 arcsecond aperture, is 1.1 +- 0.9 solar masses. For the sources fitted by power-law intensity distributions (I_nu(b)/I_nu(0) = (b/b_0) ^ m), the mean value of m is 1.52 +- 0.45 for Class 0 and I sources at 850 microns and 1.44 +- 0.25 at 450 microns. Based on a simple analysis, assuming the emission is in the Rayleigh-Jeans limit and that T_d(r) ~ r ^ {-0.4}, these values of m translate into power-law density distributions (n ~ r ^ {-p}) with p ~ 2.1. However, we show that this result may be changed by more careful consideration of effects such as beam size and shape, finite outer radii, more realistic T_d(r), and failure of the Rayleigh-Jeans approximation.