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 time-scales associated with the various stages of the star formation process remain poorly constrained. This includes the earliest phases of star formation, during which molecular clouds condense out of the atomic interstellar medium. We present the first in a series of papers with the ultimate goal of compiling the first multi-tracer timeline of star formation, through a comprehensive set of evolutionary phases from atomic gas clouds to unembedded young stellar populations. In this paper, we present an empirical determination of the lifetime of atomic clouds using the Uncertainty Principle for Star Formation formalism, based on the de-correlation of H$\alpha$ and H\,{\sc i} emission as a function of spatial scale. We find an atomic gas cloud lifetime of 48$\substack{+13\\-8}$\,Myr. This timescale is consistent with the predicted average atomic cloud lifetime in the LMC (based on galactic dynamics) that is dominated by the gravitational collapse of the mid-plane ISM. We also determine the overlap time-scale for which both H\,{\sc i} and H$\alpha$ emission are present to be very short ($t_{over}<1.7$\,Myr), consistent with zero, indicating that there is a near-to-complete phase change of the gas to a molecular form in an intermediary stage between H\,{\sc i} clouds and H\,{\sc ii} regions. We utilise the time-scales derived in this work to place empirically determined limits on the time-scale of molecular cloud formation. By performing the same analysis with and without the 30 Doradus region included, we find that the most extreme star forming environment in the LMC has little effect on the measured average atomic gas cloud lifetime. By measuring the lifetime of the atomic gas clouds, we place strong constraints on the physics that drives the formation of molecular clouds and establish a solid foundation for the development of a multi-tracer timeline of star formation in the LMC.

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 examine the role of the large-scale galactic-dynamical environment in setting the properties of giant molecular clouds in Milky Way-like galaxies. We perform three high-resolution simulations of Milky Way-like discs with the moving-mesh hydrodynamics code Arepo, yielding a statistical sample of $\sim 80,000$ giant molecular clouds and $\sim 55,000$ HI clouds. We account for the self-gravity of the gas, momentum and thermal energy injection from supernovae and HII regions, mass injection from stellar winds, and the non-equilibrium chemistry of hydrogen, carbon and oxygen. By varying the external gravitational potential, we probe galactic-dynamical environments spanning an order of magnitude in the orbital angular velocity, gravitational stability, mid-plane pressure and the gradient of the galactic rotation curve. The simulated molecular clouds are highly overdense ($\sim 100\times$) and over-pressured ($\sim 25\times$) relative to the ambient interstellar medium. Their gravo-turbulent and star-forming properties are decoupled from the dynamics of the galactic mid-plane, so that the kpc-scale star formation rate surface density is related only to the number of molecular clouds per unit area of the galactic mid-plane. Despite this, the clouds display clear, statistically-significant correlations of their rotational properties with the rates of galactic shearing and gravitational free-fall. We find that galactic rotation and gravitational instability can influence their elongation, angular momenta, and tangential velocity dispersions. The lower pressures and densities of the HI clouds allow for a greater range of significant dynamical correlations, mirroring the rotational properties of the molecular clouds, while also displaying a coupling of their gravitational and turbulent properties to the galactic-dynamical environment.

In this chapter, we review some historical understanding and recent advances on the Initial Mass Function (IMF) and the Core Mass Function (CMF), both in terms of observations and theories. We focus mostly on star formation in clustered environment since this is suggested by observations to be the dominant mode of star formation. The statistical properties and the fragmentation behaviour of turbulent gas is discussed, and we also discuss the formation of binaries and small multiple systems.

We combine data from ALMA and MUSE to study the resolved (~300 pc scale) star formation relation (star formation rate vs. molecular gas surface density) in cluster galaxies. Our sample consists of 9 Fornax cluster galaxies, including spirals, ellipticals, and dwarfs, covering a stellar mass range of ~10^8.8 - 10^11 M_Sun. CO(1-0) and extinction corrected Halpha were used as tracers for the molecular gas mass and star formation rate, respectively. We compare our results with Kennicutt (1998) and Bigiel et al. (2008). Furthermore, we create depletion time maps to reveal small-scale variations in individual galaxies. We explore these further in FCC290, using the 'uncertainty principle for star formation' (Kruijssen & Longmore, 2014a) to estimate molecular cloud lifetimes, which we find to be short (<10 Myr) in this galaxy. Galaxy-averaged depletion times are compared with other parameters such as stellar mass and cluster-centric distance. We find that the star formation relation in the Fornax cluster is close to those from Kennicutt (1998) and Bigiel et al. (2008}), but overlaps mostly with the shortest depletion times predicted by Bigiel et al. (2008). This slight decrease in depletion time is mostly driven by dwarf galaxies with disturbed molecular gas reservoirs close to the virial radius. In FCC90, a dwarf galaxy with a molecular gas tail, we find that depletion times are a factor >~10 higher in its tail than in its stellar body.

The nature and rate of (viscous) angular momentum transport in protoplanetary discs (PPDs) has important consequences for the formation process of planetary systems. While accretion rates onto the central star yield constraints on such transport in the inner regions of a PPD, empirical constraints on viscous spreading in the outer regions remain challenging to obtain. Here we demonstrate a novel method to probe the angular momentum transport at the outer edge of the disc. This method applies to PPDs that have lost a significant fraction of their mass due to thermal winds driven by UV irradiation from a neighbouring OB star. We demonstrate that this external photoevaporation can explain the observed depletion of discs in the 3-5 Myr old $\sigma$ Orionis region, and use our model to make predictions motivating future empirical investigations of disc winds. For populations of intermediate-age PPDs, in viscous models we show that the mass flux outwards due to angular momentum redistribution is balanced by the mass-loss in the photoevaporative wind. A comparison between wind mass-loss and stellar accretion rates therefore offers an independent constraint on viscous models in the outer regions of PPDs.

Star clusters are fundamental units of stellar feedback and unique tracers of their host galactic properties. In this review, we will first focus on their constituents, i.e.\ detailed insight into their stellar populations and their surrounding ionised, warm, neutral, and molecular gas. We, then, move beyond the Local Group to review star cluster populations at various evolutionary stages, and in diverse galactic environmental conditions accessible in the local Universe. At high redshift, where conditions for cluster formation and evolution are more extreme, we are only able to observe the integrated light of a handful of objects that we believe will become globular clusters. We therefore discuss how numerical and analytical methods, informed by the observed properties of cluster populations in the local Universe, are used to develop sophisticated simulations potentially capable of disentangling the genetic map of galaxy formation and assembly that is carried by globular cluster populations.

Molecular clouds are a fundamental ingredient of galaxies: they are the channels that transform the diffuse gas into stars. The detailed process of how they do it is not completely understood. We review the current knowledge of molecular clouds and their substructure from scales $\sim~$1~kpc down to the filament and core scale. We first review the mechanisms of cloud formation from the warm diffuse interstellar medium down to the cold and dense molecular clouds, the process of molecule formation and the role of the thermal and gravitational instabilities. We also discuss the main physical mechanisms through which clouds gather their mass, and note that all of them may have a role at various stages of the process. In order to understand the dynamics of clouds we then give a critical review of the widely used virial theorem, and its relation to the measurable properties of molecular clouds. Since these properties are the tools we have for understanding the dynamical state of clouds, we critically analyse them. We finally discuss the ubiquitous filamentary structure of molecular clouds and its connection to prestellar cores and star formation.

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.

Star clusters form in dense, hierarchically collapsing gas clouds. Bulk kinetic energy is transformed to turbulence with stars forming from cores fed by filaments. In the most compact regions, stellar feedback is least effective in removing the gas and stars may form very efficiently. These are also the regions where, in high-mass clusters, ejecta from some kind of high-mass stars are effectively captured during the formation phase of some of the low mass stars and effectively channeled into the latter to form multiple populations. Star formation epochs in star clusters are generally set by gas flows that determine the abundance of gas in the cluster. We argue that there is likely only one star formation epoch after which clusters remain essentially clear of gas by cluster winds. Collisional dynamics is important in this phase leading to core collapse, expansion and eventual dispersion of every cluster. We review recent developments in the field with a focus on theoretical work.

Star formation is a complex multi-scale phenomenon that is of significant importance for astrophysics in general. Stars and star formation are key pillars in observational astronomy from local star forming regions in the Milky Way up to high-redshift galaxies. From a theoretical perspective, star formation and feedback processes (radiation, winds, and supernovae) play a pivotal role in advancing our understanding of the physical processes at work, both individually and of their interactions. In this review we will give an overview of the main processes that are important for the understanding of star formation. We start with an observationally motivated view on star formation from a global perspective and outline the general paradigm of the life-cycle of molecular clouds, in which star formation is the key process to close the cycle. After that we focus on the thermal and chemical aspects in star forming regions, discuss turbulence and magnetic fields as well as gravitational forces. Finally, we review the most important stellar feedback mechanisms.

Linking globular clusters (GCs) to the assembly of their host galaxies is an overarching goal in GC studies. The inference of tight scaling relations between GC system properties and the mass of both the stellar and dark halo components of their host galaxies are indicative of an intimate physical connection, yet have also raised fundamental questions about how and when GCs form. Specifically, the inferred correlation between the mass of a GC system (M_gc) and the dark matter halo mass (M_halo) of a galaxy has been posited as a consequence of a causal relation between the formation of dark matter mini-haloes and GC formation during the early epochs of galaxy assembly. We present the first results from a new simulation of a cosmological volume (L=34.4cMpc on a side) from the E-MOSAICS suite, which includes treatments of the formation and evolution of GCs within the framework of a detailed galaxy formation model. The simulated M_gc-M_halo relation is linear for halo masses $>5\times10^{11}$ Msun, and is driven by the hierarchical assembly of galaxies, in agreement with previous studies. Below this halo mass, the simulated relation features a downturn, which we show is consistent with observations, and is driven by the underlying stellar mass (M_star)-halo mass relation of galaxies. Our fiducial model reproduces the observed M_gc-M_star relation across the full mass range, which we argue is more physically relevant than the M_gc-M_halo relation. We also explore the physical processes driving the observed constant value of M_gc / M_halo $\sim 5\times10^{-5}$ and find that it is the result of a combination of cluster formation physics and cluster disruption.

Globular clusters (GCs) are typically old, with most having formed at z >~ 2. This makes understanding their birth environments difficult, as they are typically too distant to observe with sufficient angular resolution to resolve GC birth sites. Using 25 cosmological zoom-in simulations of Milky Way-like galaxies from the E-MOSAICS project, with physically-motivated models for star formation, feedback, and the formation, evolution, and disruption of GCs, we identify the birth environments of present-day GCs. We find roughly half of GCs in these galaxies formed in-situ (52.0 +/- 1.0 per cent) between z ~ 2 - 4, in turbulent, high-pressure discs fed by gas that was accreted without ever being strongly heated through a virial shock or feedback. A minority of GCs form during mergers (12.6 +/- 0.6 per cent in major mergers, and 7.2 +/- 0.5 per cent in minor mergers), but we find that mergers are important for preserving the GCs seen today by ejecting them from their natal, high density interstellar medium (ISM), where proto-GCs are rapidly destroyed due to tidal shocks from ISM substructure. This chaotic history of hierarchical galaxy assembly acts to mix the spatial and kinematic distribution of GCs formed through different channels, making it difficult to use observable GC properties to distinguish GCs formed in mergers from ones formed by smooth accretion, and similarly GCs formed in-situ from those formed ex-situ. These results suggest a simple picture of GC formation, in which GCs are a natural outcome of normal star formation in the typical, gas-rich galaxies that are the progenitors of present-day galaxies.

Quenched galaxies are often observed to contain a strong bulge component. The key question is whether this reflects a causal connection - can star formation be quenched dynamically by bulges or the spheroids of early-type galaxies? We systematically investigate the impact of these morphological components on star formation, by performing a suite of hydrodynamical simulations of isolated galaxies containing a spheroid. We vary the bulge mass and scale radius, while the total initial stellar, halo and gas mass are kept constant, with a gas fraction of 5 per cent. In addition, we consider two different sub-grid star formation prescriptions. The first follows most simulations in the literature by assuming a constant star formation efficiency per free-fall time, whereas in the second model it depends on the gas virial parameter, following high-resolution simulations of turbulent fragmentation. Across all simulations, central spheroids increase the gas velocity dispersion towards the galactic centre. This increases the gravitational stability of the gas disc, suppresses fragmentation and star formation, and results in galaxies hosting extremely smooth and quiescent gas discs that fall below the galaxy main sequence. These effects amplify when using the more sophisticated, dynamics-dependent star formation model. Finally, we discover a pronounced relation between the central stellar surface density and star formation rate (SFR), such that the most bulge-dominated galaxies show the strongest deviation from the main sequence. We conclude that the SFR of galaxies is not only set by the balance between accretion and feedback, but carries a (sometimes dominant) dependence on the gravitational potential.

We report ALMA Band 6 continuum observations of 2000 AU resolution toward four massive molecular clouds in the Central Molecular Zone of the Galaxy. To study gas fragmentation, we use the dendrogram method to identify cores as traced by the dust continuum emission. The four clouds exhibit different fragmentation states at the observed resolution despite having similar masses at the cloud scale ($\sim$1--5 pc). Assuming a constant dust temperature of 20 K, we construct core mass functions of the clouds and find a slightly top-heavy shape as compared to the canonical initial mass function, but we note several significant uncertainties that may affect this result. The characteristic spatial separation between the cores as identified by the minimum spanning tree method, $\sim$$10^4$ AU, and the characteristic core mass, 1--7 $M_\odot$, are consistent with predictions of thermal Jeans fragmentation. The three clouds showing fragmentation may be forming OB associations (stellar mass $\sim$$10^3$ $M_\odot$). None of the four clouds under investigation seem to be currently able to form massive star clusters like the Arches and the Quintuplet ($\sim$$10^4$ $M_\odot$), but they may form such clusters by further gas accretion onto the cores.

We present a detailed comparison of the Milky Way (MW) globular cluster (GC) kinematics with the 25 Milky Way-mass cosmological simulations from the E-MOSAICS project. While the MW falls within the kinematic distribution of GCs spanned by the simulations, the relative kinematics of its metal-rich ($[\rm{Fe/H}]>-1.2$) versus metal-poor ($[\rm{Fe/H}]<-1.2$), and inner ($r<8\rm{kpc}$) versus outer ($r>8\rm{kpc}$) populations are atypical for its mass. To understand the origins of these features, we perform a comprehensive statistical analysis of the simulations, and find 18 correlations describing the assembly of $L^*$ galaxies and their dark matter haloes based on their GC population kinematics. The correlations arise because the orbital distributions of accreted and in-situ GCs depend on the masses and accretion redshifts of accreted satellites, driven by the combined effects of dynamical fraction, tidal stripping, and dynamical heating. Because the kinematics of in-situ/accreted GCs are broadly traced by the metal-rich/metal-poor and inner/outer populations, the observed GC kinematics are a sensitive probe of galaxy assembly. We predict that relative to the population of $L^*$ galaxies, the MW assembled its dark matter and stellar mass rapidly through a combination of in-situ star formation, more than a dozen low-mass mergers, and $1.4\pm1.2$ early ($z=3.1\pm1.3$) major merger. The rapid assembly period ended early, limiting the fraction of accreted stars. We conclude by providing detailed quantitative predictions for the assembly history of the MW.

Determining the efficiency with which gas is converted into stars in galaxies requires an accurate determination of the total reservoir of molecular gas mass. However, despite being the most abundant molecule in the Universe, H$_2$ is challenging to detect through direct observations and indirect methods have to be used to estimate the total molecular gas reservoir. These are often based on scaling relations from tracers such as CO or dust, and are generally calibrated in the Milky Way. Yet, evidence that these scaling relations are environmentally dependent is growing. In particular, the commonly used CO-to-H$_2$ conversion factor (X$_{\rm CO}$) is expected to be higher in metal-poor and/or strongly UV-irradiated environments. We use new SOFIA/FIFI-LS observations of far-infrared fine structure lines from the ionised and neutral gas and the Meudon photodissociation region model to constrain the physical properties and the structure of the gas in the massive star-forming region of 30 Doradus in the Large Magellanic Cloud, and determine the spatially resolved distribution of the total reservoir of molecular gas in the proximity of the young massive cluster R136. We compare this value with the molecular gas mass inferred from ground-based CO observations and dust-based estimates to quantify the impact of this extreme environment on commonly used tracers of the molecular gas. We find that the strong radiation field combined with the half-solar metallicity of the surrounding gas are responsible for a large reservoir of "CO-dark" molecular gas, leaving a large fraction of the total H$_2$ gas (> 75%) undetected when adopting a standard X$_{\rm CO}$ factor in this massive star-forming region.

Historically, it has often been asserted that most stars form in compact clusters. In this scenario, present-day gravitationally-unbound OB associations are the result of the expansion of initially gravitationally-bound star clusters. However, this paradigm is inconsistent with recent results, both theoretical and observational, that instead favour a hierarchical picture of star formation in which stars are formed across a continuous distribution of gas densities and most OB associations never were bound clusters. Instead they are formed in-situ as the low-density side of this distribution, rather than as the remnants of expanding clusters. We utilise the second $Gaia$ data release to quantify the degree to which OB associations are undergoing expansion and, therefore, whether OB associations are the product of expanding clusters, or whether they were born in-situ, as the large-scale, globally-unbound associations that we see today. We find that the observed kinematic properties of associations are consistent with highly substructured velocity fields and additionally require some degree of localised expansion from sub-clusters within the association. While most present-day OB associations do exhibit low levels of expansion, there is no significant correlation between radial velocity and radius. Therefore, the large-scale structure of associations is not set by the expansion of clusters, rather it is a relic of the molecular gas cloud from which the association was formed. This finding is inconsistent with a monolithic model of association formation and instead favours a hierarchical model, in which OB associations form in-situ, following the fractal structure of the gas from which they form.

Giant molecular clouds (GMCs) and their stellar offspring are the building blocks of galaxies. The physical characteristics of GMCs and their evolution are tightly connected to galaxy evolution. The macroscopic properties of the interstellar medium propagate into the properties of GMCs condensing out of it, with correlations between e.g. the galactic and GMC scale gas pressures, surface densities and volume densities. That way, the galactic environment sets the initial conditions for star formation within GMCs. After the onset of massive star formation, stellar feedback from e.g. photoionisation, stellar winds, and supernovae eventually contributes to dispersing the parent cloud, depositing energy, momentum and metals into the surrounding medium, thereby changing the properties of galaxies. This cycling of matter between gas and stars, governed by star formation and feedback, is therefore a major driver of galaxy evolution. Much of the recent debate has focused on the durations of the various evolutionary phases that constitute this cycle in galaxies, and what these can teach us about the physical mechanisms driving the cycle. We review results from observational, theoretical, and numerical work to build a dynamical picture of the evolutionary lifecycle of GMC evolution, star formation, and feedback in galaxies.

Feedback from core collapse supernovae (SNe), the final stage of evolution of the most massive stars, has long been a key element in simulations of galaxy formation. In this paper, we examine how simplifying assumptions made in approximating the SN rates along with underlying physical uncertainty in those rates can lead to large variations in the overall evolution of simulated Milky Way-like galaxies. We find that the clustering of star formation is strongly impacted by the delay between star formation and SN feedback. In addition, the choice to use a realistic delay time distribution or instantaneous injection for SN can have a significant impact on the galaxy. These effects appear even when identical sub-grid models are used for coupling SN energy and momentum, and the total SN energy budget is kept constant. In addition, we show that the uncertain minimum SN progenitor mass has a significant impact on the SN energy budget and injection timescale, and can completely change the overall evolution of the galaxy. These underlying uncertainties mean that despite advances in the sub-grid modelling of SN feedback, there are still serious difficulties in constraining the effects of SN feedback. This complicates the task of comparing different simulations to each other, as well as comparing simulations to observations. We conclude by providing practical limits on the parameters of subgrid models for SN feedback, which bound the uncertainty arising from SN progenitor physics for future predictions from galaxy simulations.

Non-thermal particles and high-energy radiation can play a role in the dynamical processes in star-forming regions and provide an important piece of the multiwavelength observational picture of their structure and components. Powerful stellar winds and supernovae in compact clusters of massive stars and OB associations are known to be favourable sites of high-energy particle acceleration and sources of non-thermal radiation and neutrinos. Namely, young massive stellar clusters are likely sources of the PeV (petaelectronvolt) regime cosmic rays (CRs). They can also be responsible for the cosmic ray composition, e.g., 22Ne/20Ne anomalous isotopic ratio in CRs. Efficient particle acceleration can be accompanied by super-adiabatic amplification of the fluctuating magnetic fields in the systems converting a part of kinetic power of the winds and supernovae into the magnetic energy through the CR-driven instabilities. The escape and CR propagation in the vicinity of the sources are affected by the non-linear CR feedback. These effects are expected to be important in starburst galaxies, which produce high-energy neutrinos and gamma-rays. We give a brief review of the theoretical models and observational data on high-energy particle acceleration and their radiation in star-forming regions with young stellar population.

