Globular Cluster System Research Articles

Estimating galaxy masses from kinematics of globular cluster systems: a new method based on deep learning

ABSTRACT We present a new method by which the total masses of galaxies including dark matter can be estimated from the kinematics of their globular cluster systems (GCSs). In the proposed method, we apply the convolutional neural networks (CNNs) to the 2D maps of line-of-sight velocities (V) and velocity dispersions (σ) of GCSs predicted from numerical simulations of disc and elliptical galaxies. In this method, we first train the CNN using either only a larger number ($\sim 200\, 000$) of the synthesized 2D maps of σ (‘one-channel’) or those of both σ and V (‘two-channel’). Then, we use the CNN to predict the total masses of galaxies (i.e. test the CNN) for the totally unknown data set that is not used in training the CNN. The principal results show that overall accuracy for one-channel and two-channel data is 97.6 per cent and 97.8 per cent, respectively, which suggests that the new method is promising. The mean absolute errors (MAEs) for one-channel and two-channel data are 0.288 and 0.275, respectively, and the value of root mean square errors (RMSEs) are 0.539 and 0.51 for one-channel and two-channel, respectively. These smaller MAEs and RMSEs for two-channel data (i.e. better performance) suggest that the new method can properly consider the global rotation of GCSs in the mass estimation. We also applied our proposed method to real data collected from observations of NGC 3115 to compare the total mass predicted by our proposed method and other popular methods from the literature.

Eight more low luminosity globular clusters in the Sagittarius dwarf galaxy

Context. The Sagittarius (Sgr) dwarf galaxy is merging with the Milky Way, and the study of its globular clusters (GCs) is important to understand the history and outcome of this ongoing process. Aims. Our main goal is to characterize the GC system of the Sgr dwarf galaxy. This task is hampered by high foreground stellar contamination, mostly from the Galactic bulge. Methods. We performed a GC search specifically tailored to find new GC members within the main body of this dwarf galaxy using the combined data of the VISTA Variables in the Via Lactea Extended Survey near-infrared survey and the Gaia Early Data Release 3 optical database. Results. We applied proper motion cuts to discard foreground bulge and disk stars, and we found a number of GC candidates in the main body of the Sgr dwarf galaxy. We selected the best GCs as those objects that have significant overdensities above the stellar background of the Sgr galaxy and that possess color-magnitude diagrams with well-defined red giant branches consistent with the distance and reddening of this galaxy. Conclusions. We discover eight new GC members of the Sgr galaxy, which adds up to 29 total GCs known in this dwarf galaxy. This total number of GCs shows that the Sgr dwarf galaxy hosts a rather rich GC system. Most of the new GCs appear to be predominantly metal-rich and have low luminosity. In addition, we identify ten other GC candidates that are more uncertain and need more data for proper confirmation.

Dynamical evidence for a morphology-dependent relation between the stellar and halo masses of galaxies

We derive the stellar-to-halo mass relation (SHMR), namely f⋆ ∝ M⋆/Mh versus M⋆ and Mh, for early-type galaxies from their near-infrared luminosities (for M⋆) and the position-velocity distributions of their globular cluster systems (for Mh). Our individual estimates of Mh are based on fitting a flexible dynamical model with a distribution function expressed in terms of action-angle variables and imposing a prior on Mh from the correlation between halo concentration and mass in the standard Λ cold dark matter (ΛCDM) cosmology. We find that the SHMR for early-type galaxies declines with mass beyond a peak at M⋆ ∼ 5 × 1010 M⊙ and Mh ∼ 1 × 1012 M⊙ (near the mass of the Milky Way). This result is consistent with the standard SHMR derived by abundance matching for the general population of galaxies, and also with previous, less robust derivations of the SHMR for early-type galaxies. However, it contrasts sharply with the monotonically rising SHMR for late-type galaxies derived from extended HI rotation curves and the same ΛCDM prior on Mh that we adopt for early-type galaxies. We show that the SHMR for massive galaxies varies more or less continuously with disc fraction and Hubble type between these rising and falling branches. We also show that the different SHMRs for late-type and early-type galaxies are consistent with the similar scaling relations between their stellar velocities and masses (the Tully–Fisher and the Faber–Jackson relations). As we demonstrate explicitly, differences in the relations between the stellar and halo virial velocities account for the similarity of the scaling relations. We argue that all these empirical findings are natural consequences of a picture in which galactic discs are built mainly by relatively smooth and gradual inflow, regulated by feedback from young stars, while galactic spheroids are built by a combination of merging, black-hole fuelling, and feedback from active galactic nuclei.

