Source author record

J. M. Diederik Kruijssen

J. M. Diederik Kruijssen appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.GA astro-ph.CO astro-ph.SR astro-ph.EP astro-ph.IM astro-ph.HE physics.data-an physics.ins-det

Catalog footprint

What is connected

100works

8topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2026arXiv

Unveiling the 3D structure of the central molecular zone from stellar kinematics and photometry: The 50 and 20 km/s clouds

The central molecular zone (CMZ), surrounding the Galactic centre, is the largest reservoir of dense molecular gas in the Galaxy. Despite its relative proximity, the 3D structure of the CMZ remains poorly constrained, primarily due to projection effects. We aim to constrain the line-of-sight location of two molecular clouds in the CMZ -- the 50 and 20 km/s clouds -- and to investigate their possible physical connection using stellar kinematics and photometry. This study serves as a pilot for future applications across the full CMZ. We estimated the line-of-sight position of the clouds by analysing stellar kinematics, stellar densities, and stellar populations towards the cloud regions and a control field. We find an absence of westward moving stars in the cloud regions, which indicates that they lie on the near side of the CMZ. This interpretation is supported by the stellar density distributions. The similar behaviour observed in the two clouds, as well as in the region between them (the ridge), suggests that they are located at comparable distances and are physically linked. We also identified an intermediate-age stellar population (2-7 Gyr) in both regions, consistent with that observed on the near side of the CMZ. We estimated the line-of-sight distances at which the clouds and the ridge become kinematically detectable (i.e. where the proper motion component parallel to the Galactic plane differs from that of the control field at the 3 sigma level) by converting their measured proper motions parallel to the Galactic plane using a theoretical model of the stellar distribution. We find that the 50 and 20 km/s clouds are located at $43\pm8$ pc and $56\pm11$ pc from Sgr A*, respectively, and that the ridge lies at $56\pm11$ pc; this supports the idea that the clouds are physically connected through the ridge.

preprint2025arXiv

Rotation and stability of the circumnuclear gas disk in the Galactic Center potential by the ALMA CMZ Exploration Survey (ACES)

We investigated the gravitational potential and mass distribution in the Galactic Center by examining the morphology and kinematics of the circumnuclear gaseous disk revealed by the molecular line data from the ALMA CMZ Exploration Survey (ACES). We obtain an estimate of the shape of the potential {within the central $\sim 20$ pc} to reproduce the observed properties of the circumnuclear gas disk (CND) by simulating the motion of test particles for various axial ratios and show that the potential is approximately spherical. We construct a rotation curve by applying the terminal velocity method to the position-velocity diagrams, and calculate the mass distribution in the Galactic Center. The distribution of mass density is found to be of cusp type, approximated by $ρ_{\rm mass} \sim 1.56\times 10^5(R/1 {\rm pc})^{-1.9}~M_{\odot} {\rm pc}^{-3}$, where $R$ is the distance from the nucleus. We discuss the tidal effect caused by the gravitational potential that produces the rotation curve and show that the gas disk is stable against self-gravitational contraction within a critical radius of $ R_{\rm T}\sim 14 ~(ρ_{\rm gas}/10^5 {\rm H_2~cm^{-3}})^{-1/2}~{\rm pc}$. This suggests suppression of star formation and a top-heavy IMF in the circmunuclear region.

preprint2023arXiv

In-situ or accreted? Using deep learning to infer the origin of extragalactic globular clusters from observables

Globular clusters (GCs) are powerful tracers of the galaxy assembly process, and have already been used to obtain a detailed picture of the progenitors of the Milky Way. Using the E-MOSAICS cosmological simulation of a (34.4 Mpc)$^3$ volume that follows the formation and co-evolution of galaxies and their star cluster populations, we develop a method to link the origin of GCs to their observable properties. We capture this complex link using a supervised deep learning algorithm trained on the simulations, and predict the origin of individual GCs (whether they formed in the main progenitor or were accreted from satellites) based solely on extragalactic observables. An artificial neural network classifier trained on $\sim50,000$ GCs hosted by $\sim 700$ simulated galaxies successfully predicts the origin of GCs in the test set with a mean accuracy of $89$ per cent for the objects with [Fe/H]<-0.5 that have unambiguous classifications. The network relies mostly on the alpha-element abundances, metallicities, projected positions, and projected angular momenta of the clusters to predict their origin. A real-world test using the known progenitor associations of the Milky Way GCs achieves up to $90$ per cent accuracy, and successfully identifies as accreted most of the GCs in the inner Galaxy associated to the Kraken progenitor, as well as all the Gaia-Enceladus GCs. We demonstrate that the model is robust to observational uncertainties, and develop a method to predict the classification accuracy across observed galaxies. The classifier can be optimized for available observables (e.g. to improve the accuracy by including GC ages), making it a valuable tool to reconstruct the assembly histories of galaxies in upcoming wide-field surveys.

preprint2022arXiv

A 2-3 mm high-resolution molecular line survey towards the centre of the nearby spiral galaxy NGC 6946

The complex physical, kinematic, and chemical properties of galaxy centres make them interesting environments to examine with molecular line emission. We present new $2-4$" (${\sim}75{-}150$ pc at $7.7$ Mpc) observations at 2 and 3 mm covering the central $50$" (${\sim}1.9$ kpc) of the nearby double-barred spiral galaxy NGC 6946 obtained with the IRAM Plateau de Bure Interferometer. We detect spectral lines from ten molecules: CO, HCN, HCO$^+$, HNC, CS, HC$_3$N, N$_2$H$^+$, C$_2$H, CH$_3$OH, and H$_2$CO. We complemented these with published 1mm CO observations and 33 GHz continuum observations to explore the star formation rate surface density ${Σ_{\mathrm{SFR}}}$ on 150 pc scales. In this paper, we analyse regions associated with the inner bar of NGC 6946 $-$ the nuclear region (NUC), the northern (NBE), and southern inner bar end (SBE) and we focus on short-spacing corrected bulk (CO) and dense gas tracers (HCN, HCO$^+$, and HNC). We find that HCO$^+$ correlates best with ${Σ_{\mathrm{SFR}}}$, but the dense gas fraction ($f_{\mathrm{dense}}$) and star formation efficiency of the dense gas (${\mathrm{SFE_{dense}}}$) fits show different behaviours than expected from large-scale disc observations.The SBE has a higher ${Σ_{\mathrm{SFR}}}$, $f_{\mathrm{dense}}$, and shocked gas fraction than the NBE. We examine line ratio diagnostics and find a higher CO(2-1)/CO(1-0) ratio towards NBE than for the NUC. Moreover, comparison with existing extragalactic datasets suggests that using the HCN/HNC ratio to probe kinetic temperatures is not suitable on kiloparsec and sub-kiloparsec scales in extragalactic regions. Lastly, our study shows that the HCO$^+$/HCN ratio might not be a unique indicator to diagnose AGN activity in galaxies.

preprint2022arXiv

A CO isotopologue Line Atlas within the Whirlpool galaxy Survey (CLAWS)

We present the CO isotopologue Line Atlas within the Whirpool galaxy Survey (CLAWS) based on an IRAM 30-m large programme which provides a benchmark study of numerous, faint CO isotopologues in the mm-wavelength regime across the full disc of M51 (NGC 5194). The survey's core goal is to use the low-J CO isotopologue lines to constrain CO excitation and chemistry, and therefrom the local physical conditions of the gas. In this survey paper, we describe the CLAWS observing and data reduction strategies. We map the J=1-0 and 2-1 transitions of the CO isotopologues $^{12}$CO,$^{13}$CO, C$^{18}$O and C$^{17}$O, as well as several supplementary lines within the 1 mm and 3 mm window (CN(1-0), CS(2-1), CH$_3$OH(2-1), N$_2$H$^+$(1-0), HC$_3$N(10-9)) at ~1 kpc resolution. A total observation time of 149 h offers unprecedented sensitivity. We use these data to explore several CO isotopologue line ratios in detail, study their radial (and azimuthal) trends and investigate whether changes in line ratios stem from changes in ISM properties such as gas temperatures, densities or chemical abundances. For example, we find negative radial trends for the $^{13}$CO}/$^{12}$CO, C$^{18}$O/$^{12}$CO and C$^{18}$O/$^{13}$CO line ratios in their J=1-0 transitions. We also find variations with local environment, such as higher $^{12}$CO(2-1)/(1-0) or $^{13}$CO/$^{12}$CO(1-0) line ratios in interarm regions compared to spiral arm regions. We propose that these aforementioned variations of CO line ratios are most likely due to a variation of the optical depth, while abundance variations due to selective nucleosynthesis on a galaxy-wide scale could also play a role. We also study the CO spectral line energy distribution (SLED) using archival JCMT $^{12}$CO(3-2) data and find a variation of the SLED shape with local environmental parameters further underlying changes in optical depth, gas temperatures or densities.

preprint2022arXiv

A tale of two DIGs: The relative role of HII regions and low-mass hot evolved stars in powering the diffuse ionised gas (DIG) in PHANGS-MUSE galaxies

