Showing 1–50 of 118 results for author: Stanimirovic, S

Revisiting the Vertical Distribution of HI Absorbing Clouds in the Solar Neighborhood

Authors: Trey V. Wenger, Daniel R. Rybarczyk, Snežana Stanimirović

Abstract: The vertical distribution of cold neutral hydrogen (HI) clouds is a constraint on models of the structure, dynamics, and hydrostatic balance of the interstellar medium. In 1978, Crovisier pioneered a method to infer the vertical distribution of HI absorbing clouds in the solar neighborhood. Using data from the Nançay 21-cm absorption survey, they determine the mean vertical displacement of cold HI… ▽ More The vertical distribution of cold neutral hydrogen (HI) clouds is a constraint on models of the structure, dynamics, and hydrostatic balance of the interstellar medium. In 1978, Crovisier pioneered a method to infer the vertical distribution of HI absorbing clouds in the solar neighborhood. Using data from the Nançay 21-cm absorption survey, they determine the mean vertical displacement of cold HI clouds, $\langle|z|\rangle$. We revisit Crovisier's analysis and explore the consequences of truncating the HI absorption sample in Galactic latitude. For any non-zero latitude limit, we find that the quantity inferred by Crovisier is not the mean vertical displacement but rather a ratio involving higher moments of the vertical distribution. The resultant distribution scale heights are thus ${\sim}1.5$ to ${\sim}3$ times smaller than previously determined. In light of this discovery, we develop a Bayesian Monte Carlo Markov Chain method to infer the vertical distribution of HI absorbing clouds. We fit our model to the original Nançay data and find a vertical distribution moment ratio $\langle|z|^3\rangle/\langle|z|^2\rangle = 97 \pm 15\,\text{pc}$, which corresponds to a Gaussian scale height $σ_z = 61 \pm 9\,\text{pc}$, an exponential scale height $λ_z = 32 \pm 5\,\text{pc}$, and a rectangular half-width $W_{z, 1/2} = 129 \pm 20\,\text{pc}$. Consistent with recent simulations, the vertical scale height of cold HI clouds appears to remain constant between the inner-Galaxy and the Galactocentric distance of the solar neighborhood. Local fluctuations might explain the large scale height observed at the same Galactocentric distance on the far side of the Galaxy. △ Less

Submitted 27 March, 2024; originally announced March 2024.

Comments: 13 pages, 3 figures, 3 tables. Accepted for publication in the Astrophysical Journal

Where do stars explode in the ISM? -- The distribution of dense gas around massive stars and supernova remnants in M33

Authors: Sumit K. Sarbadhicary, Jordan Wagner, Eric W. Koch, Ness Mayker Chen, Adam K. Leroy, Natalia Lahén, Erik Rosolowsky, Kathryn F. Neugent, Chang-Goo Kim, Laura Chomiuk, Julianne J. Dalcanton, Laura A. Lopez, Nickolas M. Pingel, Remy Indebetouw, Thomas G. Williams, Elizabeth Tarantino, Jennifer Donovan Meyer, Evan D. Skillman, Adam Smercina, Amanda A. Kepley, Eric J. Murphy, Jay Strader, Tony Wong, Snežana Stanimirović, Vicente Villanueva , et al. (5 additional authors not shown)

Abstract: Star formation in galaxies is regulated by turbulence, outflows, gas heating and cloud dispersal -- processes which depend sensitively on the properties of the interstellar medium (ISM) into which supernovae (SNe) explode. Unfortunately, direct measurements of ISM environments around SNe remain scarce, as SNe are rare and often distant. Here we demonstrate a new approach: mapping the ISM around th… ▽ More Star formation in galaxies is regulated by turbulence, outflows, gas heating and cloud dispersal -- processes which depend sensitively on the properties of the interstellar medium (ISM) into which supernovae (SNe) explode. Unfortunately, direct measurements of ISM environments around SNe remain scarce, as SNe are rare and often distant. Here we demonstrate a new approach: mapping the ISM around the massive stars that are soon to explode. This provides a much larger census of explosion sites than possible with only SNe, and allows comparison with sensitive, high-resolution maps of the atomic and molecular gas from the Jansky VLA and ALMA. In the well-resolved Local Group spiral M33, we specifically observe the environments of red supergiants (RSGs, progenitors of Type II SNe), Wolf-Rayet stars (WRs, tracing stars $>$30 M$_{\odot}$, and possibly future stripped-envelope SNe), and supernova remnants (SNRs, locations where SNe have exploded). We find that massive stars evolve not only in dense, molecular-dominated gas (with younger stars in denser gas), but also a substantial fraction ($\sim$45\% of WRs; higher for RSGs) evolve in lower-density, atomic-gas-dominated, inter-cloud media. We show that these measurements are consistent with expectations from different stellar-age tracer maps, and can be useful for validating SN feedback models in numerical simulations of galaxies. Along with the discovery of a 20-pc diameter molecular gas cavity around a WR, these findings re-emphasize the importance of pre-SN/correlated-SN feedback evacuating the dense gas around massive stars before explosion, and the need for high-resolution (down to pc-scale) surveys of the multi-phase ISM in nearby galaxies. △ Less

Submitted 26 October, 2023; originally announced October 2023.

Comments: 34 pages, 14 figures. Submitted to ApJ. Comments welcome! The density distributions will be made publicly available after journal acceptance of manuscript. Please feel free to contact us in the meantime if you would like to use them

Neutral Hydrogen (HI) 21 cm as a probe: Investigating Spatial Variations in Interstellar Turbulent Properties

Authors: Amit K. Mittal, Brian L Babler, Snezana Stanimirovic, Nickolas Pingel

Abstract: Interstellar turbulence shapes the HI distribution in the Milky Way (MW). How this affects large-scale statistical properties of HI column density across the MW remains largely unconstrained. We use approx 13,000 square-degree GALFA-HI survey to map statistical fluctuations of HI over the 40 km s-1 velocity range. We calculate the spatial power spectrum (SPS) of HI column density image by running… ▽ More Interstellar turbulence shapes the HI distribution in the Milky Way (MW). How this affects large-scale statistical properties of HI column density across the MW remains largely unconstrained. We use approx 13,000 square-degree GALFA-HI survey to map statistical fluctuations of HI over the 40 km s-1 velocity range. We calculate the spatial power spectrum (SPS) of HI column density image by running a 3-degree kernel and measuring SPS slope over a range of angular scales from 16 arcmin to 20 degree. Due to GALFA complex observing and calibration strategy, we construct detailed estimates of the noise contribution and account for GALFA beam effects on SPS. This allows us to systematically analyze HI images that trace a wide range of interstellar environments. We find that SPS slope varies between -2.6 at high Galactic latitudes, and -3.2 close to Galactic plane. The range of SPS slope values becomes tighter when we consider HI optical depth and line-of-sight length caused by the plane-parallel geometry of HI disk. This relatively uniform, large-scale distribution of SPS slope is suggestive of large-scale turbulent driving being a dominant mechanism for shaping HI structures in the MW and/or the stellar feedback turbulence being efficiently dissipated within dense molecular clouds. Only at latitudes above 60 degrees we find evidence for HI SPS slope being consistently more shallow. Those directions are largely within the Local Bubble, suggesting the recent history of this cavity, shaped by multiple supernovae explosions, has modified the turbulent state of HI and/or fractions of HI phases. △ Less

Submitted 21 October, 2023; originally announced October 2023.

Comments: Accepted in ApJ

A new method for spatially resolving the turbulence driving mixture in the ISM with application to the Small Magellanic Cloud

Authors: Isabella A. Gerrard, Christoph Federrath, Nickolas M. Pingel, Naomi M. McClure-Griffiths, Antoine Marchal, Gilles Joncas, Susan E. Clark, Snežana Stanimirović, Min-Young Lee, Jacco Th. van Loon, John Dickey, Helga Dénes, Yik Ki Ma, James Dempsey, Callum Lynn

Abstract: Turbulence plays a crucial role in shaping the structure of the interstellar medium. The ratio of the three-dimensional density contrast ($σ_{ρ/ρ_0}$) to the turbulent sonic Mach number ($\mathcal{M}$) of an isothermal, compressible gas describes the ratio of solenoidal to compressive modes in the turbulent acceleration field of the gas, and is parameterised by the turbulence driving parameter:… ▽ More Turbulence plays a crucial role in shaping the structure of the interstellar medium. The ratio of the three-dimensional density contrast ($σ_{ρ/ρ_0}$) to the turbulent sonic Mach number ($\mathcal{M}$) of an isothermal, compressible gas describes the ratio of solenoidal to compressive modes in the turbulent acceleration field of the gas, and is parameterised by the turbulence driving parameter: $b=σ_{ρ/ρ_0}/\mathcal{M}$. The turbulence driving parameter ranges from $b=1/3$ (purely solenoidal) to $b=1$ (purely compressive), with $b=0.38$ characterising the natural mixture (1/3~compressive, 2/3~solenoidal) of the two driving modes. Here we present a new method for recovering $σ_{ρ/ρ_0}$, $\mathcal{M}$, and $b$, from observations on galactic scales, using a roving kernel to produce maps of these quantities from column density and centroid velocity maps. We apply our method to high-resolution HI emission observations of the Small Magellanic Cloud (SMC) from the GASKAP-HI survey. We find that the turbulence driving parameter varies between $b\sim 0.3$ and $b\sim 1.0$ within the main body of the SMC, but the median value converges to $b\sim0.51$, suggesting that the turbulence is overall driven more compressively ($b>0.38$). We observe no correlation between the $b$ parameter and HI or H$α$ intensity, indicating that compressive driving of HI turbulence cannot be determined solely by observing HI or H$α$ emission density, and that velocity information must also be considered. Further investigation is required to link our findings to potential driving mechanisms such as star-formation feedback, gravitational collapse, or cloud-cloud collisions. △ Less

Submitted 19 September, 2023; originally announced September 2023.

Comments: 20 pages, 16 figures, accepted to MNRAS

Probing the Conditions for the Hı-to-H$_{2}$ Transition in the Interstellar Medium

Authors: Gyueun Park, Min-Young Lee, Shmuel Bialy, Blakesley Burkhart, J. R. Dawson, Carl Heiles, Di Li, Claire Murray, Hiep Nguyen, Anita Hafner, Daniel R. Rybarczyk, Snežana Stanimirović

Abstract: In this paper, we investigate the conditions for the HI-to-H$_{2}$ transition in the solar neighborhood by analyzing HI emission and absorption measurements toward 58 Galactic lines of sight (LOSs) along with $^{12}$CO(1$-$0) (CO) and dust data. Based on the accurate column densities of the cold and warm neutral medium (CNM and WNM), we first perform a decomposition of gas into atomic and molecula… ▽ More In this paper, we investigate the conditions for the HI-to-H$_{2}$ transition in the solar neighborhood by analyzing HI emission and absorption measurements toward 58 Galactic lines of sight (LOSs) along with $^{12}$CO(1$-$0) (CO) and dust data. Based on the accurate column densities of the cold and warm neutral medium (CNM and WNM), we first perform a decomposition of gas into atomic and molecular phases and show that the observed LOSs are mostly HI-dominated. In addition, we find that the CO-dark H$_{2}$, not the optically thick HI, is a major ingredient of the dark gas in the solar neighborhood. To examine the conditions for the formation of CO-bright molecular gas, we analyze the kinematic association between HI and CO and find that the CNM is kinematically more closely associated with CO than the WNM. When CNM components within CO line widths are isolated, we find the following characteristics: spin temperature $<$ 200 K, peak optical depth $>$ 0.1, CNM fraction of $\sim$0.6, and $V$-band dust extinction $>$ 0.5 mag. These results suggest that CO-bright molecular gas preferentially forms in environments with high column densities where the CNM becomes colder and more abundant. Finally, we confront the observed CNM properties with the steady-state H$_{2}$ formation model of Sternberg et al. and infer that the CNM must be clumpy with a small volume filling factor. Another possibility would be that missing processes in the model, such as cosmic-rays and gas dynamics, play an important role in the HI-to-H$_{2}$ transition. △ Less

Submitted 18 July, 2023; originally announced July 2023.

Comments: 25 pages, 12 figures, 6 tables, Accepted for publication in ApJ

Atomic Hydrogen in the Milky Way: A Stepping Stone in the Evolution of Galaxies

Authors: Naomi M. McClure-Griffiths, Snezana Stanimirovic, Daniel R. Rybarczyk

