In galactic disks, galactic rotation sets the bulk motion of gas, and its energy and momentum can be transferred toward small scales. Additionally, in the interstellar medium, random and noncircular motions arise from stellar feedback, cloud-cloud interactions, and instabilities, among other processes. Our aim is to comprehend to which extent small-scale gas dynamics is decoupled from galactic rotation. We study the relative contributions of galactic rotation and local noncircular motions to the circulation of gas, $\Gamma$, a macroscopic measure of local rotation, defined as the line integral of the velocity field around a closed path. We measure the circulation distribution as a function of spatial scale in a set of simulated disk galaxies and we model the velocity field as the sum of galactic rotation and a Gaussian random field. The random field is parameterized by a broken power law in Fourier space, with a break at the scale $\lambda_c$. We define the spatial scale $\lambda_{\rm eq}$ at which galactic rotation and non-circular motions contribute equally to $\Gamma$. For our simulated galaxies, the gas dynamics at the scale of molecular clouds is usually dominated by noncircular motions, but in the center of galactic disks galactic rotation is still relevant. Our model shows that the transfer of rotation from large scales breaks at the scale $\lambda_c$ and this transition is necessary to reproduce the circulation distribution. We find that $\lambda_{\rm eq}$, and therefore the structure of the gas velocity field, is set by the local conditions of gravitational stability and stellar feedback.

We present an overview of the CMZoom survey and its first data release. CMZoom is the first blind, high-resolution survey of the Central Molecular Zone (CMZ; the inner 500 pc of the Milky Way) at wavelengths sensitive to the pre-cursors of high-mass stars. CMZoom is a 500-hour Large Program on the Submillimeter Array (SMA) that mapped at 1.3 mm all of the gas and dust in the CMZ above a molecular hydrogen column density of 10^23 cm^-2 at a resolution of ~3" (0.1 pc). In this paper, we focus on the 1.3 mm dust continuum and its data release, but also describe CMZoom spectral line data which will be released in a forthcoming publication. While CMZoom detected many regions with rich and complex substructure, its key result is an overall deficit in compact substructures on 0.1 - 2 pc scales (the compact dense gas fraction: CDGF). In comparison with clouds in the Galactic disk, the CDGF in the CMZ is substantially lower, despite having much higher average column densities. CMZ clouds with high CDGFs are well-known sites of active star formation. The inability of most gas in the CMZ to form compact substructures is likely responsible for the dearth of star formation in the CMZ, surprising considering its high density. The factors responsible for the low CDGF are not yet understood but are plausibly due to the extreme environment of the CMZ, having far-reaching ramifications for our understanding of the star formation process across the cosmos.

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.

We report here Atacama Large Millimeter/submillimeter Array (ALMA) N$_2$H$^+$ (1-0) images of the Orion Molecular Cloud 2 and 3 (OMC-2/3) with high angular resolution (3'' or 1200 au) and high spatial dynamic range. Combining dataset from the ALMA main array, ALMA Compact Array (ACA), the Nobeyama 45m Telescope, and the JVLA (providing temperature measurement on matching scales), we find that most of the dense gas in OMC-2/3 is subsonic ($\rm \sigma_{NT}/c_{s}$ = 0.62) with a mean line width ($\Delta\upsilon$) of 0.39 km s$^{-1}$ FWHM. This is markedly different from the majority of previous observations of massive star-forming regions. In contrast, line widths from the Nobeyama Telescope are transonic at 0.69 km s$^{-1}$ ($\rm \sigma_{NT}/c_{s}$ = 1.08). We demonstrated that the larger line widths obtained by the single-dish telescope arose from unresolved sub-structures within their respective beams. The dispersions from larger scales $\sigma_{ls}$ (as traced by the Nobeyama Telescope) can be decomposed into three components $\rm \sigma_{ls}^2 = \sigma_{ss}^2+ \sigma_{bm}^2+ \sigma_{rd}^2$, where small-scale $\sigma_{ss}$ is the line dispersion of each ALMA beam, bulk motion $\sigma_{bm}$ is dispersion between peak velocity of each ALMA beam, and $\sigma_{rd}$ is the residual dispersion. Such decomposition, though purely empirical, appears to be robust throughout our data cubes. Apparent supersonic line widths, commonly found in massive molecular clouds, are thus likely due to the effect of poor spatial resolution. The observed non-thermal line dispersion (sometimes referred to as 'turbulence') transits from supersonic to subsonic at $\sim 0.05$ pc scales in OMC-2/3 region. Such transition could be commonly found with sufficient spatial (not just angular) resolution, even in regions with massive young clusters, such as Orion molecular clouds studied here.

Knowledge about the distribution of CO emission in the Milky Way is essential to understand the impact of Galactic environment on the formation and evolution of structures in the interstellar medium. However, currently our insight about the fraction of CO in spiral arm and interarm regions is still limited by large uncertainties in assumed rotation curve models or distance determination techniques. In this work we use a Bayesian approach to obtain the current best assessment of the distribution of 13CO from the Galactic Ring Survey. We performed two different distance estimates that either included or excluded a model for Galactic features. We also include a prior for the solution of the kinematic distance ambiguity that was determined from a compilation of literature distances and an assumed size-linewidth relationship. We find that the fraction of 13CO emission associated with spiral arm features varies from 76% to 84% between the two distance runs. The vertical distribution of the gas is concentrated around the Galactic midplane showing FWHM values of ~75 pc. We do not find any significant difference between gas emission properties associated with spiral arm and interarm features. In particular the distribution of velocity dispersion values of gas emission in spurs and spiral arms is very similar. We detect a trend of higher velocity dispersion values with increasing heliocentric distance, which we attribute to beam averaging effects caused by differences in spatial resolution. We argue that the true distribution of the gas emission is likely more similar to a combination of our two distance results, and highlight the importance of using complementary distance estimations to safeguard against the pitfalls of any single approach. We conclude that the methodology presented in this work is a promising way to determine distances to gas emission features in Galactic plane surveys.

