• We present Keck/DEIMOS spectroscopy of globular clusters (GCs) around the ultra-diffuse galaxies (UDGs) VLSB-B, VLSB-D, and VCC615 located in the central regions of the Virgo cluster. We spectroscopically identify 4, 12, and 7 GC satellites of these UDGs, respectively. We find that the three UDGs have systemic velocities ($V_{sys}$) consistent with being in the Virgo cluster, and that they span a wide range of velocity dispersions, from $\sim 16$ to $\sim 47$ km/s, and high dynamical mass-to-light ratios within the radius that contains half the number of GCs ($ 407^{+916}_{-407}$, $21^{+15}_{-11}$, $60^{+65}_{-38}$, respectively). VLSB-D shows possible evidence for rotation along the stellar major axis and its $V_{sys}$ is consistent with that of the massive galaxy M84 and the center of the Virgo cluster itself. These findings, in addition to having a dynamically and spatially ($\sim 1$ kpc) off-centered nucleus and being extremely elongated, suggest that VLSB-D could be tidally perturbed. On the contrary, VLSB-B and VCC615 show no signals of tidal deformation. Whereas the dynamics of VLSB-D suggest that it has a less massive dark matter halo than expected for its stellar mass, VLSB-B and VCC615 are consistent with a $\sim 10^{12}$ M$_{\odot}$ dark matter halo. Although our samples of galaxies and GCs are small, these results suggest that UDGs may be a diverse population, with their low surface brightnesses being the result of very early formation, tidal disruption, or a combination of the two.
  • We derive stellar population parameters for a representative sample of ultracompact dwarfs (UCDs) and a large sample of massive globular clusters (GCs) with stellar masses $\gtrsim$ 10$^{6}$ $M_{\odot}$ in the central galaxy M87 of the Virgo galaxy cluster, based on model fitting to the Lick-index measurements from both the literature and new observations. After necessary spectral stacking of the relatively faint objects in our initial sample of 40 UCDs and 118 GCs, we obtain 30 sets of Lick-index measurements for UCDs and 80 for GCs. The M87 UCDs have ages $\gtrsim$ 8 Gyr and [$\alpha$/Fe] $\simeq$ 0.4 dex, in agreement with previous studies based on smaller samples. The literature UCDs, located in lower-density environments than M87, extend to younger ages and smaller [$\alpha$/Fe] (at given metallicities) than M87 UCDs, resembling the environmental dependence of the Virgo dE nuclei. The UCDs exhibit a positive mass-metallicity relation (MZR), which flattens and connects compact ellipticals at stellar masses $\gtrsim$ 10$^{8}$ $M_{\odot}$. The Virgo dE nuclei largely follow the average MZR of UCDs, whereas most of the M87 GCs are offset towards higher metallicities for given stellar masses. The difference between the mass-metallicity distributions of UCDs and GCs may be qualitatively understood as a result of their different physical sizes at birth in a self-enrichment scenario or of galactic nuclear cluster star formation efficiency being relatively low in a tidal stripping scenario for UCD formation. The existing observations provide the necessary but not sufficient evidence for tidally stripped dE nuclei being the dominant contributors to the M87 UCDs.
  • Substructure in globular cluster (GC) populations around large galaxies is expected in galaxy formation scenarios that involve accretion or merger events, and it has been searched for using direct associations between GCs and structure in the diffuse galaxy light, or with GC kinematics. Here, we present a search for candidate substructures in the GC population around the Virgo cD galaxy M87 through the analysis of the spatial distribution of the GC colors.~The study is based on a sample of $\sim\!1800$ bright GCs with high-quality $u,g,r,i,z,K_s$ photometry, selected to ensure a low contamination by foreground stars or background galaxies.~The spectral energy distributions of the GCs are associated with formal estimates of age and metallicity, which are representative of its position in a 4-D color-space relative to standard single stellar population models.~Dividing the sample into broad bins based on the relative formal ages, we observe inhomogeneities which reveal signatures of GC substructures.~The most significant of these is a spatial overdensity of GCs with relatively young age labels, of diameter $\sim\!0.1$\,deg ($\sim\!30\,$kpc), located to the south of M87.~The significance of this detection is larger than about 5$\sigma$ after accounting for estimates of random and systematic errors.~Surprisingly, no large Virgo galaxy is present in this area, that could potentially host these GCs.~But candidate substructures in the M87 halo with equally elusive hosts have been described based on kinematic studies in the past.~The number of GC spectra available around M87 is currently insufficient to clarify the nature of the new candidate substructure.
  • We present a photometric study of the dwarf galaxy population in the core region ($< r_{\rm vir}/4$) of the Fornax galaxy cluster based on deep $u'g'i'$ photometry from the Next Generation Fornax Cluster Survey. All imaging data were obtained with the Dark Energy Camera mounted on the 4-meter Blanco telescope at the Cerro-Tololo Interamerican Observatory. We identify 258 dwarf galaxy candidates with luminosities $-17 < M_{g'} < -8$ mag, corresponding to typical stellar masses of $9.5\gtrsim \log{\cal M}_{\star}/M_\odot \gtrsim 5.5$, reaching $\sim\!3$ mag deeper in point-source luminosity and $\sim\!4$ mag deeper in surface-brightness sensitivity compared to the classic Fornax Cluster Catalog. Morphological analysis shows that surface-brightness profiles are well represented by single-component S\'ersic models with average S\'ersic indices of $\langle n\rangle_{u',g',i'}=(0.78-0.83) \pm 0.02$, and average effective radii of $\langle r_e\rangle_{u',g',i'}\!=(0.67-0.70) \pm 0.02$ kpc. Color-magnitude relations indicate a flattening of the galaxy red sequence at faint galaxy luminosities, similar to the one recently discovered in the Virgo cluster. A comparison with population synthesis models and the galaxy mass-metallicity relation reveals that the average faint dwarf galaxy is likely older than ~5 Gyr. We study galaxy scaling relations between stellar mass, effective radius, and stellar mass surface density over a stellar mass range covering six orders of magnitude. We find that over the sampled stellar mass range several distinct mechanisms of galaxy mass assembly can be identified: i) dwarf galaxies assemble mass inside the half-mass radius up to $\log{\cal M}_{\star}$ ~8.0, ii) isometric mass assembly in the range $8.0 < \log{\cal M}_{\star}/M_\odot < 10.5$, and iii) massive galaxies assemble stellar mass predominantly in their halos at $\log{\cal M}_{\star}$ ~10.5 and above.
  • We present an analysis of 39 nuclei and their early-type hosts in the Virgo Cluster using ten broadband filters: F300W, F475W, F850LP, F160W, $u^*griz$, and $K_s$. We describe the Virgo Redux program, which provides high-resolution UV and NIR imaging. Combining this data with optical and NIR imaging from the ACS Virgo Cluster Survey and the Next Generation Virgo Cluster Survey, we estimate masses, metallicities and ages using simple stellar population (SSP) models. For 19 nuclei, we compare to SSP parameters derived from Keck and Gemini spectra and find reasonable agreement between the photometric and spectroscopic metallicity: the RMS scatter is 0.3 dex. We reproduce the nucleus-galaxy mass fraction of $0.33^{+0.09}_{-0.07}$ percent for galaxy stellar masses $10^{8.4}-10^{10.3} M_\odot$ with a typical precision of $\sim$35% for the nuclei masses. Based on available model predictions, there is no single preferred formation scenario for nuclei, suggesting that nuclei are formed stochastically through a mix of processes. Nuclei metallicities are statistically identical to those of their hosts, appearing $0.07 \pm 0.3$ dex more metal-rich on average -- although, omitting galaxies with unusual origins, nuclei are $0.20\pm0.28$ dex more metal-rich. Nuclei appear to be $0.56 \pm 0.12$ dex more metal rich than ultra-compact dwarf galaxies (UCDs) at fixed mass. We find no clear age difference between nuclei and their galaxies, with nuclei displaying a broad range of ages. Interestingly, we find that the most massive nuclei may be flatter and more closely aligned with the semi-major axes of their hosts, suggesting that they formed through predominantly dissipative processes.
  • The Canada-France Imaging Survey (CFIS) will map the northern high Galactic latitude sky in the $u$-band ("CFIS-u", 10,000$\, {\rm deg^2}$) and in the $r$-band ("CFIS-r", 5,000$\, {\rm deg^2}$), enabling a host of stand-alone science investigations, and providing some of the ground-based data necessary for photometric redshift determination for the Euclid mission. In this first contribution we present the $u$-band component of the survey, describe the observational strategy, and discuss some first highlight results, based on approximately one third of the final area. We show that the Galactic anticenter structure is distributed continuously along the line of sight, out to beyond 20 kpc, and possesses a metallicity distribution that is essentially identical to that of the outer disk sampled by APOGEE. This suggests that it is probably a buckled disk of old metal-rich stars, rather than a stream or a flare. We also discuss the future potential for CFIS-u in discovering star-forming dwarf galaxies around the Local Group, the characterization of the white dwarf and blue straggler population of the Milky Way, as well as its sensitivity to low-surface brightness structures in external galaxies.
  • We present the chemical distribution of the Milky Way, based on 2,900$\, {\rm deg^2}$ of $u$-band photometry taken as part of the Canada-France Imaging Survey. When complete, this survey will cover 10,000$\, {\rm deg^2}$ of the Northern sky. By combing the CFHT $u$-band photometry together with SDSS and Pan-STARRS $g,r,$ and $i$, we demonstrate that we are able to measure reliably the metallicities of individual stars to $\sim 0.2$ dex, and hence additionally obtain good photometric distance estimates. This survey thus permits the measurement of metallicities and distances of the dominant main-sequence population out to approximately 30 kpc, and provides much higher number of stars at large extraplanar distances than have been available from previous surveys. We develop a non-parametric distance-metallicity decomposition algorithm and apply it to the sky at $30\deg < |b| < 70\deg$ and to the North Galactic Cap. We find that the metallicity-distance distribution is well-represented by three populations whose metallicity distributions do not vary significantly with vertical height above the disk. As traced in main-sequence stars, the stellar halo component shows a vertical density profile that is close to exponential, with a scale height of around 3 kpc. This may indicate that the inner halo was formed partly from disk stars ejected in an ancient minor merger.
  • We build a background cluster candidate catalog from the Next Generation Virgo Cluster Survey, using our detection algorithm RedGOLD. The NGVS covers 104$deg^2$ of the Virgo cluster in the $u^*,g,r,i,z$-bandpasses to a depth of $ g \sim 25.7$~mag (5$\sigma$). Part of the survey was not covered or has shallow observations in the $r$--band. We build two cluster catalogs: one using all bandpasses, for the fields with deep $r$--band observations ($\sim 20 \ deg^2$), and the other using four bandpasses ($u^*,g,i,z$) for the entire NGVS area. Based on our previous CFHT-LS W1 studies, we estimate that both of our catalogs are $\sim100\%$($\sim70\%$) complete and $\sim80\%$ pure, at $z\le 0.6$($z\lesssim1$), for galaxy clusters with masses of $M\gtrsim10^{14}\ M_{\odot}$. We show that when using four bandpasses, though the photometric redshift accuracy is lower, RedGOLD detects massive galaxy clusters up to $z\sim 1$ with completeness and purity similar to the five-band case. This is achieved when taking into account the bias in the richness estimation, which is $\sim40\%$ lower at $0.5\le z<0.6$ and $\sim20\%$ higher at $0.6<z< 0.8$, with respect to the five-band case. RedGOLD recovers all the X-ray clusters in the area with mass $M_{500} > 1.4 \times 10^{14} \rm M_{\odot}$ and $0.08<z<0.5$. Because of our different cluster richness limits and the NGVS depth, our catalogs reach to lower masses than the published redMaPPer cluster catalog over the area, and we recover $\sim 90-100\%$ of its detections.
  • Large samples of globular clusters (GC) with precise multi-wavelength photometry are becoming increasingly available and can be used to constrain the formation history of galaxies. We present the results of an analysis of Milky Way (MW) and Virgo core GCs based on five optical-near-infrared colors and ten synthetic stellar population models. For the MW GCs, the models tend to agree on photometric ages and metallicities, with values similar to those obtained with previous studies. When used with Virgo core GCs, for which photometry is provided by the Next Generation Virgo cluster Survey (NGVS), the same models generically return younger ages. This is a consequence of the systematic differences observed between the locus occupied by Virgo core GCs and models in panchromatic color space. Only extreme fine-tuning of the adjustable parameters available to us can make the majority of the best-fit ages old. Although we cannot exclude that the formation history of the Virgo core may lead to more conspicuous populations of relatively young GCs than in other environments, we emphasize that the intrinsic properties of the Virgo GCs are likely to differ systematically from those assumed in the models. Thus, the large wavelength coverage and photometric quality of modern GC samples, such as used here, is not by itself sufficient to better constrain the GC formation histories. Models matching the environment-dependent characteristics of GCs in multi-dimensional color space are needed to improve the situation.
  • Wide-field $u'g'r'i'z'$ Dark Energy Camera observations centred on the giant elliptical galaxy NGC5128 covering $\sim21deg^2$ are used to compile a new catalogue of $\sim3200$ globular clusters (GCs). We report 2404 new candidates, including the vast majority within $\sim140$kpc of NGC5128. We find evidence for a transition at a galactocentric radius of $R_{\rm gc}\approx55$kpc from GCs intrinsic to NGC5128 to those likely to have been accreted from dwarf galaxies or that may transition to the intra-group medium of the Centaurus A galaxy group. We fit power-law surface number density profiles of the form $\Sigma_{N, R_{\rm gc}}\propto R_{\rm gc}^\Gamma$ and find that inside the transition radius, the red GCs are more centrally concentrated than the blue, with $\Gamma_{\rm inner,red}\approx-1.78$ and $\Gamma_{\rm inner,blue}\approx-1.40$. Outside this region both profiles flatten, more dramatically for the red GCs ($\Gamma_{\rm outer,red}\approx-0.33$) compared to the blue ($\Gamma_{\rm outer,blue}\approx-0.61$), although the former is more likely to suffer contamination by background sources. The median $(g'\!-\!z')_0\!=\!1.27$mag colour of the inner red population is consistent with arising from the amalgamation of two giant galaxies each less luminous than present-day NGC5128. Both in- and out-ward of the transition radius, we find the fraction of blue GCs to dominate over the red GCs, indicating a lively history of minor-mergers. Assuming the blue GCs to originate primarily in dwarf galaxies, we model the population required to explain them, while remaining consistent with NGC5128's present-day spheroid luminosity. We find that several dozen dwarfs of luminosities $L_{dw,V}\simeq10^{6-9.3}L_{V,\odot}$, following a Schechter luminosity function with a faint-end slope of $-1.50\leq\alpha\leq-1.25$ is favoured, many of which may have already been disrupted in NGC5128's tidal field.
  • Globular clusters (GCs) are some of the most visible tracers of the merging and accretion history of galaxy halos. Metal-poor GCs, in particular, are thought to arrive in massive galaxies largely through dry, minor merging events, but it is rare to see a direct connection between GCs and visible stellar streams. NGC 474 is a post-merger early-type galaxy with dramatic fine structures made of concentric shells and radial streams that have been more clearly revealed by deep imaging. We present a study of GCs in NGC 474 to better establish the relationship between merger-induced fine structure and the GC system. We find that many GCs are superimposed on visible streams and shells, and about 35% of GCs outside $3R_{\rm e,galaxy}$ are located in regions of fine structure. The spatial correlation between the GCs and fine structure is significant at the 99.9% level, showing that this correlation is not coincidental. The colors of the GCs on the fine structures are mostly blue, and we also find an intermediate-color population that is dominant in the central region, and which will likely passively evolve to have colors consistent with a traditional metal-rich GC population. The association of the blue GCs with fine structures is direct confirmation that many metal-poor GCs are accreted onto massive galaxy halos through merging events, and that progenitors of these mergers are sub-L* galaxies.
  • We use deep optical photometry from the Next Generation Virgo Cluster Survey [NGVS] to investigate the color-magnitude diagram for the galaxies inhabiting the core of this cluster. The sensitivity of the NGVS imaging allows us to continuously probe galaxy colors over a factor of $\sim 2 \times 10^5$ in luminosity, from brightest cluster galaxies to scales overlapping classical satellites of the Milky Way [$M_{g^{\prime}}$ $\sim$ $-$9; $M_{*}$ $\sim 10^6$ M$_{\odot}$], within a single environment. Remarkably, we find the first evidence that the RS flattens in all colors at the faint-magnitude end [starting between $-$14 $\le$ $M_{g^{\prime}}$ $\le$ $-$13, around $M_{*}$ $\sim 4 \times 10^7$ M$_{\odot}$], with the slope decreasing to $\sim$60% or less of its value at brighter magnitudes. This could indicate that the stellar populations of faint dwarfs in Virgo's core share similar characteristics [e.g. constant mean age] over $\sim$3 mags in luminosity, suggesting that these galaxies were quenched coevally, likely via pre-processing in smaller hosts. We also compare our results to galaxy formation models, finding that the RS in model clusters have slopes at intermediate magnitudes that are too shallow, and in the case of semi-analytic models, do not reproduce the flattening seen at both extremes [bright/faint] of the Virgo RS. Deficiencies in the chemical evolution of model galaxies likely contribute to the model-data discrepancies at all masses, while overly efficient quenching may also be a factor at dwarf scales. Deep UV and near-IR photometry are required to unambiguously diagnose the cause of the faint-end flattening.
  • We present an analysis of high-quality photometry for globular clusters (GCs) in the Virgo cluster core region, based on data from the Next Generation Virgo Cluster Survey (NGVS) pilot field, and in the Milky Way (MW) based on VLT/X-Shooter spectrophotometry. We find significant discrepancies in color-color diagrams between sub-samples from different environments, confirming that the environment has a strong influence on the integrated colors of GCs. GC color distributions along a single color are not sufficient to capture the differences we observe in color-color space. While the average photometric colors become bluer with increasing radial distance to the cD galaxy M87, we also find a relation between the environment and the slope and intercept of the color-color relations. A denser environment seems to produce a larger dynamic range in certain color indices. We argue that these results are not due solely to differential extinction, IMF variations, calibration uncertainties, or overall age/metallicity variations. We therefore suggest that the relation between the environment and GC colors is, at least in part, due to chemical abundance variations, which affect stellar spectra and stellar evolution tracks. Our results demonstrate that stellar population diagnostics derived from model predictions which are calibrated on one particular sample of GCs may not be appropriate for all extragalactic GCs. These results advocate a more complex model of the assembly history of GC systems in massive galaxies that goes beyond the simple bimodality found in previous decades.
  • We report the discovery of a very diverse set of five low-surface brightness (LSB) dwarf galaxy candidates in Hickson Compact Group 90 (HCG 90) detected in deep U- and I-band images obtained with VLT/VIMOS. These are the first LSB dwarf galaxy candidates found in a compact group of galaxies. We measure spheroid half-light radii in the range $0.7\!\lesssim\! r_{\rm eff}/{\rm kpc}\! \lesssim\! 1.5$ with luminosities of $-11.65\!\lesssim\! M_U\! \lesssim\! -9.42$ and $-12.79\!\lesssim\! M_I\! \lesssim\! -10.58$ mag, corresponding to a color range of $(U\!-\!I)_0\!\simeq\!1.1\!-\!2.2$ mag and surface brightness levels of $\mu_U\!\simeq\!28.1\,{\rm mag/arcsec^2}$ and $\mu_I\!\simeq\!27.4\,{\rm mag/arcsec^2}$. Their colours and luminosities are consistent with a diverse set of stellar population properties. Assuming solar and 0.02 Z$_\odot$ metallicities we obtain stellar masses in the range $M_*|_{Z_\odot} \simeq 10^{5.7-6.3} M_{\odot}$ and $M_*|_{0.02\,Z_\odot}\!\simeq\!10^{6.3-8}\,M_{\odot}$. Three dwarfs are older than 1 Gyr, while the other two significantly bluer dwarfs are younger than $\sim 2$ Gyr at any mass/metallicity combination. Altogether, the new LSB dwarf galaxy candidates share properties with dwarf galaxies found throughout the Local Volume and in nearby galaxy clusters such as Fornax. We find a pair of candidates with $\sim\!2$ kpc projected separation, which may represent one of the closest dwarf galaxy pairs found. We also find a nucleated dwarf candidate, with a nucleus size of $r_{\rm eff}\!\simeq\!46\!-\!63$ pc and magnitude M$_{U,0}=-7.42$ mag and $(U\!-\!I)_0\!=\!1.51$ mag, which is consistent with a nuclear stellar disc with a stellar mass in the range $10^{4.9-6.5}\,M_\odot$.
  • The central region of the Virgo cluster of galaxies contains thousands of globular clusters (GCs), an order of magnitude more than the numbers found in the Local Group. Relics of early star formation epochs in the universe, these GCs also provide ideal targets to test our understanding of the Spectral Energy Distributions (SEDs) of old stellar populations. Based on photometric data from the Next Generation Virgo cluster Survey (NGVS) and its near-infrared counterpart NGVS-IR, we select a robust sample of 1846 GCs with excellent photometry and spanning the full range of colors present in the Virgo core. The selection exploits the well defined locus of GCs in the uiK diagram and the fact that the globular clusters are marginally resolved in the images. We show that the GCs define a narrow sequence in 5-dimensional color space, with limited but real dispersion around the mean sequence. The comparison of these SEDs with the predictions of eleven widely used population synthesis models highlights differences between models, and also shows that no single model adequately matches the data in all colors. We discuss possible causes for some of these discrepancies. Forthcoming papers of this series will examine how best to estimate photometric metallicities in this context, and compare the Virgo globular cluster colors with those in other environments.
  • We present measurements of the galaxy luminosity and stellar mass function in a 3.71 deg$^2$ (0.3 Mpc$^2$) area in the core of the Virgo cluster, based on $ugriz$ data from the Next Generation Virgo Cluster Survey (NGVS). The galaxy sample consists of 352 objects brighter than $M_g=-9.13$ mag, the 50% completeness limit of the survey. Using a Bayesian analysis, we find a best-fit faint end slope of $\alpha=-1.33 \pm 0.02$ for the g-band luminosity function; consistent results are found for the stellar mass function as well as the luminosity function in the other four NGVS bandpasses. We discuss the implications for the faint-end slope of adding 92 ultra compact dwarfs galaxies (UCDs) -- previously compiled by the NGVS in this region -- to the galaxy sample, assuming that UCDs are the stripped remnants of nucleated dwarf galaxies. Under this assumption, the slope of the luminosity function (down to the UCD faint magnitude limit, $M_g = -9.6$ mag) increases dramatically, up to $\alpha = -1.60 \pm 0.06$ when correcting for the expected number of disrupted non-nucleated galaxies. We also calculate the total number of UCDs and globular clusters that may have been deposited in the core of Virgo due to the disruption of satellites, both nucleated and non-nucleated. We estimate that ~150 objects with $M_g\lesssim-9.6$ mag and that are currently classified as globular clusters, might, in fact, be the nuclei of disrupted galaxies. We further estimate that as many as 40% of the (mostly blue) globular clusters in the core of Virgo might once have belonged to such satellites; these same disrupted satellites might have contributed ~40% of the total luminosity in galaxies observed in the core region today. Finally, we use an updated Local Group galaxy catalog to provide a new measurement of the luminosity function of Local Group satellites, $\alpha=-1.21\pm0.05$.
  • We use imaging from the Next Generation Virgo cluster Survey (NGVS) to present a comparative study of ultra-compact dwarf (UCD) galaxies associated with three prominent Virgo sub-clusters: those centered on the massive, red-sequence galaxies M87, M49 and M60. We show how UCDs can be selected with high completeness using a combination of half-light radius and location in color-color diagrams ($u^*iK_s$ or $u^*gz$). Although the central galaxies in each of these sub-clusters have nearly identical luminosities and stellar masses, we find large differences in the sizes of their UCD populations, with M87 containing ~3.5 and 7.8 times more UCDs than M49 and M60, respectively. The relative abundance of UCDs in the three regions scales in proportion to sub-cluster mass, as traced by X-ray gas mass, total gravitating mass, number of globular clusters, and number of nearby galaxies. We find that the UCDs are predominantly blue in color, with ~85% of the UCDs having colors similar to blue GCs and stellar nuclei of dwarf galaxies. We present evidence that UCDs surrounding M87 and M49 may follow a morphological sequence ordered by the prominence of their outer, low surface brightness envelope, ultimately merging with the sequence of nucleated low-mass galaxies, and that envelope prominence correlates with distance from either galaxy. Our analysis provides evidence that tidal stripping of nucleated galaxies is an important process in the formation of UCDs.
  • The origin of ultra-compact dwarfs (UCDs)--objects larger and more massive than typical globular clusters (GCs), but more compact than typical dwarf galaxies--has been hotly debated in the 15 years since their discovery. Even whether UCDs should be considered galactic in origin, or simply the most extreme GCs, is not yet settled. We present the dynamical properties of 97 spectroscopically confirmed UCDs (rh >~10 pc) and 911 GCs associated with central cD galaxy of the Virgo cluster, M87. Our UCDs, of which 89% have M_star > ~2X10^6 M_sun and 92% are as blue as the classic blue GCs, nearly triple the sample of previous confirmed Virgo UCDs, providing by far the best opportunity for studying the global dynamics of a UCD system. We found that (1) UCDs have a surface number density profile that is shallower than that of the blue GCs in the inner ~ 70 kpc and as steep as that of the red GCs at larger radii; (2) UCDs exhibit a significantly stronger rotation than the GCs, and the blue GCs seem to have a velocity field that is more consistent with that of the surrounding dwarf ellipticals than with that of UCDs; (3) UCDs have a radially increasing orbital anisotropy profile, and are tangentially-biased at radii < ~ 40 kpc and radially-biased further out. In contrast, the blue GCs become more tangentially-biased at larger radii beyond ~ 40 kpc; (4) GCs with M_star > 2X10^6 M_sun have rotational properties indistinguishable from the less massive ones, suggesting that it is the size, instead of mass, that differentiates UCDs from GCs as kinematically distinct populations. We conclude that most UCDs in M87 are not consistent with being merely the most luminous and extended examples of otherwise normal GCs. The radially-biased orbital structure of UCDs at large radii is in general agreement with the "tidally threshed dwarf galaxy" scenario.
  • We report on a large-scale study of the distribution of globular clusters (GCs) throughout the Virgo cluster, based on photometry from the Next Generation Virgo Cluster Survey, a large imaging survey covering Virgo's primary subclusters to their virial radii. Using the g', (g'-i') color-magnitude diagram of unresolved and marginally-resolved sources, we constructed 2-D maps of the GC distribution. We present the clearest evidence to date showing the difference in concentration between red and blue GCs over the extent of the cluster, where the red (metal-rich) GCs are largely located around the massive early-type galaxies, whilst the blue (metal-poor) GCs have a more extended spatial distribution, with significant populations present beyond 83' (215 kpc) along the major axes of M49 and M87. The GC distribution around M87 and M49 shows remarkable agreement with the shape, ellipticity and boxiness of the diffuse light surrounding both galaxies. We find evidence for spatial enhancements of GCs surrounding M87 that may be indicative of recent interactions or an ongoing merger history. We compare the GC map to the locations of Virgo galaxies and the intracluster X-ray gas, and find good agreement between these baryonic structures. The Virgo cluster contains a total population of 67300$\pm$14400 GCs, of which 35% are located in M87 and M49 alone. We compute a cluster-wide specific frequency S_N,CL=$2.8\pm0.7$, including Virgo's diffuse light. The GC-to-baryonic mass fraction is e_b=$5.7\pm1.1\times10^{-4} $and the GC-to-total cluster mass formation efficiency is e_t=$2.9\pm0.5\times10^{-5}$, values slightly lower than, but consistent with, those derived for individual galactic halos. Our results show that the production of the complex structures in the unrelaxed Virgo cluster core (including the diffuse intracluster light) is an ongoing process.(abridged)
  • We study the dynamics of the giant elliptical galaxy M87 from the central to the outermost regions with the made-to-measure (M2M) method. We use a new catalogue of 922 globular cluster line-of- sight velocities extending to a projected radius of 180 kpc (equivalent to 25 M87 effective radii), and SAURON integral field unit data within the central 2.4 kpc. 263 globular clusters, mainly located beyond 40 kpc, are newly observed by the Next Generation Virgo Survey (NGVS). For the M2M modelling, the gravitational potential is taken as a combination of a luminous matter potential with a constant stellar mass-to-light ratio and a dark matter potential modelled as a logarithmic potential. Our best dynamical model returns a stellar mass-to-light ratio in the I band of M/LI = 6.0(+ -0.3) M_sun/L_sun with a dark matter potential scale velocity of 591(+ -50) km/s and scale radius of 42(+ -10) kpc. We determine the total mass of M87 within 180 kpc to be (1.5 + - 0.2) 10^13 M_sun. The mass within 40 kpc is smaller than previous estimates determined using globular cluster kinematics that did not extend beyond 45 kpc. With our new globular cluster velocities at much larger radii, we see that globular clusters around 40 kpc show an anomalously large velocity dispersion which affected previous results. The mass we derived is in good agreement with that inferred from ROSAT X-ray observation out to 180 kpc. Within 30 kpc our mass is also consistent with that inferred from Chandra and XMM-Newton X-ray observations, while within 120 kpc it is about 20% smaller. The model velocity dispersion anisotropy beta parameter for the globular clusters in M87 is small, varying from -0.2 at the centre to 0.2 at 40 kpc, and gradually decreasing to zero at 120 kpc.
  • The NGVS-IR project (Next Generation Virgo Survey - Infrared) is a contiguous near-infrared imaging survey of the Virgo cluster of galaxies. It complements the optical wide-field survey of Virgo (NGVS). The current state of NGVS-IR consists of Ks-band imaging of 4 deg^2 centered on M87, and J and Ks-band imaging of 16 deg^2 covering the region between M49 and M87. In this paper, we present the observations of the central 4 deg^2 centered on Virgo's core region. The data were acquired with WIRCam on the Canada-France-Hawaii Telescope and the total integration time was 41 hours distributed in 34 contiguous tiles. A survey-specific strategy was designed to account for extended galaxies while still measuring accurate sky brightness within the survey area. The average 5\sigma limiting magnitude is Ks=24.4 AB mag and the 50% completeness limit is Ks=23.75 AB mag for point source detections, when using only images with better than 0.7" seeing (median seeing 0.54"). Star clusters are marginally resolved in these image stacks, and Virgo galaxies with \mu_Ks=24.4 AB mag arcsec^-2 are detected. Combining the Ks data with optical and ultraviolet data, we build the uiK color-color diagram which allows a very clean color-based selection of globular clusters in Virgo. This diagnostic plot will provide reliable globular cluster candidates for spectroscopic follow-up campaigns needed to continue the exploration of Virgo's photometric and kinematic sub-structures, and will help the design of future searches for globular clusters in extragalactic systems. Equipped with this powerful new tool, future NGVS-IR investigations based on the uiK diagram will address the mapping and analysis of extended structures and compact stellar systems in and around Virgo galaxies.
  • The majority of clusters in the Universe have masses well below 10^5 Msun. Hence their integrated fluxes and colors can be affected by the random presence of a few bright stars introduced by stochastic sampling of the stellar mass function. Specific methods are being developed to extend the analysis of cluster SEDs into the low-mass regime. In this paper, we apply such a method to observations of star clusters, in the nearby spiral galaxy M83. We reassess ages and masses of a sample of 1242 objects for which UBVIHalpha fluxes were obtained with the HST/WFC3 images. Synthetic clusters with known properties are used to characterize the limitations of the method. The ensemble of color predictions of the discrete cluster models are in good agreement with the distribution of observed colors. We emphasize the important role of the Halpha data in the assessment of the fraction of young objects, particularly in breaking the age-extinction degeneracy that hampers an analysis based on UBVI only. We find the mass distribution of the cluster sample to follow a power-law of index -2.1 +/-0.2, and the distribution of ages a power-law of index -1.0 +/-0.2 for M > 10^3.5 Msun and ages between 10^7 and 10^9 yr. An extension of our main method, that makes full use of the probability distributions of age and mass of the individual clusters, is explored. It produces similar power-law slopes and will deserve further investigation. Although the properties derived for individual clusters significantly differ from those obtained with traditional, non-stochastic models in ~30% of the objects, the first order aspect of the age and mass distributions are similar to those obtained previously for this M83 sample in the range of overlap of the studies. We extend the power-law description to lower masses with better mass and age resolution and without most of the artifacts produced by the classical method.
  • The evolution of AGB stars is notoriously complex. The confrontation of AGB population models with observed stellar populations is a useful alternative to the detailed study of individual stars in efforts to converge towards a reliable evolution theory. I review here the impact of studies of star clusters on AGB models and AGB population synthesis, deliberately leaving out any more complex stellar populations. Over the last 10 years, despite much effort, the absolute uncertainties in the predictions of the light emitted by intermediate age populations have not been reduced to a satisfactory level. Observational sample definitions, as well as the combination of the natural variance in AGB properties with small number statistics, are largely responsible for this situation. There is hope that the constraints may soon become strong enough, thanks to large unbiased surveys of star clusters, resolved colour-magnitude diagrams, and new analysis methods that can account for the stochastic nature of AGB populations in clusters.
  • Star clusters are studied widely both as benchmarks for stellar evolution models and in their own right. Cluster age and mass distributions within galaxies are probes of star formation histories, and of cluster formation and disruption processes. The vast majority of clusters in the Universe is small, and it is well known that the integrated fluxes and colors have broad probability distributions, due to small numbers of bright stars. This paper goes beyond the description of predicted probability distributions, and presents results of the analysis of cluster energy distributions in an explicitly stochastic context. The method developed is Bayesian. It provides posterior probability distributions in the age-mass-extinction space, using multi-wavelength photometric observations and a large collection of Monte-Carlo simulations of clusters of finite stellar masses. Both UBVI and UBVIK datasets are considered, and the study conducted in this paper is restricted to the solar metallicity. We first reassess and explain errors arising from the use of standard analysis methods, which are based on continuous population synthesis models: systematic errors on ages and random errors on masses are large, while systematic errors on masses tend to be smaller. The age-mass distributions obtained after analysis of a synthetic sample are very similar to those found for real galaxies in the literature. The Bayesian approach on the other hand, is very successful in recovering the input ages and masses. Taking stochastic effects into account is important, more important for instance than the choice of adding or removing near-IR data in many cases. We found no immediately obvious reason to reject priors inspired by previous (standard) analyses of cluster populations in galaxies, i.e. cluster distributions that scale with mass as M^-2 and are uniform on a logarithmic age scale.
  • We point out a strong time-evolution of the mass-to-light conversion factor eta commonly used to estimate masses of unresolved star clusters from observed cluster spectro-photometric measures. We present a series of gas-dynamical models coupled with the Cambridge stellar evolution tracks to compute line-of-sight velocity dispersions and half-light radii weighted by the luminosity. We explore a range of initial conditions, varying in turn the cluster mass and/or density, and the stellar population's IMF. We find that eta, and hence the estimated cluster mass, may increase by factors as large as 3 over time-scales of 50 million years. We apply these results to an hypothetic cluster mass distribution function (d.f.) and show that the d.f. shape may be strongly affected at the low-mass end by this effect. Fitting truncated isothermal (Michie-King) models to the projected light profile leads to over-estimates of the concentration parameter c of delta c ~ 0.3 compared to the same functional fit applied to the projected mass density.