We use integral field spectroscopy from the PHANGS-MUSE survey, which resolves the ionised interstellar medium at ${\sim}50$ pc resolution in 19 nearby spiral galaxies, to study the origin of the diffuse ionised gas (DIG). We examine the physical conditions of the diffuse gas by first removing morphologically-defined HII regions and then binning the low-surface-brightness areas to achieve significant detections of the key nebular lines. A simple model for the leakage and propagation of ionising radiation from HII regions is able to reproduce the observed distribution of H$α$ in the DIG. Leaking radiation from HII regions also explains the observed decrease in line ratios of low-ionisation species ([SII]/H$α$, [NII]/H$α$ and [OI]/H$α$) with increasing H$α$ surface brightness ($Σ_{Hα}$). Emission from hot low-mass evolved stars, however, is required to explain: (1) the enhanced low-ionisation line ratios observed in the central regions of some galaxies; (2) the observed trends of a flat or decreasing [OIII]/H$β$ with $Σ_{Hα}$; and (3) the offset of some DIG regions from the locus of HII regions in the Baldwin-Phillips-Terlevich (BPT) diagram, extending into the area of low-ionisation (nuclear) emission-line regions (LI[N]ERs). Hot low-mass evolved stars make a small contribution to the energy budget of the DIG (2% of the galaxy-integrated H$α$ emission), but their harder spectra make them fundamental contributors to [OIII] emission. The DIG might result from a superposition of two components, an energetically dominant contribution from young stars and a more diffuse background of harder ionising photons from old stars. This unified framework bridges observations of the Milky Way DIG with LI(N)ER-like emission observed in nearby galaxy bulges.

preprint2022arXiv

A trail of dark matter-free galaxies from a bullet dwarf collision

The ultra-diffuse galaxies DF2 and DF4 in the NGC1052 group share several unusual properties: they both have large sizes, rich populations of overluminous and large globular clusters, and very low velocity dispersions indicating little or no dark matter. It has been suggested that these galaxies were formed in the aftermath of high velocity encounters of gas rich galaxies, events that resemble the collision that created the bullet cluster but on much smaller scales. The gas separates from the dark matter in the collision and subsequent star formation leads to the formation of one or more dark matter-free galaxies. Here we show that the present-day line-of-sight distances and radial velocities of DF2 and DF4 are consistent with their joint formation in the aftermath of a single bullet-dwarf collision, around eight billion years ago. Moreover, we find that DF2 and DF4 are part of an apparent linear substructure of 7-11 large, low-luminosity objects. We propose that these all originated in the same event, forming a trail of dark matter-free galaxies that is more than 2 Mpc long and angled 7 +- 2 degrees from the line of sight. We also tentatively identify the highly dark matter-dominated remnants of the two progenitor galaxies that are expected at the leading edges of the trail.

preprint2022arXiv

Constraining the formation of NGC1052-DF2 from its unusual globular cluster population

The ultra-diffuse galaxy (UDG) NGC1052-DF2 has a low dark matter content and hosts a very unusual globular cluster (GC) population, with a median luminosity $\sim4$ times higher than in most galaxies and containing about 5~per~cent of the galaxy's stars. We apply a theoretical model that predicts the initial cluster mass function as a function of the galactic environment to investigate the origin of DF2's peculiar GC system. Using the GC mass function, the model constrains the star-forming conditions in the galaxy during the formation of its GCs, $\sim9~\rm{Gyr}$ ago. We predict that the GCs formed in an environment with very high gas surface density, $Σ_{\rm ISM} \gtrsim 10^3 M_{\odot} \rm{pc}^{-2}$, and strong centrifugal support, $Ω\gtrsim 0.7~\rm{Myr}^{-1}$, similar to nearby circum-nuclear starbursts and the central region of the Milky Way. These extreme conditions required to form the observed GC population imply a very high cluster formation efficiency of $\gtrsim 78$ per cent, and contrast strongly with the current diffuse nature of the galaxy. Since a nuclear starburst would lead to the rapid in-spiral of the GCs and is ruled out by the absence of a nuclear star cluster, we propose that the GCs plausibly formed during a major merger at $z\sim1.3$. The merger remnant must have undergone significant expansion of its stellar (and perhaps also its dark matter) component to reach its low present surface brightness, leading to the interesting possibility that it was the formation of DF2's extreme GC population that caused it to become a UDG. If true, this strong structural evolution would have important implications for understanding the origins of UDGs.

preprint2022arXiv

Constraining the shape of dark matter haloes with globular clusters

We explore how diffuse stellar light and globular clusters (GCs) can be used to trace the matter distribution of their host halo using an observational methodology. For this, we use 117 simulated dark matter (DM) haloes from the $(34.4~\rm cMpc)^3$ periodic volume of the E-MOSAICS project. For each halo, we compare the stellar surface brightness and GC projected number density maps to the surface densities of DM and total mass. We find that the dominant structures identified in the stellar light and in the GCs correspond closely with those from the DM and total mass. Our method is unaffected by the presence of satellites and its precision improves with fainter GC samples. We recover tight relations between the profiles of stellar surface brightness and GC number density to those of the DM, suggesting that the profile of DM can be accurately recovered from the stars and GCs ($σ\leq0.5~$dex). We quantify the projected morphology of DM, stars and GCs, and find that the stars and GCs are more flattened than the DM. Additionally, the semi-major axes of the distribution of stars and GCs are typically misaligned by $\sim 10~$degrees from that of DM. We demonstrate that deep imaging of diffuse stellar light and GCs can place constraints on the shape, profile and orientation of their host halo. These results extend down to haloes with central galaxies $M_{\star}\geq10^{10}~M_{\odot}$, and the analysis will be applicable to future data from the Euclid, Roman and the Rubin observatories.

preprint2022arXiv

Empirically-motivated early feedback: momentum input by stellar feedback in galaxy simulations inferred through observations

We present a novel method for including the effects of early (pre-supernova) feedback in simulations of galaxy evolution. Rather than building a model which attempts to match idealized, small-scale simulations or analytic approximations, we rely on direct observational measurements of the time-scales over which star-forming molecular clouds are disrupted by early feedback. We combine observations of the spatial de-correlation between molecular gas and star formation tracers on $\sim100$~pc scales with an analytic framework for the expansion of feedback fronts driven by arbitrary sources or mechanisms, and use these to constrain the time-scale and momentum injection rate by early feedback. This allows us to directly inform our model for feedback from these observations, sidestepping the complexity of multiple feedback mechanisms and their interaction below the resolution scale. We demonstrate that this new model has significant effects on the spatial clustering of star formation, the structure of the ISM, and the driving of outflows from the galactic plane, while preserving the overall regulation of the galaxy-integrated star formation rate. We find that this new feedback model results in galaxies that regulate star formation through the rapid disruption of star-forming clouds, rather than by highly efficient, global galactic outflows. We also demonstrate that these results are robust to stochasticity, degraded numerical resolution, changes in the star formation model parameters, and variations in the single free model parameter that is unconstrained by observations.

preprint2022arXiv

Environmental dependence of the molecular cloud lifecycle in 54 main sequence galaxies

The processes of star formation and feedback, regulating the cycle of matter between gas and stars on the scales of giant molecular clouds (GMCs; $\sim$100pc), play a major role in governing galaxy evolution. Measuring the time-scales of GMC evolution is important to identify and characterise the specific physical mechanisms that drive this transition. By applying a robust statistical method to high-resolution CO and narrow-band H$α$ imaging from the PHANGS survey, we systematically measure the evolutionary timeline from molecular clouds to exposed young stellar regions on GMC scales, across the discs of an unprecedented sample of 54 star-forming main-sequence galaxies (excluding their unresolved centres). We find that clouds live for about $1{-}3$ GMC turbulence crossing times ($5{-}30$Myr) and are efficiently dispersed by stellar feedback within $1{-}5$Myr once the star-forming region becomes partially exposed, resulting in integrated star formation efficiencies of $1{-}8$%. These ranges reflect physical galaxy-to-galaxy variation. In order to evaluate whether galactic environment influences GMC evolution, we correlate our measurements with average properties of the GMCs and their local galactic environment. We find several strong correlations that can be physically understood, revealing a quantitative link between galactic-scale environmental properties and the small-scale GMC evolution. Notably, the measured CO-visible cloud lifetimes become shorter with decreasing galaxy mass, mostly due to the increasing presence of CO-dark molecular gas in such environment. Our results represent a first step towards a comprehensive picture of cloud assembly and dispersal, which requires further extension and refinement with tracers of the atomic gas, dust, and deeply-embedded stars.

preprint2022arXiv

Introducing EMP-Pathfinder: modelling the simultaneous formation and evolution of stellar clusters in their host galaxies

