Search

We present a comprehensive analysis of globular cluster (GC) formation and evolution across the $34^3$ Mpc$^3$ volume of the E-MOSAICS galaxy formation simulations. Defining GCs as surviving, high-mass ($>10^5$ M$_\odot$) clusters, we analyse their formation histories as a function of their metallicity and host galaxy mass, also distinguishing between central and satellite galaxies. The redshift of peak GC formation rate increases weakly with galaxy mass, decreases with metallicity, and does not differ between centrals and satellites. The epoch of peak GC formation precedes that of the stars by a factor of $1.1{-}1.6$, primarily due to `downsizing', i.e. low-mass galaxies form their stars later. Consequently, this offset decreases with galaxy mass, leading to nearly coeval stellar and GC populations in massive galaxies ($>10^{11}$ M$_\odot$). GCs themselves do not exhibit strong downsizing, because they predominantly formed at early cosmic epochs conducive to the formation (through high gas pressures) and survival (through high galaxy merger and GC migration rates) of massive, compact stellar systems. The total GC formation rate in the volume peaks at $z\approx 2.5$, shortly before star formation peaks at $z\approx 2$, but well after the general cluster formation rate at $z\approx 4$, reflecting a survivor bias where surviving GCs formed more recently. We find that GC formation commenced early, at $z>10$, such that the results of this work may provide a framework for interpreting direct observations of proto-GC formation with the JWST, especially as these observations accumulate to enable statistical studies., Comment: 12 pages, 7 figures, 3 tables; submitted to The Open Journal of Astrophysics (December 6, 2024); Figures 6 and 7 show the Madau plot for GCs

This study compiles stellar populations and internal properties of ultra-diffuse galaxies (UDGs) to highlight correlations with their local environment, globular cluster (GC) richness, and star formation histories. Complementing our sample of 88 UDGs, we include 36 low-surface brightness dwarf galaxies with UDG-like properties, referred to as NUDGes (nearly-UDGs). All galaxies were studied using the same spectral energy distribution fitting methodology to explore what sets UDGs apart from other galaxies. We show that NUDGes are similar to UDGs in all properties except for being, by definition, smaller and having higher surface brightness. We find that UDGs and NUDGes show similar behaviours in their GC populations, with the most metal-poor galaxies hosting consistently more GCs on average. This suggests that GC content may provide an effective way to distinguish extreme galaxies within the low surface brightness regime alongside traditional parameters like size and surface brightness. We confirm previous results using clustering algorithms that UDGs split into two main classes, which might be associated with the formation pathways of a puffy dwarf and a failed galaxy. The clustering applied to the UDGs+NUDGes dataset yields an equivalent result. The difference in mass contained in the GC system suggests that galaxies in different environments have not simply evolved from one another but may have formed through distinct processes., Comment: Accepted for publication in MNRAS. 12 pages, 7 figures

High-resolution imaging and strong gravitational lensing of high-redshift galaxies have enabled the detection of compact sources with properties similar to nearby massive star clusters. Often found to be very young, these sources may be globular clusters detected in their earliest stages. In this work, we compare predictions of high-redshift ($z \sim 1$--$10$) star cluster properties from the E-MOSAICS simulation of galaxy and star cluster formation with those of the star cluster candidates in strongly lensed galaxies from James Webb (JWST) and Hubble Space Telescope (HST) imaging. We select galaxies in the simulation that match the luminosities of the majority of lensed galaxies with star cluster candidates observed with JWST. We find that the luminosities, ages and masses of the brightest star cluster candidates in the high-redshift galaxies are consistent with the E-MOSAICS model. In particular, the brightest cluster ages are in excellent agreement. The results suggest that star clusters in both low- and high-redshift galaxies may form via common mechanisms. However, the brightest clusters in the lensed galaxies tend to be $\approx 1$--$1.5$ mag brighter and $\approx 0.5$ dex more massive than the median E-MOSAICS predictions. We discuss the large number of effects that could explain the discrepancy, including simulation and observational limitations, stellar population models, cluster detection biases and nuclear star clusters. Understanding these limitations would enable stronger tests of globular cluster formation models., Comment: 17 pages, 9 figures, 1 table. Accepted for publication in MNRAS

We present Hubble Space Telescope ACS/WFC and WFC3/UVIS imaging for a sample of 50 low surface brightness (LSB) galaxies in the $\sim$10$^{15}$ M$_{\odot}$ Perseus cluster, which were originally identified in ground-based imaging. We measure the structural properties of these galaxies and estimate the total number of globular clusters (GCs) they host. Around half of our sample galaxies meet the strict definition of an ultra-diffuse galaxy (UDG), while the others are UDG-like but are either somewhat more compact or slightly brighter. A small number of galaxies reveal systems with many tens of GCs, rivalling some of the richest GC systems known around UDGs in the Coma cluster. We find the sizes of rich GC systems, in terms of their half-number radii, extending to $\sim$1.2 times the half-light radii of their host galaxy on average. The mean colours of the GC systems are the same, within the uncertainties, as those of their host galaxy stars. This suggests that GCs and galaxy field stars may have formed at the same epoch from the same enriched gas. It may also indicate a significant contribution from disrupted GCs to the stellar component of the host galaxy as might be expected in the 'failed galaxy' formation scenario for UDGs., Comment: 16 pages, 11 figures, 2 tables, accepted for publication in MNRAS

Globular clusters (GCs) are sensitive tracers of galaxy assembly histories but interpreting the information they encode is challenging because mergers are thought to promote both the formation and disruption of GCs. We use simulations with controlled merger histories to examine the influence of merger mass ratio on the GC population of a present-day $L^\ast$ galaxy, using the genetic modification technique to adjust the initial conditions of a galaxy that experiences major mergers at $z = 1.7$ and $z = 0.77$ (ORGANIC case), so the later merger has twice its original mass ratio (ENHANCED case), or is prevented from occurring (SUPPRESSED case). We evolve the three realizations with E-MOSAICS, which couples sub-grid star cluster formation and evolution models to the EAGLE galaxy formation model. Relative to the ORGANIC case, the mass of surviving GCs is elevated (reduced) in the ENHANCED (SUPPRESSED) case, indicating that major mergers promote a net boost to the GC population. The boost is clearly quantified by the GC specific mass, $S_{\rm M}$, because it is sensitive to the number of the most massive GCs, whose long characteristic disruption timescales enable them to survive their hostile natal environments. In contrast, the specific frequency, $T_{\rm N}$, is insensitive to assembly history because it primarily traces low-mass GCs that tend to be disrupted soon after their formation. The promotion of GC formation and disruption by major mergers imprints a lasting and potentially observable signature: an elevated mass fraction of field stars in the galaxy's stellar halo that were born in star clusters., Comment: 17 pages, 8 figures. Submitted, comments welcome

Observational surveys have found that the dynamical masses of ultra-diffuse galaxies (UDGs) correlate with the richness of their globular cluster (GC) system. This could be explained if GC-rich galaxies formed in more massive dark matter haloes. We use simulations of galaxies and their GC systems from the E-MOSAICS project to test whether the simulations reproduce such a trend. We find that GC-rich simulated galaxies in galaxy groups have enclosed masses that are consistent with the dynamical masses of observed GC-rich UDGs. However, simulated GC-poor galaxies in galaxy groups have higher enclosed masses than those observed. We argue that GC-poor UDGs with low stellar velocity dispersions are discs observed nearly face on, such that their true mass is underestimated by observations. Using the simulations, we show that galactic star-formation conditions resulting in dispersion-supported stellar systems also leads to efficient GC formation. Conversely, conditions leading to rotationally-supported discs leads to inefficient GC formation. This result may explain why early-type galaxies typically have richer GC systems than late-type galaxies. This is also supported by comparisons of stellar axis ratios and GC specific frequencies in observed dwarf galaxy samples, which show GC-rich systems are consistent with being spheroidal, while GC-poor systems are consistent with being discs. Therefore, particularly for GC-poor galaxies, rotation should be included in dynamical mass measurements from stellar dynamics., Comment: 16 pages, 8 figures. Accepted for publication in MNRAS

We use spectral energy distribution (SED) fitting to place constraints on the stellar populations of 59 ultra-diffuse galaxies (UDGs) in the low-to-moderate density fields of the MATLAS survey. We use the routine PROSPECTOR, coupled with archival data in the optical from DECaLS, and near- and mid-infrared imaging from WISE, to recover the stellar masses, ages, metallicities and star formation timescales of the UDGs. We find that a subsample of the UDGs lies within the scatter of the mass-metallicity relation (MZR) for local classical dwarfs. However, another subsample is more metal-poor, being consistent with the evolving MZR at high-redshift. We investigate UDG positioning trends in the mass-metallicity plane as a function of surface brightness, effective radius, axis ratio, local volume density, mass-weighted age, star formation timescale, globular cluster (GC) counts and GC specific frequency. We find that our sample of UDGs can be separated into two main classes. Class A: Comprised of UDGs with lower stellar masses, prolonged star formation histories (SFHs), more elongated, inhabiting less dense environments, hosting fewer GCs, younger, consistent with the classical dwarf MZR, and fainter. Class B: UDGs with higher stellar masses, rapid SFHs, rounder, inhabiting the densest of our probed environments, hosting on average the most numerous GC systems, older, consistent with the high-redshift MZR (i.e., consistent with early-quenching), and brighter. The combination of these properties suggests that UDGs of Class A are consistent with a `puffed-up dwarf' formation scenario, while UDGs of Class B seem to be better explained by `failed galaxy' scenarios., Comment: Accepted for publication in MNRAS. 19 pages (+6 of appendices), 10 figures, 5 tables

