Spectroscopy and Photometry of the Least Massive Type II Globular Clusters: NGC 1261 and NGC 6934*

Recent work has revealed two classes of globular clusters (GCs), dubbed Type I and Type II. Type II GCs are characterized by both a blue and a red red giant branch composed of stars with different metallicities, often coupled with distinct abundances in the slow neutron-capture elements (s-elements). Here we continue the chemical tagging of Type II GCs by adding the two least massive clusters of this class, NGC 1261 and NGC 6934. Based on both spectroscopy and photometry, we find red stars in NGC 1261 to be slightly enhanced in [Fe/H] by ∼0.1 dex and confirm that red stars of NGC 6934 are enhanced in iron by ∼0.2 dex. Neither NGC 1261 nor NGC 6934 show internal variations in the s-elements, which suggests a GC mass threshold for the occurrence of s-process enrichment. We found a significant correlation between the additional Fe locked in the red stars of Type II GCs and the present-day mass of the cluster. Nevertheless, most Type II GCs retained a small fraction of Fe produced by SNe II, lower than the 2%; NGC 6273, M54, and ω Centauri are remarkable exceptions. In the Appendix, we infer for the first time chemical abundances of lanthanum, assumed as representative of the s-elements, in M54, the GC located in the nucleus of the Sagittarius dwarf galaxy. Red-sequence stars are marginally enhanced in [La/Fe] by 0.10 ± 0.06 dex, in contrast with the large [La/Fe] spread of most Type II GCs. We suggest that different processes are responsible for the enrichment in iron and s-elements in Type II GCs.

RELICS: ICL Analysis of the z = 0.566 Merging Cluster WHL J013719.8–08284

We present a pilot study of the intracluster light (ICL) in massive clusters using imaging of the z = 0.566 cluster of galaxies WHL J013719.8–08284 observed by the RELICS project with the Hubble Space Telescope. We measure the ICL fraction in four optical ACS/WFC filters (F435W, F475W, F606W, and F814W) and five infrared WFC3/IR bands (F105W, F110W, F125W, F140W, and F160W). The ICL maps are calculated using the free-of-a-priori-assumptions algorithm CICLE, and the cluster membership is estimated from photometric properties. We find optical ICL fractions that range between ∼6% and 19%, in nice agreement with the values found in previous works for merging clusters. We also observe an ICL fraction excess between 3800 Å and 4800 Å, previously identified as a signature of merging clusters at 0.18 < z < 0.55. This excess suggests the presence of an enhanced population of young/low-metallicity stars in the ICL. All indicators thus point to WHL J013719.8–08284 as a disturbed cluster with a significant amount of recently injected stars, bluer than the average stars hosted by the cluster members and likely stripped out from infalling galaxies during the current merging event. Infrared ICL fractions are ∼50% higher than optical ones, which could be signatures of an older and/or higher-metallicity ICL population that can be associated with the buildup of the brightest cluster galaxy, passive evolution of previously injected young stars, or preprocessing in infalling groups. Finally, investigating the photometry of the cluster members, we tentatively conclude that WHL J013719.8–08284 fulfills the expected conditions for a fossil system progenitor.

Probing Galaxy Evolution in Massive Clusters Using ACT and DES: Splashback as a Cosmic Clock

We measure the projected number density profiles of galaxies and the splashback feature in clusters selected by the Sunyaev–Zel’dovich effect from the Advanced Atacama Cosmology Telescope (AdvACT) survey using galaxies observed by the Dark Energy Survey (DES). The splashback radius is consistent with CDM-only simulations and is located at 2.4 − 0.4 + 0.3 Mpc h − 1 . We split the galaxies on color and find significant differences in their profile shapes. Red and green-valley galaxies show a splashback-like minimum in their slope profile consistent with theory, while the bluest galaxies show a weak feature at a smaller radius. We develop a mapping of galaxies to subhalos in simulations and assign colors based on infall time onto their hosts. We find that the shift in location of the steepest slope and different profile shapes can be mapped to the average time of infall of galaxies of different colors. The steepest slope traces a discontinuity in the phase space of dark matter halos. By relating spatial profiles to infall time, we can use splashback as a clock to understand galaxy quenching. We find that red galaxies have on average been in clusters over 3.2 Gyr, green galaxies about 2.2 Gyr, while blue galaxies have been accreted most recently and have not reached apocenter. Using the full radial profiles, we fit a simple quenching model and find that the onset of galaxy quenching occurs after a delay of about a gigayear and that galaxies quench rapidly thereafter with an exponential timescale of 0.6 Gyr.

Shocks in the stacked Sunyaev–Zel’dovich profiles of clusters – I. Analysis with the Three Hundred simulations

