Search

We use the Milky Way neutral hydrogen (HI) absorption and emission spectra from the Galactic Australian Square Kilometre Array Pathfinder (GASKAP) Phase II Pilot survey along with toy models to investigate the effects of stacking multicomponent spectra on measurements of peak optical depth and spin temperature. Shifting spectra by the peak in emission, 'primary' components shifted to 0 km s$^{-1}$ are correctly averaged. Additional components on individual sightlines are averaged with non-centred velocities, producing a broader and shallower 'secondary' component in the resulting stack. Peak optical depths and brightness temperatures of the secondary components from stacks are lower limits of their true average values due to the velocity offset of each component. The spin temperature however is well correlated with the truth since the velocity offset of components affects the emission and absorption spectra equally. Stacking 462 GASKAP absorption-emission spectral pairs, we detect a component with a spin temperature of 1320 $\pm$ 263 K, consistent with gas from the unstable neutral medium and higher than any previous GASKAP detection in this region. We also stack 2240 pilot survey spectra containing no Milky Way absorption, revealing a primary narrow and secondary broad component, with spin temperatures belonging to the cold neutral medium (CNM). Spatially binning and stacking the non-detections across the plane-of-sky by their distance from CNM absorption detections, the primary component's optical depth decreases with distance from known locations of cold gas. The spin temperature however remains stable in both components, over an approximate physical plane-of-sky distance of $\sim$ 100 pc., Comment: This paper has 16 pages, 17 figures. This paper has been accepted for publication in Monthly Notices of the Royal Astronomical Society Main Journal

We present a catalog of individual stellar and dust extinction properties along close to 500,000 sight lines in the southwest bar of the Small Magellanic Cloud (SMC). The catalog is based on multiband Hubble Space Telescope photometric data spanning near-ultraviolet to near-infrared wavelengths from the Small Magellanic Cloud Investigation of Dust and Gas Evolution survey (SMIDGE) covering a 100 x 200 pc area. We use the probabilistic technique of the Bayesian Extinction And Stellar Tool (BEAST) to model the spectral energy distributions of individual stars in SMIDGE and include the effects of observational uncertainties in the data. We compare BEAST-derived dust extinction properties with tracers of the interstellar medium, such as the emission from the 12CO (2-1) transition (I(CO)), the dust mass surface density ({\Sigma}dust) from far-IR emission, the H I column density (N(HI)) from the 21cm transition, and the mass fraction of polycyclic aromatic hydrocarbons (PAHs; qPAH, derived from IR emission). We find that the dust extinction (A(V )) in the SMIDGE field is strongly correlated with {\Sigma}dust and I(CO), and less so with N(HI) and qPAH, and suggest potential explanations. Our extinction measurements are also sensitive to the presence of the 2175 {\AA} bump in the extinction curve toward UV bright stars. While most do not show evidence for the bump, we identify ~200 lines of sight that are 2175 {\AA} bump candidates. Furthermore, we find distinct structures in the dust extinction-distance distributions that provide insights into the 3D geometry of the SMC., Comment: 20 pages, 13 figures, 4 tables

The Local Ultraviolet to Infrared Treasury (LUVIT) is a Hubble Space Telescope program that combines newly acquired data in the near ultraviolet (NUV), optical, and near infrared (NIR) with archival optical and NIR imaging to produce multiband panchromatic resolved stellar catalogs for 23 pointings in 22 low-mass, star-forming galaxies ranging in distance from the outskirts of the Local Group to ~3.8 Mpc. We describe the survey design, detail the LUVIT broadband filter observations and the archival datasets included in the LUVIT reductions, and summarize the simultaneous multiband data reduction steps. The spatial distributions and color-magnitude diagrams (CMDs) from the resulting stellar catalogs are presented for each target, from the NUV to the NIR. We demonstrate in which regions of the CMDs stars with NUV and optical, optical and NIR, and NUV through NIR detections reside. For each target, we use the results from artificial star tests to measure representative completeness, bias, and total photometric uncertainty as a function of magnitude in each broadband filter. We also assess which LUVIT targets have significant spatial variation in the fraction of stars recovered at a given magnitude. The panchromatic LUVIT stellar catalogs will provide a rich legacy dataset for a host of resolved stellar population studies., Comment: 48 pages, 14 figures, 8 tables, published in ApJS

