Search

The Milky Way's Central Molecular Zone (CMZ) differs dramatically from our local solar neighbourhood, both in the extreme interstellar medium conditions it exhibits (e.g. high gas, stellar, and feedback density) and in the strong dynamics at play (e.g. due to shear and gas influx along the bar). Consequently, it is likely that there are large-scale physical structures within the CMZ that cannot form elsewhere in the Milky Way. In this paper, we present new results from the Atacama Large Millimeter/submillimeter Array (ALMA) large programme ACES (ALMA CMZ Exploration Survey) and conduct a multi-wavelength and kinematic analysis to determine the origin of the M0.8$-$0.2 ring, a molecular cloud with a distinct ring-like morphology. We estimate the projected inner and outer radii of the M0.8$-$0.2 ring to be 79" and 154", respectively (3.1 pc and 6.1 pc at an assumed Galactic Centre distance of 8.2 kpc) and calculate a mean gas density $> 10^{4}$ cm$^{-3}$, a mass of $\sim$ $10^6$ M$_\odot$, and an expansion speed of $\sim$ 20 km s$^{-1}$, resulting in a high estimated kinetic energy ($> 10^{51}$ erg) and momentum ($> 10^7$ M$_\odot$ km s$^{-1}$). We discuss several possible causes for the existence and expansion of the structure, including stellar feedback and large-scale dynamics. We propose that the most likely cause of the M0.8$-$0.2 ring is a single high-energy hypernova explosion. To viably explain the observed morphology and kinematics, such an explosion would need to have taken place inside a dense, very massive molecular cloud, the remnants of which we now see as the M0.8$-$0.2 ring. In this case, the structure provides an extreme example of how supernovae can affect molecular clouds., Comment: 17 pages, 13 figures, and 2 tables. Accepted for publication in Astronomy & Astrophysics

We present the initial data for the ($J = 3 \to 2$) transition of $^{13}$CO obtained from the Central Molecular Zone (CMZ) of the Milky Way as part of the CO Heterodyne Inner Milky Way Plane Survey 2 (CHIMPS2). Covering $359^\circ \leq l \leq 1^\circ$ and $|b| \leq 0.5^\circ$ with an angular resolution of 19 arcsec, velocity resolution of 1 km s$^{-1}$, and rms $T_A^* = 0.59$ K at these resolutions, our observations unveil the complex structure of the CMZ molecular gas in improved detail. Complemented by the $^{12}$CO CHIMPS2 data, we estimate a median optical depth of $\tau_{13} = 0.087$. The preliminary analysis yields a median $^{13}$CO column density range equal to $N(^{13}\text{CO})= 2$--$5 \times 10^{18}$ cm$^{-2}$, median H$_2$ column density equal to $N(\text{H}_2)= 4 \times 10^{22}$ cm$^{-2}$ to $1 \times 10^{23}$ cm$^{-2}$. We derive $N(\text{H}_2)$-based total mass estimates of $M(\text{H}_2)= 2$--$6 \times 10^7\, M_{\odot}$, in agreement with previous studies. We analyze the relationship between the integrated intensity of $^{13}$CO and the surface density of compact sources identified by Herschel Hi-GAL, and find that younger Hi-GAL sources detected at 500 $\mu$m but not at 70 $\mu$m follow the dense gas of the CMZ more closely than those that are bright at 70 $\mu$m. The latter, actively star-forming sources, appear to be more associated with material in the foreground spiral arms., Comment: 13 pages, 9 figures

The interstellar medium in the Milky Way's Central Molecular Zone (CMZ) is known to be strongly magnetised, but its large-scale morphology and impact on the gas dynamics are not well understood. We explore the impact and properties of magnetic fields in the CMZ using three-dimensional non-self gravitating magnetohydrodynamical simulations of gas flow in an external Milky Way barred potential. We find that: (1) The magnetic field is conveniently decomposed into a regular time-averaged component and an irregular turbulent component. The regular component aligns well with the velocity vectors of the gas everywhere, including within the bar lanes. (2) The field geometry transitions from parallel to the Galactic plane near $z=0$ to poloidal away from the plane. (3) The magneto-rotational instability (MRI) causes an in-plane inflow of matter from the CMZ gas ring towards the central few parsecs of $0.01-0.1$ M$_\odot$ yr$^{-1}$ that is absent in the unmagnetised simulations. However, the magnetic fields have no significant effect on the larger-scale bar-driven inflow that brings the gas from the Galactic disc into the CMZ. (4) A combination of bar inflow and MRI-driven turbulence can sustain a turbulent vertical velocity dispersion of $\sigma_z \simeq 5$ km s$^{-1}$ on scales of $20$ pc in the CMZ ring. The MRI alone sustains a velocity dispersion of $\sigma_z \simeq 3$ km s$^{-1}$. Both these numbers are lower than the observed velocity dispersion of gas in the CMZ, suggesting that other processes such as stellar feedback are necessary to explain the observations. (5) Dynamo action driven by differential rotation and the MRI amplifies the magnetic fields in the CMZ ring until they saturate at a value that scales with the average local density as $B \simeq 102 (n/10^3 cm^{-3})^{0.33}$ $\mu$G. Finally, we discuss the implications of our results within the observational context in the CMZ., Comment: submitted to A&A, comments welcome

The fragmentation mode of high-mass molecular clumps and the accretion processes that form the most massive stars ($M\gtrsim 8M_\odot$) are still not well understood. To this end, we have undertaken a large observational program (CORE) making use of interferometric observations from the Northern Extended Millimetre Array (NOEMA) for a sample of 20 luminous ($L>10^4L_\odot$) protostellar objects in the 1.37 mm wavelength regime in both continuum and line emission, reaching $\sim$0.4" resolution (800 au at 2 kpc). Using the dense gas tracer CH$_3$CN, we find velocity gradients across 13 cores perpendicular to the directions of bipolar molecular outflows, making them excellent disk candidates. Specific angular momentum ($j$) radial profiles are on average $\sim10^{-3}$ km /s pc and follow $j \propto r^{1.7}$, consistent with a poorly resolved rotating and infalling envelope/disk model. Fitting the velocity profiles with a Keplerian model, we find protostellar masses in the range of $\sim 10-25$ $M_\odot$. Modelling the level population of CH$_3$CN lines, we present temperature maps and find median gas temperatures in the range $70-210$ K. We create Toomre $Q$ maps to study the stability of the disks and find almost all (11 of 13) disk candidates to be prone to fragmentation due to gravitational instabilities at the scales probed by our observations. In particular, disks with masses greater than $\sim10-20\%$ of the mass of their host (proto)stars are Toomre unstable, and more luminous protostellar objects tend to have disks that are more massive and hence more prone to fragmentation. Our finings show that most disks around high-mass protostars are prone to disk fragmentation early in their formation due to their high disk to stellar mass ratio. This impacts the accretion evolution of high-mass protostars which will have significant implications for the formation of the most massive stars., Comment: 27 pages, 12 figures, 6 appendices - accepted for publication in Astronomy and Astrophysics

Context: Core accretion models of massive star formation require the existence of massive, starless cores within molecular clouds. Yet, only a small number of candidates for such truly massive, monolithic cores are currently known. Aims: Here we analyse a massive core in the well-studied infrared-dark cloud (IRDC) called the 'dragon cloud' (also known as G028.37+00.07 or 'Cloud C'). This core (C2c1) sits at the end of a chain of a roughly equally spaced actively star-forming cores near the centre of the IRDC. Methods: We present new high-angular resolution 1 mm ALMA dust continuum and molecular line observations of the massive core. Results: The high-angular resolution observations show that this region fragments into two cores C2c1a and C2c1b, which retain significant background-subtracted masses of 23 Msun and 2 Msun (31 Msun and 6 Msun without background subtraction), respectively. The cores do not appear to fragment further on the scales of our highest angular resolution images (0.200 arcsec, 0.005 pc ~ 1000 AU). We find that these cores are very dense (nH2 > 10^6 cm-3) and have only trans-sonic non-thermal motions (Ms ~ 1). Together the mass, density and internal motions imply a virial parameter of < 1, which suggests the cores are gravitationally unstable, unless supported by strong magnetic fields with strengths of ~ 1 - 10 mG. From CO line observations, we find that there is tentative evidence for a weak molecular outflow towards the lower-mass core, and yet the more massive core remains devoid of any star formation indicators. Conclusions: We present evidence for the existence of a massive, pre-stellar core, which has implications for theories of massive star formation. This source warrants follow-up higher-angular-resolution observations to further assess its monolithic and pre-stellar nature., Comment: 8+4 pages, 4+2 Figures, 2 Tables, accepted for publication in Astronomy & Astrophysics

