Context. The [CII] 158 {\mu}m far-infrared (FIR) fine-structure line is one of the most important cooling lines of the star-forming interstellar medium (ISM). High spectral resolution observations have shown complex structures in the line profiles of the [CII] emission. Aims. Our aim is to determine whether the complex profiles observed in [^{12}CII] are due to individual velocity components along the line-of-sight or to self-absorption based on a comparison of the [^{12}CII] and isotopic [^{13}CII] line profiles. Methods. Deep integrations with the SOFIA/upGREAT 7-pixel array receiver in M43, Horsehead~PDR, Monoceros~R2, and M17~SW allow for the detection of optically thin [^{13}CII] emission lines, along with the [^{12}CII] emission lines, with a high signal-to-noise ratio. We first derived the [^{12}CII] optical depth and the [CII] column density from a single component model. However, the complex line profiles observed require a double layer model with an emitting background and an absorbing foreground. A multi-component velocity fit allows us to derive the physical conditions of the [CII] gas: column density and excitation temperature. Results. We find moderate to high [^{12}CII] optical depths in all four sources and self-absorption of [^{12}CII] in Mon R2 and M17 SW. The high column density of the warm background emission corresponds to an equivalent Av of up to 41 mag. The foreground absorption requires substantial column densities of cold and dense [CII] gas, with an equivalent Av ranging up to about 13 mag. Conclusions. The column density of the warm background material requires multiple photon-dominated region (PDR) surfaces stacked along the line of sight and in velocity. The substantial column density of dense and cold foreground [CII] gas detected in absorption cannot be explained with any known scenario and we can only speculate on its origins

Molecular clouds form from the atomic phase of the interstellar medium. However, characterizing the transition between the atomic and the molecular interstellar medium (ISM) is a difficult observational task. Here we address cloud formation processes by combining HSIA with molecular line data. One scenario proposed by numerical simulations is that the column density probability density functions (N-PDF) evolves from a log-normal shape at early times to a power-law-like shape at later times. In this paper, we study the cold atomic component of the giant molecular filament GMF38a (d=3.4 kpc, length$\sim230$ pc). We identify an extended HISA feature, which is partly correlated with the 13CO emission. The peak velocities of the HISA and 13CO observations agree well on the eastern side of the filament, whereas a velocity offset of approximately 4 km s$^{-1}$ is found on the western side. The sonic Mach number we derive from the linewidth measurements shows that a large fraction of the HISA, which is ascribed to the cold neutral medium (CNM), is at subsonic and transonic velocities. The column density of the CNM is on the order of 10$^{20}$ to 10$^{21}$ cm$^{-2}$. The column density of molecular hydrogen is an order of magnitude higher. The N-PDFs from HISA (CNM), HI emission (WNM+CNM), and 13CO (molecular component) are well described by log-normal functions, which is in agreement with turbulent motions being the main driver of cloud dynamics. The N-PDF of the molecular component also shows a power law in the high column-density region, indicating self-gravity. We suggest that we are witnessing two different evolutionary stages within the filament. The eastern subregion seems to be forming a molecular cloud out of the atomic gas, whereas the western subregion already shows high column density peaks, active star formation and evidence of related feedback processes.

With the $Karl~G.~Jansky$ Very Large Array (VLA) in C-configuration, we observed a large portion of the first Galactic quadrant ($l=14.0-67.4^\circ $ and $\lvert b \rvert \leq 1.25^\circ $) achieving an angular resolution of $\leq 40^{\prime\prime}$. At $L$ Band, the WIDAR correlator at the VLA was set to cover the 21~cm HI line, four OH transitions, a series of H$n\alpha$ radio recombination lines (RRLs; $n=151$ to 186), and eight 128~MHz wide continuum spectral windows (SPWs) simultaneously. The HI emission shows clear filamentary substructures at negative velocities with low velocity crowding. The emission at positive velocities is more smeared-out likely due to higher spatial and velocity crowding of structures at the positive velocities. Comparing to the spiral arm model of the Milky Way, the atomic gas follows the Sagittarius and Perseus Arm well but with significant material in the inter-arm regions. With the C-configuration-only HI+continuum data, we produced a HI optical depth map of the THOR areal coverage from 228 absorption spectra with the nearest-neighbor method. With this $\tau$ map, we corrected the HI emission for optical depth and the derived column density is 38% higher than the column density with optically thin assumption. The total HI mass with optical depth correction in the survey region is 4.7$\times10^8~M_\odot$, 31% more than the mass derived assuming the emission is optically thin. If we apply this 31% correction to the whole Milky Way, the total atomic gas mass would be 9.4-10.5$\times 10^9~M_\odot$. Comparing the HI with existing CO data, we find a significant increase in the atomic-to-molecular gas ratio from the spiral arms to the inter-arm regions.

