Showing 1–36 of 36 results for author: Rosolowsky, E W

PHANGS-HST catalogs for $\sim$100,000 star clusters and compact associations in 38 galaxies: I. Observed properties

Authors: Daniel Maschmann, Janice C. Lee, David A. Thilker, Bradley C. Whitmore, Sinan Deger, Mederic Boquien, Rupali Chandar, Daniel A. Dale, Aida Wofford, Stephen Hannon, Kirsten L. Larson, Adam K. Leroy, Eva Schinnerer, Erik W. Rosolowsky, Leonardo Ubeda, Ashley Barnes, Eric Emsellem, Kathryn Grasha, Brent Groves, Hwihyun Kim, Ralf S. Klessen, Kathryn Kreckel, Rebecca C. Levy, Francesca Pinna, Jimena Rodriguez , et al. (2 additional authors not shown)

Abstract: We present the largest catalog to-date of star clusters and compact associations in nearby galaxies. We have performed a V-band-selected census of clusters across the 38 spiral galaxies of the PHANGS-HST Treasury Survey, and measured integrated, aperture-corrected NUV-U-B-V-I photometry. This work has resulted in uniform catalogs that contain $\sim$20,000 clusters and compact associations which ha… ▽ More We present the largest catalog to-date of star clusters and compact associations in nearby galaxies. We have performed a V-band-selected census of clusters across the 38 spiral galaxies of the PHANGS-HST Treasury Survey, and measured integrated, aperture-corrected NUV-U-B-V-I photometry. This work has resulted in uniform catalogs that contain $\sim$20,000 clusters and compact associations which have passed human inspection and morphological classification, and a larger sample of $\sim$100,000 classified by neural network models. Here, we report on the observed properties of these samples, and demonstrate that tremendous insight can be gained from just the observed properties of clusters, even in the absence of their transformation into physical quantities. In particular, we show the utility of the UBVI color-color diagram, and the three principal features revealed by the PHANGS-HST cluster sample: the young cluster locus, the middle-age plume, and the old globular cluster clump. We present an atlas of maps of the 2D spatial distribution of clusters and compact associations in the context of the molecular clouds from PHANGS-ALMA. We explore new ways of understanding this large dataset in a multi-scale context by bringing together once-separate techniques for the characterization of clusters (color-color diagrams and spatial distributions) and their parent galaxies (galaxy morphology and location relative to the galaxy main sequence). A companion paper presents the physical properties: ages, masses, and dust reddenings derived using improved spectral energy distribution (SED) fitting techniques. △ Less

Submitted 7 March, 2024; originally announced March 2024.

Comments: 48 pages, 26 figures, article in review at APJ

The EDGE-CALIFA Survey: An Extragalactic Database for Galaxy Evolution Studies

Authors: Tony Wong, Yixian Cao, Yufeng Luo, Alberto D. Bolatto, Sebastián F. Sánchez, Jorge K. Barrera-Ballesteros, Leo Blitz, Dario Colombo, Helmut Dannerbauer, Alex Green, Veselina Kalinova, Ferzem Khan, Andrew Kim, Eduardo A. D. Lacerda, Adam K. Leroy, Rebecca C. Levy, Xincheng Lin, Yuanze Luo, Erik W. Rosolowsky, Mónica Rubio, Peter Teuben, Dyas Utomo, Vicente Villanueva, Stuart N. Vogel, Xinyu Wang

Abstract: The EDGE-CALIFA survey provides spatially resolved optical integral field unit (IFU) and CO spectroscopy for 125 galaxies selected from the CALIFA Data Release 3 sample. The Extragalactic Database for Galaxy Evolution (EDGE) presents the spatially resolved products of the survey as pixel tables that reduce the oversampling in the original images and facilitate comparison of pixels from different i… ▽ More The EDGE-CALIFA survey provides spatially resolved optical integral field unit (IFU) and CO spectroscopy for 125 galaxies selected from the CALIFA Data Release 3 sample. The Extragalactic Database for Galaxy Evolution (EDGE) presents the spatially resolved products of the survey as pixel tables that reduce the oversampling in the original images and facilitate comparison of pixels from different images. By joining these pixel tables to lower dimensional tables that provide radial profiles, integrated spectra, or global properties, it is possible to investigate the dependence of local conditions on large-scale properties. The database is freely accessible and has been utilized in several publications. We illustrate the use of this database and highlight the effects of CO upper limits on the inferred slopes of the local scaling relations between stellar mass, star formation rate (SFR), and H$_2$ surface densities. We find that the correlation between H$_2$ and SFR surface density is the tightest among the three relations. △ Less

Submitted 23 January, 2024; originally announced January 2024.

Comments: 21 pages, accepted for publication in ApJS, see DOIs below for code and data access

The EDGE-CALIFA Survey: Spatially Resolved 13CO(1-0) Observations and Variations in 12CO(1-0)/13CO(1-0) in Nearby Galaxies on kpc Scales

Authors: Yixian Cao, Tony Wong, Alberto D. Bolatto, Adam Leroy, Erik W. Rosolowsky, Dyas Utomo, Sebastian Sanchez, Jorge Barrera-Ballesteros, Rebecca Levy, Dario Colombo, Leo Blitz, Stuart Vogel, Johannes Puschnig, Vicente Villanueva, Monica Rubio

Abstract: We present 13CO(1-0) observations for the EDGE-CALIFA survey, which is a mapping survey of 126 nearby galaxies at a typical spatial resolution of 1.5 kpc. Using detected 12CO(1-0) emission as a prior, we detect 13CO(1-0) in 41 galaxies via integrated line flux over the entire galaxy, and in 30 galaxies via integrated line intensity in resolved synthesized beams. Incorporating our CO observations a… ▽ More We present 13CO(1-0) observations for the EDGE-CALIFA survey, which is a mapping survey of 126 nearby galaxies at a typical spatial resolution of 1.5 kpc. Using detected 12CO(1-0) emission as a prior, we detect 13CO(1-0) in 41 galaxies via integrated line flux over the entire galaxy, and in 30 galaxies via integrated line intensity in resolved synthesized beams. Incorporating our CO observations and optical IFU spectroscopy, we perform a systematic comparison between the line ratio R12/13 and the properties of the stars and ionized gas. Higher R12/13 values are found in interacting galaxies than in non-interacting galaxies. The global R12/13 slightly increases with infrared color F60/F100, but appears insensitive to other host galaxy properties such as morphology, stellar mass, or galaxy size. We also present annulus-averaged R12/13 profiles for our sample up to a galactocentric radius of 0.4r25 (~6 kpc), taking into account the 13CO(1-0) non-detections by spectral stacking. The radial profiles of R12/13 are quite flat across our sample. Within galactocentric distances of 0.2r25, azimuthally-averaged R12/13 increases with star formation rate. However, the Spearman rank correlation tests show the azimuthally-averaged R12/13 does not strongly correlate with any other gas or stellar properties in general, especially beyond 0.2r25 from the galaxy centers. Our findings suggest that in the complex environments in galaxy disks, R12/13 is not a sensitive tracer for ISM properties. Dynamical disturbances, like galaxy interactions or the presence of a bar, also have an overall impact on R12/13, which further complicate the interpretations of R12/13 variations. △ Less

Submitted 13 June, 2023; originally announced June 2023.

Comments: 35pages, 11 figure, 6 tables, accepted for publication in ApJS

Fuelling the nuclear ring of NGC 1097

Authors: Mattia C. Sormani, Ashley T. Barnes, Jiayi Sun, Sophia K. Stuber, Eva Schinnerer, Eric Emsellem, Adam K. Leroy, Simon C. O. Glover, Jonathan D. Henshaw, Sharon E. Meidt, Justus Neumann, Miguel Querejeta, Thomas G. Williams, Frank Bigiel, Cosima Eibensteiner, Francesca Fragkoudi, Rebecca C. Levy, Kathryn Grasha, Ralf S. Klessen, J. M. Diederik Kruijssen, Nadine Neumayer, Francesca Pinna, Erik W. Rosolowsky, Rowan J. Smith, Yu-Hsuan Teng , et al. (2 additional authors not shown)

Abstract: Galactic bars can drive cold gas inflows towards the centres of galaxies. The gas transport happens primarily through the so-called bar ``dust lanes'', which connect the galactic disc at kpc scales to the nuclear rings at hundreds of pc scales much like two gigantic galactic rivers. Once in the ring, the gas can fuel star formation activity, galactic outflows, and central supermassive black holes.… ▽ More Galactic bars can drive cold gas inflows towards the centres of galaxies. The gas transport happens primarily through the so-called bar ``dust lanes'', which connect the galactic disc at kpc scales to the nuclear rings at hundreds of pc scales much like two gigantic galactic rivers. Once in the ring, the gas can fuel star formation activity, galactic outflows, and central supermassive black holes. Measuring the mass inflow rates is therefore important to understanding the mass/energy budget and evolution of galactic nuclei. In this work, we use CO datacubes from the PHANGS-ALMA survey and a simple geometrical method to measure the bar-driven mass inflow rate onto the nuclear ring of the barred galaxy NGC~1097. The method assumes that the gas velocity in the bar lanes is parallel to the lanes in the frame co-rotating with the bar, and allows one to derive the inflow rates from sufficiently sensitive and resolved position-position-velocity diagrams if the bar pattern speed and galaxy orientations are known. We find an inflow rate of $\dot{M}=(3.0 \pm 2.1)\, \rm M_\odot\, yr^{-1}$ averaged over a time span of 40 Myr, which varies by a factor of a few over timescales of $\sim$10 Myr. Most of the inflow appears to be consumed by star formation in the ring which is currently occurring at a rate of ${\rm SFR}\simeq~1.8$-$2 \rm M_\odot\, yr^{-1}$, suggesting that the inflow is causally controlling the star formation rate in the ring as a function of time. △ Less