We present the spectroscopic confirmation of the brightest known gravitationally lensed Lyman break galaxy in the Epoch of Reionisation, A1703-zD1, through the detection of [C II] at a redshift of z = 6.8269 +/- 0.0004. This source was selected behind the strong lensing cluster Abell 1703, with an intrinsic L$_{UV}$ ~ L$^*$$_{z=7}$ luminosity and a very blue Spitzer/IRAC [3.6]-[4.5] colour, implying high equivalent width line emission of [O III]+H$\beta$. [C II] is reliably detected at 6.1$\sigma$ co-spatial with the rest-frame UV counterpart, showing similar spatial extent. Correcting for the lensing magnification, the [C II] luminosity in A1703-zD1 is broadly consistent with the local L$_{[CII]}$ - SFR relation. We find a clear velocity gradient of 103 +/- 22 km/s across the source which possibly indicates rotation or an ongoing merger. We furthermore present spectral scans with no detected [C II] above 4.6$\sigma$ in two unlensed Lyman break galaxies in the EGS-CANDELS field at z ~ 6.6 - 6.9. This is the first time that NOEMA has been successfully used to observe [C II] in a 'normal' star-forming galaxy at z > 6, and our results demonstrate its capability to complement ALMA in confirming galaxies in the Epoch of Reionisation., Comment: 9 pages, 15 figures. Accepted for publication in MNRAS

The atomic phase of the interstellar medium plays a key role in the formation process of molecular clouds. Due to the line-of-sight confusion in the Galactic plane that is associated with its ubiquity, atomic hydrogen emission has been challenging to study. Employing the high-angular resolution data from the THOR survey, we identify one of the largest, coherent, mostly atomic HI filaments in the Milky Way at the line-of-sight velocities around -54 km/s. The giant atomic filament "Maggie", with a total length of 1.2 kpc, is not detected in most other tracers, and does not show signs of active star formation. At a kinematic distance of 17 kpc, Maggie is situated below (by 500 pc) but parallel to the Galactic HI disk and is trailing the predicted location of the Outer Arm by 5-10 km/s in longitude-velocity space. The centroid velocity exhibits a smooth gradient of less than $\pm$3 km/s /10 pc and a coherent structure to within $\pm$6 km/s. The line widths of 10 km/s along the spine of the filament are dominated by non-thermal effects. After correcting for optical depth effects, the mass of Maggie's dense spine is estimated to be $7.2\times10^5\,M_{\odot}$. The mean number density of the filament is 4$\rm\,cm^{-3}$, which is best explained by the filament being a mix of cold and warm neutral gas. In contrast to molecular filaments, the turbulent Mach number and velocity structure function suggest that Maggie is driven by transonic to moderately supersonic velocities that are likely associated with the Galactic potential rather than being subject to the effects of self-gravity or stellar feedback. The column density PDF displays a log-normal shape around a mean of $N_{\rm HI} = 4.8\times 10^{20}\rm\,cm^{-2}$, thus reflecting the absence of dominating effects of gravitational contraction., Comment: 19 pages, 17 figures, accepted for publication in A&A

G0.253+0.016, commonly referred to as "the Brick" and located within the Central Molecular Zone, is one of the densest ($\approx10^{3-4}$ cm$^{-3}$) molecular clouds in the Galaxy to lack signatures of widespread star formation. We set out to constrain the origins of an arc-shaped molecular line emission feature located within the cloud. We determine that the arc, centred on $\{l_{0},b_{0}\}=\{0.248^{\circ}, 0.18^{\circ}\}$, has a radius of $1.3$ pc and kinematics indicative of the presence of a shell expanding at $5.2^{+2.7}_{-1.9}$ km s$^{-1}$. Extended radio continuum emission fills the arc cavity and recombination line emission peaks at a similar velocity to the arc, implying that the molecular and ionised gas are physically related. The inferred Lyman continuum photon rate is $N_{\rm LyC}=10^{46.0}-10^{47.9}$ photons s$^{-1}$, consistent with a star of spectral type B1-O8.5, corresponding to a mass of $\approx12-20$ M$_{\odot}$. We explore two scenarios for the origin of the arc: i) a partial shell swept up by the wind of an interloper high-mass star; ii) a partial shell swept up by stellar feedback resulting from in-situ star formation. We favour the latter scenario, finding reasonable (factor of a few) agreement between its morphology, dynamics, and energetics and those predicted for an expanding bubble driven by the wind from a high-mass star. The immediate implication is that G0.253+0.016 may not be as quiescent as is commonly accepted. We speculate that the cloud may have produced a $\lesssim10^{3}$ M$_{\odot}$ star cluster $\gtrsim0.4$ Myr ago, and demonstrate that the high-extinction and stellar crowding observed towards G0.253+0.016 may help to obscure such a star cluster from detection., Comment: 17 pages, 7 figures. Accepted for publication in MNRAS (October 15, 2021)

The feedback from young stars (i.e. pre-supernova) is thought to play a crucial role in molecular cloud destruction. In this paper, we assess the feedback mechanisms acting within a sample of 5810 HII regions identified from the PHANGS-MUSE survey of 19 nearby ($<$ 20 Mpc) star-forming, main sequence spiral galaxies (log($M_\star$/M$_\odot$)= 9.4 $-$ 11). These optical spectroscopic maps are essential to constrain the physical properties of the HII regions, which we use to investigate their internal pressure terms. We estimate the photoionised gas ($P_\mathrm{therm}$), direct radiation ($P_\mathrm{rad}$), and mechanical wind pressure ($P_\mathrm{wind}$), which we compare to the confining pressure of their host environment ($P_\mathrm{de}$). The HII regions remain unresolved within our ${\sim}50{-}100$ pc resolution observations, so we place upper ($P_\mathrm{max}$) and lower ($P_\mathrm{min}$) limits on each of the pressures by using a minimum (i.e. clumpy structure) and maximum (i.e. smooth structure) size, respectively. We find that the $P_\mathrm{max}$ measurements are broadly similar, and for $P_\mathrm{min}$ the $P_\mathrm{therm}$ is mildly dominant. We find that the majority of HII regions are over-pressured, $P_\mathrm{tot}/P_\mathrm{de} = (P_\mathrm{therm}+P_\mathrm{wind}+P_\mathrm{rad})/P_\mathrm{de} > 1$, and expanding, yet there is a small sample of compact HII regions with $P_\mathrm{tot,max}/P_\mathrm{de} < 1$ ($\sim$1% of the sample). These mostly reside in galaxy centres ($R_\mathrm{gal}<1$kpc), or, specifically, environments of high gas surface density; log($\Sigma_\mathrm{gas}/\mathrm{M_\odot} \mathrm{pc}^{-2}$)$\sim$2.5 (measured on kpc-scales). Lastly, we compare to a sample of literature measurements for $P_\mathrm{therm}$ and $P_\mathrm{rad}$ to investigate how dominant pressure term transitions over around 5dex in spatial dynamic range and 10 dex in pressure., Comment: 28 pages total, 16+3 figures, and 3 tables. Accepted for publication in MNRAS

We aim to characterize the physical and chemical properties of fragmented cores during the earliest evolutionary stages in the very young star-forming regions ISOSS J22478+6357 and ISOSS J23053+5953. NOEMA 1.3 mm data are used in combination with archival mid- and far-infrared observations to construct and fit the SEDs of individual fragmented cores. The radial density profiles are inferred from the 1.3 mm continuum visibility profiles and the radial temperature profiles are estimated from H2CO rotation temperature maps. Molecular column densities are derived with the line fitting tool XCLASS. The physical and chemical properties are combined by applying the physical-chemical model MUSCLE in order to constrain the chemical timescales of a few line-rich cores. The morphology and spatial correlations of the molecular emission are analyzed using the HOG method. The mid-infrared data show that both regions contain a cluster of young stellar objects. Bipolar molecular outflows are observed in the CO 2-1 transition toward the strong mm cores indicating protostellar activity. We find strong molecular emission of SO, SiO, H2CO, and CH3OH in locations which are not associated with the mm cores. These shocked knots can be either associated with the bipolar outflows or, in the case of ISOSS J23053+5953, with a colliding flow that creates a large shocked region between the mm cores. The mean chemical timescale of the cores is lower (20 000 yr) compared to that of the sources of the more evolved CORE sample (60 000 yr). With the HOG method, we find that the spatial emission of species tracing the extended emission and of shock-tracing molecules are well correlated within transitions of these groups., Comment: 34 pages, 17 figures, accepted for publication in A&A