By analyzing the spectral energy distributions (SEDs) of resolved stars in nearby galaxies, we can constrain their stellar properties and line-of-sight dust extinction. From the Scylla survey, we obtain ultraviolet to near-infrared photometry from Wide Field Camera 3 onboard the {\it Hubble Space Telescope} for more than 1.5 million stars in the SMC and LMC. We use the Bayesian Extinction and Stellar Tool (BEAST) to analyze the multi-band SEDs of these sources and characterize their initial masses, ages, metallicities, distances, and line-of-sight extinction properties (e.g.~$A_V$, $R_V$). We apply quality cuts and perform validation simulations to construct a catalog of over 550,000 stars with high-reliability SED fits, which we use to analyze the stellar content and extinction properties of the SMC and LMC. We detect stars with masses as low as 0.6 $M_{\odot}$. Observed stellar age distributions show a jump in stars around 6 Gyrs ago, which is in agreement with other star-formation histories. Extinctions ($A_V$) in both galaxies follow a log-normal distribution. We compare $A_V$ with ancillary gas and dust tracers like $HI$, $H_\alpha$, and far infrared (FIR) dust emission and find positive correlations on a field-by-field basis. We convert observed $A_V$ to predicted dust surface densities using the Draine et. al. (2014) model and find $A_V$-based dust surface densities are a factor of $\sim$2.5 lower than observed FIR-based dust surface densities, a correction factor similar to other studies., Comment: Submitted to ApJ, 31 pages

The proximity of the Large and Small Magellanic Clouds (LMC and SMC) provides the opportunity to study the impact of dwarf-dwarf interactions on their mass assembly with a unique level of detail. To this end, we analyze two-filter broadband imaging of 83 Hubble Space Telescope (HST) pointings covering 0.203 deg$^2$ towards the SMC, extending out to $\sim$3.5 kpc in projection from its optical center. Lifetime star formation histories (SFHs) fit to each pointing independently reveal an outside-in age gradient such that fields in the SMC outskirts are older on average. We measure radial gradients of the lookback time to form 90%, 75% and 50% of the cumulative stellar mass for the first time, finding $\delta$($\tau_{90}$, $\tau_{75}$, $\tau_{50}$)/$\delta$R = (0.61$^{+0.08}_{-0.07}$, 0.65$^{+0.09}_{-0.08}$, 0.82$^{+0.12}_{-0.16}$) Gyr/kpc assuming PARSEC evolutionary models and a commonly used elliptical geometry of the SMC, although our results are robust to these assumptions. The wing of the SMC deviates from this trend, forming 25\% of its cumulative mass over the most recent 3 Gyr due to a best-fit star formation rate that remains approximately constant. Our results are consistent with chemodynamical evidence of a tidally stripped SMC component in the foreground, and imply contributions to the observed SFH from multiple previous LMC-SMC interactions. We also compare our SMC SFH with results from a companion study of the LMC, finding that while the two galaxies present different internal, spatially resolved SFH trends, both the LMC and SMC have similar near-constant lifetime SFHs when viewed globally., Comment: ApJ in press. 40 pages, 18 figures, 5 tables including Appendices