Globular clusters (GCs) formed when the Milky Way experienced a phase of rapid assembly. We use the wealth of information contained in the Galactic GC population to quantify the properties of the satellite galaxies from which the Milky Way assembled. To achieve this, we train an artificial neural network on the E-MOSAICS cosmological simulations of the co-formation and co-evolution of GCs and their host galaxies. The network uses the ages, metallicities, and orbital properties of GCs that formed in the same progenitor galaxies to predict the stellar masses and accretion redshifts of these progenitors. We apply the network to Galactic GCs associated with five progenitors: {\it Gaia}-Enceladus, the Helmi streams, Sequoia, Sagittarius, and the recently discovered, `low-energy' GCs, which provide an excellent match to the predicted properties of the enigmatic galaxy `Kraken'. The five galaxies cover a narrow stellar mass range [$M_\star=(0.6{-}4.6)\times10^8~{\rm M}_\odot$], but have widely different accretion redshifts ($z_{\rm acc}=0.57{-}2.65$). All accretion events represent minor mergers, but Kraken likely represents the most major merger ever experienced by the Milky Way, with stellar and virial mass ratios of $r_{M_\star}=1$:$31^{+34}_{-16}$ and $r_{M_{\rm h}}=1$:$7^{+4}_{-2}$, respectively. The progenitors match the $z=0$ relation between GC number and halo virial mass, but have elevated specific frequencies, suggesting an evolution with redshift. Even though these progenitors likely were the Milky Way's most massive accretion events, they contributed a total mass of only $\log{(M_{\rm \star,tot}/{\rm M}_\odot)}=9.0\pm0.1$, similar to the stellar halo. This implies that the Milky Way grew its stellar mass mostly by in-situ star formation. We conclude by organising these accretion events into the most detailed reconstruction to date of the Milky Way's merger tree.

The ages and metallicities of globular clusters (GCs) are known to be powerful tracers of the properties of their progenitor galaxies, enabling their use in determining the merger histories of galaxies. However, while useful in separating GCs into individual accretion events, the orbits of GC groups themselves have received less attention as probes of their progenitor galaxy properties. In this work, we use simulations of galaxies and their GC systems from the E-MOSAICS project to explore how the present-day orbital properties of GCs are related to the properties of their progenitor galaxies. We find that the orbits of GCs deposited by accretion events are sensitive to the mass and merger redshift of the satellite galaxy. Earlier mergers and larger galaxy masses deposit GCs at smaller median apocentres and lower total orbital energy. The orbital properties of accreted groups of GCs can therefore be used to infer the properties of their progenitor galaxy, though there exists a degeneracy between galaxy mass and accretion time. Combining GC orbits with other tracers (GC ages, metallicities) will help to break the galaxy mass/accretion time degeneracy, enabling stronger constraints on the properties of their progenitor galaxy. In situ GCs generally orbit at lower energies (small apocentres) than accreted GCs, however they exhibit a large tail to high energies and even retrograde orbits (relative to the present-day disc), showing significant overlap with accreted GCs. Applying the results to Milky Way GCs groups suggests a merger redshift $z \sim 1.5$ for the Gaia Sausage/Enceladus and $z>2$ for the `low-energy'/Kraken group, adding further evidence that the Milky Way had two significant mergers in its past.

In recent years, exciting developments have taken place in the identification of the role of cosmic rays in star-forming environments. Observations from radio to infrared wavelengths and theoretical modelling have shown that low-energy cosmic rays (<1 TeV) play a fundamental role in shaping the chemical richness of the interstellar medium, determining the dynamical evolution of molecular clouds. In this review we summarise in a coherent picture the main results obtained by observations and by theoretical models of propagation and generation of cosmic rays, from the smallest scales of protostars and circumstellar discs, to young stellar clusters, up to Galactic and extragalactic scales. We also discuss the new fields that will be explored in the near future thanks to new generation instruments, such as: CTA, for the $\gamma$-ray emission from high-mass protostars; SKA and precursors, for the synchrotron emission at different scales; and ELT/HIRES, JWST, and ARIEL, for the impact of cosmic rays on exoplanetary atmospheres and habitability.

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.

We present a new theoretical framework for using entropy to understand how outflows driven by supernovae are launched from disc galaxies: via continuous, buoyant acceleration through the circumgalactic medium (CGM). When young star clusters detonate supernovae in the interstellar medium (ISM) of a galaxy, they generate hot, diffuse bubbles that push on the surrounding ISM and evaporate that ISM into their interiors. As these bubbles reach the scale height of the ISM, they break out of the disc, rising into the CGM. Once these bubbles break out, if they have sufficiently high entropy, they will feel an upward acceleration, owing to a local buoyant force. This upward force will accelerate these bubbles, driving them to high galactocentric radii, keeping them in the CGM for $>\Gyr$, even if their initial velocity is much lower than the local escape velocity. We derive an equation of motion for these entropy-driven winds that connects the ISM properties, halo mass, and CGM profile of galaxies to the ultimate evolution of feedback-driven winds. We explore the parameter space of these equations, and show how this novel framework can explain both self-consistent simulations of star formation and galactic outflows as well as the new wealth of observations of CGM kinematics. We show that these entropy-driven winds can produce long wind recycling times, while still carrying a significant amount of mass. Comparisons to simulations and observations show entropy-driven winds convincingly explain the kinematics of galactic outflows.

In Meidt et al. (2018), we showed that gas kinematics on the scale of individual molecular clouds are not dominated by self-gravity but also track a component that originates with orbital motion in the potential of the host galaxy. This agrees with observed cloud line widths, which show systematic variations from virial motions with environment, pointing at the influence of the galaxy potential. In this paper, we hypothesize that these motions act to slow down the collapse of gas and so help regulate star formation. Extending the results of Meidt et al. (2018), we derive a dynamical collapse timescale that approaches the free-fall time only once the gas has fully decoupled from the galactic potential. Using this timescale we make predictions for how the fraction of free-falling, strongly self-gravitating gas varies throughout the disks of star-forming galaxies. We also use this collapse timescale to predict variations in the molecular gas star formation efficiency, which is lowered from a maximum, feedback-regulated level in the presence of strong coupling to the galactic potential. Our model implies that gas can only decouple from the galaxy to collapse and efficiently form stars deep within clouds. We show that this naturally explains the observed drop in star formation rate per unit gas mass in the Milky Way's CMZ and other galaxy centers. The model for a galactic bottleneck to star formation also agrees well with resolved observations of dense gas and star formation in galaxy disks and the properties of local clouds.

We present MUSE Integral Field Unit (IFU) observations of five individual HII regions in two giant (>100 pc in radius) star-forming complexes in the low-metallicity ($Z$~0.33 $Z_{\odot}$) nearby (D ~ 2 Mpc) dwarf spiral galaxy NGC 300. We combine the IFU data with high spatial resolution HST photometry to demonstrate the extraction of stellar spectra and the classification of individual stars from ground-based data at the distance of 2 Mpc. For the two star-forming complexes, in which no O-type stars had previously been identified, we find a total of 13 newly identified O-type stars in the mass range 15-50 M$_{\odot}$, as well as 4 Wolf-Rayet stars. We use the derived massive stellar content to analyze the impact of stellar feedback on the HII regions. As already found for HII regions in the Magellanic Clouds, the dynamics of the analyzed NGC 300 HII regions are dominated by a combination of the pressure of the ionized gas and stellar winds. By comparing the derived ionized gas mass loading factors to the total gas mass loading factor across the NGC 300 disk, we find that the latter is an order of magnitude higher, either indicating very early evolutionary stages for these HII regions, or being a direct result of the multi-phase nature of feedback-driven bubbles. Moreover, we analyze the relation between the star formation rate and the pressure of the ionized gas as derived from small (<100 pc) scales, as both quantities are systematically overestimated when derived on galactic scales. With the wealth of ongoing and upcoming IFU instruments and programs, this study serves as a pathfinder for the systematic investigation of resolved stellar feedback in nearby galaxies, and it delivers the necessary analysis tools to enable massive stellar content and feedback studies sampling an unprecedented range of HII region properties across entire galaxies in the nearby Universe.

Young stellar clusters across nearly five orders of magnitude in mass appear to follow a simple mass-radius relationship (MRR), $R_{\star} \propto M_{\star}^{\alpha}$, with $\alpha \approx 0.2 - 0.33$. Here, we develop a simple analytic model to explain this observation. We begin by considering giant molecular clouds near virial equilibrium and subsequently relate the properties of the cluster to those of the cloud. In turn, we relate the cloud properties to those of the large-scale galactic environment. The model predicts an \textit{initial} mass-radius relation of constant surface density, $R_{\star} \propto M_{\star}^{1/2}$. It also predicts the initial cluster radius depends on the large-scale gas density $\Sigma_g$ of the ambient ISM, scaling as $R_{\star} \propto \Sigma_g^{-1/2}$. We argue that the tendency of observed clusters to fall along lines of shallower MRRs than our initial $R_{\star} \propto M_{\star}^{1/2}$ is in fact a combination of two effects. The fact that massive clusters can \textit{only} form in high gas-density environments, when combined with the $R_{\star} \propto \Sigma_g^{-1/2}$ scaling we find here, ultimately shallows the global slope at high masses to nearly $R_{\star} \propto M_{\star}^{1/3}$. Meanwhile, at low masses relaxation-driven expansion quickly shallows the MRR. We combine our predicted MRR with a simple population synthesis model and apply it to a range of star-forming environments, from nearby disc galaxies to nuclear starbursts, and find good agreement throughout. We provide quantitative predictions for the radii of proto-globular clusters (GCs) and discuss the implications of the model for GC evolution and survival across cosmic time as well as the dynamical assembly of black hole binaries in stellar clusters.

The {\alpha}-element abundances of the globular cluster (GC) and field star populations of galaxies encode information about the formation of each of these components. We use the E-MOSAICS cosmological simulations of ~L* galaxies and their GCs to investigate the [{\alpha}/Fe]-[Fe/H] distribution of field stars and GCs in 25 Milky Way-mass galaxies. The [{\alpha}/Fe]-[Fe/H] distribution go GCs largely follows that of the field stars and can also therefore be used as tracers of the [{\alpha}/Fe]-[Fe/H] evolution of the galaxy. Due to the difference in their star formation histories, GCs associated with stellar streams (i.e. which have recently been accreted) have systematically lower [{\alpha}/Fe] at fixed [Fe/H]. Therefore, if a GC is observed to have low [{\alpha}/Fe] for its [Fe/H] there is an increased probability that this GC was accreted recently alongside a dwarf galaxy. There is a wide range of shapes for the field star [{\alpha}/Fe]-[Fe/H] distribution, with a notable subset of galaxies exhibiting bimodal distributions, in which the high [{\alpha}/Fe] sequence is mostly comprised of stars in the bulge, a high fraction of which are from disrupted GCs. We calculate the contribution of disrupted GCs to the bulge component of the 25 simulated galaxies and find values between 0.3-14 per cent, where this fraction correlates with the galaxy's formation time. The upper range of these fractions is compatible with observationally-inferred measurements for the Milky Way, suggesting that in this respect the Milky Way is not typical of L* galaxies, having experienced a phase of unusually rapid growth at early times.

The Central Molecular Zone (CMZ) is the most extreme star-forming environment in the Milky Way in terms of gas pressures, densities, temperatures, and dynamics. It acts as a critical test bed for developing star formation theories applicable to the (high-redshift-like) conditions under which most stars in the Universe formed. We present a set of numerical simulations of molecular clouds orbiting on the 100-pc stream that dominates the molecular gas reservoir of the CMZ, with the goal of characterising their morphological and kinematic evolution in response to the external gravitational potential and their eccentric orbital motion. These simulations capture the evolution of single clouds in a strong and plausibly dominant background potential. We find that the evolution of the clouds is closely coupled to the orbital dynamics and their arrival on the 100-pc stream marks a transformative event in their lifecycle. The clouds' sizes, aspect ratios, position angles, filamentary structure, column densities, velocity dispersions, line-of-sight velocity gradients, angular momenta, and overall kinematic complexity are controlled by the background potential and their passage through the orbit's pericentre. We compare these predictions of our simulations to observations of clouds on the Galactic Centre `dust ridge' and find that the inclusion of galactic dynamics naturally reproduces a surprisingly wide variety of key observed morphological and kinematic features. We argue that the accretion of gas clouds onto the central regions of galaxies, where the rotation curve turns over and the tidal field becomes fully compressive, is likely to lead to their collapse and associated star formation. This can generate an evolutionary progression of cloud collapse with a common zero point. Together, these processes may naturally give rise to the synchronised starbursts observed in numerous galactic nuclei.

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.

Cloud-scale surveys of molecular gas reveal the link between molecular clouds properties and star formation (SF) across a range of galactic environments. Cloud populations in galaxy disks are considered to be representative of the `normal' SF. At high resolution, however, clouds with exceptional gas properties and SF activity may also be observed in normal disk environments. In this paper, we study the brightest cloud traced in CO emission in the disk of NGC628. The cloud is spatially coincident with an extremely bright HII region. We characterize its molecular gas properties and investigate how feedback and large-scale processes influence the properties of the molecular gas. High resolution CO ALMA observations are used to characterize its mass and dynamical state, which are compared to other clouds in NGC628. A LVG analysis is used to constrain the beam-diluted density and temperature of the molecular gas. We analyze the MUSE spectrum using Starburst99 to characterize the young stellar population associated with the HII region. The cloud is massive ($1-2\times10^7$M$_{\odot}$), with a beam-diluted density of $n_{\rm H_2}=5\times10^4$ cm$^{-3}$. It has a low virial parameter, suggesting that its CO emission may be overluminous due to heating by the HII region. A young ($2-4$ Myr), massive $3\times10^{5}$ M$_{\odot}$ stellar population is associated. We argue that the cloud is currently being destroyed by feedback from young massive stars. Due to the cloud's large mass, this phase of the cloud's evolution is long enough for the impact of feedback on the excitation of the gas to be observed. Its high mass may be related to its location at a spiral co-rotation radius, where gas experiences reduced galactic shear compared to other regions of the disk, and receives a sustained inflow of gas that can promote the cloud's mass growth.

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 report one of the first extragalactic observations of electron temperature variations across a spiral arm. Using MUSE mosaic observations of the nearby galaxy NGC 1672, we measure the [N II]5755 auroral line in a sample of 80 HII regions in the eastern spiral arm of NGC1672. We discover systematic temperature variations as a function of distance perpendicular to the spiral arm. The electron temperature is lowest on the spiral arm itself and highest on the downstream side. Photoionization models of different metallicity, pressure, and age of the ionizing source are explored to understand what properties of the interstellar medium drive the observed temperature variations. An azimuthally varying metallicity appears to be the most likely cause of the temperature variations. The electron temperature measurements solidify recent discoveries of azimuthal variations of oxygen abundance based on strong lines, and rule out the possibility that the abundance variations are artefacts of the strong-line calibrations.

The distribution of metals within a galaxy traces the baryon cycle and the buildup of galactic disks, but the detailed gas phase metallicity distribution remains poorly sampled. We have determined the gas phase oxygen abundances for 7,138 HII regions across the disks of eight nearby galaxies using VLT/MUSE optical integral field spectroscopy as part of the PHANGS-MUSE survey. After removing the first order radial gradients present in each galaxy, we look at the statistics of the metallicity offset (Delta O/H) and explore azimuthal variations. Across each galaxy, we find low (sigma=0.03-0.05 dex) scatter at any given radius, indicative of efficient mixing. We compare physical parameters for those HII regions that are 1 sigma outliers towards both enhanced and reduced abundances. Regions with enhanced abundances have high ionization parameter, higher Halpha luminosity, lower Halpha velocity dispersion, younger star clusters and associated molecular gas clouds show higher molecular gas densities. This indicates recent star formation has locally enriched the material. Regions with reduced abundances show increased Halpha velocity dispersions, suggestive of mixing introducing more pristine material. We observe subtle azimuthal variations in half of the sample, but can not always cleanly associate this with the spiral pattern. Regions with enhanced and reduced abundances are found distributed throughout the disk, and in half of our galaxies we can identify subsections of spiral arms with clearly associated metallicity gradients. This suggests spiral arms play a role in organizing and mixing the ISM.

Globular clusters (GCs) have been posited, alongside dwarf galaxies, as significant contributors to the field stellar population of the Galactic halo. In order to quantify their contribution, we examine the fraction of halo stars formed in stellar clusters in the suite of 25 present-day Milky Way-mass cosmological zoom simulations from the E-MOSAICS project. We find that a median of $2.3$ and $0.3$ per cent of the mass in halo field stars formed in clusters and GCs, defined as clusters more massive than $5\times 10^3$ and $10^5~M_{\odot}$, respectively, with the $25$-$75$th percentiles spanning $1.9$-$3.0$ and $0.2$-$0.5$ per cent being caused by differences in the assembly histories of the host galaxies. Under the extreme assumption that no stellar cluster survives to the present day, the mass fractions increase to a median of $5.9$ and $1.8$ per cent. These small fractions indicate that the disruption of GCs plays a sub-dominant role in the build-up of the stellar halo. We also determine the contributed halo mass fraction that would present signatures of light-element abundance variations considered to be unique to GCs, and find that clusters and GCs would contribute a median of $1.1$ and $0.2$ per cent, respectively. We estimate the contributed fraction of GC stars to the Milky Way halo, based on recent surveys, and find upper limits of $2$-$5$ per cent (significantly lower than previous estimates), suggesting that models other than those invoking strong mass loss are required to describe the formation of chemically enriched stellar populations in GCs.

The stellar birth environment can significantly shorten protoplanetary disc (PPD) lifetimes due to the influence of stellar feedback mechanisms. The degree to which these mechanisms suppress the time and mass available for planet formation is dependent on the local far-ultraviolet (FUV) field strength, stellar density, and ISM properties. In this work, we present the first theoretical framework quantifying the distribution of PPD dispersal time-scales as a function of parameters that describe the galactic environment. We calculate the probability density function for FUV flux and stellar density in the solar neighbourhood. In agreement with previous studies, we find that external photoevaporation is the dominant environment-related factor influencing local stellar populations after the embedded phase. Applying our general prescription to the Central Molecular Zone of the Milky Way (i.e. the central ~250 pc), we predict that 90% of PPDs in the region are destroyed within 1 Myr of the dispersal of the parent molecular cloud. Even in such dense environments, we find that external photoevaporation is the dominant disc depletion mechanism over dynamical encounters between stars. PPDs around low-mass stars are particularly sensitive to FUV-induced mass loss, due to a shallower gravitational potential. For stars of mass ~1 $M_\odot$, the solar neighbourhood lies at approximately the highest gas surface density for which PPD dispersal is still relatively unaffected by external FUV photons, with a median PPD dispersal timescale of ~4 Myr. We highlight the key questions to be addressed to further contextualise the significance of the local galactic environment for planet formation.

The cluster M54 lies at the centre of the Sagittarius dwarf spheroidal galaxy, and therefore may be the closest example of a nuclear star cluster. Either in-situ star formation, inspiralling globular clusters, or a combination have been invoked to explain the wide variety of stellar sub-populations in nuclear star clusters. Globular clusters are known to exhibit light element variations, which can be identified using the photometric construct called a chromosome map. In this letter, we create chromosome maps for three distinct age-metallicity sub-populations in the vicinity of M54. We find that the old, metal-poor population shows the signature of light element variations, while the young and intermediate-age metal rich populations do not. We conclude that the nucleus of Sagittarius formed through a combination of in-situ star formation and globular cluster accretion. This letter demonstrates that properly constructed chromosome maps of iron-complex globular clusters can provide insight into the formation locations of the different stellar populations.

We present an analysis of young star clusters (YSCs) that form in the E-MOSAICS cosmological, hydrodynamical simulations of galaxies and their star cluster populations. Through comparisons with observed YSC populations, this work aims to test models for YSC formation and obtain an insight into the formation processes at work in part of the local galaxy population. We find that the models used in E-MOSAICS for the cluster formation efficiency and high-mass truncation of the initial cluster mass function ($M_\mathrm{c,\ast}$) both quantitatively reproduce the observed values of cluster populations in nearby galaxies. At higher redshifts ($z \geq 2$, near the peak of globular cluster formation) we find that, at a constant star formation rate (SFR) surface density, $M_\mathrm{c,\ast}$ is larger than at $z=0$ by a factor of four due to the higher gas fractions in the simulated high-redshift galaxies. Similar processes should be at work in local galaxies, offering a new way to test the models. We find that cluster age distributions may be sensitive to variations in the cluster formation rate (but not SFR) with time, which may significantly affect their use in tests of cluster mass loss. By comparing simulations with different implementations of cluster formation physics, we find that (even partially) environmentally-independent cluster formation is inconsistent with the brightest cluster-SFR and specific luminosity-$\Sigma_\mathrm{SFR}$ relations, whereas these observables are reproduced by the fiducial, environmentally-varying model. This shows that models in which a constant fraction of stars form in clusters are inconsistent with observations.

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.

We present the results of a proof-of-concept experiment which demonstrates that deep learning can successfully be used for production-scale classification of compact star clusters detected in HST UV-optical imaging of nearby spiral galaxies in the PHANGS-HST survey. Given the relatively small and unbalanced nature of existing, human-labelled star cluster datasets, we transfer the knowledge of neural network models for real-object recognition to classify star clusters candidates into four morphological classes. We show that human classification is at the 66%:37%:40%:61% agreement level for the four classes considered. Our findings indicate that deep learning algorithms achieve 76%:63%:59%:70% for a star cluster sample within 4Mpc < D <10Mpc. We tested the robustness of our deep learning algorithms to generalize to different cluster images using the first data obtained by PHANGS-HST of NGC1559, which is more distant at D = 19Mpc, and found that deep learning produces classification accuracies 73%:42%:52%:67%. We furnish evidence for the robustness of these analyses by using two different neural network models for image classification, trained multiple times from the ground up to assess the variance and stability of our results. We quantified the importance of the NUV, U, B, V and I images for morphological classification with our deep learning models, and find that the V-band is the key contributor as human classifications are based on images taken in that filter. This work lays the foundations to automate classification for these objects at scale, and the creation of a standardized dataset.

We investigate how the dynamical state of molecular clouds relates to host galaxy environment, and how this impacts the star formation efficiency in the Milky Way and seven nearby galaxies. We compile measurements of molecular cloud and host galaxy properties and determine mass-weighted mean cloud properties for entire galaxies and distinct subregions within. We find molecular clouds to be in ambient pressure-balanced virial equilibrium, where clouds in gas-rich, molecular-dominated, high-pressure regions are close to self-virialization, whereas clouds in gas-poor, atomic-dominated, low-pressure environments achieve a balance between their internal kinetic pressure and external pressure from the ambient medium. The star formation efficiency per free-fall time of molecular clouds is low ~0.1%-1% and shows systematic variations of 2 dex as a function of the virial parameter and host galactic environment. The trend observed for clouds in low-pressure environments--as the solar neighborhood--is well matched by state-of-the-art turbulence-regulated models of star formation. However, these models substantially overpredict the low observed star formation efficiencies of clouds in high-pressure environments, which suggests the importance of additional physical parameters not yet considered by these models.

We present a simple physical model for the minimum mass of bound stellar clusters as a function of the galactic environment. The model evaluates which parts of a hierarchically-clustered star-forming region remain bound given the time-scales for gravitational collapse, star formation, and stellar feedback. We predict the initial cluster mass functions (ICMFs) for a variety of galaxies and we show that these predictions are consistent with observations of the solar neighbourhood and nearby galaxies, including the Large Magellanic Cloud and M31. In these galaxies, the low minimum cluster mass of $\sim10^2~\rm{M}_{\odot}$ is caused by sampling statistics, representing the lowest mass at which massive (feedback-generating) stars are expected to form. At the high gas density and shear found in the Milky Way's Central Molecular Zone and the nucleus of M82, the model predicts that a mass $>10^2~\rm{M}_{\odot}$ must collapse into a single cluster prior to feedback-driven dispersal, resulting in narrow ICMFs with elevated characteristic masses. We find that the minimum cluster mass is a sensitive probe of star formation physics due to its steep dependence on the star formation efficiency per free-fall time. Finally, we provide predictions for globular cluster (GC) populations, finding a narrow ICMF for dwarf galaxy progenitors at high redshift, which can explain the high specific frequency of GCs at low metallicities observed in Local Group dwarfs like Fornax and WLM. The predicted ICMFs in high-redshift galaxies constitute a critical test of the model, ideally-suited for the upcoming generation of telescopes.

We present the evolution of the rest-frame ultraviolet (UV) properties of the globular cluster (GC) populations and their host galaxies formed in the E-MOSAICS suite of cosmological hydrodynamical simulations. We compute the luminosities of all clusters associated with 25 simulated Milky Way-mass galaxies, discussed in previous works, in the rest-frame UV and optical bands by combining instantaneous cluster properties (age, mass, metallicity) with simple stellar population models, from redshifts $z=0$ to 10. Due to the rapid fading of young stellar populations in the UV, most of the simulated galaxies do not host GCs bright enough to be individually identified in deep Hubble Space Telescope (HST) observations, even in highly magnified systems. The median age of the most UV-luminous GCs is $<10$ Myr (assuming no extinction), increasing to $\gtrsim 100$ Myr for red optical filters. We estimate that these GCs typically only contribute a few per cent of the total UV luminosity of their host galaxies at any epoch. We predict that the number density of UV-bright proto-GCs (or cluster clumps) will peak between redshifts $z=1-3$. In the main progenitors of Milky Way-mass galaxies, $10$-$20$ per cent of the galaxies at redshifts $1 \lesssim z \lesssim 3$ have clusters brighter than $M_{\rm UV} < -15$, and less than $10$ per cent at other epochs. The brightest cluster in the galaxy sample at $z>2$ is typically $M_{\rm UV} \sim -16$, consistent with the luminosities of observed compact, high-redshift sources.

Direct comparisons between galaxy simulations and observations that both reach scales < 100 pc are strong tools to investigate the cloud-scale physics of star formation and feedback in nearby galaxies. Here we carry out such a comparison for hydrodynamical simulations of a Milky Way-like galaxy, including stochastic star formation, HII region and supernova feedback, and chemical post-processing at 8 pc resolution. Our simulation shows excellent agreement with almost all kpc-scale and larger observables, including total star formation rates, radial profiles of CO, HI, and star formation through the galactic disc, mass ratios of the ISM components, both whole-galaxy and resolved Kennicutt-Schmidt relations, and giant molecular cloud properties. However, we find that our simulation does not reproduce the observed de-correlation between tracers of gas and star formation on < 100 pc scales, known as the star formation 'uncertainty principle', which indicates that observed clouds undergo rapid evolutionary lifecycles. We conclude that the discrepancy is driven by insufficiently-strong pre-supernova feedback in our simulation, which does not disperse the surrounding gas completely, leaving star formation tracer emission too strongly associated with molecular gas tracer emission, inconsistent with observations. This result implies that the cloud-scale de-correlation of gas and star formation is a fundamental test for feedback prescriptions in galaxy simulations, one that can fail even in simulations that reproduce all other macroscopic properties of star-forming galaxies.

The physics of star formation and the deposition of mass, momentum, and energy into the interstellar medium by massive stars (`feedback') are the main uncertainties in modern cosmological simulations of galaxy formation and evolution. These processes determine the properties of galaxies, but are poorly understood on the $\lesssim$100 pc scale of individual giant molecular clouds (GMCs) resolved in modern galaxy formation simulations. The key question is why the timescale for depleting molecular gas through star formation in galaxies ($t_{\rm dep}\approx 2$ Gyr) exceeds the dynamical timescale of GMCs by two orders of magnitude. Either most of a GMC's mass is converted into stars over many dynamical times, or only a small fraction turns into stars before the GMC is dispersed on a dynamical timescale. Here we report our observation that molecular gas and star formation are spatially de-correlated on GMC scales in the nearby flocculent spiral galaxy NGC300, contrary to their tight correlation on galactic scales. We demonstrate that this de-correlation implies rapid evolutionary cycling between GMCs, star formation, and feedback. We apply a novel statistical method to quantify the evolutionary timeline and find that star formation is regulated by efficient stellar feedback, driving GMC dispersal on short timescales (~1.5 Myr) due to radiation and stellar winds, prior to supernova explosions. This feedback limits GMC lifetimes to about one dynamical timescale (~10 Myr), with integrated star formation efficiencies of only 2-3%. Our findings reveal that galaxies consist of building blocks undergoing vigorous, feedback-driven lifecycles, that vary with the galactic environment and collectively define how galaxies form stars. Systematic applications of this multi-scale analysis to large galaxy samples will provide key input for a predictive, bottom-up theory of galaxy formation and evolution.