Increasing evidence suggests that, similar to their low-mass counterparts, high-mass stars form through a disk-mediated accretion process. At the same time, formation of high-mass stars still necessitates high accretion rates, and hence, high gas densities, which in turn can cause disks to become unstable against gravitational fragmentation. We study the kinematics and fragmentation of the disk around the high-mass star forming region AFGL 2591-VLA 3 which was hypothesized to be fragmenting based on the observations that show multiple outflow directions. We use a new set of high-resolution (0.19 arcsec) IRAM/NOEMA observations at 843 micron towards VLA 3 which allow us to resolve its disk, characterize the fragmentation, and study its kinematics. In addition to the 843 micron continuum emission, our spectral setup targets warm dense gas and outflow tracers such as HCN, HC$_3$N and SO$_2$, as well as vibrationally excited HCN lines. The high resolution continuum and line emission maps reveal multiple fragments with subsolar masses within the inner 1000 AU of VLA 3. Furthermore, the velocity field of the inner disk observed at 843 micron shows a similar behavior to that of the larger scale velocity field studied in the CORE project at 1.37 mm. We present the first observational evidence for disk fragmentation towards AFGL 2591-VLA 3, a source that was thought to be a single high-mass core. While the fragments themselves are low-mass, the rotation of the disk is dominated by the protostar with a mass of 10.3$\pm 1.8~M_{\odot}$. These data also show that NOEMA Band 4 can obtain the highest currently achievable spatial resolution at (sub-)mm wavelengths in observations of strong northern sources., Comment: 16 pages, 11 figures, Accepted for publication in A&A

Context: The formation of high-mass star-forming regions from their parental gas cloud and the subsequent fragmentation processes lie at the heart of star formation research. Aims: We aim to study the dynamical and fragmentation properties at very early evolutionary stages of high-mass star formation. Methods: Employing the NOrthern Extended Millimeter Array (NOEMA) and the IRAM 30m telescope, we observed two young high-mass star-forming regions, ISOSS22478 and ISOSS23053, in the 1.3mm continuum and spectral line emission at a high angular resolution (~0.8''). Results: We resolved 29 cores that are mostly located along filament-like structures. Depending on the temperature assumption, these cores follow a mass-size relation of approximately M~r^2.0, corresponding to constant mean column densities. However, with different temperature assumptions, a steeper mass-size relation up to M~r^3.0, which would be more likely to correspond to constant mean volume densities, cannot be ruled out. The correlation of the core masses with their nearest neighbor separations is consistent with thermal Jeans fragmentation. We found hardly any core separations at the spatial resolution limit, indicating that the data resolve the large-scale fragmentation well. Although the kinematics of the two regions appear very different at first sight - multiple velocity components along filaments in ISOSS22478 versus a steep velocity gradient of more than 50km/s/pc in ISOSS23053 - the findings can be explained within the framework of a dynamical cloud collapse scenario. Conclusions: While our data are consistent with a dynamical cloud collapse scenario and subsequent thermal Jeans fragmentation, the importance of additional environmental properties, such as the magnetization of the gas or external shocks triggering converging gas flows, is nonetheless not as well constrained and would require future investigation., Comment: 16 pages, 16 figures, accepted for Astronomy and Astrophysics, a higher-resolution version can also be found at https://www.mpia.de/homes/beuther/papers.html

Infrared dark clouds (IRDCs) are potential hosts of the elusive early phases of high-mass star formation (HMSF). Here we conduct an in-depth analysis of the fragmentation properties of a sample of 10 IRDCs, which have been highlighted as some of the best candidates to study HMSF within the Milky Way. To do so, we have obtained a set of large mosaics covering these IRDCs with ALMA at band 3 (or 3mm). These observations have a high angular resolution (~3arcsec or ~0.05pc), and high continuum and spectral line sensitivity (~0.15mJy/beam and ~0.2K per 0.1km/s channel at the N2H+(1-0) transition). From the dust continuum emission, we identify 96 cores ranging from low- to high-mass (M = 3.4 to 50.9Msun) that are gravitationally bound (alpha_vir = 0.3 to 1.3) and which would require magnetic field strengths of B = 0.3 to 1.0mG to be in virial equilibrium. We combine these results with a homogenised catalogue of literature cores to recover the hierarchical structure within these clouds over four orders of magnitude in spatial scale (0.01pc to 10pc). Using supplementary observations at an even higher angular resolution, we find that the smallest fragments (<0.02pc) within this hierarchy do not currently have the mass and/or the density required to form high-mass stars. Nonetheless, the new ALMA observations presented in this paper have facilitated the identification of 19 (6 quiescent and 13 star-forming) cores that retain >16Msun without further fragmentation. These high-mass cores contain trans-sonic non-thermal motions, are kinematically sub-virial, and require moderate magnetic field strengths for support against collapse. The identification of these potential sites of high-mass star formation represents a key step in allowing us to test the predictions from high-mass star and cluster formation theories., Comment: 21(+6 in appendix) pages, 12(+3) figures, 3(+3) tables. Machine-readable versions of Table A1, A2 and A3 are available online. Accepted for publication in MNRAS

For centuries researchers have used sound to monitor and study wildlife. Traditionally, conservationists have identified species by ear; however, it is now common to deploy audio recording technology to monitor animal and ecosystem sounds. Animals use sound for communication, mating, navigation and territorial defence. Animal sounds provide valuable information and help conservationists to quantify biodiversity. Acoustic monitoring has grown in popularity due to the availability of diverse sensor types which include camera traps, portable acoustic sensors, passive acoustic sensors, and even smartphones. Passive acoustic sensors are easy to deploy and can be left running for long durations to provide insights on habitat and the sounds made by animals and illegal activity. While this technology brings enormous benefits, the amount of data that is generated makes processing a time-consuming process for conservationists. Consequently, there is interest among conservationists to automatically process acoustic data to help speed up biodiversity assessments. Processing these large data sources and extracting relevant sounds from background noise introduces significant challenges. In this paper we outline an approach for achieving this using state of the art in machine learning to automatically extract features from time-series audio signals and modelling deep learning models to classify different bird species based on the sounds they make. The acquired bird songs are processed using mel-frequency cepstrum (MFC) to extract features which are later classified using a multilayer perceptron (MLP). Our proposed method achieved promising results with 0.74 sensitivity, 0.92 specificity and an accuracy of 0.74.

We use sub-arcsecond resolution ($\sim$0.4$''$) observations with NOEMA at 1.37 mm to study the dust emission and molecular gas of 18 high-mass star-forming regions. We combine the derived physical and chemical properties of individual cores in these regions to estimate their ages. The temperature structure of these regions are determined by fitting H2CO and CH3CN line emission. The density profiles are inferred from the 1.37 mm continuum visibilities. The column densities of 11 different species are determined by fitting the emission lines with XCLASS. Within the 18 observed regions, we identify 22 individual cores with associated 1.37 mm continuum emission and with a radially decreasing temperature profile. We find an average temperature power-law index of q = 0.4$\pm$0.1 and an average density power-law index of p = 2.0$\pm$0.2 on scales on the order of several 1 000 au. Comparing these results with values of p derived in the literature suggest that the density profiles remain unchanged from clump to core scales. The column densities relative to N(C18O) between pairs of dense gas tracers show tight correlations. We apply the physical-chemical model MUSCLE to the derived column densities of each core and find a mean chemical age of $\sim$60 000 yrs and an age spread of 20 000-100 000 yrs. With this paper we release all data products of the CORE project available at https://www.mpia.de/core. The CORE sample reveals well constrained density and temperature power-law distributions. Furthermore, we characterize a large variety in molecular richness that can be explained by an age spread confirmed by our physical-chemical modeling. The hot molecular cores show the most emission lines, but we also find evolved cores at an evolutionary stage, in which most molecules are destroyed and thus the spectra appear line-poor again., Comment: 29 pages, 14 figures, accepted for publication in A&A