Confirmation and physical characterization of the new bulge globular cluster Patchick 99 from the VVV and Gaia surveys

Context. Globular clusters (GCs) are recognised as important tools for understanding the formation and evolution of the Milky Way (MW) because they are the oldest objects in our Galaxy. Unfortunately, the known sample in our MW is still incomplete, especially towards the innermost regions, because of the high differential reddening, extinction, and stellar crowding. Therefore, the discovery of new GC candidates and the confirmation of their true nature are crucial for the census of the MW GC system. Aims. Our main goal is to confirm the physical nature of two GC candidates: Patchick 99 and TBJ 3. Both are located towards the Galactic bulge. We use public data in the near-infrared(IR) passband from the VISTA Variables in the Via Láctea Survey (VVV), the VVV eXtended Survey, and the Two Micron All Sky Survey along with deep optical data from the Gaia Mission DR 2 in order to estimate their main astrophysical parameters, such as reddening and extinction, distance, total luminosity, mean cluster proper motions, size, metallicity, and age. Methods. We investigated both candidates at different wavelengths, allowing us to discard TBJ3 as a possible GC. We use near-IR (Ks vs. (J − Ks)) and optical (G vs. (BP − RP)) colour–magnitude diagrams (CMDs) in order to analyse Patchick 99. First, we decontaminated CMDs, following a statistical procedure, as well as selecting only stars which have similar proper motions (PMs) and are situated within 3′ of the centre. Mean PMs were measured from Gaia DR 2 data. Reddening and extinction were derived by adopting optical and near-IR reddening maps, and were used to estimate the distance modulus and the heliocentric distance. Metallicity and age were evaluated by fitting theoretical stellar isochrones. Results. Reddening and extinction values for Patchick 99 are E(J − Ks) = (0.12 ± 0.02) mag and AKs = (0.09 ± 0.01) mag from the VVV data, whereas we calculate E(BP − RP) = (0.21 ± 0.03) mag and AG = (0.68 ± 0.08) mag from Gaia DR 2 data. We use those values and the magnitude of the RC to estimate the distance, finding good agreement between the near-IR and optical measurements. In fact, we obtain (m − M)0 = (14.02 ± 0.01) mag, equivalent to a distance D = (6.4 ± 0.2) kpc in near-IR and (m − M)0 = (14.23 ± 0.1) mag and so D = (7.0 ± 0.2) kpc in optical. In addition, we derive the metallicity and age for Patchick 99 using our distance and extinction values and fitting PARSEC isochrones. We find [Fe/H]=(−0.2 ± 0.2) dex and t = (10 ± 2) Gyr. The mean PMs for Patchick 99 are μα = ( − 2.98 ± 1.74) mas yr−1 and μδ = ( − 5.49 ± 2.02) mas yr−1 using the Gaia DR 2 data. These are consistent with the bulge kinematics. We also calculate the total luminosity of our cluster and confirm that it is a low-luminosity GC, with MKs = ( − 7.0 ± 0.6) mag. The radius estimation is performed building the radial density profile and we find its angular radius rP99 ∼ 10′. We also recognise seven RR Lyrae star members within 8.2 arcmin from the Patchick 99 centre, but only three of them have PMs matching the mean GC PM, confirming the distance found by other methods. Conclusions. We find that TBJ 3 shows mid-IR emissions that are not present in GCs. We therefore discard TBJ 3 as a GC candidate and focus our work on Patchick 99. We conclude that Patchick 99 is an old metal-rich GC situated in the Galactic bulge. TBJ 3 is a background galaxy.