Understanding the evolution of stellar clusters in an evolving tidal field is critical for studying the disruption of stellar clusters in a cosmological context. We systematically characterise the response of stellar clusters to tidal shocks using controlled N-body simulations of clusters with various properties that are subjected to different types of shocks. We find that the strength of the shock and the density of the cluster within the half-mass radius are the dominant properties that drive the amount of mass lost by the cluster, with the shape of the cluster profile being of minor influence. When the shock is applied as two separate sub-shocks, the amount of mass loss during the second sub-shock is sensitive to the gap time between them. Clusters that experience successive sub-shocks separated by less than their crossing time attain the same masses and sizes at the end of the simulation. However, clusters subjected to sub-shocks separated by more than a crossing time experience different evolutionary histories. The amount of mass lost in the N-body models and its scaling with shock and cluster properties differs from that predicted by classical tidal disruption theory. We demonstrate that the discrepancy is alleviated by including a dependence on the escape time-scale of unbound stars, analogously to mass loss driven by two-body relaxation. With our new theoretical model for shock-driven mass loss, the predicted relative amounts of mass loss agree with the results of the N-body simulations to ~0.3 dex across the full suite of simulations.

The formation histories of globular clusters (GCs) are a key diagnostic for understanding their relation to the evolution of the Universe through cosmic time. We use the suite of 25 cosmological zoom-in simulations of present-day Milky Way-mass galaxies from the E-MOSAICS project to study the formation histories of stars, clusters, and GCs, and how these are affected by the environmental dependence of the cluster formation physics. We find that the median lookback time of GC formation in these galaxies is ${\sim}10.73~$Gyr ($z=2.1$), roughly $2.5~$Gyr earlier than that of the field stars (${\sim}8.34~$Gyr or $z=1.1$). The epoch of peak GC formation is mainly determined by the time evolution of the maximum cluster mass, which depends on the galactic environment and largely increases with the gas pressure. Different metallicity subpopulations of stars, clusters and GCs present overlapping formation histories, implying that star and cluster formation represent continuous processes. The metal-poor GCs ($-2.5<[\rm Fe/H]<-1.5$) of our galaxies are older than the metal-rich GC subpopulation ($-1.0<[\rm Fe/H]<-0.5$), forming $12.13~$Gyr and $10.15~$Gyr ago ($z=3.7$ and $z=1.8$), respectively. The median ages of GCs are found to decrease gradually with increasing metallicity, which suggests different GC metallicity subpopulations do not form independently and their spatial and kinematic distributions are the result of their evolution in the context of hierarchical galaxy formation and evolution. We predict that proto-GC formation is most prevalent at $2\lesssim z \lesssim 3$, which could be tested with observations of lensed galaxies using JWST.

In the local Universe, globular clusters (GCs) with metallicities $[{\rm Fe}/{\rm H}]<-2.5$ are extremely rare. In this Letter, the close connection between GC formation and galaxy evolution is used to show that this GC metallicity `floor' results from the galaxy mass-metallicity relation of ultra low-luminosity galaxies (ULLGs) at high redshift, where the most metal-poor GCs must have formed. Galaxies with metallicities $[{\rm Fe}/{\rm H}]\lesssim-2.5$ have too low masses to form GCs with initial masses $M_{\rm i}\gtrsim10^5~{\rm M}_\odot$, needed to survive for a Hubble time. This translates the galaxy mass-metallicity relation into a maximum initial cluster mass-metallicity relation for $[{\rm Fe}/{\rm H}]\lesssim-1.8$, which naturally leads to the observed colour-magnitude relation of metal-poor GCs at $z=0$ (the `blue tilt'). Its strength traces the slope of the gas phase mass-metallicity relation of ULLGs. Based on the observed blue tilt of GCs in the Virgo and Fornax Clusters, the galaxy mass-metallicity relation is predicted to have a slope of $\alpha=0.4\pm0.1$ for $10^5\lesssim M_\star/{\rm M}_\odot\lesssim10^7$ at $z\gtrsim2$. The GC metallicity floor implies a minimum host galaxy mass and a maximum redshift for GC formation. Any proto-GCs that may be detected at $z>9$ are most likely to end up in galaxies presently more massive than the Milky Way, whereas GCs in low-mass galaxies such as the Fornax dSph ($M_\star\approx4\times10^7~{\rm M}_\odot$) formed at $z\lesssim3$.

It remains unclear what sets the efficiency with which molecular gas transforms into stars. Here we present a new VLA map of the spiral galaxy M51 in 33GHz radio continuum, an extinction-free tracer of star formation, at 3" scales (~100pc). We combined this map with interferometric PdBI/NOEMA observations of CO(1-0) and HCN(1-0) at matched resolution for three regions in M51 (central molecular ring, northern and southern spiral arm segments). While our measurements roughly fall on the well-known correlation between total infrared and HCN luminosity, bridging the gap between Galactic and extragalactic observations, we find systematic offsets from that relation for different dynamical environments probed in M51, e.g. the southern arm segment is more quiescent due to low star formation efficiency (SFE) of the dense gas, despite having a high dense gas fraction. Combining our results with measurements from the literature at 100pc scales, we find that the SFE of the dense gas and the dense gas fraction anti-correlate and correlate, respectively, with the local stellar mass surface density. This is consistent with previous kpc-scale studies. In addition, we find a significant anti-correlation between the SFE and velocity dispersion of the dense gas. Finally, we confirm that a correlation also holds between star formation rate surface density and the dense gas fraction, but it is not stronger than the correlation with dense gas surface density. Our results are hard to reconcile with models relying on a universal gas density threshold for star formation and suggest that turbulence and galactic dynamics play a major role in setting how efficiently dense gas converts into stars.

We present 25 cosmological zoom-in simulations of Milky Way-mass galaxies in the `MOdelling Star cluster population Assembly In Cosmological Simulations within EAGLE' (E-MOSAICS) project. E-MOSAICS couples a detailed physical model for the formation, evolution, and disruption of star clusters to the EAGLE galaxy formation simulations. This enables following the co-formation and co-evolution of galaxies and their star cluster populations, thus realising the long-standing promise of using globular clusters (GCs) as tracers of galaxy formation and assembly. The simulations show that the age-metallicity distributions of GC populations exhibit strong galaxy-to-galaxy variations, resulting from differences in their evolutionary histories. We develop a formalism for systematically constraining the assembly histories of galaxies using GC age-metallicity distributions. These distributions are characterised through 13 metrics that we correlate with 30 quantities describing galaxy formation and assembly (e.g. halo properties, formation/assembly redshifts, stellar mass assembly time-scales, galaxy merger statistics), resulting in 20 statistically (highly) significant correlations. The GC age-metallicity distribution is a sensitive probe of the mass growth, metal enrichment, and minor merger history of the host galaxy. No such relation is found between GCs and major mergers, which play a sub-dominant role in GC formation for Milky Way-mass galaxies. Finally, we show how the GC age-metallicity distribution enables the reconstruction of the host galaxy's merger tree, allowing us to identify all progenitors with masses $M_*\gtrsim10^8$ M$_\odot$ for redshifts $1\leq z\leq2.5$. These results demonstrate that cosmological simulations of the co-formation and co-evolution of GCs and their host galaxies successfully unlock the potential of GCs as quantitative tracers of galaxy formation and assembly.

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.

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 closest galaxy center, our own Central Molecular Zone (CMZ; the central 500 pc of the Milky Way), is a powerful laboratory for studying the secular processes that shape galaxies across cosmic time, from large-scale gas flows and star formation to stellar feedback and interaction with a central supermassive black hole. Research over the last decade has revealed that the process of converting gas into stars in galaxy centers differs from that in galaxy disks. The CMZ is the only galaxy center in which we can identify and weigh individual forming stars, so it is the key location to establish the physical laws governing star formation and feedback under the conditions that dominate star formation across cosmic history. Large-scale surveys of molecular and atomic gas within the inner kiloparsec of the Milky Way (~10 degrees) will require efficient mapping capabilities on single-dish radio telescopes. Characterizing the detailed star formation process will require large-scale, high-resolution surveys of the protostellar populations and small-scale gas structure with dedicated surveys on the Atacama Large Millimeter/submillimeter Array, and eventually with the James Webb Space Telescope, the Next Generation Very Large Array, and the Origins Space Telescope.

We recently reported the discovery of a candidate jet-driving microquasar (S10) in the nearby spiral galaxy NGC 300. However, in the absence of kinematic information, we could not reliably determine the jet power or the dynamical age of the jet cavity. Here, we present optical MUSE integral field unit (IFU) observations of S10, which reveal a bipolar line-emitting jet structure surrounding a continuum-emitting central source. The optical jet lobes of S10 have a total extent of $\sim$ 40 pc and a shock velocity of $\sim$ 150 km s$^{-1}$. Together with the jet kinematics, we exploit the MUSE coverage of the Balmer H$\beta$ line to estimate the density of the surrounding matter and therefore compute the jet power to be $P_{jet}\approx$ 6.3 $\times$ 10$^{38}$ erg s$^{-1}$. An optical analysis of a microquasar jet bubble and a consequent robust derivation of the jet power have been possible only in a handful of similar sources. This study therefore adds valuable insight into microquasar jets, and demonstrates the power of optical integral field spectroscopy in identifying and analysing these objects.

We investigate star formation at very early evolutionary phases in five massive clouds in the inner 500 pc of the Galaxy, the Central Molecular Zone. Using interferometer observations of H$_2$O masers and ultra-compact H II regions, we find evidence of ongoing star formation embedded in cores of 0.2 pc scales and $\gtrsim$10$^5$ cm$^{-3}$ densities. Among the five clouds, Sgr C possesses a high (9%) fraction of gas mass in gravitationally bound and/or protostellar cores, and follows the dense ($\gtrsim$10$^4$ cm$^{-3}$) gas star formation relation that is extrapolated from nearby clouds. The other four clouds have less than 1% of their cloud masses in gravitationally bound and/or protostellar cores, and star formation rates 10 times lower than predicted by the dense gas star formation relation. At the spatial scale of these cores, the star formation efficiency is comparable to that in Galactic disk sources. We suggest that the overall inactive star formation in these Central Molecular Zone clouds could be because there is much less gas confined in gravitationally bound cores, which may be a result of the strong turbulence in this region and/or the very early evolutionary stage of the clouds when collapse has only recently started.

The fraction of stars forming in compact, gravitationally bound clusters (the `cluster formation efficiency' or CFE) is an important quantity for deriving the spatial clustering of stellar feedback and for tracing star formation using stellar clusters across the Universe. Observations of clusters in nearby galaxies have revealed a strong dependence of the CFE on the local gas density, indicating that more stars form in star clusters when the star formation rate surface density is higher. Previously, it has not been possible to test this relation at very young ages and in clusters with individual stars resolved due to the universally-low densities in the cluster-forming regions in the Local Group. This has even led to the suggestion that the CFE increases with distance from the Sun, which would suggest an observational bias. However, the Central Molecular Zone of the Milky Way hosts clouds with densities that are orders of magnitude higher than anywhere else in the Local Group. We report a measurement of the CFE in the highest-density region in the Galaxy, Sgr~B2, based on ALMA observations of high-mass young stellar objects. We find that over a third of the stars ($37\pm7\%$) in Sgr~B2 are forming in bound clusters. This value is consistent with the predictions of environmentally-dependent models for the CFE and is inconsistent with a constant CFE in the Galaxy.

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.

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.

Star formation is spatially clustered across a range of environments, from dense stellar clusters to unbound associations. As a result, radiative or dynamical interactions with neighbouring stars disrupt (proto)planetary systems and limit their radii, leaving a lasting impact on their potential habitability. In the solar neighbourhood, we find that the vast majority of stars form in unbound associations, such that the interaction of (proto)planetary systems with neighbouring stars is limited to the densest sub-regions. However, the fraction of star formation occurring in compact clusters was considerably higher in the past, peaking at ~50% in the young Milky Way at redshift z~2. These results demonstrate that the large-scale star formation environment affects the demographics of planetary systems and the occupation of the habitable zone. We show that planet formation is governed by multi-scale physics, in which Mpc-scale events such as galaxy mergers affect the AU-scale properties of (proto)planetary systems.

We use ALMA and IRAM 30-m telescope data to investigate the relationship between the spectroscopically-traced dense gas fraction and the cloud-scale (120 pc) molecular gas surface density in five nearby, star-forming galaxies. We estimate the dense gas mass fraction at 650 pc and 2800 pc scales using the ratio of HCN (1-0) to CO (1-0) emission. We then use high resolution (120 pc) CO (2-1) maps to calculate the mass-weighted average molecular gas surface density within 650 pc or 2770 pc beam where the dense gas fraction is estimated. On average, the dense gas fraction correlates with the mass-weighted average molecular gas surface density. Thus, parts of a galaxy with higher mean cloud-scale gas surface density also appear to have a larger fraction of dense gas. The normalization and slope of the correlation do vary from galaxy to galaxy and with the size of the regions studied. This correlation is consistent with a scenario where the large-scale environment sets the gas volume density distribution, and this distribution manifests in both the cloud-scale surface density and the dense gas mass fraction.