In the massive star-forming region IRAS 21078+5211, a highly fragmented cluster (0.1~pc in size) of molecular cores is observed, located at the density peak of an elongated (1~pc in size) molecular cloud. A small (1~km/s per 0.1~pc) LSR velocity (Vlsr) gradient is detected across the axis of the molecular cloud. Assuming we are observing a mass flow from the harboring cloud to the cluster, we derive a mass infall rate of about 10^{-4}~M_{sun}~yr^{-1}. The most massive cores (labeled 1, 2, and 3) are found at the center of the cluster, and these are the only ones that present a signature of protostellar activity in terms of emission from high-excitation molecular lines or a molecular outflow. We reveal an extended (size about 0.1~pc), bipolar collimated molecular outflow emerging from core 1. We believe this is powered by a (previously discovered) compact (size <= 1000~au) radio jet, ejected by a YSO embedded in core 1 (named YSO-1), since the molecular outflow and the radio jet are almost parallel and have a comparable momentum rate. By means of high-excitation lines, we find a large (14~km/s over 500~au) Vlsr gradient at the position of YSO-1, oriented approximately perpendicular to the radio jet. Assuming this is an edge-on, rotating disk and fitting a Keplerian rotation pattern, we determine the YSO-1 mass to be 5.6+/-2.0~M_{sun}. The water masers (previously observed with VLBI) emerge within 100-300~au from YSO-1 and are unique tracers of the jet kinematics. Their three-dimensional (3D) velocity pattern reveals that the gas flows along, and rotates about, the jet axis. We show that the 3D maser velocities are fully consistent with the magneto-centrifugal disk-wind models predicting a cylindrical rotating jet. Under this hypothesis, we determine the jet radius to be about 16~au and the corresponding launching radius and terminal velocity to be about 2.2~au and 200~km/s, respectively., Comment: 25 pages, 16 figures

The latest generation of Galactic-plane surveys is enhancing our ability to study the effects of galactic environment upon the process of star formation. We present the first data from CO Heterodyne Inner Milky Way Plane Survey 2 (CHIMPS2). CHIMPS2 is a survey that will observe the Inner Galaxy, the Central Molecular Zone (CMZ), and a section of the Outer Galaxy in $^{12}$CO, $^{13}$CO, and C$^{18}$O $(J = 3\rightarrow2)$ emission with the Heterodyne Array Receiver Program on the James Clerk Maxwell Telescope (JCMT). The first CHIMPS2 data presented here are a first look towards the CMZ in $^{12}$CO J = 3$\rightarrow$2 and cover $-3^{\circ}\leq\,\ell\,\leq\,5^{\circ}$ and $\mid$b$\mid \leq 0.5^{\circ}$ with angular resolution of 15 arcsec, velocity resolution of 1 km s$^{-1}$, and rms $\Delta T_A ^\ast =$ 0.58 K at these resolutions. Such high-resolution observations of the CMZ will be a valuable data set for future studies, whilst complementing the existing Galactic Plane surveys, such as SEDIGISM, the Herschel infrared Galactic Plane Survey, and ATLASGAL. In this paper, we discuss the survey plan, the current observations and data, as well as presenting position-position maps of the region. The position-velocity maps detect foreground spiral arms in both absorption and emission., Comment: 18 pages, 14 figures, accepted for publication in MNRAS

The density structure of the interstellar medium (ISM) determines where stars form and release energy, momentum, and heavy elements, driving galaxy evolution. Density variations are seeded and amplified by gas motion, but the exact nature of this motion is unknown across spatial scale and galactic environment. Although dense star-forming gas likely emerges from a combination of instabilities, convergent flows, and turbulence, establishing the precise origin is challenging because it requires quantifying gas motion over many orders of magnitude in spatial scale. Here we measure the motion of molecular gas in the Milky Way and in nearby galaxy NGC 4321, assembling observations that span an unprecedented spatial dynamic range ($10^{-1}{-}10^3$ pc). We detect ubiquitous velocity fluctuations across all spatial scales and galactic environments. Statistical analysis of these fluctuations indicates how star-forming gas is assembled. We discover oscillatory gas flows with wavelengths ranging from $0.3{-}400$ pc. These flows are coupled to regularly-spaced density enhancements that likely form via gravitational instabilities. We also identify stochastic and scale-free velocity and density fluctuations, consistent with the structure generated in turbulent flows. Our results demonstrate that ISM structure cannot be considered in isolation. Instead, its formation and evolution is controlled by nested, interdependent flows of matter covering many orders of magnitude in spatial scale., Comment: Published in Nature Astronomy on July 6th 2020. This is the authors' version before final edits. Includes methods and supplementary information. Link to the NA publication: https://www.nature.com/articles/s41550-020-1126-z

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., Comment: 22 pages, 25 figures, accepted for publication by A&A

We present a detailed observational and modeling study of the hot core VLA 3 in the high-mass star-forming region AFGL 2591, which is a target region of the NOrthern Extended Millimeter Array (NOEMA) large program CORE. Using NOEMA observations at 1.37 mm with an angular resolution of ~0."42 (1 400 au at 3.33 kpc), we derived the physical and chemical structure of the source. We modeled the observed molecular abundances with the chemical evolution code MUSCLE (MUlti Stage ChemicaL codE). Results. With the kinetic temperature tracers CH3CN and H2CO we observe a temperature distribution with a power-law index of q = 0.41+-0.08. Using the visibilities of the continuum emission we derive a density structure with a power-law index of p = 1.7+-0.1. The hot core spectra reveal high molecular abundances and a rich diversity in complex molecules. The majority of the molecules have an asymmetric spatial distribution around the forming protostar(s), which indicates a complex physical structure on scales < 1 400 au. Using MUSCLE, we are able to explain the observed molecular abundance of 10 out of 14 modeled species at an estimated hot core chemical age of ~21 100 years. In contrast to the observational analysis, our chemical modeling predicts a lower density power-law index of p < 1.4. Reasons for this discrepancy are discussed. Conclusions. Combining high spatial resolution observations with detailed chemical modeling allows us to derive a concise picture of the physical and chemical structure of the famous AFGL 2591 hot core. The next steps are to conduct a similar analysis for the whole CORE sample, and then use this analysis to constrain the chemical diversity in high-mass star formation to a much greater depth., Comment: 22 pages, 13 figures, accepted for publication in A&A

Context: OH masers trace diverse physical processes, from the expanding envelopes around evolved stars to star-forming regions or supernovae remnants. Aims: We identify the ground-state OH masers at 18cm wavelength in the area covered by ``The HI/OH/Recombination line survey of the Milky Way (THOR)''. We present a catalogue of all OH maser features and their possible associated environments. Methods: The THOR survey covers longitude and latitude ranges of 14.3

The formation process of high-mass stars (>8M$_\odot$) is poorly constrained, particularly, the effects of clump fragmentation creating multiple systems and the mechanism of mass accretion onto the cores. We study the fragmentation of dense gas clumps, and trace the circumstellar rotation and outflows by analyzing observations of the high-mass (~500M$_\odot$) star-forming region IRAS 23033+5951. Using the Northern Extended Millimeter Array (NOEMA) in three configurations and the IRAM 30-m single-dish telescope at 220GHz, we probe the gas and dust emission at an angular resolution of ~0.45arcsec, corresponding to 1900au. In the mm continuum emission, we identify a protostellar cluster with at least four mm-sources, where three of them show a significantly higher peak intensity well above a signal-to-noise ratio of 100. Hierarchical fragmentation from large to small spatial scales is discussed. Two fragments are embedded in rotating structures and drive molecular outflows, traced by $^{13}$CO (2-1) emission. The velocity profiles across two of the cores are similar to Keplerian but are missing the highest velocity components close to the center of rotation, which is a common phenomena from observations like these, and other rotation scenarios are not excluded entirely. Position-velocity diagrams suggest protostellar masses of ~6 and 19M$_\sun$. Rotational temperatures from fitting CH$_3$CN ($12_K-11_K$) spectra are used for estimating the gas temperature and by that the disk stability against gravitational fragmentation, utilizing Toomre's $Q$ parameter. [We] identify only one candidate disk to be unstable against gravitational instability caused by axisymmetric perturbations. The dominant sources cover different evolutionary stages within the same maternal gas clump. The appearance of rotation and outflows of the cores are similar to those found in low-mass star-forming regions., Comment: 23 pages, 16 figures, 5 tables, accepted for publication in A&A

This study is part of the project ``CORE'', an IRAM/NOEMA large program consisting of observations of the millimeter continuum and molecular line emission towards 20 selected high-mass star forming regions. We focus on IRAS23385+6053, which is believed to be the least evolved source of the CORE sample. The observations were performed at ~1.4 mm and employed three configurations of NOEMA and additional single-dish maps, merged with the interferometric data to recover the extended emission. Our correlator setup covered a number of lines from well-known hot core tracers and a few outflow tracers. The angular (~0.45"$-$0.9") and spectral (0.5 km/s) resolutions were sufficient to resolve the clump in IRAS23385+6053 and investigate the existence of large-scale motions due to rotation, infall, or expansion. We find that the clump splits into six distinct cores when observed at sub-arcsecond resolution. These are identified through their 1.4 mm continuum and molecular line emission. We produce maps of the velocity, line width, and rotational temperature from the methanol and methyl cyanide lines, which allow us to investigate the cores and reveal a velocity and temperature gradient in the most massive core. We also find evidence of a bipolar outflow, possibly powered by a low-mass star. We present the tentative detection of a circumstellar self-gravitating disk lying in the most massive core and powering a large-scale outflow previously known in the literature. In our scenario, the star powering the flow is responsible for most of the luminosity of IRAS23385+6053 (~$3000~L_\odot$). The other cores, albeit with masses below the corresponding virial masses, appear to be accreting material from their molecular surroundings and are possibly collapsing or on the verge of collapse. We conclude that we are observing a sample of star-forming cores that is bound to turn into a cluster of massive stars.