The 21-cm spectral line widths, $w_{50}$, of galaxies are an approximate tracer of their dynamical masses, such that the dark matter halo mass function is imprinted in the number density of galaxies as a function of $w_{50}$. Correcting observed number counts for survey incompleteness at the level of accuracy needed to place competitive constraints on warm dark matter (WDM) cosmological models is very challenging, but forward-modelling the results of cosmological hydrodynamical galaxy formation simulations into observational data space is more straightforward. We take this approach to make predictions for an ALFALFA-like survey from simulations using the EAGLE galaxy formation model in both cold (CDM) and WDM cosmogonies. We find that for WDM cosmogonies more galaxies are detected at the low-$w_{50}$ end of the 21-cm velocity width function than in the CDM cosmogony, contrary to what might na\"ively be expected from the suppression of power on small scales in such models. This is because low-mass galaxies form later and retain more gas in WDM cosmogonies (with EAGLE). While some shortcomings in the treatment of cold gas in the EAGLE model preclude placing definitive constraints on WDM scenarios, our analysis illustrates that near-future simulations with more accurate modelling of cold gas will likely make strong constraints possible, especially in conjunction with new 21-cm surveys such as WALLABY., Comment: 12 pages, 7 figures. MNRAS accepted version

We study the statistics of dwarf galaxy populations as a function of environment in the cold dark matter (CDM) and warm dark matter (WDM) cosmogonies, using hydrodynamical simulations starting from initial conditions with matched phases but differing power spectra, and evolved with the EAGLE galaxy formation model. We measure the abundance of dwarf galaxies within 3~Mpc of DM haloes with a present-day halo mass similar to that of the Milky Way (MW), and find that the radial distribution of galaxies $M_{*}>10^7$\Msun is nearly identical for WDM and CDM. However, the cumulative mass function becomes shallower for WDM at lower masses, yielding 50~per~cent fewer dwarf galaxies of $M_{*}\gtrsim10^{5}$~\Msun than CDM. The suppression of low-mass halo counts in WDM relative to CDM increases significantly from high-density regions to low-density regions for haloes in the region of the half-mode mass, $M_\rm{hm}$. The luminous fraction in the two models also diverges from the overdense to the underdense regions for $M>2M_\rm{hm}$, as the increased collapse delay at small densities pushes the collapse to after the reionization threshold. However, the stellar mass--halo mass relation of WDM haloes relative to CDM increases towards lower-density regions. Finally, we conclude that the suppression of galaxies with $M_{*}\gtrsim10^5$\Msun between WDM and CDM is independent of density: the suppression of halo counts and the luminous fraction is balanced by an enhancement in stellar mass--halo mass relation., Comment: 10 pages, 6 figures, manuscript to be submitted to MNRAS

