We present observations of velocity-resolved [CII] 158 $\mu$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.

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.

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.

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 $\sigma = 7$ km s$^{-1}$ and all other methods giving higher values ($\sigma\sim10$ km s$^{-1}$). All measurements of the line width show a shallow radial trend, with $\sigma$ 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.

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.

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.

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.

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.

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.

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.

We present results of mid-infrared spectroscopic mapping observations of six star-forming regions in the Small Magellanic Cloud from the Spitzer Spectroscopic Survey of the SMC (S4MC). We detect the mid-IR emission from polycyclic aromatic hydrocarbons (PAHs) in all of the mapped regions, greatly increasing the range of environments where PAHs have been spectroscopically detected in the SMC. We investigate the variations of the mid-IR bands in each region and compare our results to studies of the PAH bands in the SINGS sample and in a sample of low-metallicity starburst galaxies. PAH emission in the SMC is characterized by low ratios of the 6-9 micron features relative to the 11.3 micron feature and weak 8.6 and 17.0 micron features. Interpreting these band ratios in the light of laboratory and theoretical studies, we find that PAHs in the SMC tend to be smaller and less ionized than those in higher metallicity galaxies. Based on studies of PAH destruction, we argue that a size distribution shifted towards smaller PAHs cannot be the result of processing in the interstellar medium, but instead reflects differences in the formation of PAHs at low metallicity. Finally, we discuss the implications of our observations for our understanding of the PAH life-cycle in low-metallicity galaxies---namely that the observed deficit of PAHs may be a consequence of PAHs forming with smaller average sizes and therefore being more susceptible to destruction under typical interstellar medium conditions.

The Magellanic Stream (MS) is the nearest example of a gaseous trail formed by interacting galaxies. While the substantial gas masses in these kinds of circumgalactic structures are postulated to represent important sources of fuel for future star formation, the mechanisms whereby this material might be accreted back into galaxies remain unclear. Recent neutral hydrogen (HI) observations have demonstrated that the northern portion of the MS, which probably has been interacting with the Milky Way's hot gaseous halo for close to 1000~Myr, has a larger spatial extent than previously recognized, while also containing significant amounts of small-scale structure. After a brief consideration of the large-scale kinematics of the MS as traced by the recently-discovered extension of the MS, we explore the aging process of the MS gas through the operation of various hydrodynamic instabilities and interstellar turbulence. This in turn leads to consideration of processes whereby MS material survives as cool gas, and yet also evidently fails to form stars. Parallels between the MS and extragalactic tidal features are briefly discussed with an emphasis on steps toward establishing what the MS reveals about the critical role of local processes in determining the evolution of these kinds of systems.

We have applied the unsharp-masking technique to the 24 $\mu$m image of the Small Magellanic Cloud (SMC), obtained with the Spitzer Space Telescope, to search for high-extinction regions. This technique has been used to locate very dense and cold interstellar clouds in the Galaxy, particularly infrared dark clouds (IRDCs). Fifty five candidate regions of high-extinction, namely high-contrast regions (HCRs), have been identified from the generated decremental contrast image of the SMC. Most HCRs are located in the southern bar region and mainly distributed in the outskirts of CO clouds, but most likely contain a significant amount of H2. HCRs have a peak-contrast at 24 $\mu$m of 2 - 2.5 % and a size of 8 - 14 pc. This corresponds to the size of typical and large Galactic IRDCs, but Galactic IRDCs are 2 - 3 times darker at 24 $\mu$m than our HCRs. To constrain the physical properties of the HCRs, we have performed NH3, N2H+, HNC, HCO+, and HCN observations toward one of the HCRs, HCR LIRS36-EAST, using the Australia Telescope Compact Array and the Mopra single-dish radio telescope. We did not detect any molecular line emission, however, our upper limits to the column densities of molecular species suggest that HCRs are most likely moderately dense with n ~ 10^{3} cm-3. This volume density is in agreement with predictions for the cool atomic phase in low metallicity environments. We suggest that HCRs may be tracing clouds at the transition from atomic to molecule-dominated medium, and could be a powerful way to study early stages of gas condensation in low metallicity galaxies. Alternatively, if made up of dense molecular clumps < 0.5 pc in size, HCRs could be counterparts of Galactic IRDCs, and/or regions with highly unusual abundance of very small dust grains.

To understand the impact of low metallicities on giant molecular cloud (GMC) structure, we compare far infrared dust emission, CO emission, and dynamics in the star-forming complex N83 in the Wing of the Small Magellanic Cloud. Dust emission (measured by Spitzer as part of the S3MC and SAGE-SMC surveys) probes the total gas column independent of molecular line emission and traces shielding from photodissociating radiation. We calibrate a method to estimate the dust column using only the high-resolution Spitzer data and verify that dust traces the ISM in the HI-dominated region around N83. This allows us to resolve the relative structures of H2, dust, and CO within a giant molecular cloud complex, one of the first times such a measurement has been made in a low-metallicity galaxy. Our results support the hypothesis that CO is photodissociated while H2 self-shields in the outer parts of low-metallicity GMCs, so that dust/self shielding is the primary factor determining the distribution of CO emission. Four pieces of evidence support this view. First, the CO-to-H2 conversion factor averaged over the whole cloud is very high 4-11 \times 10^21 cm^-2/(K km/s), or 20-55 times the Galactic value. Second, the CO-to-H2 conversion factor varies across the complex, with its lowest (most nearly Galactic) values near the CO peaks. Third, bright CO emission is largely confined to regions of relatively high line-of-sight extinction, A_V >~ 2 mag, in agreement with PDR models and Galactic observations. Fourth, a simple model in which CO emerges from a smaller sphere nested inside a larger cloud can roughly relate the H2 masses measured from CO kinematics and dust.

Radio observations of gas in the Milky Way and Local Group are vital for understanding how galaxies function as systems. The unique sensitivity of Arecibo's 305m dish, coupled with the 7-beam Arecibo L-Band Feed Array (ALFA), provides an unparalleled tool for investigating the full range of interstellar phenomena traced by the HI 21cm line. The GALFA (Galactic ALFA) HI Survey is mapping the entire Arecibo sky over a velocity range of -700 to +700 km/s with 0.2 km/s velocity channels and an angular resolution of 3.4 arcminutes. We present highlights from the TOGS (Turn on GALFA Survey) portion of GALFA-HI, which is covering thousands of square degrees in commensal drift scan observations with the ALFALFA and AGES extragalactic ALFA surveys. This work is supported in part by the National Astronomy and Ionosphere Center, operated by Cornell University under cooperative agreement with the National Science Foundation.

This paper reports results on the statistical analysis of HI turbulence in the Small Magellanic Cloud (SMC). We use 21 cm channel maps, obtained with the Australia Telescope Compact Array and the Parkes telescope, and analyze the spectrum of observed intensity fluctuations as a function of the velocity slice thickness. We confirm predictions by Lazarian & Pogosyan (2000) on the change of the power law index and establish the spectra of 3-D density and velocity. The obtained spectral indices, -3.3 and -3.4, are slightly more shallow than the predictions for the Kolmogorov spectrum. This contrasts to the predictions for the shock-type spectra that are steeper than the Kolmogorov one. The nature of the energy injection in the SMC is unclear as no distinct energy injection scales are observed up to the entire scale of the SMC.