The formation and evolution of stellar clusters is intimately linked to that of their host galaxies. To study this connection, we present the EMP-Pathfinder suite of cosmological zoom-in Milky Way-mass simulations. These simulations contain a sub-grid description for stellar cluster formation and evolution, allowing us to study the simultaneous formation and evolution of stellar clusters alongside their host galaxies across cosmic time. As a key ingredient in these simulations, we include the physics of the multi-phase nature of the interstellar medium (ISM), which enables studies of how the presence of a cold, dense ISM affects cluster formation and evolution. We consider two different star formation prescriptions: a constant star formation efficiency per free-fall time, as well as an environmentally-dependent, turbulence-based prescription. We identify two key results drawn from these simulations. Firstly, we find that tidal shock-driven disruption caused by the graininess of the cold ISM produces old ($τ>10~$Gyr) stellar cluster populations with properties that are in excellent agreement with the observed populations in the Milky Way and M31. Importantly, the addition of the cold ISM addresses the areas of disagreement found in previous simulations that lacked the cold gas phase. Secondly, the formation of stellar clusters is extremely sensitive to the baryonic physics that govern the properties of the cold, dense gas reservoir in the galaxy. This implies that the demographics of stellar cluster populations represent an important diagnostic tool for constraining baryonic physics models in upcoming galaxy formation simulations that also include a description of the cold ISM.

preprint2022arXiv

Planetary Nebula Luminosity Function distances for 19 galaxies observed by PHANGS-MUSE

We provide new planetary nebula luminosity function (PNLF) distances to 19 nearby spiral galaxies that were observed with VLT/MUSE by the PHANGS collaboration. Emission line ratios are used to separate planetary nebulae (PNe) from other bright [OIII] emitting sources like compact supernovae remnants (SNRs) or HII regions. While many studies have used narrowband imaging for this purpose, the detailed spectral line information provided by integral field unit (IFU) spectroscopy grants a more robust way of categorising different [OIII] emitters. We investigate the effects of SNR contamination on the PNLF and find that we would fail to classify all objects correctly, when limited to the same data narrowband imaging provides. However, the few misclassified objects usually do not fall on the bright end of the luminosity function, and only in three cases does the distance change by more than $1σ$. We find generally good agreement with literature values from other methods. Using metallicity constraints that have also been derived from the same IFU data, we revisit the PNLF zero point calibration. Over a range of $8.34 < 12 + \log(\mathrm{O}/\mathrm{H}) < 8.59$, our sample is consistent with a constant zero point and yields $M^*=-4.542^{+0.103}_{-0.059}\, \mathrm{mag}$, within $1σ$ of other literature values. MUSE pushes the limits of PNLF studies and makes galaxies beyond $20\, \mathrm{Mpc}$ accessible for this kind of analysis. This approach to the PNLF shows great promise for leveraging existing archival IFU data on nearby galaxies.

preprint2022arXiv

The emergence of dark matter-deficient ultra-diffuse galaxies driven by scatter in the stellar mass-halo mass relation and feedback from globular clusters

In addition to their low stellar densities, ultra-diffuse galaxies (UDGs) have a broad variety of dynamical mass-to-light ratios, ranging from dark matter (DM) dominated systems to objects nearly devoid of DM. To investigate the origin of this diversity, we develop a simple, semi-empirical model that predicts the structural evolution of galaxies, driven by feedback from massive star clusters, as a function of their departure from the mean SMHM relation. The model predicts that a galaxy located $\gtrsim 0.5$ dex above the mean relation at $M_{\rm halo}=10^{10}M_{\odot}$ will host a factor of $\sim 10-100$ larger globular cluster (GC) populations, and that feedback from these GCs drives a significant expansion of the stellar component and loss of DM compared to galaxies on the SMHM relation. This effect is stronger in haloes that collapse earlier and have enhanced star formation rates at $z\gtrsim2$, which leads to increased gas pressures, stellar clustering, and mean cluster masses, and significantly enhances the energy loading of galactic winds and its impact on the DM and stellar orbits. The impact on galaxy size and DM content can be large enough to explain observed galaxies that contain nearly the universal baryon fraction, as well as NGC1052-DF2 and DF4 and other isolated UDGs that contain almost no DM. The trend of increasing galaxy size with GC specific frequency observed in galaxy clusters also emerges naturally in the model. Our predictions can be tested with large and deep surveys of the stellar and GC populations in dwarfs and UDGs. Because stellar clustering drives the efficiency of galactic winds, it may be a dominant factor in the structural evolution of galaxies and should be included as an essential ingredient in galaxy formation models.

preprint2022arXiv

The Gas-Star Formation Cycle in Nearby Star-forming Galaxies II. Resolved Distributions of CO and H$α$ Emission for 49 PHANGS Galaxies

The relative distribution of molecular gas and star formation in galaxies gives insight into the physical processes and timescales of the cycle between gas and stars. In this work, we track the relative spatial configuration of CO and H$α$ emission at high resolution in each of our galaxy targets, and use these measurements to quantify the distributions of regions in different evolutionary stages of star formation: from molecular gas without star formation traced by H$α$ to star-forming gas, and to HII regions. The large sample, drawn from the Physics at High Angular resolution in Nearby GalaxieS ALMA and narrowband H$α$ (PHANGS-ALMA and PHANGS-H$α$) surveys, spans a wide range of stellar mass and morphological types, allowing us to investigate the dependencies of the gas-star formation cycle on global galaxy properties. At a resolution of 150 pc, the incidence of regions in different stages shows a dependence on stellar mass and Hubble type of galaxies over the radial range probed. Massive and/or earlier-type galaxies exhibit a significant reservoir of molecular gas without star formation traced by H$α$, while lower-mass galaxies harbor substantial HII regions that may have dispersed their birth clouds or formed from low-mass, more isolated clouds. Galactic structures add a further layer of complexity to relative distribution of CO and H$α$ emission. Trends between galaxy properties and distributions of gas traced by CO and H$α$ are visible only when the observed spatial scale is $\ll$ 500 pc, reflecting the critical resolution requirement to distinguish stages of star formation process.

preprint2022arXiv

The initial conditions for young massive cluster formation in the Galactic Centre: convergence of large-scale gas flows

Young massive clusters (YMCs) are compact ($\lesssim$1 pc), high-mass (>10${}^4$ M${}_{\odot}$) stellar systems of significant scientific interest. Due to their rarity and rapid formation, we have very few examples of YMC progenitor gas clouds before star formation has begun. As a result, the initial conditions required for YMC formation are uncertain. We present high-resolution (0.13$^{\prime\prime}$, $\sim$1000 au) ALMA observations and Mopra single-dish data, showing that Galactic Centre dust ridge `Cloud d' (G0.412$+$0.052, mass$\sim 7.6 \times 10^4$ M$_{\odot}$, radius$\sim 3.2$ pc) has the potential to become an Arches-like YMC (10$^4$ M$_{\odot}$, r$\sim$1 pc), but is not yet forming stars. This would mean it is the youngest known pre-star forming massive cluster and therefore could be an ideal laboratory for studying the initial conditions of YMC formation. We find 96 sources in the dust continuum, with masses $\lesssim$3 M$_{\odot}$ and radii of $\sim$10${}^3$ au. The source masses and separations are more consistent with thermal rather than turbulent fragmentation. It is not possible to unambiguously determine the dynamical state of most of the sources, as the uncertainty on virial parameter estimates is large. We find evidence for large-scale ($\sim$1 pc) converging gas flows, which could cause the cloud to grow rapidly, gaining 10$^4$ M$_{\odot}$ within 10$^5$ yr. The highest density gas is found at the convergent point of the large-scale flows. We expect this cloud to form many high-mass stars, but find no high-mass starless cores. If the sources represent the initial conditions for star formation, the resulting IMF will be bottom-heavy.

preprint2022arXiv

The PHANGS-MUSE survey -- Probing the chemo-dynamical evolution of disc galaxies

We present the PHANGS-MUSE survey, a programme using the MUSE IFS at the ESO VLT to map 19 massive $(9.4 < \log(M_{*}/M_\odot) < 11.0)$ nearby (D < 20 Mpc) star-forming disc galaxies. The survey consists of 168 MUSE pointings (1'x1' each), a total of nearly 15 Million spectra, covering ~1.5 Million independent spectra. PHANGS-MUSE provides the first IFS view of star formation across different local environments (including galaxy centres, bars, spiral arms) in external galaxies at a median resolution of 50~pc, better than the mean inter-cloud distance in the ionised interstellar medium. This `cloud-scale' resolution allows detailed demographics and characterisations of HII regions and other ionised nebulae. PHANGS-MUSE further delivers a unique view on the associated gas and stellar kinematics, and provides constraints on the star formation history. The PHANGS-MUSE survey is complemented by dedicated ALMA CO(2-1) and multi-band HST observations, therefore allowing us to probe the key stages of the star formation process from molecular clouds to HII regions and star clusters. This paper describes the scientific motivation, sample selection, observational strategy, data reduction and analysis process of the PHANGS-MUSE survey. We present our bespoke automated data-reduction framework, which is built on the reduction recipes provided by ESO, but additionally allows for mosaicking and homogenisation of the point spread function. We further present a detailed quality assessment and a brief illustration of the potential scientific applications of the large set of PHANGS-MUSE data products generated by our data analysis framework. The data cubes and analysis data products described in this paper represent the basis for the first PHANGS-MUSE public data release and are available in the ESO archive and via the Canadian Astronomy Data Centre.

preprint2022arXiv