We present images of radio recombination lines (RRLs) at wavelengths around 17 cm from the star-forming region W49A to determine the kinematics of ionized gas in the THOR survey (The HI/OH/Recombination line survey of the inner Milky Way) at an angular resolution of 16.8"x13.8". The distribution of ionized gas appears to be affected by feedback processes from the star clusters in W49A. The velocity structure of the RRLs shows a complex behavior with respect to the molecular gas. We find a shell-like distribution of ionized gas as traced by RRL emission surrounding the central cluster of OB stars in W49A. We describe the evolution of the shell with the recent feedback model code WARPFIELD that includes the important physical processes and has previously been applied to the 30 Doradus region in the Large Magellanic Cloud. The cloud structure and dynamics of W49A are in agreement with a feedback-driven shell that is re-collapsing. The shell may have triggered star formation in other parts of W49A. We suggest that W49A is a potential candidate for star formation regulated by feedback-driven and re-collapsing shells.

We have carried out the largest and most unbiased search for hypercompact (HC) HII regions. Our method combines four interferometric radio continuum surveys (THOR, CORNISH, MAGPIS and White2005) with far-infrared and sub-mm Galactic Plane surveys to identify embedded HII regions with positive spectral indices. 120 positive spectrum HII regions have been identified from a total sample of 534 positive spectral index radio sources. None of these HII regions, including the known HCHII regions recovered in our search, fulfills the canonical definition of an HCHII region at 5 GHz. We suggest that the current canonical definition of HCHII regions is not accurate and should be revised to include a hierarchical structure of ionized gas that results in an extended morphology at 5 GHz. Correlating our search with known ultracompact (UC) HII region surveys, we find that roughly half of detected UCHII regions have positive spectral indices, instead of more commonly assumed flat and optically thin spectra. This implies a mix of optically thin and thick emission and has important implications for previous analyses which have so far assumed optically thin emission for these objects. Positive spectrum HII regions are statistically more luminous and possess higher Lyman continuum fluxes than HII regions with flat or negative indices. Positive spectrum HII regions are thus more likely to be associated with more luminous and massive stars. No differences are found in clump mass, linear diameter or luminosity-to-mass ratio between positive spectrum and non-positive spectrum HII regions.

The hydroxyl radical (OH) is present in the diffuse molecular and partially atomic phases of the interstellar medium (ISM), but its abundance relative to hydrogen is not clear. We aim to evaluate the abundance of OH with respect to molecular hydrogen using OH absorption against cm-continuum sources over the first Galactic quadrant. This OH study is part of the HI/OH/Recombination line survey (THOR). THOR is a Karl G. Jansky Very Large Array large program of atomic, molecular and ionized gas in the range 15{\deg}$\leq$l$\leq$67{\deg} and |b|$\leq$1{\deg}. It is the highest-resolution unbiased OH absorption survey to date towards this region. We combine the derived optical depths with literature 13CO(1-0) and HI observations to determine the OH abundance. We detect absorption in the 1665 and 1667 MHz transitions for continuum sources stronger than $F_{\rm cont}\geq$0.1 Jy/beam. OH absorption is found against $\sim$15% of these continuum sources with increasing fractions for stronger sources. Most of the absorption is associated with Galactic HII regions. We find OH and 13CO gas to have similar kinematic properties. The OH abundance decreases with increasing hydrogen column density. The OH abundance with respect to the total hydrogen nuclei column density (atomic and molecular phase) is in agreement with a constant abundance for $A_V$ < 10-20. Towards the lowest column densities, we find sources that exhibit OH absorption but no 13CO emission, indicating that OH is a well suited tracer of the low column density molecular gas. We present spatially resolved OH absorption towards W43. The unbiased nature of the THOR survey opens a new window onto the gas properties of the ISM. The characterization of the OH abundance over a large range of hydrogen gas column densities contributes to the understanding of OH as a molecular gas tracer and provides a starting point for future investigations.