The evolution of galaxies is governed by equations with chaotic solutions: gravity and compressible hydrodynamics. While this micro-scale chaos and stochasticity has been well studied, it is poorly understood how it couples to macro-scale properties examined in simulations of galaxy formation. In this paper, we show how perturbations introduced by floating-point roundoff, random number generators, and seemingly trivial differences in algorithmic behaviour can produce non-trivial differences in star formation histories, circumgalactic medium (CGM) properties, and the distribution of stellar mass. We examine the importance of stochasticity due to discreteness noise, variations in merger timings and how self-regulation moderates the effects of this stochasticity. We show that chaotic variations in stellar mass can grow until halted by feedback-driven self-regulation or gas exhaustion. We also find that galaxy mergers are critical points from which large (as much as a factor of 2) variations in quantities such as the galaxy stellar mass can grow. These variations can grow and persist for more than a Gyr before regressing towards the mean. These results show that detailed comparisons of simulations require serious consideration of the magnitude of effects compared to run-to-run chaotic variation, and may significantly complicate interpreting the impact of different physical models. Understanding the results of simulations requires us to understand that the process of simulation is not a mapping of an infinitesimal point in configuration space to another, final infinitesimal point. Instead, simulations map a point in a space of possible initial conditions points to a volume of possible final states.

Diffuse emission is observed in galaxies in many tracers across the electromagnetic spectrum, including tracers of star formation, such as H$\alpha$ and ultraviolet (UV), and tracers of gas mass, such as carbon monoxide (CO) transition lines and the 21-cm line of atomic hydrogen (HI). Its treatment is key to extracting meaningful information from observations such as cloud-scale star formation rates. Finally, studying diffuse emission can reveal information about the physical processes taking place in the ISM, such as chemical transitions and the nature of stellar feedback (through the photon escape fraction). We present a physically-motivated method for decomposing astronomical images containing both diffuse emission and compact regions of interest, such as HII regions or molecular clouds, into diffuse and compact component images through filtering in Fourier space. We have previously presented a statistical method for constraining the evolutionary timeline of star formation and mean separation length between compact star forming regions with galaxy-scale observations. We demonstrate how these measurements are biased by the presence of diffuse emission in tracer maps and that by using the mean separation length as a critical length scale to separate diffuse emission from compact emission, we are able to filter out this diffuse emission, thus removing its biasing effect. Furthermore, this method provides, without the need for interferometry or ancillary spectral data, a measurement of the diffuse emission fraction in input tracer maps and decomposed diffuse and compact emission maps for further analysis.

We recently presented a new statistical method to constrain the physics of star formation and feedback on the cloud scale. Fundamentally, this new method is only able to recover the relative durations of different evolutionary phases. It therefore requires a "reference time-scale" to convert the relative time-scales into absolute values. The phase over which the star formation rate (SFR) tracer is visible is the logical choice to assign this reference time-scale to, since it can be characterised using stellar population synthesis models. In this paper, we use hydrodynamical disc galaxy simulations to produce synthetic emission maps of several SFR tracers and apply our statistical method to measure the associated characteristic time-scale of each tracer. These cover 12 ultraviolet (UV) filters (from GALEX, Swift, and HST) covering a wavelength range 150-350 nm, as well as H$\alpha$. For solar-metallicity environments, we find the characteristic time-scales for H$\alpha$ with (without) continuum subtraction to be ~4.3 (6-16) Myr and for the UV filters to be in the range 17-33 Myr, monotonically increasing with wavelength. We find that the characteristic time-scale decreases towards higher metallicities, as well as to lower star formation rate surface densities, if the stellar initial mass function is not well-sampled. We provide fitting functions to the resulting reference time-scales to facilitate observational applications of our statistical method across a wide range of galactic environments. More generally, our results predict the time-scales over which photoionisation and UV heating take place around star-forming regions.

Stellar haloes encode a fossil record of a galaxy's accretion history, generally in the form of structures of low surface brightness, such as stellar streams. While their low surface brightness makes it challenging to determine their age, metallicity, kinematics and spatial structure, the infalling galaxies also deposit globular clusters (GCs) in the halo, which are bright and therefore easier to observe and characterise. To understand how GCs associated with stellar streams can be used to estimate the stellar mass and the infall time of their parent galaxy, we examine a subset of 15 simulations of galaxies and their star clusters from the E-MOSAICS project. E-MOSAICS is a suite of hydrodynamical simulations incorporating a sub-grid model for GC formation and evolution. We find that more massive accreted galaxies typically contribute younger and more metal rich GCs. This lower age results from a more extended cluster formation history in more massive galaxies. In addition, at fixed stellar mass, galaxies that are accreted later host younger clusters, because they can continue to form GCs without being subjected to environmental influences for longer. This explains the large range of ages observed for clusters associated with the Sagittarius dwarf galaxy in the halo of the Milky Way compared to clusters which are thought to have formed in satellites accreted early in the Milky Way's formation history. Using the ages of the GCs associated with the Sagittarius dwarf, we estimate a virial radius crossing lookback time (infall time) of $9.3 \pm 1.8 Gyr$.

We present MUSE integral field data of two HII region complexes in the Large Magellanic Cloud (LMC), N44 and N180. Both regions consist of a main superbubble and a number of smaller, more compact HII regions that formed on the edge of the superbubble. For a total of 11 HII regions, we systematically analyse the radiative and mechanical feedback from the massive O-type stars on the surrounding gas. We exploit the integral field property of the data and the coverage of the HeII$\lambda$5412 line to identify and classify the feedback-driving massive stars, and from the estimated spectral types and luminosity classes we determine the stellar radiative output in terms of the ionising photon flux $Q_{0}$. We characterise the HII regions in terms of their sizes, morphologies, ionisation structure, luminosity and kinematics, and derive oxygen abundances via emission line ratios. We analyse the role of different stellar feedback mechanisms for each region by measuring the direct radiation pressure, the pressure of the ionised gas, and the pressure of the shock-heated winds. We find that stellar winds and ionised gas are the main drivers of HII region expansion in our sample, while the direct radiation pressure is up to three orders of magnitude lower than the other terms. We relate the total pressure to the star formation rate per unit area, $\Sigma_{SFR}$, for each region and find that stellar feedback has a negative effect on star formation, and sets an upper limit to $\Sigma_{SFR}$ as a function of increasing pressure.

Several models have been advanced to explain the multiple stellar populations observed in globular clusters (GCs). Most models necessitate a large initial population of unenriched stars that provide the pollution for an enriched population, and which are subsequently lost from the cluster. This scenario generally requires clusters to lose $>90$ per cent of their birth mass. We use a suite of 25 cosmological zoom-in simulations of present-day Milky Way-mass galaxies from the \emosaics project to study whether dynamical disruption by evaporation and tidal shocking provides the necessary mass loss. We find that GCs with present-day masses $M>10^5~M_{\odot}$ were only $2$-$4$ times more massive at birth, in conflict with the requirements of the proposed models. This factor correlates weakly with metallicity, gas pressure at birth, or galactocentric radius, but increases towards lower GC masses. To reconcile our results with observational data, either an unphysically steep cluster mass-size relation must be assumed, or the initial enriched fractions must be similar to their present values. We provide the required relation between the initial enriched fraction and cluster mass. Dynamical cluster mass loss cannot reproduce the high observed enriched fractions nor their trend with cluster mass.

We use the age-metallicity distribution of 96 Galactic globular clusters (GCs) to infer the formation and assembly history of the Milky Way (MW), culminating in the reconstruction of its merger tree. Based on a quantitative comparison of the Galactic GC population to the 25 cosmological zoom-in simulations of MW-mass galaxies in the E-MOSAICS project, which self-consistently model the formation and evolution of GC populations in a cosmological context, we find that the MW assembled quickly for its mass, reaching $\{25,50\}\%$ of its present-day halo mass already at $z=\{3,1.5\}$ and half of its present-day stellar mass at $z=1.2$. We reconstruct the MW's merger tree from its GC age-metallicity distribution, inferring the number of mergers as a function of mass ratio and redshift. These statistics place the MW's assembly $\textit{rate}$ among the 72th-94th percentile of the E-MOSAICS galaxies, whereas its $\textit{integrated}$ properties (e.g. number of mergers, halo concentration) match the median of the simulations. We conclude that the MW has experienced no major mergers (mass ratios $>$1:4) since $z\sim4$, sharpening previous limits of $z\sim2$. We identify three massive satellite progenitors and constrain their mass growth and enrichment histories. Two are proposed to correspond to Sagittarius (few $10^8~{\rm M}_\odot$) and the GCs formerly associated with Canis Major ($\sim10^9~{\rm M}_\odot$). The third satellite has no known associated relic and was likely accreted between $z=0.6$-$1.3$. We name this enigmatic galaxy $\textit{Kraken}$ and propose that it is the most massive satellite ($M_*\sim2\times10^9~{\rm M}_\odot$) ever accreted by the MW. We predict that $\sim40\%$ of the Galactic GCs formed ex-situ (in galaxies with masses $M_*=2\times10^7$-$2\times10^9~{\rm M}_\odot$), with $6\pm1$ being former nuclear clusters.

The cloud-scale physics of star formation and feedback represent the main uncertainty in galaxy formation studies. Progress is hampered by the limited empirical constraints outside the restricted environment of the Local Group. In particular, the poorly-quantified time evolution of the molecular cloud lifecycle, star formation, and feedback obstructs robust predictions on the scales smaller than the disc scale height that are resolved in modern galaxy formation simulations. We present a new statistical method to derive the evolutionary timeline of molecular clouds and star-forming regions. By quantifying the excess or deficit of the gas-to-stellar flux ratio around peaks of gas or star formation tracer emission, we directly measure the relative rarity of these peaks, which allows us to derive their lifetimes. We present a step-by-step, quantitative description of the method and demonstrate its practical application. The method's accuracy is tested in nearly 300 experiments using simulated galaxy maps, showing that it is capable of constraining the molecular cloud lifetime and feedback time-scale to $<0.1$ dex precision. Access to the evolutionary timeline provides a variety of additional physical quantities, such as the cloud-scale star formation efficiency, the feedback outflow velocity, the mass loading factor, and the feedback energy or momentum coupling efficiencies to the ambient medium. We show that the results are robust for a wide variety of gas and star formation tracers, spatial resolutions, galaxy inclinations, and galaxy sizes. Finally, we demonstrate that our method can be applied out to high redshift ($z\lesssim4$) with a feasible time investment on current large-scale observatories. This is a major shift from previous studies that constrained the physics of star formation and feedback in the immediate vicinity of the Sun.

We present Atacama Large Millimeter/submillimeter Array (ALMA) observations of 12CO(1-0) and 12CO(2-1) in the central 40" (680 pc) of the nuclear starburst galaxy NGC 253, including its molecular outflow. We measure the ratio of brightness temperature for CO(2-1)/CO(1-0), r_21, in the central starburst and outflow-related features. We discuss how r_21 can be used to constrain the optical depth of the CO emission, which impacts the inferred mass of the outflow and consequently the molecular mass outflow rate. We find r_21 less than or equal to 1 throughout, consistent with a majority of the CO emission being optically-thick in the outflow, as it is in the starburst. This suggests that the molecular outflow mass is 3-6 times larger than the lower limit reported for optically thin CO emission from warm molecular gas. The implied molecular mass outflow rate is 25-50 solar masses per year, assuming that conversion factor for the outflowing gas is similar to our best estimates for the bulk of the starburst. This is a factor of 9-19 times larger than the star formation rate in NGC 253. We see tentative evidence for an extended, diffuse CO(2-1) component.

Star formation is a multi-scale process that requires tracing cloud formation and stellar feedback within the local (<kpc) and global galaxy environment. We present first results from two large observing programs on ALMA and VLT/MUSE, mapping cloud scales (1arcsec = 47pc) in both molecular gas and star forming tracers across 90 kpc^2 of the central disk of NGC 628 to probe the physics of star formation. Systematic spatial offsets between molecular clouds and HII regions illustrate the time evolution of star-forming regions. Using uniform sampling of both maps on 50-500 pc scales, we infer molecular gas depletion times of 1-3 Gyr, but also find the increase of scatter in the star formation relation on small scales is consistent with gas and HII regions being only weakly correlated at the cloud (50 pc) scale. This implies a short overlap phase for molecular clouds and HII regions, which we test by directly matching our catalog of 1502 HII regions and 738 GMCs. We uncover only 74 objects in the overlap phase, and we find depletion times >1 Gyr, significantly longer than previously reported for individual star-forming clouds in the Milky Way. Finally, we find no clear trends that relate variations in the depletion time observed on 500 pc scales to physical drivers (metallicity, molecular and stellar mass surface density, molecular gas boundedness) on 50 pc scales.

We present HST/ACS imaging of twenty-three very low surface brightness ($\mu_{e,V}$=25-27.5) galaxies detected in the fields of four nearby galaxy groups. These objects were selected from deep optical imaging obtained with the Dragonfly Telephoto Array. Seven are new discoveries, while most of the others had been identified previously in visual surveys of deep photographic plates and more recent surveys. Few have previously been studied in detail. From the ACS images, we measure distances to the galaxies using both the tip of the red giant branch method and the surface brightness fluctuations method. We demonstrate that the two methods are consistent with each other in the regime where both can be applied. The distances to 15 out of 20 galaxies with stable measurements are consistent with that of the targeted group within errors. This suggests that assuming group membership based solely on projected proximity is ~75% successful in this regime. The galaxies are nearly round, with a median axis ratio of 0.85, and visually resemble dwarf spheroidal galaxies. The objects have a range of sizes, from $R_e$=0.4 kpc to $R_e$=1.8 kpc, with a median of 1.0 kpc. They range in luminosity from $M_V$=-11.4 to $M_V$=-15.6, with a median of -12.4. Galaxies with $R_e$~1 kpc and $M_V$~-12 are fairly rare in the Local Group but we find many of them in this relatively small sample. Four of the objects fall in the class of ultra diffuse galaxies (UDGs), with $R_e$>1.5 kpc and $\mu_{0,V}$>24 mag/arcsec^2, including the recently identified dark matter deficient galaxy NGC1052-DF2.

The metal-poor sub-population of globular cluster (GC) systems exhibits a correlation between the GC average colour and luminosity, especially in those systems associated with massive elliptical galaxies. More luminous (more massive) GCs are typically redder and hence more metal-rich. This 'blue tilt' is often interpreted as a mass-metallicity relation stemming from GC self-enrichment, whereby more massive GCs retain a greater fraction of the enriched gas ejected by their evolving stars, fostering the formation of more metal-rich secondary generations. We examine the E-MOSAICS simulations of the formation and evolution of galaxies and their GC populations, and find that their GCs exhibit a colour-luminosity relation similar to that observed in local galaxies, without the need to invoke mass-dependent self-enrichment. We find that the blue tilt is most appropriately interpreted as a dearth of massive, metal-poor GCs: the formation of massive GCs requires high interstellar gas surface densities, conditions that are most commonly fostered by the most massive, and hence most metal rich, galaxies, at the peak epoch of GC formation. The blue tilt is therefore a consequence of the intimate coupling between the small-scale physics of GC formation and the evolving properties of interstellar gas hosted by hierarchically-assembling galaxies.

We estimate the star formation efficiency per gravitational free fall time, $\epsilon_{\rm ff}$, from observations of nearby galaxies with resolution matched to the typical size of a Giant Molecular Cloud. This quantity, $\epsilon_{\rm ff}$, is theoretically important but so far has only been measured for Milky Way clouds or inferred indirectly in a few other galaxies. Using new, high resolution CO imaging from the PHANGS-ALMA survey, we estimate the gravitational free-fall time at 60 to 120 pc resolution, and contrast this with the local molecular gas depletion time to estimate $\epsilon_{\rm ff}$. Assuming a constant thickness of the molecular gas layer ($H = 100$ pc) across the whole sample, the median value of $\epsilon_{\rm ff}$ in our sample is $0.7\%$. We find a mild scale-dependence, with higher $\epsilon_{\rm ff}$ measured at coarser resolution. Individual galaxies show different values of $\epsilon_{\rm ff}$, with the median $\epsilon_{\rm ff}$ ranging from $0.3\%$ to $2.6\%$. We find the highest $\epsilon_{\rm ff}$ in our lowest mass targets, reflecting both long free-fall times and short depletion times, though we caution that both measurements are subject to biases in low mass galaxies. We estimate the key systematic uncertainties, and show the dominant uncertainty to be the estimated line-of-sight depth through the molecular gas layer and the choice of star formation tracers.

We recently published velocity measurements of luminous globular clusters in the galaxy NGC1052-DF2, concluding that it lies far off the canonical stellar mass - halo mass relation. Here we present a revised velocity for one of the globular clusters, GC-98, and a revised velocity dispersion measurement for the galaxy. We find that the intrinsic dispersion $\sigma=5.6^{+5.2}_{-3.8}$ km/s using Approximate Bayesian Computation, or $\sigma=7.8^{+5.2}_{-2.2}$ km/s using the likelihood. The expected dispersion from the stars alone is ~7 km/s. Responding to a request from the Editors of ApJ Letters and RNAAS, we also briefly comment on the recent analysis of our measurements by Martin et al. (2018).

Gamma ray observations have found evidence of an extremely energetic outflow emanating from the Galactic Centre, and an `excess' of emission at GeV energies towards the Galactic Centre over that expected from current models. Determining whether the outflow is AGN- or star formation-driven, and whether the `excess' is astrophysical in origin or requires new physics (e.g. self-annihilation of dark matter), requires the accurate modelling of the expected energy injection from astrophysical sources and the subsequent interaction with the surrounding environment. We briefly summarise current constraints on the distribution of gas and young stars in the inner few hundred parsec of the Galaxy that can be included in future 2D and 3D modelling of the astrophysical gamma ray emission. The key points to highlight with respect to predominantly axisymmetric models currently in use are: (i) the distribution of dense gas, young stars and interstellar radiation field is highly asymmetric around the Galactic Centre; (ii) star formation is almost exclusively constrained to a Galactocentric radius of ~100pc; (iii) the star formation rate in this region has been constant at <0.1Msun/yr to within a factor of 2 over the last ~5 Myr.

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.

We measure the velocity dispersion, $\sigma$, and surface density, $\Sigma$, of the molecular gas in nearby galaxies from CO spectral line cubes with spatial resolution $45$-$120$ pc, matched to the size of individual giant molecular clouds. Combining $11$ galaxies from the PHANGS-ALMA survey with $4$ targets from the literature, we characterize ${\sim}30,000$ independent sightlines where CO is detected at good significance. $\Sigma$ and $\sigma$ show a strong positive correlation, with the best-fit power law slope close to the expected value for resolved, self-gravitating clouds. This indicates only weak variation in the virial parameter $\alpha_\mathrm{vir}\propto\sigma^2/\Sigma$, which is ${\sim}1.5$-$3.0$ for most galaxies. We do, however, observe enormous variation in the internal turbulent pressure $P_\mathrm{turb}\propto\Sigma\,\sigma^2$, which spans ${\sim}5\rm\;dex$ across our sample. We find $\Sigma$, $\sigma$, and $P_\mathrm{turb}$ to be systematically larger in more massive galaxies. The same quantities appear enhanced in the central kpc of strongly barred galaxies relative to their disks. Based on sensitive maps of M31 and M33, the slope of the $\sigma$-$\Sigma$ relation flattens at $\Sigma\lesssim10\rm\;M_\odot\,pc^{-2}$, leading to high $\sigma$ for a given $\Sigma$ and high apparent $\alpha_\mathrm{vir}$. This echoes results found in the Milky Way, and likely originates from a combination of lower beam filling factors and a stronger influence of local environment on the dynamical state of molecular gas in the low density regime.

We report tentative evidence for a cold stellar stream in the ultra-diffuse galaxy NGC1052-DF2. If confirmed, this stream (which we refer to as "The Maybe Stream") would be the first cold stellar stream detected outside of the Local Group. The candidate stream is very narrow and has an unusual and highly curved shape.

We recently found an ultra diffuse galaxy (UDG) with a half-light radius of R_e = 2.2 kpc and little or no dark matter. The total mass of NGC1052-DF2 was measured from the radial velocities of bright compact objects that are associated with the galaxy. Here we analyze these objects using a combination of HST imaging and Keck spectroscopy. Their average size is <r_h> = 6.2+-0.5 pc and their average ellipticity is <{\epsilon}> = 0.18+-0.02. From a stacked Keck spectrum we derive an age >9 Gyr and a metallicity of [Fe/H] = -1.35+-0.12. Their properties are similar to {\omega} Centauri, the brightest and largest globular cluster in the Milky Way, and our results demonstrate that the luminosity function of metal-poor globular clusters is not universal. The fraction of the total stellar mass that is in the globular cluster system is similar to that in other UDGs, and consistent with "failed galaxy" scenarios where star formation terminated shortly after the clusters were formed. However, the galaxy is a factor of ~1000 removed from the relation between globular cluster mass and total galaxy mass that has been found for other galaxies, including other UDGs. We infer that a dark matter halo is not a prerequisite for the formation of metal-poor globular cluster-like objects in high redshift galaxies.

