Search

Understanding the chemical processes during starless core and prestellar core evolution is an important step in understanding the initial stages of star and disk formation. This project is a study of deuterated ammonia, o-NH$_2$D, in the L1251 star-forming region toward Cepheus. Twenty-two dense cores (twenty of which are starless or prestellar, and two of which have a protostar), previously identified by p-NH$_3$ (1,1) observations, were targeted with the 12m Arizona Radio Observatory telescope on Kitt Peak. o-NH$_2$D J$_{\rm{K_a} \rm{K_c}}^{\pm} =$ $1_{11}^{+} \rightarrow 1_{01}^{-}$ was detected in 13 (59\%) of the NH$_3$-detected cores with a median sensitivity of $\sigma_{T_{mb}} = 17$ mK. All cores detected in o-NH$_2$D at this sensitivity have p-NH$_3$ column densities $> 10^{14}$ cm$^{-2}$. The o-NH$_2$D column densities were calculated using the constant excitation temperature (CTEX) approximation while correcting for the filling fraction of the NH$_3$ source size. The median deuterium fraction was found to be 0.11 (including 3$\sigma$ upper limits). However, there are no strong, discernible trends in plots of deuterium fraction with any physical or evolutionary variables. If the cores in L1251 have similar initial chemical conditions, then this result is evidence of the cores physically evolving at different rates., Comment: 18 pages, 14 figures, MNRAS, in press

Chemistry along the star- and planet-formation sequence regulates how prebiotic building blocks -- carriers of the elements CHNOPS -- are incorporated into nascent planetesimals and planets. Spectral line observations across the electromagnetic spectrum are needed to fully characterize interstellar CHNOPS chemistry, yet to date there are only limited astrochemical constraints at THz frequencies. Here, we highlight advances to the study of CHNOPS astrochemistry that will be possible with the Orbiting Astronomical Satellite for Investigating Stellar Systems (OASIS). OASIS is a NASA mission concept for a space-based observatory that will utilize an inflatable 14-m reflector along with a heterodyne receiver system to observe at THz frequencies with unprecedented sensitivity and angular resolution. As part of a survey of H2O and HD towards ~100 protostellar and protoplanetary disk systems, OASIS will also obtain statistical constraints on the inventories of light hydrides including NH3 and H2S towards protoplanetary disks, as well as complex organics in protostellar hot corinos and envelopes. Line surveys of additional star-forming regions, including high-mass hot cores, protostellar outflow shocks, and prestellar cores, will also leverage the unique capabilities of OASIS to probe high-excitation organics and small hydrides, as is needed to fully understand the chemistry of these objects., Comment: Accepted to Frontiers in Astronomy and Space Sciences

The physical processes behind the transfer of mass from parsec-scale clumps to massive-star-forming cores remain elusive. We investigate the relation between the clump morphology and the mass fraction that ends up in its most massive core (MMC) as a function of infrared brightness, i.e. a clump evolutionary tracer. Using ALMA 12 m and ACA we surveyed 6 infrared-dark hubs in 2.9mm continuum at $\sim$3" resolution. To put our sample into context, we also re-analysed published ALMA data from a sample of 29 high mass-surface density ATLASGAL sources. We characterise the size, mass, morphology, and infrared brightness of the clumps using Herschel and Spitzer data. Within the 6 newly observed hubs, we identify 67 cores, and find that the MMCs have masses between 15-911 $\mathrm{M}_{\odot}$ within a radius of 0.018-0.156 pc. The MMC of each hub contains 3-24% of the clump mass ($f_\mathrm{MMC}$), becoming 5-36% once core masses are normalised to the median core radius. Across the 35 clumps, we find no significant difference in the median $f_\mathrm{MMC}$ values of hub and non-hub systems, likely the consequence of a sample bias. However, we find that $f_\mathrm{MMC}$ is $\sim$7.9 times larger for infrared-dark clumps compared to infrared-bright ones. This factor increases up to $\sim$14.5 when comparing our sample of 6 infrared-dark hubs to infrared-bright clumps. We speculate that hub-filament systems efficiently concentrate mass within their MMC early on during its evolution. As clumps evolve, they grow in mass, but such growth does not lead to the formation of more massive MMCs., Comment: 15 pages, 9 figures

Recent observations indicate that organic molecules are prevalent towards starless and prestellar cores. Deuteration of these molecules has not been well-studied during the starless phase. Published observations of singly-deuterated methanol, CH$_2$DOH, have only been observed in a couple of well-studied, dense and evolved prestellar cores (e.g. L1544, L183). Since the formation of gas-phase methanol during this cold phase is believed to occur via desorption from the icy grain surfaces, observations of CH$_2$DOH may be useful as a probe of the deuterium fraction in the ice mantles of dust grains. We present a systematic survey of CH$_2$DOH towards 12 starless and prestellar cores in the B10 region of the Taurus Molecular Cloud. Nine of the twelve cores are detected with [CH$_2$DOH]/[CH$_3$OH] ranging from $< 0.04$ to $0.23^{+0.12}_{-0.06}$ with a median value of $0.11$. Sources not detected tend to have larger virial parameters and larger methanol linewidths than detected sources. The results of this survey indicate that deuterium fractionation of organic molecules, such as methanol, during the starless phase may be more easily detectable than previously thought., Comment: To be published in MNRAS

