Globular Cluster Systems in Brightest Cluster Galaxies. III: Beyond Bimodality

We present new deep photometry of the rich globular cluster (GC) systems around the Brightest Cluster Galaxies UGC 9799 (Abell 2052) and UGC 10143 (Abell 2147), obtained with the HST ACS and WFC3 cameras. For comparison, we also present new reductions of similar HST/ACS data for the Coma supergiants NGC 4874 and 4889. All four of these galaxies have huge cluster populations (to the radial limits of our data, comprising from 12000 to 23000 clusters per galaxy). The metallicity distribution functions (MDFs) of the GCs can still be matched by a bimodal-Gaussian form where the metal-rich and metal-poor modes are separated by ~0.8 dex, but the internal dispersions of each mode are so large that the total MDF becomes very broad and nearly continuous from [Fe/H] = -2.4 to Solar. There are, however, significant differences between galaxies in the relative numbers of \emph{metal-rich} clusters, suggesting that they underwent significantly different histories of mergers with massive, gas-rich halos. Lastly, the proportion of metal-poor GCs rises especially rapidly outside projected radii R > 4 R_eff, suggesting the importance of accreted dwarf satellites in the outer halo. Comprehensive models for the formation of GCs as part of the hierarchical formation of their parent galaxies will be needed to trace the systematic change in structure of the MDF with galaxy mass, from the distinctly bimodal form in smaller galaxies up to the broad continuum that we see in the very largest systems.