It has recently been established that the properties of young star clusters (YSCs) can vary as a function of the galactic environment in which they are found. We use the cluster catalogue produced by the Legacy Extragalactic UV Survey (LEGUS) collaboration to investigate cluster properties in the spiral galaxy M51. We analyse the cluster population as a function of galactocentric distance and in arm and inter-arm regions. The cluster mass function exhibits a similar shape at all radial bins, described by a power law with a slope close to $-2$ and an exponential truncation around $10^5\ \rm{M}_{\odot}$ . While the mass functions of the YSCs in the spiral arm and inter-arm regions have similar truncation masses, the inter-arm region mass function has a significantly steeper slope than the one in the arm region; a trend that is also observed in the giant molecular cloud mass function and predicted by simulations. The age distribution of clusters is dependent on the region considered, and is consistent with rapid disruption only in dense regions, while little disruption is observed at large galactocentric distances and in the inter-arm region. The fraction of stars forming in clusters does not show radial variations, despite the drop in the $H_2$ surface density measured as function of galactocentric distance. We suggest that the higher disruption rate observed in the inner part of the galaxy is likely at the origin of the observed flat cluster formation efficiency radial profile.

We propose a simple analytic theory for environmentally-dependent molecular cloud lifetimes, based on the large-scale (galactic) dynamics of the interstellar medium. Within this theory, the cloud lifetime is set by the time-scales for gravitational collapse, galactic shear, spiral arm interactions, epicyclic perturbations and cloud-cloud collisions. It is dependent on five observable quantities, accessible through measurements of the galactic rotation curve, the gas and stellar surface densities, and the gas and stellar velocity dispersions of the host galaxy. We determine how the relative importance of each dynamical mechanism varies throughout the space of observable galactic properties, and conclude that gravitational collapse and galactic shear play the greatest role in setting the cloud lifetime for the considered range of galaxy properties, while cloud-cloud collisions exert a much lesser influence. All five environmental mechanisms are nevertheless required to obtain a complete picture of cloud evolution. We apply our theory to the galaxies M31, M51, M83, and the Milky Way, and find a strong dependence of the cloud lifetime upon galactocentric radius in each case, with a typical cloud lifetime between 10 and 50 Myr. Our theory is ideally-suited for systematic observational tests with the Atacama Large Millimetre/submillimetre array.

We discuss some of the key open questions regarding the formation and evolution of globular clusters (GCs) during galaxy formation and assembly within a cosmological framework. The current state-of-the-art for both observations and simulations is described, and we briefly mention directions for future research. The oldest GCs have ages $\ge$ 12.5 Gyr and formed around the time of reionisation. Resolved colour-magnitude diagrams of Milky Way GCs and direct imaging of lensed proto-GCs at z $\sim$ 6 with JWST promise further insight. Globular clusters are known to host multiple populations of stars with variations in their chemical abundances. Recently, such multiple populations have been detected in $\sim$2 Gyr old compact, massive star clusters. This suggests a common, single pathway for the formation of GCs at high and low redshift. The shape of the initial mass function for GCs remains unknown, however for massive galaxies a power-law mass function is favoured. Significant progress has been made recently modelling GC formation in the context of galaxy formation, with success in reproducing many of the observed GC-galaxy scaling relations.

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.

It is often stated that star clusters are the fundamental units of star formation and that most (if not all) stars form in dense stellar clusters. In this monolithic formation scenario, low density OB associations are formed from the expansion of gravitationally bound clusters following gas expulsion due to stellar feedback. $N$-body simulations of this process show that OB associations formed this way retain signs of expansion and elevated radial anisotropy over tens of Myr. However, recent theoretical and observational studies suggest that star formation is a hierarchical process, following the fractal nature of natal molecular clouds and allowing the formation of large-scale associations in-situ. We distinguish between these two scenarios by characterising the kinematics of OB associations using the Tycho-$Gaia$ Astrometric Solution catalogue. To this end, we quantify four key kinematic diagnostics: the number ratio of stars with positive radial velocities to those with negative radial velocities, the median radial velocity, the median radial velocity normalised by the tangential velocity, and the radial anisotropy parameter. Each quantity presents a useful diagnostic of whether the association was more compact in the past. We compare these diagnostics to models representing random motion and the expanding products of monolithic cluster formation. None of these diagnostics show evidence of expansion, either from a single cluster or multiple clusters, and the observed kinematics are better represented by a random velocity distribution. This result favours the hierarchical star formation model in which a minority of stars forms in bound clusters and large-scale, hierarchically-structured associations are formed in-situ.

Modern extragalactic molecular gas surveys now reach the scales of star-forming giant molecular clouds (GMCs, 20-50 pc). Systematic variations in GMC properties with galaxy environment imply that clouds are not universally self-gravitating objects, decoupled from their surroundings. Here we reexamine the coupling of clouds to their environment and develop a model for 3D gas motions generated by forces arising with the galaxy gravitational potential defined by the background disk of stars and dark matter. We show that these motions can resemble or even exceed the motions needed to support gas against its own self-gravity throughout typical galaxy disks. The importance of the galactic potential in spiral arms and galaxy centers suggests that the response to self-gravity does not always dominate the motions of gas at GMC scales, with implications for observed gas kinematics, virial equilibrium and cloud morphology. We describe how a uniform treatment of gas motions in the plane and in the vertical direction synthesizes the two main mechanisms proposed to regulate star formation: vertical pressure equilibrium and shear/Coriolis forces as parameterized by Toomre Q~1. As the modeled motions are coherent and continually driven by the external potential, they represent support for the gas that is distinct from that conventionally attributed to turbulence, which decays rapidly and requires thus maintenance, e.g. via feedback from star formation. Thus our model suggests the galaxy itself can impose an important limit to star formation, as we explore in a second paper in this series.

We introduce the MOdelling Star cluster population Assembly In Cosmological Simulations within EAGLE (E-MOSAICS) project. E-MOSAICS incorporates models describing the formation, evolution and disruption of star clusters into the EAGLE galaxy formation simulations, enabling the examination of the co-evolution of star clusters and their host galaxies in a fully cosmological context. A fraction of the star formation rate of dense gas is assumed to yield a cluster population; this fraction, and the population's initial properties, are governed by the physical properties of the natal gas. The subsequent evolution and disruption of the entire cluster population is followed accounting for two-body relaxation, stellar evolution, and gravitational shocks induced by the local tidal field. This introductory paper presents a detailed description of the model and initial results from a suite of 10 simulations of $\sim L^\star$ galaxies with disc-like morphologies at $z=0$. The simulations broadly reproduce key observed characteristics of young star clusters and globular clusters (GCs), without invoking separate formation mechanisms for each population. The simulated GCs are the surviving population of massive clusters formed at early epochs ($z\gtrsim1-2$), when the characteristic pressures and surface densities of star-forming gas were significantly higher than observed in local galaxies. We examine the influence of the star formation and assembly histories of galaxies on their cluster populations, finding that (at similar present-day mass) earlier-forming galaxies foster a more massive and disruption-resilient cluster population, while galaxies with late mergers are capable of forming massive clusters even at late cosmic epochs. (Abridged)

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.

We compare the structure of molecular gas at $40$ pc resolution to the ability of gas to form stars across the disk of the spiral galaxy M51. We break the PAWS survey into $370$ pc and $1.1$ kpc resolution elements, and within each we estimate the molecular gas depletion time ($\tau_{\rm Dep}^{\rm mol}$), the star formation efficiency per free fall time ($\epsilon_{\rm ff}$), and the mass-weighted cloud-scale (40 pc) properties of the molecular gas: surface density, $\Sigma$, line width, $\sigma$, and $b\equiv\Sigma/\sigma^2\propto\alpha_{\rm vir}^{-1}$, a parameter that traces the boundedness of the gas. We show that the cloud-scale surface density appears to be a reasonable proxy for mean volume density. Applying this, we find a typical star formation efficiency per free-fall time, $\epsilon_{ff} \left( \left< \Sigma_{40pc} \right> \right) \sim 0.3{-}0.36\%$, lower than adopted in many models and found for local clouds. More, the efficiency per free fall time anti-correlates with both $\Sigma$ and $\sigma$, in some tension with turbulent star formation models. The best predictor of the rate of star formation per unit gas mass in our analysis is $b \equiv \Sigma / \sigma^2$, tracing the strength of self gravity, with $\tau_{\rm Dep}^{\rm mol} \propto b^{-0.9}$. The sense of the correlation is that gas with stronger self-gravity (higher $b$) forms stars at a higher rate (low $\tau_{\rm Dep}^{\rm mol}$). The different regions of the galaxy mostly overlap in $\tau_{\rm Dep}^{\rm mol}$ as a function of $b$, so that low $b$ explains the surprisingly high $\tau_{\rm Dep}^{\rm mol}$ found towards the inner spiral arms found by by Meidt et al. (2013).

Observations of globular clusters (GCs) and field stars in the halos of the giant elliptical galaxy Cen A and the spiral galaxy M31 show a large range of cluster-to-star number ratios ('specific frequencies'). The cluster-to-star ratio decreases with increasing metallicity by a factor of 100-1000, at all galactocentric radii and with a slope that does not seem to depend on radius. In dwarf galaxies, the GCs are also more metal-poor than the field stars on average. These observations indicate a strong dependence of either the cluster formation efficiency or the cluster destruction rate on metallicity and environment. We aim to explain these trends by considering various effects that may influence the observed cluster-to-star ratio as a function of metallicity, environment and cosmological history. We show that both the cluster formation efficiency and the maximum cluster mass increase with metallicity, so they cannot explain the observed trend. Destruction of GCs by tidal stripping and dynamical friction destroy clusters over too small a range of galactocentric radii. We show that cluster destruction by tidal shocks from giant molecular clouds in the high-density formation environments of GCs becomes increasingly efficient towards high galaxy masses and, hence, towards high metallicities. The predicted cluster-to-star ratio decreases by a factor 100-1000 towards high metallicities and should only weakly depend on galactocentric radius due to orbital mixing during hierarchical galaxy merging, consistent with the observations. The observed, strong dependence of the cluster-to-star ratio on metallicity and the independence of its slope on galactocentric radius can be explained by cluster destruction and hierarchical galaxy growth. As a result, we find that the metallicity-dependence of the cluster-to-star ratio does not reflect a GC formation efficiency, but a survival fraction. (Abridged)