Radio continuum surveys of the Galactic plane can find and characterize HII regions, supernova remnants (SNRs), planetary nebulae (PNe), and extragalactic sources. A number of surveys at high angular resolution (<25") at different wavelengths exist to study the interstellar medium (ISM), but no comparable high-resolution and high-sensitivity survey exists at long radio wavelengths around 21cm. We observed a large fraction of the Galactic plane in the first quadrant of the Milky Way (l=14.0-67.4deg and |b| < 1.25deg) with the Karl G. Jansky Very Large Array (VLA) in the C-configuration covering six continuum spectral windows. These data provide a detailed view on the compact as well as extended radio emission of our Galaxy and thousands of extragalactic background sources. We used the BLOBCAT software and extracted 10916 sources. After removing spurious source detections caused by the sidelobes of the synthesised beam, we classified 10387 sources as reliable detections. We smoothed the images to a common resolution of 25" and extracted the peak flux density of each source in each spectral window (SPW) to determine the spectral indices $\alpha$ (assuming $I(\nu)\propto\nu^\alpha$). By cross-matching with catalogs of HII regions, SNRs, PNe, and pulsars, we found radio counterparts for 840 HII regions, 52 SNRs, 164 PNe, and 38 pulsars. We found 79 continuum sources that are associated with X-ray sources. We identified 699 ultra-steep spectral sources ($\alpha < -1.3$) that could be high-redshift galaxies. Around 9000 of the sources we extracted are not classified specifically, but based on their spatial and spectral distribution, a large fraction of them is likely to be extragalactic background sources. More than 7750 sources do not have counterparts in the SIMBAD database, and more than 3760 sources do not have counterparts in the NED database.

There is a considerable deficiency in the number of known supernova remnants (SNRs) in the Galaxy compared to that expected. Searches for extended low-surface brightness radio sources may find new Galactic SNRs, but confusion with the much larger population of HII regions makes identifying such features challenging. SNRs can, however, be separated from HII regions using their significantly lower mid-infrared (MIR) to radio continuum intensity ratios. We use the combination of high-resolution 1-2 GHz continuum data from The HI, OH, Recombination line survey of the Milky Way (THOR) and lower-resolution VLA 1.4 GHz Galactic Plane Survey (VGPS) continuum data, together with MIR data from the Spitzer GLIMPSE, Spitzer MIPSGAL, and WISE surveys to identify SNR candidates. To ensure that the candidates are not being confused with HII regions, we exclude radio continuum sources from the WISE Catalog of Galactic HII Regions, which contains all known and candidate H II regions in the Galaxy. We locate 76 new Galactic SNR candidates in the THOR and VGPS combined survey area of 67.4deg>l>17.5deg, |b|<1.25deg and measure the radio flux density for 52 previously-known SNRs. The candidate SNRs have a similar spatial distribution to the known SNRs, although we note a large number of new candidates near l=30deg, the tangent point of the Scutum spiral arm. The candidates are on average smaller in angle compared to the known regions, 6.4'+/-4.7' versus 11.0'+/-7.8', and have lower integrated flux densities. If the 76 candidates are confirmed as true SNRs, for example using radio polarization measurements or by deriving radio spectral indices, this would more than double the number of known Galactic SNRs in the survey area. This large increase would still, however, leave a discrepancy between the known and expected SNR populations of about a factor of two.

The process of atomic-to-molecular (HI-to-H$_2$) gas conversion is fundamental for molecular-cloud formation and star formation. 21 cm observations of the star-forming region W43 revealed extremely high HI column densities, of 120-180 M$_{\odot}$ pc$^{-2}$, a factor of 10-20 larger than predicted by HI-to-H$_2$ transition theories. We analyze the observed HI with an HI-to-H$_2$ transition theoretical model, and show that the theory-observation discrepancy cannot be explained by the intense radiation in W43, nor by variations of the assumed volume density or H$_2$ formation-rate coefficient. We show that the large observed HI columns are naturally explained by several ($9-22$) HI-to-H$_2$ transition layers, superimposed along the sightlines of W43. We discuss other possible interpretations such as a non-steady-state scenario, and inefficient dust absorption. The case of W43 suggests that HI thresholds reported in extra-galactic observations are probably not associated with a single HI-to-H$_2$ transition, but are rather a result of several transition layers (clouds) along the sightlines, beam-diluted with diffuse inter-cloud gas.

Context: The past decade has witnessed a large number of Galactic plane surveys at angular resolutions below 20". However, no comparable high-resolution survey exists at long radio wavelengths around 21cm in line and continuum emission. Methods: Employing the Very Large Array (VLA) in the C-array configuration and a large program, we observe the HI 21cm line, four OH lines, nineteen Halpha radio recombination lines as well as the continuum emission from 1 to 2GHz in full polarization over a large part of the first Galactic quadrant. Results: Covering Galactic longitudes from 14.5 to 67.4deg and latitudes between +-1.25deg, we image all of these lines and the continuum at ~20" resolution. These data allow us to study the various components of the interstellar medium (ISM): from the atomic phase, traced by the HI line, to the molecular phase, observed by the OH transitions, to the ionized medium, revealed by the cm continuum and the Halpha radio recombination lines. Furthermore, the polarized continuum emission enables magnetic field studies. In this overview paper, we discuss the survey outline and present the first data release as well as early results from the different datasets. We now release the first half of the survey; the second half will follow later after the ongoing data processing has been completed. The data in fits format (continuum images and line data cubes) can be accessed through the project web-page http://www.mpia.de/thor. Conclusions: The HI/OH/Recombination line survey of the Milky Way (THOR) opens a new window to the different parts of the ISM. It enables detailed studies of molecular cloud formation, conversion of atomic to molecular gas, and feedback from HII regions as well as the magnetic field in the Milky Way. It is highly complementary to other surveys of our Galaxy, and comparing different datasets allows us to address many open questions.

We carried out a large program with the Karl G. Jansky Very Large Array (VLA): "THOR: The HI, OH, Recombination line survey of the Milky Way". We observed a significant portion of the Galactic plane in the first quadrant of the Milky Way in the 21cm HI line, 4 OH transitions, 19 radio recombination lines, and continuum from 1 to 2 GHz. In this paper we present a catalog of the continuum sources in the first half of the survey (l=14.0-37.9deg and l=47.1-51.2deg, |b|<1.1deg) at a spatial resolution of 10-25", with a spatially varying noise level of ~0.3-1 mJy/beam. The catalog contains ~4400 sources. Around 1200 of these are spatially resolved, and ~1000 are possible artifacts, given their low signal-to-noise ratios. Since the spatial distribution of the unresolved objects is evenly distributed and not confined to the Galactic plane, most of them are extragalactic. Thanks to the broad bandwidth of the observations from 1 to 2 GHz, we are able to determine a reliable spectral index for ~1800 sources. The spectral index distribution reveals a double-peaked profile with maxima at spectral indices of alpha = -1 and alpha = 0 , corresponding to steep declining and flat spectra, respectively. This allows us to distinguish between thermal and non-thermal emission, which can be used to determine the nature of each source. We examine the spectral index of ~300 known HII regions, for which we find thermal emission with spectral indices around alpha = 0. In contrast, supernova remnants (SNR) show non-thermal emission with alpha = -0.5 and extragalactic objects generally have a steeper spectral index of alpha = -1. Using the spectral index information of the THOR survey, we investigate potential SNR candidates. We classify the radiation of four SNR candidates as non-thermal, and for the first time, we provide strong evidence for the SNR origin of these candidates.

We report on observations of the hydroxyl radical (OH) within The H{\sc I}, OH Recombination line survey (THOR) pilot region. The region is bounded approximately between Galactic coordinates l=29.2 to 31.5$^\circ$ and b=-1.0 to +1.0$^\circ$ and includes the high-mass star forming region W43. We identify 103 maser sites, including 72 with 1612\,MHz masers, 42 showing masers in either of the main line transitions at 1665 and 1667\,MHz and four showing 1720\,MHz masers. Most maser sites with either main-line or 1720\,MHz emission are associated with star formation, whereas most of the 1612\,MHz masers are associated with evolved stars. We find that nearly all of the main-line maser sites are co-spatial with an infrared source, detected by GLIMPSE. We also find diffuse OH emission, as well as OH in absorption towards selected unresolved or partially resolved sites. Extended OH absorption is found towards the well known star forming complex W43 Main.

Aims: Understanding the fragmentation and collapse properties of the dense gas during the onset of high-mass star formation. Methods: We observed the massive (~800M_sun) starless gas clump IRDC18310-4 with the Plateau de Bure Interferometer (PdBI) at sub-arcsecond resolution in the 1.07mm continuum andN2H+(3-2) line emission. Results: Zooming from a single-dish low-resolution map to previous 3mm PdBI data, and now the new 1.07mm continuum observations, the sub-structures hierarchically fragment on the increasingly smaller spatial scales. While the fragment separations may still be roughly consistent with pure thermal Jeans fragmentation, the derived core masses are almost two orders of magnitude larger than the typical Jeans mass at the given densities and temperatures. However, the data can be reconciled with models using non-homogeneous initial density structures, turbulence and/or magnetic fields. While most sub-cores remain (far-)infrared dark even at 70mum, we identify weak 70mum emission toward one core with a comparably low luminosity of ~16L_sun, re-enforcing the general youth of the region. The spectral line data always exhibit multiple spectral components toward each core with comparably small line widths for the individual components (in the 0.3 to 1.0km/s regime). Based on single-dish C18O(2-1) data we estimate a low virial-to-gas-mass ratio <=0.25. We discuss that the likely origin of these spectral properties may be the global collapse of the original gas clump that results in multiple spectral components along each line of sight. Even within this dynamic picture the individual collapsing gas cores appear to have very low levels of internal turbulence.

To study the atomic, molecular and ionized emission of Giant Molecular Clouds (GMCs), we have initiated a Large Program with the VLA: 'THOR - The HI, OH, Recombination Line survey of the Milky Way'. We map the 21cm HI line, 4 OH lines, 19 H_alpha recombination lines and the continuum from 1 to 2 GHz of a significant fraction of the Milky Way (l=15-67deg, |b|<1deg) at ~20" resolution. In this paper, we focus on the HI emission from the W43 star-formation complex. Classically, the HI 21cm line is treated as optically thin with column densities calculated under this assumption. This might give reasonable results for regions of low-mass star-formation, however, it is not sufficient to describe GMCs. We analyzed strong continuum sources to measure the optical depth, and thus correct the HI 21cm emission for optical depth effects and weak diffuse continuum emission. Hence, we are able to measure the HI mass of W43 more accurately and our analysis reveals a lower limit of M=6.6x10^6 M_sun, which is a factor of 2.4 larger than the mass estimated with the assumption of optically thin emission. The HI column densities are as high as N(HI)~150 M_sun/pc^2 ~ 1.9x10^22 cm^-2, which is an order of magnitude higher than for low mass star formation regions. This result challenges theoretical models that predict a threshold for the HI column density of ~10 M_sun/pc^2, at which the formation of molecular hydrogen should set in. By assuming an elliptical layered structure for W43, we estimate the particle density profiles. While at the cloud edge atomic and molecular hydrogen are well mixed, the center of the cloud is strongly dominated by molecular hydrogen. We do not identify a sharp transition between hydrogen in atomic and molecular form. Our results are an important characterization of the atomic to molecular hydrogen transition in an extreme environment and challenges current theoretical models.

We want to understand the kinematic and thermal properties of young massive gas clumps prior to and at the earliest evolutionary stages of high-mass star formation. Do we find signatures of gravitational collapse? Do we find temperature gradients in the vicinity or absence of infrared emission sources? Do we find coherent velocity structures toward the center of the dense and cold gas clumps? To determine kinematics and gas temperatures, we used ammonia, because it is known to be a good tracer and thermometer of dense gas. We observed the NH$_3$(1,1) and (2,2) lines within seven very young high-mass star-forming regions with the VLA and the Effelsberg 100m telescope. This allows us to study velocity structures, linewidths, and gas temperatures at high spatial resolution of 3-5$"$, corresponding to $\sim$0.05 pc. We find on average cold gas clumps with temperatures in the range between 10 K and 30 K. The observations do not reveal a clear correlation between infrared emission peaks and ammonia temperature peaks. We report an upper limit for the linewidth of $\sim$1.3 km s$^{-1}$, at the spectral resolution limit of our VLA observation. This indicates a relatively low level of turbulence on the scale of the observations. Velocity gradients are present in almost all regions with typical velocity differences of 1 to 2 km s$^{-1}$ and gradients of 5 to 10 km s$^{-1}$ pc$^{-1}$. These velocity gradients are smooth in most cases, but there is one exceptional source (ISOSS23053), for which we find several velocity components with a steep velocity gradient toward the clump centers that is larger than 30 km s$^{-1}$ pc$^{-1}$. This steep velocity gradient is consistent with recent models of cloud collapse. Furthermore, we report a spatial correlation of ammonia and cold dust, but we also find decreasing ammonia emission close to infrared emission sources.

Understanding the chemical evolution of young (high-mass) star-forming regions is a central topic in star formation research. Chemistry is employed as a unique tool 1) to investigate the underlying physical processes and 2) to characterize the evolution of the chemical composition. We observed a sample of 59 high-mass star-forming regions at different evolutionary stages varying from the early starless phase of infrared dark clouds to high-mass protostellar objects to hot molecular cores and, finally, ultra-compact HII regions at 1mm and 3mm with the IRAM 30m telescope. We determined their large-scale chemical abundances and found that the chemical composition evolves along with the evolutionary stages. On average, the molecular abundances increase with time. We modeled the chemical evolution, using a 1D physical model where density and temperature vary from stage to stage coupled with an advanced gas-grain chemical model and derived the best-fit chi^2 values of all relevant parameters. A satisfying overall agreement between observed and modeled column densities for most of the molecules was obtained. With the best-fit model we also derived a chemical age for each stage, which gives the timescales for the transformation between two consecutive stages. The best-fit chemical ages are ~10,000 years for the IRDC stage, ~60,000 years for the HMPO stage, ~40,000 years for the HMC stage, and ~10,000 years for the UCHII stage. The total chemical timescale for the entire evolutionary sequence of the high-mass star formation process is on the order of 10^5 years, which is consistent with theoretical estimates. Furthermore, based on the approach of a multiple-line survey of unresolved data, we were able to constrain an intuitive and reasonable physical and chemical model. The results of this study can be used as chemical templates for the different evolutionary stages in high-mass star formation.