A 2-d orbital period for a redback millisecond pulsar candidate in the globular cluster NGC 6397

ABSTRACT We report optical modulation of the companion to the X-ray source U18 in the globular cluster NGC 6397. U18, with combined evidence from radio and X-ray measurements, is a strong candidate as the second redback in this cluster, initially missed in pulsar searches. This object is a bright variable star with an anomalous red colour and optical variability (∼0.2 mag in amplitude) with a periodicity ∼1.96 d that can be interpreted as the orbital period. This value corresponds to the longest orbital period for known redback candidates and confirmed systems in Galactic globular clusters and one of the few with a period longer than 1 d.

Globular cluster systems of relic galaxies

ABSTRACT We analyse the globular cluster (GC) systems of a sample of 15 massive, compact early-type galaxies (ETGs), 13 of which have already been identified as good relic galaxy candidates on the basis of their compact morphologies, old stellar populations and stellar kinematics. These relic galaxy candidates are likely the nearby counterparts of high-redshift red nugget galaxies. Using F814W (≈I) and F160W (≈H) data from the WFC3 camara onboard the Hubble Space Telescope, we determine the total number, luminosity function, specific frequency, colour, and spatial distribution of the GC systems (GCSs). We find lower specific frequencies (SN < 2.5 with a median of SN = 1) than ETGs of comparable mass. This is consistent with a scenario of rapid, early dissipative formation, with relatively low levels of accretion of low-mass, high-SN satellites. The GC half-number radii are compact, but follow the relations found in normal ETGs. We identify an anticorrelation between the specific angular momentum (λR) of the host galaxy and the (I − H) colour distribution width of their GCSs. Assuming that λR provides a measure of the degree of dissipation in massive ETGs, we suggest that the (I − H) colour distribution width can be used as a proxy for the degree of complexity of the accretion histories in these systems.

The likelihood of undiscovered globular clusters in the outskirts of the Milky Way

ABSTRACT The currently known Galactic globular cluster population extends out to a maximum galactocentric distance of ∼145 kpc, with the peculiarity that the outermost clusters predominantly have an inward velocity. Orbit averaging finds that this configuration occurs by chance about $6{{\ \rm per\ cent}}$ of the time, suggesting that several globular clusters with positive radial velocities remain undiscovered. We evaluate the expected number of undiscovered clusters at large distances under the assumption that the cluster population has a smooth radial distribution and is in equilibrium within the Milky Way’s virial radius. By comparing the present day kinematic properties of outer clusters to random orbital configurations of the Galactic globular cluster system through orbit averaging, we estimate a likelihood of $73{{\ \rm per\ cent}}$ of there being at least one undiscovered globular cluster within the Milky Way. This estimate assumes the current population is complete out to 50 kpc, and increases to $91{{\ \rm per\ cent}}$ if the population is complete out to 150 kpc. The likelihood of there being two undiscovered clusters is between $60$ and $70{{\ \rm per\ cent}}$, with the likelihood of there being three undiscovered clusters being on the order of $50{{\ \rm per\ cent}}$. The most likely scenario is that the undiscovered clusters are moving outwards, which results in the outer cluster population being consistent with an equilibrium state. Searches for distant and possibly quite low concentration and very low metallicity globular clusters will be enabled with upcoming deep imaging surveys.

Mass Segregation as a New Indicator of Binary, IMBH, and Stellar-mass Black Hole Systems in Globular Clusters

Abstract This study presents the mass segregation of 35 Galactic globular clusters (GCs) through the radial profile of the mean mass of main-sequence stars. Most of them are found to be mass segregated due to two-body relaxation effects. The GCs with a larger core to half-mass–radius ratio tend to be less mass segregated. Recent simulations show that energy sources can quench mass segregation in relaxed GCs. Combining mass segregation and core binary fractions (f c ), we find a moderate anticorrelation between them when f c is below 0.15. Four GCs that have a small amount of mass segregation are selected as possible intermediate-mass black hole holders. However, none of them are strong candidates considering their half-mass relaxation times.

The Milky Way globular cluster system within the context of a closed box