The age distributions of stellar cluster populations have long been proposed to probe the recent formation history of the host galaxy. However, progress is hampered by the limited understanding of cluster disruption by evaporation and tidal shocks. We study the age distributions of clusters in smoothed particle hydrodynamics simulations of isolated disc galaxies, which include a self-consistent, physical model for the formation and dynamical evolution of the cluster population and account for the variation of cluster disruption in time and space. We show that the downward slope of the cluster age distribution due to disruption cannot be reproduced with a single functional form, because the disruption rate exhibits systematic trends with cluster age (the `cruel cradle effect'). This problem is resolved by using the median cluster age to trace cluster disruption. Across 120 independent galaxy snapshots and simulated cluster populations, we perform two-dimensional power law fits of the median cluster age to various macroscopic physical quantities and find that it scales as $t_{\rm med}\propto \Sigma^{-0.51\pm0.03}\sigma_{\rm 1D}^{-0.85\pm0.10}M_{\rm min}^\gamma$, for the gas surface density $\Sigma$, gas velocity dispersion $\sigma_{\rm 1D}$, and minimum cluster mass $M_{\rm min}$. This scaling accurately describes observed cluster populations and indicates disruption by impulsive tidal shocks from the interstellar medium. The term $M_{\rm min}^\gamma$ provides a model-independent way to measure the mass dependence of the cluster disruption time $\gamma$. Finally, the ensemble-average cluster lifetime depends on the gas density less strongly than the instantaneous disruption time of single clusters. These results reflect the variation of cluster disruption in time and space. We provide quantitative ways of accounting for these physics in cluster population studies.

High-mass stars have generally been assumed to accrete most of their mass while already contracted onto the main sequence, but this hypothesis has not been observationally tested. We present ALMA observations of a 3 x 1.5 pc area in the W51 high-mass star-forming complex. We identify dust continuum sources and measure the gas and dust temperature through both rotational diagram modeling of CH3OH and brightness-temperature-based limits. The observed region contains three high-mass YSOs that appear to be at the earliest stages of their formation, with no signs of ionizing radiation from their central sources. The data reveal high gas and dust temperatures (T > 100 K) extending out to about 5000 AU from each of these sources. There are no clear signs of disks or rotating structures down to our 1000 AU resolution. The extended warm gas provides evidence that, during the process of forming, these high-mass stars heat a large volume and correspondingly large mass of gas in their surroundings, inhibiting fragmentation and therefore keeping a large reservoir available to feed from. By contrast, the more mature massive stars that illuminate compact H ii regions have little effect on their surrounding dense gas, suggesting that these main sequence stars have completed most or all of their accretion. The high luminosity of the massive protostars (L > 10^4 solar luminosities), combined with a lack of centimeter continuum emission from these sources, implies that they are not on the main sequence while they accrete the majority of their mass; instead, they may be bloated and cool.

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.

We present a simple, self-consistent model to predict the maximum masses of giant molecular clouds (GMCs), stellar clusters and high-redshift clumps as a function of the galactic environment. Recent works have proposed that these maximum masses are set by shearing motions and centrifugal forces, but we show that this idea is inconsistent with the low masses observed across an important range of local-Universe environments, such as low-surface density galaxies and galaxy outskirts. Instead, we propose that feedback from young stars can disrupt clouds before the global collapse of the shear-limited area is completed. We develop a shear-feedback hybrid model that depends on three observable quantities: the gas surface density, the epicylic frequency, and the Toomre parameter. The model is tested in four galactic environments: the Milky Way, the Local Group galaxy M31, the spiral galaxy M83, and the high-redshift galaxy zC406690. We demonstrate that our model simultaneously reproduces the observed maximum masses of GMCs, clumps and clusters in each of these environments. We find that clouds and clusters in M31 and in the Milky Way are feedback-limited beyond radii of 8.4 and 4 kpc, respectively, whereas the masses in M83 and zC406690 are shear-limited at all radii. In zC406690, the maximum cluster masses decrease further due to their inspiral by dynamical friction. These results illustrate that the maximum masses change from being shear-limited to being feedback-limited as galaxies become less gas-rich and evolve towards low shear. This explains why high-redshift clumps are more massive than GMCs in the Local Universe.

We recently reported a population of protostellar candidates in the 20 km s$^{-1}$ cloud in the Central Molecular Zone of the Milky Way, traced by H$_2$O masers in gravitationally bound dense cores. In this paper, we report high-angular-resolution ($\sim$3'') molecular line studies of the environment of star formation in this cloud. Maps of various molecular line transitions as well as the continuum at 1.3 mm are obtained using the Submillimeter Array. Five NH$_3$ inversion lines and the 1.3 cm continuum are observed with the Karl G. Jansky Very Large Array. The interferometric observations are complemented with single-dish data. We find that the CH$_3$OH, SO, and HNCO lines, which are usually shock tracers, are better correlated spatially with the compact dust emission from dense cores among the detected lines. These lines also show enhancement in intensities with respect to SiO intensities toward the compact dust emission, suggesting the presence of slow shocks or hot cores in these regions. We find gas temperatures of $\gtrsim$100 K at 0.1-pc scales based on RADEX modelling of the H$_2$CO and NH$_3$ lines. Although no strong correlations between temperatures and linewidths/H$_2$O maser luminosities are found, in high-angular-resolution maps we notice several candidate shock heated regions offset from any dense cores, as well as signatures of localized heating by protostars in several dense cores. Our findings suggest that at 0.1-pc scales in this cloud star formation and strong turbulence may together affect the chemistry and temperature of the molecular gas.

The work of Adamo et al. (2015) showed that the mass distributions of young massive stellar clusters were truncated above a maximum-mass scale in the nearby galaxy M83 and that this truncation mass varies with galactocentric radius. Here, we present a cloud-based analysis of ALMA CO($1\to 0$) observations of M83 to search for such a truncation mass in the molecular cloud population. We identify a population of 873 molecular clouds in M83 that is largely similar to those found in the Milky Way and Local Group galaxies, though clouds in the centre of the galaxy show high surface densities and enhanced turbulence, as is common for clouds in high-density nuclear environments. Like the young massive clusters, we find a maximum-mass scale for the molecular clouds that decreases radially in the galaxy. We find the most massive young massive cluster tracks the most massive molecular cloud with the cluster mass being $10^{-2}$ times that of the most massive molecular cloud. Outside the nuclear region of M83 ($R_{g}>0.5$ kpc), there is no evidence for changing internal conditions in the population of molecular clouds, with the average internal pressures, densities, and free-fall times remaining constant for the cloud population over the galaxy. This result is consistent with the bound cluster formation efficiency depending only on the large-scale properties of the ISM, rather than the internal conditions of individual clouds.

We present 1.05 mm ALMA observations of the deeply embedded high-mass protocluster G11.92-0.61, designed to search for low-mass cores within the accretion reservoir of the massive protostars. Our ALMA mosaic, which covers an extent of ~0.7 pc at sub-arcsecond (~1400 au) resolution, reveals a rich population of 16 new millimetre continuum sources surrounding the three previously-known millimetre cores. Most of the new sources are located in the outer reaches of the accretion reservoir: the median projected separation from the central, massive (proto)star MM1 is ~0.17 pc. The derived physical properties of the new millimetre continuum sources are consistent with those of low-mass prestellar and protostellar cores in nearby star-forming regions: the median mass, radius, and density of the new sources are 1.3 Msun, 1600 au, and n(H2)~10^7 cm^-3. At least three of the low-mass cores in G11.92-0.61 drive molecular outflows, traced by high-velocity 12CO(3-2) (observed with the SMA) and/or by H2CO and CH3OH emission (observed with ALMA). This finding, combined with the known outflow/accretion activity of MM1, indicates that high- and low-mass stars are forming (accreting) simultaneously within this protocluster. Our ALMA results are consistent with the predictions of competitive-accretion-type models in which high-mass stars form along with their surrounding clusters.

We present the ALMA survey of CO(2-1) emission from the 1/5 solar metallicity, Local Group dwarf galaxy NGC 6822. We achieve high (0.9 arcsec ~ 2 pc) spatial resolution while covering large area: four 250 pc x 250 pc regions that encompass ~2/3 of NGC 6822's star formation. In these regions, we resolve ~150 compact CO clumps that have small radii (~2-3 pc), narrow line width (~1 km/s), and low filling factor across the galaxy. This is consistent with other recent studies of low metallicity galaxies, but here shown with a 15 times larger sample. At parsec scales, CO emission correlates with 8 micron emission better than with 24 micron emission and anti-correlates with Halpha, so that PAH emission may be an effective tracer of molecular gas at low metallicity. The properties of the CO clumps resemble those of similar-size structures in Galactic clouds except of slightly lower surface brightness and CO-to-H2 ratio ~1-2 times the Galactic value. The clumps exist inside larger atomic-molecular complexes with masses typical for giant molecular cloud. Using dust to trace H2 for the entire complex, we find CO-to-H2 to be ~20-25 times the Galactic value, but with strong dependence on spatial scale and variations between complexes that may track their evolutionary state. The H2-to-HI ratio is low globally and only mildly above unity within the complexes. The SFR-to-H2 ratio is ~3-5 times higher in the complexes than in massive disk galaxies, but after accounting for the bias from targeting star-forming regions, we conclude that the global molecular gas depletion time may be as long as in massive disk galaxies.

We present a dynamical model for gas transport, star formation, and winds in the nuclear regions of galaxies, focusing on the Milky Way's Central Molecular Zone (CMZ). In our model angular momentum and mass are transported by a combination of gravitational and bar-driven acoustic instabilities. In gravitationally-unstable regions the gas can form stars, and the resulting feedback drives both turbulence and a wind that ejects mass from the CMZ. We show that the CMZ is in a quasi-steady state where mass deposited at large radii by the bar is transported inward to a star-forming, ring-shaped region at $\sim 100$ pc from the Galactic Centre, where the shear reaches a minimum. This ring undergoes episodic starbursts, with bursts lasting $\sim 5-10$ Myr occurring at $\sim 20-40$ Myr intervals. During quiescence the gas in the ring is not fully cleared, but is driven out of a self-gravitating state by the momentum injected by expanding supernova remnants. Starbursts also drive a wind off the star-forming ring, with a time-averaged mass flux comparable to the star formation rate. We show that our model agrees well with the observed properties of the CMZ, and places it near a star formation minimum within the evolutionary cycle. We argue that such cycles of bursty star formation and winds should be ubiquitous in the nuclei of barred spiral galaxies, and show that the resulting distribution of galactic nuclei on the Kennicutt-Schmidt relation is in good agreement with that observed in nearby galaxies.

We explore the use of mm-wave emission line ratios to trace molecular gas density when observations integrate over a wide range of volume densities within a single telescope beam. For observations targeting external galaxies, this case is unavoidable. Using a framework similar to that of Krumholz and Thompson (2007), we model emission for a set of common extragalactic lines from lognormal and power law density distributions. We consider the median density of gas producing emission and the ability to predict density variations from observed line ratios. We emphasize line ratio variations, because these do not require knowing the absolute abundance of our tracers. Patterns of line ratio variations have the prospect to illuminate the high-end shape of the density distribution, and to capture changes in the dense gas fraction and median volume density. Our results with and without a high density power law tail differ appreciably; we highlight better knowledge of the PDF shape as an important area. We also show the implications of sub-beam density distributions for isotopologue studies targeting dense gas tracers. Differential excitation often implies a significant correction to the naive case. We provide tabulated versions of many of our results, which can be used to interpret changes in mm-wave line ratios in terms of changes in the underlying density distributions.

Distance uncertainties plague our understanding of the physical scales relevant to the physics of star formation in extragalactic studies. The planetary nebulae luminosity function (PNLF) is one of very few techniques that can provide distance estimates to within ~10%, however it requires a planetary nebula (PN) sample that is uncontaminated by other ionizing sources. We employ optical IFU spectroscopy using MUSE on the VLT to measure [OIII] line fluxes for sources unresolved on 50 pc scales within the central star-forming galaxy disk of NGC 628. We use diagnostic line ratios to identify 62 PNe, 30 supernova remnants and 87 HII regions within our fields. Using the 36 brightest PNe we determine a new PNLF distance modulus of 29.91^{+0.08}_{-0.13} mag (9.59^{+0.35}_{-0.57} Mpc), in good agreement with literature values but significantly larger than the previously reported PNLF distance. We are able to explain the discrepancy and recover the previous result when we reintroduce SNR contaminants to our sample. This demonstrates the power of full spectral information over narrowband imaging in isolating PNe. Given our limited spatial coverage within the galaxy, we show that this technique can be used to refine distance estimates even when IFU observations cover only a fraction of a galaxy disk.

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.

We infer the absolute time dependence of kinematic gas temperature along a proposed orbit of molecular clouds in the Central Molecular Zone (CMZ) of the Galactic Center (GC). Ammonia gas temperature maps are one of the results of the "Survey of Water and Ammonia in the Galactic Center" (SWAG, PI: J. Ott); the dynamical model of molecular clouds in the CMZ was taken from Kruijssen et al. (2015). We find that gas temperatures increase as a function of time in both regimes before and after the cloud passes pericenter on its orbit in the GC potential. This is consistent with the recent proposal that pericenter passage triggers gravitational collapse. Other investigated quantities (line width, column density, opacity) show no strong sign of time dependence but are likely dominated by cloud-to-cloud variations.

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.

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.

The environment within the inner few hundred parsecs of the Milky Way, known as the "Central Molecular Zone" (CMZ), harbours densities and pressures orders of magnitude higher than the Galactic Disc; akin to that at the peak of cosmic star formation (Kruijssen & Longmore 2013). Previous studies have shown that current theoretical star-formation models under-predict the observed level of star-formation (SF) in the CMZ by an order of magnitude given the large reservoir of dense gas it contains. Here we explore potential reasons for this apparent dearth of star formation activity.

The Central Molecular Zone (CMZ, the central 500 pc of the Milky Way) contains the largest reservoir of high-density molecular gas in the Galaxy, but forms stars at a rate 10-100 times below commonly-used star formation relations. We discuss recent efforts in understanding how the nearest galactic nucleus forms its stars. The latest models of the gas inflow, star formation, and feedback duty cycle reproduce the main observable features of the CMZ, showing that star formation is episodic and that the CMZ currently resides at a star formation minimum. Using orbital modelling, we derive the three-dimensional geometry of the CMZ and show how the orbital dynamics and the star formation potential of the gas are closely coupled. We discuss how this coupling reveals the physics of star formation and feedback under the conditions seen in high-redshift galaxies, and promotes the formation of the densest stellar clusters in the Galaxy.

The Central Molecular Zone (CMZ), usually referring to the inner 500 pc of the Galaxy, contains a dozen of massive ($\sim10^5$ $M_\odot$) molecular clouds. Are these clouds going to actively form stars like Sgr B2? How are they affected by the extreme physical conditions in the CMZ, such as strong turbulence? Here we present a first step towards answering these questions. Using high-sensitivity, high angular resolution radio and (sub)millimeter observations, we studied deeply embedded star formation in six massive clouds in the CMZ, including the 20 and 50 km s$^{-1}$ clouds, Sgr B1 off (as known as dust ridge clouds e/f), Sgr C, Sgr D, and G0.253-0.016. The VLA water maser observations suggest a population of deeply embedded protostellar candidates, many of which are new detections. The SMA 1.3 mm continuum observations reveal peaks in dust emission associated with the masers, suggesting the existence of dense cores. While our findings confirm that clouds such as G0.253-0.016 lack internal compact substructures and are quiescent in terms of star formation, two clouds (the 20 km s$^{-1}$ cloud and Sgr C) stand out with clusters of water masers with associated dense cores which may suggest a population of deeply embedded protostars at early evolutionary phases. Follow-up observations with VLA and ALMA are necessary to confirm their protostellar nature.

We present ALMA observations of the inner 1' (1.2 kpc) of the Circinus galaxy, the nearest Seyfert. We target CO (1-0) in the region associated with a well-known multiphase outflow driven by the central active galactic nucleus (AGN). While the geometry of Circinus and its outflow make disentangling the latter difficult, we see indications of outflowing molecular gas at velocities consistent with the ionized outflow. We constrain the mass of the outflowing molecular gas to be 1.5e5 -5.1e6 solar masses, yielding a molecular outflow rate of 0.35-12.3 solar masses per year. The values within this range are comparable to the star formation rate in Circinus, indicating that the outflow indeed regulates star formation to some degree. The molecular outflow in Circinus is considerably lower in mass and energetics than previously-studied AGN-driven outflows, especially given its high ratio of AGN luminosity to bolometric luminosity. The molecular outflow in Circinus is, however, consistent with some trends put forth in Cicone et al. (2014), including a linear relation between kinetic power and AGN luminosity, as well as its momentum rate vs. bolometric luminosity (although the latter places Circinus among the starburst galaxies in that sample). We detect additional molecular species including CN and C17O.

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.

The cloud-scale density, velocity dispersion, and gravitational boundedness of the interstellar medium (ISM) vary within and among galaxies. In turbulent models, these properties play key roles in the ability of gas to form stars. New high fidelity, high resolution surveys offer the prospect to measure these quantities across galaxies. We present a simple approach to make such measurements and to test hypotheses that link small-scale gas structure to star formation and galactic environment. Our calculations capture the key physics of the Larson scaling relations, and we show good correspondence between our approach and a traditional "cloud properties" treatment. However, we argue that our method is preferable in many cases because of its simple, reproducible characterization of all emission. Using, low-J 12CO data from recent surveys, we characterize the molecular ISM at 60pc resolution in the Antennae, the Large Magellanic Cloud, M31, M33, M51, and M74. We report the distributions of surface density, velocity dispersion, and gravitational boundedness at 60pc scales and show galaxy-to-galaxy and intra-galaxy variations in each. The distribution of flux as a function of surface density appears roughly lognormal with a 1sigma width of ~0.3 dex, though the center of this distribution varies from galaxy to galaxy. The 60pc resolution line width and molecular gas surface density correlate well, which is a fundamental behavior expected for virialized or free-falling gas. Varying the measurement scale for the LMC and M31, we show that the molecular ISM has higher surface densities, lower line widths, and more self-gravity at smaller scales.

We present new JVLA observations of the high-mass cluster-forming region W51A from 2 to 16 GHz with resolution ${\theta}_{fwhm} \approx$ 0.3 - 0.5". The data reveal a wealth of observational results: (1) Currently-forming, very massive (proto-O) stars are traced by o-H2CO $2_{1,1}-2_{1,2}$ emission, suggesting that this line can be used efficiently as a massive protostar tracer. (2) There is a spatially distributed population of $\sim$mJy continuum sources, including hypercompact H ii regions and candidate colliding wind binaries, in and around the W51 proto-clusters. (3) There are two clearly detected protoclusters, W51e and W51 IRS2, that are gas-rich but may have most of their mass in stars within their inner $\sim$ 0.05 pc. The majority of the bolometric luminosity in W51 most likely comes from a third population of OB stars between these clusters. The presence of a substantial population of exposed O-stars coincident with a population of still-forming massive stars, along with a direct measurement of the low mass loss rate via ionized gas outflow from W51 IRS2, together imply that feedback is ineffective at halting star formation in massive protoclusters. Instead, feedback may shut off the large-scale accretion of diffuse gas onto the W51 protoclusters, implying that they are evolving towards a state of gas exhaustion rather than gas expulsion. Recent theoretical models predict gas exhaustion to be a necessary step in the formation of gravitationally bound stellar clusters, and our results provide an observational validation of this process.

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.

Young massive clusters (YMCs) have central stellar mass surface densities exceeding $10^{4} M_{\odot} pc^{-2}$. It is currently unknown whether the stars formed at such high (proto)stellar densities. We compile a sample of gas clouds in the Galaxy which have sufficient gas mass within a radius of a few parsecs to form a YMC, and compare their radial gas mass distributions to the stellar mass distribution of Galactic YMCs. We find that the gas in the progenitor clouds is distributed differently than the stars in YMCs. The mass surface density profiles of the gas clouds are generally shallower than the stellar mass surface density profiles of the YMCs, which are characterised by prominent dense core regions with radii ~ 0.1 pc, followed by a power-law tail. On the scale of YMC core radii, we find that there are no known clouds with significantly more mass in their central regions when compared to Galactic YMCs. Additionally, we find that models in which stars form from very dense initial conditions require surface densities that are generally higher than those seen in the known candidate YMC progenitor clouds. Our results show that the quiescent, less evolved clouds contain less mass in their central regions than in the highly star-forming clouds. This suggests an evolutionary trend in which clouds continue to accumulate mass towards their centres after the onset of star formation. We conclude that a conveyor-belt scenario for YMC formation is consistent with the current sample of Galactic YMCs and their progenitor clouds.

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.

The colour-magnitude diagrams of some intermediate-age clusters (1-2 Gyr) star clusters show unexpectedly broad main-sequence turnoffs, raising the possibility that these clusters have experienced more than one episode of star formation. Such a scenario predicts the existence of an extended main sequence turn off (eMSTO) only in clusters with escape velocities above a certain threshold ($>15$ km s$^{-1}$), which would allow them to retain or accrete gas that eventually would fuel a secondary extended star-formation episode. This paper presents a test of this scenario based on the study of the young and massive cluster NGC 7252: W3. We use the HST photometry from WFPC2 and WFC3 images obtained with UV and optical filters, as well as MagE echellette spectrograph data from the Las Campanas Clay 6.5m telescope, in order to construct the observed UV/optical SED of NGC 7252: W3. The observations are then compared with synthetic spectra based on different star formation histories consistent with those of the eMSTO clusters. We find that the SED of this cluster is best fitted by a synthetic spectrum with a single stellar population of age $570^{+70}_{-62}$ Myr and mass $1.13^{+0.14}_{-0.13}\times 10^8$ M$_\odot$, confirming earlier works on NGC 7252: W3. We also estimate the lower limit on the central escape velocity of 193 km s$^{-1}$. We rule out extended star-formation histories, like those inferred for the eMSTO clusters in the Magellanic Clouds, at high confidence. We conclude that the escape velocity of a cluster does not dictate whether a cluster can undergo extended periods of star formation.

We have discovered a new H$_2$CO (formaldehyde) $1_{1,0}-1_{1,1}$ 4.82966 GHz maser in Galactic Center Cloud C, G0.38+0.04. At the time of acceptance, this is the eighth region containing an H$_2$CO maser detected in the Galaxy. Cloud C is one of only two sites of confirmed high-mass star formation along the Galactic Center Ridge, affirming that H$_2$CO masers are exclusively associated with high-mass star formation. This discovery led us to search for other masers, among which we found new SiO vibrationally excited masers, making this the fourth star-forming region in the Galaxy to exhibit SiO maser emission. Cloud C is also a known source of CH$_3$OH Class-II and OH maser emission. There are now two known SiO and H$_2$CO maser containing regions in the CMZ, compared to two and six respectively in the Galactic disk, while there is a relative dearth of H$_2$O and CH$_3$OH Class-II masers in the CMZ. SiO and H$_2$CO masers may be preferentially excited in the CMZ, perhaps due to higher gas-phase abundances from grain destruction and heating, or alternatively H$_2$O and CH$_3$OH maser formation may be suppressed in the CMZ. In any case, Cloud C is a new testing ground for understanding maser excitation conditions.

Stars generally form in aggregates, some of which are bound ('clusters') while others are unbound and disperse on short ($\sim10$ Myr) timescales ('associations'). The fraction of stars forming in bound clusters ($\Gamma$) is a fundamental outcome of the star formation process. Recent observational and theoretical work has suggested that $\Gamma$ increases with the gas surface density ($\Sigma$) or star formation rate (SFR) surface density ($\Sigma_{\rm SFR}$), both within galaxies and between different ones. However, a recent paper by Chandar et al. has challenged these results, showing that the $total$ number of stellar aggregates per unit SFR does not vary systematically with the host galaxy's absolute SFR. In this Letter, we show that no variations are expected when no distinction is made between bound and unbound aggregates, because the sum of these two fractions should be close to unity. We also demonstrate that any scaling of $\Gamma$ with the absolute SFR is much weaker than with $\Sigma_{\rm SFR}$, due to the mass-radius-SFR relation of star-forming 'main sequence' galaxies. The environmental variation of $\Gamma$ should therefore be probed as a function of area-normalised quantities, such as $\Sigma$ or $\Sigma_{\rm SFR}$. We present a set of guidelines for meaningful observational tests of cluster formation theories and show that these resolve the reported discrepancy.

We report the discovery of a population of deeply embedded protostellar candidates in the 20 km s$^{-1}$ cloud, one of the massive molecular clouds in the Central Molecular Zone (CMZ) of the Milky Way, using interferometric submillimeter continuum and H$_2$O maser observations. The submillimeter continuum emission shows five 1-pc scale clumps, each of which further fragments into several 0.1-pc scale cores. We identify 17 dense cores, among which 12 are gravitationally bound. Among the 18 H$_2$O masers detected, 13 coincide with the cores and probably trace outflows emanating from the protostars. There are also 5 gravitationally bound dense cores without H$_2$O maser detection. In total the 13 masers and 5 cores may represent 18 protostars with spectral types later than B1 or potential growing more massive stars at earlier evolutionary stage, given the non-detection in the centimeter radio continuum. In combination with previous studies of CH$_3$OH masers, we conclude that the star formation in this cloud is at an early evolutionary phase, before the presence of any significant ionizing or heating sources. Our findings indicate that star formation in this cloud may be triggered by a tidal compression as it approaches pericenter, similar to the case of G0.253+0.016 but with a higher star formation rate, and demonstrate that high angular resolution, high sensitivity maser and submillimeter observations are a promising technique to unveil deeply embedded star formation in the CMZ.

We summarise some of the recent progress in understanding the formation and evolution of globular clusters (GCs) in the context of galaxy formation and evolution. It is discussed that an end-to-end model for GC formation and evolution should capture four different phases: (1) star and cluster formation in the high-pressure interstellar medium of high-redshift galaxies, (2) cluster disruption by tidal shocks in the gas-rich host galaxy disc, (3) cluster migration into the galaxy halo, and (4) the final evaporation-dominated evolution of GCs until the present day. Previous models have mainly focussed on phase 4. We present and discuss a simple model that includes each of these four steps - its key difference with respect to previous work is the simultaneous addition of the high-redshift formation and early evolution of young GCs, as well as their migration into galaxy haloes. The new model provides an excellent match to the observed GC mass spectrum and specific frequency, as well as the relations of GCs to the host dark matter halo mass and supermassive black hole mass. These results show (1) that the properties of present-day GCs are reproduced by assuming that they are the natural outcome of regular high-redshift star formation (i.e. they form according to same physical processes that govern massive cluster formation in the local Universe), and (2) that models only including GC evaporation strongly underestimate their integrated mass loss over a Hubble time.

The formation of globular clusters (GCs) remains one of the main unsolved problems in star and galaxy formation. The past decades have seen important progress in constraining the physics of GC formation from a variety of directions. In this review, we discuss the latest constraints obtained from studies of present-day GC populations, the formation of young massive clusters (YMCs) in the local Universe, and the observed, large-scale conditions for star and cluster formation in high-redshift galaxies. The main conclusion is that the formation of massive, GC progenitor clusters is restricted to high-pressure environments similar to those observed at high redshift and at the sites of YMC formation in the local Universe. However, the correspondingly high gas densities also lead to efficient cluster disruption by impulsive tidal shocks, which limits the survival of GCs progenitor clusters. As a result, the long-term survival of GC progenitor clusters requires them to migrate into the host galaxy halo on a short time-scale. It is proposed that the necessary cluster migration is facilitated by the frequent galaxy mergers occurring at high redshift. We use the available observational and theoretical constraints to condense the current state of the field into a coherent picture of GC formation, in which regular star and cluster formation in high-redshift galaxies naturally leads to the GC populations observed today.

We present an end-to-end, two-phase model for the origin of globular clusters (GCs). In the model, populations of stellar clusters form in the high-pressure discs of high-redshift ($z>2$) galaxies (a rapid-disruption phase due to tidal perturbations from the dense interstellar medium), after which the galaxy mergers associated with hierarchical galaxy formation redistribute the surviving, massive clusters into the galaxy haloes, where they remain until the present day (a slow-disruption phase due to tidal evaporation). The high galaxy merger rates of $z>2$ galaxies allow these clusters to be `liberated' into the galaxy haloes before they are disrupted within the high-density discs. This physically-motivated toy model is the first to include the rapid-disruption phase, which is shown to be essential for simultaneously reproducing the wide variety of properties of observed GC systems, such as their universal characteristic mass-scale, the dependence of the specific frequency on metallicity and galaxy mass, the GC system mass-halo mass relation, the constant number of GCs per unit supermassive black hole mass, and the colour bimodality of GC systems. The model predicts that most of these observables were already in place at $z=1$-$2$, although under rare circumstances GCs may still form in present-day galaxies. In addition, the model provides important constraints on models for multiple stellar populations in GCs by putting limits on initial GC masses and the amount of pristine gas accretion. The paper is concluded with a discussion of these and several other predictions and implications, as well as the main open questions in the field.

The Galactic center is the closest region in which we can study star formation under extreme physical conditions like those in high-redshift galaxies. We measure the temperature of the dense gas in the central molecular zone (CMZ) and examine what drives it. We mapped the inner 300 pc of the CMZ in the temperature-sensitive J = 3-2 para-formaldehyde (p-H$_2$CO) transitions. We used the $3_{2,1} - 2_{2,0} / 3_{0,3} - 2_{0,2}$ line ratio to determine the gas temperature in $n \sim 10^4 - 10^5 $cm$^{-3}$ gas. We have produced temperature maps and cubes with 30" and 1 km/s resolution and published all data in FITS form. Dense gas temperatures in the Galactic center range from ~60 K to > 100 K in selected regions. The highest gas temperatures T_G > 100 K are observed around the Sgr B2 cores, in the extended Sgr B2 cloud, the 20 km/s and 50 km/s clouds, and in "The Brick" (G0.253+0.016). We infer an upper limit on the cosmic ray ionization rate ${\zeta}_{CR} < 10^{-14}$ 1/s. The dense molecular gas temperature of the region around our Galactic center is similar to values found in the central regions of other galaxies, in particular starburst systems. The gas temperature is uniformly higher than the dust temperature, confirming that dust is a coolant in the dense gas. Turbulent heating can readily explain the observed temperatures given the observed line widths. Cosmic rays cannot explain the observed variation in gas temperatures, so CMZ dense gas temperatures are not dominated by cosmic ray heating. The gas temperatures previously observed to be high in the inner ~75 pc are confirmed to be high in the entire CMZ.

We develop a simple dynamical model for the evolution of gas in the centres of barred spiral galaxies, using the Milky Way's Central Molecular Zone (CMZ, i.e., the central few hundred pc) as a case study. We show that, in the presence of a galactic bar, gas in a disc in the central regions of a galaxy will be driven inwards by angular momentum transport induced by acoustic instabilities within the bar's inner Lindblad resonance. This transport process drives turbulence within the gas that temporarily keeps it strongly gravitationally stable and prevents the onset of rapid star formation. However, at some point the rotation curve must transition from approximately flat to approximately solid body, and the resulting reduction in shear reduces the transport rates and causes gas to build up, eventually producing a gravitationally-unstable region that is subject to rapid and violent star formation. For the observed rotation curve of the Milky Way, the accumulation happens $\sim 100$ pc from the centre of the Galaxy, in good agreement with the observed location of gas clouds and young star clusters in the CMZ. The characteristic timescale for gas accumulation and star formation is of order $10-20$ Myr. We argue that similar phenomena should be ubiquitous in other barred spiral galaxies.

Whilst young massive clusters (YMCs; $M$ $\gtrsim$ 10$^{4}$ M$_{\odot}$, age $\lesssim$ 100 Myr) have been identified in significant numbers, their progenitor gas clouds have eluded detection. Recently, four extreme molecular clouds residing within 200 pc of the Galactic centre have been identified as having the properties thought necessary to form YMCs. Here we utilise far-IR continuum data from the Herschel Infrared Galactic Plane Survey (HiGAL) and millimetre spectral line data from the Millimetre Astronomy Legacy Team 90 GHz Survey (MALT90) to determine their global physical and kinematic structure. We derive their masses, dust temperatures and radii and use virial analysis to conclude that they are all likely gravitationally bound -- confirming that they are likely YMC progenitors. We then compare the density profiles of these clouds to those of the gas and stellar components of the Sagittarius B2 Main and North proto-clusters and the stellar distribution of the Arches YMC. We find that even in these clouds -- the most massive and dense quiescent clouds in the Galaxy -- the gas is not compact enough to form an Arches-like ($M$ = 2x10$^{4}$ M$_{\odot}$, R$_{eff}$ = 0.4 pc) stellar distribution. Further dynamical processes would be required to condense the resultant population, indicating that the mass becomes more centrally concentrated as the (proto)-cluster evolves. These results suggest that YMC formation may proceed hierarchically rather than through monolithic collapse.

G0.253+0.016 is a molecular clump that appears to be on the verge of forming a high mass, Arches-like cluster. Here we present new ALMA observations of its small-scale (~0.07 pc) 3mm dust continuum and molecular line emission. The data reveal a complex network of emission features, the morphology of which ranges from small, compact regions to extended, filamentary structures that are seen in both emission and absorption. The dust column density is well traced by molecules with higher excitation energies and critical densities, consistent with a clump that has a denser interior. A statistical analysis supports the idea that turbulence shapes the observed gas structure within G0.253+0.016. We find a clear break in the turbulent power spectrum derived from the optically thin dust continuum emission at a spatial scale of ~0.1 pc, which may correspond to the spatial scale at which gravity has overcome the thermal pressure. We suggest that G0.253+0.016 is on the verge of forming a cluster from hierarchical, filamentary structures that arise from a highly turbulent medium. Although the stellar distribution within Arches-like clusters is compact, centrally condensed and smooth, the observed gas distribution within G0.253+0.016 is extended, with no high-mass central concentration, and has a complex, hierarchical structure. If this clump gives rise to a high-mass cluster and its stars are formed from this initially hierarchical gas structure, then the resulting cluster must evolve into a centrally condensed structure via a dynamical process.

We recently proposed that the star-forming potential of dense molecular clouds in the Central Molecular Zone (CMZ, i.e. the central few 100 pc) of the Milky Way is linked to their orbital dynamics, potentially giving rise to an absolute-time sequence of star-forming clouds. In this paper, we present an orbital model for the gas stream(s) observed in the CMZ. The model is obtained by integrating orbits in the observed gravitational potential and represents a good fit to the distribution of dense gas, reproducing all of its key properties. The orbit is also consistent with observational constraints not included in the fitting process, such as the velocities of Sgr B2 and the Arches and Quintuplet clusters. It differs from previous models: (1) the orbit is open rather than closed due to the extended mass distribution in the CMZ, (2) its orbital velocity is twice as high as in previous models, and (3) Sgr A$^*$ coincides with the focus of the (eccentric) orbit rather than being offset. Our orbital solution supports the scenario in which the dust ridge between G0.253+0.016 ('the Brick') and Sgr B2 represents an absolute-time sequence of star-forming clouds, triggered by the tidal compression during their recent pericentre passage. We position the clouds on a common timeline and find that their pericentre passages occurred 0.30-0.74 Myr ago. Given their short free-fall times (0.3-0.4 Myr), the quiescent cloud G0.253+0.016 and the vigorously star-forming complex Sgr B2 are separated by a single free-fall time of evolution, implying that star formation proceeds rapidly once collapse has been initiated. We provide several quantitative predictions of our model and conclude with a discussion of the model in the Galactic context, highlighting its relation to large-scale gas accretion, the dynamics of the bar, the $x_2$ orbital family, and the origin of the Arches and Quintuplet clusters. (Abridged)

ALMA HCO+ observations of the infrared dark cloud G0.253+0.016 located in the Central Molecular Zone of the Galaxy are presented. The 89 GHz emission is area-filling, optically thick, and sub-thermally excited. Two types of filaments are seen in absorption against the HCO+ emission. Broad-line absorption filaments (BLAs) have widths of less than a few arcseconds (0.07 - 0.14 pc), lengths of 30 to 50 arcseconds (1.2 - 1.8 pc), and absorption profiles extending over a velocity range larger than 20 km/sec. The BLAs are nearly parallel to the nearby G0.18 non-thermal filaments and may trace HCO+ molecules gyrating about highly ordered magnetic fields located in front of G0.253+0.016 or edge-on sheets formed behind supersonic shocks propagating orthogonal to our line-of-sight in the foreground. Narrow-line absorption filaments (NLAs) have line-widths less than 20 km/sec. Some NLAs are also seen in absorption in other species with high optical depth such as HCN and occasionally in emission where the background is faint. The NLAs, which also trace low-density, sub-thermally excited HCO+ molecules, are mostly seen on the blueshifted side of the emission from G0.253+0.016. If associated with the surface of G0.253+0.016, the kinematics of the NLAs indicate that the cloud surface is expanding. The decompression of entrained, milli-Gauss magnetic fields may be responsible for the re-expansion of the surface layers of G0.253+0.016 as it recedes from the Galactic center following a close encounter with Sgr A.

Despite the simplicity of theoretical models of supersonically turbulent, isothermal media, their predictions successfully match the observed gas structure and star formation activity within low-pressure (P/k < 10^5 K cm^-3) molecular clouds in the solar neighbourhood. However, it is unknown if these theories extend to clouds in high-pressure (P/k > 10^7 K cm^-3) environments, like those in the Galaxy's inner 200 pc Central Molecular Zone (CMZ) and in the early Universe. Here we present ALMA 3mm dust continuum emission within a cloud, G0.253+0.016, which is immersed in the high-pressure environment of the CMZ. While the log-normal shape and dispersion of its column density PDF is strikingly similar to those of solar neighbourhood clouds, there is one important quantitative difference: its mean column density is 1--2 orders of magnitude higher. Both the similarity and difference in the PDF compared to those derived from solar neighbourhood clouds match predictions of turbulent cloud models given the high-pressure environment of the CMZ. The PDF shows a small deviation from log-normal at high column densities confirming the youth of G0.253+0.016. Its lack of star formation is consistent with the theoretically predicted, environmentally dependent volume density threshold for star formation which is orders of magnitude higher than that derived for solar neighbourhood clouds. Our results provide the first empirical evidence that the current theoretical understanding of molecular cloud structure derived from the solar neighbourhood also holds in high-pressure environments. We therefore suggest that these theories may be applicable to understand star formation in the early Universe.

Most stars form in a clustered environment. Therefore, it is important to assess how this environment influences the evolution of protoplanetary discs around young stars. In turn, this affects their ability to produce planets and ultimately life. We present here for the first time 3D SPH/N-body simulations that include both the hydrodynamical evolution of the discs around their natal stars, as well as the dynamics of the stars themselves. The discs are viscously evolving, accreting mass onto the central star and spreading. We find penetrating encounters to be very destructive for the discs as in previous studies, although the frequency of such encounters is low. We also find, however, that encounter influence the disc radii more strongly than other disc properties such as the disc mass. The disc sizes are set by the competition between viscous spreading and the disruptive effect of encounters. As discs spread, encounters become more and more important. In the regime of rapid spreading encounters simply truncate the discs, stripping the outer portions. In the opposite regime, we find that the effect of many distant encounters is able to limit the disc size. Finally, we predict from our simulations that disc sizes are limited by encounters at stellar densities exceeding $\sim 2-3 \times 10^3 \ \mathrm{pc}^{-2}$.

The star formation rate (SFR) in the Central Molecular Zone (CMZ, i.e. the central 500 pc) of the Milky Way is lower by a factor of >10 than expected for the substantial amount of dense gas it contains, which challenges current star formation theories. In this paper, we quantify which physical mechanisms could be responsible. On scales larger than the disc scale height, the low SFR is found to be consistent with episodic star formation due to secular instabilities or possibly variations of the gas inflow along the Galactic bar. The CMZ is marginally Toomre-stable when including gas and stars, but highly Toomre-stable when only accounting for the gas, indicating a low condensation rate of self-gravitating clouds. On small scales, we find that the SFR in the CMZ may be caused by an elevated critical density for star formation due to the high turbulent pressure. The existence of a universal density threshold for star formation is ruled out. The HI-H$_2$ phase transition of hydrogen, the tidal field, a possible underproduction of massive stars due to a bottom-heavy initial mass function, magnetic fields, and cosmic ray or radiation pressure feedback also cannot individually explain the low SFR. We propose a self-consistent cycle of star formation in the CMZ, in which the effects of several different processes combine to inhibit star formation. The rate-limiting factor is the slow evolution of the gas towards collapse - once star formation is initiated it proceeds at a normal rate. The ubiquity of star formation inhibitors suggests that a lowered central SFR should be a common phenomenon in other galaxies. We discuss the implications for galactic-scale star formation and supermassive black hole growth, and relate our results to the star formation conditions in other extreme environments.

Despite their importance as stellar nurseries and the building blocks of galaxies, very little is known about the formation of the highest mass clusters. The dense clump G0.253+0.016 represents an example of a clump that may form an Arches-like, high-mass cluster. Here we present molecular line maps toward G0.253+0.016 taken as part of the MALT90 molecular line survey, complemented with APEX observations. Combined, these data reveal the global physical properties and kinematics of G0.253+0.016. Recent Herschel data show that while the dust temperature is low (~19 K) toward its centre, the dust temperature on the exterior is higher (~27 K) due to external heating. Our new molecular line data reveal that, overall, the morphology of dense gas detected toward G0.253+0.016 matches very well its IR extinction and dust continuum emission. An anti-correlation between the dust and gas column densities toward its centre indicates that the clump is centrally condensed with a cold, dense interior in which the molecular gas is chemically depleted. The velocity field shows a strong gradient along the clump's major axis, with the blue-shifted side at higher Galactic longitude. The optically thick gas tracers are systematically red-shifted with respect to the optically thin and hot gas tracers, indicating radial motions. The gas kinematics and line ratios support the recently proposed scenario in which G0.253+0.016 results from a tidal compression during a recent pericentre passage near SgrA*. Because G0.253+0.016 represents an excellent example of a clump that may form a high-mass cluster, its detailed study should reveal a wealth of knowledge about the early stages of cluster formation.

Galactic scaling relations between the (surface densities of) the gas mass and the star formation (SF) rate are known to develop substantial scatter or even change form when considered below a certain spatial scale. We quantify how this behaviour should be expected due to the incomplete statistical sampling of independent star-forming regions. Other included limiting factors are the incomplete sampling of SF tracers from the stellar initial mass function and the spatial drift between gas and stars. We present a simple uncertainty principle for SF, which can be used to predict and interpret the failure of galactic SF relations on small spatial scales. This uncertainty principle explains how the scatter of SF relations depends on the spatial scale and predicts a scale-dependent bias of the gas depletion time-scale when centering an aperture on gas or SF tracer peaks. We show how the scatter and bias are sensitive to the physical size and time-scales involved in the SF process (such as its duration or the molecular cloud lifetime), and illustrate how our formalism provides a powerful tool to constrain these largely unknown quantities. Thanks to its general form, the uncertainty principle can also be applied to other astrophysical systems, e.g. addressing the time-evolution of star-forming cores, protoplanetary discs, or galaxies and their nuclei.

We review the formation and early evolution of the most massive and dense young stellar clusters, focusing on the role they can play in our understanding of star and planet formation as a whole. Young massive cluster (YMC) progenitor clouds in the Galactic Center can accumulate to a high enough density without forming stars that the initial protostellar densities are close to the final stellar density. For this to hold in the disk, the time scale to accumulate the gas to such high densities must be much shorter than the star formation timescale. Otherwise the gas begins forming stars while it is being accumulated to high density. The distinction between the formation regimes in the two environments is consistent with the predictions of environmentally-dependent density thresholds for star formation. This implies that stars in YMCs of similar total mass and radius can have formed at widely different initial protostellar densities. The fact that no systematic variations in fundamental properties are observed between YMCs in the disk and Galactic Center suggests stellar mass assembly is not strongly affected by the initial protostellar density. We review recent theoretical advances and summarize the debate on three key open questions: the initial (proto)stellar distribution, infant (im)mortality and age spreads within YMCs. We conclude: the initial protostellar distribution is likely hierarchical; YMCs likely experienced a formation history that was dominated by gas exhaustion rather than gas expulsion; YMCs are dynamically stable from a young age; and YMCs have age spreads much smaller than their mean age. Finally, we show that it is plausible that metal-rich globular clusters may have formed in a similar way to YMCs in nearby galaxies. In summary, the study of YMC formation bridges star/planet formation in the solar neighborhood to the oldest structures in the local Universe. [abridged]

At a distance of 8 kpc, the center of our Galaxy is the nearest galactic nucleus, and has been the subject of numerous key projects undertaken by great observatories such as Chandra, Spitzer, and Herschel. However, there are still no surveys of molecular gas properties in the Galactic center with less than 30" (1 pc) resolution. There is also no sensitive polarization survey of this region, despite numerous nonthermal magnetic features apparently unique to the central 300 parsecs. In this paper, we outline the potential the VLASS has to fill this gap. We assess multiple considerations in observing the Galactic center, and recommend a C-band survey with 10 micro-Jy continuum RMS and sensitive to molecular gas with densities greater than 10^4 cm^{-3}, covering 17 square degrees in both DnC and CnB configurations ( resolution ~5"), totaling 750 hours of observing time. Ultimately, we wish to note that the upgraded VLA is not just optimized for fast continuum surveys, but has a powerful correlator capable of simultaneously observing continuum emission and dozens of molecular and recombination lines. This is an enormous strength that should be fully exploited and highlighted by the VLASS, and which is ideally suited for surveying the center of our Galaxy.