As part of a series of papers aimed at understanding the evolution of the Milky Way's Central Molecular Zone (CMZ), we present hydrodynamical simulations of turbulent molecular clouds orbiting in an accurate model of the gravitational potential extant there. We consider two sets of model clouds differing in the energy content of their velocity fields. In the first, self--virialised set, the turbulent kinetic energies are chosen to be close in magnitude to the clouds' self--gravitational potential energies. Comparison with isolated clouds evolving without an external potential shows that the self--virialised clouds are unable to withstand the compressive tidal field of the CMZ and rapidly collapse, forming stars much faster and reaching gas exhaustion after a small fraction of a Galactocentric orbit. In the second, tidally--virialised, set of simulations, the clouds' turbulent kinetic energies are in equilibrium with the external tidal field. These models are better supported against the field and the stronger turbulence suppresses star formation. Our results strongly support the inference that anomalously low star formation rates in the CMZ are due primarily to high velocity dispersions in the molecular gas. The clouds follow open, eccentric orbits oscillating in all three spatial coordinates. We examine the consequences of the orbital dynamics, particularly pericentre passage, by performing companion simulations of clouds on circular orbits. The increased tidal forces at pericentre produce transient accelerations in star formation rates of at most a factor of 2.7. Our results demonstrate that modelling star formation in galactic centres requires the inclusion of tidal forces., Comment: 21 pages, 18 figures, accepted by MNRAS

Young massive clusters (YMCs) are the most compact, high-mass stellar systems still forming at the present day. The precursor clouds to such systems are, however, rare due to their large initial gas mass reservoirs and rapid dispersal timescales due to stellar feedback. Nonetheless, unlike their high-z counterparts, these precursors are resolvable down to the sites of individually forming stars, and hence represent the ideal environments in which to test the current theories of star and cluster formation. Using high angular resolution (1$^{\prime\prime}$ / 0.05pc) and sensitivity ALMA observations of two YMC progenitor clouds in the Galactic Centre, we have identified a suite of molecular line transitions -- e.g. c-C$_{3}$H$_{2} $($7-6$) -- that are believed to be optically thin, and reliably trace the gas structure in the highest density gas on star-forming core scales. We conduct a virial analysis of the identified core and proto-cluster regions, and show that half of the cores (5/10) and both proto-clusters are unstable to gravitational collapse. This is the first kinematic evidence of global gravitational collapse in YMC precursor clouds at such an early evolutionary stage. The implications are that if these clouds are to form YMCs, then they likely do so via the "conveyor-belt" mode, whereby stars continually form within dispersed dense gas cores as the cloud undergoes global gravitational collapse. The concurrent contraction of both the cluster-scale gas and embedded (proto)stars ultimately leads to the high (proto)stellar density in YMCs., Comment: Accepted for publication in MNRAS. 22 pages (+22 pages of appendix), 10 (+21 in appendix) figures, 5 (+1 in appendix) tables

In this paper we provide a comprehensive description of the internal dynamics of G0.253+0.016 (a.k.a. 'the Brick'); one of the most massive and dense molecular clouds in the Galaxy to lack signatures of widespread star formation. As a potential host to a future generation of high-mass stars, understanding largely quiescent molecular clouds like G0.253+0.016 is of critical importance. In this paper, we reanalyse Atacama Large Millimeter Array cycle 0 HNCO $J=4(0,4)-3(0,3)$ data at 3 mm, using two new pieces of software which we make available to the community. First, scousepy, a Python implementation of the spectral line fitting algorithm scouse. Secondly, acorns (Agglomerative Clustering for ORganising Nested Structures), a hierarchical n-dimensional clustering algorithm designed for use with discrete spectroscopic data. Together, these tools provide an unbiased measurement of the line of sight velocity dispersion in this cloud, $\sigma_{v_{los}, {\rm 1D}}=4.4\pm2.1$ kms$^{-1}$, which is somewhat larger than predicted by velocity dispersion-size relations for the Central Molecular Zone (CMZ). The dispersion of centroid velocities in the plane of the sky are comparable, yielding $\sigma_{v_{los}, {\rm 1D}}/\sigma_{v_{pos}, {\rm 1D}}\sim1.2\pm0.3$. This isotropy may indicate that the line-of-sight extent of the cloud is approximately equivalent to that in the plane of the sky. Combining our kinematic decomposition with radiative transfer modelling we conclude that G0.253+0.016 is not a single, coherent, and centrally-condensed molecular cloud; 'the Brick' is not a \emph{brick}. Instead, G0.253+0.016 is a dynamically complex and hierarchically-structured molecular cloud whose morphology is consistent with the influence of the orbital dynamics and shear in the CMZ., Comment: 29 pages, 19 figures, 1 table (including appendix). Accepted for publication in MNRAS (February 4, 2019). Scousepy is available here: https://github.com/jdhenshaw/scousepy. Acorns is available here: https://github.com/jdhenshaw/acorns

The evolution of molecular clouds in galactic centres is thought to differ from that in galactic discs due to a significant influence of the external gravitational potential. We present a set of numerical simulations of molecular clouds orbiting on the 100-pc stream of the Central Molecular Zone (the central $\sim500$ pc of the Galaxy) and characterise their morphological and kinematic evolution in response to the background potential and eccentric orbital motion. We find that the clouds are shaped by strong shear and torques, by tidal and geometric deformation, and by their passage through the orbital pericentre. Within our simulations, these mechanisms control cloud sizes, aspect ratios, position angles, filamentary structure, column densities, velocity dispersions, line-of-sight velocity gradients, spin angular momenta, and kinematic complexity. By comparing these predictions to observations of clouds on the Galactic Centre 'dust ridge', we find that our simulations naturally reproduce a broad range of key observed morphological and kinematic features, which can be explained in terms of well-understood physical mechanisms. We argue that the accretion of gas clouds onto the central regions of galaxies, where the rotation curve turns over and the tidal field is fully compressive, is accompanied by transformative dynamical changes to the clouds, leading to collapse and star formation. This can generate an evolutionary progression of cloud collapse with a common starting point, which either marks the time of accretion onto the tidally-compressive region or of the most recent pericentre passage. Together, these processes may naturally produce the synchronised starbursts observed in numerous (extra)galactic nuclei., Comment: 21 pages, 8 figures, 3 tables; accepted by MNRAS (February 5, 2019); Figures 2, 3, and 4 show the main results of the paper. An animated version of Figure 2 is available as an ancillary file in the arXiv source and will be included as online Supporting Information with the MNRAS publication

Methods: We observed the high-mass hot core region G351.77-0.54 with ALMA and more than 16km baselines. Results: At a spatial resolution of 18/40au (depending on the distance), we identify twelve sub-structures within the inner few thousand au of the region. The brightness temperatures are high, reaching values greater 1000K, signposting high optical depth toward the peak positions. Core separations vary between sub-100au to several 100 and 1000au. The core separations and approximate masses are largely consistent with thermal Jeans fragmentation of a dense gas core. Due to the high continuum optical depth, most spectral lines are seen in absorption. However, a few exceptional emission lines are found that most likely stem from transitions with excitation conditions above1000K. Toward the main continuum source, these emission lines exhibit a velocity gradient across scales of 100-200au aligned with the molecular outflow and perpendicular to the previously inferred disk orientation. While we cannot exclude that these observational features stem from an inner hot accretion disk, the alignment with the outflow rather suggests that it stems from the inner jet and outflow region. The highest-velocity features are found toward the peak position, and no Hubble-like velocity structure can be identified. Therefore, these data are consistent with steady-state turbulent entrainment of the hot molecular gas via Kelvin-Helmholtz instabilities at the interface between the jet and the outflow. Conclusions: Resolving this high-mass star-forming region at sub-50au scales indicates that the hierarchical fragmentation process in the framework of thermal Jeans fragmentation can continue down to the smallest accessible spatial scales. Velocity gradients on these small scales have to be treated cautiously and do not necessarily stem from disks, but may be better explained with outflow emission., Comment: Accepted by Astronomy & Astrophysics, high-resolution version available at http://www.mpia.de/homes/beuther/papers.html