Metallicities of the Milky Way globular cluster (GC) system are interpreted within the framework of a closed box model of chemical evolution. Model age-metallicity relations normalised to the Sun are compared with an observed age-metallicity relation constructed from literature data for the Galactic globular clusters and solar neighbourhood dwarf stars. One branch of the GC age-metallicity relation is seen to form a natural extension of the solar neighbourhood relation to low metallicities. Such clusters may fit in situ into the same Galactic closed-box model as many disk stars. Timescales for halo and disk GC formation are computed based on this assumption. However, a single closed box cannot account for a second branch in the GC age-metallicity data. This problem is addressed by assuming that some GCs did not form in situ within the closed box of the Galaxy, but rather were formed within separate parent objects whose chemical evolution was characterised by longer star formation timescales than the Milky Way box. These parent objects were eventually acquired by the Milky Way. Estimates are made of the initial mass and star formation timescale of those objects that contributed GCs to the accreted component of the Milky Way.

Nonlinear Color–Metallicity Relations of Globular Clusters. IX. Different Radial Number Density Profiles between Blue and Red Clusters

Abstract The optical colors of globular clusters (GCs) in most large early-type galaxies are bimodal. Blue and red GCs show a sharp difference in the radial profile of their surface number density in the sense that red GCs are more centrally concentrated than blue GCs. An instant interpretation is that there exist two distinct GC subsystems having different radial distributions. This view, however, was challenged by a scenario in which, due to the nonlinear nature of the GC metallicity-to-color transformation for old (≳10 Gyr) GCs, a broad unimodal metallicity spread can exhibit a bimodal color distribution. Here we show, by simulating the radial trends in the GC color distributions of the four nearby giant elliptical galaxies (M87, M49, M60, and NGC 1399), that the difference in the radial profile between blue and red GCs stems naturally from the metallicity-to-color nonlinearity plus the well-known radial metallicity gradient of GC systems. The model suggests no or little radial variation in GC age even out to ∼20 R eff. Our results provide a simpler solution to the distinct radial profiles of blue and red GCs that does not necessarily invoke the presence of two GC subsystems and further fortify the nonlinearity scenario for the GC color bimodality phenomenon.

Mysterious Globular Cluster System of the Peculiar Massive Galaxy M85

Abstract We present a study on the stellar population and kinematics of globular clusters (GCs) in the peculiar galaxy M85. We obtain optical spectra of 89 GCs at 8 kpc < R < 160 kpc using the MMT/Hectospec. We divide them into three groups, blue/green/red GCs (B/G/RGCs), with their (g − i)0 colors. All GC subpopulations have mean ages of about 10 Gyr, but showing differences in metallicities. The BGCs and RGCs are the most metal-poor ([Z/H] ∼ −1.49) and metal-rich ([Z/H] ∼ −0.45), respectively, and the GGCs are in between. We find that the inner GC system exhibits a strong overall rotation that is entirely due to a disklike rotation of the RGC system. The BGC system shows little rotation. The GGCs show kinematic properties clearly distinct among the GC subpopulations, having higher mean velocities than the BGCs and RGCs and being aligned along the major axis of M85. This implies that the GGCs have an origin different from the other GC subpopulations. The rotation-corrected velocity dispersion of the RGC system is much lower than that of the BGC system, indicating the truncation of the red halo of M85. The BGCs have a flat velocity dispersion profile out to R = 67 kpc, reflecting the dark matter extent of M85. Using the velocity dispersion of the BGC system, we estimate the dynamical mass of M85 to be 3.8 × 1012 M ⊙. We infer that M85 has undergone merging events lately, resulting in the peculiar kinematics of the GC system.

On the Evolution of the Globular Cluster System in NGC 1052-DF2: Dynamical Friction, Globular–Globular Interactions, and Galactic Tides