Several recent studies have provided evidence for a `bottom-heavy' stellar initial mass function (IMF) in massive elliptical galaxies. Here we investigate the influence of the IMF shape on the recently discovered color-magnitude relation (CMR) among globular clusters (GCs) in such galaxies. To this end we use calculations of GC mass loss due to stellar and dynamical evolution to evaluate (i) the shapes of stellar mass functions in GCs after 12 Gyr of evolution as a function of current GC mass along with their effects on integrated-light colors and mass-to-light ratios, and (ii) their impact on the effects of GC self-enrichment using the 2009 `reference' model of Bailin & Harris. As to the class of metal-poor GCs, we find the observed shape of the CMR (often referred to as the `blue tilt') to be very well reproduced by Bailin & Harris' reference self-enrichment model once 12 Gyr of GC mass loss is taken into account. The influence of the IMF on this result is found to be insignificant. However, we find that the observed CMR among the class of metal-rich GCs (the `red tilt') can only be adequately reproduced if the IMF was bottom-heavy (e.g., $-3.0 <~ \alpha <~ -2.3$ in $dN/dM \propto M^{\alpha}$), which causes the stellar mass function at subsolar masses to depend relatively strongly on GC mass. This constitutes additional evidence that the metal-rich stellar populations in giant elliptical galaxies were formed with a bottom-heavy IMF.

Detailed observations of the nearest star-forming regions in the Milky Way (MW) provide the ultimate benchmark for studying star formation. The extent to which the results of these Galaxy-based studies can be extrapolated to extragalactic systems depends on the overlap of the environmental conditions probed. In this paper, we compare the properties of clouds and star-forming regions in the MW with those in nearby galaxies and in the high-redshift Universe. We find that in terms of their baryonic composition, kinematics, and densities, the clouds in the solar neighbourhood are similar to those in nearby galaxies. The clouds and regions in the Central Molecular Zone (CMZ, i.e. the inner 250 pc) of the MW are indistinguishable from high-redshift clouds and galaxies. The presently low star formation rate in the CMZ therefore implies that either (1) its gas represents the initial conditions for high-redshift starbursts or (2) some yet unidentified process consistently suppresses star formation over >100 Myr time-scales. We conclude that the MW contains large reservoirs of gas with properties directly comparable to most of the known range of star formation environments and is therefore an excellent template for studying star formation across cosmological time-scales.

Intermediate-mass black holes (IMBHs) may provide the missing link to understanding the growth of supermassive black holes in the early Universe. Some formation scenarios predict that IMBHs could have formed by runaway collisions in globular clusters (GCs). Understanding the influence of an IMBH in the center of a GC on its environment might provide indirect detection methods. Our goal is to test the effects of different initial compositions of GCs on their evolution in a tidal field. We pin down the crucial observables that indicate the presence of an IMBH at the center of the cluster. In addition to central IMBHs, we also consider the effects of different stellar-mass black hole retention and stellar binary fractions. We performed a set of 22 N-body simulations and varied particle numbers, IMBH masses, stellar-mass black-hole retention fractions, and stellar binary fractions. These models are all run in an external tidal field to study the effect of black holes on the cluster mass loss, mass function, and life times. Finally, we compared our results with observational data. We found that a central massive black hole increases the escape rate of high-mass stars from a star cluster, implying that the relative depletion of the mass function at the low-mass end proceeds less rapidly. Furthermore, we found a similar behavior for a cluster hosting a high number of stellar-mass black holes instead of one massive central IMBH. The presence of an IMBH also weakly affects the fraction of the cluster mass that is constituted by stellar remnants, as does the presence of primordial binaries. We compared our simulations with observational data from the literature and found good agreement between our models and observed mass functions and structural parameters of GCs. We exploited this agreement to identify GCs that could potentially host IMBHs.

The masses of supermassive black holes (SMBHs, M_BH=10^6-10^11 Msun) in the centres of galaxies are related to the host stellar spheroid mass and velocity dispersion. A key question is how these relations originate, and over which range of black hole masses they hold. It has been speculated that intermediate-mass black holes (IMBHs, M_BH=10^2-10^5 Msun) could play a fundamental role in the growth of SMBHs. A handful of IMBHs has recently been detected in Galactic globular clusters (GCs), but their masses are inconsistent with the galactic scaling relations of SMBHs. In this Letter, we derive the initial properties of the GCs using a standard analytical evolutionary model, of which the free parameters are fixed by independent constraints. We find that the observed IMBH masses initially followed the galactic SMBH scaling relations, and subsequently moved off these relations due to the dynamical evolution of their host GCs. This work is concluded with a brief discussion of the uncertainties and the implications of our results for the possible universality of massive black hole growth.

Numerical simulations of star cluster formation have advanced greatly during the past decade, covering increasingly massive gas clouds while accounting for more and more complex physics. In this review, I discuss the present state of the field, paying particular attention to the key physics that need to be included in cluster formation simulations. The main numerical techniques are summarized for a broad audience, before evaluating their application to the problem of cluster formation. A faithful reproduction of the observed characteristics of cluster formation can presently be achieved in numerical simulations. Ideally, this requires turbulent initial conditions to be combined with radiative feedback, protostellar outflows, and magnetic fields. With the exciting prospect in mind that our understanding of cluster formation will soon be revolutionized by facilities like ALMA, JWST, and the EVLA, this review also identifies a number of areas that would particularly benefit from a joint observational and theoretical effort.

Super star clusters are the end product of star formation under the most extreme conditions. As such, studying how their final stellar populations are assembled from their natal progenitor gas clouds can provide strong constraints on star formation theories. An obvious place to look for the initial conditions of such extreme stellar clusters are gas clouds of comparable mass and density, with no star formation activity. We present a method to identify such progenitor gas clouds and demonstrate the technique for the gas in the inner few hundred pc of our Galaxy. The method highlights three clouds in the region with similar global physical properties to the previously identified extreme cloud, G0.253+0.016, as potential young massive cluster (YMC) precursors. The fact that four potential YMC progenitor clouds have been identified in the inner 100 pc of the Galaxy, but no clouds with similar properties have been found in the whole first quadrant despite extensive observational efforts, has implications for cluster formation/destruction rates across the Galaxy. We put forward a scenario to explain how such dense gas clouds can arise in the Galactic centre environment, in which YMC formation is triggered by gas streams passing close to the minimum of the global Galactic gravitational potential at the location of the central supermassive black hole, Sgr A*. If this triggering mechanism can be verified, we can use the known time interval since closest approach to Sgr A* to study the physics of stellar mass assembly in an extreme environment as a function of absolute time.