The observed scaling relations between supermassive black hole masses and their host galaxy properties indicate that supermassive black holes influence the evolution of galaxies. However, the scaling relations may be affected by selection biases. We propose to measure black hole masses in a mass-limited galaxy sample including all non-detections to inprove constraints on galaxy mass - black hole mass scaling relations and test for selection bias. We use high spatial resolution spectroscopy from the Keck and Gemini telescopes, and the Jeans Anisotropic Modelling method to measure black hole masses in early type galaxies from the Virgo Cluster. We present four new black hole masses and one upper limit in our mass-selected sample of galaxies of galaxy mass (1.0-3.2) x $10^{10} M_\odot$. This brings the total measured to 11 galaxies out of a full sample of 18 galaxies, allowing us to constrain scaling relations. We calculate a lower limit for the average black hole mass in our sample of 3.7 x $10^{7} M_\odot$. This is at an average galaxy stellar mass of (1.81 +/- 0.14) x $10^{10} M_\odot $ and an average bulge mass of (1.31 +/- 0.15) x $10^{10} M_\odot$. This lower limit shows that black hole masses in early type galaxies are not strongly affected by selection biases.

We examine the stellar [$\alpha$/Fe]-[Fe/H] distribution of $\simeq1000$ present-day galaxies in a high-resolution EAGLE simulation. A slight majority of galaxies exhibit the canonical distribution, characterised by a sequence of low-metallicity stars with high [$\alpha$/Fe] that transitions at a "knee" to a sequence of declining [$\alpha$/Fe] with increasing metallicity. This population yields a knee metallicity - galaxy-mass relation similar to that observed in Local Group galaxies, both in slope and scatter. However, many simulated galaxies lack a knee or exhibit more complicated distributions. Knees are found only in galaxies with star formation histories (SFHs) featuring a sustained decline from an early peak ($t\simeq7~{\rm Gyr}$), which enables enrichment by Type Ia supernovae to dominate that due to Type II supernovae (SN II), reducing [$\alpha$/Fe] in the interstellar gas. The simulation thus indicates that, contrary to the common interpretation implied by analytic galactic chemical evolution (GCE) models, knee formation is not a consequence of the onset of enrichment by SN Ia. We use the SFH of a simulated galaxy exhibiting a knee as input to the VICE GCE model, finding it yields an $\alpha$-rich plateau enriched only by SN II, but the plateau comprises little stellar mass and the galaxy forms few metal-poor ([Fe/H]$\lesssim - 1$) stars. This follows from the short, constant gas consumption timescale typically assumed by GCEs, which implies the presence of a readily-enriched, low-mass gas reservoir. When an initially longer, evolving consumption timescale is adopted, VICE reproduces the simulated galaxy's track through the [$\alpha$/Fe]-[Fe/H] plane and its metallicity distribution function., Comment: 17 pages, 12 figures, 1 table, submitted to MNRAS