ABSTRACT Gas infalling into the gravitational potential wells of massive galaxy clusters is expected to experience one or more shocks on its journey to becoming part of the intracluster medium (ICM). These shocks are important for setting the thermodynamic properties of the ICM and can therefore impact cluster observables such as X-ray emission and the Sunyaev–Zel’dovich (SZ) effect. We investigate the possibility of detecting signals from cluster shocks in the averaged thermal SZ profiles of galaxy clusters. Using zoom-in hydrodynamic simulations of massive clusters from the Three Hundred Project, we show that if cluster SZ profiles are stacked as a function of R/R200m, shock-induced features appear in the averaged SZ profile. These features are not accounted for in standard fitting formulae for the SZ profiles of galaxy clusters. We show that the shock features should be detectable with samples of clusters from ongoing and future SZ surveys. We also demonstrate that the location of these features is correlated with the cluster accretion rate, as well as the location of the cluster splashback radius. Analyses of ongoing and future surveys, such as SPT-3g, AdvACT, Simons Observatory, and CMB-S4, which include gas shocks will gain a new handle on the properties and dynamics of the outskirts of massive haloes, both in gas and in mass.

On the initial mass-radius relation of stellar clusters

Young stellar clusters across nearly five orders of magnitude in mass appear to follow a power-law mass-radius relationship (MRR), $R_{\star }\propto M_{\star }^{\alpha }$, with α ≈ 0.2 − 0.33. We develop a simple analytic model for the cluster mass-radius relation. We consider a galaxy disc in hydrostatic equilibrium, which hosts a population of molecular clouds that fragment into clumps undergoing cluster formation and feedback-driven expansion. The model predicts a mass-radius relation of $R_{\star }\propto M_{\star }^{1/2}$ and a dependence on the kpc-scale gas surface density $R_{\star }\propto \Sigma _{\rm g}^{-1/2}$, which results from the formation of more compact clouds (and cluster-forming clumps within) at higher gas surface densities. This environmental dependence implies that the high-pressure environments in which the most massive clusters can form also induce the formation of clusters with the smallest radii, thereby shallowing the observed MRR at high-masses towards the observed $R_{\star }\propto M_{\star }^{1/3}$. At low cluster masses, relaxation-driven expansion induces a similar shallowing of the MRR. We combine our predicted MRR with a simple population synthesis model and apply it to a variety of star-forming environments, finding good agreement. Our model predicts that the high-pressure formation environments of globular clusters at high redshift naturally led to the formation of clusters that are considerably more compact than those in the local Universe, thereby increasing their resilience to tidal shock-driven disruption and contributing to their survival until the present day.

The history of metal enrichment traced by X-ray observations of high redshift galaxy clusters

We present the analysis of deep X-ray observations of 10 massive galaxy clusters at redshifts 1.05 &lt; z &lt; 1.71, with the primary goal of measuring the metallicity of the intracluster medium (ICM) at intermediate radii, to better constrain models of the metal enrichment of the intergalactic medium. The targets were selected from X-ray and Sunyaev-Zel’dovich (SZ) effect surveys, and observed with both the XMM-Newton and Chandra satellites. For each cluster, a precise gas mass profile was extracted, from which the value of r500 could be estimated. This allows us to define consistent radial ranges over which the metallicity measurements can be compared. In general, the data are of sufficient quality to extract meaningful metallicity measurements in two radial bins, r &lt; 0.3r500 and 0.3 &lt; r/r500 &lt; 1.0. For the outer bin, the combined measurement for all ten clusters, Z/Z⊙ = 0.21 ± 0.09, represents a substantial improvement in precision over previous results. This measurement is consistent with, but slightly lower than, the average metallicity of 0.315 Solar measured at intermediate-to-large radii in low-redshift clusters. Combining our new high-redshift data with the previous low-redshift results allows us to place the tightest constraints to date on models of the evolution of cluster metallicity at intermediate radii. Adopting a power law model of the form Z∝(1 + z)γ, we measure a slope $\gamma = -0.5^{+0.4}_{-0.3}$, consistent with the majority of the enrichment of the ICM having occurred at very early times and before massive clusters formed, but leaving open the possibility that some additional enrichment in these regions may have occurred since a redshift of 2.

Exploring the role of cosmological shock waves in the Dianoga simulations of galaxy clusters

Cosmological shock waves are ubiquitous to cosmic structure formation and evolution. As a consequence, they play a major role in the energy distribution and thermalization of the intergalactic medium (IGM). We analyse the Mach number distribution in the Dianoga simulations of galaxy clusters performed with the SPH code GADGET-3. The simulations include the effects of radiative cooling, star formation, metal enrichment, supernova and active galactic nuclei feedback. A grid-based shock-finding algorithm is applied in post-processing to the outputs of the simulations. This procedure allows us to explore in detail the distribution of shocked cells and their strengths as a function of cluster mass, redshift and baryonic physics. We also pay special attention to the connection between shock waves and the cool-core/non-cool core (CC/NCC) state and the global dynamical status of the simulated clusters. In terms of general shock statistics, we obtain a broad agreement with previous works, with weak (low-Mach number) shocks filling most of the volume and processing most of the total thermal energy flux. As a function of cluster mass, we find that massive clusters seem more efficient in thermalising the IGM and tend to show larger external accretion shocks than less massive systems. We do not find any relevant difference between CC and NCC clusters. However, we find a mild dependence of the radial distribution of the shock Mach number on the cluster dynamical state, with disturbed systems showing stronger shocks than regular ones throughout the cluster volume.