(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., Comment: 32 pages, 17 figures, and 6 tables. Accepted for publication in ApJ (August 26, 2019)

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)., Comment: ApJ accepted

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., Comment: accepted for publication in ApJ. 34 pages

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., Comment: Accepted to The Astrophysical Journal (ApJ), 45 pages, 24 figures (some with reduced resolution in this preprint); Project website is at http://carma.astro.umd.edu/classy

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., Comment: to appear in Astronomy & Astrophysics

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., Comment: Accepted for publication in ApJ; 33 pages; 22 figures; 7 tables

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%., Comment: Accepted for publication in ApJ; 32 pages, 21 figures

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., Comment: ApJ Accepted

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., Comment: 76 pages, 26 figures, 6 tables. Fig. 25 shows the complete sample of 107 maps. Accepted for publication in ApJS. Full tables 3 and 5 will be available online through ApJS. Typos corrected and other small changes on text after journal proof-review

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)., Comment: Tutorial article accepted for publication in PASP. Table 1 lists critical densities and effective excitation densities for 72 transitions of 12 species. 24 pages

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., Comment: 35 pages with 15 figures and 20 tables. Appeared in the March 2015 issue of PASP. Latest version corrects a factor-of-two error in Eq. 100 (corrigendum to appear in PASP) and previous corrigendum (February 2016 issue of PASP): (1) Fix error in dipole moment units; (2) Fixed errors in constants for worked column density equations for C18O 1-0 (Eq. 90), C17O 1-0 (Eq. 93), and N2H+ 1-0 (Eq. 96)

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., Comment: 28 pages, 19 figures. Accepted for publication in ApJ. Distance-Omnibus code available at https://github.com/BGPS/distance-omnibus

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., Comment: Accepted to ApJ. 38 pages, 16 figures, 5 tables

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., Comment: 67 pages, 22 figures, 4 table. Accepted by The Astronomical Journal

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%)., Comment: Accepted for publication in ApJS. The full velocity catalog and data tables (complete Tables 2,3, and 5) may be downloaded from http://eldora.as.arizona.edu/~yshirley/Arizona/BGPS/ Kinematic finder charts (0.5 deg x 0.5 deg) are also available at the same address

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., Comment: Accepted for publication in ApJ; 29 pages, 16 figures. Data available prior to publication by contacting the authors

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 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., Comment: 42 pages, 19 figures, Accepted for publication in the Astronomical Journal

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., Comment: 13 pages

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., Comment: 38 pages (with full table), 3 figures

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., Comment: 15 pages, 4 figures, accepted for publication in ApJ

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., Comment: 48 pages, 25 figures, Accepted for publication

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, Comment: 25 pages, 16 Figures, Accepted to The Astrophysical Journal

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., Comment: 47 pages, 19 figures, accepted in ApJ on Jan 4th 2010

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., Comment: 28 pages, 7 figure. Accepted for publication in ApJ

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., Comment: 19 pages, 3 figures, AJ accepted, in press

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], Comment: 76 pages, 22 figures, published in ApJ: http://iopscience.iop.org/0004-637X/721/1/137/

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., Comment: 33 pages, 26 figures, 6 tables, accepted for publication in ApJ. Full resolution version available here: http://casa.colorado.edu/~battersb/Publications.html

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., Comment: 26 pages, 19 figures, 8 tables. Accepted for publication by ApJ

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., Comment: 41 pages, 28 figures, 5 tables; accepted by ApJ

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., Comment: accepted for publication in ApJ

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., Comment: 5 pages, 3 figures, emulateapj format, accepted for publication in ApJL

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., Comment: 15 pages, 17 figures, emulateapj format, accepted for publication in ApJ

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., Comment: 7 pages. To appear in Massive Star Formation: Observations Confront Theory, ed. H. Beuther, ASP Conference Series

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}., Comment: To appear in "Molecules in Space & Laboratory", Paris 2007, eds. J.L. Lemaire & F. Combes. Corrected typo, otherwise "more or less" the same

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., Comment: 16 pages, 13 figures, accepted for publication in ApJ

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., Comment: 14 pages. Accepted for publication in ApJ

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., Comment: 10 pages with 3 figures, accepted by ApJ Letters

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., Comment: 10 pages, 2 figures. In press for ApJ Letters

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., Comment: 60 pages, 22 figures, accepted for publication by ApJ

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), Comment: 26 pages, 10 figures (2 color figs). Accepted to ApJ

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., Comment: 9 pages; 8 figures; accepted for publication in MNRAS

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)., Comment: Figure converted to .gif format for astro-ph. Please download high-res version from http://www.aoc.nrao.edu/~yshirley/papers.html

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., Comment: 3 pages, to appear in the proceedings of the conference, "Star Formation as a Diagnostic of Star Formation"

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 = 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., Comment: 45 pages, 20 figures, accepted to ApJ Supplement

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., Comment: 5 pages, 3 figures. To appear in Galactic Star Formation Across the Stellar Mass Spectrum, ASP conf. series, ed. James M. De Buizer