Abstract: Atomic hydrogen (HI) is a critical stepping stone in the gas evolution cycle of the interstellar medium (ISM) of the Milky Way. Hi traces both the cold, premolecular state before star formation and the warm, diffuse ISM before and after star formation. This review describes new, sensitive HI absorption and emission surveys, which, together with high angular and spectral resolution Hi emission data… ▽ More Atomic hydrogen (HI) is a critical stepping stone in the gas evolution cycle of the interstellar medium (ISM) of the Milky Way. Hi traces both the cold, premolecular state before star formation and the warm, diffuse ISM before and after star formation. This review describes new, sensitive HI absorption and emission surveys, which, together with high angular and spectral resolution Hi emission data, have revealed the physical properties of HI, its structure, and its association with magnetic fields. We give an overview of the HI phases and discuss how Hi properties depend on the environment and what its structure can tell us about feedback in the ISM. Key findings include the following: - The mass fraction of the cold neutral medium is $\lesssim 40$\% on average, increasing with $A_V$ due to the increase of mean gas density. - The cold disk extends to at least $R\sim 25$ kpc. - Approximately 40% of the HI is warm, with structural characteristics that derive from feedback events. - Cold HI is highly filamentary, whereas warm HI is more smoothly distributed. We summarize future observational and simulation opportunities that can be used to unravel the 3D structure of the atomic ISM and the effects of heating and cooling on HI properties. △ Less

Submitted 17 July, 2023; originally announced July 2023.

Comments: Manuscript (48 pages, 11 figures) accepted to Annual Reviews of Astronomy and Astrophysics for publication in Volume 61. This is the authors' own version. The final version, and associated supplementary material, will become available from: https://doi.org/10.1146/annurev-astro-052920-104851

Shaken or stirred: the diffuse interstellar medium with exceptionally high SiO abundance

Authors: Daniel R. Rybarczyk, Snezana Stanimirovic, Antoine Gusdorf

Abstract: Interstellar shocks, a key element of stellar feedback processes, shape the structure of the interstellar medium (ISM) and are essential for the chemistry, thermodynamics, and kinematics of interstellar gas. Powerful, high-velocity shocks are driven by stellar winds, young supernova explosions, more evolved supernova remnants, cloud-cloud collisions, and protostellar outflows, whereas the existenc… ▽ More Interstellar shocks, a key element of stellar feedback processes, shape the structure of the interstellar medium (ISM) and are essential for the chemistry, thermodynamics, and kinematics of interstellar gas. Powerful, high-velocity shocks are driven by stellar winds, young supernova explosions, more evolved supernova remnants, cloud-cloud collisions, and protostellar outflows, whereas the existence and origin of much-lower-velocity shocks ($\lesssim $ 10 km$~$s$^{-1}$) are not understood. Direct observational evidence for interstellar shocks in diffuse and translucent ISM environments have been especially lacking. We present the most sensitive survey to date of SiO -- often considered an unambiguous tracer of interstellar shocks -- in absorption, obtained with the Northern Extended Millimeter Array interferometer. We detect SiO in 5/8 directions probing diffuse and translucent environments without ongoing star formation. Our results demonstrate that SiO formation in the diffuse ISM (i.e., in the absence of significant star formation and stellar feedback) is more widespread and effective than previously reported. The observed SiO linewidths are all $\lesssim$ 4 km$~$s$^{-1}$, excluding high-velocity shocks as a formation mechanism. Yet, the SiO abundances we detect are mostly 1 to 2 orders of magnitude higher than typically assumed in quiescent environments and are often accompanied with other molecular transitions whose column densities cannot be explained with UV-dominated chemical models. Our results challenge the traditional view of SiO production via stellar-feedback sources and emphasize the need for observational constraints on the distribution of Si in the gas phase and grain mantles, which are crucial for understanding the physics of grain processing and diffuse interstellar chemistry. △ Less

Submitted 13 April, 2023; originally announced April 2023.

Comments: 19 pages, 6 figures. Accepted to ApJ

GNOMES II: Analysis of the Galactic diffuse molecular ISM in all four ground state hydroxyl transitions using Amoeba

Authors: Anita Petzler, J. R. Dawson, Hiep Nguyen, Carl Heiles, M. Wardle, M. -Y. Lee, Claire E. Murray, K. L. Thompson, Snezana Stanimirovic

Abstract: We present observations of the four 2 Pi 3/2 J = 3/2 ground-rotational state transitions of the hydroxyl molecule (OH) along 107 lines of sight both in and out of the Galactic plane: 92 sets of observations from the Arecibo telescope and 15 sets of observations from the Australia Telescope Compact Array (ATCA). Our Arecibo observations included off-source pointings, allowing us to measure excitati… ▽ More We present observations of the four 2 Pi 3/2 J = 3/2 ground-rotational state transitions of the hydroxyl molecule (OH) along 107 lines of sight both in and out of the Galactic plane: 92 sets of observations from the Arecibo telescope and 15 sets of observations from the Australia Telescope Compact Array (ATCA). Our Arecibo observations included off-source pointings, allowing us to measure excitation temperature (Tex) and optical depth, while our ATCA observations give optical depth only. We perform Gaussian decomposition using the Automated Molecular Excitation Bayesian line-fitting Algorithm 'AMOEBA' (Petzler, Dawson, and Wardle 2021) fitting all four transitions simultaneously with shared centroid velocity and width. We identify 109 features across 38 sightlines (including 58 detections along 27 sightlines with excitation temperature measurements). While the main lines at 1665 and 1667 MHz tend to have similar excitation temperatures (median Tex(main) difference = 0.6 K, 84% show Tex(main) difference < 2 K), large differences in the 1612 and 1720 MHz satellite line excitation temperatures show that the gas is generally not in LTE. For a selection of sightlines we compare our OH features to associated (on-sky and in velocity) HI cold gas components (CNM) identified by Nguyen et al. (2019) and find no strong correlations. We speculate that this may indicate an effective decoupling of the molecular gas from the CNM once it accumulates. △ Less

Submitted 21 February, 2023; originally announced February 2023.

Comments: Accepted for publication to PASA. 41 pages, 27 figures

GASKAP-HI Pilot Survey Science III: An unbiased view of cold gas in the Small Magellanic Cloud

Authors: James Dempsey, N. M. McClure-Griffiths, Claire Murray, John M. Dickey, Nickolas M. Pingel, Katherine Jameson, Helga Dénes, Jacco Th. van Loon, D. Leahy, Min-Young Lee, S. Stanimirović, Shari Breen, Frances Buckland-Willis, Steven J. Gibson, Hiroshi Imai, Callum Lynn, C. D. Tremblay

Abstract: We present the first unbiased survey of neutral hydrogen (HI) absorption in the Small Magellanic Cloud (SMC). The survey utilises pilot HI observations with the Australian Square Kilometre Array Pathfinder (ASKAP) telescope as part of the Galactic ASKAP HI (GASKAP-HI) project whose dataset has been processed with the GASKAP-HI absorption pipeline, also described here. This dataset provides absorpt… ▽ More We present the first unbiased survey of neutral hydrogen (HI) absorption in the Small Magellanic Cloud (SMC). The survey utilises pilot HI observations with the Australian Square Kilometre Array Pathfinder (ASKAP) telescope as part of the Galactic ASKAP HI (GASKAP-HI) project whose dataset has been processed with the GASKAP-HI absorption pipeline, also described here. This dataset provides absorption spectra towards 229 continuum sources, a 275% increase in the number of continuum sources previously published in the SMC region, as well as an improvement in the quality of absorption spectra over previous surveys of the SMC. Our unbiased view, combined with the closely matched beam size between emission and absorption, reveals a lower cold gas faction (11%) than the 2019 ATCA survey of the SMC and is more representative of the SMC as a whole. We also find that the optical depth varies greatly between the SMC's bar and wing regions. In the bar we find that the optical depth is generally low (correction factor to the optically thin column density assumption of $\mathcal{R}_{\rm HI} \sim 1.04$) but increases linearly with column density. In the wing however, there is a wide scatter in optical depth despite a tighter range of column densities. △ Less

Submitted 13 April, 2022; originally announced April 2022.

Comments: Accepted for publication in PASA, 19 pages, 17 figures, 5 tables

The role of neutral hydrogen in setting the abundances of molecular species in the Milky Way's diffuse interstellar medium. II. Comparison between observations and theoretical models

Authors: Daniel R. Rybarczyk, Munan Gong, Snezana Stanimirovic, Brian Babler, Claire E. Murray, Jan Martin Winters, Gan Luo, T. M. Dame, Lucille Steffes

Abstract: We compare observations of HI from the Very Large Array (VLA) and the Arecibo Observatory and observations of HCO$^+$ from the Atacama Large Millimeter/submillimeter Array (ALMA) and the Northern Extended Millimeter Array (NOEMA) in the diffuse ($A_V\lesssim1$) interstellar medium (ISM) to predictions from a photodissociation region (PDR) chemical model and multi-phase ISM simulations. Using a coa… ▽ More We compare observations of HI from the Very Large Array (VLA) and the Arecibo Observatory and observations of HCO$^+$ from the Atacama Large Millimeter/submillimeter Array (ALMA) and the Northern Extended Millimeter Array (NOEMA) in the diffuse ($A_V\lesssim1$) interstellar medium (ISM) to predictions from a photodissociation region (PDR) chemical model and multi-phase ISM simulations. Using a coarse grid of PDR models, we estimate the density, FUV radiation field, and cosmic ray ionization rate (CRIR) for each structure identified in HCO$^+$ and HI absorption. These structures fall into two categories. Structures with $T_s<40~\mathrm{K}$, mostly with $N(\mathrm{HCO^+})\lesssim10^{12}~\mathrm{cm^{-2}}$, are consistent with modest density, FUV radiation field, and CRIR models, typical of the diffuse molecular ISM. Structures with spin temperature $T_s>40~\mathrm{K}$, mostly with $N(\mathrm{HCO^+})\gtrsim10^{12}~\mathrm{cm^{-2}}$, are consistent with high density, FUV radiation field, and CRIR models, characteristic of environments close to massive star formation. The latter are also found in directions with a significant fraction of thermally unstable HI. In at least one case, we rule out the PDR model parameters, suggesting that alternative mechanisms (e.g., non-equilibrium processes like turbulent dissipation and/or shocks) are required to explain the observed HCO$^+$ in this direction. Similarly, while our observations and simulations of the turbulent, multi-phase ISM agree that HCO$^+$ formation occurs along sightlines with $N(\mathrm{HI})\gtrsim10^{21}~\mathrm{cm^{-2}}$, the simulated data fail to explain HCO$^+$ column densities $\gtrsim\rm{few}\times10^{12}~\mathrm{cm^{-2}}$. Since a majority of our sightlines with HCO$^+$ had such high column densities, this likely indicates that non-equilibrium chemistry is important for these lines of sight. △ Less

Submitted 10 December, 2021; originally announced December 2021.

Comments: 15 pages, 3 figures, 1 table. Accepted to ApJ

GASKAP-HI Pilot Survey Science I: ASKAP Zoom Observations of HI Emission in the Small Magellanic Cloud

Authors: N. M. Pingel, J. Dempsey, N. M. McClure-Griffiths, J. M. Dickey, K. E. Jameson, H. Arce, G. Anglada, J. Bland-Hawthorn, S. L. Breen, F. Buckland-Willis, S. E. Clark, J. R. Dawson, H. Dénes, E. M. Di Teodoro, B. -Q. For, Tyler J. Foster, J. F. Gómez, H. Imai, G. Joncas, C. -G. Kim, M. -Y. Lee, C. Lynn, D. Leahy, Y. K. Ma, A. Marchal , et al. (31 additional authors not shown)

Abstract: We present the most sensitive and detailed view of the neutral hydrogen (HI) emission associated with the Small Magellanic Cloud (SMC), through the combination of data from the Australian Square Kilometre Array Pathfinder (ASKAP) and Parkes (Murriyang), as part of the Galactic Australian Square Kilometre Array Pathfinder (GASKAP) pilot survey. These GASKAP-HI pilot observations, for the first time… ▽ More We present the most sensitive and detailed view of the neutral hydrogen (HI) emission associated with the Small Magellanic Cloud (SMC), through the combination of data from the Australian Square Kilometre Array Pathfinder (ASKAP) and Parkes (Murriyang), as part of the Galactic Australian Square Kilometre Array Pathfinder (GASKAP) pilot survey. These GASKAP-HI pilot observations, for the first time, reveal HI in the SMC on similar physical scales as other important tracers of the interstellar medium, such as molecular gas and dust. The resultant image cube possesses an rms noise level of 1.1 K (1.6 mJy/beam) per 0.98 km s$^{-1}$ spectral channel with an angular resolution of 30$''$ ($\sim$10 pc). We discuss the calibration scheme and the custom imaging pipeline that utilizes a joint deconvolution approach, efficiently distributed across a computing cluster, to accurately recover the emission extending across the entire $\sim$25 deg$^2$ field-of-view. We provide an overview of the data products and characterize several aspects including the noise properties as a function of angular resolution and the represented spatial scales by deriving the global transfer function over the full spectral range. A preliminary spatial power spectrum analysis on individual spectral channels reveals that the power-law nature of the density distribution extends down to scales of 10 pc. We highlight the scientific potential of these data by comparing the properties of an outflowing high velocity cloud with previous ASKAP+Parkes HI test observations. △ Less

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