We present kinematic orientations and high resolution (150 pc) rotation curves for 67 main sequence star-forming galaxies surveyed in CO (2-1) emission by PHANGS-ALMA. Our measurements are based on the application of a new fitting method tailored to CO velocity fields. Our approach identifies an optimal global orientation as a way to reduce the impact of non-axisymmetric (bar and spiral) features and the uneven spatial sampling characteristic of CO emission in the inner regions of nearby galaxies. The method performs especially well when applied to the large number of independent lines-of-sight contained in the PHANGS CO velocity fields mapped at 1'' resolution. The high resolution rotation curves fitted to these data are sensitive probes of mass distribution in the inner regions of these galaxies. We use the inner slope as well as the amplitude of our fitted rotation curves to demonstrate that CO is a reliable global dynamical mass tracer. From the consistency between photometric orientations from the literature and kinematic orientations determined with our method, we infer that the shapes of stellar disks in the mass range of log($\rm M_{\star}(M_{\odot})$)=9.0-10.9 probed by our sample are very close to circular and have uniform thickness.

Context: Atomic and molecular cloud formation is a dynamical process. However, kinematic signatures of these processes are still observationally poorly constrained. Methods: Targeting the cloud-scale environment of the prototypical infrared dark cloud G28.3, we employ spectral line imaging observations of the two atomic lines HI and [CI] as well as molecular lines observations in 13CO in the 1--0 and 3--2 transitions. The analysis comprises investigations of the kinematic properties of the different tracers, estimates of the mass flow rates, velocity structure functions, a Histogram of Oriented Gradients (HOG) study as well as comparisons to simulations. Results: The central IRDC is embedded in a more diffuse envelope of cold neutral medium (CNM) traced by HI self-absorption (HISA) and molecular gas. The spectral line data as well as the HOG and structure function analysis indicate a possible kinematic decoupling of the HI from the other gas compounds. Spectral analysis and position-velocity diagrams reveal two velocity components that converge at the position of the IRDC. Estimated mass flow rates appear rather constant from the cloud edge toward the center. The velocity structure function analysis is consistent with gas flows being dominated by the formation of hierarchical structures. Conclusions: The observations and analysis are consistent with a picture where the IRDC G28 is formed at the center of two converging gas flows. While the approximately constant mass flow rates are consistent with a self-similar, gravitationally driven collapse of the cloud, external compression by, e.g., spiral arm shocks or supernovae explosions cannot be excluded yet. Future investigations should aim at differentiating the origin of such converging gas flows.

We compare the observed turbulent pressure in molecular gas, $P_\mathrm{turb}$, to the required pressure for the interstellar gas to stay in equilibrium in the gravitational potential of a galaxy, $P_\mathrm{DE}$. To do this, we combine arcsecond resolution CO data from PHANGS-ALMA with multi-wavelength data that traces the atomic gas, stellar structure, and star formation rate (SFR) for 28 nearby star-forming galaxies. We find that $P_\mathrm{turb}$ correlates with, but almost always exceeds the estimated $P_\mathrm{DE}$ on kiloparsec scales. This indicates that the molecular gas is over-pressurized relative to the large-scale environment. We show that this over-pressurization can be explained by the clumpy nature of molecular gas; a revised estimate of $P_\mathrm{DE}$ on cloud scales, which accounts for molecular gas self-gravity, external gravity, and ambient pressure, agrees well with the observed $P_\mathrm{turb}$ in galaxy disks. We also find that molecular gas with cloud-scale ${P_\mathrm{turb}}\approx{P_\mathrm{DE}}\gtrsim{10^5\,k_\mathrm{B}\,\mathrm{K\,cm^{-3}}}$ in our sample is more likely to be self-gravitating, whereas gas at lower pressure appears more influenced by ambient pressure and/or external gravity. Furthermore, we show that the ratio between $P_\mathrm{turb}$ and the observed SFR surface density, $\Sigma_\mathrm{SFR}$, is compatible with stellar feedback-driven momentum injection in most cases, while a subset of the regions may show evidence of turbulence driven by additional sources. The correlation between $\Sigma_\mathrm{SFR}$ and kpc-scale $P_\mathrm{DE}$ in galaxy disks is consistent with the expectation from self-regulated star formation models. Finally, we confirm the empirical correlation between molecular-to-atomic gas ratio and kpc-scale $P_\mathrm{DE}$ reported in previous works.

The analysis of large molecular line surveys of the Galactic plane is essential for our understanding of the gas kinematics on Galactic scales, in particular its link with the formation and evolution of dense structures in the interstellar medium. An approximation of the emission peaks with Gaussian functions allows for an efficient and straightforward extraction of useful physical information contained in the shape and Doppler-shifted frequency of the emission lines contained in these enormous data sets. In this work we present an overview and first results of a Gaussian decomposition of the entire Galactic Ring Survey (GRS) 13CO (1-0) data that consists of about 2.3 million spectra. We performed the decomposition with the fully automated GaussPy+ algorithm and fitted about 4.6 million Gaussian components to the GRS spectra. We discuss the statistics of the fit components and relations between the fitted intensities, velocity centroids, and velocity dispersions. We find that the magnitude of the velocity dispersion values increase toward the inner Galaxy and around the Galactic midplane, which we speculate is partly due to the influence of the Galactic bar and regions with higher non-thermal motions located in the midplane, respectively. We also use our decomposition results to infer global properties of the gas emission and find that the number of fit components used per spectrum is indicative for the amount of structure along the line of sight. We find that the emission lines from regions located on the far side of the Galaxy show increased velocity dispersion values, likely due to beam averaging effects. We demonstrate how this trend has the potential to aid in characterising Galactic structure by disentangling emission that is belonging to the nearby Aquila Rift molecular cloud from emission that is more likely associated with the Perseus and Outer spiral arms.