The present-day globular cluster kinematics of lenticular galaxies from the E-MOSAICS simulations and their relation to the galaxy assembly histories

We study the present-day rotational velocity ($V_{rot}$) and velocity dispersion ($σ$) profiles of the globular cluster (GC) systems in a sample of 50 lenticular (S0) galaxies from the E-MOSAICS galaxy formation simulations. We find that 82% of the galaxies have GCs that are rotating along the photometric major axis of the galaxy ($aligned$), while the remaining 18% of the galaxies do not ($misaligned$). This is generally consistent with the observations from the SLUGGS survey. For the $aligned$ galaxies, classified as $peaked$ $and$ $outwardly$ $decreasing$ (49%), $flat$ (24%) and $increasing$ (27%) based on the $V_{rot}/σ$ profiles out to large radii, we do not find any clear correlation between these present-day $V_{rot}/σ$ profiles of the GCs and the past merger histories of the S0 galaxies, unlike in previous simulations of galaxy stars. For just over half of the $misaligned$ galaxies, we find that the GC misalignment is the result of a major merger within the last 10 Gyr so that the $ex$-$situ$ GCs are misaligned by an angle between 0° (co-rotation) to 180° (counter-rotation) with respect to the $in$-$situ$ GCs, depending on the orbital configuration of the merging galaxies. For the remaining $misaligned$ galaxies, we suggest that the $in$-$situ$ metal-poor GCs, formed at early times, have undergone more frequent kinematic perturbations than the $in$-$situ$ metal-rich GCs. We also find that the GCs accreted early and the $in$-$situ$ GCs are predominantly located within 0.2 virial radii ($R_{200}$) from the centre of galaxies in 3D phase-space diagrams.

preprint2022arXiv

Uncertainties in supernova input rates drive qualitative differences in simulations of galaxy evolution

Feedback from core collapse supernovae (SNe), the final stage of evolution of massive stars, is a key element in galaxy formation theory. The energy budget of SN feedback, as well as the duration over which SNe occur, are constrained by stellar lifetime models and the minimum mass star that ends its life as a SN. Simplifying approximations for this SN rate are ubiquitous in simulation studies. We show here how the choice of SN budget and timings ($t_0$ for the delay between star formation and the first SN, $τ_{\rm SN}$ for the duration of SN injection, and the minimum SN progenitor mass) drive changes in the regulation of star formation and outflow launching. Extremely long delays for instantaneous injection of SN energy $(t_0 << 20\;\rm{Myr})$ reduces star formation and drive stronger outflows compared smaller delays. This effect is primarily driven by enhanced clustering of young stars. With continuous injection of energy, longer SN durations results in a larger fraction of SN energy deposited in low ambient gas densities, where cooling losses are lower. This is effect is particularly when driven by the choice of the minimum SN progenitor mass, which also sets the total SN energy budget. These underlying uncertainties mean that despite advances in the sub-grid modeling of SN feedback, serious difficulties in constraining the strength of SN feedback remain. We recommend future simulations use realistic SN injection durations, and bound their results using SN energy budgets and durations for minimum SN progenitors of $7M_\odot$ and $9M_\odot$.

preprint2022arXiv

WISDOM Project -- X. The morphology of the molecular ISM in galaxy centres and its dependence on galaxy structure

We use high-resolution maps of the molecular interstellar medium (ISM) in the centres of eighty-six nearby galaxies from the millimetre-Wave Interferometric Survey of Dark Object Masses (WISDOM) and Physics at High Angular Resolution in Nearby GalaxieS (PHANGS) surveys to investigate the physical mechanisms setting the morphology of the ISM at molecular cloud scales. We show that early-type galaxies tend to have smooth, regular molecular gas morphologies, while the ISM in spiral galaxy bulges is much more asymmetric and clumpy when observed at the same spatial scales. We quantify these differences using non-parametric morphology measures (Asymmetry, Smoothness and Gini), and compare these measurements with those extracted from idealised galaxy simulations. We show that the morphology of the molecular ISM changes systematically as a function of various large-scale galaxy parameters, including galaxy morphological type, stellar mass, stellar velocity dispersion, effective stellar mass surface density, molecular gas surface density, star formation efficiency and the presence of a bar. We perform a statistical analysis to determine which of these correlated parameters best predicts the morphology of the ISM. We find the effective stellar mass surface (or volume) density to be the strongest predictor of the morphology of the molecular gas, while star formation and bars maybe be important secondary drivers. We find that gas self-gravity is not the dominant process shaping the morphology of the molecular gas in galaxy centres. Instead effects caused by the depth of the potential well such as shear, suppression of stellar spiral density waves and/or inflow affect the ability of the gas to fragment.

preprint2022arXiv

WISDOM project -- XI. Star Formation Efficiency in the Bulge of the AGN-host Galaxy NGC 3169 with SITELLE and ALMA

The star formation efficiency (SFE) has been shown to vary across different environments, particularly within galactic starbursts and deep within the bulges of galaxies. Various quenching mechanisms may be responsible, ranging from galactic dynamics to feedback from active galactic nuclei (AGN). Here, we use spatially-resolved observations of warm ionised gas emission lines from the imaging Fourier transform spectrograph SITELLE at the Canada-France-Hawaii Telescope (CFHT) and cold molecular gas (CO(2-1)) from the Atacama Large Millimeter/sub-millimeter Array (ALMA) to study the SFE in the bulge of the AGN-host galaxy NGC 3169. After distinguishing star-forming regions from AGN-ionised regions using emission-line ratio diagnostics, we measure spatially-resolved molecular gas depletion times (τ_dep = 1/SFE) with a spatial resolution of \approx 100 pc within a galactocentric radius of 1.8 kpc. We identify a star-forming ring located at radii 1.25 \pm 0.6 kpc with an average τ_dep of 0.3 Gyr. At radii < 0.9 kpc, however, the molecular gas surface densities and depletion times increase with decreasing radius, the latter reaching approximately 2.3 Gyr at a radius \approx 500 pc. Based on analyses of the gas kinematics and comparisons with simulations, we identify AGN feedback, bulge morphology and dynamics as the possible causes of the radial profile of SFE observed in the central region of NGC 3169.

preprint2021arXiv

Age Determination of Galaxy Merger Remnant Stars using Asteroseismology

The Milky Way was shaped by the mergers with several galaxies in the past. We search for remnant stars that were born in these foreign galaxies and assess their ages in an effort to put upper limits on the merger times and thereby better understand the evolutionary history of our Galaxy. Using 6D-phase space information from Gaia eDR3 and chemical information from APOGEE DR16, we kinematically and chemically select $23$ red giant stars belonging to former dwarf galaxies that merged with the Milky Way. With added asteroseismology from Kepler and K2, we determine the ages of the $23$ ex-situ stars and $55$ in-situ stars with great precision. We find that all the ex-situ stars are consistent with being older than $8$ Gyr. While it is not possible to associate all the stars with a specific dwarf galaxy we classify eight of them as Gaia-Enceladus/Sausage stars, which is one of the most massive mergers in our Galaxy's history. We determine their mean age to be $9.5^{+1.2}_{-1.3}$ Gyr consistent with a merger time of $8$-$10$ Gyr ago. The rest of the stars are possibly associated with Kraken, Thamnos, Sequoia, or another extragalactic progenitor. The age determination of ex-situ stars paves the way to more accurately pinning down when the merger events occurred and hence provide tight constraints useful for simulating how these events unfolded.

preprint2021arXiv

Giant Molecular Cloud Catalogues for PHANGS-ALMA: Methods and Initial Results

We present improved methods for segmenting CO emission from galaxies into individual molecular clouds, providing an update to the CPROPS algorithms presented by Rosolowsky & Leroy (2006; arXiv:astro-ph/0601706 ). The new code enables both homogenization of the noise and spatial resolution among data, which allows for rigorous comparative analysis. The code also models the completeness of the data via false source injection and includes an updated segmentation approach to better deal with blended emission. These improved algorithms are implemented in a publicly available python package, PYCPROPS. We apply these methods to ten of the nearest galaxies in the PHANGS-ALMA survey, cataloguing CO emission at a common 90 pc resolution and a matched noise level. We measure the properties of 4986 individual clouds identified in these targets. We investigate the scaling relations among cloud properties and the cloud mass distributions in each galaxy. The physical properties of clouds vary among galaxies, both as a function of galactocentric radius and as a function of dynamical environment. Overall, the clouds in our target galaxies are well-described by approximate energy equipartition, although clouds in stellar bars and galaxy centres show elevated line widths and virial parameters. The mass distribution of clouds in spiral arms has a typical mass scale that is 2.5x larger than interarm clouds and spiral arms clouds show slightly lower median virial parameters compared to interarm clouds (1.2 versus 1.4).

preprint2021arXiv

Globular clusters as tracers of the dark matter halo: insights from the E-MOSAICS simulations