JWST observations of the strongly lensed galaxy The Sparkler have revealed a population of gravitationally bound globular cluster (GC) candidates. Different analyses have resulted in broadly similar ages but significantly different metallicities, questioning the assembly history that has led to the formation of such a population. In this letter, we re-analyse the two sets of photometry available in the literature with the code MCMAME especially tailored to fit physical properties of GCs. We find the ages and metallicities from both datasets are consistent within 1 $\sigma$ uncertainties. A significant group of GCs is consistent with being old and metal poor ([Fe/H] $\sim -1.7$). For this group, the ages do not converge, hence, we conclude that they are definitively older than 1 Gyr and can be as old as the age of the Universe. The remaining GCs have younger ages and a metallicity spread. The ages and metallicities distribution of GCs in the Sparkler are consistent with those observed in Local Group's galaxies at similar lookback times. Comparing with predictions from E-MOSAICS simulations we confirm that the Sparkler GC population traces the self-enrichment history of a galaxy which might become a few times $10^9$ M$_{\odot}$ massive system at redshift $z = 0$, Comment: MNRAS letter, accepted for publication

We use the hydrodynamic EAGLE simulation to predict the numbers and masses of supermassive black holes in remnant nuclei of disrupted galaxies (stripped nuclei) and compare these to confirmed measurements of black holes in observed ultra-compact dwarf galaxies (UCDs). We find that black holes in stripped nuclei are consistent with the numbers and masses of those in observed UCDs. Approximately 50 per cent of stripped nuclei with $M > 2 \times 10^6 M_\odot$ should contain supermassive black holes. We further calculate how the presence of a black hole increases the dynamical mass of a stripped nucleus via the mass elevation ratio $\Psi$ defined as the ratio of the kinematically derived mass to the expected mass from stellar population synthesis. We find $\Psi_{sim} = 1.51^{+0.06}_{-0.04}$ for $M > 10^7 M_\odot$ stripped nuclei, consistent with that of observed UCDs which have $\Psi_{obs} = 1.7 \pm 0.2$ above $M > 10^7 M_\odot$. We also find that the mass elevation ratios of stripped nuclei with supermassive black holes can explain the observed number of UCDs with elevated mass-to-light ratios. Finally, we predict the relative number of massive black holes in stripped nuclei and galaxy nuclei and find that stripped nuclei should increase the number of black holes in galaxy clusters by 30 - 100 per cent, depending on the black hole occupation fraction of low-mass galaxies. We conclude that the population of supermassive black holes in UCDs represents a large and unaccounted-for portion of supermassive black holes in galaxy clusters., Comment: 15 pages, revised after submission to MNRAS

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., Comment: 21 pages, 22 figures. Submitted to MNRAS. Comments welcome