We report a systematic and statistically significant offset between the optical (g-z or B-I) colors of seven massive elliptical galaxies and the mean colors of their associated massive metal-rich globular clusters (GCs) in the sense that the parent galaxies are redder by 0.12-0.20 mag at a given galactocentric distance. However, spectroscopic indices in the blue indicate that the luminosity-weighted ages and metallicities of such galaxies are equal to that of their averaged massive metal-rich GCs at a given galactocentric distance, to within small uncertainties. The observed color differences between the red GC systems and their parent galaxies cannot be explained by the presence of multiple stellar generations in massive metal-rich GCs, as the impact of the latter to the populations' integrated g-z or B-I colors is found to be negligible. However, we show that this paradox can be explained if the stellar initial mass function (IMF) in these massive elliptical galaxies was significantly steeper at subsolar masses than canonical IMFs derived from star counts in the solar neighborhood, with the GC colors having become bluer due to dynamical evolution, causing a significant flattening of the stellar MF of the average surviving GC.

We present observational evidence of environmental effects on the formation and evolution of planetary systems. We combine catalogues of resolved protoplanetary discs (PPDs) and young stellar objects in the solar neighbourhood to analyse the PPD size distribution as a function of ambient stellar density. By running Kolmogorov-Smirnov tests between the PPD radii at different densities, we find empirical evidence, at the >97% confidence level, for a change in the PPD radius distribution at ambient stellar densities \Sigma > 10^{3.5} pc^{-2}. This coincides with a simple theoretical estimate for the truncation of PPDs or planetary systems by dynamical encounters. If this agreement is causal, the ongoing disruption of PPDs and planetary systems limits the possible existence of planets in the habitable zone, with shorter lifetimes at higher host stellar masses and ambient densities. Therefore, habitable planets are not likely to be present in long-lived stellar clusters, and may have been ejected altogether to form a population of unbound, free-floating planets. We conclude that, while highly suggestive, our results should be verified through other methods. Our simple model shows that truncations should lead to a measurable depletion of the PPD mass function that can be detected with ALMA observations of the densest nearby and young clusters.

We recently presented a model for the cluster formation efficiency (CFE), i.e. the fraction of star formation occurring in bound stellar clusters. It utilizes the idea that the formation of stars and stellar clusters occurs across a continuous spectrum of ISM densities. Bound stellar clusters naturally arise from the high-density end of this density spectrum. Due to short free-fall times, these high-density regions can achieve high star formation efficiencies (SFEs) and can be unaffected by gas expulsion. Lower-density regions remain gas-rich and substructured, and are unbound upon gas expulsion. The model enables the CFE to be calculated using galactic-scale observables. I present a brief summary of the model physics, assumptions and caveats, and show that it agrees well with observations. Fortran and IDL routines for calculating the CFE are publicly available at http://www.mpa-garching.mpg.de/cfe.

We present a theoretical framework in which bound stellar clusters arise naturally at the high-density end of the hierarchy of the interstellar medium (ISM). Due to short free-fall times, these high-density regions achieve high local star formation efficiencies, enabling them to form bound clusters. Star-forming regions of lower density remain substructured and gas-rich, ending up unbound when the residual gas is expelled. Additionally, the tidal perturbation of star-forming regions by nearby, dense giant molecular clouds imposes a minimum density contrast required for the collapse to a bound cluster. The fraction of all star formation that occurs in bound stellar clusters (the cluster formation efficiency or CFE) follows by integration of these local clustering and survival properties over the full density spectrum of the ISM, and hence is set by galaxy-scale physics. We derive the CFE as a function of observable galaxy properties, and find that it increases with the gas surface density, from ~1% in low-density galaxies to a peak value of ~70% at densities of ~10^3 Msun pc^-2. This explains the observation that the CFE increases with the star formation rate density in nearby dwarf, spiral, and starburst galaxies. Indeed, comparing our model results with observed galaxies yields excellent agreement. The model is applied further by calculating the spatial variation of the CFE within single galaxies. We also consider the variation of the CFE with cosmic time and show that it increases with redshift, peaking in high-redshift, gas-rich disc galaxies. It is estimated that up to 30-35% of all stars in the Universe once formed in bound stellar clusters. We discuss how our theory can be verified with Gaia and ALMA, and provide implementations for future theoretical work and for simulations of galaxy formation and evolution.

Until the early 2000s, the research portfolio of the Astronomical Institute in Utrecht (SIU) did not include galaxy evolution. Somewhat serendipitously, this changed with the advent of the star cluster group. In only a few years, a simple framework was developed to describe and quantify the properties of dynamically evolving star cluster populations. Since then, the `Utrecht cluster disruption model' has shown that the galactic environment plays an important role in setting the evolution of stellar clusters. From this simple result, it follows that cluster populations bear some imprint of the characteristics and histories of their host galaxies, and that star clusters can be used to trace galaxy evolution -- an aim for which the Utrecht star cluster models were never designed, but which they are well-capable of fulfilling. I review some of the work in this direction, with a strong emphasis on the contributions from the SIU.

We report the detection of five high-velocity stars in the core of the globular cluster NGC 2808. The stars lie on the the red giant branch and show total velocities between 40 and 45 km/s. For a core velocity dispersion sigma_c = 13.4 km/s, this corresponds to up to 3.4 sigma_c. These velocities are close to the estimated escape velocity (~ 50 km/s) and suggest an ejection from the core. Two of these stars have been confirmed in our recent integral field spectroscopy data and we will discuss them in more detail here. These two red giants are located at a projected distance of ~ 0.3 pc from the center. According to their positions on the color magnitude diagram, both stars are cluster members. We investigate several possible origins for the high velocities of the stars and conceivable ejection mechanisms. Since the velocities are close to the escape velocity, it is not obvious whether the stars are bound or unbound to the cluster. We therefore consider both cases in our analysis. We perform numerical simulations of three-body dynamical encounters between binaries and single stars and compare the resulting velocity distributions of escapers with the velocities of our stars. We compare the predictions for a single dynamical encounter with a compact object with those of a sequence of two-body encounters due to relaxation. If the stars are unbound, the encounter must have taken place recently, when the stars were already in the giant phase. After including binary fractions and black-hole retention fractions, projection effects, and detection probabilities from Monte-Carlo simulations, we estimate the expected numbers of detections for all the different scenarios. Based on these numbers, we conclude that the most likely scenario is that the stars are bound and were accelerated by a single encounter between a binary of main-sequence stars and a ~ 10 M_sun black hole.

Whether massive stars can occasionally form in relative isolation or if they require a large cluster of lower-mass stars around them is a key test in the differentiation of star formation theories as well as how the initial mass function of stars is sampled. Previous attempts to find O-type stars that formed in isolation were hindered by the possibility that such stars are merely runaways from clusters, i.e., their current isolation does not reflect their birth conditions. We introduce a new method to find O-type stars that are not affected by such a degeneracy. Using the VLT-FLAMES Tarantula Survey and additional high resolution imaging we have identified stars that satisfy the following constraints: 1) they are O-type stars that are not detected to be part of a binary system based on RV time series analysis; 2) they are designated spectral type O7 or earlier ; 3) their velocities are within 1\sigma of the mean of OB-type stars in the 30 Doradus region, i.e. they are not runaways along our line-of-sight; 4) the projected surface density of stars does not increase within 3 pc towards the O-star (no evidence for clusters); 5) their sight lines are associated with gaseous and/or dusty filaments in the ISM, and 6) if a second candidate is found in the direction of the same filament with which the target is associated, both are required to have similar velocities. With these criteria, we have identified 15 stars in the 30 Doradus region, which are strong candidates for being high-mass stars that have formed in isolation. Additionally, we employed extensive MC stellar cluster simulations to confirm that our results rule out the presence of clusters around the candidates. Eleven of these are classified as Vz stars, possibly associated with the zero-age main sequence. We include a newly discovered W-R star as a candidate, although it does not meet all of the above criteria.

(Abridged) Interacting galaxies are well-known for their high star formation rates and rich star cluster populations, but the rapidly changing tidal field can also efficiently destroy clusters. We use numerical simulations of merging disc galaxies to investigate which mechanism dominates. The simulations include a model for the formation and dynamical disruption of the entire star cluster population. We find that the dynamical heating of clusters by tidal shocks is about an order of magnitude higher in interacting galaxies than in isolated galaxies. This is driven by the increased gas density, and is sufficient to destroy star clusters at a higher rate than new clusters are formed: the total number of clusters in the merger remnant is 2-50% of the amount in the progenitor discs, with low-mass clusters being disrupted preferentially. By adopting observationally motivated selection criteria, we find that the observed surplus of star clusters in nearby merging galaxies is caused by the bias to detect young, massive clusters. We provide a general expression for the survival fraction of clusters, which increases with the gas depletion time-scale. Due to the preferential disruption of low-mass clusters, the mass distribution of the surviving star clusters in a merger remnant develops a peak at a mass of about 10^3 Msun, which evolves to higher masses at a rate of 0.3-0.4 dex per Gyr. The peak mass initially depends weakly on the galactocentric radius, but this correlation disappears as the system ages. We discuss the similarities between the cluster populations of the simulated merger remnants and (young) globular cluster systems. Our results suggest that the combination of cluster formation and destruction should be widespread in the dense star-forming environments at high redshifts, which could provide a natural origin to present-day globular cluster systems.

We test whether the masses of old globular clusters (GCs) in dwarf galaxies are consistent with the same initial mass spectrum as young massive clusters (YMCs) in nearby star-forming galaxies. The most massive GCs of dwarf galaxies are compared to their expected masses when drawing from the Schechter-type ICMF of YMCs. It is found that the most massive GCs of galaxies in the stellar mass range M_*,gal=10^7-10^9 Msun are consistent with the same initial mass spectrum as YMCs in about 90% of the cases, suggesting that their formation mechanisms were the same. For the remaining 10%, the most massive clusters are nuclear GCs, which have been able to grow to higher masses through further merging after their initial formation ended. Because the effects of cluster disruption are weaker for more massive clusters, we estimate that up to one third of the metal-poor GCs in the Milky Way may have a nuclear origin, while the remaining two thirds formed through the same process as YMCs in the local universe. A log-normal ICMF is inconsistent with observed GCs at a 99.6% confidence level.

The fraction of star formation that results in bound star clusters is influenced by the density spectrum in which stars are formed and by the response of the stellar structure to gas expulsion. We analyse hydrodynamical simulations of turbulent fragmentation in star-forming regions to assess the dynamical properties of the resulting population of stars and (sub)clusters. Stellar subclusters are identified using a minimum spanning tree algorithm. When considering only the gravitational potential of the stars and ignoring the gas, we find that the identified subclusters are close to virial equilibrium (the typical virial ratio Q_vir~0.59, where virial equilibrium would be Q_vir~0.5). This virial state is a consequence of the low gas fractions within the subclusters, caused by the accretion of gas onto the stars and the accretion-induced shrinkage of the subclusters. Because the subclusters are gas-poor, up to a length scale of 0.1-0.2 pc at the end of the simulation, they are only weakly affected by gas expulsion. The fraction of subclusters that reaches the high density required to evolve to a gas-poor state increases with the density of the star-forming region. We extend this argument to star cluster scales, and suggest that the absence of gas indicates that the early disruption of star clusters due to gas expulsion (infant mortality) plays a smaller role than anticipated, and is potentially restricted to star-forming regions with low ambient gas densities. We propose that in dense star-forming regions, the tidal shocking of young star clusters by the surrounding gas clouds could be responsible for the early disruption. This `cruel cradle effect' would work in addition to disruption by gas expulsion. We suggest possible methods to quantify the relative contributions of both mechanisms.

The fraction of star formation that results in bound stellar clusters (cluster formation efficiency or CFE) is a central quantity in many studies of star formation, star clusters and galaxies. Recent results suggest that contrary to popular assumption, the CFE is not (solely) set by gas expulsion, but is also influenced by the primordial environment, although its precise behaviour remains unknown. Here it is discussed which mechanisms set the CFE, which recent advancements have been made to disentangle their contributions, and which studies are needed in the near future to achieve a quantitative understanding of the CFE.

(Abridged) The formation and evolution of star cluster populations are related to the galactic environment. Cluster formation is governed by processes acting on galactic scales, and star cluster disruption is driven by the tidal field. In this paper, we present a self-consistent model for the formation and evolution of star cluster populations, for which we combine an N-body/SPH galaxy evolution code with semi-analytic models for star cluster evolution. The model includes star formation, feedback, stellar evolution, and star cluster disruption by two-body relaxation and tidal shocks. We apply the model by simulating a suite of 9 isolated disc galaxies and 24 galaxy mergers. The evolutionary histories of individual clusters in these simulations are discussed to illustrate how the environment of clusters changes in time and space. The resulting variability of the disruption rate with time and space affects the properties of star cluster populations. The combined effect of clusters escaping their dense formation sites (`cluster migration') and the preferential disruption of clusters residing in dense environments (`natural selection') implies that the mean disruption rate of the population decreases with cluster age. This affects the slope of the cluster age distribution, which becomes a function of the star formation rate density. The evolutionary histories of clusters in a galaxy merger vary widely and determine which clusters survive the merger. This impacts the age distributions and the locations of the surviving clusters at all times during a merger. We conclude that accounting for the interplay between the formation, disruption, and orbital histories of clusters enables a more sophisticated interpretation of observed properties of cluster populations, thereby extending the role of cluster populations as tracers of galaxy evolution.

The relation between the globular cluster luminosity function (GCLF, dN/dlogL) and globular cluster mass function (GCMF, dN/dlogM) is considered. Due to low-mass star depletion, dissolving GCs have mass-to-light (M/L) ratios that are lower than expected from their metallicities. This has been shown to lead to an M/L ratio that increases with GC mass and luminosity. We model the GCLF and GCMF and show that the power law slopes inherently differ (1.0 versus 0.7, respectively) when accounting for the variability of M/L. The observed GCLF is found to be consistent with a Schechter-type initial cluster mass function and a mass-dependent mass-loss rate.

The observed mass-to-light (M/L) ratios of globular clusters (GCs) are on average ~20% lower than expected from Simple Stellar Population (SSP) models, which only account for the effects of stellar evolution. We study the M/L ratio evolution of a sample of 24 Galactic GCs using parameterised cluster models. The dynamical evolution of GCs is included by accounting for their dissolution and by using a detailed description of the evolution of the stellar mass function. The ejection of low-mass stars leads to a decrease of M/L, which is found to explain the discrepancy between the observations and SSP models.

(Abridged) The dynamical ejection of stars from star clusters affects the shape of the stellar mass function (MF) in these clusters, because the escape probability of a star depends on its mass. The objective of this paper is to provide and to apply a simple physical model for the evolution of the MF in star clusters for a large range of the parameter space. The model is derived from the basic principles of two-body encounters and energy considerations. It is independent of the adopted mass loss rate or initial mass function (IMF), and contains stellar evolution, stellar remnant retention, dynamical dissolution in a tidal field, and mass segregation. It is found that the MF evolution in star clusters depends on the disruption time, remnant retention fraction, initial-final stellar mass relation, and IMF. Low-mass stars are preferentially ejected after t~400 Myr. Before that time, masses around 15-20% of the maximum stellar mass are lost. The degree of low-mass star depletion grows for increasing disruption times, but can be quenched when the retained fraction of massive remnants is large. The highly depleted MFs of certain Galactic globular clusters are explained by the enhanced low-mass star depletion that occurs for low remnant retention fractions. Unless the retention fraction is exceptionally large, dynamical evolution always decreases the mass-to-light ratio. The modeled evolution of the MF is consistent with N-body simulations when adopting identical boundary conditions. However, it is found that the results from N-body simulations only hold for their specific boundary conditions and should not be generalised to all clusters. It is concluded that the model provides an efficient method do understand the evolution of the stellar MF in star clusters under widely varying conditions.

The conversion of the globular cluster luminosity function (GCLF, dN/dlogL) to the globular cluster mass function (GCMF, dN/dlogM) is addressed. Dissolving globular clusters (GCs) become preferentially depleted in low-mass stars, which have a high mass-to-light ratio. This has been shown to result in a mass-to-light ratio (M/L) that increases with GC luminosity or mass, because more massive GCs have lost a smaller fraction of their stars than low-mass GCs. Using GC models, we study the influence of the luminosity dependency of M/L on the inferred GCMF. The observed GCLF is consistent with a powerlaw or Schechter type GC initial mass function in combination with a cluster mass-dependent mass loss rate. Below the peak, the logarithmic slope of the GCMF is shallower than that of the GCLF (0.7 versus 1.0), whereas the peak mass is 0.1-0.3 dex lower when accounting for the variability of M/L than in the case where a constant M/L is adopted.

Differences between the inferred star formation histories (SFHs) of star clusters and field stars seem to suggest distinct star formation processes for the two. The Small Magellanic Cloud (SMC) is an example of a galaxy where such a discrepancy is observed. We model the observed age distributions of the SMC clusters and field stars using a new population synthesis code, SPACE, that includes stellar evolution, infant mortality and cluster dissolution. We find that the two observed age distributions can be explained by a single SFH, thus eliminating the need to assume two separate mechanisms for star formation.

The observed dynamical mass-to-light (M/L) ratios of globular clusters (GCs) are systematically lower than those expected from `canonical' simple stellar population models, which do not account for the preferential loss of low-mass stars due to energy equipartition. It was recently shown that low-mass star depletion can qualitatively explain the M/L discrepancy. To verify whether it is indeed the driving mechanism, we derive dissolution timescales and use these to predict the M/L_V ratios of the 24 Galactic GCs for which orbital parameters and dynamical M/L_V are known. We also predict the slopes of their low-mass stellar mass functions (MFs). We use the SPACE cluster models, which include dynamical dissolution, low-mass star depletion, stellar evolution, stellar remnants and various metallicities. The predicted M/L_V are in 1 sigma agreement with the observations for 12 out of 24 GCs. The discrepancy for the other GCs probably arises because our predictions give global M/L ratios, while the observations represent extrapolated central values that are different from global ones in case of mass segregation and a long dissolution timescale. GCs in our sample which likely have dissimilar global and central M/L ratios can be excluded by imposing limits on the dissolution timescale and King parameter. For the remaining GCs, the observed and predicted average M/L_V are 78^+9_-11% and 78+/-2% of the canonically expected values, while for the entire sample the values are 74^+6_-7% and 85+/-1%. The predicted correlation between the slope of the low-mass stellar MF and M/L_V is qualitatively consistent with observed MF slopes. It is concluded that the variation of M/L ratio due to dissolution and low-mass star depletion is a plausible explanation for the discrepancy between the observed and canonically expected M/L ratios of GCs. (Abridged)

To obtain an accurate description of broad-band photometric star cluster evolution, certain effects should be accounted for, such as the preferential loss of low-mass stars and the retain of stellar remnants. Moreover, the IMF and metallicity affect photometry as well. Due to their advanced evolutionary state, globular clusters (GCs) are interesting test cases for cluster models describing these effects. In this paper we describe cluster models in which the above effects are included. The photometric evolution of clusters is predicted, and the results are applied to Galactic GCs. The presented cluster models represent an analytical description of the evolution of the underlying stellar mass function due to stellar evolution and dynamical cluster dissolution. Stellar remnants are included by using initial-remnant mass relations, while cluster photometry is computed from the Padova 1999 isochrones. Our study shows that the preferential loss of low-mass stars, the retain of stellar remnants, as well as IMF and metallicity variations strongly influence the magnitude, colour and mass-to-light ratio evolution of clusters. The different effects can be clearly separated with our models. The models are applied to the Galactic GC population, of which the magnitude, colour and mass-to-light ratio ranges are well reproduced. We also show that the mass-to-light ratios of clusters of similar ages and metallicities cannot be assumed to be constant for all cluster luminosities. Instead, mass-to-light ratio increases with cluster luminosity and mass. These models underline the importance of detailed cluster models when considering cluster photometry. With the analytic framework provided in this paper, observed cluster properties can be interpreted in a more complete perspective. (Abridged)

The majority of observed mass-to-light ratios of globular clusters are too low to be explained by `canonical' cluster models, in which dynamical effects are not accounted for. Moreover, these models do not reproduce a recently reported trend of increasing M/L with cluster mass, but instead predict mass-to-light ratios that are independent of cluster mass for a fixed age and metallicity. This study aims to explain the M/L of globular clusters in four galaxies by including stellar evolution, stellar remnants, and the preferential loss of low-mass stars due to energy equipartition. Analytical cluster models are applied that account for stellar evolution and dynamical cluster dissolution to samples of globular clusters in Cen A, the Milky Way, M31 and the LMC. The models include stellar remnants and cover metallicities in the range Z=0.0004-0.05. Both the low observed mass-to-light ratios and the trend of increasing M/L with cluster mass can be reproduced by including the preferential loss of low-mass stars, assuming reasonable values for the dissolution timescale. This leads to a mass-dependent M/L evolution and increases the explained percentage of the observations from 39% to 92%. This study shows that the hitherto unexplained discrepancy between observations and models of the mass-to-light ratios of globular clusters can be explained by dynamical effects, provided that the globular clusters exhibiting low M/L have dissolution timescales within the ranges assumed in this Letter. Furthermore, it substantiates that M/L cannot be assumed to be constant with mass at fixed age and metallicity.