Globular clusters (GCs) are bright objects that span a wide range of galactocentric distances, and are thus probes of the structure of dark matter (DM) haloes. In this work, we explore whether the projected radial profiles of GCs can be used to infer the structural properties of their host DM haloes. We use the simulated GC populations in a sample of 166 central galaxies from the $(34.4~\rm cMpc)^3$ periodic volume of the E-MOSAICS project. We find that more massive galaxies host stellar and GC populations with shallower density profiles that are more radially extended. In addition, the metal-poor GC subpopulations tend to have shallower and more extended profiles than the metal-rich subsamples, which we relate to the preferentially accreted origin of the metal-poor GCs. We find strong correlations between the slopes and effective radii of the radial profiles of the GC populations and the structural properties of the DM haloes, such as their power-law slopes, scale radii, and concentration parameters. Accounting for a dependence on the galaxy stellar mass decreases the scatter of the two-dimensional relations. This suggests that the projected number counts of GCs, combined with their galaxy mass, trace the density profile of the DM halo of their host galaxy. When applied to extragalactic GC systems, we recover the scale radii and the extent of the DM haloes of a sample of ETGs with uncertainties smaller than $0.2~\rm dex$. Thus, extragalactic GC systems provide a novel avenue to explore the structure of DM haloes beyond the Local Group.

preprint2021arXiv

Low-J CO Line Ratios From Single Dish CO Mapping Surveys and PHANGS-ALMA

We measure the low-J CO line ratio R21=CO(2-1)/CO(1-0), R32=CO(3-2)/CO(2-1), and R31 = CO(3-2)/CO(1-0) using whole-disk CO maps of nearby galaxies. We draw CO(2-1) from PHANGS--ALMA, HERACLES, and follow-up IRAM surveys; CO(1-0) from COMING and the Nobeyama CO Atlas of Nearby Spiral Galaxies; and CO(3-2) from the JCMT NGLS and APEX LASMA mapping. Altogether this yields 76, 47, and 29 maps of R21, R32, and R31 at 20" \sim 1.3 kpc resolution, covering 43, 34, and 20 galaxies. Disk galaxies with high stellar mass, log10 M_* [Msun]=10.25-11 and star formation rate, SFR=1-5 Msun/yr, dominate the sample. We find galaxy-integrated mean values and 16%-84% range of R21 = 0.65 (0.50-0.83), R32=0.50 (0.23-0.59), and R31=0.31 (0.20-0.42). We identify weak trends relating galaxy-integrated line ratios to properties expected to correlate with excitation, including SFR/M_* and SFR/L_CO. Within galaxies, we measure central enhancements with respect to the galaxy-averaged value of \sim 0.18^{+0.09}_{-0.14} dex for R21, 0.27^{+0.13}_{-0.15} dex for R31, and 0.08^{+0.11}_{-0.09} dex for R32. All three line ratios anti-correlate with galactocentric radius and positively correlate with the local star formation rate surface density and specific star formation rate, and we provide approximate fits to these relations. The observed ratios can be reasonably reproduced by models with low temperature, moderate opacity, and moderate densities, in good agreement with expectations for the cold ISM. Because the line ratios are expected to anti-correlate with the CO(1-0)-to-H_2 conversion factor, alphaCO^(1-0), these results have general implications for the interpretation of CO emission from galaxies.

preprint2021arXiv

Molecular Gas Properties and CO-to-H2 Conversion Factors in the Central Kiloparsec of NGC 3351

The CO-to-H$_2$ conversion factor ($α_\rm{CO}$) is critical to studying molecular gas and star formation in galaxies. The value of $α_\rm{CO}$ has been found to vary within and between galaxies, but the specific environmental conditions that cause these variations are not fully understood. Previous observations on $\sim$kpc scales revealed low values of $α_\rm{CO}$ in the centers of some barred spiral galaxies, including NGC 3351. We present new ALMA Band 3, 6, and 7 observations of $^{12}$CO, $^{13}$CO, and C$^{18}$O lines on 100 pc scales in the inner $\sim$2 kpc of NGC 3351. Using multi-line radiative transfer modeling and a Bayesian likelihood analysis, we infer the H$_2$ density, kinetic temperature, CO column density per line width, and CO isotopologue abundances on a pixel-by-pixel basis. Our modeling implies the existence of a dominant gas component with a density of $2{-}3\times10^3$ $\rm{cm^{-3}}$ in the central ${\sim}$1 kpc and a high temperature of 30$-$60 K near the nucleus and near the contact points that connect to the bar-driven inflows. Assuming a CO/H$_2$ abundance of $3\times10^{-4}$, our analysis yields $α_\rm{CO}{\sim}0.5{-}2.0$ $\rm{M_\odot\,(K~km~s^{-1}~pc^2)^{-1}}$ with a decreasing trend with galactocentric radius in the central $\sim$1 kpc. The inflows show a substantially lower $α_\rm{CO} < 0.1$ $\rm{M_\odot\,(K~km~s^{-1}~pc^2)^{-1}}$, likely due to lower optical depths caused by turbulence or shear in the inflows. Over the whole region, this gives an intensity-weighted $α_\rm{CO}$ of ${\sim}1.5$ $\rm{M_\odot\,(K~km~s^{-1}~pc^2)^{-1}}$, which is similar to previous dust modeling based results at kpc scales. This suggests that low $α_\rm{CO}$ on kpc scales in the centers of some barred galaxies may be due to the contribution of low optical depth CO emission in bar-driven inflows.

preprint2021arXiv

NGC5846-UDG1: A galaxy formed mostly by star formation in massive, extremely dense clumps of gas

It has been shown that ultra-diffuse galaxies (UDGs) have higher specific frequencies of globular clusters on average than other dwarf galaxies with similar luminosities. The UDG NGC5846-UDG1 is among the most extreme examples of globular cluster-rich galaxies found so far. Here we present new Hubble Space Telescope (HST) observations and analysis of this galaxy and its globular cluster system. We find that NGC5846-UDG1 hosts $54 \pm 9$ globular clusters, three to four times more than any previously known galaxy with a similar luminosity, and higher than reported in previous studies. With a galaxy luminosity of $L_{V,\mathrm{gal}} \approx 6 \times 10^7\,{\rm L}_{\odot}$ ($M_\star \approx 1.2 \times 10^8\,{\rm M}_\odot$) and a total globular cluster luminosity of $L_{V,\mathrm{GCs}} \approx 7.6 \times 10^6\,{\rm L}_{\odot}$, we find that the clusters currently comprise $\sim 13 \%$ of the total light. Taking into account the effects of mass loss from clusters during their formation and throughout their lifetime, we infer that most of the stars in the galaxy likely formed in globular clusters, and very little to no "normal" low-density star formation occurred. This result implies that the most extreme conditions during early galaxy formation promoted star formation in massive and dense clumps, in contrast to the dispersed star formation observed in galaxies today.

preprint2021arXiv

PHANGS-HST: Star Cluster Spectral Energy Distribution Fitting with CIGALE

The sensitivity and angular resolution of photometric surveys executed by the Hubble Space Telescope (HST) enable studies of individual star clusters in galaxies out to a few tens of megaparsecs. The fitting of spectral energy distributions (SEDs) of star clusters is essential for measuring their physical properties and studying their evolution. We report on the use of the publicly available Code Investigating GALaxy Emission (CIGALE) SED fitting package to derive ages, stellar masses, and reddenings for star clusters identified in the Physics at High Angular resolution in Nearby GalaxieS-HST (PHANGS-HST) survey. Using samples of star clusters in the galaxy NGC 3351, we present results of benchmark analyses performed to validate the code and a comparison to SED fitting results from the Legacy ExtraGalactic Ultraviolet Survey (LEGUS). We consider procedures for the PHANGS-HST SED fitting pipeline, e.g., the choice of single stellar population models, the treatment of nebular emission and dust, and the use of fluxes versus magnitudes for the SED fitting. We report on the properties of clusters in NGC 3351 and find, on average, the clusters residing in the inner star-forming ring of NGC 3351 are young ($< 10$ Myr) and massive ($10^{5} M_{\odot}$) while clusters in the stellar bulge are significantly older. Cluster mass function fits yield $β$ values around -2, consistent with prior results with a tendency to be shallower at the youngest ages. Finally, we explore a Bayesian analysis with additional physically-motivated priors for the distribution of ages and masses and analyze the resulting cluster distributions.

preprint2021arXiv

PHANGS-MUSE: the HII region luminosity function of local star-forming galaxies

We use an unprecedented sample of about 23,000 HII regions detected at an average physical resolution of 67pc in the PHANGS-MUSE sample to study the extragalactic HII region Ha luminosity function (LF). Our observations probe the star-forming disk of 19 nearby spiral galaxies with low inclination and located close to the star formation main sequence at z=0. The mean LF slope $α$ in our sample is =1.73 with a $σ$ of 0.15. We find that $α$ decreases with the galaxy's star formation rate surface density and argue that this is driven by an enhanced clustering of young stars at high gas surface densities. Looking at the HII regions within single galaxies we find that no significant variations occur between the LF of the inner and outer part of the star-forming disk, whereas the LF in the spiral arm areas is shallower than in the inter-arm areas for six out of the 13 galaxies with clearly visible spiral arms. We attribute these variations to the spiral arms increasing the molecular clouds arm--inter-arm mass contrast and find suggestive evidence that they are more evident for galaxies with stronger spiral arms. Furthermore, we find systematic variations in $α$ between samples of HII regions with high and low ionization parameter q and argue that they are driven by the aging of HII regions.