Submitted 23 May, 2023; originally announced May 2023.

Comments: Accepted in MNRAS

Spectral Stacking of Radio-Interferometric Data

Authors: Lukas Neumann, Jakob S. den Brok, Frank Bigiel, Adam Leroy, Antonio Usero, Ashley T. Barnes, Ivana Bešlić, Cosima Eibensteiner, Malena Held, María J. Jiménez-Donaire, Jérôme Pety, Erik W. Rosolowsky, Eva Schinnerer, Thomas G. Williams

Abstract: Mapping molecular line emission beyond the bright low-J CO transitions is still challenging in extragalactic studies, even with the latest generation of (sub-)mm interferometers, such as ALMA and NOEMA. We summarise and test a spectral stacking method that has been used in the literature to recover low-intensity molecular line emission, such as HCN(1-0), HCO+(1-0), and even fainter lines in extern… ▽ More Mapping molecular line emission beyond the bright low-J CO transitions is still challenging in extragalactic studies, even with the latest generation of (sub-)mm interferometers, such as ALMA and NOEMA. We summarise and test a spectral stacking method that has been used in the literature to recover low-intensity molecular line emission, such as HCN(1-0), HCO+(1-0), and even fainter lines in external galaxies. The goal is to study the capabilities and limitations of the stacking technique when applied to imaged interferometric observations. The core idea of spectral stacking is to align spectra of the low S/N spectral lines to a known velocity field calculated from a higher S/N line expected to share the kinematics of the fainter line, e.g., CO(1-0) or 21-cm emission. Then these aligned spectra can be coherently averaged to produce potentially high S/N spectral stacks. Here, we use imaged simulated interferometric and total power observations at different signal-to-noise levels, based on real CO observations. For the combined interferometric and total power data, we find that the spectral stacking technique is capable of recovering the integrated intensities even at low S/N levels across most of the region where the high S/N prior is detected. However, when stacking interferometer-only data for low S/N emission, the stacks can miss up to 50% of the emission from the fainter line. A key result of this analysis is that the spectral stacking method is able to recover the true mean line intensities in low S/N cubes and to accurately measure the statistical significance of the recovered lines. To facilitate the application of this technique we provide a public Python package, called PyStacker. △ Less

Submitted 17 May, 2023; originally announced May 2023.

Comments: 10 pages, 10 figures, accepted for pub in A&A, Apr 28, 2023

Journal ref: A&A 675, A104 (2023)

Wide-field CO isotopologue emission and the CO-to-H$_2$ factor across the nearby spiral galaxy M101

Authors: Jakob S. den Brok, Frank Bigiel, Jérémy Chastenet, Karin Sandstrom, Adam Leroy, Antonio Usero, Eva Schinnerer, Erik W. Rosolowsky, Eric W. Koch, I-Da Chiang, Ashley T. Barnes, Johannes Puschnig, Toshiki Saito, Ivana Bešlić, Melanie Chevance, Daniel A. Dale, Cosima Eibensteiner, Simon Glover, María J. Jiménez-Donaire, Yu-Hsuan Teng, Thomas G. Williams

Abstract: Carbon monoxide (CO) emission is the most widely used tracer of the bulk molecular gas in the interstellar medium (ISM) in extragalactic studies. The CO-to-H$_2$ conversion factor, $α_{\rm CO}$, links the observed CO emission to the total molecular gas mass. However, no single prescription perfectly describes the variation of $α_{\rm CO}$ across all environments across galaxies as a function of me… ▽ More Carbon monoxide (CO) emission is the most widely used tracer of the bulk molecular gas in the interstellar medium (ISM) in extragalactic studies. The CO-to-H$_2$ conversion factor, $α_{\rm CO}$, links the observed CO emission to the total molecular gas mass. However, no single prescription perfectly describes the variation of $α_{\rm CO}$ across all environments across galaxies as a function of metallicity, molecular gas opacity, line excitation, and other factors. Using resolved spectral line observations of CO and its isotopologues, we can constrain the molecular gas conditions and link them to a variation in the conversion factor. We present new IRAM 30-m 1mm and 3mm line observations of $^{12}$CO, $^{13}$CO, and C$^{18}$O} across the nearby galaxy M101. Based on the CO isotopologue line ratios, we find that selective nucleosynthesis and opacity changes are the main drivers of the variation in the line emission across the galaxy. Furthermore, we estimated $α_{\rm CO(1-0)}$ using different approaches, including (i) the dust mass surface density derived from far-IR emission as an independent tracer of the total gas surface density and (ii) LTE-based measurements using the optically thin $^{13}$CO(1-0) intensity. We find an average value of $α_{\rm CO}=4.4{\pm}0.9\rm\,M_\odot\,pc^{-2}(K\,km\,s^{-1})^{-1}$ across the galaxy, with a decrease by a factor of 10 toward the 2 kpc central region. In contrast, we find LTE-based values are lower by a factor of 2-3 across the disk relative to the dust-based result. Accounting for $α_{\rm CO}$ variations, we found significantly reduced molecular gas depletion time by a factor 10 in the galaxy's center. In conclusion, our result suggests implications for commonly derived scaling relations, such as an underestimation of the slope of the Kennicutt Schmidt law, if $α_{\rm CO}$ variations are not accounted for. △ Less

Submitted 6 February, 2023; originally announced February 2023.

Comments: Accepted for publication in A&A, 25 pages, 15 figures

Journal ref: A&A 676, A93 (2023)

PHANGS-JWST First Results: Multi-wavelength view of feedback-driven bubbles (The Phantom Voids) across NGC 628

Authors: Ashley T. Barnes, Elizabeth J. Watkins, Sharon E. Meidt, Kathryn Kreckel, Mattia C. Sormani, Robin G. Tress, Simon C. O. Glover, Frank Bigiel, Rupali Chandar, Eric Emsellem, Janice C. Lee, Adam K. Leroy, Karin M. Sandstrom, Eva Schinnerer, Erik W. Rosolowsky, Francesco Belfiore, Guillermo Blanc, Mederic Boquien, Jakob S. den Brok, Yixian Cao, Mélanie Chevance, Daniel A. Dale, Oleg Egorov, Cosima Eibensteiner, Kathryn Grasha , et al. (29 additional authors not shown)