Using molecular line data from the Millimetre Astronomy Legacy Team 90 GHz Survey (MALT90), we have searched the optically thick \hcop\, line for the "blue asymmetry" spectroscopic signature of infall motion in a large sample of high-mass, dense molecular clumps observed to be at different evolutionary stages of star cluster formation according to their mid-infrared appearance. To quantify the degree of the line asymmetry, we measure the asymmetry parameter $A = {{I_{blue}-I_{red}}\over{I_{blue}+I_{red}}}$, the fraction of the integrated intensity that lies to the blueshifted side of the systemic velocity determined from the optically thin tracer \nthp. For a sample of 1,093 sources, both the mean and median of $A$ are positive ($A = 0.083\pm0.010$ and $0.065\pm0.009$, respectively) with high statistical significance, and a majority of sources (a fraction of $0.607 \pm 0.015$ of the sample) show positive values of A, indicating a preponderance of blue-asymmetric profiles over red-asymmetric profiles. Two other measures, the local slope of the line at the systemic velocity and the $\delta v$ parameter of \citet{Mardones1997}, also show an overall blue asymmetry for the sample, but with smaller statistical significance. This blue asymmetry indicates that these high-mass clumps are predominantly undergoing gravitational collapse. The blue asymmetry is larger ($A \sim 0.12$) for the earliest evolutionary stages (quiescent, protostellar and compact H II region) than for the later H II region ($A \sim 0.06$) and PDR ($A \sim 0$) classifications., Comment: 35 pages, 16 figures

We introduce the histogram of oriented gradients (HOG), a tool developed for machine vision that we propose as a new metric for the systematic characterization of observations of atomic and molecular gas and the study of molecular cloud formation models. In essence, the HOG technique takes as input extended spectral-line observations from two tracers and provides an estimate of their spatial correlation across velocity channels. We characterize HOG using synthetic observations of HI and $^{13}$CO(J=1-0) emission from numerical simulations of MHD turbulence leading to the formation of molecular gas after the collision of two atomic clouds. We find a significant spatial correlation between the two tracers in velocity channels where $v_{HI}\approx v_{^{13}CO}$, independent of the orientation of the collision with respect to the line of sight. We use HOG to investigate the spatial correlation of the HI, from the THOR survey, and the $^{13}$CO(J=1-0) emission, from the GRS, toward the portion of the Galactic plane 33.75$\lt l\lt$35.25$^{o}$ and $|b|\lt$1.25$^{o}$. We find a significant spatial correlation between the tracers in extended portions of the studied region. Although some of the regions with high spatial correlation are associated with HI self-absorption features, suggesting that it is produced by the cold atomic gas, the correlation is not exclusive to this kind of region. The HOG results also indicate significant differences between individual regions: some show spatial correlation in channels around $v_{HI}\approx v_{^{13}CO}$ while others present this correlation in velocity channels separated by a few km/s. We associate these velocity offsets to the effect of feedback and to the presence of physical conditions that are not included in the atomic-cloud-collision simulations, such as more general magnetic field configurations, shear, and global gas infall., Comment: 32 pages, 36 figures. Accepted for publication at A&A (28DEC2018)

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., Comment: 10 pages, 21 figures. Accepted to be published in Astronomy and Astrophysics. The full continuum catalog, full table of Table 4 and all the fits files of the continuum data are available at the project website http://www2.mpia-hd.mpg.de/thor/DATA/www

The fragmentation mode of high-mass molecular clumps and the properties of the central rotating structures surrounding the most luminous objects have yet to be comprehensively characterised. Using the IRAM NOrthern Extended Millimeter Array (NOEMA) and the IRAM 30-m telescope, the CORE survey has obtained high-resolution observations of 20 well-known highly luminous star-forming regions in the 1.37 mm wavelength regime in both line and dust continuum emission. We present the spectral line setup of the CORE survey and a case study for W3(H2O). At ~0.35" (700 AU at 2 kpc) resolution, the W3(H2O) clump fragments into two cores (West and East), separated by ~2300 AU. Velocity shifts of a few km/s are observed in the dense-gas tracer, CH3CN, across both cores, consistent with rotation and perpendicular to the directions of two bipolar outflows, one emanating from each core. The kinematics of the rotating structure about W3(H2O) W shows signs of differential rotation of material, possibly in a disk-like object. The observed rotational signature around W3(H2O) E may be due to a disk-like object, an unresolved binary (or multiple) system, or a combination of both. We fit the emission of CH3CN (12-11) K = 4-6 and derive a gas temperature map with a median temperature of ~165 K across W3(H2O). We create a Toomre Q map to study the stability of the rotating structures against gravitational instability. The rotating structures appear to be Toomre unstable close to their outer boundaries, with a possibility of further fragmentation in the differentially-rotating core W3(H2O) W. Rapid cooling in the Toomre-unstable regions supports the fragmentation scenario. Combining millimeter dust continuum and spectral line data toward the famous high-mass star-forming region W3(H2O), we identify core fragmentation on large scales, and indications for possible disk fragmentation on smaller spatial scales., Comment: 23 pages, 26 figures, accepted for publication in Astronomy and Astrophysics

We apply an analytic theory for environmentally-dependent molecular cloud lifetimes to the Central Molecular Zone of the Milky Way. Within this theory, the cloud lifetime in the Galactic centre is obtained by combining the time-scales for gravitational instability, galactic shear, epicyclic perturbations and cloud-cloud collisions. We find that at galactocentric radii $\sim 45$-$120$ pc, corresponding to the location of the '100-pc stream', cloud evolution is primarily dominated by gravitational collapse, with median cloud lifetimes between 1.4 and 3.9 Myr. At all other galactocentric radii, galactic shear dominates the cloud lifecycle, and we predict that molecular clouds are dispersed on time-scales between 3 and 9 Myr, without a significant degree of star formation. Along the outer edge of the 100-pc stream, between radii of 100 and 120 pc, the time-scales for epicyclic perturbations and gravitational free-fall are similar. This similarity of time-scales lends support to the hypothesis that, depending on the orbital geometry and timing of the orbital phase, cloud collapse and star formation in the 100-pc stream may be triggered by a tidal compression at pericentre. Based on the derived time-scales, this should happen in approximately 20 per cent of all accretion events onto the 100-pc stream., Comment: 6 pages, 2 figures, 1 table; accepted by MNRAS Letters (1 May 2018)

Aims: We aim to understand the fragmentation as well as the disk formation, outflow generation and chemical processes during high-mass star formation on spatial scales of individual cores. Methods: Using the IRAM Northern Extended Millimeter Array (NOEMA) in combination with the 30m telescope, we have observed in the IRAM large program CORE the 1.37mm continuum and spectral line emission at high angular resolution (~0.4'') for a sample of 20 well-known high-mass star-forming regions with distances below 5.5kpc and luminosities larger than 10^4Lsun. Results: We present the overall survey scope, the selected sample, the observational setup and the main goals of CORE. Scientifically, we concentrate on the mm continuum emission on scales on the order of 1000AU. We detect strong mm continuum emission from all regions, mostly due to the emission from cold dust. The fragmentation properties of the sample are diverse. We see extremes where some regions are dominated by a single high-mass core whereas others fragment into as many as 20 cores. A minimum-spanning-tree analysis finds fragmentation at scales on the order of the thermal Jeans length or smaller suggesting that turbulent fragmentation is less important than thermal gravitational fragmentation. The diversity of highly fragmented versus singular regions can be explained by varying initial density structures and/or different initial magnetic field strengths. Conclusions: The smallest observed separations between cores are found around the angular resolution limit which indicates that further fragmentation likely takes place on even smaller spatial scales. The CORE project with its numerous spectral line detections will address a diverse set of important physical and chemical questions in the field of high-mass star formation., Comment: 29 pages, 17 figures, Astronomy & Astrophysics in press, for a higher-resolution version of the paper see http://www.mpia.de/homes/beuther/papers.html