preprint2021arXiv

Star formation in 'the Brick': ALMA reveals an active proto-cluster in the Galactic centre cloud G0.253+0.016

G0.253+0.016, aka 'the Brick', is one of the most massive (> 10^5 Msun) and dense (> 10^4 cm-3) molecular clouds in the Milky Way's Central Molecular Zone. Previous observations have detected tentative signs of active star formation, most notably a water maser that is associated with a dust continuum source. We present ALMA Band 6 observations with an angular resolution of 0.13" (1000 AU) towards this 'maser core', and report unambiguous evidence of active star formation within G0.253+0.016. We detect a population of eighteen continuum sources (median mass ~ 2 Msun), nine of which are driving bi-polar molecular outflows as seen via SiO (5-4) emission. At the location of the water maser, we find evidence for a protostellar binary/multiple with multi-directional outflow emission. Despite the high density of G0.253+0.016, we find no evidence for high-mass protostars in our ALMA field. The observed sources are instead consistent with a cluster of low-to-intermediate-mass protostars. However, the measured outflow properties are consistent with those expected for intermediate-to-high-mass star formation. We conclude that the sources are young and rapidly accreting, and may potentially form intermediate and high-mass stars in the future. The masses and projected spatial distribution of the cores are generally consistent with thermal fragmentation, suggesting that the large-scale turbulence and strong magnetic field in the cloud do not dominate on these scales, and that star formation on the scale of individual protostars is similar to that in Galactic disc environments.

preprint2021arXiv

The PHANGS-HST Survey: Physics at High Angular resolution in Nearby GalaxieS with the Hubble Space Telescope

The PHANGS program is building the first dataset to enable the multi-phase, multi-scale study of star formation across the nearby spiral galaxy population. This effort is enabled by large survey programs with ALMA, VLT/MUSE, and HST, with which we have obtained CO(2-1) imaging, optical spectroscopic mapping, and high resolution UV-optical imaging, respectively. Here, we present PHANGS-HST, which has obtained five band NUV-U-B-V-I imaging of the disks of 38 spiral galaxies at distances of 4-23 Mpc, and parallel V and I band imaging of their halos, to provide a census of tens of thousands of compact star clusters, and multi-scale stellar associations. The combination of HST, ALMA, and VLT/MUSE observations will yield an unprecedented joint catalog of the observed and physical properties of ~100,000 star clusters, associations, HII regions, and molecular clouds. With these basic units of star formation, PHANGS will systematically chart the evolutionary cycling between gas and stars, across a diversity of galactic environments found in nearby galaxies. We discuss the design of the PHANGS-HST survey, and provide an overview of the HST data processing pipeline and first results, highlighting new methods for selecting star cluster candidates, morphological classification of candidates with convolutional neural networks, and identification of stellar associations over a range of physical scales with a watershed algorithm. We describe the cross-observatory imaging, catalogs, and software products to be released. These high-level science products will seed a broad range of investigations, in particular, the study of embedded stellar populations and dust with JWST, for which a PHANGS Cycle 1 Treasury program to obtain eight band 2-21 $μ$m imaging has been approved.

preprint2021arXiv

The physics governing the upper truncation mass of the globular cluster mass function

The mass function of globular cluster (GC) populations is a fundamental observable that encodes the physical conditions under which these massive stellar clusters formed and evolved. The high-mass end of star cluster mass functions are commonly described using a Schechter function, with an exponential truncation mass $M_{c,*}$. For the GC mass functions in the Virgo galaxy cluster, this truncation mass increases with galaxy mass ($M_{*}$). In this paper we fit Schechter mass functions to the GCs in the most massive galaxy group ($M_{\mathrm{200}} = 5.14 \times 10^{13} M_{\odot}$) in the E-MOSAICS simulations. The fiducial cluster formation model in E-MOSAICS reproduces the observed trend of $M_{c,*}$ with $M_{*}$ for the Virgo cluster. We therefore examine the origin of the relation by fitting $M_{c,*}$ as a function of galaxy mass, with and without accounting for mass loss by two-body relaxation, tidal shocks and/or dynamical friction. In the absence of these mass-loss mechanisms, the $M_{c,*}$-$M_{*}$ relation is flat above $M_* > 10^{10} M_{\odot}$. It is therefore the disruption of high-mass GCs in galaxies with $M_{*}\sim 10^{10} M_{\odot}$ that lowers the $M_{c,*}$ in these galaxies. High-mass GCs are able to survive in more massive galaxies, since there are more mergers to facilitate their redistribution to less-dense environments. The $M_{c,*}-M_*$ relation is therefore a consequence of both the formation conditions of massive star clusters and their environmentally-dependent disruption mechanisms.

preprint2021arXiv

What to expect when using globular clusters as tracers of the total mass distribution in Milky Way-mass galaxies

Dynamical models allow us to connect the motion of a set of tracers to the underlying gravitational potential, and thus to the total (luminous and dark) matter distribution. They are particularly useful for understanding the mass and spatial distribution of dark matter (DM) in a galaxy. Globular clusters (GCs) are an ideal tracer population in dynamical models, since they are bright and can be found far out into the halo of galaxies. We aim to test how well Jeans-Anisotropic-MGE (JAM) models using GCs (positions and line-of-sight velocities) as tracers can constrain the mass and radial distribution of DM halos. For this, we use the E-MOSAICS suite of 25 zoom-in simulations of L* galaxies. We find that the DM halo properties are reasonably well recovered by the JAM models. There is, however, a strong correlation between how well we recover the mass and the radial distribution of the DM and the number of GCs in the galaxy: the constraints get exponentially worse with fewer GCs, and at least 150 GCs are needed in order to guarantee that the JAM model will perform well. We find that while the data quality (uncertainty on the radial velocities) can be important, the number of GCs is the dominant factor in terms of the accuracy and precision of the measurements. This work shows promising results for these models to be used in extragalactic systems with a sample of more than 150 GCs.

preprint2020arXiv

A model for the onset of self-gravitation and star formation in molecular gas governed by galactic forces: II. the bottleneck to collapse set by cloud-environment decoupling

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.

preprint2020arXiv

AlFoCS + Fornax3D: resolved star formation in the Fornax cluster with ALMA and MUSE

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.

preprint2020arXiv

ALMA Observations of Massive Clouds in the Central Molecular Zone: Jeans Fragmentation and Cluster Formation

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.

preprint2020arXiv

An uncertainty principle for star formation -- V. The influence of dust extinction on star formation rate tracer lifetimes and the inferred molecular cloud lifecycle

Recent observational studies aiming to quantify the molecular cloud lifecycle require the use of known 'reference time-scales' to turn the relative durations of different phases of the star formation process into absolute time-scales. We previously constrained the characteristic emission time-scales of different star formation rate (SFR) tracers, as a function of the SFR surface density and metallicity. However, we omitted the effects of dust extinction. Here, we extend our suite of SFR tracer emission time-scales by accounting for extinction, using synthetic emission maps of a high-resolution hydrodynamical simulation of an isolated, Milky-Way-like disc galaxy. The stellar feedback included in the simulation is inefficient compared to observations, implying that it represents a limiting case in which the duration of embedded star formation (and the corresponding effect of extinction) is overestimated. Across our experiments, we find that extinction mostly decreases the SFR tracer emission time-scale, changing the time-scales by factors of 0.04-1.74, depending on the gas column density. UV filters are more strongly affected than H$α$ filters. We provide the limiting correction factors as a function of the gas column density and flux sensitivity limit for a wide variety of SFR tracers. Applying these factors to observational characterisations of the molecular cloud lifecycle produces changes that broadly fall within the quoted uncertainties, except at high kpc-scale gas surface densities ($Σ_{\rm g}\gtrsim20~{\mathrm{M_{\odot}\,pc^{-2}}}$). Under those conditions, correcting for extinction may decrease the measured molecular cloud lifetimes and feedback time-scales, which further strengthens previous conclusions that molecular clouds live for a dynamical time and are dispersed by early, pre-supernova feedback.

preprint2020arXiv

CMZoom: Survey Overview and First Data Release

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.

preprint2020arXiv

Deep Transfer Learning for Star Cluster Classification: I. Application to the PHANGS-HST Survey

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 (D < 20 Mpc) in the PHANGS-HST survey. Given the relatively small nature of existing, human-labelled star cluster samples, we transfer the knowledge of state-of-the-art neural network models for real-object recognition to classify star clusters candidates into four morphological classes. We perform a series of experiments to determine the dependence of classification performance on: neural network architecture (ResNet18 and VGG19-BN); training data sets curated by either a single expert or three astronomers; and the size of the images used for training. We find that the overall classification accuracies are not significantly affected by these choices. The networks are used to classify star cluster candidates in the PHANGS-HST galaxy NGC 1559, which was not included in the training samples. The resulting prediction accuracies are 70%, 40%, 40-50%, 50-70% for class 1, 2, 3 star clusters, and class 4 non-clusters respectively. This performance is competitive with consistency achieved in previously published human and automated quantitative classification of star cluster candidate samples (70-80%, 40-50%, 40-50%, and 60-70%). The methods introduced herein lay the foundations to automate classification for star clusters at scale, and exhibit the need to prepare a standardized dataset of human-labelled star cluster classifications, agreed upon by a full range of experts in the field, to further improve the performance of the networks introduced in this study.

