A dynamical study of Galactic globular clusters under different relaxation conditions

It has long been argued that the globular clusters (GCs) in the Fornax dwarf galaxy indicate that its dark matter halo is likely to have a shallow density profile with a core of size ∼1 kpc. We revisit this argument by investigating analogues of Fornax formed in MOdelling Star cluster population Assembly In Cosmological Simulations within eagle (E-MOSAICS), a cosmological hydrodynamical simulation that follows the formation and evolution of GCs in the Evolution and Assembly of GaLaxies and their Environments (EAGLE) galaxy formation model. In eagle, Fornax-mass haloes are cuspy and well described by the Navarro–Frenk–White profile. We post-process the E-MOSAICS to account for GC orbital decay by dynamical friction, which is not included in the original model. Dynamical friction causes 33 per cent of GCs with masses $M_{\rm GC}\ge 4\times 10^4{~\rm M_\odot }$ to sink to the centre of their host with the majority being tidally disrupted before forming a nuclear star cluster. Fornax has a total of five GCs, an exceptionally large number compared to other galaxies of similar stellar mass. In the simulations, we find that only 3 per cent of the Fornax analogues have five or more GCs, while 30 per cent have only one and 35 per cent have none. We find that GC systems in satellites are more centrally concentrated than in field dwarfs, and that those formed in situ (45 per cent) are more concentrated than those that were accreted. The present-day radial distribution of GCs in E-MOSAICS Fornax analogues is indistinguishable from that in Fornax, demonstrating that the presence of five GCs in the central kiloparsec of Fornax is consistent with a cuspy dark matter halo.