Aims: We study the fragmentation and dynamical properties of a massive starless gas clump at the onset of high-mass star formation. Methods: Based on Herschel continuum data we identify a massive gas clump that remains far-infrared dark up to 100mum wavelengths. The fragmentation and dynamical properties are investigated by means of Plateau de Bure Interferometer and Nobeyama 45m single-dish spectral line and continuum observations. Results: The massive gas reservoir fragments at spatial scales of ~18000AU in four cores. Comparing the spatial extent of this high-mass region with intermediate- to low-mass starless cores from the literature, we find that linear sizes do not vary significantly over the whole mass regime. However, the high-mass regions squeeze much more gas into these similar volumes and hence have orders of magnitude larger densities. The fragmentation properties of the presented low-to high-mass regions are consistent with gravitational instable Jeans fragmentation. Furthermore, we find multiple velocity components associated with the resolved cores. Recent radiative transfer hydrodynamic simulations of the dynamic collapse of massive gas clumps also result in multiple velocity components along the line of sight because of the clumpy structure of the regions. This result is supported by a ratio between viral and total gas mass for the whole region <1. Conclusions: This apparently still starless high-mass gas clump exhibits clear signatures of early fragmentation and dynamic collapse prior to the formation of an embedded heating source. A comparison with regions of lower mass reveals that the linear size of star-forming regions does not necessarily have to vary much for different masses, however, the mass reservoirs and gas densities are orders of magnitude enhanced for high-mass regions compared to their lower-mass siblings.