preprint2020arXiv

Distances to PHANGS Galaxies: New Tip of the Red Giant Branch Measurements and Adopted Distances

PHANGS-HST is an ultraviolet-optical imaging survey of 38 spiral galaxies within ~20 Mpc. Combined with the PHANGS-ALMA, PHANGS-MUSE surveys and other multiwavelength data, the dataset will provide an unprecedented look into the connections between young stars, HII regions, and cold molecular gas in these nearby star-forming galaxies. Accurate distances are needed to transform measured observables into physical parameters (e.g., brightness to luminosity, angular to physical sizes of molecular clouds, star clusters and associations). PHANGS-HST has obtained parallel ACS imaging of the galaxy halos in the F606W and F814W bands. Where possible, we use these parallel fields to derive tip of the red giant branch (TRGB) distances to these galaxies. In this paper, we present TRGB distances for 11 galaxies from ~4 to ~15 Mpc, based on the first year of PHANGS-HST observations. Five of these represent the first published TRGB distance measurements (IC 5332, NGC 2835, NGC 4298, NGC 4321, and NGC 4328), and eight of which are the best available distances to these targets. We also provide a compilation of distances for the 118 galaxies in the full PHANGS sample, which have been adopted for the first PHANGS-ALMA public data release.

preprint2020arXiv

Dynamical Equilibrium in the Molecular ISM in 28 Nearby Star-Forming Galaxies

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, $Σ_\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 $Σ_\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.

preprint2020arXiv

Entropy Driven Winds: Outflows and Fountains Lifted Gently by Buoyancy

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.

preprint2020arXiv

From diffuse gas to dense molecular cloud cores

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.

preprint2020arXiv

Heart of Darkness: the influence of galactic dynamics on quenching star formation in galaxy spheroids

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.

preprint2020arXiv

Impact of low-energy cosmic rays on star formation

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 $γ$-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.

preprint2020arXiv

Not all stars form in clusters -- $Gaia$-DR2 uncovers the origin of OB associations

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.

preprint2020arXiv

PHANGS CO kinematics: disk orientations and rotation curves at 150 pc resolution

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.

preprint2020arXiv

Physical Processes in Star Formation

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.

preprint2020arXiv

Star clusters near and far; tracing star formation across cosmic time

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.

preprint2020arXiv

Stellar Feedback and Resolved Stellar IFU Spectroscopy in the nearby Spiral Galaxy NGC 300

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.

preprint2020arXiv

Testing viscous disc theory using the balance between stellar accretion and external photoevaporation of protoplanetary discs

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 $σ$ 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.

preprint2020arXiv

The CO-dark molecular gas mass in 30 Doradus

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.

preprint2020arXiv

The globular cluster system mass-halo mass relation in the E-MOSAICS simulations

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 (Mgc) and the dark matter halo mass (Mhalo) 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.4$~cMpc 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 Mgc-Mhalo relation is linear for halo masses $>5\times10^{11}~Msun$, and is driven by the hierarchical assembly of galaxies. 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-halo mass relation of galaxies. Our fiducial model reproduces the observed Mgc-Mstar relation across the full mass range, which we argue is more physically relevant than the Mgc-Mhalo relation. We also explore the physical processes driving the observed constant value of $Mgc / Mhalo \sim 5\times10^{-5}$ and find that it is the result of a combination of cluster formation physics and cluster disruption.

preprint2020arXiv

The molecular cloud lifecycle

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.

preprint2020arXiv

The Origin of the Stellar Mass Distribution and Multiplicity

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.

preprint2020arXiv

The Physics of Star Cluster Formation and Evolution

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.

preprint2020arXiv

The role of galactic dynamics in shaping the physical properties of giant molecular clouds in Milky Way-like galaxies

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.

preprint2020arXiv

Towards a multi-tracer timeline of star formation in the LMC -- I.\ Deriving the lifetimes of H\,{\sc i} clouds

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$α$ 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$α$ 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.

preprint2020arXiv

When Gas Dynamics Decouples from Galactic Rotation: Characterizing ISM Circulation in Disk Galaxies

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, $Γ$, 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 $λ_c$. We define the spatial scale $λ_{\rm eq}$ at which galactic rotation and non-circular motions contribute equally to $Γ$. 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 $λ_c$ and this transition is necessary to reproduce the circulation distribution. We find that $λ_{\rm eq}$, and therefore the structure of the gas velocity field, is set by the local conditions of gravitational stability and stellar feedback.

preprint2020arXiv

Where did the globular clusters of the Milky Way form? Insights from the E-MOSAICS simulations

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.

preprint2019arXiv

Mapping electron temperature variations across a spiral arm in NGC 1672

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.

preprint2019arXiv

The Gas-Star Formation Cycle in Nearby Star-Forming Galaxies I. Assessment of Multi-scale Variations

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.

preprint2019arXiv

The headlight cloud in NGC 628: An extreme giant molecular cloud in a typical galaxy disk

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.

preprint2019arXiv

The implications of clustered star formation for (proto)planetary systems and habitability

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.

preprint2019arXiv