Abstract The ultradiffuse galaxy NGC 1052-DF2 has an overabundance of luminous globular clusters (GCs), and its kinematics is consistent with the presence of little to no dark matter. As the velocity dispersion among the GCs is comparable to the expected internal dispersions of the individual GCs, the galaxy might be highly conducive to GC–GC merging. If true, this could explain the puzzling luminosity function of its GCs. Here we examine this possibility by resimulating three of our earlier simulations of the GC system, where the GCs were modeled as single particles, with live GCs. Somewhat surprisingly, we infer a low merger rate of ∼0.03 Gyr−1. The main reason is that the GCs are too dense for tidal shock capture, caused by impulsive encounters among them, to operate efficiently (we infer a tidal capture rate of only ∼0.002 Gyr−1). Therefore, whatever mergers occur are driven by other mechanisms, which we find to be captures induced by dynamical friction and compressive tides from other GCs. The low merger rate inferred here makes it unlikely that the unusually large luminosities of the GCs can be explained as a result of past GC–GC mergers. Our simulations also indicate that, if NGC 1052-DF2 is indeed largely devoid of dark matter, its tidal field is too weak to induce any significant mass loss from the GCs. Therefore, in such a scenario, we predict that it is improbable for the GCs to reveal tidal features, something that can be tested with future deep observations.

Predicting accreted satellite galaxy masses and accretion redshifts based on globular cluster orbits in the E-MOSAICS simulations

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

The number of globular clusters around the iconic UDG DF44 is as expected for dwarf galaxies

ABSTRACT There is a growing consensus that the vast majority of ultradiffuse galaxies (UDGs) are dwarf galaxies. However, there remain a few UDGs that seem to be special in terms of their globular cluster (GC) systems. In particular, according to some authors, certain UDGs exhibit large GC populations when compared to expectations from their stellar (or total) mass. Among these special UDGs, DF44 in the Coma cluster is one of the better-known examples. DF44 has been claimed to have a relatively high number of GCs, $N_{\mathrm{ GC}}=74^{+18}_{-18}$, for a stellar mass of only $3\times 10^8\, \mathrm{ M}_{ \odot}$ which would indicate a much larger dark halo mass than dwarfs of similar stellar mass. In this paper, we revisit this number and, contrary to previous results, find $N_{\mathrm{ GC}}=21^{+7}_{-9}$ assuming that the distribution of the GCs follows the same geometry as the galaxy. If we assume that the GCs around DF44 are distributed in a (projected) circularly symmetric way and, if we use a less strict criterion for the selection of the GCs, we find $N_{\mathrm{ GC}}=18^{+23}_{-12}$. Making use of the MGC–Mhalo relation, this number of GCs suggests a dark matter halo mass of $M_{\mathrm{ halo}}=1.1^{+0.4}_{-0.5} \times 10^{11} \mathrm{ M}_{\odot}$, a value which is consistent with the expected total mass for DF44 based on its velocity dispersion, $\sigma =33^{+3}_{-3}$ km s−1. We conclude that the number of GCs around DF44 is as expected for regular dwarf galaxies of similar stellar mass and DF44 is not extraordinary in this respect.

The first detection of ultra-diffuse galaxies in the Hydra I cluster from the VEGAS survey

In this paper, we report the discovery of 27 low-surface brightness galaxies, of which 12 are candidates for ultra-diffuse galaxies (UDG) in the Hydra I cluster, based on deep observations taken as part of the VST Early-type Galaxy Survey (VEGAS). This first sample of UDG candidates in the Hydra I cluster represents an important step in our project that aims to enlarge the number of confirmed UDGs and, through study of statistically relevant samples, constrain the nature and formation of UDGs. This study presents the main properties of this class of galaxies in the Hydra I cluster. For all UDGs, we analysed the light and colour distribution, and we provide a census of the globular cluster (GC) systems around them. Given the limitations of a reliable GC selection based on two relatively close optical bands only, we find that half of the UDG candidates have a total GC population consistent with zero. Of the other half, two galaxies have a total population larger than zero at 2σ level. We estimate the stellar mass, the total number of GCs, and the GC specific frequency (SN). Most of the candidates span a range of stellar masses of 107 − 108 M⊙. Based on the GC population of these newly discovered UDGs, we conclude that most of these galaxies have a standard or low dark matter content, with a halo mass of ≤1010 M⊙.

Wide-field study of the peculiar globular cluster system hosted by the field lenticular NGC 1172

