We present the first results from our HST Brightest Cluster Galaxy (BCG) survey of seven central supergiant cluster galaxies and their globular cluster (GC) systems. We measure a total of 48000 GCs in all seven galaxies, representing the largest single GC database. We find that a log-normal shape accurately matches the observed luminosity function (LF) of the GCs down to the GCLF turnover point, which is near our photometric limit. In addition, the LF has a virtually identical shape in all seven galaxies. Our data underscore the similarity in the formation mechanism of massive star clusters in diverse galactic environments. At the highest luminosities (log L > 10^7 L_Sun) we find small numbers of "superluminous" objects in five of the galaxies; their luminosity and color ranges are at least partly consistent with those of UCDs (Ultra-Compact Dwarfs). Lastly, we find preliminary evidence that in the outer halo (R > 20 kpc), the LF turnover point shows a weak dependence on projected distance, scaling as L_0 ~ R^-0.2, while the LF dispersion remains nearly constant.

We use data from the Pan-Andromeda Archaeological Survey (PAndAS) to search for evidence of an extended halo component belonging to M33 (the Triangulum Galaxy). We identify a population of red giant branch (RGB) stars at large radii from M33's disk whose connection to the recently discovered extended "disk substructure" is ambiguous, and which may represent a "bona-fide" halo component. After first correcting for contamination from the Milky Way foreground population and misidentified background galaxies, we average the radial density of RGB candidate stars over circular annuli centered on the galaxy and away from the disk substructure. We find evidence of a low-luminosity, centrally concentrated component that is everywhere in our data fainter than mu_V ~ 33 mag arcsec^(-2). The scale length of this feature is not well constrained by our data, but it appears to be of order r_exp ~ 20 kpc; there is weak evidence to suggest it is not azimuthally symmetric. Inspection of the overall CMD for this region that specifically clips out the disk substructure reveals that this residual RGB population is consistent with an old population with a photometric metallicity of around [Fe/H] ~ -2 dex, but some residual contamination from the disk substructure appears to remain. We discuss the likelihood that our findings represent a bona-fide halo in M33, rather than extended emission from the disk substructure. We interpret our findings in terms of an upper limit to M33's halo that is a few percent of its total luminosity, although its actual luminosity is likely much less.

We use CFHT/MegaCam data to search for outer halo star clusters in M33 as part of the Pan-Andromeda Archaeological Survey (PAndAS). This work extends previous studies out to a projected radius of 50 kpc and covers over 40 square degrees. We find only one new unambiguous star cluster in addition to the five previously known in the M33 outer halo (10 kpc <= r <= 50 kpc). Although we identify 2440 cluster candidates of various degrees of confidence from our objective image search procedure, almost all of these are likely background contaminants, mostly faint unresolved galaxies. We measure the luminosity, color and structural parameters of the new cluster in addition to the five previously-known outer halo clusters. At a projected radius of 22 kpc, the new cluster is slightly smaller, fainter and redder than all but one of the other outer halo clusters, and has g' ~ 19.9, (g'-i') ~ 0.6, concentration parameter c ~ 1.0, a core radius r_c ~ 3.5 pc, and a half-light radius r_h ~ 5.5 pc. For M33 to have so few outer halo clusters compared to M31 suggests either tidal stripping of M33's outer halo clusters by M31, or a very different, much calmer accretion history of M33.

In hierarchical cosmological models, galaxies grow in mass through the continual accretion of smaller ones. The tidal disruption of these systems is expected to result in loosely bound stars surrounding the galaxy, at distances that reach $10 - 100$ times the radius of the central disk. The number, luminosity and morphology of the relics of this process provide significant clues to galaxy formation history, but obtaining a comprehensive survey of these components is difficult because of their intrinsic faintness and vast extent. Here we report a panoramic survey of the Andromeda galaxy (M31). We detect stars and coherent structures that are almost certainly remnants of dwarf galaxies destroyed by the tidal field of M31. An improved census of their surviving counterparts implies that three-quarters of M31's satellites brighter than $M_V < -6$ await discovery. The brightest companion, Triangulum (M33), is surrounded by a stellar structure that provides persuasive evidence for a recent encounter with M31. This panorama of galaxy structure directly confirms the basic tenets of the hierarchical galaxy formation model and reveals the shared history of M31 and M33 in the unceasing build-up of galaxies.

We combine the globular cluster data for fifteen Brightest Cluster Galaxies and use this material to trace the mass-metallicity relations (MMR) in their globular cluster systems (GCSs). This work extends previous studies which correlate the properties of the MMR with those of the host galaxy. Our combined data sets show a mean trend for the metal-poor (MP) subpopulation which corresponds to a scaling of heavy-element abundance with cluster mass Z ~ M^(0.30+/-0.05). No trend is seen for the metal-rich (MR) subpopulation which has a scaling relation that is consistent with zero. We also find that the scaling exponent is independent of the GCS specific frequency and host galaxy luminosity, except perhaps for dwarf galaxies. We present new photometry in (g',i') obtained with Gemini/GMOS for the globular cluster populations around the southern giant ellipticals NGC 5193 and IC 4329. Both galaxies have rich cluster populations which show up as normal, bimodal sequences in the colour-magnitude diagram. We test the observed MMRs and argue that they are statistically real, and not an artifact caused by the method we used. We also argue against asymmetric contamination causing the observed MMR as our mean results are no different from other contamination-free studies. Finally, we compare our method to the standard bimodal fitting method (KMM or RMIX) and find our results are consistent. Interpretation of these results is consistent with recent models for globular cluster formation in which the MMR is determined by GC self-enrichment during their brief formation period.