It remains a major challenge to derive a theory of cloud-scale ($\lesssim100$ pc) star formation and feedback, describing how galaxies convert gas into stars as a function of the galactic environment. Progress has been hampered by a lack of robust empirical constraints on the giant molecular cloud (GMC) lifecycle. We address this problem by systematically applying a new statistical method for measuring the evolutionary timeline of the GMC lifecycle, star formation, and feedback to a sample of nine nearby disc galaxies, observed as part of the PHANGS-ALMA survey. We measure the spatially-resolved ($\sim100$ pc) CO-to-H$\alpha$ flux ratio and find a universal de-correlation between molecular gas and young stars on GMC scales, allowing us to quantify the underlying evolutionary timeline. GMC lifetimes are short, typically 10-30 Myr, and exhibit environmental variation, between and within galaxies. At kpc-scale molecular gas surface densities $\Sigma_{\rm H_2}\geqslant8$M$_{\odot}$pc$^{-2}$, the GMC lifetime correlates with time-scales for galactic dynamical processes, whereas at $\Sigma_{\rm H_2}\leqslant8$M$_{\odot}$pc$^{-2}$ GMCs decouple from galactic dynamics and live for an internal dynamical time-scale. After a long inert phase without massive star formation traced by H$\alpha$ (75-90% of the cloud lifetime), GMCs disperse within just 1-5 Myr once massive stars emerge. The dispersal is most likely due to early stellar feedback, causing GMCs to achieve integrated star formation efficiencies of 4-10% These results show that galactic star formation is governed by cloud-scale, environmentally-dependent, dynamical processes driving rapid evolutionary cycling. GMCs and HII regions are the fundamental units undergoing these lifecycles, with mean separations of 100-300 pc in star-forming discs. Future work should characterise the multi-scale physics and mass flows driving these lifecycles.

The processes regulating star formation in galaxies are thought to act across a hierarchy of spatial scales. To connect extragalactic star formation relations from global and kpc-scale measurements to recent cloud-scale resolution studies, we have developed a simple, robust method that quantifies the scale dependence of the relative spatial distributions of molecular gas and recent star formation. In this paper, we apply this method to eight galaxies with roughly 1 arcsec resolution molecular gas imaging from the PHANGS-ALMA and PAWS surveys that have matched resolution, high quality narrowband Halpha imaging. At a common scale of 140pc, our massive (log(Mstar/Msun)=9.3-10.7), normally star-forming (SFR/Msun/yr=0.3-5.9) galaxies exhibit a significant reservoir of quiescent molecular gas not associated with star formation as traced by Halpha emission. Galactic structures act as backbones for both molecular and HII region distributions. As we degrade the spatial resolution, the quiescent molecular gas disappears, with the most rapid changes occurring for resolutions up to about 0.5kpc. As the resolution becomes poorer, the morphological features become indistinct for spatial scales larger than about 1kpc. The method is a promising tool to search for relationships between the quiescent or star-forming molecular reservoir and galaxy properties, but requires a larger sample size to identify robust correlations between the star-forming molecular gas fraction and global galaxy parameters.

We present new observations of C-band continuum emission and masers to assess high-mass ($>$8 $M_\odot$) star formation at early evolutionary phases in the inner 200 pc of the Central Molecular Zone (CMZ) of the Galaxy. The continuum observation is complete to free-free emission from stars above 10-11 $M_\odot$ in 91% of the covered area. We identify 104 compact sources in the continuum emission, among which five are confirmed ultracompact H II regions, 12 are candidates of ultra-compact H II regions, and the remaining 87 sources are mostly massive stars in clusters, field stars, evolved stars, pulsars, extragalactic sources, or of unknown nature that is to be investigated. We detect class II CH$_3$OH masers at 23 positions, among which six are new detections. We confirm six known H$_2$CO masers in two high-mass star forming regions, and detect two new H$_2$CO masers toward the Sgr C cloud, making it the ninth region in the Galaxy that contains masers of this type. In spite of these detections, we find that current high-mass star formation in the inner CMZ is only taking place in seven isolated clouds. The results suggest that star formation at early evolutionary phases in the CMZ is about 10 times less efficient than expected by the dense gas star formation relation, which is in line with previous studies that focus on more evolved phases of star formation. This means that if there will be any impending, next burst of star formation in the CMZ, it has not yet begun.

Magneto-hydrodynamic (MHD) shocks are violent events that inject large amounts of energy in the interstellar medium (ISM) dramatically modifying its physical properties and chemical composition. Indirect evidence for the presence of such shocks has been reported from the especial chemistry detected toward a variety of astrophysical shocked environments. However, the internal physical structure of these shocks remains unresolved since their expected spatial scales are too small to be measured with current instrumentation. Here we report the first detection of a fully spatially resolved, MHD shock toward the Infrared Dark Cloud (IRDC) G034.77-00.55. The shock, probed by Silicon Monoxide (SiO) and observed with the Atacama Large Millimetre/sub-millimetre Array (ALMA), is associated with the collision between the dense molecular gas of the cloud and a molecular gas flow pushed toward the IRDC by the nearby supernova remnant (SNR) W44. The interaction is occurring on sub-parsec spatial scales thanks to the enhanced magnetic field of the SNR, making the dissipation region of the MHD shock large enough to be resolved with ALMA. Our observations suggest that molecular flow-flow collisions can be triggered by stellar feedback, inducing shocked molecular gas densities compatible with those required for massive star formation.

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

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.

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.

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.