ABSTRACT We present a wide-field study of the globular cluster system (GCS) of the field lenticular galaxy NGC 1172, based on observations from Gemini Multi-Object Spectrograph/Gemini (optical), FourStar/Magellan (NIR), and archival data from ACS/Hubble Space Telescope (optical). This analysis covers the full extension of the GCS, and results in a value of specific frequency (SN = 8.6 ± 1.5) peculiarly high for an intermediate-mass galaxy in a low-density environment such as this one. We find that the GCS appears to be bimodal, although the colour distribution is narrow and does not allow for an accurate separation of the subpopulations. However, the combination of optical and NIR filters allows us to obtain an estimation of the metallicity distribution based on the photometry, which supports bimodality. We conclude that the presence of a large fraction of metal-poor globular clusters (GCs) and the high specific frequency point to NGC 1172 having accreted a significant amount of GCs from low-mass satellites in the past.

The Assembly History of M87 through Radial Variations in Chemical Abundances of Its Field Star and Globular Cluster Populations

Abstract We present an extensive study of spectroscopically derived chemical abundances for M87 and its globular cluster (GC) system. Using observations from the Mitchell spectrograph at McDonald, LRIS at Keck, and Hectospec on the MMT, we derive new metallicity gradients from ∼2 to 140 kpc. We use a novel hierarchical statistical framework to simultaneously separate the GC system into subpopulations while measuring the metallicity gradients of those subpopulations. We create physically motivated spectral stacks of the GC subpopulations by leveraging the output of this statistical framework to perform the first application of abundance tagging in a massive early-type galaxy to better constrain the origins of the GC subpopulations and thus the assembly history of M87. We find a metal-poor, α-enhanced population of GCs in both the inner and outer halos unanticipated by current cosmological simulations of galaxy evolution. We use the remarkably flat metallicity gradients we find for both the metal-rich and metal-poor GC subpopulations in the inner halo as tentative evidence that some amount of the metal-poor GCs formed directly in the halo of M87 at high redshift.

Substructure in the Globular Cluster Populations of the Virgo Cluster Elliptical Galaxies M84 and M86

Abstract We have carried out a search for substructure within the globular cluster (GC) systems of M84 (NGC 4374) and M86 (NGC 4406), two giant elliptical galaxies in the Virgo Cluster. We use wide-field (36′ × 36′), multicolor broadband imaging to identify GC candidates in these two galaxies, as well as several other nearby lower-mass galaxies. Our analysis of the spatial locations of the GC candidates reveals several substructures, including a peak in the projected number density of GCs in M86 that is offset from the system center and may be at least partly due to the presence of the dwarf elliptical galaxy NGC 4406B, a bridge that connects the M84 and M86 GC systems, and a boxy isodensity contour along the southeast side of the M86 GC system. We divide our sample into red (metal-rich) and blue (metal-poor) GC candidates to look for differences in the spatial distributions of the two populations and find that the blue cluster candidates are the dominant population in each of the substructures we identify. We also incorporate the measurements from two radial velocity surveys of the GCs in the region and find that the bridge substructure is populated by GCs with a mix of velocities that are consistent with either M86 and M84, possibly providing further evidence for interaction signatures between the two galaxies.

К ТЕОРИИ ПРОИСХОЖДЕНИЯ СИСТЕМ ШАРОВЫХ СКОПЛЕНИЙ ВОКРУГ ГАЛАКТИК

The origin of poor and moderate globular cluster systems around galaxies against the background of a nonlinear non-stationary model of collapsing galaxies is studied. The stability of the model with respect to four perturbation modes is studied, and the degree of the mode determines on average the number of clusters in the system. A critical dependence of the initial virial relation on the degree of rotation is constructed for each mode. The dependences of the instability increments on the physical parameters of the model are found.

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

ABSTRACT 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 × 1011 M⊙, 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–M⋆ 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 $\hbox{$M_{\rm GC}$}/ \hbox{$M_{\rm halo}$}\sim 5\times 10^{-5}$ and find that it is the result of a combination of cluster formation physics and cluster disruption.