The proximity of the Magellanic Clouds provides the opportunity to study interacting dwarf galaxies near a massive host, and spatial trends in their stellar population properties in particular, with a unique level of detail. The Scylla pure parallel program has obtained deep (80% complete to >1 mag below the ancient main sequence turnoff), homogeneous two-filter Hubble Space Telescope (HST) imaging sampling the inner star-forming disk of the Large Magellanic Cloud (LMC), the perfect complement to shallower, contiguous ground-based surveys. We harness this imaging together with extant archival data and fit lifetime star formation histories (SFHs) to resolved color-magnitude diagrams (CMDs) of 111 individual fields, using three different stellar evolutionary libraries. We validate per-field recovered distances and extinctions as well as the combined global LMC age-metallicity relation and SFH against independent estimates. We find that the present-day radial age gradient reverses from an inside-out gradient in the inner disk to an outside-in gradient beyond $\sim$2 disk scalelengths, supported by ground-based measurements. The gradients become relatively flatter at earlier lookback times, while the location of the inversion remains constant over an order of magnitude in lookback time, from $\sim$1$-$10 Gyr. This suggests at least one mechanism that predates the recent intense LMC-SMC interaction. We compare observed radial age trends to other late-type galaxies at fixed stellar mass and discuss similarities and differences in the context of potential drivers, implying strong radial migration in the LMC., Comment: ApJ in press. 45 pages, 17 figures, 9 tables including Appendices

Scylla is a deep Hubble Space Telescope survey of the stellar populations, interstellar medium and star formation in the LMC and SMC. As a pure-parallel complement to the Ultraviolet Legacy Library of Young Stars as Essential Standards (ULLYSES) survey, Scylla obtained 342 orbits of ultraviolet (UV) through near-infrared (IR) imaging of the LMC and SMC with Wide Field Camera 3. In this paper, we describe the science objectives, observing strategy, data reduction procedure, and initial results from our photometric analysis of 96 observed fields. Although our observations were constrained by ULYSSES primary exposures, we imaged all fields in at least two filters (F475W and F814W), and 64% of fields in at least three and as many as seven WFC3 filters spanning the UV to IR. Overall, we reach average 50% completeness of $m_{\rm F225W}=26.0$, $m_{\rm F275W}=26.2$, $m_{\rm F336W}=26.9$, $m_{\rm F475W}=27.8$, $m_{\rm F814W}=25.5$, $m_{\rm F110W}=24.7$, and $m_{\rm F160W}=24.0$ Vega magnitudes in our photometric catalogs, which is faintward of the ancient main sequence turnoff in all filters. The primary science goals of Scylla include characterizing the structure and properties of dust in the MCs, as well as their spatially-resolved star formation and chemical enrichment histories. Our images and photometric catalogs, which represent the widest-area coverage of MCs with HST photometry to date, are available as a high-level science product at the Barbara A. Mikulski Archive for Space Telescopes., Comment: 25 pages, 16 figures, 8 tables. Accepted for publication in ApJS

We present the largest Galactic neutral hydrogen HI absorption survey to date, utilizing the Australian SKA Pathfinder Telescope at an unprecedented spatial resolution of 30''. This survey, GASKAP-HI, unbiasedly targets 2,714 continuum background sources over 250 square degrees in the direction of the Magellanic Clouds, a significant increase compared to a total of 373 sources observed by previous Galactic absorption surveys across the entire Milky Way. We aim to investigate the physical properties of cold (CNM) and warm (WNM) neutral atomic gas in the Milky Way foreground, characterized by two prominent filaments at high Galactic latitudes (between $-45^{\circ}$ and $-25^{\circ}$). We detected strong HI absorption along 462 lines of sight above the 3$\sigma$ threshold, achieving an absorption detection rate of 17%. GASKAP-HI's unprecedented angular resolution allows for simultaneous absorption and emission measurements to sample almost the same gas clouds along a line of sight. A joint Gaussian decomposition is then applied to absorption-emission spectra to provide direct estimates of HI optical depths, temperatures, and column densities for the CNM and WNM components. The thermal properties of CNM components are consistent with those previously observed along a wide range of Solar neighborhood environments, indicating that cold HI properties are widely prevalent throughout the local interstellar medium. Across our region of interest, CNM accounts for ~30% of the total HI gas, with the CNM fraction increasing with column density toward the two filaments. Our analysis reveals an anti-correlation between CNM temperature and its optical depth, which implies that CNM with lower optical depth leads to a higher temperature., Comment: Largest Galactic HI Absorption Survey To Date (GASKAP-HI): Cold Atomic Gas in the Magellanic Cloud foreground using Australian SKA Pathfinder. This paper has 19 pages, 17 figures. This paper has been accepted for publication in Monthly Notices of the Royal Astronomical Society Main Journal