To probe the initial conditions for high-mass star and cluster formation, we investigate the properties of dense filaments within the infrared dark cloud G035.39-00.33 (IRDC G035.39) in a combined Very Large Array (VLA) and the Green Bank Telescope (GBT) mosaic tracing the NH3 (1,1) and (2,2) emission down to 0.08 pc scales. Using agglomerative hierarchical clustering on multiple line-of-sight velocity component fitting results, we identify seven extended velocity-coherent components in our data, likely representing spatially coherent physical structures, some exhibiting complex gas motions. The velocity gradient magnitude distribution peaks at its mode of 0.35 km/s/pc and has a long tail extending into higher values of 1.5 - 2 km/s/pc, and is generally consistent with those found toward the same cloud in other molecular tracers and with the values found towards nearby low-mass dense cloud cores at the same scales. Contrary to observational and theoretical expectations, we find the non-thermal ammonia line widths to be systematically narrower (by about 20%) than those of N2H+ (1-0) line transition observed with similar resolution. If the observed ordered velocity gradients represent the core envelope solid-body rotation, we estimate the specific angular momentum to be about 2 x 10^21 cm^2/s, similar to the low-mass star-forming cores. Together with the previous finding of subsonic motions in G035.39, our results demonstrate high levels of similarity between kinematics of a high-mass star-forming IRDC and the low-mass star formation regime.

Stellar feedback plays a fundamental role in shaping the evolution of galaxies. Here we explore the use of ionised gas kinematics in young, bipolar H II regions as a probe of early feedback in these star-forming environments. We have undertaken a multiwavelength study of a young, bipolar H II region in the Galactic disc, G$316.81-0.06$, which lies at the centre of a massive ($\sim10^3$ M$_{\odot}$) infrared-dark cloud filament. It is still accreting molecular gas as well as driving a $\sim 0.2$ pc ionised gas outflow perpendicular to the filament. Intriguingly, we observe a large velocity gradient ($47.81 \pm 3.21$ km s$^{-1}$ pc$^{-1}$) across the ionised gas in a direction perpendicular to the outflow. This kinematic signature of the ionised gas shows a reasonable correspondence with the simulations of young H II regions. Based on a qualitative comparison between our observations and these simulations, we put forward a possible explanation for the velocity gradients observed in G$316.81-0.06$. If the velocity gradient perpendicular to the outflow is caused by rotation of the ionised gas, then we infer that this rotation is a direct result of the initial net angular momentum in the natal molecular cloud. If this explanation is correct, this kinematic signature should be common in other young (bipolar) H II regions. We suggest that further quantitative analysis of the ionised gas kinematics of young H II regions, combined with additional simulations, should improve our understanding of feedback at these early stages.

High-mass star forming regions are typically thought to be dominated by supersonic motions. We present combined Very Large Array and Green Bank Telescope (VLA+GBT) observations of NH$_3$ (1,1) and (2,2) in the infrared dark cloud (IRDC) G035.39-00.33, tracing cold and dense gas down to scales of 0.07 pc. We find that, in contrast to previous similar studies of IRDCs, more than a third of the fitted ammonia spectra show subsonic non-thermal motions (mean line width of 0.71 $\mathrm{km~s^{-1}}$), and the sonic Mach number distribution peaks around $\mathcal{M} = 1$. As possible observational and instrumental biases would only broaden the line profiles, our results provide strong upper limits to the actual value of $\mathcal{M}$, further strengthening our findings of narrow line widths. This finding calls for a reevaluation of the role of turbulent dissipation and subsonic regions in massive-star and cluster formation. Based on our findings in G035.39, we further speculate that the coarser spectral resolution used in the previous VLA NH$_3$ studies may have inhibited the detection of subsonic turbulence in IRDCs. The reduced turbulent support suggests that dynamically important magnetic fields of the 1 mG order would be required to support against possible gravitational collapse. Our results offer valuable input into the theories and simulations that aim to recreate the initial conditions of high-mass star and cluster formation.

Infrared dark clouds (IRDCs) are thought to be potential hosts of the elusive early phases of high-mass star formation. Here we conduct an in-depth kinematic analysis of one such IRDC, G034.43+00.24 (Cloud F), using high sensitivity and high spectral resolution IRAM-30m N$_2$H$^+$ ($1-0$) and C$^{18}$O ($1-0$) observations. To disentangle the complex velocity structure within this cloud we use Gaussian decomposition and hierarchical clustering algorithms. We find that four distinct coherent velocity components are present within Cloud F. The properties of these components are compared to those found in a similar IRDC, G035.39-00.33 (Cloud H). We find that the components in both clouds have: high densities (inferred by their identification in N$_2$H$^+$), trans-to-supersonic non-thermal velocity dispersions with Mach numbers of $\sim$ $1.5-4$, a separation in velocity of $\sim$3 km s$^{-1}$, and a mean red-shift of $\sim$ 0.3 km s$^{-1}$ between the N$_2$H$^+$ (dense gas) and C$^{18}$O emission (envelope gas). The latter of these could suggest that these clouds share a common formation scenario. We investigate the kinematics of the larger-scale Cloud F structures, using lower-density-tracing $^{13}$CO ($1-0$) observations. A good correspondence is found between the components identified in the IRAM-30m observations and the most prominent component in the$^{13}$CO data. We find that the IRDC Cloud F is only a small part of a much larger structure, which appears to be an inter-arm filament of the Milky Way.