The formation of the most massive O-type stars is poorly understood. We present a case study of a young massive clump from the ATLASGAL survey, G328.2551-0.5321. It exhibits a bolometric luminosity of 1.3$\times$10$^4$ L$_{\odot}$ corresponding to a current protostellar mass of $\sim$11 and 16 M$_{\odot}$. We analyze high angular-resolution observations with ALMA at $\sim$0.17" corresponding a physical scale of $\sim$400 au in dust continuum and molecular lines. The dust continuum emission reveals a single high-mass protostellar envelope and shows evidence for a marginally resolved continuum source. We detect a rotational line of CH$_3$OH within its $v_{\rm t}$=1 torsionally excited state revealing two bright peaks of emission spatially offset from the dust continuum peak, and exhibiting a distinct velocity component $\pm$4.5 km s$^{-1}$ offset compared to the source $v_{\rm lsr}$. Local thermodynamic equilibrium analysis suggests N(CH$_3$OH)=1.2$-$2$\times$10$^{19}$ cm$^{-2}$, and kinetic temperatures of 160$-$200 K at the position of these peaks. Their velocity shifts correspond well to the expected Keplerian velocity around a central object with 15M$_{\odot}$ consistent with the mass estimate based on the source's bolometric luminosity. We propose a picture where the CH$_3$OH emission peaks trace the accretion shocks around the centrifugal barrier, pinpointing the interaction region between the collapsing envelope and an accretion disk. Because the HC$_3$N $v_{\rm 7}$=1e ($J$=38-37) line shows compact emission, and a velocity pattern consistent with models of Keplerian rotation, we suggest that this could be a new tracer for compact accretion disks around high-mass protostars. The estimated physical properties of the accretion disk suggest a specific angular momentum several times larger than typically observed towards low-mass protostars., Comment: Manuscript submitted to A&A, and revised after first referee report (abstract abbreviated)

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., Comment: accepted for publication in Astronomy & Astrophysics; corrections before publication

We have compared the occurrence of 6.7-GHz and 12.2-GHz methanol masers with 22-GHz water masers and 6035-MHz excited-state OH masers in the 100 square degree region of the southern Galactic plane common to the Methanol Multibeam (MMB) and H2O southern Galactic Plane surveys (HOPS). We find the most populous star formation species to be 6.7-GHz methanol, followed by water, then 12.2-GHz and, finally, excited-state OH masers. We present association statistics, flux density (and luminosity where appropriate) and velocity range distributions across the largest, fully surveyed portion of the Galactic plane for four of the most common types of masers found in the vicinity of star formation regions. Comparison of the occurrence of the four maser types with far-infrared dust temperatures shows that sources exhibiting excited-state OH maser emission are warmer than sources showing any of the other three maser types. We further find that sources exhibiting both 6.7-GHz and 12.2-GHz methanol masers are warmer than sources exhibiting just 6.7-GHz methanol maser emission. These findings are consistent with previously made suggestions that both OH and 12.2-GHz methanol masers generally trace a later stage of star formation compared to other common maser types., Comment: 16 pages, 12 figures, accepted MNRAS

The star formation rate in the Central Molecular Zone (CMZ) is an order of magnitude lower than predicted according to star formation relations that have been calibrated in the disc of our own and nearby galaxies. Understanding how and why star formation appears to be different in this region is crucial if we are to understand the environmental dependence of the star formation process. Here, we present the detection of a sample of high-mass cores in the CMZ's "dust ridge" that have been discovered with the Submillimeter Array as part of the CMZoom survey. These cores range in mass from ~ 50 - 2150 Msun within radii of 0.1 - 0.25 pc. All appear to be young (pre-UCHII), meaning that they are prime candidates for representing the initial conditions of high-mass stars and sub-clusters. We report that at least two of these cores ('c1' and 'e1') contain young, high-mass protostars. We compare all of the detected cores with high-mass cores in the Galactic disc and find that they are broadly similar in terms of their masses and sizes, despite being subjected to external pressures that are several orders of magnitude greater - ~ 10^8 K/cm^3, as opposed to ~ 10^5 K/cm^3. The fact that > 80% of these cores do not show any signs of star-forming activity in such a high-pressure environment leads us to conclude that this is further evidence for an increased critical density threshold for star formation in the CMZ due to turbulence., Comment: 19 pages, 10 figures, 3 tables. Accepted for publication in MNRAS

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., Comment: Accepted for publication in A&A

The inner few hundred parsecs of the Milky Way harbours gas densities, pressures, velocity dispersions, an interstellar radiation field and a cosmic ray ionisation rate orders of magnitude higher than the disc; akin to the environment found in star-forming galaxies at high-redshift. Previous studies have shown that this region is forming stars at a rate per unit mass of dense gas which is at least an order of magnitude lower than in the disc, potentially violating theoretical predictions. We show that all observational star formation rate diagnostics - both direct counting of young stellar objects and integrated light measurements - are in agreement within a factor two, hence the low star formation rate is not the result of the systematic uncertainties that affect any one method. As these methods trace the star formation over different timescales, from $0.1 - 5$ Myr, we conclude that the star formation rate has been constant to within a factor of a few within this time period. We investigate the progression of star formation within gravitationally bound clouds on $\sim$ parsec scales and find $1 - 4$ per cent of the cloud masses are converted into stars per free-fall time, consistent with a subset of the considered "volumetric" star formation models. However, discriminating between these models is obstructed by the current uncertainties on the input observables and, most importantly and urgently, by their dependence on ill-constrained free parameters. The lack of empirical constraints on these parameters therefore represents a key challenge in the further verification or falsification of current star formation theories., Comment: 26 pages, 10 figures, 7 tables. Accepted for publication in MNRAS

The H2O Southern Galactic Plane Survey (HOPS) has mapped 100 square degrees of the Galactic plane for water masers and thermal molecular line emission using the 22-m Mopra telescope. We describe the automated spectral-line fitting pipelines used to determine the properties of emission detected in HOPS datacubes, and use these to derive the physical and kinematic properties of gas in the survey. A combination of the angular resolution, sensitivity, velocity resolution and high critical density of lines targeted make the HOPS data cubes ideally suited to finding precursor clouds to the most massive and dense stellar clusters in the Galaxy. We compile a list of the most massive HOPS ammonia regions and investigate whether any may be young massive cluster progenitor gas clouds. HOPS is also ideally suited to trace the flows of dense gas in the Galactic Centre. We find the kinematic structure of gas within the inner 500pc of the Galaxy is consistent with recent predictions for the dynamical evolution of gas flows in the centre of the Milky Way. We confirm a recent finding that the dense gas in the inner 100pc has an oscillatory kinematic structure with characteristic length scale of ~20pc, and also identify similar oscillatory kinematic structure in the gas at radii larger than 100pc. Finally, we make all of the above fits and the remaining HOPS data cubes across the 100 square degrees of the survey available to the community., Comment: 53 pages, 34 figures, 4 tables. Submitted MNRAS. Paper with full resolution images can be downloaded here: http://www.astro.ljmu.ac.uk/~snl/hops_p3.pdf