The Small Magellanic Cloud (SMC) shows a large variation in ultraviolet (UV) dust extinction curves, ranging from Milky Way-like (MW) to significantly steeper curves with no detectable 2175 A bump. This result is based on a sample of only nine sightlines. From HST/STIS and IUE spectra of OB stars, we have measured UV extinction curves along 32 SMC sightlines where eight of these curves were published previously. We find 16 sightlines with steep extinction with no detectable 2175 A bump, four sightlines with MW-like extinction with a detectable 2175 A bump, two sightlines with fairly flat UV extinction and weak/absent 2175 A bumps, and 10 sightlines with unreliable curves due to low SMC dust columns. Our expanded sample shows that the sightlines with and without the 2175 A bump are located throughout the SMC and not limited to specific regions. The average extinction curve of the 16 bumpless sightlines is very similar to the previous average based on four sightlines. We find no correlation between dust column and the strength of the 2175 A bump. We test the hypothesis that the 2175 A bump is due to the same dust grains that are responsible for the mid-infrared carbonaceous (PAH) emission features and find they are correlated, confirming recent work in the MW. Overall, the slope of the UV extinction increases as the amplitudes of the 2175 A bump and far-UV curvature decrease. Finally, the UV slope is correlated with $N(HI)/A(V)$ and the 2175 A bump and nonlinear far-UV rise amplitudes are anti-correlated with $N(HI)/A(V)$., Comment: 24 pages, 12 figures, ApJ, accepted

We leverage recent numerical simulations of highly resolved star-forming regions in a Milky Way-like Galaxy to explore the nature of extended gaseous envelopes around molecular clouds. We extract a sample of two dozen star-forming clouds from the feedback-dominated suite of the "Cloud Factory'' simulations. With the goal of exploring the 3D thermal and chemical structure of the gas, we measure and fit the clouds' radial profiles in multiple tracers, including $n_{H_I}$, $n_{H_2}$, $n_{H_{tot}}$, $n_{CO}$, and gas temperature. We find that while solar neighborhood clouds recently detected via 3D dust mapping have radially symmetric, low-density envelopes that extend $\sim$ 10-15 pc, the simulated cloud envelopes are primarily radially asymmetric with low-density envelopes that extend only $\sim$ 2-3 pc. One potential explanation for the absence of extended envelopes in the simulated clouds may be the lack of magnetic fields, while a stronger local feedback prescription compared to solar neighborhood conditions may drive the radially asymmetric cloud morphologies. We make the pipeline to extract and characterize the radial profile of clouds publicly available, which can be used in complementary and future simulations to shed additional light on the key physics shaping the formation and evolution of star-forming structures in the Milky Way., Comment: 22 pages, 13 figures, 2 tables. Accepted for publication in ApJ (March 2024)

The structure and dynamics of the star-forming disk of the Small Magellanic Cloud (SMC) have long confounded us. The SMC is widely used as a prototype for galactic physics at low metallicity, and yet we fundamentally lack an understanding of the structure of its interstellar medium (ISM). In this work, we present a new model for the SMC by comparing the kinematics of young, massive stars with the structure of the ISM traced by high-resolution observations of neutral atomic hydrogen (HI) from the Galactic Australian Square Kilometer Array Pathfinder survey (GASKAP-HI). Specifically, we identify thousands of young, massive stars with precise radial velocity constraints from the Gaia and APOGEE surveys and match these stars to the ISM structures in which they likely formed. By comparing the average dust extinction towards these stars, we find evidence that the SMC is composed of two structures with distinct stellar and gaseous chemical compositions. We construct a simple model that successfully reproduces the observations and shows that the ISM of the SMC is arranged into two, superimposed, star-forming systems with similar gas mass separated by ~5 kpc along the line of sight., Comment: ApJ accepted. 20 pages, 18 figures

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

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., Comment: Accepted for publication to PASA. 41 pages, 27 figures

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., Comment: 15 pages, 3 figures, 1 table. Accepted to ApJ