Galaxy clusters represent one of the largest-scale structures in the Universe, with their environments of varying density being ideal for inducing an increase in activity of the cluster galaxy members and the pervading Intracluster Medium (ICM) between them.<br> However the precise drivers, thereby the inferred consequences, of differing states of galaxy cluster activity have not been adequately disentangled.<br> If one can procure a sufficient dichotomic sample of cluster galaxies of differing sub-populations from relaxed and unrelaxed galaxy cluster dynamical states, one could analyse how these sub-populations `move' as a function of both radius and dynamical state.<br> Specifically, using velocity dispersion profiles (VDPs) that weight against cluster galaxies across all radii with a Gaussian window function, one can determine how cluster galaxies of properties that pertain to galaxy evolutionary indicators respond to the cluster environment as a function of the clustocentric radius. <br> In addition, analysing the rotational profiles of these sub-populations as a function of the clustocentric radius can help identify where the bulk of the kinetic energy lies.<br> In view of these points, we have pooled data of galaxies from the Sloan Digital Sky Survey (SDSS) two composite galaxy clusters are assembled into unrelaxed and relaxed states, which are compared with one another through splitting their cluster galaxy populations by their evolutionary proxies (i.e. morphology, colour, mass and AGN), thus analysing their resultant VDPs and rotational profiles.<br> It is therefore found that unrelaxed galaxy clusters possess increased interaction between cluster galaxies as the radius increases, driven by blue and red elliptical galaxies, suggesting an induction of activity by pre-processed groups.<br> AGN-hosting cluster galaxy sub-populations within these unrelaxed environs however present increased activity as \(R\rightarrow0\), indicating these are recent infallers, potentially triggered due to interactions with the ICM.<br> Rotational profiles of unrelaxed clusters build upon the VDPs, showing significant bulk rotation of cluster galaxies with prominent core rotation driven by redder sub-populations.

We study the evolution of globular clusters with mass spectra under the influence of the steady Galactic tidal field, including the effects of velocity anisotropy. Similarly to single-mass models, velocity anisotropy develops as the cluster evolves, but the degree of anisotropy is much smaller than in isolated clusters. Except for very early epochs of the cluster evolution, the velocity distributions of nearly all mass components become tangentially anisotropic at the outer parts. We examine how the mass function (MF) changes in time. Specifically, we find that the power-law index of the MF decreases monotonically with the total mass of the cluster, in agreement with previous findings based on isotropic models or N-body studies. This is also consistent with the behaviour of the observed slopes of MFs for a limited number of clusters. We attempt to compare our results with multimass King models, although it is almost impossible to fit the entire density profiles for all mass components. When the MF is fixed, the central densities of individual components show significant differences between Fokker-Planck and King models. We obtain ‘best-fitting’ multimass King models, for which the central density of individual components as well as the total density distribution agrees with the Fokker-Planck models by adjusting the MF. The MFs obtained in this way closely resemble the MF within the half-mass radius of the Fokker-Planck result. Also, we find that the local MFs predicted by Fokker-Planck calculations vary more rapidly with radius than best-fitting multimass King models. The projected velocity profiles for anisotropic models show significant flattening toward the tidal radius compared with the isotropic model. This is caused by the fact that the tangential velocity dispersion becomes dominant at the outer parts. Such a behaviour of velocity profile appears to be consistent with the observed profiles of the collapsed cluster M15.

Nearly all globular clusters (GCs) studied to date show evidence for multiple stellar populations, in stark contrast to the conventional view that GCs are a mono-metallic, coeval population of stars. Building on earlier work, we propose a simple physical model for the early evolution (several 10^8 yr) of GCs. We consider the effects of stellar mass-loss, type II and prompt type Ia supernovae, ram pressure, and accretion from the ambient ISM on the development of a young GC's own gas reservoir. In our model, type II SNe from a first generation of star formation clears the GC of its initial gas reservoir. Over the next several 10^8 yr, mass lost from AGB stars and matter accreted from the ambient ISM collect at the center of the GC. This material must remain quite cool (T~10^2K), but does not catastrophically cool on a crossing time because of the high Lyman-Werner flux density in young GCs. The collection of gas within the GC must compete with ram pressure from the ambient ISM. After several 10^8 yr, the Lyman-Werner photon flux density drops by more than three orders of magnitude, allowing molecular hydrogen and then stars to form. After this second generation of star formation, type II SNe from the second generation and then prompt type Ia SNe associated with the first generation maintain a gas-free GC, thereby ending the cycle of star formation events. Our model makes clear predictions for the presence or absence of multiple stellar populations within GCs as a function of GC mass and formation environment. Analyzing intermediate-age LMC clusters, we find evidence for a mass threshold of ~10^4 Msun below which LMC clusters appear to be truly coeval. This threshold mass is consistent with our predictions for the mass at which ram pressure is capable of clearing gas from clusters in the LMC at the present epoch. (ABRIDGED)

We study the rich globular cluster (GC) system in the center of the massive cluster of galaxies Abell 1689 (z=0.18), one of the most powerful gravitational lenses known. With 28 HST/ACS orbits in the F814W bandpass, we reach magnitude I_814=29 with >90% completeness and sample the brightest ~5% of the GC system. Assuming the well-known Gaussian form of the GC luminosity function (GCLF), we estimate a total population of N(GC_total) = 162,850 GCs within a projected radius of 400kpc. As many as half may comprise an intracluster component. Even with the sizable uncertainties, which mainly result from the uncertain GCLF parameters, this is by far the largest GC system studied to date. The specific frequency S_N is high, but not uncommon for central galaxies in massive clusters, rising from S_N~5 near the center to ~12 at large radii. Passive galaxy fading would increase S_N by ~20% at z=0. We construct the radial mass profiles of the GCs, stars, intracluster gas, and lensing-derived total mass, and we compare the mass fractions as a function of radius. The estimated mass in GCs, M(GC_total)=3.9x10^10 Msun, is comparable to ~80% of the total stellar mass of the Milky Way. The shape of the GC mass profile appears intermediate between those of the stellar light and total cluster mass. Despite the extreme nature of this system, the ratios of the GC mass to the baryonic and total masses, and thus the GC formation efficiency, are typical of those in other rich clusters when comparing at the same physical radii. The GC formation efficiency is not constant, but varies with radius, in a manner that appears similar for different clusters; we speculate on the reasons for this similarity in profile.

Aims. The Fornax dwarf spheroidal galaxy (dSph) represents a challenge for some globular cluster (GC) formation models, because an exceptionally high fraction of its stellar mass is locked in its GC system. In order to shed light on our understanding of GC formation, we aim to constrain the amount of stellar mass that Fornax has lost via tidal interaction with the Milky Way (MW). Methods. Exploiting the flexibility of effective multi-component N-body simulations and relying on state-of-the-art estimates of Fornax’s orbital parameters, we study the evolution of the mass distribution of the Fornax dSph in observationally justified orbits in the gravitational potential of the MW over 12 Gyr. Results. We find that, though the dark-matter mass loss can be substantial, the fraction of stellar mass lost by Fornax to the MW is always negligible, even in the most eccentric orbit considered. Conclusions. We conclude that stellar-mass loss due to tidal stripping is not a plausible explanation for the unexpectedly high stellar mass of the GC system of the Fornax dSph and we discuss quantitatively the implications for GC formation scenarios.

XMM-Newton observations of faint X-ray sources in globular clusters

The results from high-resolution spectroscopy and accurate photometry have challenged the traditional paradigm that stars in globular clusters (GC) are simple stellar populations, rather suggesting that these structures harbor distinct groups of stars, differing in the chemical composition, particularly in the abundances of the light elements, from helium to silicon. Because this behavior is not shared by field stars, it is generally believed that some self-enrichment mechanism must have acted in GC, such that new stellar generations formed from the ashes of stars belonging to the original population. In this review, after presenting the state-of-the-art of the observations of GC stars, we discuss the possibility that the pollution of the intra-cluster medium was provided by the winds of AGB stars of initial mass above ∼3 M⊙. These objects evolve with time scales of 40 − 100 Myr and contaminate their surroundings with gas processed by p-capture nucleosynthesis, in agreement with the chemical patterns traced by GC stars.

Heavy elements, those produced by neutron-capture reactions, have traditionally shown no star-to-star dispersion in all but a handful of metal-poor globular clusters (GCs). Recent detections of low [Pb/Eu] ratios or upper limits in several metal-poor GCs indicate that the heavy elements in these GCs were produced exclusively by an r-process. Reexamining GC heavy element abundances from the literature, we find unmistakable correlations between the [La/Fe] and [Eu/Fe] ratios in 4 metal-poor GCs (M5, M15, M92, and NGC 3201), only 2 of which were known previously. This indicates that the total r-process abundances vary star-to-star (by factors of 2-6) relative to Fe within each GC. We also identify potential dispersion in two other GCs (M3 and M13). Several GCs (M12, M80, and NGC 6752) show no evidence of r-process dispersion. The r-process dispersion is not correlated with the well-known light element dispersion, indicating it was present in the gas throughout the duration of star formation. The observations available at present suggest that star-to-star r-process dispersion within metal-poor GCs may be a common but not ubiquitous phenomenon that is neither predicted by nor accounted for in current models of GC formation and evolution.

We have derived the mean proper motions and space velocities of 154 Galactic globular clusters and the velocity dispersion profiles of 141 globular clusters based on a combination of Gaia DR2 proper motions with ground-based line-of-sight velocities. Combining the velocity dispersion profiles derived here with new measurements of the internal mass functions allows us to model the internal kinematics of 144 clusters, more than 90% of the currently known Galactic globular cluster population. We also derive the initial cluster masses by calculating the cluster orbits backwards in time applying suitable recipes to account for mass loss and dynamical friction. We find a correlation between the stellar mass function of a globular cluster and the amount of mass lost from the cluster, pointing to dynamical evolution as one of the mechanisms shaping the mass function of stars in clusters. The mass functions also show strong evidence that globular clusters started with a bottom-light initial mass function. Our simulations show that the currently surviving globular cluster population has lost about 80% of its mass since the time of formation. If globular clusters started from a log-normal mass function, we estimate that the Milky Way contained about 500 globular clusters initially, with a combined mass of about $2.5 \cdot 10^8$ $M_\odot$. For a power-law initial mass function, the initial mass in globular clusters could have been a factor of three higher.

Using a recent homogeneous sample of 40 high-quality velocity dispersion profiles for Galactic globular clusters, we study the regime of low gravitational acceleration relevant to the outskirts of these systems. We find that a simple empirical profile having a central Gaussian component and a constant large-radius asymptote, σ∞, accurately describes the variety of observed velocity dispersion profiles. We use published population synthesis models, carefully tailored to each individual cluster, to estimate mass-to-light ratios from which total stellar masses, M, are inferred. We obtain a clear scaling, reminiscent of the galactic Tully–Fisher relation of $\sigma _{\infty }[\, \mathrm{km \, s}^{-1}]= 0.084^{+0.075}_{-0.040} (\mathrm{{\it M}/M}_{\odot })^{0.3 \pm 0.051}$, which is interesting to compare to the deep modified Newtonian dynamics (MOND) limit of $\sigma _{\infty } [\mathrm{km \, s}^{-1}]=0.2(\mathrm{{\it M}/M}_{\odot })^{0.25}$. Under a Newtonian interpretation, our results constitute a further restriction on models where initial conditions are crafted to yield the outer flattening observed today. Within a modified gravity scheme, because the globular clusters studied are not isolated objects in the deep MOND regime, the results obtained point towards a modified gravity where the external field effect of MOND does not appear, or is significantly suppressed.

In the last decade, the Kepler and CoRoT space-photometry missions have demonstrated the potential of asteroseismology as a novel, versatile and powerful tool to perform exquisite tests of stellar physics, and to enable precise and accurate characterisations of stellar properties, with impact on both exoplanetary and Galactic astrophysics. Based on our improved understanding of the strengths and limitations of such a tool, we argue for a new small/medium space mission dedicated to gathering high-precision, high-cadence, long photometric series in dense stellar fields. Such a mission will lead to breakthroughs in stellar astrophysics, especially in the metal poor regime, will elucidate the evolution and formation of open and globular clusters, and aid our understanding of the assembly history and chemodynamics of the Milky Way’s bulge and a few nearby dwarf galaxies.

We present an empirical model for the number of globular clusters (GCs) in galaxies based on recent data showing a tight relationship between dark matter halo virial masses and GC numbers. While a simple base model forming GCs in low-mass haloes reproduces this relation, we show that a second formation pathway for GCs is needed to account for observed younger GC populations. We confirm previous works that reported the observed linear correlation as being a consequence of hierarchical merging and its insensitivity to the exact GC formation processes at higher virial masses, even for a dual formation scenario. We find that the scatter of the linear relation is strongly correlated with the relative amount of smooth accretion: the more dark matter is smoothly accreted, the fewer GCs a halo has compared to other haloes of the same mass. This scatter is smaller than that introduced by halo mass measurements, indicating that the number of GCs in a galaxy is a good tracer for its dark matter mass. Smooth accretion is also the reason for a lower average dark matter mass per GC in low-mass haloes. Finally, we successfully reproduce the observed general trend of GCs being old and the tendency of more massive haloes hosting older GC systems. Including the second GC formation mechanism through gas-rich mergers leads to a more realistic variety of GC age distributions and also introduces an age inversion in the halo virial mass range log Mvir/M⊙ = 11–13.

The discovery of young globular clusters in merging galaxies and other environments provides an opportunity to study directly the process of globular cluster formation. Empirically it appears that globular cluster formation occurs preferentially in regions in which star formation occurs at a high rate and efficiency. Further, the interstellar medium in such regions is likely to be at a higher pressure than less active star-forming environments. An additional observational clue to the globular cluster formation process is that young globular clusters have little or no mass-radius relationship. In this paper I argue that high pressure and high star-formation efficiency are responsible for current globular cluster formation. I suggest that the precursors to globular clusters are molecular clouds and that the mass-radius relationship exhibited by such clouds is wiped out by a variable star formation efficiency.

We present detailed radiative transfer simulations of the reionization of the Milky Way by metal-poor globular clusters. We identify potential metal-poor globular cluster candidates within the Aquarius simulation using dark matter halo velocity dispersions. We calculate the local ionization fields via a photon-conserving, three dimensional non-equilibrium chemistry code. The key feature of the model is that globular cluster formation is suppressed if the local gas is ionized. We assume that at these early times, the ionization field is dominated by the flux from metal-poor globular clusters. Our spatial treatment of the ionization field leads to drastically different numbers and spatial distributions when compared to models where globular cluster formation is simply truncated at early redshifts (z ~ 13). The spatial distributions are more extended and more globular clusters are produced. We find that additional sources of ionization are required at later epochs (z ~ 10) to ionize the remaining gas and recover radial distributions statistically consistent with that of the Milky Way metal-poor globular clusters. We investigate a range of plausible ionization efficiencies to determine the effect photon-rich and photon-poor models have on present-day globular cluster properties. If globular clusters do indeed form within high-redshift dark matter haloes, they produce enough photons to ionize 98 and 90 per cent local (i.e. 23 h-3 Mpc3 centred on the host galaxy) volume and mass by redshift 10, respectively. In our photon-poorest model, this contribution drops to 60 and 50 per cent. Our model therefore implies that globular clusters are important contributors to the reionization process on local scales at high-redshift until more photon-rich sources dominate the photon budget at later times. The surviving clusters in all models have a narrow average age range (mean = 13.34 Gyr, s = 0.04 Gyr) consistent with current age estimates of the Milky Way metal-poor globular clusters. We also test a simple dynamical destruction model and estimate that ~60 per cent of all metal-poor globular clusters formed at high redshift have since been destroyed via tidal interactions with the host galaxy.