Comments: Accepted for publication in PASA, 34 pages, 18 figures, 5 tables

GASKAP Pilot Survey Science II: ASKAP Zoom Observations of Galactic 21-cm Absorption

Authors: J. M. Dickey, J. M. Dempsey, N. M. Pingel, N. M. McClure-Griffiths, K. Jameson, J. R. Dawson, H. Dénes, S. E. Clark, G. Joncas, D. Leahy, Min-Young Lee, M. -A. Miville-Deschênes, S. Stanimirović, C. D. Tremblay, J. Th. van Loon

Abstract: Using the Australian Square Kilometre Array Pathfinder to measure 21-cm absorption spectra toward continuum background sources, we study the cool phase of the neutral atomic gas in the far outer disk, and in the inner Galaxy near the end of the Galactic bar at longitude 340 degrees. In the inner Galaxy the cool atomic gas has a smaller scale height than in the solar neighborhood, similar to the mo… ▽ More Using the Australian Square Kilometre Array Pathfinder to measure 21-cm absorption spectra toward continuum background sources, we study the cool phase of the neutral atomic gas in the far outer disk, and in the inner Galaxy near the end of the Galactic bar at longitude 340 degrees. In the inner Galaxy the cool atomic gas has a smaller scale height than in the solar neighborhood, similar to the molecular gas and the superthin stellar population in the bar. In the outer Galaxy the cool atomic gas is mixed with the warm, neutral medium, with the cool fraction staying roughly constant with Galactic radius. The mean spin temperature, i.e. the ratio of the emission brightness temperature to the absorption, is roughly constant for velocities corresponding to Galactic radius greater than about twice the solar circle radius. The ratio has a value of about 300 K, but this does not correspond to a physical temperature in the gas. If the gas causing the absorption has kinetic temperature of about 100 K, as in the solar neighborhood, then the value 300 K indicates that the fraction of the gas mass in this phase is one-third of the total HI mass. △ Less

Submitted 26 February, 2022; v1 submitted 8 November, 2021; originally announced November 2021.

Comments: 29 pages, 12 figures, submitted to Astrophysical Journal

The role of neutral hydrogen in setting the abundances of molecular species in the Milky Way's diffuse interstellar medium. I. Observational constraints from ALMA and NOEMA

Authors: Daniel R. Rybarczyk, Snezana Stanimirovic, Munan Gong, Brian Babler, Claire E. Murray, Maryvonne Gerin, Jan Martin Winters, Gan Luo, T. M. Dame, Lucille Steffes

Abstract: We have complemented existing observations of HI absorption with new observations of HCO$^+$, C$_2$H, HCN, and HNC absorption from the Atacama Large Millimeter/submillimeter Array (ALMA) and the Northern Extended Millimeter Array (NOEMA) in the direction of 20 background radio continuum sources with $4^\circ \leq |b| \leq 81^\circ$ to constrain the atomic gas conditions that are suitable for the f… ▽ More We have complemented existing observations of HI absorption with new observations of HCO$^+$, C$_2$H, HCN, and HNC absorption from the Atacama Large Millimeter/submillimeter Array (ALMA) and the Northern Extended Millimeter Array (NOEMA) in the direction of 20 background radio continuum sources with $4^\circ \leq |b| \leq 81^\circ$ to constrain the atomic gas conditions that are suitable for the formation of diffuse molecular gas. We find that these molecular species form along sightlines where $A_V \gtrsim 0.25$, consistent with the threshold for the HI-to-H$_2$ transition at solar metallicity. Moreover, we find that molecular gas is associated only with structures that have an HI optical depth $> 0.1$, a spin temperature $< 80$ K, and a turbulent Mach number $\gtrsim 2$. We also identify a broad, faint component to the HCO$^+$ absorption in a majority of sightlines. Compared to the velocities where strong, narrow HCO$^+$ absorption is observed, the HI at these velocities has a lower cold neutral medium (CNM) fraction and negligible CO emission. The relative column densities and linewidths of the different molecular species observed here are similar to those observed in previous experiments over a range of Galactic latitudes, suggesting that gas in the solar neighborhood and gas in the Galactic plane are chemically similar. For a select sample of previously-observed sightlines, we show that the absorption line profiles of HCO$^+$, HCN, HNC, and C$_2$H are stable over periods of $\sim 3$ years and $\sim 25$ years, likely indicating that molecular gas structures in these directions are at least $\gtrsim 100$ AU in size △ Less

Submitted 13 September, 2021; originally announced September 2021.

Comments: 31 pages, 12 figures, 8 tables. Accepted to ApJ

The MACH HI absorption survey I: Physical conditions of cold atomic gas outside of the Galactic plane

Authors: Claire E. Murray, Snežana Stanimirović, Carl Heiles, John M. Dickey, N. M. McClure-Griffiths, M. -Y. Lee, W. M. Goss, Nicholas Killerby-Smith

Abstract: Tracing the transition between the diffuse atomic interstellar medium (ISM) and cold, dense gas is crucial for deciphering the star formation cycle in galaxies. Here we present MACH, a new survey of cold neutral hydrogen (HI) absorption at $21\rm\,cm$ by the Karl G. Jansky Very Large Array. We target 42 bright background sources with $60<l<110^{\circ}$, $30<b<62^{\circ}$, significantly expanding t… ▽ More Tracing the transition between the diffuse atomic interstellar medium (ISM) and cold, dense gas is crucial for deciphering the star formation cycle in galaxies. Here we present MACH, a new survey of cold neutral hydrogen (HI) absorption at $21\rm\,cm$ by the Karl G. Jansky Very Large Array. We target 42 bright background sources with $60<l<110^{\circ}$, $30<b<62^{\circ}$, significantly expanding the sample of publicly-available, sensitive $21\rm\,cm$ absorption outside the Galactic plane. With matching $21\rm\,cm$ emission data from the EBHIS survey, we measure the total column density and cold HI fraction, and quantify the properties of individual HI structures along each sightline via autonomous Gaussian decomposition. Combining the MACH sample with results from recent HI absorption surveys, we produce a robust characterization of the cool atomic medium at high and intermediate Galactic latitudes. We find that MACH HI has significantly smaller column density relative to samples at similar latitudes, and the detected cold HI structures have smaller line widths, temperatures and turbulent Mach numbers, suggesting that MACH probes a particularly quiescent region. Using all available observations, we compute the cumulative covering fraction ($c$) of cold HI at local velocities outside the disk: structures with $τ>0.001$ are ubiquitous ($c\sim100\%$), whereas high optical depths ($τ>1$) are extremely rare ($c\sim0\%$). △ Less

Submitted 29 June, 2021; originally announced June 2021.

Comments: 33 pages, 16 figures. Accepted for publication in ApJS

Tiny-scale Structure Discovered toward PSR B1557$-$50

Authors: Mengting Liu, Marko Krčo, Di Li, George Hobbs, J. R. Dawson, Carl Heiles, Andrew Jameson, Snežana Stanimirović, Simon Johnston, John M. Dickey