Using new large area maps of the cold neutral medium (CNM) fraction, $f_{\rm CNM}$, we investigate the relationship between the CNM, the abundance of polycyclic aromatic hydrocarbons (PAHs), and the anomalous microwave emission (AME). We first present our $f_{\rm CNM}$ map based on full-sky HI4PI data, using a convolutional neural network to covert the spectroscopic HI data to $f_{\rm CNM}$. We demonstrate that $f_{\rm CNM}$ is strongly correlated with the fraction of dust in PAHs as estimated from mid- and far-infrared dust emission. In contrast, we find no correlation between $f_{\rm CNM}$ and the amount of AME per dust emission, nor between $f_{\rm CNM}$ and the AME peak frequency. These results suggest PAHs preferentially reside in cold, relatively dense gas, perhaps owing to enhanced destruction in more diffuse media. The lack of correlation between $f_{\rm CNM}$ and AME peak frequency is in tension with expectations from theoretical models positing different spectral energy distributions of AME in the cold versus warm neutral medium. We suggest that different PAH abundances and emission physics in different interstellar environments may explain the weaker than expected correlation between 12$\mu$m PAH emission and AME even if PAHs are the AME carriers., Comment: Submitted to ApJ. 12 pages, 4 figures, comments welcome

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, Comment: 31 pages, 12 figures, 8 tables. Accepted to ApJ

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 $600.001$ are ubiquitous ($c\sim100\%$), whereas high optical depths ($\tau>1$) are extremely rare ($c\sim0\%$)., Comment: 33 pages, 16 figures. Accepted for publication in ApJS

Resolving the phase structure of neutral hydrogen (HI) is crucial for understanding the life cycle of the interstellar medium (ISM). However, accurate measurements of HI temperature and density are limited by the availability of background continuum sources for measuring HI absorption. Here we test the use of deep learning for extracting HI properties over large areas without optical depth information. We train a 1D convolutional neural network using synthetic observations of 3D numerical simulations of the ISM to predict the fraction of cold neutral medium (f_CNM) and the correction to the optically-thin HI column density for optical depth (R_HI) from $21\rm\,cm$ emission alone. We restrict our analysis to high Galactic latitudes ($|b|>30 deg), where the complexity of spectral line profiles is minimized. We verify that the network accurately predicts f_CNM and R_HI by comparing the results with direct constraints from 21 cm absorption. By applying the network to the GALFA-HI survey, we generate large-area maps of f_CNM and R_HI. Although the overall contribution to the total HI column of cold neutral medium (CNM)-rich structures is small (~5%), we find that these structures are ubiquitous. Our results are consistent with the picture that small-scale structures observed in 21 cm emission aligned with the magnetic field are dominated by CNM. Finally, we demonstrate that the observed correlation between HI column density and dust reddening (E(B-V)) declines with increasing R_HI, indicating that future efforts to quantify foreground Galactic E(B-V) using HI, even at high latitudes, should increase fidelity by accounting for HI phase structure., Comment: ApJ accepted

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\sigma$ 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., Comment: 29 pages, 19 figures, 4 tables. Accepted to ApJ. Section 4.3.1 updated from previous version

We investigate the kinematics of neutral gas in the Small Magellanic Cloud (SMC) and test the hypothesis that it is rotating in a disk. To trace the 3D motions of the neutral gas distribution, we identify a sample of young, massive stars embedded within it. These are stars with radial velocity measurements from spectroscopic surveys and proper motion measurements from Gaia, whose radial velocities match with dominant HI components. We compare the observed radial and tangential velocities of these stars with predictions from the state-of-the-art rotating disk model based on high-resolution 21 cm observations of the SMC from the Australian Square Kilometer Array Pathfinder telescope. We find that the observed kinematics of gas-tracing stars are inconsistent with disk rotation. We conclude that the kinematics of gas in the SMC are more complex than can be inferred from the integrated radial velocity field. As a result of violent tidal interactions with the LMC, non-rotational motions are prevalent throughout the SMC, and it is likely composed of distinct sub-structures overlapping along the line of sight., Comment: 9 pages, 5 figures, 1 Appendix; ApJ accepted