We report ALMA observations with resolution $\approx0.5$" at 3 mm of the extended Sgr B2 cloud in the Central Molecular Zone (CMZ). We detect 271 compact sources, most of which are smaller than 5000 AU. By ruling out alternative possibilities, we conclude that these sources consist of a mix of hypercompact HII regions and young stellar objects (YSOs). Most of the newly-detected sources are YSOs with gas envelopes which, based on their luminosities, must contain objects with stellar masses $M_*\gtrsim8$ M$_\odot$. Their spatial distribution spread over a $\sim12\times3$ pc region demonstrates that Sgr B2 is experiencing an extended star formation event, not just an isolated `starburst' within the protocluster regions. Using this new sample, we examine star formation thresholds and surface density relations in Sgr B2. While all of the YSOs reside in regions of high column density ($N(H_2)\gtrsim2\times10^{23}$ cm$^{-2}$), not all regions of high column density contain YSOs. The observed column density threshold for star formation is substantially higher than that in solar vicinity clouds, implying either that high-mass star formation requires a higher column density or that any star formation threshold in the CMZ must be higher than in nearby clouds. The relation between the surface density of gas and stars is incompatible with extrapolations from local clouds, and instead stellar densities in Sgr B2 follow a linear $\Sigma_*-\Sigma_{gas}$ relation, shallower than that observed in local clouds. Together, these points suggest that a higher volume density threshold is required to explain star formation in CMZ clouds.

Infrared-Dark Clouds (IRDCs) are cold, dense regions of high (optical and infrared) extinction, believed to be the birthplace of high-mass stars and stellar clusters. The physical mechanisms leading to the formation of these IRDCs are not completely understood and it is thus important to study their molecular gas kinematics and chemical content to search for any signature of the IRDCs formation process. Using the 30m-diameter antenna at the Instituto de Radioastronom\'ia Milim\'etrica, we have obtained emission maps of dense gas tracers (H$^{13}$CO$^{+}$ and HN$^{13}$C) and typical shock tracers (SiO and CH$_3$OH) toward three IRDCs, G028.37+00.07, G034.43+00.24 and G034.77-00.55 (clouds C, F and G, respectively). We have studied the molecular gas kinematics in these clouds and, consistent with previous works toward other IRDCs, the clouds show complex gas kinematics with several velocity-coherent sub-structures separated in velocity space by a few km s$^{-1}$. Correlated with these complex kinematic structures, widespread (parsec-scale) emission of SiO and CH$_3$OH is present in all the three clouds. For clouds C and F, known to be actively forming stars, widespread SiO and CH$_3$OH is likely associated with on-going star formation activity. However, for cloud G, which lacks either 8 $\mu$m or 24 $\mu$m sources and 4.5 $\mu$m H$_2$ shock-excited emission, the detected widespread SiO and CH$_3$OH emission may have originated in a large-scale shock interaction, although a scenario involving a population of low-mass stars driving molecular outflows cannot be fully ruled out.

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.

The Survey of Water and Ammonia in the Galactic Center (SWAG) covers the Central Molecular Zone (CMZ) of the Milky Way at frequencies between 21.2 and 25.4 GHz obtained at the Australia Telescope Compact Array at $\sim 0.9$ pc spatial and $\sim 2.0$ km s$^{-1}$ spectral resolution. In this paper, we present data on the inner $\sim 250$ pc ($1.4^\circ$) between Sgr C and Sgr B2. We focus on the hyperfine structure of the metastable ammonia inversion lines (J,K) = (1,1) - (6,6) to derive column density, kinematics, opacity and kinetic gas temperature. In the CMZ molecular clouds, we find typical line widths of $8-16$ km s$^{-1}$ and extended regions of optically thick ($\tau > 1$) emission. Two components in kinetic temperature are detected at $25-50$ K and $60-100$ K, both being significantly hotter than dust temperatures throughout the CMZ. We discuss the physical state of the CMZ gas as traced by ammonia in the context of the orbital model by Kruijssen et al. (2015) that interprets the observed distribution as a stream of molecular clouds following an open eccentric orbit. This allows us to statistically investigate the time dependencies of gas temperature, column density and line width. We find heating rates between $\sim 50$ and $\sim 100$ K Myr$^{-1}$ along the stream orbit. No strong signs of time dependence are found for column density or line width. These quantities are likely dominated by cloud-to-cloud variations. Our results qualitatively match the predictions of the current model of tidal triggering of cloud collapse, orbital kinematics and the observation of an evolutionary sequence of increasing star formation activity with orbital phase.

Aims. Infrared dark clouds represent the earliest stages of high-mass star formation. Detailed observations of their physical conditions on all physical scales are required to improve our understanding of their role in fueling star formation. Methods. We investigate the large-scale structure of the IRDC G035.39-00.33, probing the dense gas with the classical ammonia thermometer. This allows us to put reliable constraints on the temperature of the extended, pc-scale dense gas reservoir and to probe the magnitude of its non-thermal motions. Available far-infrared observations can be used in tandem with the observed ammonia emission to estimate the total gas mass contained in G035.39-00.33. Results. We identify a main velocity component as a prominent filament, manifested as an ammonia emission intensity ridge spanning more than 6 pc, consistent with the previous studies on the Northern part of the cloud. A number of additional line-of-sight components are found, and a large scale, linear velocity gradient of ~0.2 km s$^{-1}$ pc$^{-1}$ is found along the ridge of the IRDC. In contrast to the dust temperature map, an ammonia-derived kinetic temperature map, presented for the entirety of the cloud, reveals local temperature enhancements towards the massive protostellar cores. We show that without properly accounting for the line of sight contamination, the dust temperature is 2-3 K larger than the gas temperature measured with NH$_3$. Conclusions. While both the large scale kinematics and temperature structure are consistent with that of starless dark filaments, the kinetic gas temperature profile on smaller scales is suggestive of tracing the heating mechanism coincident with the locations of massive protostellar cores.