Abstract: We present a high-resolution view of bubbles within The Phantom Galaxy (NGC 628); a nearby (~10Mpc), star-forming (~2Msun/yr), face-on (i~9deg) grand-design spiral galaxy. With new data obtained as part of the PHANGS-JWST treasury program, we perform a detailed case-study of two regions of interest, one of which contains the largest and most prominent bubble in the galaxy (The Phantom Void; over 1… ▽ More We present a high-resolution view of bubbles within The Phantom Galaxy (NGC 628); a nearby (~10Mpc), star-forming (~2Msun/yr), face-on (i~9deg) grand-design spiral galaxy. With new data obtained as part of the PHANGS-JWST treasury program, we perform a detailed case-study of two regions of interest, one of which contains the largest and most prominent bubble in the galaxy (The Phantom Void; over 1kpc in diameter), and the other being a smaller region that may be the precursor to such a large bubble (The Precursor Phantom Void). When comparing to matched resolution Halpha observations from the Hubble Space Telescope (HST), we see that the ionized gas is brightest in the shells of both bubbles, and is coincident with the youngest (~1Myr) and most massive (~100,000Msun) stellar associations. We also find an older generation (~20Myr) of stellar associations is present within the bubble of The Phantom Void. From our kinematic analysis of the HI, H2 (CO) and HII gas across The Phantom Void, we infer a high expansion speed of around 15 to 50km/s. The large size and high expansion speed of The Phantom Void suggest that the driving mechanism is sustained stellar feedback due to multiple mechanisms, where early feedback first cleared a bubble (as we observe now in The Precursor Phantom Void), and since then SNe have been exploding within the cavity, and have accelerated the shell. Finally, comparison to simulations shows a striking resemblance to our JWST observations, and suggests that such large-scale stellar feedback-driven bubbles should be common within other galaxies. △ Less

Submitted 1 December, 2022; originally announced December 2022.

Comments: 12 pages total, 8 figures, and 1 table. Accepted as part of a PHANGS-JWST Focus Issue to appear in ApJ

PHANGS-JWST First Results: A statistical view on bubble evolution in NGC628

Authors: Elizabeth J. Watkins, Ashley Barnes, Kiana F. Henny, Hwihyun Kim, Kathryn Kreckel, Sharon E. Meidt, Ralf S. Klessen, Simon C. O. Glover, Thomas G. Williams, B. W. Keller, Adam K. Leroy, Erik W. Rosolowsky, Mederic Boquien, Gagandeep S. Anand, Francesco Belfiore, Frank Bigiel, Guillermo Blanc, Yixian Cao, Rupali Chandar, Ness Mayker Chen, Mélanie Chevance, Enrico Congiu, Daniel A. Dale, Sinan Deger, Oleg Egorov , et al. (27 additional authors not shown)

Abstract: The first JWST observations of nearby galaxies have unveiled a rich population of bubbles that trace the stellar feedback mechanisms responsible for their creation. Studying these bubbles therefore allows us to chart the interaction between stellar feedback and the interstellar medium, and the larger galactic flows needed to regulate star formation processes globally. We present the first catalog… ▽ More The first JWST observations of nearby galaxies have unveiled a rich population of bubbles that trace the stellar feedback mechanisms responsible for their creation. Studying these bubbles therefore allows us to chart the interaction between stellar feedback and the interstellar medium, and the larger galactic flows needed to regulate star formation processes globally. We present the first catalog of bubbles in NGC628, visually identified using MIRI F770W PHANGS-JWST observations, and use them to statistically evaluate bubble characteristics. We classify 1694 structures as bubbles with radii between 6-552 pc. Of these, 31% contain at least one smaller bubble at their edge, indicating that previous generations of star formation have a local impact on where new stars form. On large scales, most bubbles lie near a spiral arm, and their radii increase downstream compared to upstream. Furthermore, bubbles are elongated in a similar direction to the spiral arm ridge-line. These azimuthal trends demonstrate that star formation is intimately connected to the spiral arm passage. Finally, the bubble size distribution follows a power-law of index $p=-2.2\pm0.1$, which is slightly shallower than the theoretical value by 1-3.5$σ$ that did not include bubble mergers. The fraction of bubbles identified within the shells of larger bubbles suggests that bubble merging is a common process. Our analysis therefore allows us to quantify the number of star-forming regions that are influenced by an earlier generation, and the role feedback processes have in setting the global star formation rate. With the full PHANGS-JWST sample, we can do this for more galaxies. △ Less

Submitted 1 December, 2022; originally announced December 2022.

Comments: 25 pages total, 13 Figures and 1 Table. Accepted for publication in ApJL as part of a PHANGS-JWST First Results Focus issue

PHANGS-JWST First Results: Dust embedded star clusters in NGC 7496 selected via 3.3 $μ$m PAH emission

Authors: Jimena Rodriguez, Janice Lee, Bradley Whitmore, David Thilker, Daniel Maschmann, Rupali Chandar, Daniel Dale, Diederik Kruijssen, Mederic Boquien, Kathryn Grasha, Elizabeth Watkins, Ashley Barnes, Mattia Sormani, Thomas Williams, Jaeyeon Kim, Gagandeep Anand, Mélanie Chevance, Frank Bigiel, Adam Leroy, Ralf Klessen, Erik W. Rosolowsky, Karin Sandstrom, Hamid Hassani, Hwihyun Kim, Eva Schinnerer , et al. (16 additional authors not shown)

Abstract: The earliest stages of star formation occur enshrouded in dust and are not observable in the optical. Here we leverage the extraordinary new high-resolution infrared imaging from JWST to begin the study of dust-embedded star clusters in nearby galaxies throughout the local volume. We present a technique for identifying dust-embedded clusters in NGC 7496 (18.7 Mpc), the first galaxy to be observed… ▽ More The earliest stages of star formation occur enshrouded in dust and are not observable in the optical. Here we leverage the extraordinary new high-resolution infrared imaging from JWST to begin the study of dust-embedded star clusters in nearby galaxies throughout the local volume. We present a technique for identifying dust-embedded clusters in NGC 7496 (18.7 Mpc), the first galaxy to be observed by the PHANGS-JWST Cycle 1 Treasury Survey. We select sources that have strong 3.3$μ$m PAH emission based on a $\rm F300M-F335M$ color excess, and identify 67 candidate embedded clusters. Only eight of these are found in the PHANGS-HST optically-selected cluster catalog and all are young (six have SED-fit ages of $\sim1$ Myr). We find that this sample of embedded cluster candidates may significantly increase the census of young clusters in NGC 7496 from the PHANGS-HST catalog -- the number of clusters younger than $\sim$2 Myr could be increased by a factor of two. Candidates are preferentially located in dust lanes, and are coincident with peaks in PHANGS-ALMA CO (2-1) maps. We take a first look at concentration indices, luminosity functions, SEDs spanning from 2700A to 21$μ$m, and stellar masses (estimated to be between $\sim10^4-10^5 M_{\odot}$). The methods tested here provide a basis for future work to derive accurate constraints on the physical properties of embedded clusters, characterize the completeness of cluster samples, and expand analysis to all 19 galaxies in the PHANGS-JWST sample, which will enable basic unsolved problems in star formation and cluster evolution to be addressed. △ Less

Submitted 28 November, 2022; v1 submitted 24 November, 2022; originally announced November 2022.

Comments: 12 pages, 6 figures; accepted for publication in ApJL as part of PHANGS-JWST First Results Special Issue

The Life and Times of Giant Molecular Clouds

Authors: Mélanie Chevance, Mark R. Krumholz, Anna F. McLeod, Eve C. Ostriker, Erik W. Rosolowsky, Amiel Sternberg

Abstract: Giant molecular clouds (GMCs) are the sites of star formation and stellar feedback in galaxies. Their properties set the initial conditions for star formation and their lifecycles determine how feedback regulates galaxy evolution. In recent years, the advent of high-resolution telescopes has enabled systematic GMC-scale studies of the molecular interstellar medium in nearby galaxies, now covering… ▽ More Giant molecular clouds (GMCs) are the sites of star formation and stellar feedback in galaxies. Their properties set the initial conditions for star formation and their lifecycles determine how feedback regulates galaxy evolution. In recent years, the advent of high-resolution telescopes has enabled systematic GMC-scale studies of the molecular interstellar medium in nearby galaxies, now covering a wide range of physical conditions and allowing for the first studies of how GMC properties depend on galactic environment. These observational developments have been accompanied by numerical simulations of improving resolution that are increasingly accurately accounting for the effects of the galactic-scale environment on GMCs, while simultaneously improving the treatment of the small-scale processes of star-formation and stellar feedback within them. The combination of these recent developments has greatly improved our understanding of the formation, evolution, and destruction of GMCs. We review the current state of the field, highlight current open questions, and discuss promising avenues for future studies. △ Less

Submitted 17 March, 2022; originally announced March 2022.

Comments: 37 pages, 5 figures; To appear in Protostars and Planets VII; Editors: Shu-ichiro Inutsuka, Yuri Aikawa, Takayuki Muto, Kengo Tomida, and Motohide Tamura

Are massive dense clumps truly sub-virial? A new analysis using Gould Belt ammonia data

Authors: Ayushi Singh, Christopher D. Matzner, Rachel K. Friesen, Peter G. Martin, Jaime E. Pineda, Erik W. Rosolowsky, Felipe Alves, Ana Chacón-Tanarro, Hope How-Huan Chen, Michael Chun-Yuan Chen, Spandan Choudhury, James Di Francesco, Jared Keown, Helen Kirk, Anna Punanova, Youngmin Seo, Yancy Shirley, Adam Ginsburg, Stella S. R. Offner, Héctor G. Arce, Paola Caselli, Alyssa A. Goodman, Philip C. Myers, Elena Redaelli

Abstract: Dynamical studies of dense structures within molecular clouds often conclude that the most massive clumps contain too little kinetic energy for virial equilibrium, unless they are magnetized to an unexpected degree. This raises questions about how such a state might arise, and how it might persist long enough to represent the population of massive clumps. In an effort to re-examine the origins of… ▽ More Dynamical studies of dense structures within molecular clouds often conclude that the most massive clumps contain too little kinetic energy for virial equilibrium, unless they are magnetized to an unexpected degree. This raises questions about how such a state might arise, and how it might persist long enough to represent the population of massive clumps. In an effort to re-examine the origins of this conclusion, we use ammonia line data from the Green Bank Ammonia Survey and Planck-calibrated dust emission data from Herschel to estimate the masses and kinetic and gravitational energies for dense clumps in the Gould Belt clouds. We show that several types of systematic error can enhance the appearance of low kinetic-to-gravitational energy ratios: insufficient removal of foreground and background material; ignoring the kinetic energy associated with velocity differences across a resolved cloud; and over-correcting for stratification when evaluating the gravitational energy. Using an analysis designed to avoid these errors, we find that the most massive Gould Belt clumps harbor virial motions, rather than sub-virial ones. As a byproduct, we present a catalog of masses, energies, and virial energy ratios for 85 Gould Belt clumps. △ Less

Submitted 11 August, 2021; originally announced August 2021.

Comments: Submitted to ApJ

A lack of constraints on the cold opaque HI mass: HI spectra in M31 and M33 prefer multi-component models over a single cold opaque component

Authors: Eric W. Koch, Erik W. Rosolowsky, Adam K. Leroy, Jeremy Chastenet, I-Da Chiang, Julianne Dalcanton, Amanda A. Kepley, Karin M. Sandstrom, Andreas Schruba, Snezana Stanimirovic, Dyas Utomo, Thomas G. Williams

Abstract: Previous work has argued that atomic gas mass estimates of galaxies from 21 cm HI emission are systematically low due to a cold opaque atomic gas component. If true, this opaque component necessitates a ~35% correction factor relative to the mass from assuming optically-thin HI emission. These mass corrections are based on fitting HI spectra with a single opaque component model that produces a dis… ▽ More Previous work has argued that atomic gas mass estimates of galaxies from 21 cm HI emission are systematically low due to a cold opaque atomic gas component. If true, this opaque component necessitates a ~35% correction factor relative to the mass from assuming optically-thin HI emission. These mass corrections are based on fitting HI spectra with a single opaque component model that produces a distinct "top-hat" shaped line profile. Here, we investigate this issue using deep, high spectral resolution HI VLA observations of M31 and M33 to test if these top-hat profiles are instead superpositions of multiple HI components along the line-of-sight. We fit both models and find that >80% of the spectra strongly prefer a multi-component Gaussian model while <2% prefer the single opacity-corrected component model. This strong preference for multiple components argues against previous findings of lines-of-sight dominated by only cold HI. Our findings are enabled by the improved spectral resolution (0.42 km/s), whereas coarser spectral resolution blends multiple components together. We also show that the inferred opaque atomic ISM mass strongly depends on the goodness-of-fit definition and is highly uncertain when the inferred spin temperature has a large uncertainty. Finally, we find that the relation of the HI surface density with the dust surface density and extinction has significantly more scatter when the inferred HI opacity correction is applied. These variations are difficult to explain without additionally requiring large variations in the dust properties. Based on these findings, we suggest that the opaque HI mass is best constrained by HI absorption studies. △ Less

Submitted 6 April, 2021; originally announced April 2021.

Comments: Accepted in MNRAS. 25 pg, 13 figures

Jet-ISM Interactions near the Microquasars GRS 1758-258 and 1E 1740.7-2942

Authors: A. J. Tetarenko, E. W. Rosolowsky, J. C. A Miller-Jones, G. R. Sivakoff

Abstract: We present Atacama Large Millimeter/Sub-millimeter Array observations of the candidate jet-ISM interaction zones near the black hole X-ray binaries GRS 1758$-$258 and 1E 1740.7$-$2942. Using these data, we map the molecular line emission in the regions, detecting emission from the HCN [$J=1-0$], HCO$^+$ [$J=1-0$], SiO [$J=2-1$], CS [$J=2-1$], $^{13}$CO [$J=1-0$], C$^{18}$O [$J=1-0$], HNCO [… ▽ More We present Atacama Large Millimeter/Sub-millimeter Array observations of the candidate jet-ISM interaction zones near the black hole X-ray binaries GRS 1758$-$258 and 1E 1740.7$-$2942. Using these data, we map the molecular line emission in the regions, detecting emission from the HCN [$J=1-0$], HCO$^+$ [$J=1-0$], SiO [$J=2-1$], CS [$J=2-1$], $^{13}$CO [$J=1-0$], C$^{18}$O [$J=1-0$], HNCO [$J=4_{0,4}-3_{0,3}$], HNCO [$J=5_{0,5}-4_{0,4}$], and CH$_3$OH [$J=2_{1,1}-1_{1,0}$] molecular transitions. Through examining the morphological, spectral, and kinematic properties of this emission, we identify molecular structures that may trace jet-driven cavities in the gas surrounding these systems. Our results from the GRS 1758$-$258 region in particular, are consistent with recent work, which postulated the presence of a jet-blown cocoon structure in deep radio continuum maps of the region. Using these newly discovered molecular structures as calorimeters, we estimate the time averaged jet power from these systems, finding $(1.1-5.7)\times10^{36}{\rm erg\,s}^{-1}$ over $0.12-0.31$ Myr for GRS 1758$-$258 and $(0.7-3.5)\times10^{37}{\rm erg\,s}^{-1}$ over $0.10-0.26$ Myr for 1E 1740.7$-$2942. Additionally, the spectral line characteristics of the detected emission place these molecular structures in the central molecular zone of our Galaxy, thereby constraining the distances to the black hole X-ray binaries to be $8.0\pm1.0$ kpc. Overall, our analysis solidifies the diagnostic capacity of molecular lines, and highlights how astro-chemistry can both identify jet-ISM interaction zones and probe jet feedback from Galactic X-ray binaries. △ Less

Submitted 21 July, 2020; originally announced July 2020.

Comments: 24 pages, 21 figures, accepted to MNRAS

The X-ray Emissivity of Low-Density Stellar Populations

Authors: C. O. Heinke, M. G. Ivanov, E. W. Koch, R. Andrews, L. Chomiuk, H. N. Cohn, S. Crothers, T. de Boer, N. Ivanova, A. K. H. Kong, N. W. C. Leigh, P. M. Lugger, L. Nelson, C. J. Parr, E. W. Rosolowsky, A. J. Ruiter, C. L. Sarazin, A. W. Shaw, G. R. Sivakoff, M. van den Berg

Abstract: The dynamical production of low-mass X-ray binaries and brighter cataclysmic variables (CVs) in dense globular clusters is well-established. We investigate how the X-ray emissivity of fainter X-ray binaries (principally CVs and coronally active binaries) varies between different environments. We compile calculations (largely from the literature) of the X-ray emissivity of old stellar populations,… ▽ More The dynamical production of low-mass X-ray binaries and brighter cataclysmic variables (CVs) in dense globular clusters is well-established. We investigate how the X-ray emissivity of fainter X-ray binaries (principally CVs and coronally active binaries) varies between different environments. We compile calculations (largely from the literature) of the X-ray emissivity of old stellar populations, including open and globular clusters and several galaxies. We investigate three literature claims of unusual X-ray sources in low-density stellar populations. We show that a suggested quiescent neutron star in the open cluster NGC 6819 is a foreground M dwarf. We show that the suggested diffuse X-ray emission from an old nova shell in the globular cluster NGC 6366 is actually a background galaxy cluster. And we show that a suggested population of quiescent X-ray binaries in the Sculptor Dwarf Galaxy is mostly (perhaps entirely) background galaxies. We find that above densities of $10^4$ M$_{\odot}$/pc$^3$, the X-ray emissivity of globular clusters increases, due to dynamical production of X-ray emitting systems. Below this density, globular clusters have lower X-ray emissivity than the other populations, and we do not see a strong dependence of X-ray emissivity due to density effects. We find significant correlations between X-ray emissivity and binary fraction, metallicity, and density. Sampling these fits via bootstrap techniques gives less significant correlations, but confirms the effect of metallicity on low-density populations, and that of density on the full globular cluster sample. △ Less

Submitted 26 January, 2020; originally announced January 2020.

Comments: 28 pages, 23 figures. MNRAS, in press

Spatial Power Spectra of Dust across the Local Group: No Constraint on Disc Scale Height

Authors: Eric W. Koch, I-Da Chiang, Dyas Utomo, Jérémy Chastenet, Adam K. Leroy, Erik W. Rosolowsky, Karin M. Sandstrom

Abstract: We analyze the 1D spatial power spectra of dust surface density and mid to far-infrared emission at $24-500\,μ$m in the LMC, SMC, M31, and M33. By forward-modelling the point-spread-function (PSF) on the power spectrum, we find that nearly all power spectra have a single power-law and point source component. A broken power-law model is only favoured for the LMC 24 $μ$m MIPS power spectrum and is d… ▽ More We analyze the 1D spatial power spectra of dust surface density and mid to far-infrared emission at $24-500\,μ$m in the LMC, SMC, M31, and M33. By forward-modelling the point-spread-function (PSF) on the power spectrum, we find that nearly all power spectra have a single power-law and point source component. A broken power-law model is only favoured for the LMC 24 $μ$m MIPS power spectrum and is due to intense dust heating in 30 Doradus. We also test for local power spectrum variations by splitting the LMC and SMC maps into $820$ pc boxes. We find significant variations in the power-law index with no strong evidence for breaks. The lack of a ubiquitous break suggests that the spatial power spectrum does not constrain the disc scale height. This contradicts claims of a break where the turbulent motion changes from 3D to 2D. The power spectrum indices in the LMC, SMC, and M31 are similar ($2.0-2.5$). M33 has a flatter power spectrum ($1.3$), similar to more distant spiral galaxies with a centrally-concentrated H$_{2}$ distribution. We compare the power spectra of HI, CO, and dust in M31 and M33 and find that HI power spectra are consistently flatter than CO power spectra. These results cast doubt on the idea that the spatial power spectrum traces large scale turbulent motion in nearby galaxies. Instead, we find that the spatial power spectrum is influenced by (1) the PSF on scales below $\sim3$ times the FWHM, (2) bright compact regions (30 Doradus), and (3) the global morphology of the tracer (an exponential CO disc). △ Less

Submitted 10 January, 2020; v1 submitted 19 December, 2019; originally announced December 2019.

Comments: Accepted by MNRAS; 19 pgs, 10 figures. Supplementary figures of the power spectrum fits for each band can be viewed here: https://doi.org/10.7910/DVN/RYRXK7

TurbuStat: Turbulence Statistics in Python

Authors: Eric W. Koch, Erik W. Rosolowsky, Ryan D. Boyden, Blakesley Burkhart, Adam Ginsburg, Jason L. Loeppky, Stella S. R. Offner

Abstract: We present TurbuStat (v1.0): a Python package for computing turbulence statistics in spectral-line data cubes. TurbuStat includes implementations of fourteen methods for recovering turbulent properties from observational data. Additional features of the software include: distance metrics for comparing two data sets; a segmented linear model for fitting lines with a break-point; a two-dimensional e… ▽ More We present TurbuStat (v1.0): a Python package for computing turbulence statistics in spectral-line data cubes. TurbuStat includes implementations of fourteen methods for recovering turbulent properties from observational data. Additional features of the software include: distance metrics for comparing two data sets; a segmented linear model for fitting lines with a break-point; a two-dimensional elliptical power-law model; multi-core fast-fourier-transform support; a suite for producing simulated observations of fractional Brownian Motion fields, including two-dimensional images and optically-thin HI data cubes; and functions for creating realistic world coordinate system information for synthetic observations. This paper summarizes the TurbuStat package and provides representative examples using several different methods. TurbuStat is an open-source package and we welcome community feedback and contributions. △ Less

Submitted 25 April, 2019; v1 submitted 23 April, 2019; originally announced April 2019.

Comments: Accepted in AJ. 21 pages, 8 figures

Relationship between the Line Width of the Atomic and Molecular ISM in M33

Authors: Eric W. Koch, Erik W. Rosolowsky, Andreas Schruba, Adam Leroy, Amanda Kepley, Jonathan Braine, Julianne Dalcanton, Megan C. Johnson

Abstract: We investigate how the spectral properties of atomic (HI) and molecular (H$_2$) gas, traced by CO(2-1), are related in M33 on $80$ pc scales. We find the HI and CO(2-1) velocity at peak intensity to be highly correlated, consistent with previous studies. By stacking spectra aligned to the velocity of HI peak intensity, we find that the CO line width ($σ_{\rm HWHM}=4.6\pm0.9$ km s$^{-1}$;… ▽ More We investigate how the spectral properties of atomic (HI) and molecular (H$_2$) gas, traced by CO(2-1), are related in M33 on $80$ pc scales. We find the HI and CO(2-1) velocity at peak intensity to be highly correlated, consistent with previous studies. By stacking spectra aligned to the velocity of HI peak intensity, we find that the CO line width ($σ_{\rm HWHM}=4.6\pm0.9$ km s$^{-1}$; $σ_{\rm HWHM}$ is the effective Gaussian width) is consistently smaller than the HI line width ($σ_{\rm HWHM}=6.6\pm0.1$ km s$^{-1}$), with a ratio of ${\sim}0.7$, in agreement with Druard et al. (2014). The ratio of the line widths remains less than unity when the data are smoothed to a coarser spatial resolution. In other nearby galaxies, this line width ratio is close to unity which has been used as evidence for a thick, diffuse molecular disk that is distinct from the thin molecular disk dominated by molecular clouds. The smaller line width ratio found here suggests that M33 has a marginal thick molecular disk. From modelling individual lines-of-sight, we recover a strong correlation between HI and CO line widths when only the HI located closest to the CO component is considered. The median line width ratio of the line-of-sight line widths is $0.56\pm0.01$. There is substantial scatter in the HI--CO(2-1) line width relation, larger than the uncertainties, that results from regional variations on $<500$ pc scales, and there is no significant trend in the line widths, or their ratios, with galactocentric radius. These regional line width variations may be a useful probe of changes in the local cloud environment or the evolutionary state of molecular clouds. △ Less

Submitted 13 February, 2019; originally announced February 2019.

Comments: Accepted for publication in MNRAS; 20 pages, 15 figures

Kinematics of the Atomic ISM in M33 on 80 pc scales

Authors: Eric W. Koch, Erik W. Rosolowsky, Felix J. Lockman, Amanda A. Kepley, Adam Leroy, Andreas Schruba, Jonathan Braine, Julianne Dalcanton, Megan C. Johnson, Snežana Stanimirović

Abstract: We present new L-band ($1\mbox{--}2$ GHz) observations of the nearby spiral galaxy M33 with 80 pc resolution obtained with the Karl G. Jansky Very Large Array. The HI observations, combined with HI measurements from the Green Bank Telescope, improve the spectral resolution and sensitivity ($2.8$ K rms noise in a $0.2$ km s$^{-1}$ channel) compared to previous observations. We find individual profi… ▽ More We present new L-band ($1\mbox{--}2$ GHz) observations of the nearby spiral galaxy M33 with 80 pc resolution obtained with the Karl G. Jansky Very Large Array. The HI observations, combined with HI measurements from the Green Bank Telescope, improve the spectral resolution and sensitivity ($2.8$ K rms noise in a $0.2$ km s$^{-1}$ channel) compared to previous observations. We find individual profiles are usually non-Gaussian, harbouring line wings, multiple components, and asymmetries. Given this spectral complexity, we quantify the motions in the atomic ISM through moment analysis of the spectra and fits to aligned, stacked profiles. The measured value of the HI line width depends strongly on the method used, with the velocity stacked profiles aligned to the peak velocity giving the minimum value of $σ= 7$ km s$^{-1}$ and all other methods giving higher values ($σ\sim10$ km s$^{-1}$). All measurements of the line width show a shallow radial trend, with $σ$ decreasing by $\sim2$ km s$^{-1}$ from $R_{\rm gal}=0$ to $R_{\rm gal}=8$ kpc. We consider a number of energy sources that might maintain the line width against turbulent dissipation, but no single source is adequate. We find excess emission relative to a Gaussian in the stacked profile line wings, ranging from 9% to 26% depending on how the spectra are aligned. By splitting the line wings into symmetric and asymmetric components, we find that the lagging rotational disk accounts for one-third of the line wing flux. We also find emission far from the rotation axis of the galaxy in multiple discrete HI clouds, including a filament with a projected length of $\sim8$ kpc. △ Less

Submitted 8 August, 2018; v1 submitted 28 June, 2018; originally announced June 2018.

Comments: MNRAS, in press. 30 pages, 24 figures

Assessing the Impact of Astrochemistry on Molecular Cloud Turbulence Statistics

Authors: Ryan D. Boyden, Stella S. R. Offner, Eric W. Koch, Erik W. Rosolowsky

Abstract: We analyze hydrodynamic simulations of turbulent, star-forming molecular clouds that are post-processed with the photo-dissociation region astrochemistry code 3D-PDR. We investigate the sensitivity of 15 commonly applied turbulence statistics to post-processing assumptions, namely variations in gas temperature, abundance and external radiation field. We produce synthetic $^{12}$CO(1-0) and CI(… ▽ More We analyze hydrodynamic simulations of turbulent, star-forming molecular clouds that are post-processed with the photo-dissociation region astrochemistry code 3D-PDR. We investigate the sensitivity of 15 commonly applied turbulence statistics to post-processing assumptions, namely variations in gas temperature, abundance and external radiation field. We produce synthetic $^{12}$CO(1-0) and CI($^{3}$P$_{1}$-$^{3}$P$_{0}$) observations and examine how the variations influence the resulting emission distributions. To characterize differences between the datasets, we perform statistical measurements, identify diagnostics sensitive to our chemistry parameters, and quantify the statistic responses by using a variety of distance metrics. We find that multiple turbulent statistics are sensitive not only to the chemical complexity but also to the strength of the background radiation field. The statistics with meaningful responses include principal component analysis, spatial power spectrum and bicoherence. A few of the statistics, such as the velocity coordinate spectrum, are primarily sensitive to the type of tracer being utilized, while others, like the delta-variance, strongly respond to the background radiation field. Collectively, these findings indicate that more realistic chemistry impacts the responses of turbulent statistics and is necessary for accurate statistical comparisons between models and observed molecular clouds. △ Less

Submitted 24 May, 2018; originally announced May 2018.

Comments: 27 pages, 21 figures, accepted to ApJ

Mapping Jet-ISM Interactions in X-ray Binaries with ALMA: A GRS 1915$+$105 Case Study

Authors: A. J. Tetarenko, P. Freeman, E. W. Rosolowsky, J. C. A Miller-Jones, G. R. Sivakoff

Abstract: We present Atacama Large Millimetre/Sub-Millimetre Array (ALMA) observations of IRAS 19132+1035, a candidate jet-ISM interaction zone near the black hole X-ray binary (BHXB) GRS 1915+105. With these ALMA observations (combining data from the 12 m array and the Atacama Compact Array), we map the molecular line emission across the IRAS 19132+1035 region. We detect emission from the $^{12}$CO [… ▽ More We present Atacama Large Millimetre/Sub-Millimetre Array (ALMA) observations of IRAS 19132+1035, a candidate jet-ISM interaction zone near the black hole X-ray binary (BHXB) GRS 1915+105. With these ALMA observations (combining data from the 12 m array and the Atacama Compact Array), we map the molecular line emission across the IRAS 19132+1035 region. We detect emission from the $^{12}$CO [$J=2-1$], $^{13}$CO [$ν=0$, $J=2-1$], C$^{18}$O [$J=2-1$], ${\rm H}_{2}{\rm CO}$ [$J=3_{0,3}-2_{0,2}$], ${\rm H}_{2}{\rm CO}$ [$J=3_{2,2}-2_{2,1}$], ${\rm H}_{2}{\rm CO}$ [$J=3_{2,1}-2_{2,0}$], SiO [$ν=0$, $J=5-4$], CH$_3$OH [$J=4_{2,2}-3_{1,2}$], and CS [$ν=0$, $J=5-4$] transitions. Given the morphological, spectral, and kinematic properties of this molecular emission, we present several lines of evidence that support the presence of a jet-ISM interaction at this site, including a jet-blown cavity in the molecular gas. This compelling new evidence identifies this site as a jet-ISM interaction zone, making GRS 1915$+$105 the third Galactic BHXB with at least one conclusive jet-ISM interaction zone. However, we find that this interaction occurs on much smaller scales than was postulated by previous work, where the BHXB jet does not appear to be dominantly powering the entire IRAS 19132+1035 region. Using estimates of the ISM conditions in the region, we utilize the detected cavity as a calorimeter to estimate the time-averaged power carried in the GRS 1915+105 jets of $(8.4^{+7.7}_{-8.1})\times10^{32}\,{\rm erg\,s}^{-1}$. Overall, our analysis demonstrates that molecular lines are excellent diagnostic tools to identify and probe jet-ISM interaction zones near Galactic BHXBs. △ Less

Submitted 1 December, 2017; originally announced December 2017.

Comments: 23 pages, 16 Figures, Accepted to MNRAS

Identifying Tools for Comparing Simulations and Observations of Spectral-line Data Cubes

Authors: Eric W. Koch, Caleb G. Ward, Stella Offner, Jason L. Loeppky, Erik W. Rosolowsky

Abstract: We present a statistical framework to compare spectral-line data cubes of molecular clouds and use the framework to perform an analysis of various statistical tools developed from methods proposed in the literature. We test whether our methods are sensitive to changes in the underlying physical properties of the clouds or whether their behaviour is governed by random fluctuations. We perform a set… ▽ More We present a statistical framework to compare spectral-line data cubes of molecular clouds and use the framework to perform an analysis of various statistical tools developed from methods proposed in the literature. We test whether our methods are sensitive to changes in the underlying physical properties of the clouds or whether their behaviour is governed by random fluctuations. We perform a set of 32 self-gravitating magnetohydrodynamic simulations that test all combinations of five physical parameters -- Mach number, plasma parameter, virial parameter, driving scales, and solenoidal driving fraction -- each of which can be set to a low or high value. We create mock observational data sets of ${\rm ^{13}CO}$(1-0) emission from each simulation. We compare these mock data to a those generated from a set of baseline simulations using pseudo-distance metrics based on 18 different statistical techniques that have previously been used to study molecular clouds. We analyze these results using methods from the statistical field of experimental design and find that several of the statistics can reliably track changes in the underlying physics. Our analysis shows that the interactions between parameters are often among the most significant effects. A small fraction of statistics are also sensitive to changes in magnetic field properties. We use this framework to compare the set of simulations to observations of three nearby star-forming regions: NGC 1333, Oph A, and IC 348. We find that no one simulation agrees significantly better with the observations, although there is evidence that the high Mach number simulations are more consistent with the observations. △ Less

Submitted 17 July, 2017; originally announced July 2017.

Comments: 26 pages, 12 figures, accepted to MNRAS. If the data DOI link is not yet active, go here for the data: http://apps.canfar.net/storage/list/ekoch/Koch2017_SimulatedObs

Extreme Jet Ejections from the Black Hole X-ray Binary V404 Cygni

Authors: A. J. Tetarenko, G. R. Sivakoff, J. C. A Miller-Jones, E. W. Rosolowsky, G. Petitpas, M. Gurwell, J. Wouterloot, R. Fender, S. Heinz, D. Maitra, S. B. Markoff, S. Migliari, M. P. Rupen, A. P. Rushton, D. M. Russell, T. D. Russell, C. L. Sarazin

Abstract: We present simultaneous radio through sub-mm observations of the black hole X-ray binary (BHXB) V404 Cygni during the most active phase of its June 2015 outburst. Our $4$ hour long set of overlapping observations with the Very Large Array, the Sub-millimeter Array, and the James Clerk Maxwell Telescope (SCUBA-2), covers 8 different frequency bands (including the first detection of a BHXB jet at… ▽ More We present simultaneous radio through sub-mm observations of the black hole X-ray binary (BHXB) V404 Cygni during the most active phase of its June 2015 outburst. Our $4$ hour long set of overlapping observations with the Very Large Array, the Sub-millimeter Array, and the James Clerk Maxwell Telescope (SCUBA-2), covers 8 different frequency bands (including the first detection of a BHXB jet at $666 \,{\rm GHz}/450μm$), providing an unprecedented multi-frequency view of the extraordinary flaring activity seen during this period of the outburst. In particular, we detect multiple rapidly evolving flares, which reach Jy-level fluxes across all of our frequency bands. With this rich data set we performed detailed MCMC modeling of the repeated flaring events. Our custom model adapts the van der Laan synchrotron bubble model to include twin bi-polar ejections, propagating away from the black hole at bulk relativistic velocities, along a jet axis that is inclined to the line of sight. The emission predicted by our model accounts for projection effects, relativistic beaming, and the geometric time delay between the approaching and receding ejecta in each ejection event. We find that a total of 8 bi-polar, discrete jet ejection events can reproduce the emission that we observe in all of our frequency bands remarkably well. With our best fit model, we provide detailed probes of jet speed, structure, energetics, and geometry. Our analysis demonstrates the paramount importance of the mm/sub-mm bands, which offer a unique, more detailed view of the jet than can be provided by radio frequencies alone. △ Less

Submitted 27 April, 2017; originally announced April 2017.

Comments: 25 pages, 11 figures, Accepted to MNRAS

An Exploration of the Statistical Signatures of Stellar Feedback

Authors: Ryan D. Boyden, Eric W. Koch, Erik W. Rosolowsky, Stella S. R. Offner

Abstract: All molecular clouds are observed to be turbulent, but the origin, means of sustenance, and evolution of the turbulence remain debated. One possibility is that stellar feedback injects enough energy into the cloud to drive observed motions on parsec scales. Recent numerical studies of molecular clouds have found that feedback from stars, such as protostellar outflows and winds, injects energy and… ▽ More All molecular clouds are observed to be turbulent, but the origin, means of sustenance, and evolution of the turbulence remain debated. One possibility is that stellar feedback injects enough energy into the cloud to drive observed motions on parsec scales. Recent numerical studies of molecular clouds have found that feedback from stars, such as protostellar outflows and winds, injects energy and impacts turbulence. We expand upon these studies by analyzing magnetohydrodynamic simulations of molecular clouds, including stellar winds, with a range of stellar mass-loss rates and magnetic field strengths. We generate synthetic $^{12}$CO(1-0) maps assuming that the simulations are at the distance of the nearby Perseus molecular cloud. By comparing the outputs from different initial conditions and evolutionary times, we identify differences in the synthetic observations and characterize these using common astrostatistics. We quantify the different statistical responses using a variety of metrics proposed in the literature. We find that multiple astrostatistics, including principal component analysis, the spectral correlation function and the velocity coordinate spectrum, are sensitive to changes in stellar mass-loss rates and/or time evolution. A few statistics, including the Cramer statistic and velocity coordinate spectrum, are sensitive to the magnetic field strength. These findings demonstrate that stellar feedback influences molecular cloud turbulence and can be identified and quantified observationally using such statistics. △ Less

Submitted 2 November, 2016; originally announced November 2016.

Comments: 22 pages, 23 figures, accepted to ApJ

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

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

Filament Identification through Mathematical Morphology

Authors: Eric W. Koch, Erik W. Rosolowsky

Abstract: We present a new algorithm for detecting filamentary structure FilFinder. The algorithm uses the techniques of mathematical morphology for filament identification, presenting a complementary approach to current algorithms which use matched filtering or critical manifolds. Unlike other methods, FilFinder identifies filaments over a wide dynamic range in brightness. We apply the new algorithm to far… ▽ More We present a new algorithm for detecting filamentary structure FilFinder. The algorithm uses the techniques of mathematical morphology for filament identification, presenting a complementary approach to current algorithms which use matched filtering or critical manifolds. Unlike other methods, FilFinder identifies filaments over a wide dynamic range in brightness. We apply the new algorithm to far infrared imaging data of dust emission released by the Herschel Gould Belt Survey team. Our preliminary analysis characterizes both filaments and fainter striations. We find a typical filament width of 0.09 pc across the sample, but the brightness varies from cloud to cloud. Several regions show a bimodal filament brightness distribution, with the bright mode (filaments) being an order of magnitude brighter than the faint mode (striations). Using the Rolling Hough Transform, we characterize the orientations of the striations in the data, finding preferred directions that agree with magnetic field direction where data are available. There is a suggestive but noisy correlation between typical filament brightness and literature values of the star formation rates for clouds in the Gould Belt. △ Less

Submitted 8 July, 2015; originally announced July 2015.

Comments: 18 pages, 14 figures, Accepted for publication in MNRAS. The FilFinder algorithm is available at http://github.com/e-koch/FilFinder and documentation at http://fil-finder.readthedocs.org/en/latest/. Maps of the identified filaments are available at http://e-koch.github.io/FilFinder_Results/hgbs.html and as online figures

Astrochemical Correlations in Molecular Clouds

Authors: Brandt A. L. Gaches, Stella S. R. Offner, Erik W. Rosolowsky, Thomas G. Bisbas

Abstract: We investigate the spectral correlations between different species used to observe molecular clouds. We use hydrodynamic simulations and a full chemical network to study the abundances of over 150 species in typical Milky Way molecular clouds. We perform synthetic observations in order to produce emission maps of a subset of these tracers. We study the effects of different lines of sight and spati… ▽ More We investigate the spectral correlations between different species used to observe molecular clouds. We use hydrodynamic simulations and a full chemical network to study the abundances of over 150 species in typical Milky Way molecular clouds. We perform synthetic observations in order to produce emission maps of a subset of these tracers. We study the effects of different lines of sight and spatial resolution on the emission distribution and perform a robust quantitative comparison of the species to each other. We use the Spectral Correlation Function (SCF), which quantifies the root mean squared difference between spectra separated by some length scale, to characterize the structure of the simulated cloud in position-position-velocity (PPV) space. We predict the observed SCF for a broad range of observational tracers, and thus, identify homologous species. In particular, we show that the pairs C and CO, C$^{+}$ and CN, NH$_3$ and H$_2$CS have very similar SCFs. We measure the SCF slope variation as a function of beam size for all species and demonstrate that the beam size has a distinct effect on different species emission. However, for beams of up to 10", placing the cloud at 1 kpc, the change is not large enough to move the SCF slopes into different regions of parameter space. The results from this study provide observational guidance for choosing the best tracer to probe various cloud length scales. △ Less

Submitted 8 December, 2014; originally announced December 2014.

Comments: Accepted by ApJ. 19 pages, 13 figures

Power Spectrum Analysis of Polarized Emission from the Canadian Galactic Plane Survey

Authors: R. A. Stutz, E. W. Rosolowsky, R. Kothes, T. Landecker

Abstract: Angular power spectra are calculated and presented for the entirety of the Canadian Galactic Plane Survey polarization dataset at 1.4 GHz covering an area of 1060 deg$^2$. The data analyzed are a combination of data from the 100-m Effelsberg Telescope, the 26-m Telescope at the Dominion Radio Astrophysical Observatory, and the Synthesis Telescope at the Dominion Radio Astrophysical Observatory, al… ▽ More Angular power spectra are calculated and presented for the entirety of the Canadian Galactic Plane Survey polarization dataset at 1.4 GHz covering an area of 1060 deg$^2$. The data analyzed are a combination of data from the 100-m Effelsberg Telescope, the 26-m Telescope at the Dominion Radio Astrophysical Observatory, and the Synthesis Telescope at the Dominion Radio Astrophysical Observatory, allowing all scales to be sampled down to arcminute resolution. The resulting power spectra cover multipoles from $\ell \approx 60$ to $\ell \approx 10^4$ and display both a power-law component at low multipoles and a flattening at high multipoles from point sources. We fit the power spectrum with a model that accounts for these components and instrumental effects. The resulting power-law indices are found to have a mode of 2.3, similar to previous results. However, there are significant regional variations in the index, defying attempts to characterize the emission with a single value. The power-law index is found to increase away from the Galactic plane. A transition from small-scale to large-scale structure is evident at $b= 9^{\circ}$, associated with the disk-halo transition in a 15$^{\circ}$ region around $l=108^{\circ}$. Localized variations in the index are found toward HII regions and supernova remnants, but the interpretation of these variations is inconclusive. The power in the polarized emission is anticorrelated with bright thermal emission (traced by H$α$ emission) indicating that the thermal emission depolarizes background synchrotron emission. △ Less

Submitted 10 April, 2014; originally announced April 2014.

Comments: Accepted to ApJ; 17 pages

Direct observation of a sharp transition to coherence in Dense Cores

Authors: Jaime E. Pineda, Alyssa A. Goodman, Héctor G. Arce, Paola Caselli, Jonathan B. Foster, Philip C. Myers, Erik W. Rosolowsky

Abstract: We present NH3 observations of the B5 region in Perseus obtained with the Green Bank Telescope (GBT). The map covers a region large enough (~11'x14') that it contains the entire dense core observed in previous dust continuum surveys. The dense gas traced by NH3(1,1) covers a much larger area than the dust continuum features found in bolometer observations. The velocity dispersion in the central… ▽ More We present NH3 observations of the B5 region in Perseus obtained with the Green Bank Telescope (GBT). The map covers a region large enough (~11'x14') that it contains the entire dense core observed in previous dust continuum surveys. The dense gas traced by NH3(1,1) covers a much larger area than the dust continuum features found in bolometer observations. The velocity dispersion in the central region of the core is small, presenting subsonic non-thermal motions which are independent of scale. However, it is thanks to the coverage and high sensitivity of the observations that we present the detection, for the first time, of the transition between the coherent core and the dense but more turbulent gas surrounding it. This transition is sharp, increasing the velocity dispersion by a factor of 2 in less than 0.04 pc (the 31" beam size at the distance of Perseus, ~250 pc). The change in velocity dispersion at the transition is ~3 km/s/pc. The existence of the transition provides a natural definition of dense core: the region with nearly-constant subsonic non-thermal velocity dispersion. From the analysis presented here we can not confirm nor rule out a corresponding sharp density transition. △ Less

Submitted 15 February, 2010; originally announced February 2010.

Comments: Accepted for publication in ApJ Letters. 6 pages, 5 figures

The Effect of Projection on Derived Mass-Size and Linewidth-Size Relationships

Authors: Rahul Shetty, David C. Collins, Jens Kauffmann, Alyssa A. Goodman, Erik W. Rosolowsky, Michael L. Norman

Abstract: Power law mass-size and linewidth-size correlations, two of "Larson's laws," are often studied to assess the dynamical state of clumps within molecular clouds. Using the result of a hydrodynamic simulation of a molecular cloud, we investigate how geometric projection may affect the derived Larson relationships. We find that large scale structures in the column density map have similar masses and… ▽ More Power law mass-size and linewidth-size correlations, two of "Larson's laws," are often studied to assess the dynamical state of clumps within molecular clouds. Using the result of a hydrodynamic simulation of a molecular cloud, we investigate how geometric projection may affect the derived Larson relationships. We find that large scale structures in the column density map have similar masses and sizes to those in the 3D simulation (PPP). Smaller scale clumps in the column density map are measured to be more massive than the PPP clumps, due to the projection of all emitting gas along lines of sight. Further, due to projection effects, structures in a synthetic spectral observation (PPV) may not necessarily correlate with physical structures in the simulation. In considering the turbulent velocities only, the linewidth-size relationship in the PPV cube is appreciably different from that measured from the simulation. Including thermal pressure in the simulated linewidths imposes a minimum linewidth, which results in a better agreement in the slopes of the linewidth-size relationships, though there are still discrepancies in the offsets, as well as considerable scatter. Employing commonly used assumptions in a virial analysis, we find similarities in the computed virial parameters of the structures in the PPV and PPP cubes. However, due to the discrepancies in the linewidth- and mass- size relationships in the PPP and PPV cubes, we caution that applying a virial analysis to observed clouds may be misleading due to geometric projection effects. We speculate that consideration of physical processes beyond kinetic and gravitational pressure would be required for accurately assessing whether complex clouds, such as those with highly filamentary structure, are bound. △ Less

Submitted 10 February, 2010; v1 submitted 25 January, 2010; originally announced January 2010.

Comments: 25 pages, including 7 Figures; Accepted for publication in ApJ

The Perils of Clumpfind: The Mass Spectrum of Sub-structures in Molecular Clouds

Authors: Jaime E. Pineda, Erik W. Rosolowsky, Alyssa A. Goodman

Abstract: We study the mass spectrum of sub-structures in the Perseus Molecular Cloud Complex traced by 13CO (1-0), finding that $dN/dM\propto M^{-2.4}$ for the standard Clumpfind parameters. This result does not agree with the classical $dN/dM\propto M^{-1.6}$. To understand this discrepancy we study the robustness of the mass spectrum derived using the Clumpfind algorithm. Both 2D and 3D Clumpfind versi… ▽ More We study the mass spectrum of sub-structures in the Perseus Molecular Cloud Complex traced by 13CO (1-0), finding that $dN/dM\propto M^{-2.4}$ for the standard Clumpfind parameters. This result does not agree with the classical $dN/dM\propto M^{-1.6}$. To understand this discrepancy we study the robustness of the mass spectrum derived using the Clumpfind algorithm. Both 2D and 3D Clumpfind versions are tested, using 850 $μ$m dust emission and 13CO spectral-line observations of Perseus, respectively. The effect of varying threshold is not important, but varying stepsize produces a different effect for 2D and 3D cases. In the 2D case, where emission is relatively isolated (associated with only the densest peaks in the cloud), the mass spectrum variability is negligible compared to the mass function fit uncertainties. In the 3D case, however, where the 13CO emission traces the bulk of the molecular cloud, the number of clumps and the derived mass spectrum are highly correlated with the stepsize used. The distinction between "2D" and "3D" here is more importantly also a distinction between "sparse" and "crowded" emission. In any "crowded" case, Clumpfind should not be used blindly to derive mass functions. Clumpfind's output in the "crowded" case can still offer a statistical description of emission useful in inter-comparisons, but the clump-list should not be treated as a robust region decomposition suitable to generate a physically-meaningful mass function. We conclude that the 13CO mass spectrum depends on the observations resolution, due to the hierarchical structure of MC. △ Less

Submitted 1 June, 2009; originally announced June 2009.

Comments: 5 pages, 3 figures. Accepted for publication in ApJ Letters

Journal ref: Astrophys.J.699:L134-L138,2009

Dense Cores in Perseus: The Influence of Stellar Content and Cluster Environment

Authors: Jonathan B. Foster, Erik W. Rosolowsky, Jens Kauffman, Jaime E. Pineda, Michelle A. Borkin, Paola Caselli, Phil C. Myers, Alyssa A. Goodman

Abstract: We present the chemistry, temperature, and dynamical state of a sample of 193 dense cores or core candidates in the Perseus Molecular cloud and compare the properties of cores associated with young stars and clusters with those which are not. The combination of our NH3 and CCS observations with previous millimeter, sub-millimeter, and Spitzer data available for this cloud enable us both to deter… ▽ More We present the chemistry, temperature, and dynamical state of a sample of 193 dense cores or core candidates in the Perseus Molecular cloud and compare the properties of cores associated with young stars and clusters with those which are not. The combination of our NH3 and CCS observations with previous millimeter, sub-millimeter, and Spitzer data available for this cloud enable us both to determine core properties precisely and to accurately classify cores as starless or protostellar. The properties of cores in different cluster environments and before-and-after star formation provide important constraints on simulations of star-formation, particularly under the paradigm that the essence of star formation is set by the turbulent formation of prestellar cores. We separate the influence of stellar content from that of cluster environment and find that cores within clusters have (1) higher kinetic temperatures and (2) lower fractional abundances of CCS and NH3. Cores associated with protostars have (1) slightly higher kinetic temperatures (2) higher NH3 excitation temperatures), (3) are at higher column density, have (4) slightly more non-thermal/turbulent NH3 linewidths, have (5) higher masses and have (6) lower fractional abundance of CCS. We find that neither cluster environment nor protostellar content makes a significant difference to the dynamical state of cores as estimated by the virial parameter -- most cores in each category are gravitationally bound. Overall, cluster environment and protostellar content have a smaller influence on the properties of the cores than is typically assumed, and the variation within categories is larger than the differences between categories. △ Less

Submitted 3 February, 2009; originally announced February 2009.

Comments: 28 pages, 17 figures. Accepted to ApJ

Journal ref: Astrophys.J.696:298-319,2009

Structural Analysis of Molecular Clouds: Dendrograms

Authors: E. W. Rosolowsky, J. E. Pineda, J. Kauffmann, A. A. Goodman

Abstract: We demonstrate the utility of dendrograms at representing the essential features of the hierarchical structure of the isosurfaces for molecular line data cubes. The dendrogram of a data cube is an abstraction of the changing topology of the isosurfaces as a function of contour level. The ability to track hierarchical structure over a range of scales makes this analysis philosophically different… ▽ More We demonstrate the utility of dendrograms at representing the essential features of the hierarchical structure of the isosurfaces for molecular line data cubes. The dendrogram of a data cube is an abstraction of the changing topology of the isosurfaces as a function of contour level. The ability to track hierarchical structure over a range of scales makes this analysis philosophically different from local segmentation algorithms like CLUMPFIND. Points in the dendrogram structure correspond to specific volumes in data cubes defined by their bounding isosurfaces. We further refine the technique by measuring the properties associated with each isosurface in the analysis allowing for a multiscale calculation of molecular gas properties. Using COMPLETE 13CO(1-0) data from the L1448 region in Perseus and mock observations of a simulated data cube, we identify regions that have a significant contribution by self-gravity to their energetics on a range of scales. We find evidence for self-gravitation on all spatial scales in L1448 though not in all regions. In the simulated observations, nearly all of the emission is found in objects that would be self-gravitating if gravity were included in the simulation. We reconstruct the size-line width relationship within the data cube using the dendrogram-derived properties and find it follows the standard relation: s_v ~ R^0.58. Finally, we show that constructing the dendrogram of CO J=1-0 emission from the Orion-Monoceros region allows for the identification of giant molecular clouds in a blended molecular line data set using only a physically motivated definition (self-gravitating clouds with masses 5x10^4 Msun. △ Less

Submitted 20 February, 2008; originally announced February 2008.

Comments: 15 pages, 16 figures. Accepted to ApJ. Paper will full resolution figures available at http://people.ok.ubc.ca/erosolo/dendrograms.pdf

An Ammonia Spectral Atlas of Dense Cores in Perseus

Authors: E. W. Rosolowsky, J. E. Pineda, J. B. Foster, M. A. Borkin, J. Kauffmann, P. Caselli, P. C. Myers, A. A. Goodman

Abstract: We present ammonia observations of 193 dense cores and core candidates in the Perseus molecular cloud made using the Robert F. Byrd Green Bank Telescope. We simultaneously observed the NH3(1,1), NH3(2,2), CCS (2_1 -> 1_0) and CC34S (2_1 -> 1_0) transitions near 23 GHz for each of the targets with a spectral resolution of dv ~ 0.024 km/s. We find ammonia emission associated with nearly all of the… ▽ More We present ammonia observations of 193 dense cores and core candidates in the Perseus molecular cloud made using the Robert F. Byrd Green Bank Telescope. We simultaneously observed the NH3(1,1), NH3(2,2), CCS (2_1 -> 1_0) and CC34S (2_1 -> 1_0) transitions near 23 GHz for each of the targets with a spectral resolution of dv ~ 0.024 km/s. We find ammonia emission associated with nearly all of the (sub)millimeter sources as well as at several positions with no associated continuum emission. For each detection, we have measured physical properties by fitting a simple model to every spectral line simultaneously. Where appropriate, we have refined the model by accounting for low optical depths, multiple components along the line of sight and imperfect coupling to the GBT beam. For the cores in Perseus, we find a typical kinetic temperature of T=11 K, a typical column density of N(NH3)~ 10^14.5 /cm^2 and velocity dispersions ranging from sigma_v = 0.07 km/s to 0.7 km/s. However, many cores with velocity dispersions > 0.2 km/s show evidence for multiple velocity components along the line of sight. △ Less

Submitted 2 November, 2007; originally announced November 2007.

Comments: 19 pages; Accepted to ApJS; version with high resolution figures available at http://www.cfa.harvard.edu/COMPLETE/papers/nh3-paper1.pdf ; online data at http://www.cfa.harvard.edu/COMPLETE/data_html_pages/GBT_NH3.html

The Effects of Noise and Sampling on the Spectral Correlation Function

Authors: Paolo Padoan, Erik W. Rosolowsky, Alyssa A. Goodman

Abstract: The effects of noise and sampling on the ``Spectral Correlation Function'' (SCF) introduced by Rosolowsky et al. 1999 are studied using observational data, numerical simulations of magneto-hydrodynamic turbulence, and simple models of Gaussian spectral line profiles. The most significant innovations of this paper are: i) the normalization of the SCF based on an analytic model for the effect of n… ▽ More The effects of noise and sampling on the ``Spectral Correlation Function'' (SCF) introduced by Rosolowsky et al. 1999 are studied using observational data, numerical simulations of magneto-hydrodynamic turbulence, and simple models of Gaussian spectral line profiles. The most significant innovations of this paper are: i) the normalization of the SCF based on an analytic model for the effect of noise; ii) the computation of the SCF as a function of the spatial lag between spectra within a map. By computing of the dependence of the SCF on the spatial lag, S_o(dr), we have been able to conclude that: 1) S_o(dr) is a power law, with slope in the range of scales l_i-l_o. 2) The correlation outer scale, l_o, is determined by the size of the map, and no evidence for a true departure from self-similarity on large scales has been found. 3) The correlation inner scale, l_i, is a true estimate of the smallest self-similar scale in a map. 4) The spectral slope in a given region, is independent of velocity resolution (above a minimum resolution threshold), spatial resolution, and average spectrum quality. 5) Molecular transitions which trace higher gas density yield larger values of the spectral slope than transitions tracing lower gas density. 6) Nyquist sampling, bad pixels in detector arrays, and reference sharing data acquisition need to be taken into account for a correct determination of the SCF at lag dr=1. The value of the spectral slope, however, can be computed correctly without a detailed knowledge of observational procedures. △ Less

