Showing 1–50 of 73 results for author: Shirley, Y L

A Survey of Deuterated Ammonia in the Cepheus Star-Forming Region L1251

Authors: Maria Galloway-Sprietsma, Yancy L. Shirley, James Di Francesco, Jared Keown, Samantha Scibelli, Olli Sipilä, Rachel Smullen

Abstract: 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 ident… ▽ More 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 $σ_{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$σ$ 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. △ Less

Submitted 21 July, 2022; originally announced July 2022.

Comments: 18 pages, 14 figures, MNRAS, in press

Astrochemistry with the Orbiting Astronomical Satellite for Investigating Stellar Systems (OASIS)

Authors: Jennifer B. Bergner, Yancy L. Shirley, Jes K. Jorgensen, Brett McGuire, Susanne Aalto, Carrie M. Anderson, Gordon Chin, Maryvonne Gerin, Paul Hartogh, Daewook Kim, David Leisawitz, Joan Najita, Kamber R. Schwarz, Alexander G. G. M. Tielens, Christopher K. Walker, David J. Wilner, Edward J. Wollack

Abstract: 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. H… ▽ More 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. △ Less

Submitted 9 December, 2021; v1 submitted 14 November, 2021; originally announced November 2021.

Comments: Accepted to Frontiers in Astronomy and Space Sciences

An ALMA study of hub-filament systems I. On the clump mass concentration within the most massive cores

Authors: Michael Anderson, Nicolas Peretto, Sarah E. Ragan, Andrew J. Rigby, Adam Avison, Ana Duarte-Cabral, Gary A. Fuller, Yancy L. Shirley, Alessio Traficante, Gwenllian M. Williams

Abstract: 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… ▽ More 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. △ Less

Submitted 15 September, 2021; originally announced September 2021.

Comments: 15 pages, 9 figures

A Survey of CH2DOH Towards Starless and Prestellar Cores in the Taurus Molecular Cloud

Authors: Hannah E. Ambrose, Yancy L. Shirley, Samantha Scibelli

Abstract: 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 metha… ▽ More 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. △ Less

Submitted 17 November, 2020; originally announced November 2020.

Comments: To be published in MNRAS

ALMA observations of fragmentation, sub-structure, and protostars in high-mass starless clump candidates

Authors: Brian E. Svoboda, Yancy L. Shirley, Alessio Traficante, Cara Battersby, Gary A. Fuller, Qizhou Zhang, Henrik Beuther, Nicolas Peretto, Crystal Brogan, Todd Hunter

Abstract: (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… ▽ More (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σ$ (or $1σ$ 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. △ Less

Submitted 27 August, 2019; originally announced August 2019.

Comments: 32 pages, 17 figures, and 6 tables. Accepted for publication in ApJ (August 26, 2019)

An Ammonia Spectral Map of the L1495-B218 Filaments in the Taurus Molecular Cloud: II CCS & HC$_7$N Chemistry and Three Modes of Star Formation in the Filaments

Authors: Young Min Seo, Liton Majumdar, Paul F. Goldsmith, Yancy L. Shirley, Karen Willacy, Derek Ward-Thompson, Rachel Friesen, David Frayer, Sarah E. Church, Dongwoo Chung, Kieran Cleary, Nichol Cunningham, Kiruthika Devaraj, Dennis Egan, Todd Gaier, Rohit Gawande, Joshua O. Gundersen, Andrew I. Harris, Pekka Kangaslahti, Anthony C. S. Readhead, Lorene Samoska, Matthew Sieth, Michael Stennes, Patricia Voll, Steve White

Abstract: 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 ev… ▽ More 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). △ Less

Submitted 14 December, 2018; originally announced December 2018.

Comments: ApJ accepted

ALMA Observations of Starless Core Substructure in Ophiuchus

Authors: Helen Kirk, Michael M. Dunham, James Di Francesco, Doug Johnstone, Stella S. R. Offner, Sarah I. Sadavoy, John J. Tobin, Hector G. Arce, Tyler L. Bourke, Steve Mairs, Philip C. Myers, Jaime E. Pineda, Scott Schnee, Yancy L. Shirley

Abstract: 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.… ▽ More 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. △ Less

Submitted 1 March, 2017; originally announced March 2017.

Comments: accepted for publication in ApJ. 34 pages

CARMA Large Area Star Formation Survey: Dense Gas in the Young L1451 Region of Perseus

Authors: Shaye Storm, Lee G. Mundy, Katherine I. Lee, Manuel Fernández-López, Leslie W. Looney, Peter Teuben, Héctor G. Arce, Erik W. Rosolowsky, Aaron M. Meisner, Andrea Isella, Jens Kauffmann, Yancy L. Shirley, Woojin Kwon, Adele L. Plunkett, Marc W. Pound, Dominique M. Segura-Cox, Konstantinos Tassis, John J. Tobin, Nikolaus H. Volgenau, Richard M. Crutcher, Leonardo Testi

Abstract: 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 det… ▽ More 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. △ Less

Submitted 28 June, 2016; originally announced June 2016.

Comments: 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

Deuteration of ammonia in the starless core Ophiuchus/H-MM1

Authors: Jorma Harju, Fabien Daniel, Olli Sipilä, Paola Caselli, Jaime E. Pineda, Rachel K. Friesen, Anna Punanova, Rolf Güsten, Laurent Wiesenfeld, Philip C. Myers, Alexandre Faure, Pierre Hily-Blant, Claire Rist, Erik Rosolowsky, Stephan Schlemmer, Yancy L. Shirley

Abstract: 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… ▽ More 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. △ Less

Submitted 22 December, 2016; v1 submitted 19 April, 2016; originally announced April 2016.

Comments: to appear in Astronomy & Astrophysics

Journal ref: A&A 600, A61 (2017)

The Bolocam Galactic Plane Survey. XIV. Physical Properties of Massive Starless and Star Forming Clumps

Authors: Brian E Svoboda, Yancy L Shirley, Cara Battersby, Erik W Rosolowsky, Adam G Ginsburg, Timothy P Ellsworth-Bowers, Michele R Pestalozzi, Miranda K Dunham, Neal J Evans II, John Bally, Jason Glenn

Abstract: 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$ $μ{\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… ▽ More 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$ $μ{\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 ($α\sim 0.7$) with $>75\%$ of clumps with known distance being gravitationally bound ($α< 2$). These samples show a statistically significant increase in the median clump mass of $Δ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. △ Less

Submitted 27 November, 2015; originally announced November 2015.

Comments: Accepted for publication in ApJ; 33 pages; 22 figures; 7 tables

The Bolocam Galactic Plane Survey. XIII. Physical Properties and Mass Functions of Dense Molecular Cloud Structures

Authors: Timothy P. Ellsworth-Bowers, Jason Glenn, Allyssa Riley, Erik Rosolowsky, Adam Ginsburg, Neal J. Evans II, John Bally, Cara Battersby, Yancy L. Shirley, Manuel Merello

Abstract: 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 abo… ▽ More 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 $α= 2.0\pm0.1$ for sources nearer than 6.5 kpc and $α= 1.9\pm0.1$ for objects between 2 kpc and 10 kpc). The high-mass power-law indices are generally $1.85 \leq α\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α = 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%. △ Less

Submitted 6 April, 2015; originally announced April 2015.

Comments: Accepted for publication in ApJ; 32 pages, 21 figures

An Ammonia Spectral Map of the L1495-B218 Filaments in the Taurus Molecular Cloud : I. Physical Properties of Filaments and Dense cores

Authors: Young Min Seo, Yancy L. Shirley, Paul Goldsmith, Derek Ward-Thompson, Jason M. Kirk, Markus Schmalzl, Jeong-Eun Lee, Rachel Friesen, Glen Langston, Joe Masters, Robert W. Garwood

Abstract: 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… ▽ More 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 $μ$m and 500 $μ$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. △ Less

Submitted 17 March, 2015; v1 submitted 17 March, 2015; originally announced March 2015.

Comments: ApJ Accepted

The Bolocam Galactic Plane Survey. XI. Temperatures and Substructure of Galactic Clumps Based on 350 micron Observations

Authors: Manuel Merello, Neal J. Evans II, Yancy L. Shirley, Erik Rosolowsky, Adam Ginsburg, John Bally, Cara Battersby, Michael M. Dunham

Abstract: 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 pop… ▽ More 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. △ Less

Submitted 2 April, 2015; v1 submitted 23 January, 2015; originally announced January 2015.

Comments: 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

How to Calculate Molecular Column Density

Authors: Jeffrey G. Mangum, Yancy L. Shirley

Abstract: 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 expres… ▽ More 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. △ Less

Submitted 7 April, 2017; v1 submitted 7 January, 2015; originally announced January 2015.

Comments: 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)

The Critical Density and the Effective Excitation Density of Commonly Observed Molecular Dense Gas Tracers

Authors: Yancy L. Shirley

Abstract: The optically thin critical densities and the effective excitation densities to produce a 1 K km/s (or 0.818 Jy km/s $(\frac{ν_{jk}}{100 \rm{GHz}})^2 \, (\frac{θ_{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 temperatur… ▽ More The optically thin critical densities and the effective excitation densities to produce a 1 K km/s (or 0.818 Jy km/s $(\frac{ν_{jk}}{100 \rm{GHz}})^2 \, (\frac{θ_{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). △ Less

Submitted 7 January, 2015; originally announced January 2015.

Comments: 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 Bolocam Galactic Plane Survey. XII. Distance Catalog Expansion Using Kinematic Isolation of Dense Molecular Cloud Structures With 13CO(1-0)

Authors: Timothy P. Ellsworth-Bowers, Erik Rosolowsky, Jason Glenn, Adam Ginsburg, Neal J. Evans II, Cara Battersby, Yancy L. Shirley, Brian Svoboda

Abstract: 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 tr… ▽ More 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. △ Less

Submitted 10 November, 2014; originally announced November 2014.

Comments: 28 pages, 19 figures. Accepted for publication in ApJ. Distance-Omnibus code available at https://github.com/BGPS/distance-omnibus

Journal ref: ApJ, 2015, 799, 29

CARMA Large Area Star Formation Survey: Structure and Kinematics of Dense Gas in Serpens Main

Authors: Katherine I. Lee, Manuel Fernandez-Lopez, Shaye Storm, Leslie W. Looney, Lee G. Mundy, Dominique Segura-Cox, Peter Teuben, Erik Rosolowsky, Hector G. Arce, Eve C. Ostriker, Yancy L. Shirley, Woojin Kwon, Jens Kauffmann, John J. Tobin, Adele L. Plunkett, Marc W. Pound, Demerese M. Salter, N. H. Volgenau, Che-Yu Chen, Konstantinos Tassis, Andrea Isella, Richard M. Crutcher, Charles F. Gammie, Leonardo Testi

Abstract: 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 c… ▽ More 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. △ Less

Submitted 13 October, 2014; originally announced October 2014.

Comments: Accepted to ApJ. 38 pages, 16 figures, 5 tables

Outflows, Dusty Cores, and a Burst of Star Formation in the North America and Pelican Nebulae

Authors: John Bally, Adam Ginsburg, Ron Probst, Bo Reipurth, Yancy L. Shirley, Guy S. Stringfellow

Abstract: 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 /… ▽ More 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. △ Less

Submitted 23 September, 2014; originally announced September 2014.

Comments: 67 pages, 22 figures, 4 table. Accepted by The Astronomical Journal

CARMA Large Area Star Formation Survey: Project Overview with Analysis of Dense Gas Structure and Kinematics in Barnard 1

Authors: S. Storm, L. G. Mundy, M. Fernández-López, K. I. Lee, L. W. Looney, P. J. Teuben, E. Rosolowsky, H. G. Arce, E. C. Ostriker, D. Segura-Cox, M. W. Pound, D. M. Salter, N. H. Volgenau, Y. L. Shirley, C. Chen, H. Gong, A. L. Plunkett, J. J. Tobin, W. Kwon, A. Isella, J. Kauffmann, K. Tassis, R. M. Crutcher, C. F. Gammie, L. Testi

Abstract: 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 squa… ▽ More 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. △ Less

Submitted 3 September, 2014; originally announced September 2014.

Comments: Accepted to The Astrophysical Journal (ApJ), 51 pages, 27 figures (some with reduced resolution in this preprint); Project website is at http://carma.astro.umd.edu/classy

CARMA Large Area Star Formation Survey: Observational Analysis of Filaments in the Serpens South Molecular Cloud

Authors: M. Fernández-López, H. G. Arce, L. Looney, L. G. Mundy, S. Storm, P. J. Teuben, K. Lee, D. Segura-Cox, A. Isella, J. J. Tobin, E. Rosolowsky, A. Plunkett, W. Kwon, J. Kauffmann, E. Ostriker, K. Tassis, Y. L. Shirley, M. Pound

Abstract: 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… ▽ More 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. △ Less

Submitted 17 July, 2014; v1 submitted 2 July, 2014; originally announced July 2014.

Comments: 12 pages, 4 figures, published in ApJL (July 2014)

Journal ref: 2014, ApJ, 790, L19

Chemical evolution in the early phases of massive star formation. I

Authors: T. Gerner, H. Beuther, D. Semenov, H. Linz, T. Vasyunina, S. Bihr, Y. L. Shirley, Th. Henning

Abstract: 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 sta… ▽ More 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. △ Less

Submitted 18 March, 2014; v1 submitted 24 January, 2014; originally announced January 2014.

Comments: 31 pages, 11 figures, 21 tables, accepted by A&A; typos added

Journal ref: A&A 563, A97 (2014)

The Bolocam Galactic Plane Survey. X. A Complete Spectroscopic Catalog of Dense Molecular Gas Observed toward 1.1 mm Dust Continuum Sources with 7.5 <= l <= 194 degrees

Authors: Yancy L. Shirley, Timothy P. Ellsworth-Bowers, Brian Svoboda, Wayne M. Schlingman, Adam Ginsburg, Erik Rosolowsky, Thomas Gerner, Steven Mairs, Cara Battersby, Guy Stringfellow, Miranda K. Dunham, Jason Glenn, John Bally

Abstract: 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 large… ▽ More 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%). △ Less

Submitted 19 August, 2013; originally announced August 2013.

Comments: 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

The Bolocam Galactic Plane Survey. VIII. A Mid-Infrared Kinematic Distance Discrimination Method

Authors: Timothy P. Ellsworth-Bowers, Jason Glenn, Erik Rosolowsky, Steven Mairs, Neal J. Evans II, Cara Battersby, Adam Ginsburg, Yancy L. Shirley, John Bally

Abstract: 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 S… ▽ More 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. △ Less

Submitted 16 April, 2013; originally announced April 2013.

Comments: Accepted for publication in ApJ; 29 pages, 16 figures. Data available prior to publication by contacting the authors

Journal ref: ApJ, 2013, 770, 39

On the internal dynamics of starless cores: stability of starless cores with internal motions and collapse dynamics

Authors: Young Min Seo, Seung Soo Hong, Yancy L. Shirley

Abstract: 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 condit… ▽ More 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. △ Less

Submitted 10 April, 2013; originally announced April 2013.

The earliest phases of star formation - A Herschel key project. The thermal structure of low-mass molecular cloud cores

Authors: R. Launhardt, A. M. Stutz, A. Schmiedeke, Th. Henning, O. Krause, Z. Balog, H. Beuther, S. Birkmann, M. Hennemann, J. Kainulainen, T. Khanzadyan, H. Linz, N. Lippok, M. Nielbock, J. Pitann, S. Ragan, C. Risacher, M. Schmalzl, Y. L. Shirley, B. Stecklum, J. Steinacker, J. Tackenberg

Abstract: 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 constrai… ▽ More 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. △ Less

Submitted 8 January, 2013; originally announced January 2013.

Comments: 35 pages, 47 figures, Accepted for publication in Astronomy and Astrophysics

A Systematic Search for Molecular Outflows Toward Candidate Low-Luminosity Protostars and Very Low Luminosity Objects

Authors: Kamber R. Schwarz, Yancy L. Shirley, Michael M. Dunham

Abstract: 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 simul… ▽ More 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. △ Less

Submitted 11 September, 2012; originally announced September 2012.

Comments: 42 pages, 19 figures, Accepted for publication in the Astronomical Journal

Tentative Detection of Deuterated Methane toward the Low-Mass Protostar IRAS 04368+2557 in L1527

Authors: Nami Sakai, Yancy L. Shirley, Takeshi Sakai, Tomoya Hirota, Yoshimasa Watanabe, Satoshi Yamamoto

Abstract: 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… ▽ More 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. △ Less

Submitted 7 September, 2012; originally announced September 2012.

Comments: 13 pages

Testing the global star formation relation: An HCO+ (3-2) mapping study of Red MSX sources in the Bolocam Galactic Plane Survey

Authors: David E. Schenck, Yancy L. Shirley, Megan Reiter, Stephanie Juneau

Abstract: 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 H… ▽ More 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. △ Less

Submitted 6 September, 2011; originally announced September 2011.

Comments: 38 pages (with full table), 3 figures

Journal ref: The Astronomical Journal, 142, 3, 94 (2011)

Evidence for Inflow in High-Mass Star-Forming Clumps

Authors: Megan Reiter, Yancy L. Shirley, Jingwen Wu, Crystal Brogan, Alwyn Wootten, Ken'ichi Tatematsu

Abstract: 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 l… ▽ More 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. △ Less

Submitted 26 July, 2011; originally announced July 2011.

Comments: 15 pages, 4 figures, accepted for publication in ApJ

The Bolocam Galactic Plane Survey V: HCO+ and N2H+ Spectroscopy of 1.1 mm Dust Continuum Sources

Authors: Wayne M. Schlingman, Yancy L. Shirley, David E. Bolin, Erik Rosolowsky, John Bally, Cara Battersby, Miranda K. Dunham, Timothy P. Ellsworth-Bowers, Neal J. Evans II, Adam Ginsburg, Guy Stringfellow

Abstract: 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$σ$. We find a strong correlation between the integrated intensity of both dense gas tracers and the 1.1 mm dust emission of BGPS sourc… ▽ More 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$σ$. 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. △ Less

Submitted 20 May, 2011; originally announced May 2011.

Comments: 48 pages, 25 figures, Accepted for publication

C18O Depletion in Starless Cores in Taurus

Authors: Amanda Brady Ford, Yancy L. Shirley

Abstract: 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 c… ▽ More 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 △ Less

Submitted 12 January, 2011; originally announced January 2011.

Comments: 25 pages, 16 Figures, Accepted to The Astrophysical Journal

The Spitzer c2d Survey of Nearby Dense Cores XI : Infrared and Submillimeter Observations of CB130

Authors: Hyo Jeong Kim, Neal J. Evans II, Michael M. Dunham, Jo-Hsin Chen, Jeong-Eun Lee, Tyler L. Bourke, Tracy L. Huard, Yancy L. Shirley, Christopher De Vries

Abstract: 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… ▽ More 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. △ Less

Submitted 5 January, 2011; originally announced January 2011.

Comments: 47 pages, 19 figures, accepted in ApJ on Jan 4th 2010

Observational Constraints on Submillimeter Dust Opacity

Authors: Yancy L. Shirley, Tracy L. Huard, Klaus M. Pontoppidan, David J. Wilner, Amelia M. Stutz, John H. Bieging, Neal J. Evans II

Abstract: 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 (… ▽ More 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 ($κ\propto λ^{-β}$). 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{κ_{850}}{κ_{2.2}} = (3.21 - 4.80)^{+0.44}_{-0.30} \times 10^{-4}$ and $\frac{κ_{450}}{κ_{2.2}} = (12.8 - 24.8)^{+2.4}_{-1.3} \times 10^{-4}$ with a submillimeter opacity power-law index of $β_{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 $κ_{850} = (1.18 - 1.77)^{+0.36}_{-0.24}$ and $κ_{450} = (4.72 - 9.13)^{+1.9}_{-0.98}$ cm$^2$g$^{-1}$ of dust. These opacities are from $(0.65 - 0.97) κ^{\rm{OH}5}_{850}$ of the widely used theoretical opacities of Ossenkopf and Henning for coagulated ice grains with thin mantles at 850\micron. △ Less

Submitted 15 December, 2010; originally announced December 2010.

Comments: 28 pages, 7 figure. Accepted for publication in ApJ

MUSTANG 3.3 Millimeter Continuum Observations of Class 0 Protostars

Authors: Yancy L. Shirley, Brian S. Mason, Jeffrey G. Mangum, David E. Bolin, Mark J. Devlin, Simon R. Dicker, Phillip M. Korngut

Abstract: 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… ▽ More 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. △ Less

Submitted 16 November, 2010; originally announced November 2010.

Comments: 19 pages, 3 figures, AJ accepted, in press

The Bolocam Galactic Plane Survey IV: 1.1 and 0.35 mm Dust Continuum Emission in the Galactic Center Region

Authors: John Bally, James Aguirre, Cara Battersby, Eric Todd Bradley, Claudia Cyganowski, Darren Dowell, Meredith Drosback, Miranda K Dunham, Neal J. Evans II, Adam Ginsburg, Jason Glenn, Paul Harvey, Elisabeth Mills, Manuel Merello, Erik Rosolowsky, Wayne Schlingman, Yancy L. Shirley, Guy S. Stringfellow, Josh Walawender, Jonathan Williams

Abstract: 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… ▽ More 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 μ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 μ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 μ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] △ Less

Submitted 3 November, 2010; originally announced November 2010.

Comments: 76 pages, 22 figures, published in ApJ: http://iopscience.iop.org/0004-637X/721/1/137/

Journal ref: Astrophys.J.721:137-163,2010

An Infrared through Radio Study of the Properties and Evolution of IRDC Clumps

Authors: Cara Battersby, John Bally, James M. Jackson, Adam Ginsburg, Yancy L. Shirley, Wayne Schlingman, Jason Glenn

Abstract: 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 wi… ▽ More 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. △ Less

Submitted 4 August, 2010; originally announced August 2010.

Comments: 33 pages, 26 figures, 6 tables, accepted for publication in ApJ. Full resolution version available here: http://casa.colorado.edu/~battersb/Publications.html

The Initial Conditions of Clustered Star Formation III. The Deuterium Fractionation of the Ophiuchus B2 Core

Authors: R. K. Friesen, J. Di Francesco, P. C. Myers, A. Belloche, Y. L. Shirley, T. L. Bourke, P. André

Abstract: 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 seve… ▽ More 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. △ Less

Submitted 23 June, 2010; originally announced June 2010.

Comments: ApJ accepted

The Bolocam Galactic Plane Survey -- III. Characterizing Physical Properties of Massive Star-Forming Regions in the Gemini OB1 Molecular Cloud

Authors: Miranda K. Dunham, Erik Rosolowsky, Neal J. Evans II, Claudia J. Cyganowski, James Aguirre, John Bally, Cara Battersby, Eric Todd Bradley, Darren Dowell, Meredith Drosback, Adam Ginsburg, Jason Glenn, Paul Harvey, Manuel Merello, Wayne Schlingman, Yancy L. Shirley, Guy S. Stringfellow, Josh Walawender, Jonathan P. Williams

Abstract: 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 sourc… ▽ More 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. △ Less

Submitted 26 May, 2010; originally announced May 2010.

Comments: 26 pages, 19 figures, 8 tables. Accepted for publication by ApJ

Enhanced Dense Gas Fraction in Ultra-Luminous Infrared Galaxies

Authors: S. Juneau, D. T. Narayanan, J. Moustakas, Y. L. Shirley, R. S. Bussmann, R. C. Kennicutt Jr, P. A. Vanden Bout

Abstract: 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… ▽ More 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. △ Less

Submitted 3 December, 2009; v1 submitted 2 November, 2009; originally announced November 2009.

Comments: 20 pages, 10 figures. Published in The Astrophysical Journal (ApJ 707, 1217). Fixed a reference

Journal ref: Astrophys.J.707:1217-1232,2009

Spitzer and HHT observations of starless cores: masses and environments

Authors: Amelia M. Stutz, George H. Rieke, John H. Bieging, Zoltan Balog, Fabian Heitsch, Miju Kang, William L. Peters, Yancy L. Shirley, Michael W. Werner

Abstract: 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)… ▽ More 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. △ Less

Submitted 3 November, 2009; v1 submitted 29 October, 2009; originally announced October 2009.

Comments: 41 pages, 28 figures, 5 tables; accepted by ApJ

Journal ref: Astrophys.J.707:137-166,2009

An Observational Determination of the Proton to Electron Mass Ratio in the Early Universe

Authors: Rodger I. Thompson, Jill Bechtold, John H. Black, Daniel Eisenstein, Xiaohui Fan, Robert C. Kennicutt, Carlos Martins, J. Xavier Prochaska, Yancey L. Shirley

Abstract: 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 measur… ▽ More 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. △ Less

Submitted 24 July, 2009; originally announced July 2009.

Comments: Accepted for publication in the Astrophysical Journal

The Spitzer c2d Survey of Nearby Dense Cores: VI. The Protostars of Lynds Dark Nebula 1221

Authors: Chadwick H. Young, Tyler L. Bourke, Michael M. Dunham, Neal J. Evans II, Jes K. Jorgensen, Yancy L. Shirley, Kaisa E. Young, Christopher De Vries, Mark J. Claussen, Victor Popa

Abstract: 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 Col… ▽ More 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. △ Less

Submitted 9 July, 2009; originally announced July 2009.

Comments: accepted for publication in ApJ

Journal ref: Astrophys.J.702:340-351,2009

The Initial Conditions of Clustered Star Formation I: NH3 Observations of Dense Cores in Ophiuchus

Authors: R. K. Friesen, J. Di Francesco, Y. L. Shirley, P. C. Myers

Abstract: 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 micr… ▽ More 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 (Δ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] △ Less

Submitted 3 March, 2009; originally announced March 2009.

Comments: ApJ, accepted

Journal ref: Astrophys.J.697:1457-1480,2009

Spitzer Observations of L429: A Near-collapse or Collapsing Starless Core

Authors: Amelia M. Stutz, Tyler L. Bourke, George H. Rieke, John H. Bieging, Karl A. Misselt, Philip C. Myers, Yancy L. Shirley

Abstract: 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 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. △ Less

Submitted 10 November, 2008; originally announced November 2008.

Comments: 5 pages, 3 figures, emulateapj format, accepted for publication in ApJL

Journal ref: Astrophys.J.Lett.690:L35-L38,2009

Spitzer Observations of Bok Globule B335: Isolated Star Formation Efficiency and Cloud Structure

Authors: Amelia M. Stutz, Mark Rubin, Michael W. Werner, George H. Rieke, John H. Bieging, Jocelyn Keene, Miju Kang, Yancy L. Shirley, K. Y. L. Su, Thangasamy Velusamy, David J. Wilner

Abstract: 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… ▽ More 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. △ Less

Submitted 14 July, 2008; originally announced July 2008.

Comments: 15 pages, 17 figures, emulateapj format, accepted for publication in ApJ

The Star Formation Rate - Dense Gas Relation in Galaxies as Measured by HCN (3-2) Emission

Authors: R. S. Bussmann, D. Narayanan, Y. L. Shirley, S. Juneau, J. Wu, P. M. Solomon, P. A. Vanden Bout, J. Moustakas, C. K. Walker

Abstract: 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 (… ▽ More 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. △ Less

Submitted 28 May, 2008; originally announced May 2008.

Comments: 4 pages, 1 table, 2 figures; accepted to the ApJL

The Properties of Dense Molecular Gas in the Milky Way and Galaxies

Authors: Yancy L. Shirley, Jingwen Wu, R. Shane Bussmann, Al Wootten

Abstract: 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 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. △ Less

Submitted 28 November, 2007; originally announced November 2007.

Comments: 7 pages. To appear in Massive Star Formation: Observations Confront Theory, ed. H. Beuther, ASP Conference Series

Developing A Chemical Evolutionary Sequence for Low-mass Starless Cores

Authors: Yancy L. Shirley

Abstract: 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 $μ$m e… ▽ More 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 $μ$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}. △ Less

Submitted 29 November, 2007; v1 submitted 18 September, 2007; originally announced September 2007.

Comments: To appear in "Molecules in Space & Laboratory", Paris 2007, eds. J.L. Lemaire & F. Combes. Corrected typo, otherwise "more or less" the same

Spitzer observations of a 24 micron shadow: Bok Globule CB190

Authors: Amelia M. Stutz, John H. Bieging, George H. Rieke, Yancy L. Shirley, Zoltan Balog, Karl D. Gordon, Elizabeth M. Green, Jocelyn Keene, Brandon C. Kelly, Mark Rubin, Michael W. Werner

Abstract: 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… ▽ More 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. △ Less

Submitted 15 May, 2007; originally announced May 2007.

Comments: 16 pages, 13 figures, accepted for publication in ApJ

Journal ref: Astrophys.J.665:466-477,2007

The Detection and Characterization of cm Radio Continuum Emission from the Low-mass Protostar L1014-IRS

Authors: Yancy L. Shirley, Mark J. Claussen, Tyler M. Bourke, Chadwick H. Young, Geoffrey A. Blake

Abstract: 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 ca… ▽ More 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 $α= 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σ$ 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. △ Less

Submitted 11 May, 2007; originally announced May 2007.

Comments: 14 pages. Accepted for publication in ApJ

Journal ref: Astrophys.J.667:329-339,2007