Nitrogen is one of the most abundant elements in the Universe and its 14N/15N isotopic ratio has the potential to provide information about the initial environment in which our Sun formed. Recent findings suggest that the Solar System may have formed in a massive cluster since the presence of short-lived radioisotopes in meteorites can only be explained by the influence of a supernova. The aim of this project is to determine the 14N/15N ratio towards a sample of cold, massive dense cores at the initial stages in their evolution. We have observed the J=1-0 transitions of HCN, H13CN, HC15N, HN13C and H15NC toward a sample of 22 cores in 4 Infrared-Dark Clouds (IRDCs). IRDCs are believed to be the precursors of high-mass stars and star clusters. Assuming LTE and a temperature of 15K, the column densities of HCN, H13CN, HC15N, HN13C and H15NC are calculated and their 14N/15N ratio is determined for each core. The 14N/15N ratio measured in our sample of IRDC cores range between ~70 and >763 in HCN and between ~161 and ~541 in HNC. They are consistent with the terrestrial atmosphere (TA) and protosolar nebula (PSN) values, and with the ratios measured in low-mass pre-stellar cores. However, the 14N/15N ratios measured in cores C1, C3, F1, F2 and G2 do not agree with the results from similar studies toward the same massive cores using nitrogen bearing molecules with nitrile functional group (-CN) and nitrogen hydrides (-NH) although the ratio spread covers a similar range. Amongst the 4 IRDCs we measured relatively low 14N/15N ratios towards IRDC G which are comparable to those measured in small cosmomaterials and protoplanetary disks. The low average gas density of this cloud suggests that the gas density, rather than the gas temperature, may be the dominant parameter influencing the initial nitrogen isotopic composition in young PSN.

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.

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.

The inner few hundred parsecs of the Milky Way, the Central Molecular Zone (CMZ), is our closest laboratory for understanding star formation in the extreme environments (hot, dense, turbulent gas) that once dominated the universe. We present an update on the first large-area survey to expose the sites of star formation across the CMZ at high-resolution in submillimeter wavelengths: the CMZoom survey with the Submillimeter Array (SMA). We identify the locations of dense cores and search for signatures of embedded star formation. CMZoom is a three-year survey in its final year and is mapping out the highest column density regions of the CMZ in dust continuum and a variety of spectral lines around 1.3 mm. CMZoom combines SMA compact and subcompact configurations with single-dish data from BGPS and the APEX telescope, achieving an angular resolution of about 4" (0.2 pc) and good image fidelity up to large spatial scales.

Star formation within the Central Molecular Zone (CMZ) may be intimately linked to the orbital dynamics of the gas. Recent models suggest that star formation within the dust ridge molecular clouds (from G0.253+0.016 to Sgr B2) follows an evolutionary time sequence, triggered by tidal compression during their preceding pericentre passage. Given that these clouds are the most likely precursors to a generation of massive stars and extreme star clusters, this scenario would have profound implications for constraining the time-evolution of star formation. In this Letter, we search for the initial conditions of the protocluster clouds, focusing on the kinematics of gas situated upstream from pericentre. We observe a highly-regular corrugated velocity field in $\{l,\,v_{\rm LSR}\}$ space, with amplitude and wavelength $A=3.7\,\pm\,0.1$ kms$^{-1}$ and $\lambda_{\rm vel, i}=22.5\,\pm\,0.1$ pc, respectively. The extremes in velocity correlate with a series of massive ($\sim10^{4}$M$_{\odot}$) and compact ($R_{\rm eq}\sim2$ pc), quasi-regularly spaced ($\sim8$ pc), molecular clouds. The corrugation wavelength and cloud separation closely agree with the predicted Toomre ($\sim17$ pc) and Jeans ($\sim6$ pc) lengths, respectively. We conclude that gravitational instabilities are driving the condensation of molecular clouds within the Galactic Centre gas stream. Furthermore, we speculate these seeds are the historical analogue of the dust-ridge molecular clouds, representing the initial conditions of star and cluster formation in the CMZ.

The central molecular zone (CMZ) hosts some of the most massive and dense molecular clouds and star clusters in the Galaxy, offering an important window into star formation under extreme conditions. Star formation in this extreme environment may be closely linked to the 3-D distribution and orbital dynamics of the gas. Here I discuss how our new, accurate description of the $\{l,b,v\}$ structure of the CMZ is helping to constrain its 3-D geometry. I also present the discovery of a highly-regular, corrugated velocity field located just upstream from the dust ridge molecular clouds (which include G0.253+0.016 and Sgr B2). The extremes in this velocity field correlate with a series of massive ($\sim10^{4}$ M$_{\odot}$) cloud condensations. The corrugation wavelength ($\sim23$ pc) and cloud separation ($\sim8$ pc) closely agree with the predicted Toomre ($\sim17$ pc) and Jeans ($\sim6$ pc) lengths, respectively. I conclude that gravitational instabilities are driving the formation of molecular clouds within the Galactic Centre gas stream. Furthermore, I suggest that these seeds are the historical analogues of the dust ridge molecular clouds - possible progenitors of some of the most massive and dense molecular clouds in the Galaxy. If our current best understanding for the 3-D geometry of this system is confirmed, these clouds may pinpoint the beginning of an evolutionary sequence that can be followed, in time, from cloud condensation to star formation.