Our understanding of the dynamics of the interstellar medium is informed by the study of the detailed velocity structure of emission line observations. One approach to study the velocity structure is to decompose the spectra into individual velocity components; this leads to a description of the dataset that is significantly reduced in complexity. However, this decomposition requires full automation lest it becomes prohibitive for large datasets, such as Galactic plane surveys. We developed GaussPy+, a fully automated Gaussian decomposition package that can be applied to emission line datasets, especially large surveys of HI and isotopologues of CO. We built our package upon the existing GaussPy algorithm and significantly improved its performance for noisy data. New functionalities of GaussPy+ include: i) automated preparatory steps, such as an accurate noise estimation, which can also be used as standalone applications; ii) an improved fitting routine; iii) an automated spatial refitting routine that can add spatial coherence to the decomposition results by refitting spectra based on neighbouring fit solutions. We thoroughly tested the performance of GaussPy+ on synthetic spectra and a test field from the Galactic Ring Survey. We found that GaussPy+ can deal with cases of complex emission and even low to moderate signal-to-noise values., Comment: Accepted for publication in A&A. 36 pages, 33 figures (main text: 23 pages, 18 figures; Appendix: 13 pages, 15 figures). The GaussPy+ package is written in Python 3 and is freely available at GitHub under https://github.com/mriener/gausspyplus

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., Comment: This article (22 pages, 23 figures) has been accepted for publication in The Astrophysical Journal (ApJ)

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 ($\sigma_{\tau} < 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

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., Comment: 10 pages, 7 figures; ApJ accepted

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 $\tau_{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 $\sigma_{353}$=$\tau_{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 $\sigma_{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., Comment: The revised manuscript is accepted for publication in The Astrophysical Journal

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\pi}\sigma } \rm{exp}\left[-\frac{[ln(\textit{T}_{ex})-ln(3.4\ K)]^2}{2\sigma^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., Comment: 24 pages, 23 figures. Accepted for publication in ApJS

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., Comment: Accepted for publication in ApJ. Revised to incorporate referee comments. 18 pages, 18 figures, 1 table

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 $\Delta v=4.5\pm 1.0\rm\,km\,s^{-1}$ and optical depth of $\tau(\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., Comment: Accepted for publication in ApJ. 6 pages, 2 figures, 2 tables

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

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 $\sigma_{\tau}=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., Comment: Accepted for publication in ApJ; 24 pages, 14 figures

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 $\lambda$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., Comment: 25 pages, from "Advancing Astrophysics with the Square Kilometre Array", to appear in Proceedings of Science

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., Comment: 12 pages, 8 figures, submitted to AJ

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

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 = 7200 (+1800,-1200) K (68% confidence). This temperature lies above theoretical predictions based on collisional excitation alone, implying that Ly-{\alpha} 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-{\alpha} 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., Comment: 7 pages, 3 figures, 1 table. Accepted for publication in ApJ Letters

We present the first detection of a 6.7 GHz Class II methanol (CH3OH) maser in the Andromeda galaxy (M31). The CH3OH maser was found in a Very Large Array (VLA) survey during the fall of 2009. We have confirmed the methanol maser with the new Expanded VLA (EVLA), in operation since March 2010, but were unsuccessful in detecting a water maser at this location. A direct application for this methanol maser is the determination of the proper motion of M31, such as was obtained with water masers in M33 and IC10 previously. Unraveling the three-dimensional velocity of M31 would solve for the biggest unknown in the modeling of the dynamics and evolution of the Local Group of galaxies., Comment: Letter, accepted by ApJL