We present the first data release of the James Clerk Maxwell Telescope (JCMT) Plane Survey (JPS), the JPS Public Release 1 (JPSPR1). JPS is an 850-um continuum survey of six fields in the northern inner Galactic Plane in a longitude range of l=7-63, made with the Sub-millimetre Common-User Bolometer Array 2 (SCUBA-2). This first data release consists of emission maps of the six JPS regions with an average pixel-to-pixel noise of 7.19 mJy beam^-1, when smoothed over the beam, and a compact-source catalogue containing 7,813 sources. The 95 per cent completeness limits of the catalogue are estimated at 0.04 Jy beam^-1 and 0.3 Jy for the peak and integrated flux densities, respectively. The emission contained in the compact-source catalogue is 42 +- 5 per cent of the total and, apart from the large-scale (greater than 8') emission, there is excellent correspondence with features in the 500-um Herschel maps. We find that, with two-dimensional matching, 98 +- 2 per cent of sources within the fields centred at l=20, 30, 40 and 50 are associated with molecular clouds, with 91 +- 3 per cent of the l=30 and 40 sources associated with dense molecular clumps. Matching the JPS catalogue to Herschel 70-um sources, we find that 38 +- 1 per cent of sources show evidence of ongoing star formation. The images and catalogue will be a valuable resource for studies of star formation in the Galaxy and the role of environment and spiral arms in the star formation process., Comment: 22 pages, 18 figures, 4 tables. Full version of Table 3 available from http://www.canfar.phys.uvic.ca/vosui/#/JPSPR1 Accepted for publication in Monthly Notices of the Royal Astronomical Society

Aims: We resolve the small-scale structure around the high-mass hot core region G351.77-0.54 to investigate its disk and fragmentation properties. Methods: Using ALMA at 690GHz with baselines exceeding 1.5km, we study the dense gas, dust and outflow emission at an unprecedented spatial resolution of 0.06" (130AU@2.2kpc). Results: Within the inner few 1000AU, G351.77 fragments into at least four cores (brightness temperatures between 58 and 197K). The central structure around the main submm source #1 with a diameter of ~0.5" does not show additional fragmentation. While the CO(6-5) line wing emission shows an outflow lobe in the north-western direction emanating from source #1, the dense gas tracer CH3CN shows a velocity gradient perpendicular to the outflow that is indicative of rotational motions. Absorption profile measurements against the submm source #2 indicate infall rates on the order of 10^{-4} to 10^{-3}M_sun/yr which can be considered as an upper limit of the mean accretion rates. The position-velocity diagrams are consistent with a central rotating disk-like structure embedded in an infalling envelope, but they may also be influenced by the outflow. Using the CH_3CN(37_k-36_k) k-ladder with excitation temperatures up to 1300K, we derive a gas temperature map of source #1 exhibiting temperatures often in excess of 1000K. Brightness temperatures of the submm continuum never exceed 200K. This discrepancy between gas temperatures and submm dust brightness temperatures (in the optically thick limit) indicates that the dust may trace the disk mid-plane whereas the gas could be tracing a hotter gaseous disk surface layer. In addition, we conduct a pixel-by-pixel Toomre gravitational stability analysis of the central rotating structure. The derived high Q values throughout the structure confirm that this central region appears stable against gravitational instability., Comment: 13 pages, 12 figures, accepted for Astronomy & astrophysics, a high-resolution version of the draft can be found at http://www.mpia-hd.mpg.de/homes/beuther/papers.html

The early evolution of massive cluster progenitors is poorly understood. We investigate the fragmentation properties from 0.3 pc to 0.06 pc scales of a homogenous sample of infrared-quiet massive clumps within 4.5 kpc selected from the ATLASGAL survey. Using the ALMA 7m array we detect compact dust continuum emission towards all targets, and find that fragmentation, at these scales, is limited. The mass distribution of the fragments uncovers a large fraction of cores above 40 $M_\odot$, corresponding to massive dense cores (MDCs) with masses up to ~400 $M_\odot$. 77 % of the clumps contain at most 3 MDCs per clump, and we also reveal single clumps/MDCs. The most massive cores are formed within the more massive clumps, and a high concentration of mass on small scales reveals a high core formation efficiency. The mass of MDCs highly exceeds the local thermal Jeans-mass, and observational evidence is lacking for a sufficiently high level of turbulence or strong enough magnetic fields to keep the most massive MDCs in equilibrium. If already collapsing, the observed fragmentation properties with a high core formation efficiency are consistent with the collapse setting in at parsec scales., Comment: accepted by A&A Letters

In this paper we describe an unmanned aerial system equipped with a thermal-infrared camera and software pipeline that we have developed to monitor animal populations for conservation purposes. Taking a multi-disciplinary approach to tackle this problem, we use freely available astronomical source detection software and the associated expertise of astronomers, to efficiently and reliably detect humans and animals in aerial thermal-infrared footage. Combining this astronomical detection software with existing machine learning algorithms into a single, automated, end-to-end pipeline, we test the software using aerial video footage taken in a controlled, field-like environment. We demonstrate that the pipeline works reliably and describe how it can be used to estimate the completeness of different observational datasets to objects of a given type as a function of height, observing conditions etc. -- a crucial step in converting video footage to scientifically useful information such as the spatial distribution and density of different animal species. Finally, having demonstrated the potential utility of the system, we describe the steps we are taking to adapt the system for work in the field, in particular systematic monitoring of endangered species at National Parks around the world., Comment: 16 pages, 6 figures, accepted to International Journal of Remote Sensing

Yellow hypergiants are rare and represent a fast evolutionary stage of massive evolved stars. That evolutionary phase is characterised by a very intense mass loss, the understanding of which is still very limited. Here we report ALMA Compact Array observations of a 50$"$-mosaic toward the Fried Egg nebula, around one of the few Galactic yellow hypergiants IRAS 17163-3907. The emission from the $^{12}$CO J=2-1 line, H30$\alpha$ recombination line, and continuum is imaged at a resolution of $\sim$8$"$, revealing the morphology of the molecular environment around the star. The continuum emission is unresolved and peaks at the position of the star. The radio recombination line H30$\alpha$ shows unresolved emission at the star, with an approximately gaussian spectrum centered on a velocity of 21$\pm$3~km/s with a width of 57$\pm$6~km/s. In contrast, the CO 2-1 emission is complex and decomposes into several components beyond the contamination from interstellar gas in the line of sight. The CO spectrum toward the star is a broad plateau, centered at the systemic velocity of +18 km/s and with an expansion velocity of 100$\pm$10 km/s. Assuming isotropic and constant mass-loss, we estimate a mass-loss rate of 8$\pm$1.5 $\times10^{-5}$~M$_\odot$ yr$^{-1}$. At a radius of 25$"$ from the star, we detect CO emission associated with the dust ring previously imaged by {\it Herschel}. The kinematics of this ring, however, is not consistent with an expanding shell, but show a velocity gradient of $v_{sys} \pm$20 km/s. In addition, we find a puzzling bright feature radially connecting the star to the CO ring, at a velocity of +40 km/s relative to the star. This spur feature may trace a unidirectional ejection event from the star. Our ACA observations reveal the complex morphology around IRAS 17163 and illustrate the breakthroughs that ALMA will bring to the field of massive stellar evolution., Comment: Accepted for publication in A&A. 8 pages, 7 figures, plus appendix

Star formation in the Galactic disc is primarily controlled by gravity, turbulence, and magnetic fields. It is not clear that this also applies to star formation near the Galactic Centre. Here we determine the turbulence and star formation in the CMZ cloud G0.253+0.016. Using maps of 3mm dust emission and HNCO intensity-weighted velocity obtained with ALMA, we measure the volume-density variance $\sigma_{\rho/\rho_0} = 1.3 \pm 0.5$ and turbulent Mach number $\mathcal{M} = 11 \pm 3$. Combining these with turbulence simulations to constrain the plasma $\beta = 0.34 \pm 0.35$, we reconstruct the turbulence driving parameter $b = 0.22 \pm 0.12$ in G0.253+0.016. This low value of $b$ indicates solenoidal (divergence-free) driving of the turbulence in G0.253+0.016. By contrast, typical clouds in the Milky Way disc and spiral arms have a significant compressive (curl-free) driving component ($b > 0.4$). We speculate that shear causes the solenoidal driving in G0.253+0.016 and show that this may reduce the star formation rate by a factor of 7 compared to nearby clouds., Comment: IAU Symposium 322 "The Multi-Messenger Astrophysics of the Galactic Centre", eds. R. M. Crocker, S. N. Longmore, G. V. Bicknell

Star formation is primarily controlled by the interplay between gravity, turbulence, and magnetic fields. However, the turbulence and magnetic fields in molecular clouds near the Galactic Center may differ substantially from spiral-arm clouds. Here we determine the physical parameters of the central molecular zone (CMZ) cloud G0.253+0.016, its turbulence, magnetic field and filamentary structure. Using column-density maps based on dust-continuum emission observations with ALMA+Herschel, we identify filaments and show that at least one dense core is located along them. We measure the filament width W_fil=0.17$\pm$0.08pc and the sonic scale {\lambda}_sonic=0.15$\pm$0.11pc of the turbulence, and find W_fil~{\lambda}_sonic. A strong velocity gradient is seen in the HNCO intensity-weighted velocity maps obtained with ALMA+Mopra, which is likely caused by large-scale shearing of G0.253+0.016, producing a wide double-peaked velocity PDF. After subtracting the gradient to isolate the turbulent motions, we find a nearly Gaussian velocity PDF typical for turbulence. We measure the total and turbulent velocity dispersion, 8.8$\pm$0.2km/s and 3.9$\pm$0.1km/s, respectively. Using magnetohydrodynamical simulations, we find that G0.253+0.016's turbulent magnetic field B_turb=130$\pm$50$\mu$G is only ~1/10 of the ordered field component. Combining these measurements, we reconstruct the dominant turbulence driving mode in G0.253+0.016 and find a driving parameter b=0.22$\pm$0.12, indicating solenoidal (divergence-free) driving. We compare this to spiral-arm clouds, which typically have a significant compressive (curl-free) driving component (b>0.4). Motivated by previous reports of strong shearing motions in the CMZ, we speculate that shear causes the solenoidal driving in G0.253+0.016 and show that this reduces the star formation rate (SFR) by a factor of 6.9 compared to typical nearby clouds., Comment: 18 pages, 6 figures, 1 table containing new Brick parameters; accepted by ApJ

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., Comment: Accepted for Astronomy & Astrophysics, 23 pages, 17 figures, a high-resolution version can be found at http://www.mpia.de/thor