High-mass stars shape the interstellar medium in galaxies, and yet, largely because the initial conditions are poorly constrained, we do not know how they form. One possibility is that high-mass stars and star clusters form at the junction of filamentary networks, referred to as 'hubs'. In this letter we present the complex anatomy of a protocluster hub within an Infrared Dark Cloud (IRDC), G035.39-00.33, believed to be in an early phase of its evolution. We use high-angular resolution ($\{\theta_{\rm maj}, \theta_{\rm min}\}=\{1.4{\rm arcsec}, 0.8{\rm arcsec}\}\sim\{0.02{\rm pc}, 0.01{\rm pc}\}$) and high-sensitivity ($0.2$mJy beam$^{-1}$; $\sim0.2$M$_{\odot}$) 1.07 mm dust continuum observations from the Atacama Large Millimeter Array (ALMA) to identify a network of narrow, $0.028\,\pm\,0.005$pc wide, filamentary structures. These are a factor of $\gtrsim3$ narrower than the proposed 'quasi-universal' $\sim0.1$pc width of interstellar filaments. Additionally, 28 compact objects are reported, spanning a mass range $0.3{\rm M_{\odot}}<M_{\rm c}<10.4{\rm M_{\odot}}$. This indicates that at least some low-mass objects are forming coevally with more massive counterparts. Comparing to the popular 'bead-on-a-string' analogy, the protocluster hub is poorly represented by a monolithic clump embedded within a single filament. Instead, it comprises multiple intra-hub filaments, each of which retains its integrity as an independent structure and possesses its own embedded core population.

We present 3.7 arcsec (~0.05 pc) resolution 3.2 mm dust continuum observations from the IRAM PdBI, with the aim of studying the structure and fragmentation of the filamentary Infrared Dark Cloud G035.39-00.33. The continuum emission is segmented into a series of 13 quasi-regularly spaced (~0.18pc) cores, following the major axis of the IRDC. We compare the spatial distribution of the cores with that predicted by theoretical work describing the fragmentation of hydrodynamic fluid cylinders, finding a significant (a factor of ~8) discrepancy between the two. Our observations are consistent with the picture emerging from kinematic studies of molecular clouds suggesting that the cores are harboured within a complex network of independent sub-filaments. This result emphasizes the importance of considering the underlying physical structure, and potentially, dynamically important magnetic fields, in any fragmentation analysis. The identified cores exhibit a range in (peak) beam-averaged column density ($3.6{\rm x}10^{23}{\rm cm}^{-2}<N_{H,c}<8.0{\rm x}10^{23}{\rm cm}^{-2}$), mass ($8.1M_{\odot}<M_{c}<26.1M_{\odot}$), and number density ($6.1{\rm x}10^{5}{\rm cm}^{-3}<n_{H, c, eq}<14.7{\rm x}10^{5}{\rm cm}^{-3}$). Two of these cores, dark in the mid-infrared, centrally-concentrated, monolithic (with no traceable substructure at our PdBI resolution), and with estimated masses of the order ~20-25$M_{\odot}$, are good candidates for the progenitors of intermediate-to-high-mass stars. Virial parameters span a range $0.2<\alpha_{\rm vir}<1.3$. Without additional support, possibly from dynamically important magnetic fields with strengths of the order 230$\mu$G<B<670$\mu$G, the cores are susceptible to gravitational collapse. These results may imply a multilayered fragmentation process, which incorporates the formation of sub-filaments, embedded cores, and the possibility of further fragmentation.

Infrared Dark Clouds (IRDCs) are cold, dense regions that are usually found within Giant Molecular Clouds (GMCs). Ongoing star formation within IRDCs is typically still deeply embedded within the surrounding molecular gas. Characterising the properties of relatively quiescent IRDCs may therefore help us to understand the earliest phases of the star formation process. Studies of local molecular clouds have revealed that deuterated species are enhanced in the earliest phases of star formation. In this paper we test this towards IRDC G035.39-00.33. We present an 80 arcsec by 140 arcsec map of the J=2-1 transition of N2D+, obtained with the IRAM-30m telescope. We find that N2D+ is widespread throughout G035.39-00.33. Complementary observations of N2H+(1-0) are used to estimate the deuterium fraction, N(N2D+)/N(N2H+). We report a mean deuterium fraction of 0.04+-0.01, with a maximum of 0.09+-0.02. The mean deuterium fraction is ~3 orders of magnitude greater than the interstellar [D]/[H] ratio. High angular resolution observations are required to exclude beam dilution effects of compact deuterated cores. Using chemical modelling, we find that the average observed values of the deuterium fraction are in agreement with an equilibrium deuterium fraction, given the general properties of the cloud. This implies that the IRDC is at least ~3Myr old, which is ~8 times longer than the mean free-fall time of the observed deuterated region.

Using spectral-line observations of HNCO, N2H+, and HNC, we investigate the kinematics of dense gas in the central ~250 pc of the Galaxy. We present SCOUSE (Semi-automated multi-COmponent Universal Spectral-line fitting Engine), a line fitting algorithm designed to analyse large volumes of spectral-line data efficiently and systematically. Unlike techniques which do not account for complex line profiles, SCOUSE accurately describes the {l, b, v_LSR} distribution of CMZ gas, which is asymmetric about Sgr A* in both position and velocity. Velocity dispersions range from 2.6 km/s<\sigma<53.1 km/s. A median dispersion of 9.8 km/s, translates to a Mach number, M_3D>28. The gas is distributed throughout several "streams", with projected lengths ~100-250 pc. We link the streams to individual clouds and sub-regions, including Sgr C, the 20 and 50 km/s clouds, the dust ridge, and Sgr B2. Shell-like emission features can be explained by the projection of independent molecular clouds in Sgr C and the newly identified conical profile of Sgr B2 in {l ,b, v_LSR} space. These features have previously invoked supernova-driven shells and cloud-cloud collisions as explanations. We instead caution against structure identification in velocity-integrated emission maps. Three geometries describing the 3-D structure of the CMZ are investigated: i) two spiral arms; ii) a closed elliptical orbit; iii) an open stream. While two spiral arms and an open stream qualitatively reproduce the gas distribution, the most recent parameterisation of the closed elliptical orbit does not. Finally, we discuss how proper motion measurements of masers can distinguish between these geometries, and suggest that this effort should be focused on the 20 km/s and 50 km/s clouds and Sgr C.