The lifecycle of molecular clouds in nearby star-forming disc galaxies

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$α$ 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 $Σ_{\rm H_2}\geqslant8$M$_{\odot}$pc$^{-2}$, the GMC lifetime correlates with time-scales for galactic dynamical processes, whereas at $Σ_{\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$α$ (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.

preprint2019arXiv

The mass fraction of halo stars contributed by the disruption of globular clusters in the E-MOSAICS simulations

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.

preprint2018arXiv

Cloud-Scale Molecular Gas Properties in 15 Nearby Galaxies

We measure the velocity dispersion, $σ$, and surface density, $Σ$, 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. $Σ$ and $σ$ 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 $α_\mathrm{vir}\proptoσ^2/Σ$, 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Σ\,σ^2$, which spans ${\sim}5\rm\;dex$ across our sample. We find $Σ$, $σ$, 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 $σ$-$Σ$ relation flattens at $Σ\lesssim10\rm\;M_\odot\,pc^{-2}$, leading to high $σ$ for a given $Σ$ and high apparent $α_\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.

preprint2016arXiv

A Portrait of Cold Gas in Galaxies at 60pc Resolution and a Simple Method to Test Hypotheses That Link Small-Scale ISM Structure to Galaxy-Scale Processes

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.

preprint2016arXiv

Toward gas exhaustion in the W51 high-mass protoclusters

We present new JVLA observations of the high-mass cluster-forming region W51A from 2 to 16 GHz with resolution $θ_{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.

preprint2016arXiv

Using young massive star clusters to understand star formation and feedback in high-redshift-like environments

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.

preprint2015arXiv

A Dynamical Model for the Formation of Gas Rings and Episodic Starbursts Near Galactic Centres

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.

preprint2015arXiv

Deeply Embedded Protostellar Population in the 20 km s-1 Cloud of the Central Molecular Zone

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.

preprint2015arXiv

Dense gas in the Galactic central molecular zone is warm and heated by turbulence

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 $ζ_{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.

preprint2015arXiv

Globular cluster formation in the context of galaxy formation and evolution

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.

preprint2015arXiv

Globular clusters as the relics of regular star formation in 'normal' high-redshift galaxies

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.

preprint2015arXiv

High-mass star-forming cloud G0.38+0.04 in the Galactic Center Dust Ridge contains H2CO and SiO masers

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.

preprint2015arXiv

Pitfalls when observationally characterizing the relative formation rates of stars and stellar clusters in galaxies

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 ($Γ$) is a fundamental outcome of the star formation process. Recent observational and theoretical work has suggested that $Γ$ increases with the gas surface density ($Σ$) or star formation rate (SFR) surface density ($Σ_{\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 $Γ$ with the absolute SFR is much weaker than with $Σ_{\rm SFR}$, due to the mass-radius-SFR relation of star-forming 'main sequence' galaxies. The environmental variation of $Γ$ should therefore be probed as a function of area-normalised quantities, such as $Σ$ or $Σ_{\rm SFR}$. We present a set of guidelines for meaningful observational tests of cluster formation theories and show that these resolve the reported discrepancy.

preprint2014arXiv

An uncertainty principle for star formation. I. Why galactic star formation relations break down below a certain spatial scale

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.

preprint2014arXiv

Protoplanetary disc evolution affected by star-disc interactions in young stellar clusters

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}$.

preprint2014arXiv

The dynamical evolution of molecular clouds near the Galactic Centre - I. Orbital structure and evolutionary timeline

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)

preprint2014arXiv

The Formation and Early Evolution of Young Massive Clusters

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]

preprint2014arXiv

What controls star formation in the central 500 pc of the Galaxy?

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.

preprint2013arXiv

Color-Magnitude Relations within Globular Cluster Systems of Giant Elliptical Galaxies: The Effects of Globular Cluster Mass Loss and the Stellar Initial Mass Function

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 <~ α<~ -2.3$ in $dN/dM \propto M^α$), 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.

preprint2013arXiv

Comparing molecular gas across cosmic time-scales: the Milky Way as both a typical spiral galaxy and a high-redshift galaxy analogue

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.

preprint2013arXiv

Did massive black holes in globular clusters initially satisfy galactic scaling relations?

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.

preprint2013arXiv

N-body simulations of globular clusters in tidal fields: Effects of intermediate-mass black holes

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.

preprint2013arXiv

The Current State of Cluster Formation Simulations

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.

preprint2013arXiv

The Optical Colors of Giant Elliptical Galaxies and their Metal-Rich Globular Clusters Indicate a Bottom-Heavy Initial Mass Function

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.

preprint2012arXiv

A galactic-scale origin for stellar clustering

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.

preprint2012arXiv

Collateral Damage: the Implications of Utrecht Star Cluster Astrophysics for Galaxy 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.

preprint2012arXiv

High-velocity stars in the cores of globular clusters: The illustrative case of NGC 2808

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.

preprint2012arXiv

On the fraction of star formation occurring in bound stellar clusters

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.

preprint2011arXiv

Formation versus destruction: the evolution of the star cluster population in galaxy mergers

(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.

preprint2011arXiv

Modeling the formation and evolution of star cluster populations in galaxy simulations

(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.

preprint2011arXiv

The dynamical state of stellar structure in star-forming regions

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.

preprint2011arXiv

The dynamics of star-forming regions - which mechanisms set the cluster formation efficiency?

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.

preprint2011arXiv

The initial mass spectrum of old globular clusters in dwarf galaxies

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.

preprint2009arXiv

Dissolution is the solution: on the reduced mass-to-light ratios of Galactic globular clusters

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)

preprint2009arXiv

On the Interpretation of the Globular Cluster Luminosity Function

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.

preprint2009arXiv

The evolution of the stellar mass function in star clusters

(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.

J. M. Diederik Kruijssen

What is connected

Connect this record

See the researcher in context

Building this map preview

100 published item(s)

Unveiling the 3D structure of the central molecular zone from stellar kinematics and photometry: The 50 and 20 km/s clouds

Rotation and stability of the circumnuclear gas disk in the Galactic Center potential by the ALMA CMZ Exploration Survey (ACES)

In-situ or accreted? Using deep learning to infer the origin of extragalactic globular clusters from observables

A 2-3 mm high-resolution molecular line survey towards the centre of the nearby spiral galaxy NGC 6946

A CO isotopologue Line Atlas within the Whirlpool galaxy Survey (CLAWS)

A tale of two DIGs: The relative role of HII regions and low-mass hot evolved stars in powering the diffuse ionised gas (DIG) in PHANGS-MUSE galaxies

A trail of dark matter-free galaxies from a bullet dwarf collision

Constraining the formation of NGC1052-DF2 from its unusual globular cluster population

Constraining the shape of dark matter haloes with globular clusters

Empirically-motivated early feedback: momentum input by stellar feedback in galaxy simulations inferred through observations

Environmental dependence of the molecular cloud lifecycle in 54 main sequence galaxies

Introducing EMP-Pathfinder: modelling the simultaneous formation and evolution of stellar clusters in their host galaxies

Planetary Nebula Luminosity Function distances for 19 galaxies observed by PHANGS-MUSE

The emergence of dark matter-deficient ultra-diffuse galaxies driven by scatter in the stellar mass-halo mass relation and feedback from globular clusters

The Gas-Star Formation Cycle in Nearby Star-forming Galaxies II. Resolved Distributions of CO and H$α$ Emission for 49 PHANGS Galaxies

The initial conditions for young massive cluster formation in the Galactic Centre: convergence of large-scale gas flows

The PHANGS-MUSE survey -- Probing the chemo-dynamical evolution of disc galaxies

The present-day globular cluster kinematics of lenticular galaxies from the E-MOSAICS simulations and their relation to the galaxy assembly histories

Uncertainties in supernova input rates drive qualitative differences in simulations of galaxy evolution

WISDOM Project -- X. The morphology of the molecular ISM in galaxy centres and its dependence on galaxy structure

WISDOM project -- XI. Star Formation Efficiency in the Bulge of the AGN-host Galaxy NGC 3169 with SITELLE and ALMA

Age Determination of Galaxy Merger Remnant Stars using Asteroseismology

Giant Molecular Cloud Catalogues for PHANGS-ALMA: Methods and Initial Results

Globular clusters as tracers of the dark matter halo: insights from the E-MOSAICS simulations

Low-J CO Line Ratios From Single Dish CO Mapping Surveys and PHANGS-ALMA

Molecular Gas Properties and CO-to-H2 Conversion Factors in the Central Kiloparsec of NGC 3351

NGC5846-UDG1: A galaxy formed mostly by star formation in massive, extremely dense clumps of gas

PHANGS-HST: Star Cluster Spectral Energy Distribution Fitting with CIGALE

PHANGS-MUSE: the HII region luminosity function of local star-forming galaxies

Star formation in 'the Brick': ALMA reveals an active proto-cluster in the Galactic centre cloud G0.253+0.016

The PHANGS-HST Survey: Physics at High Angular resolution in Nearby GalaxieS with the Hubble Space Telescope

The physics governing the upper truncation mass of the globular cluster mass function

What to expect when using globular clusters as tracers of the total mass distribution in Milky Way-mass galaxies

A model for the onset of self-gravitation and star formation in molecular gas governed by galactic forces: II. the bottleneck to collapse set by cloud-environment decoupling

AlFoCS + Fornax3D: resolved star formation in the Fornax cluster with ALMA and MUSE

ALMA Observations of Massive Clouds in the Central Molecular Zone: Jeans Fragmentation and Cluster Formation

An uncertainty principle for star formation -- V. The influence of dust extinction on star formation rate tracer lifetimes and the inferred molecular cloud lifecycle

CMZoom: Survey Overview and First Data Release

Deep Transfer Learning for Star Cluster Classification: I. Application to the PHANGS-HST Survey

Distances to PHANGS Galaxies: New Tip of the Red Giant Branch Measurements and Adopted Distances

Dynamical Equilibrium in the Molecular ISM in 28 Nearby Star-Forming Galaxies

Entropy Driven Winds: Outflows and Fountains Lifted Gently by Buoyancy

From diffuse gas to dense molecular cloud cores

Heart of Darkness: the influence of galactic dynamics on quenching star formation in galaxy spheroids

Impact of low-energy cosmic rays on star formation

Not all stars form in clusters -- $Gaia$-DR2 uncovers the origin of OB associations

PHANGS CO kinematics: disk orientations and rotation curves at 150 pc resolution

Physical Processes in Star Formation

Star clusters near and far; tracing star formation across cosmic time

Stellar Feedback and Resolved Stellar IFU Spectroscopy in the nearby Spiral Galaxy NGC 300

Testing viscous disc theory using the balance between stellar accretion and external photoevaporation of protoplanetary discs

The CO-dark molecular gas mass in 30 Doradus

The globular cluster system mass-halo mass relation in the E-MOSAICS simulations

The molecular cloud lifecycle

The Origin of the Stellar Mass Distribution and Multiplicity

The Physics of Star Cluster Formation and Evolution

The role of galactic dynamics in shaping the physical properties of giant molecular clouds in Milky Way-like galaxies

Towards a multi-tracer timeline of star formation in the LMC -- I.\ Deriving the lifetimes of H\,{\sc i} clouds

When Gas Dynamics Decouples from Galactic Rotation: Characterizing ISM Circulation in Disk Galaxies

Where did the globular clusters of the Milky Way form? Insights from the E-MOSAICS simulations

Mapping electron temperature variations across a spiral arm in NGC 1672

The Gas-Star Formation Cycle in Nearby Star-Forming Galaxies I. Assessment of Multi-scale Variations

The headlight cloud in NGC 628: An extreme giant molecular cloud in a typical galaxy disk

The implications of clustered star formation for (proto)planetary systems and habitability

The lifecycle of molecular clouds in nearby star-forming disc galaxies

The mass fraction of halo stars contributed by the disruption of globular clusters in the E-MOSAICS simulations

Cloud-Scale Molecular Gas Properties in 15 Nearby Galaxies

A Portrait of Cold Gas in Galaxies at 60pc Resolution and a Simple Method to Test Hypotheses That Link Small-Scale ISM Structure to Galaxy-Scale Processes

Toward gas exhaustion in the W51 high-mass protoclusters

Using young massive star clusters to understand star formation and feedback in high-redshift-like environments

A Dynamical Model for the Formation of Gas Rings and Episodic Starbursts Near Galactic Centres

Deeply Embedded Protostellar Population in the 20 km s-1 Cloud of the Central Molecular Zone

Dense gas in the Galactic central molecular zone is warm and heated by turbulence

Globular cluster formation in the context of galaxy formation and evolution