The optical morphology of galaxies is strongly related to galactic environment, with the fraction of early-type galaxies increasing with local galaxy density. In this work we present the first analysis of the galaxy morphology-density relation in a cosmological hydrodynamical simulation. We use a convolutional neural network, trained on observed galaxies, to perform visual morphological classification of galaxies with stellar masses $M_\ast > 10^{10} \, \mathrm{M}_\odot$ in the EAGLE simulation into elliptical, lenticular and late-type (spiral/irregular) classes. We find that EAGLE reproduces both the galaxy morphology-density and morphology-mass relations. Using the simulations, we find three key processes that result in the observed morphology-density relation: (i) transformation of disc-dominated galaxies from late-type (spiral) to lenticular galaxies through gas stripping in high-density environments, (ii) formation of lenticular galaxies by merger-induced black hole feedback in low-density environments, and (iii) an increasing fraction of high-mass galaxies, which are more often elliptical galaxies, at higher galactic densities., Comment: 19 pages, 17 figures, published in MNRAS

The metallicity distributions of globular cluster (GC) systems in galaxies are a critical test of any GC formation scenario. In this work, we investigate the predicted GC metallicity distributions of galaxies in the MOdelling Star cluster population Assembly In Cosmological Simulations within EAGLE (E-MOSAICS) simulation of a representative cosmological volume ($L = 34.4$ comoving Mpc). We find that the predicted GC metallicity distributions and median metallicities from the fiducial E-MOSAICS GC formation model agree well the observed distributions, except for galaxies with masses $M_\ast \sim 2 \times 10^{10}$ M$_\odot$, which contain an overabundance of metal-rich GCs. The predicted fraction of galaxies with bimodal GC metallicity distributions ($37 \pm 2$ per cent in total; $45 \pm 7$ per cent for $M_\ast > 10^{10.5}$ M$_\odot$) is in good agreement with observed fractions ($44^{+10}_{-9}$ per cent), as are the mean metallicities of the metal-poor and metal-rich peaks. We show that, for massive galaxies ($M_\ast > 10^{10}$ M$_\odot$), bimodal GC distributions primarily occur as a result of cluster disruption from initially-unimodal distributions, rather than as a result of cluster formation processes. Based on the distribution of field stars with GC-like abundances in the Milky Way, we suggest that the bimodal GC metallicity distribution of Milky Way GCs also occurred as a result of cluster disruption, rather than formation processes. We conclude that separate formation processes are not required to explain metal-poor and metal-rich GCs, and that GCs can be considered as the surviving analogues of young massive star clusters that are readily observed to form in the local Universe today., Comment: 19 pages, 16 figures. Published in MNRAS

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 ($\sigma\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., Comment: 17 pages, 15 figures; submitted to MNRAS (25/04/2022) - comments and suggestions welcome!

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 ($\tau>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., Comment: 30 pages + appendices, 12 figures. MNRAS submitted, comments and suggestions welcome!