The formation environment of stars in massive stellar clusters is similar to the environment of stars forming in galaxies at a redshift of 1 - 3, at the peak star formation rate density of the Universe. As massive clusters are still forming at the present day at a fraction of the distance to high-redshift galaxies they offer an opportunity to understand the processes controlling star formation and feedback in conditions similar to those in which most stars in the Universe formed. Here we describe a system of massive clusters and their progenitor gas clouds in the centre of the Milky Way, and outline how detailed observations of this system may be able to: (i) help answer some of the fundamental open questions in star formation and (ii) quantify how stellar feedback couples to the surrounding interstellar medium in this high-pressure, high-redshift analogue environment.

Some formation scenarios that have been put forward to explain multiple populations within Globular Clusters (GCs) require that the young massive cluster have large reservoirs of cold gas within them, which is necessary to form future generations of stars. In this paper we use deep observations taken with Atacama Large Millimeter/sub-millimeter Array (ALMA) to assess the amount of molecular gas within 3 young (50-200 Myr) massive (~10^6 Msun) clusters in the Antennae galaxies. No significant CO(3--2) emission was found associated with any of the three clusters. We place upper limits for the molecular gas within these clusters of ~1x10^5 Msun (or <9 % of the current stellar mass). We briefly review different scenarios that propose multiple episodes of star formation and discuss some of their assumptions and implications. Our results are in tension with the predictions of GC formation scenarios that expect large reservoirs of cool gas within young massive clusters at these ages.

Infrared Dark Clouds (IRDCs) are unique laboratories to study the initial conditions of high-mass star and star cluster formation. We present high-sensitivity and high-angular resolution IRAM PdBI observations of N2H+ (1-0) towards IRDC G035.39-00.33. It is found that G035.39-00.33 is a highly complex environment, consisting of several mildly supersonic filaments (sigma_NT/c_s ~1.5), separated in velocity by <1 km s^-1 . Where multiple spectral components are evident, moment analysis overestimates the non-thermal contribution to the line-width by a factor ~2. Large-scale velocity gradients evident in previous single-dish maps may be explained by the presence of substructure now evident in the interferometric maps. Whilst global velocity gradients are small (<0.7 km s^-1 pc^-1), there is evidence for dynamic processes on local scales (~1.5-2.5 km s^-1 pc^-1 ). Systematic trends in velocity gradient are observed towards several continuum peaks. This suggests that the kinematics are influenced by dense (and in some cases, starless) cores. These trends are interpreted as either infalling material, with accretion rates ~(7 \pm 4)x10^-5 M_sun yr^-1 , or expanding shells with momentum ~24 \pm 12 M_sun km s^-1 . These observations highlight the importance of high-sensitivity and high-spectral resolution data in disentangling the complex kinematic and physical structure of massive star forming regions.

Some theories of dense molecular cloud formation involve dynamical environments driven by converging atomic flows or collisions between preexisting molecular clouds. The determination of the dynamics and physical conditions of the gas in clouds at the early stages of their evolution is essential to establish the dynamical imprints of such collisions, and to infer the processes involved in their formation. We present multi-transition 13CO and C18O maps toward the IRDC G035.39-00.33, believed to be at the earliest stages of evolution. The 13CO and C18O gas is distributed in three filaments (Filaments 1, 2 and 3), where the most massive cores are preferentially found at the intersecting regions between them. The filaments have a similar kinematic structure with smooth velocity gradients of ~0.4-0.8 km s-1 pc-1. Several scenarios are proposed to explain these gradients, including cloud rotation, gas accretion along the filaments, global gravitational collapse, and unresolved sub-filament structures. These results are complemented by HCO+, HNC, H13CO+ and HN13C single-pointing data to search for gas infall signatures. The 13CO and C18O gas motions are supersonic across G035.39-00.33, with the emission showing broader linewidths toward the edges of the IRDC. This could be due to energy dissipation at the densest regions in the cloud. The average H2 densities are ~5000-7000 cm-3, with Filaments 2 and 3 being denser and more massive than Filament 1. The C18O data unveils three regions with high CO depletion factors (f_D~5-12), similar to those found in massive starless cores.

Infrared Dark Clouds (IRDCs) host the initial conditions under which massive stars and stellar clusters form. We have obtained high sensitivity and high spectral resolution observations with the IRAM 30m antenna, which allowed us to perform detailed analysis of the kinematics within one IRDC, G035.39-00.33. We focus on the 1-0 and 3-2 transitions of N2H+, C18O (1-0), and make comparison with SiO (2-1) observations and extinction mapping. Three interacting filaments of gas are found. We report large-scale velocity coherence throughout the cloud, evidenced through small velocity gradients and relatively narrow line widths. This suggests that the merging of these filaments is somewhat "gentle", possibly regulated by magnetic fields. This merging of filaments may be responsible for the weak parsec-scale SiO emission detected by Jimenez-Serra et al. 2010, via grain mantle vaporization. A systematic velocity shift between the N2H+ (1-0) and C18O (1-0) gas throughout the cloud of 0.18 +/- 0.04 kms^{-1} is also found, consistent with a scenario of collisions between filaments which is still ongoing. The N2H+ (1-0) is extended throughout the IRDC and it does not only trace dense cores, as found in nearby low-mass star-forming regions. The average H2 number density across the IRDC is ~ 5 x 10^4 cm^{-3}, at least one order of magnitude larger than in nearby molecular clouds where low-mass stars are forming. A temperature gradient perpendicular to the filament is found. From our study, we conclude that G035.39-00.33 (clearly seen in the extinction map and in N2H+) has been formed via the collision between two relatively quiescent filaments with average densities of ~ 5 x 10^3 cm^{-3}, moving with relative velocities of ~ 5 kms^{-1}. The accumulation of material at the merging points started > 1 Myr ago and it is still ongoing.