Abstract: Optical depth variations in the Galactic neutral interstellar medium (ISM) with spatial scales from hundreds to thousands of astronomical units have been observed through HI absorption against pulsars and continuum sources, while extremely small structures with spatial scales of tens of astronomical units remain largely unexplored. The nature and formation of such tiny-scale atomic structures (TSA… ▽ More Optical depth variations in the Galactic neutral interstellar medium (ISM) with spatial scales from hundreds to thousands of astronomical units have been observed through HI absorption against pulsars and continuum sources, while extremely small structures with spatial scales of tens of astronomical units remain largely unexplored. The nature and formation of such tiny-scale atomic structures (TSAS) need to be better understood. We report a tentative detection of TSAS with a signal-to-noise ratio of 3.2 toward PSR B1557$-$50 in the second epoch of two Parkes sessions just 0.36 yr apart, which are the closest-spaced spectral observations toward this pulsar. One absorption component showing marginal variations has been identified. Based on the pulsar's proper motion of 14 mas $\rm yr^{-1}$ and the component's kinematic distance of 3.3 kpc, the corresponding characteristic spatial scale is 17 au, which is among the smallest sizes of known TSAS. Assuming a similar line-of-sight (LOS) depth, the tentative TSAS cloud detected here is overdense and overpressured relative to the cold neutral medium (CNM), and can radiatively cool fast enough to be in thermal equilibrium with the ambient environment. We find that turbulence is not sufficient to confine the overpressured TSAS. We explore the LOS elongation that would be required for the tentative TSAS to be at the canonical CNM pressure, and find that it is $\sim5000$ -- much larger than filaments observed in the ISM. We see some evidence of line width and temperature variations in the CNM components observed at the two epochs, as predicted by models of TSAS-like cloud formation colliding warm neutral medium flows. △ Less

Submitted 16 April, 2021; originally announced April 2021.

Comments: 11 pages, 2 figures, 2 tables, accepted for publication in the ApJL

Journal ref: The Astrophysical Journal Letters, 2021, Volume 911, Issue 1, id.L13, 9 pp

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

Revealing Gravitational Collapse in Serpens G3-G6 Molecular Cloud using Velocity Gradients

Authors: Yue Hu, A. Lazarian, Snezana Stanimirovic

Abstract: The relative role of turbulence, magnetic fields, self-gravity in star formation is a subject of intensive debate. We present IRAM 30m telescope observations of the $^{13}$CO (1-0) emission in the Serpens G3-G6 molecular cloud and apply to the data a set of statistical methods. Those include the probability density functions (PDFs) of column density and the Velocity Gradients Technique (VGT). We c… ▽ More The relative role of turbulence, magnetic fields, self-gravity in star formation is a subject of intensive debate. We present IRAM 30m telescope observations of the $^{13}$CO (1-0) emission in the Serpens G3-G6 molecular cloud and apply to the data a set of statistical methods. Those include the probability density functions (PDFs) of column density and the Velocity Gradients Technique (VGT). We combine our data with the Planck 353 GHz polarized dust emission observations, Hershel H$_2$ column density. We suggest that the Serpens G3-G6 south clump is undergoing a gravitational collapse. Our analysis reveals that the gravitational collapse happens at volume density $n\ge10^3$ $\rm cm^{-3}$. We estimate the plane-of-the-sky magnetic field strength of approximately 120 $μG$ using the traditional Davis-Chandrasekhar-Fermi method and 100 $μG$ using a new technique proposed in Lazarian et al.(2020). We find the Serpens G3-G6 south clump's total magnetic field energy significantly surpasses kinetic energy and gravitational energy. We conclude that the gravitational collapse could be successfully triggered in a supersonic and sub-Alfvénic cloud. △ Less

Submitted 9 March, 2021; v1 submitted 11 February, 2021; originally announced February 2021.

Comments: 17 pages,11 figures, accepted for publication in ApJ

Herschel 158$ μ$m [CII] Observations of "CO-dark" Gas in the Perseus Giant Molecular Cloud

Authors: K. P. Hall, Snežana Stanimirović, Min-Young Lee, Mark Wolfire, Paul Goldsmith

Abstract: We present observations of velocity-resolved [CII] 158 $μ$m emission from both a dense and a more diffuse photodissociation region (PDR) in the Perseus giant molecular cloud (GMC) using the Heterodyne Instrument for the Far-Infrared onboard the Herschel Space Telescope. We detect [CII] emission from 80\% of the total positions, with a 95\% detection rate from the dense boundary region. The integra… ▽ More We present observations of velocity-resolved [CII] 158 $μ$m emission from both a dense and a more diffuse photodissociation region (PDR) in the Perseus giant molecular cloud (GMC) using the Heterodyne Instrument for the Far-Infrared onboard the Herschel Space Telescope. We detect [CII] emission from 80\% of the total positions, with a 95\% detection rate from the dense boundary region. The integrated intensity of the [CII] emission remains relatively constant across each boundary, despite the observed range in optical extinction between 1 mag and 10 mag. This flat profile indicates a constant heating and cooling rate within both regions observed. The integrated intensity of [CII] emission is reasonably well correlated with the neutral hydrogen (HI) column density, as well as total gas column density. This, in addition to the 80$'$ (7 pc) extent of the [CII] emission from cloud center, suggests that the HI envelope plays a dominant role in explaining the [CII] emission emanating from Perseus. We compare the [CII] and $^{12}$CO integrated intensities with predictions from a 1-D, two-sided slab PDR model and show that a simple core $+$ envelope, equilibrium model without an additional "CO-dark" H$_2$ component can reproduce observations well. Additional observations are needed to disentangle how much of the [CII] emission is associated with the "CO-dark" H$_2$ gas, as well as constrain spatial variations of the dust-to-gas ratio across Perseus. △ Less

Submitted 17 June, 2020; originally announced June 2020.

Small-scale structure traced by neutral hydrogen absorption in the direction of multiple-component radio continuum sources

Authors: Daniel R. Rybarczyk, Snezana Stanimirovic, Ellen G. Zweibel, Claire E. Murray, John M. Dickey, Brian Babler, Carl Heiles

Abstract: We have studied the small scale distribution of atomic hydrogen (HI) using 21-cm absorption spectra against multiple-component background radio continuum sources from the 21-SPONGE survey and the Millennium Arecibo Absorption Line Survey. We have found $>5σ$ optical depth variations at a level of $\sim0.03-0.5$ between 13 out of 14 adjacent sightlines separated by a few arcseconds to a few arcminu… ▽ More We have studied the small scale distribution of atomic hydrogen (HI) using 21-cm absorption spectra against multiple-component background radio continuum sources from the 21-SPONGE survey and the Millennium Arecibo Absorption Line Survey. We have found $>5σ$ optical depth variations at a level of $\sim0.03-0.5$ between 13 out of 14 adjacent sightlines separated by a few arcseconds to a few arcminutes, suggesting the presence of neutral structures on spatial scales from a few to thousands of AU (which we refer to as tiny scale atomic structure, TSAS). The optical depth variations are strongest in directions where the HI column density and the fraction of HI in the cold neutral medium (CNM) are highest, which tend to be at low Galactic latitudes. By measuring changes in the properties of Gaussian components fitted to the absorption spectra, we find that changes in both the peak optical depth and the linewidth of TSAS absorption features contribute to the observed optical depth variations, while changes in the central velocity do not appear to strongly impact the observed variations. Both thermal and turbulent motions contribute appreciably to the linewidths, but the turbulence does not appear strong enough to confine overpressured TSAS. In a majority of cases, the TSAS column densities are sufficiently high that these structures can radiatively cool fast enough to maintain thermal equilibrium with their surroundings, even if they are overpressured. We also find that a majority of TSAS is associated with the CNM. For TSAS in the direction of the Taurus molecular cloud and the local Leo cold cloud, we estimate densities over an order of magnitude higher than typical CNM densities. △ Less

Submitted 27 April, 2020; v1 submitted 24 February, 2020; originally announced February 2020.

Comments: 29 pages, 19 figures, 4 tables. Accepted to ApJ. Section 4.3.1 updated from previous version

Mapping Spatial Variations of HI Turbulent Properties in the Small and Large Magellanic Cloud

Authors: Samuel Szotkowski, Delano Yoder, Snežana Stanimirović, Brian Babler, N. M. McClure-Griffiths, Helga Dénes, Alberto Bolatto, Katherine Jameson, Lister Staveley-Smith

Abstract: We developed methods for mapping spatial variations of the spatial power spectrum (SPS) and structure function (SF) slopes, with a goal of connecting neutral hydrogen (HI) statistical properties with the turbulent drivers. The new methods were applied on the HI observations of the Small and Large Magellanic Clouds (SMC and LMC). In the case of the SMC, we find highly uniform turbulent properties o… ▽ More We developed methods for mapping spatial variations of the spatial power spectrum (SPS) and structure function (SF) slopes, with a goal of connecting neutral hydrogen (HI) statistical properties with the turbulent drivers. The new methods were applied on the HI observations of the Small and Large Magellanic Clouds (SMC and LMC). In the case of the SMC, we find highly uniform turbulent properties of HI, with no evidence for local enhancements of turbulence due to stellar feedback. Such properties could be caused by a significant turbulent driving on large-scales. Alternatively, a significant line-of-sight depth of the SMC could be masking out localized regions with a steeper SPS slope caused by stellar feedback. In contrast to the SMC, the LMC HI shows a large diversity in terms of its turbulent properties. Across most of the LMC, the small-scale SPS slope is steeper than the large-scale slope due to the presence of the HI disk. On small spatial scales, we find several areas of localized steepening of the SPS slope around major HII regions, with the 30 Doradus region being the most prominent. This is in agreement with predictions from numerical simulations which suggest steepening of the SPS slope due to stellar feedback eroding and destroying interstellar clouds. We also find localized steepening of the large-scale SPS slope in the outskirts of the LMC. This is likely caused by the flaring of the HI disk, or alternatively ram-pressure stripping of the LMC disk due to the interactions with the surrounding halo gas. △ Less

Submitted 11 November, 2019; originally announced November 2019.

Comments: 17 pages, 14 figures, accepted to the ApJ

Molecular gas in the outflow of the Small Magellanic Cloud

Authors: Enrico M. Di Teodoro, N. M. McClure-Griffiths, C. De Breuck, L. Armillotta, N. M. Pingel, K. E. Jameson, J. M. Dickey, M. Rubio, S. Stanimirovic, L. Staveley-Smith

Abstract: We report the first evidence of molecular gas in two atomic hydrogen (HI) clouds associated with gas outflowing from the Small Magellanic Cloud (SMC). We used the Atacama Pathfinder Experiment (APEX) to detect and spatially resolve individual clumps of CO(2-1) emission in both clouds. CO clumps are compact (~ 10 pc) and dynamically cold (linewidths < 1 km/s). Most CO emission appears to be offset… ▽ More We report the first evidence of molecular gas in two atomic hydrogen (HI) clouds associated with gas outflowing from the Small Magellanic Cloud (SMC). We used the Atacama Pathfinder Experiment (APEX) to detect and spatially resolve individual clumps of CO(2-1) emission in both clouds. CO clumps are compact (~ 10 pc) and dynamically cold (linewidths < 1 km/s). Most CO emission appears to be offset from the peaks of the HI emission, some molecular gas lies in regions without a clear HI counterpart. We estimate a total molecular gas mass of 10^3-10^4 Msun in each cloud and molecular gas fractions up to 30% of the total cold gas mass (molecular + neutral). Under the assumption that this gas is escaping the galaxy, we calculated a cold gas outflow rate of 0.3-1.8 Msun/yr and mass loading factors of 3 -12 at a distance larger than 1 kpc. These results show that relatively weak star-formation-driven winds in dwarf galaxies like the SMC are able to accelerate significant amounts of cold and dense matter and inject it into the surrounding environment. △ Less

Submitted 23 October, 2019; originally announced October 2019.

Comments: Accepted for publication in ApJ Letters

An ATCA Survey of HI Absorption in the Magellanic Clouds I: HI Gas Temperature Measurements in the Small Magellanic Cloud

Authors: Katherine Jameson, Naomi McClure-Griffiths, Boyang Liu, John Dickey, Lister Staveley-Smith, Snezana Stanimirovic, James Dempsey, Joanne Dawson, Helga Denes, Alberto Bolatto, Di Li, Tony Wong

Abstract: We present the first results from the Small Magellanic Cloud portion of a new Australia Telescope Compact Array (ATCA) HI absorption survey of both of the Magellanic Clouds, comprising over 800 hours of observations. Our new HI absorption line data allow us to measure the temperature and fraction of cold neutral gas in a low metallicity environment. We observed 22 separate fields, targeting a tota… ▽ More We present the first results from the Small Magellanic Cloud portion of a new Australia Telescope Compact Array (ATCA) HI absorption survey of both of the Magellanic Clouds, comprising over 800 hours of observations. Our new HI absorption line data allow us to measure the temperature and fraction of cold neutral gas in a low metallicity environment. We observed 22 separate fields, targeting a total of 55 continuum sources against 37 of which we detected HI absorption; from this we measure a column density weighted mean average spin temperature of $<T_{s}>=150$ K. Splitting the spectra into individual absorption line features, we estimate the temperatures of different gas components and find an average cold gas temperature of $\sim{30}$ K for this sample, lower than the average of $\sim{40}$ K in the Milky Way. The HI appears to be evenly distributed throughout the SMC and we detect absorption in $67\%$ of the lines of sight in our sample, including some outside the main body of the galaxy ($N_{\text{HI}}>2\times{10^{21}}$ cm$^{-2}$). The optical depth and temperature of the cold neutral atomic gas shows no strong trend with location spatially or in velocity. Despite the low metallicity environment, we find an average cold gas fraction of $\sim{20\%}$, not dissimilar from that of the Milky Way. △ Less

Submitted 13 August, 2019; originally announced August 2019.

Comments: 22 pages, 16 figures, accepted for publication in ApJS

Exploring the properties of warm and cold atomic hydrogen in the Taurus and Gemini regions

Authors: Hiep Nguyen, J. R. Dawson, Min-Young Lee, Claire E. Murray, Snezana Stanimirovic, Carl Heiles, M. -A. Miville-Deschenes, Anita Petzler

Abstract: We report Arecibo 21 cm absorption-emission observations to characterise the physical properties of neutral hydrogen (HI) in the proximity of five giant molecular clouds (GMCs): Taurus, California, Rosette, Mon OB1, and NGC 2264. Strong HI absorption was detected toward all 79 background continuum sources in the ~60x20 square degree region. Gaussian decompositions were performed to estimate temper… ▽ More We report Arecibo 21 cm absorption-emission observations to characterise the physical properties of neutral hydrogen (HI) in the proximity of five giant molecular clouds (GMCs): Taurus, California, Rosette, Mon OB1, and NGC 2264. Strong HI absorption was detected toward all 79 background continuum sources in the ~60x20 square degree region. Gaussian decompositions were performed to estimate temperatures, optical depths and column densities of the cold and warm neutral medium (CNM, WNM). The properties of individual CNM components are similar to those previously observed along random Galactic sightlines and in the vicinity of GMCs, suggesting a universality of cold HI properties. The CNM spin temperature (Ts) histogram peaks at ~50K. The turbulent Mach numbers of CNM vary widely, with a typical value of ~4, indicating that their motions are supersonic. About 60% of the total HI gas is WNM, and nearly 40% of the WNM lies in thermally unstable regime 500-5000K. The observed CNM fraction is higher around GMCs than in diffuse regions, and increases with increasing column density (NHI) to a maximum of ~75%. On average, the optically thin approximation (N*(HI)) underestimates the total N(HI) by ~21%, but we find large regional differences in the relationship between N(HI) and the required correction factor, f=N(HI)/N*(HI). We examine two different methods (linear fit of f vs log10(N*(HI)) and uniform Ts) to correct for opacity effects using emission data from the GALFA-HI survey. We prefer the uniform Ts method, since the linear relationship does not produce convincing fits for all subregions. △ Less

Submitted 17 June, 2019; originally announced June 2019.

Comments: This article (22 pages, 23 figures) has been accepted for publication in The Astrophysical Journal (ApJ)

The Magellanic Stream as a Probe of Astrophysics

Authors: Andrew J. Fox, Kathleen A. Barger, Joss Bland-Hawthorn, Dana Casetti-Dinescu, Elena D'Onghia, Felix J. Lockman, Naomi McClure-Griffiths, David Nidever, Mary Putman, Philipp Richter, Snezana Stanimirovic, Thorsten Tepper-Garcia

Abstract: Extending for over 200 degrees across the sky, the Magellanic Stream together with its Leading Arm is the most spectacular example of a gaseous stream in the local Universe. The Stream is an interwoven tail of filaments trailing the Magellanic Clouds as they orbit the Milky Way. Thought to be created by tidal forces, ram pressure, and halo interactions, the Stream is a benchmark for dynamical mode… ▽ More Extending for over 200 degrees across the sky, the Magellanic Stream together with its Leading Arm is the most spectacular example of a gaseous stream in the local Universe. The Stream is an interwoven tail of filaments trailing the Magellanic Clouds as they orbit the Milky Way. Thought to be created by tidal forces, ram pressure, and halo interactions, the Stream is a benchmark for dynamical models of the Magellanic System, a case study for gas accretion and dwarf-galaxy accretion onto galaxies, a probe of the outer halo, and the bearer of more gas mass than all other Galactic high velocity clouds combined. If it survives to reach the Galactic disk, it may maintain or even elevate the Galactic star-formation rate. In this white paper, we emphasize the Stream's importance for many areas of Galactic astronomy, summarize key unanswered questions, and identify future observations and simulations needed to resolve them. We stress the importance of ultraviolet, optical, and radio spectroscopy, and the need for computational models that capture full particle and radiation treatments within an MHD environment. △ Less

Submitted 12 March, 2019; originally announced March 2019.

Comments: Science White Paper submitted to Astro2020 Decadal Survey. 7 pages, 2 figures

On the dynamics of the Small Magellanic Cloud through high-resolution ASKAP HI observations

Authors: E. M. Di Teodoro, N. M. McClure-Griffiths, K. E. Jameson, H. Denes, John M. Dickey, S. Stanimirovic, L. Staveley-Smith, C. Anderson, J. D. Bunton, A. Chippendale, K. Lee-Waddell, A. MacLeod, M. A Voronkov

Abstract: We use new high-resolution HI data from the Australian Square Kilometre Array Pathfinder (ASKAP) to investigate the dynamics of the Small Magellanic Cloud (SMC). We model the HI gas component as a rotating disc of non-negligible angular size, moving into the plane of the sky and undergoing nutation/precession motions. We derive a high-resolution (~ 10 pc) rotation curve of the SMC out to R ~ 4 kpc… ▽ More We use new high-resolution HI data from the Australian Square Kilometre Array Pathfinder (ASKAP) to investigate the dynamics of the Small Magellanic Cloud (SMC). We model the HI gas component as a rotating disc of non-negligible angular size, moving into the plane of the sky and undergoing nutation/precession motions. We derive a high-resolution (~ 10 pc) rotation curve of the SMC out to R ~ 4 kpc. After correcting for asymmetric drift, the circular velocity slowly rises to a maximum value of Vc ~ 55 km/s at R ~ 2.8 kpc and possibly flattens outwards. In spite of the SMC undergoing strong gravitational interactions with its neighbours, its HI rotation curve is akin to that of many isolated gas-rich dwarf galaxies. We decompose the rotation curve and explore different dynamical models to deal with the unknown three-dimensional shape of the mass components (gas, stars and dark matter). We find that, for reasonable mass-to-light ratios, a dominant dark matter halo with mass M(R<4 kpc) = 1-1.5 x 10^9 solar masses is always required to successfully reproduce the observed rotation curve, implying a large baryon fraction of 30%-40%. We discuss the impact of our assumptions and the limitations of deriving the SMC kinematics and dynamics from HI observations. △ Less

Submitted 23 November, 2018; originally announced November 2018.

Comments: Accepted for publication in MNRAS

Cold gas outflows from the Small Magellanic Cloud traced with ASKAP

Authors: N. M. McClure-Griffiths, H. Dénes, J. M. Dickey, S. Stanimirović, L. Staveley-Smith, Katherine Jameson, Enrico Di Teodoro, James R. Allison, J. D. Collier, A. P. Chippendale, T. Franzen, Gülay Gürkan, G. Heald, A. Hotan, D. Kleiner, K. Lee-Waddell, D. McConnell, A. Popping, Jonghwan Rhee, C. J. Riseley, M. A. Voronkov, M. Whiting

Abstract: Feedback from massive stars plays a critical role in the evolution of the Universe by driving powerful outflows from galaxies that enrich the intergalactic medium and regulate star formation. An important source of outflows may be the most numerous galaxies in the Universe: dwarf galaxies. With small gravitational potential wells, these galaxies easily lose their star-forming material in the prese… ▽ More Feedback from massive stars plays a critical role in the evolution of the Universe by driving powerful outflows from galaxies that enrich the intergalactic medium and regulate star formation. An important source of outflows may be the most numerous galaxies in the Universe: dwarf galaxies. With small gravitational potential wells, these galaxies easily lose their star-forming material in the presence of intense stellar feedback. Here, we show that the nearby dwarf galaxy, the Small Magellanic Cloud (SMC), has atomic hydrogen outflows extending at least 2 kiloparsecs (kpc) from the star-forming bar of the galaxy. The outflows are cold, $T<400~{\rm K}$, and may have formed during a period of active star formation $25 - 60$ million years (Myr) ago. The total mass of atomic gas in the outflow is $\sim 10^7$ solar masses, ${\rm M_{\odot}}$, or $\sim 3$% of the total atomic gas of the galaxy. The inferred mass flux in atomic gas alone, $\dot{M}_{HI}\sim 0.2 - 1.0~{\rm M_{\odot}~yr^{-1}}$, is up to an order of magnitude greater than the star formation rate. We suggest that most of the observed outflow will be stripped from the SMC through its interaction with its companion, the Large Magellanic Cloud (LMC), and the Milky Way, feeding the Magellanic Stream of hydrogen encircling the Milky Way. △ Less

Submitted 5 November, 2018; originally announced November 2018.

Comments: Published in Nature Astronomy, 29 October 2018, http://dx.doi.org/10.1038/s41550-018-0608-8

Journal ref: Nature Astronomy 2, pg 901 - 906 (2018)

Science with an ngVLA: Neutral Atomic Hydrogen in the Local Universe

Authors: D. J. Pisano, F. Walter, S. Stanimirovic

Abstract: One of the outstanding questions in astronomy today is how gas flows from the circumgalactic medium (CGM) onto the disks of galaxies and then transitions from the diffuse atomic medium into molecular star-forming cores. For studies of the CGM, the Next Generation Very Large Array (ngVLA) will have the sensitivity and resolution to measure the sizes of the neutral atomic hydrogen (HI) disks of gala… ▽ More One of the outstanding questions in astronomy today is how gas flows from the circumgalactic medium (CGM) onto the disks of galaxies and then transitions from the diffuse atomic medium into molecular star-forming cores. For studies of the CGM, the Next Generation Very Large Array (ngVLA) will have the sensitivity and resolution to measure the sizes of the neutral atomic hydrogen (HI) disks of galaxies and complete a census of the HI content around galaxies. Within galaxies, the ngVLA will be able to resolve HI clouds in large numbers of galaxies beyond the Local Group providing measurements of the physical conditions of gas across a wide range of galaxy types. Finally, within our own Milky Way, the ngVLA will provide a dense grid of HI absorption spectra in the cold and warm neutral medium constraining the temperature and density of atomic gas as it transitions into molecular gas. Combined with radio continuum and molecular line data from the ngVLA plus multi-wavelength data from other planned facilities, ngVLA will have a key role in understanding star-formation in the local universe while complementing future studies with the Square Kilometer Array. △ Less

Submitted 15 October, 2018; originally announced October 2018.

Comments: To be published in the ASP Monograph Series, Science with a Next-Generation VLA, ed. E. J. Murphy (ASP, San Francisco, CA); 12 pages, 3 figures

Atomic and Ionized Microstructures in the Diffuse Interstellar Medium

Authors: Snezana Stanimirovic, Ellen G. Zweibel

Abstract: It has been known for half a century that the interstellar medium (ISM) of our Galaxy is structured on scales as small as a few hundred km, more than 10 orders of magnitude smaller than typical ISM structures and energy input scales. In this review we focus on neutral and ionized structures on spatial scales of a few to ~10^4 Astronomical Units (AU) which appear to be highly overpressured, as thes… ▽ More It has been known for half a century that the interstellar medium (ISM) of our Galaxy is structured on scales as small as a few hundred km, more than 10 orders of magnitude smaller than typical ISM structures and energy input scales. In this review we focus on neutral and ionized structures on spatial scales of a few to ~10^4 Astronomical Units (AU) which appear to be highly overpressured, as these have the most important role in the dynamics and energy balance of interstellar gas: the Tiny Scale Atomic Structure (TSAS) and Extreme Scattering Events (ESEs) as the most over-pressured example of the Tiny Scale Ionized Structures (TSIS). We review observational results and highlight key physical processes at AU scales. We present evidence for and against microstructures as part of a universal turbulent cascade and as discrete structures, and review their association with supernova remnants, the Local Bubble, and bright stars. We suggest a number of observational and theoretical programs that could clarify the nature of AU structures. TSAS and TSIS probe spatial scales in the range of what is expected for turbulent dissipation scales, therefore are of key importance for constraining exotic and not-well understood physical processes which have implications for many areas of astrophysics. The emerging picture is one in which a magnetized, turbulent cascade, driven hard by a local energy source and acting jointly with phenomena such as thermal instability, is the source of these microstructures. △ Less

Submitted 1 October, 2018; originally announced October 2018.

Comments: Paper published in 2018 Annual Review of Astronomy & Astrophysics. 51 pages. Please use the final PDF file: http://www.annualreviews.org/eprint/JsAxq3hEFMUHm6NEBdum/full/10.1146/annurev-astro-081817-051810

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

The 21-SPONGE HI Absorption Line Survey II: The temperature of Galactic HI

Authors: Claire E. Murray, Snezana Stanimirovic, W. M. Goss, Carl Heiles, John M. Dickey, Brian Babler, Chang-Goo Kim

Abstract: We present 21-cm Spectral Line Observations of Neutral Gas with the VLA (21-SPONGE), a Karl G. Jansky Very Large Array (VLA) large project (~600 hours) for measuring the physical properties of Galactic neutral hydrogen (HI). 21-SPONGE is distinguished among previous Galactic HI studies as a result of: (1) exceptional optical depth sensitivity ($σ_τ < 10^{-3}$ per $0.42\rm\,km\,s^{-1}$ channels ove… ▽ More We present 21-cm Spectral Line Observations of Neutral Gas with the VLA (21-SPONGE), a Karl G. Jansky Very Large Array (VLA) large project (~600 hours) for measuring the physical properties of Galactic neutral hydrogen (HI). 21-SPONGE is distinguished among previous Galactic HI studies as a result of: (1) exceptional optical depth sensitivity ($σ_τ < 10^{-3}$ per $0.42\rm\,km\,s^{-1}$ channels over 57 lines of sight); (2) matching 21 cm emission spectra with highest-possible angular resolution (~4') from the Arecibo Observatory; (3) detailed comparisons with numerical simulations for assessing observational biases. We autonomously decompose 21 cm spectra and derive the physical properties (i.e., spin temperature, $T_s$, column density) of the cold neutral medium (CNM; $T_s<250\rm\,K$), thermally unstable medium (UNM; $250< T_s < 1000\rm\,K$) and warm neutral medium (WNM; $T_s > 1000\rm\,K$) simultaneously. We detect 50% of the total HI mass in absorption, the majority of which is CNM (56 +/- 10%, corresponding to 28% of the total HI mass). Although CNM is detected ubiquitously, the CNM fraction along most lines of sight is <50%. We find that 20% of the total HI mass is thermally unstable (41 +/- 10% of HI detected in absorption), with no significant variation with Galactic environment. Finally, although the WNM comprises 52% of the total HI mass, we detect little evidence for WNM absorption with $1000<T_s<4000\rm\,K$. Following spectral modeling, we detect a stacked residual absorption feature corresponding to WNM with $T_s\sim10^4\rm\,K$. We conclude that excitation in excess of collisions likely produces significantly higher WNM $T_s$ than predicted by steady-state models. △ Less

Submitted 15 June, 2018; originally announced June 2018.

Comments: 33 pages, 22 figures, 5 tables; ApJ submitted

Optically thick HI does not dominate dark gas in the local ISM

Authors: Claire E. Murray, J. E. G. Peek, Min-Young Lee, Snezana Stanimirovic

Abstract: The local interstellar medium (ISM) is suffused with "dark" gas, identified by excess infrared and gamma ray emission, yet undetected by standard ISM tracers such as neutral hydrogen (HI) or carbon monoxide emission. Based on observed dust properties from Planck, recent studies have argued that HI mixed with dust is strongly saturated and that dark gas is dominated by optically-thick HI. We test t… ▽ More The local interstellar medium (ISM) is suffused with "dark" gas, identified by excess infrared and gamma ray emission, yet undetected by standard ISM tracers such as neutral hydrogen (HI) or carbon monoxide emission. Based on observed dust properties from Planck, recent studies have argued that HI mixed with dust is strongly saturated and that dark gas is dominated by optically-thick HI. We test this hypothesis by reproducing this model using data from Planck and new 21 cm emission maps from GALFA-HI -- the first large-area 21cm emission survey with comparable angular resolution to Planck. We compare the results with those from a large sample of HI column densities based on direct observations of HI optical depth, and find that the inferred column density corrections are significantly lower than those inferred by the Planck-based model. Further, we rule out the hypothesis that the pencil-beam HI absorption sight lines preferentially miss opaque "blobs" with small covering fraction, as these structures require densities and pressures which are incompatible with ISM conditions. Our results support the picture that excess dust emission in the local ISM is not dominated by optically-thick HI, but is rather a combination of intrinsic changes in dust grain emissivities and H2 missed by CO observations. △ Less

Submitted 14 June, 2018; v1 submitted 4 June, 2018; originally announced June 2018.

Comments: 10 pages, 7 figures; ApJ accepted

Dust-Gas Scaling Relations and OH Abundance in the Galactic ISM

Authors: Hiep Nguyen, J. R. Dawson, M. -A. Miville-Deschênes, Ningyu Tang, Di Li, Carl Heiles, Claire E. Murray, Snežana Stanimirović, Steven J. Gibson, N. M. McClure-Griffiths, Thomas Troland, L. Bronfman, R. Finger

Abstract: Observations of interstellar dust are often used as a proxy for total gas column density $N_\mathrm{H}$. By comparing $\textit{Planck}$ thermal dust data (Release 1.2) and new dust reddening maps from Pan-STARRS 1 and 2MASS (Green et al. 2018), with accurate (opacity-corrected) HI column densities and newly-published OH data from the Arecibo Millennium survey and 21-SPONGE, we confirm linear corre… ▽ More Observations of interstellar dust are often used as a proxy for total gas column density $N_\mathrm{H}$. By comparing $\textit{Planck}$ thermal dust data (Release 1.2) and new dust reddening maps from Pan-STARRS 1 and 2MASS (Green et al. 2018), with accurate (opacity-corrected) HI column densities and newly-published OH data from the Arecibo Millennium survey and 21-SPONGE, we confirm linear correlations between dust optical depth $τ_{353}$, reddening $E(B{-}V)$ and the total proton column density $N_\mathrm{H}$ in the range (1$-$30)$\times$10$^{20}$cm$^{-2}$, along sightlines with no molecular gas detections in emission. We derive an $N_\mathrm{H}$/$E(B{-}V)$ ratio of (9.4$\pm$1.6)$\times$10$^{21}$cm$^{-2}$mag$^{-1}$ for purely atomic sightlines at $|b|$$>$5$^{\circ}$, which is 60$\%$ higher than the canonical value of Bohlin et al. (1978). We report a $\sim$40$\%$ increase in opacity $σ_{353}$=$τ_{353}$/$N_\mathrm{H}$, when moving from the low column density ($N_\mathrm{H}$$<$5$\times$10$^{20}$cm$^{-2}$) to moderate column density ($N_\mathrm{H}$$>$5$\times$10$^{20}$cm$^{-2}$) regime, and suggest that this rise is due to the evolution of dust grains in the atomic ISM. Failure to account for HI opacity can cause an additional apparent rise in $σ_{353}$, of the order of a further $\sim$20$\%$. We estimate molecular hydrogen column densities $N_{\mathrm{H}_{2}}$ from our derived linear relations, and hence derive the OH/H$_2$ abundance ratio of $X_\mathrm{OH}$$\sim$1$\times$10$^{-7}$ for all molecular sightlines. Our results show no evidence of systematic trends in OH abundance with $N_{\mathrm{H}_{2}}$ in the range $N_{\mathrm{H}_{2}}$$\sim$(0.1$-$10)$\times$10$^{21}$cm$^{-2}$. This suggests that OH may be used as a reliable proxy for H$_2$ in this range, which includes sightlines with both CO-dark and CO-bright gas. △ Less

Submitted 29 May, 2018; originally announced May 2018.

Comments: The revised manuscript is accepted for publication in The Astrophysical Journal

The Multi-phase Turbulence Density Power Spectra in the Perseus Molecular Cloud

Authors: Nickolas Pingel, Min-Young Lee, Blakesley Burkhart, Snezana Stanimirović

Abstract: We derive two-dimensional spatial power spectra of four distinct interstellar medium tracers, HI, $^{12}$CO($J$=1--0), $^{13}$CO($J$=1--0), and dust, in the Perseus molecular cloud, covering linear scales ranging from $\sim$0.1 pc to $\sim$90 pc. Among the four tracers, we find the steepest slopes of $-3.23\pm0.05$ and $-3.22\pm0.05$ for the uncorrected and opacity-corrected HI column density imag… ▽ More We derive two-dimensional spatial power spectra of four distinct interstellar medium tracers, HI, $^{12}$CO($J$=1--0), $^{13}$CO($J$=1--0), and dust, in the Perseus molecular cloud, covering linear scales ranging from $\sim$0.1 pc to $\sim$90 pc. Among the four tracers, we find the steepest slopes of $-3.23\pm0.05$ and $-3.22\pm0.05$ for the uncorrected and opacity-corrected HI column density images. This result suggests that the HI in and around Perseus traces a non-gravitating, transonic medium on average, with a negligible effect from opacity. On the other hand, we measure the shallowest slope of $-2.72\pm0.12$ for the 2MASS dust extinction data and interpret this as the signature of a self-gravitating, supersonic medium. Possible variations in the dust-to-gas ratio likely do not change our conclusion. Finally, we derive slopes of $-3.08\pm0.08$ and $-2.88\pm0.07$ for the ${}^{12}$CO(1--0) and ${}^{13}$CO(1--0) integrated intensity images. Based on theoretical predictions for an optically thick medium, we interpret these slopes of roughly $-3$ as implying that both CO lines are susceptible to the opacity effect. While simple tests for the impact of CO formation and depletion indicate that the measured slopes of ${}^{12}$CO(1--0) and ${}^{13}$CO(1--0) are not likely affected by these chemical effects, our results generally suggest that chemically more complex and/or fully optically thick media may not be a reliable observational tracer for characterizing turbulence. △ Less

Submitted 27 February, 2018; originally announced February 2018.

Comments: Accepted for publication in The Astrophysical Journal

Where is OH and Does It Trace the Dark Molecular Gas (DMG)?

Authors: Di Li, Ningyu Tang, Hiep Nguyen, J. R. Dawson, Carl Heiles, Duo Xu, Zhichen Pan, Paul F. Goldsmith, Steven J. Gibson, Claire E. Murray, Tim Robishaw, N. M. McClure-Griffiths, John Dickey, Jorge Pineda, Snežana Stanimirović, L. Bronfman, Thomas Troland, the PRIMO collaboration

Abstract: Hydroxyl (OH) is expected to be abundant in diffuse interstellar molecular gas as it forms along with $H_2$ under similar conditions and within a similar extinction range. We have analyzed absorption measurements of OH at 1665 MHz and 1667 MHz toward 44 extragalactic continuum sources, together with the J=1-0 transitions of $^{12}$CO, $^{13}$CO , and C$^{18}$O, and the J=2-1 of $^{12}$CO. The exci… ▽ More Hydroxyl (OH) is expected to be abundant in diffuse interstellar molecular gas as it forms along with $H_2$ under similar conditions and within a similar extinction range. We have analyzed absorption measurements of OH at 1665 MHz and 1667 MHz toward 44 extragalactic continuum sources, together with the J=1-0 transitions of $^{12}$CO, $^{13}$CO , and C$^{18}$O, and the J=2-1 of $^{12}$CO. The excitation temperature of OH were found to follow a modified log-normal distribution, $ f(T{\rm_{ex}}) \propto \frac{1}{ \sqrt{2π}σ} \rm{exp}\left[-\frac{[ln(\textit{T}_{ex})-ln(3.4\ K)]^2}{2σ^2}\right] $, the peak of which is close to the temperature of the Galactic emission background (CMB+synchron). In fact, 90% of the OH has excitation temperature within 2 K of the Galactic background at the same location, providing a plausible explanation for the apparent difficulty to map this abundant molecule in emission. The opacities of OH were found to be small and peak around 0.01. For gas at intermediate extinctions (A$\rm_V$ $\sim$ 0.05-2 mag), the detection rate of OH with detection limit $N(\mathrm{OH})\simeq 10^{12}$ cm$^{-2}$ is approximately independent of $A\rm_V$. We conclude that OH is abundant in the diffuse molecular gas and OH absorption is a good tracer of `dark molecular gas (DMG)'. The measured fraction of DMG depends on assumed detection threshold of the CO data set. The next generation of highly sensitive low frequency radio telescopes, FAST and SKA, will make feasible the systematic inventory of diffuse molecular gas, through decomposing in velocity the molecular (e.g. OH and CH) absorption profiles toward background continuum sources with numbers exceeding what is currently available by orders of magnitude. △ Less

Submitted 12 January, 2018; originally announced January 2018.

Comments: 24 pages, 23 figures. Accepted for publication in ApJS

A new approach for short-spacing correction of radio interferometric data sets

Authors: S. Faridani, F. Bigiel, L. Floeer, J. Kerp, S. Stanimirovic

Abstract: The short-spacing problem describes the inherent inability of radio-interferometric arrays to measure the integrated flux and structure of diffuse emission associated with extended sources. New interferometric arrays, such as SKA, require solutions to efficiently combine interferometer and single-dish data. We present a new and open source approach for merging single-dish and cleaned interferome… ▽ More The short-spacing problem describes the inherent inability of radio-interferometric arrays to measure the integrated flux and structure of diffuse emission associated with extended sources. New interferometric arrays, such as SKA, require solutions to efficiently combine interferometer and single-dish data. We present a new and open source approach for merging single-dish and cleaned interferometric data sets requiring a minimum of data manipulation while offering a rigid flux determination and full high angular resolution. Our approach combines single-dish and cleaned interferometric data in the image domain. This approach is tested for both Galactic and extragalactic HI data sets. Furthermore, a quantitative comparison of our results to commonly used methods is provided. Additionally, for the interferometric data sets of NGC4214 and NGC5055, we study the impact of different imaging parameters as well as their influence on the combination for NGC4214. The approach does not require the raw data (visibilities) or any additional special information such as antenna patterns. This is advantageous especially in the light of upcoming radio surveys with heterogeneous antenna designs. △ Less

Submitted 27 September, 2017; originally announced September 2017.

Comments: 15 pages, accepted for publication by Astronomical Notes

An Overview of the MHONGOOSE Survey: Observing Nearby Galaxies with MeerKAT

Authors: W. J. G. de Blok, E. A. K. Adams, P. Amram, E. Athanassoula, I. Bagetakos, C. Balkowski, M. A. Bershady, R. Beswick, F. Bigiel, S. -L. Blyth, A. Bosma, R. S. Booth, A. Bouchard, E. Brinks, C. Carignan, L. Chemin, F. Combes, J. Conway, E. C. Elson, J. English, B. Epinat, B. S. Frank, J. Fiege, F. Fraternali, J. S. Gallagher , et al. (43 additional authors not shown)

Abstract: MHONGOOSE is a deep survey of the neutral hydrogen distribution in a representative sample of 30 nearby disk and dwarf galaxies with HI masses from 10^6 to ~10^{11} M_sun, and luminosities from M_R ~ -12 to M_R ~ -22. The sample is selected to uniformly cover the available range in log(M_HI). Our extremely deep observations, down to HI column density limits of well below 10^{18} cm^{-2} - or a few… ▽ More MHONGOOSE is a deep survey of the neutral hydrogen distribution in a representative sample of 30 nearby disk and dwarf galaxies with HI masses from 10^6 to ~10^{11} M_sun, and luminosities from M_R ~ -12 to M_R ~ -22. The sample is selected to uniformly cover the available range in log(M_HI). Our extremely deep observations, down to HI column density limits of well below 10^{18} cm^{-2} - or a few hundred times fainter than the typical HI disks in galaxies - will directly detect the effects of cold accretion from the intergalactic medium and the links with the cosmic web. These observations will be the first ever to probe the very low-column density neutral gas in galaxies at these high resolutions. Combination with data at other wavelengths, most of it already available, will enable accurate modelling of the properties and evolution of the mass components in these galaxies and link these with the effects of environment, dark matter distribution, and other fundamental properties such as halo mass and angular momentum. MHONGOOSE can already start addressing some of the SKA-1 science goals and will provide a comprehensive inventory of the processes driving the transformation and evolution of galaxies in the nearby universe at high resolution and over 5 orders of magnitude in column density. It will be a Nearby Galaxies Legacy Survey that will be unsurpassed until the advent of the SKA, and can serve as a highly visible, lasting statement of MeerKAT's capabilities. △ Less

Submitted 25 September, 2017; originally announced September 2017.

Comments: To be published in "MeerKAT Science: On the Pathway to the SKA". Proceedings of Science. Workshop held 25-27 May, 2016 Stellenbosch, South Africa

Spatial variations of turbulent properties of neutral hydrogen gas in the Small Magellanic Cloud using structure function analysis

Authors: David Nestingen-Palm, Snezana Stanimirovic, Diego F. Gonzalez-Casanova, Brian Babler, Katherine Jameson, Alberto Bolatto

Abstract: We investigate spatial variations of turbulent properties in the Small Magellanic Cloud (SMC) by using neutral hydrogen HI observations. With the goal of testing the importance of stellar feedback on HI turbulence, we define central and outer SMC regions based on the star formation rate (SFR) surface density, as well as the HI integrated intensity. We use the structure function and the Velocity Ch… ▽ More We investigate spatial variations of turbulent properties in the Small Magellanic Cloud (SMC) by using neutral hydrogen HI observations. With the goal of testing the importance of stellar feedback on HI turbulence, we define central and outer SMC regions based on the star formation rate (SFR) surface density, as well as the HI integrated intensity. We use the structure function and the Velocity Channel Analysis (VCA) to calculate the power-law index (gamma) for both underlying density and velocity fields in these regions. In all cases, our results show essentially no difference in gamma between the central and outer regions. This suggests that HI turbulent properties are surprisingly homogeneous across the SMC when probed at a resolution of 30 pc. Contrary to recent suggestions from numerical simulations, we do not find a significant change in gamma due to stellar feedback as traced by the SFR surface density. This could be due to the stellar feedback being widespread over the whole of the SMC, but more likely due to a large-scale gravitational driving of turbulence. We show that the lack of difference between central and outer SMC regions can not be explained by the high optical depth HI. △ Less

Submitted 10 July, 2017; originally announced July 2017.

Comments: 10 pages, 7 figures, accepted for publication in ApJ

Recovering interstellar gas properties with HI spectral lines: A comparison between synthetic spectra and 21-SPONGE

Authors: Claire E. Murray, Snezana Stanimirovic, Chang-Goo Kim, Eve C. Ostriker, Robert R. Lindner, Carl Heiles, John M. Dickey, Brian Babler

Abstract: We analyze synthetic neutral hydrogen (HI) absorption and emission spectral lines from a high- resolution, three-dimensional hydrodynamical simulation to quantify how well observational methods recover the physical properties of interstellar gas. We present a new method for uniformly decomposing HI spectral lines and estimating the properties of associated gas using the Autonomous Gaussian Decompo… ▽ More We analyze synthetic neutral hydrogen (HI) absorption and emission spectral lines from a high- resolution, three-dimensional hydrodynamical simulation to quantify how well observational methods recover the physical properties of interstellar gas. We present a new method for uniformly decomposing HI spectral lines and estimating the properties of associated gas using the Autonomous Gaussian Decomposition (AGD) algorithm. We find that Hi spectral lines recover physical structures in the simulation with excellent completeness at high Galactic latitude, and this completeness declines with decreasing latitude due to strong velocity-blending of spectral lines. The temperature and column density inferred from our decomposition and radiative transfer method agree with the simulated values within a factor of < 2 for the majority of gas structures. We next compare synthetic spectra with observations from the 21-SPONGE survey at the Karl G. Jansky Very Large Array using AGD. We find more components per line of sight in 21-SPONGE than in synthetic spectra, which reflects insufficient simulated gas scale heights and the limitations of local box simulations. In addition, we find a significant population of low-optical depth, broad absorption components in the synthetic data which are not seen in 21-SPONGE. This population is not obvious in integrated or per-channel diagnostics, and reflects the benefit of studying velocity-resolved components. The discrepant components correspond to the highest spin temperatures (1000 < Ts < 4000 K), which are not seen in 21-SPONGE despite sufficient observational sensitivity. We demonstrate that our analysis method is a powerful tool for diagnosing neutral ISM conditions, and future work is needed to improve observational statistics and implementation of simulated physics. △ Less

Submitted 27 January, 2017; v1 submitted 6 December, 2016; originally announced December 2016.

Comments: Accepted for publication in ApJ. Revised to incorporate referee comments. 18 pages, 18 figures, 1 table

The Lognormal Probability Distribution Function of the Perseus Molecular Cloud: A Comparison of HI and Dust

Authors: Blakesley Burkhart, Min-Young Lee, Claire Murray, Snezana Stanimirovic

Abstract: The shape of the probability distribution function (PDF) of molecular clouds is an important ingredient for modern theories of star formation and turbulence. Recently, several studies have pointed out observational difficulties with constraining the low column density (i.e. Av <1) PDF using dust tracers. In order to constrain the shape and properties of the low column density probability distribut… ▽ More The shape of the probability distribution function (PDF) of molecular clouds is an important ingredient for modern theories of star formation and turbulence. Recently, several studies have pointed out observational difficulties with constraining the low column density (i.e. Av <1) PDF using dust tracers. In order to constrain the shape and properties of the low column density probability distribution function, we investigate the PDF of multiphase atomic gas in the Perseus molecular cloud using opacity-corrected GALFA-HI data and compare the PDF shape and properties to the total gas PDF and the N(H2) PDF. We find that the shape of the PDF in the atomic medium of Perseus is well described by a lognormal distribution, and not by a power-law or bimodal distribution. The peak of the atomic gas PDF in and around Perseus lies at the HI-H2 transition column density for this cloud, past which the N(H2) PDF takes on a powerlaw form. We find that the PDF of the atomic gas is narrow and at column densities larger than the HI-H2 transition the HI rapidly depletes, suggesting that the HI PDF may be used to find the HI-H2 transition column density. We also calculate the sonic Mach number of the atomic gas by using HI absorption line data, which yields a median value of Ms=4.0 for the CNM, while the HI emission PDF, which traces both the WNM and CNM, has a width more consistent with transonic turbulence. △ Less

Submitted 9 September, 2015; originally announced September 2015.

Comments: Accepted to ApJL; comments welcome

The Turbulence Velocity Power Spectrum of Neutral Hydrogen in the Small Magellanic Cloud

Authors: Alexey Chepurnov, Blakesley Burkhart, Alex Lazarian, Snezana Stanimirovic

Abstract: We present the results of the Velocity Coordinate Spectrum (VCS) technique to calculate the velocity power spectrum of turbulence in the Small Magellanic Cloud (SMC) in 21cm emission. We have obtained a velocity spectral index of -3.85 and an injection scale of 2.3 kpc. The spectral index is steeper than the Kolmogorov index which is expected for shock-dominated turbulence which is in agreement wi… ▽ More We present the results of the Velocity Coordinate Spectrum (VCS) technique to calculate the velocity power spectrum of turbulence in the Small Magellanic Cloud (SMC) in 21cm emission. We have obtained a velocity spectral index of -3.85 and an injection scale of 2.3 kpc. The spectral index is steeper than the Kolmogorov index which is expected for shock-dominated turbulence which is in agreement with past works on the SMC gas dynamics. The injection scale of 2.3 kpc suggests that tidal interactions with the Large Magellanic Cloud are the dominate driver of turbulence in this dwarf galaxy. This implies turbulence maybe driven by multiple mechanisms in galaxies in addition to supernova injection and that galaxy-galaxy interactions may play an important role. △ Less

Submitted 10 June, 2015; originally announced June 2015.

Comments: Accepted to ApJ

First Detection of HCO$^+$ Absorption in the Magellanic System

Authors: Claire E. Murray, Snežana Stanimirović, N. M. McClure-Griffiths, M. E. Putman, H. S. Liszt, Tony Wong, P. Richter, J. R. Dawson, John M. Dickey, Robert R. Lindner, Brian L. Babler, J. R. Allison

Abstract: We present the first detection of HCO$^+$ absorption in the Magellanic System. Using the Australia Telescope Compact Array (ATCA), we observed 9 extragalactic radio continuum sources behind the Magellanic System and detected HCO$^+$ absorption towards one source located behind the leading edge of the Magellanic Bridge. The detection is located at LSR velocity of $v=214.0 \pm 0.4\rm\,km\,s^{-1}$, w… ▽ More We present the first detection of HCO$^+$ absorption in the Magellanic System. Using the Australia Telescope Compact Array (ATCA), we observed 9 extragalactic radio continuum sources behind the Magellanic System and detected HCO$^+$ absorption towards one source located behind the leading edge of the Magellanic Bridge. The detection is located at LSR velocity of $v=214.0 \pm 0.4\rm\,km\,s^{-1}$, with a full width at half maximum of $Δv=4.5\pm 1.0\rm\,km\,s^{-1}$ and optical depth of $τ(\rm HCO^+)=0.10\pm 0.02$. Although there is abundant neutral hydrogen (HI) surrounding the sightline in position-velocity space, at the exact location of the absorber the HI column density is low, $<10^{20}\rm\,cm^{-2}$, and there is little evidence for dust or CO emission from Planck observations. While the origin and survival of molecules in such a diffuse environment remains unclear, dynamical events such as HI flows and cloud collisions in this interacting system likely play an important role. △ Less

Submitted 10 June, 2015; originally announced June 2015.

Comments: Accepted for publication in ApJ. 6 pages, 2 figures, 2 tables

Cold and Warm Atomic Gas around the Perseus Molecular Cloud. II. The Impact of High Optical Depth on the HI Column Density Distribution and Its Implication for the HI-to-H2 Transition

Authors: Min-Young Lee, Snezana Stanimirovic, Claire E. Murray, Carl Heiles, Jesse Miller

Abstract: We investigate the impact of high optical depth on the HI saturation observed in the Perseus molecular cloud by using Arecibo HI emission and absorption measurements toward 26 radio continuum sources. The spin temperature and optical depth of individual HI components are derived along each line-of-sight, enabling us to estimate the correction for high optical depth. We examine two different method… ▽ More We investigate the impact of high optical depth on the HI saturation observed in the Perseus molecular cloud by using Arecibo HI emission and absorption measurements toward 26 radio continuum sources. The spin temperature and optical depth of individual HI components are derived along each line-of-sight, enabling us to estimate the correction for high optical depth. We examine two different methods for the correction, Gaussian decomposition and isothermal methods, and find that they are consistent (maximum correction factor ~ 1.2) likely due to the relatively low optical depth and insignificant contribution from the diffuse radio continuum emission for Perseus. We apply the correction to the optically thin HI column density on a pixel-by-pixel basis, and find that the total HI mass increases by ~10%. Using the corrected HI column density image and far-infrared data from the IRIS Survey, we then derive the H2 column density on ~0.4 pc scales. For five dark and star-forming sub-regions, the HI surface density is uniform with Sigma_HI ~ 7-9 solar mass/pc2, in agreement with the minimum HI surface density required for shielding H2 against photodissociation. As a result, Sigma_H2/Sigma_HI and Sigma_HI+Sigma_H2 show a tight relation. Our results are consistent with predictions for H2 formation in steady state and chemical equilibrium, and suggest that H2 formation is mainly responsible for the Sigma_HI saturation in Perseus. We also compare the optically thick HI with the observed "CO-dark" gas, and find that the optically thick HI only accounts for ~20% of the "CO-dark" gas in Perseus. △ Less

Submitted 23 May, 2015; v1 submitted 28 April, 2015; originally announced April 2015.

Comments: accepted for publication in the Astrophysical Journal

The 21-SPONGE HI Absorption Survey I: Techniques and Initial Results

Authors: Claire E. Murray, Snežana Stanimirović, W. M. Goss, John M. Dickey, Carl Heiles, Robert R. Lindner, Brian Babler, Nickolas M. Pingel, Allen Lawrence, Jacob Jencson, Patrick Hennebelle

Abstract: We present methods and results from "21-cm Spectral Line Observations of Neutral Gas with the EVLA" (21-SPONGE), a large survey for Galactic neutral hydrogen (HI) absorption with the Karl G. Jansky Very Large Array (VLA). With the upgraded capabilities of the VLA, we reach median root-mean-square (RMS) noise in optical depth of $σ_τ=9\times 10^{-4}$ per $0.42\rm\,km\,s^{-1}$ channel for the 31 sou… ▽ More We present methods and results from "21-cm Spectral Line Observations of Neutral Gas with the EVLA" (21-SPONGE), a large survey for Galactic neutral hydrogen (HI) absorption with the Karl G. Jansky Very Large Array (VLA). With the upgraded capabilities of the VLA, we reach median root-mean-square (RMS) noise in optical depth of $σ_τ=9\times 10^{-4}$ per $0.42\rm\,km\,s^{-1}$ channel for the 31 sources presented here. Upon completion, 21-SPONGE will be the largest HI absorption survey with this high sensitivity. We discuss the observations and data reduction strategies, as well as line fitting techniques. We prove that the VLA bandpass is stable enough to detect broad, shallow lines associated with warm HI, and show that bandpass observations can be combined in time to reduce spectral noise. In combination with matching HI emission profiles from the Arecibo Observatory ($\sim3.5'$ angular resolution), we estimate excitation (or spin) temperatures ($\rm T_s$) and column densities for Gaussian components fitted to sightlines along which we detect HI absorption (30/31). We measure temperatures up to $\rm T_s\sim1500\rm\,K$ for individual lines, showing that we can probe the thermally unstable interstellar medium (ISM) directly. However, we detect fewer of these thermally unstable components than expected from previous observational studies. We probe a wide range in column density between $\sim10^{16}$ and $>10^{21}\rm\,cm^{-2}$ for individual HI clouds. In addition, we reproduce the trend between cold gas fraction and average $\rm T_s$ found by synthetic observations of a hydrodynamic ISM simulation by Kim et al. (2014). Finally, we investigate methods for estimating HI $\rm T_s$ and discuss their biases. △ Less

Submitted 3 March, 2015; originally announced March 2015.

Comments: Accepted for publication in ApJ; 24 pages, 14 figures

Galactic and Magellanic Evolution with the SKA

Authors: Naomi M. McClure-Griffiths, Snezana Stanimirovic, Claire E. Murray, Di Li, John M. Dickey, Enrique Vazquez-Semadeni, Josh E. G. Peek, Mary Putman, Susan E. Clark, Marc-Antoine Miville-Deschenes, Joss Bland-Hawthorn, Lister Staveley-Smith

Abstract: As we strive to understand how galaxies evolve it is crucial that we resolve physical processes and test emerging theories in nearby systems that we can observe in great detail. Our own Galaxy, the Milky Way, and the nearby Magellanic Clouds provide unique windows into the evolution of galaxies, each with its own metallicity and star formation rate. These laboratories allow us to study with more d… ▽ More As we strive to understand how galaxies evolve it is crucial that we resolve physical processes and test emerging theories in nearby systems that we can observe in great detail. Our own Galaxy, the Milky Way, and the nearby Magellanic Clouds provide unique windows into the evolution of galaxies, each with its own metallicity and star formation rate. These laboratories allow us to study with more detail than anywhere else in the Universe how galaxies acquire fresh gas to fuel their continuing star formation, how they exchange gas with the surrounding intergalactic medium, and turn warm, diffuse gas into molecular clouds and ultimately stars. The $λ$21-cm line of atomic hydrogen (HI) is an excellent tracer of these physical processes. With the SKA we will finally have the combination of surface brightness sensitivity, point source sensitivity and angular resolution to transform our understanding of the evolution of gas in the Milky Way, all the way from the halo down to the formation of individual molecular clouds. △ Less

Submitted 6 January, 2015; originally announced January 2015.

Comments: 25 pages, from "Advancing Astrophysics with the Square Kilometre Array", to appear in Proceedings of Science

Autonomous Gaussian Decomposition

Authors: Robert R. Lindner, Carlos Vera-Ciro, Claire E. Murray, Snežana Stanimirović, Brian L. Babler, Carl Heiles, Patrick Hennebelle, W. M. Goss, John Dickey

Abstract: We present a new algorithm, named Autonomous Gaussian Decomposition (AGD), for automatically decomposing spectra into Gaussian components. AGD uses derivative spectroscopy and machine learning to provide optimized guesses for the number of Gaussian components in the data, and also their locations, widths, and amplitudes. We test AGD and find that it produces results comparable to human-derived sol… ▽ More We present a new algorithm, named Autonomous Gaussian Decomposition (AGD), for automatically decomposing spectra into Gaussian components. AGD uses derivative spectroscopy and machine learning to provide optimized guesses for the number of Gaussian components in the data, and also their locations, widths, and amplitudes. We test AGD and find that it produces results comparable to human-derived solutions on 21cm absorption spectra from the 21cm SPectral line Observations of Neutral Gas with the EVLA (21-SPONGE) survey. We use AGD with Monte Carlo methods to derive the HI line completeness as a function of peak optical depth and velocity width for the 21-SPONGE data, and also show that the results of AGD are stable against varying observational noise intensity. The autonomy and computational efficiency of the method over traditional manual Gaussian fits allow for truly unbiased comparisons between observations and simulations, and for the ability to scale up and interpret the very large data volumes from the upcoming Square Kilometer Array and pathfinder telescopes. △ Less

Submitted 9 September, 2014; originally announced September 2014.

Comments: 12 pages, 8 figures, submitted to AJ

Cold and warm atomic gas around the Perseus molecular cloud I: Basic Properties

Authors: Snezana Stanimirovic, Claire E. Murray, Min-Young Lee, Carl Heiles, Jesse Miller

Abstract: (Abridged) Using the Arecibo Observatory we have obtained neutral hydrogen (HI) absorption and emission spectral pairs in the direction of 26 background radio continuum sources in the vicinity of the Perseus molecular cloud. Strong absorption lines were detected in all cases allowing us to estimate spin temperature (T_s) and optical depth for 107 individual Gaussian components along these lines of… ▽ More (Abridged) Using the Arecibo Observatory we have obtained neutral hydrogen (HI) absorption and emission spectral pairs in the direction of 26 background radio continuum sources in the vicinity of the Perseus molecular cloud. Strong absorption lines were detected in all cases allowing us to estimate spin temperature (T_s) and optical depth for 107 individual Gaussian components along these lines of sight. Basic properties of individual HI clouds (spin temperature, optical depth, and the column density of the cold and warm neutral medium, CNM and WNM) in and around Perseus are very similar to those found for random interstellar lines of sight sampled by the Millennium HI survey. This suggests that the neutral gas found in and around molecular clouds is not atypical. However, lines of sight in the vicinity of Perseus have on average a higher total HI column density and the CNM fraction, suggesting an enhanced amount of cold HI relative to an average interstellar field. Our estimated optical depth and spin temperature are in stark contrast with the recent attempt at using Planck data to estimate properties of the optically thick HI. Only ~15% of lines of sight in our study have a column density weighted average spin temperature lower than 50 K, in comparison with >85% of Planck's sky coverage. The observed CNM fraction is inversely proportional to the optical-depth weighted average spin temperature, in excellent agreement with the recent numerical simulations by Kim et al. While the CNM fraction is on average higher around Perseus relative to a random interstellar field, it is generally low, 10-50%. This suggests that extended WNM envelopes around molecular clouds, and/or significant mixing of CNM and WNM throughout molecular clouds, are present and should be considered in the models of molecule and star formation. △ Less

Submitted 29 July, 2014; originally announced July 2014.

Comments: 23 pages; accepted for publication in ApJ

The Role of Feedback in Shaping the Structure of the Interstellar Medium

Authors: A. P. Walker, B. K. Gibson, K. Pilkington, C. B. Brook, P. Dutta, S. Stanimirovic, G. S. Stinson, J. Bailin

Abstract: We present an analysis of the role of feedback in shaping the neutral hydrogen (HI) content of simulated disc galaxies. For our analysis, we have used two realisations of two separate Milky Way-like (~L*) discs - one employing a conservative feedback scheme (MUGS), the other significantly more energetic (MaGICC). To quantify the impact of these schemes, we generate zeroth moment (surface density)… ▽ More We present an analysis of the role of feedback in shaping the neutral hydrogen (HI) content of simulated disc galaxies. For our analysis, we have used two realisations of two separate Milky Way-like (~L*) discs - one employing a conservative feedback scheme (MUGS), the other significantly more energetic (MaGICC). To quantify the impact of these schemes, we generate zeroth moment (surface density) maps of the inferred HI distribution; construct power spectra associated with the underlying structure of the simulated cold ISM, in addition to their radial surface density and velocity dispersion profiles. Our results are compared with a parallel, self-consistent, analysis of empirical data from THINGS (The HI Nearby Galaxy Survey). Single power-law fits (P~k^gamma) to the power spectra of the stronger-feedback (MaGICC) runs (over spatial scales corresponding to 0.5 kpc to 20 kpc) result in slopes consistent with those seen in the THINGS sample (gamma = -2.5). The weaker-feedback (MUGS) runs exhibit shallower power law slopes (gamma = -1.2). The power spectra of the MaGICC simulations are more consistent though with a two-component fit, with a flatter distribution of power on larger scales (i.e., gamma = -1.4 for scales in excess of 2 kpc) and a steeper slope on scales below 1 kpc (gamma = -5), qualitatively consistent with empirical claims, as well as our earlier work on dwarf discs. The radial HI surface density profiles of the MaGICC discs show a clear exponential behaviour, while those of the MUGS suite are essentially flat; both behaviours are encountered in nature, although the THINGS sample is more consistent with our stronger (MaGICC) feedback runs. △ Less

Submitted 3 March, 2014; originally announced March 2014.

Comments: Accepted for publication in MNRAS; 9 pages; 7 figures

The CO-to-H2 Conversion Factor across the Perseus Molecular Cloud

Authors: Min-Young Lee, Snezana Stanimirovic, Mark G. Wolfire, Rahul Shetty, Simon C. O. Glover, Faviola Z. Molina, Ralf S. Klessen

Abstract: We derive the CO-to-H2 conversion factor, X_CO = N(H2)/I_CO, across the Perseus molecular cloud on sub-parsec scales by combining the dust-based N(H2) data with the I_CO data from the COMPLETE Survey. We estimate an average X_CO ~ 3 x 10^19 cm^-2 K^-1 km^-1 s and find a factor of ~3 variations in X_CO between the five sub-regions in Perseus. Within the individual regions, X_CO varies by a factor o… ▽ More We derive the CO-to-H2 conversion factor, X_CO = N(H2)/I_CO, across the Perseus molecular cloud on sub-parsec scales by combining the dust-based N(H2) data with the I_CO data from the COMPLETE Survey. We estimate an average X_CO ~ 3 x 10^19 cm^-2 K^-1 km^-1 s and find a factor of ~3 variations in X_CO between the five sub-regions in Perseus. Within the individual regions, X_CO varies by a factor of ~100, suggesting that X_CO strongly depends on local conditions in the interstellar medium. We find that X_CO sharply decreases at Av < 3 mag but gradually increases at Av > 3 mag, with the transition occurring at Av where I_CO becomes optically thick. We compare the N(HI), N(H2), I_CO, and X_CO distributions with two models of the formation of molecular gas, a one-dimensional photodissociation region (PDR) model and a three-dimensional magnetohydrodynamic (MHD) model tracking both the dynamical and chemical evolution of gas. The PDR model based on the steady state and equilibrium chemistry reproduces our data very well but requires a diffuse halo to match the observed N(HI) and I_CO distributions. The MHD model generally matches our data well, suggesting that time-dependent effects on H2 and CO formation are insignificant for an evolved molecular cloud like Perseus. However, we find interesting discrepancies, including a broader range of N(HI), likely underestimated I_CO, and a large scatter of I_CO at small Av. These discrepancies likely result from strong compressions/rarefactions and density fluctuations in the MHD model. △ Less

Submitted 20 January, 2014; originally announced January 2014.

Comments: accepted for publication in the Astrophysical Journal

Excitation temperature of the warm neutral medium as a new probe of the Lyman-α radiation field

Authors: Claire E. Murray, Robert R. Lindner, Snezana Stanimirovic, W. M. Goss, Carl Heiles, John Dickey, Nickolas M. Pingel, Allen Lawrence, Jacob Jencson, Brian L. Babler, Patrick Hennebelle

Abstract: We use the Karl G. Jansky Very Large Array (VLA) to conduct a high-sensitivity survey of neutral hydrogen (HI) absorption in the Milky Way. In combination with corresponding HI emission spectra obtained mostly with the Arecibo Observatory, we detect a widespread warm neutral medium (WNM) component with excitation temperature <Ts>= 7200 (+1800,-1200) K (68% confidence). This temperature lies above… ▽ More We use the Karl G. Jansky Very Large Array (VLA) to conduct a high-sensitivity survey of neutral hydrogen (HI) absorption in the Milky Way. In combination with corresponding HI emission spectra obtained mostly with the Arecibo Observatory, we detect a widespread warm neutral medium (WNM) component with excitation temperature <Ts>= 7200 (+1800,-1200) K (68% confidence). This temperature lies above theoretical predictions based on collisional excitation alone, implying that Ly-α scattering, the most probable additional source of excitation, is more important in the interstellar medium (ISM) than previously assumed. Our results demonstrate that HI absorption can be used to constrain the Ly-α radiation field, a critical quantity for studying the energy balance in the ISM and intergalactic medium yet notoriously difficult to model because of its complicated radiative transfer, in and around galaxies nearby and at high redshift. △ Less

Submitted 6 January, 2014; originally announced January 2014.

Comments: 7 pages, 3 figures, 1 table. Accepted for publication in ApJ Letters

Evidence for an Interaction in the Nearest Starbursting Dwarf Irregular Galaxy IC 10

Authors: David L. Nidever, Trisha Ashley, Colin T. Slater, Juergen Ott, Megan Johnson, Eric F. Bell, Snezana Stanimirovic, Mary Putman, Steven R. Majewski, Caroline E. Simpson, Eva Juette, Tom A. Oosterloo, W. Butler Burton

Abstract: Using deep 21-cm HI data from the Green Bank Telescope we have detected an ~18.3 kpc-long gaseous extension associated with the starbursting dwarf galaxy IC 10. The newly-found feature stretches 1.3 deg to the northwest and has a large radial velocity gradient reaching to ~65 km/s lower than the IC 10 systemic velocity. A region of higher column density at the end of the extension that possesses a… ▽ More Using deep 21-cm HI data from the Green Bank Telescope we have detected an ~18.3 kpc-long gaseous extension associated with the starbursting dwarf galaxy IC 10. The newly-found feature stretches 1.3 deg to the northwest and has a large radial velocity gradient reaching to ~65 km/s lower than the IC 10 systemic velocity. A region of higher column density at the end of the extension that possesses a coherent velocity gradient (~10 km/s across ~26 arcmin) transverse to the extension suggests rotation and may be a satellite galaxy of IC 10. The HI mass of IC 10 is 9.5x10^7 (d/805 kpc)^2 Msun and the mass of the new extension is 7.1x10^5 (d/805 kpc)^2 Msun. An IC 10-M31 orbit using known radial velocity and proper motion values for IC 10 show that the HI extension is inconsistent with the trailing portion of the orbit so that an M31-tidal or ram pressure origin seems unlikely. We argue that the most plausible explanation for the new feature is that it is the result of a recent interaction (and possible late merger) with another dwarf galaxy. This interaction could not only have triggered the origin of the recent starburst in IC 10, but could also explain the existence of previously-found counter-rotating HI gas in the periphery of the IC 10 which was interpreted as originating from primordial gas infall. △ Less

Submitted 28 October, 2013; originally announced October 2013.

Comments: 7 pages, 4 figures, 1 table. Accepted for publication in ApJ Letters