Submitted 17 October, 2000; originally announced October 2000.

Comments: 28 pages, 12 figures included. To be published in The Astrophysical Journal

The Spectral Correlation Function -- A New Tool for Analyzing Spectral-Line Maps

Authors: Erik W. Rosolowsky, Alyssa A. Goodman, David J. Wilner, Jonathan P. Williams

Abstract: The "spectral correlation function" analysis we introduce in this paper is a new tool for analyzing spectral-line data cubes. Our initial tests, carried out on a suite of observed and simulated data cubes, indicate that the spectral correlation function [SCF] is likely to be a more discriminating statistic than other statistical methods normally applied. The SCF is a measure of similarity betwee… ▽ More The "spectral correlation function" analysis we introduce in this paper is a new tool for analyzing spectral-line data cubes. Our initial tests, carried out on a suite of observed and simulated data cubes, indicate that the spectral correlation function [SCF] is likely to be a more discriminating statistic than other statistical methods normally applied. The SCF is a measure of similarity between neighboring spectra in the data cube. When the SCF is used to compare a data cube consisting of spectral-line observations of the ISM with a data cube derived from MHD simulations of molecular clouds, it can find differences that are not found by other analyses. The initial results presented here suggest that the inclusion of self-gravity in numerical simulations is critical for reproducing the correlation behavior of spectra in star-forming molecular clouds. △ Less

Submitted 20 April, 1999; v1 submitted 30 March, 1999; originally announced March 1999.

Comments: 29 pages, including 4 figures (tar file submitted as source) See also: http://cfa-www.harvard.edu/~agoodman/scf/velocity_methods.html