We use the IllustrisTNG cosmological hydrodynamical simulation to study the formation of ultra-diffuse galaxies (UDGs) in galaxy clusters. We supplement the simulations with a realistic mass-size relation for galaxies at the time of infall into the cluster, as well as an analytical model to describe the tidally-induced evolution of their stellar mass, velocity dispersion and size. The model assumes "cuspy" NFW halos and, contrary to recent claims, has no difficulty reproducing the observed number of UDGs in clusters. Our results further suggest that the UDG population consists of a mixture of "normal" low surface brightness galaxies such as those found in the field ("born" UDGs, or B-UDGs), as well as a distinct population that owe their large size and low surface brightness to the effects of cluster tides ("tidal", or T-UDGs). The simulations indicate that T-UDGs entered the cluster earlier and should be more prevalent than B-UDGs near the cluster centres. T-UDGs should also have, at given stellar mass, lower velocity dispersion, higher metallicities, and lower dark matter content than B-UDGs. Our results suggest that systems like DF-44 are consistent with having been born as UDGs, while others such as DF2, DF4 and VLSB-D are possibly extreme T-UDG examples.

We study the mass distribution and kinematics of the giant elliptical galaxy M87 (NGC 4486) using discrete chemo-dynamical, axisymmetric Jeans equation modelling. Our catalogue comprises 894 globular clusters (GCs) extending to a projected radius of $\sim 430$ kpc with line-of-sight velocities and colours, and Multi Unit Spectroscopic Explorer (MUSE) integral field unit data within the central $2.4$ kpc of the main galaxy. The gravitational potential for our models is a combination of a luminous matter potential with a varying mass-to-light ratio for the main galaxy, a supermassive black hole and a dark matter (DM) potential with a cusped or cored DM halo. The best-fitting models with either a cusped or a cored DM halo show no significant differences and both are acceptable. We obtain a total mass of $(2.16 \pm 0.38) \times 10^{13} M_{\odot}$ within $\sim$ 400 kpc. By including the stellar mass-to-light ratio gradient, the DM fraction increases from $\sim$ 26 percent (with no gradient) to $\sim$ 73 percent within $1\,R_e^{\rm maj}$ (major axis of half-light isophote, 14.2 kpc), and from $\sim$ 84 percent to $\sim$ 94 percent within $5\,R_e^{\rm maj}$ (71.2 kpc). Red GCs have moderate rotation with $V_{\rm max}/\sigma \sim$ 0.4, and blue GCs have weak rotation with $V_{\rm max}/\sigma \sim$ 0.1. Red GCs have tangential velocity dispersion anisotropy, while blue GCs are consistent with being nearly isotropic. Our results suggest that red GCs are more likely to be born in-situ, while blue GCs are more likely to be accreted.

It has long been speculated that many starburst or compact dwarf galaxies are resulted from dwarf-dwarf galaxy merging, but unequivocal evidence for this possibility has rarely been reported in the literature. We present the first study of deep optical broadband images of a gas-dominated blue compact dwarf galaxy (BCD) VCC848 (Mstar=2e8Msun) which hosts extended stellar shells and thus is confirmed to be a dwarf-dwarf merger. VCC848 is located in the outskirts of the Virgo Cluster. By analyzing the stellar light distribution, we found that VCC848 is the result of a merging between two dwarf galaxies with a primary-to-secondary mass ratio < ~ 5 for the stellar components and < ~ 2 for the presumed dark matter halos. The secondary progenitor galaxy has been almost entirely disrupted. The age-mass distribution of photometrically selected star cluster candidates in VCC848 implies that the cluster formation rate (CFR, proportional to star formation rate) was enhanced by a factor of ~ 7 - 10 during the past 1 Gyr. The merging-induced enhancement of CFR peaked near the galactic center a few hundred Myr ago and has started declining in the last few tens of Myr. The current star formation activities, as traced by the youngest clusters, mainly occur at large galactocentric distances (> ~ 1 kpc). The fact that VCC848 is still (atomic) gas-dominated after the period of most violent collision suggests that gas-rich dwarf galaxy merging can result in BCD-like remnants with extended atomic gas distribution surrounding a blue compact center, in general agreement with previous numerical simulations.

We study the assembly of globular clusters (GCs) in 9 galaxy clusters using the cosmological simulation Illustris. GCs are tagged to individual galaxies at their infall time. The tidal removal of GCs from their galaxies and the distribution of the GCs within the cluster is later followed self-consistently by the simulation. The method relies on the simple assumption of a single power-law relation between halo mass (M_vir) and mass in GCs (M_GC) as found in observations. We find that the GCs specific frequency $S_N$ as a function of V-band magnitude naturally reproduces the observed "U"-shape, due to the combination of a power law M_GC-M_vir relation and the non-linear M_*-M_vir relation from the simulation. Additional scatter in the $S_N$ values are traced back to galaxies with early infall times due to the evolution in the M_*-M_vir relation with redshift. GCs that have been tidally removed from their galaxies form today the intra-cluster component from which about ~60% were brought in by galaxies that orbit today within the cluster potential. The remaining "orphan" GCs are contributed by satellite galaxies with a wide range of stellar masses that are fully tidally disrupted at z=0. This intra-cluster component is a good dynamical tracer of the dark matter potential. As a consequence of the accreted nature of most intra-cluster GCs, their orbits are fairly radial with a predicted orbital anisotropy \beta >= 0.5. However, local tangential motions may appear as a consequence of localized substructure, providing a possible interpretation to the \beta<0 values suggested in observations of M87.

We present a detailed method to simulating sensor distortions using a photon and electron Monte Carlo method. We use three dimensional electrostatic simulations to parameterize the perturbed electric field profile for non-ideal sensor details. We follow the conversion of simulated photons, and the subsequent response of the converted electrons to the electric field pattern. These non-ideal sensor details can be implemented efficiently in a Monte Carlo approach. We demonstrate that the non-ideal sensor distortions have a variety of observable consequence including the modification of the astrometric pattern, the distortion of the electron diffusion size and shape, and the distortion of flats. We show analytic validation of the diffusion physics, reproduce two kinds of edge distortion, and show qualitative validation of field-free regions, lithography errors, and fringing. We also demonstrate that there are two related effects of doping variation having different observable consequences. We show that field distortions from accumulated electrons lead to intensity-dependent point-spread-functions and the sub-linear variance in flats. The method is implemented in the Photon Simulator (PhoSim) and the code is publically available.

We report observations of the hydrogen-deficient supernova (SN) 2019bkc/ATLAS19dqr. With B- and r-band decline between peak and 10 days post peak of Delta m_10(B)=5.24+/-0.07 mag and Delta m_10(r)=3.85+/-0.10$ mag, respectively, SN 2019bkc is the most rapidly declining SN I discovered so far. While its closest matches are the rapidly declining SN 2005ek and SN 2010X, the light curves and spectra of SN 2019bkc show some unprecedented characteristics. SN 2019bkc appears "hostless," with no identifiable host galaxy near its location, although it may be associated with the galaxy cluster MKW1 at z = 0.02. We evaluate a number of existing models of fast-evolving SNe, and we find that none of them can satisfactorily explain all aspects of SN 2019bkc observations.

We present the photometric properties of 2210 newly identified dwarf galaxy candidates in the MATLAS fields. The Mass Assembly of early Type gaLAxies with their fine Structures (MATLAS) deep imaging survey mapped $\sim$142 deg$^2$ of the sky around nearby isolated early type galaxies using MegaCam on the Canada-France-Hawaii Telescope, reaching surface brightnesses of $\sim$ 28.5 - 29 in the g-band. The dwarf candidates were identified through a direct visual inspection of the images and by visually cleaning a sample selected using a partially automated approach, and were morphologically classified at the time of identification. Approximately 75% of our candidates are dEs, indicating that a large number of early type dwarfs also populate low density environments, and 23.2% are nucleated. Distances were determined for 13.5% of our sample using pre-existing $z_{spec}$ measurements and HI detections. We confirm the dwarf nature for 99% of this sub-sample based on a magnitude cut $M_g$ = -18. Additionally, most of these ($\sim$90%) have relative velocities suggesting that they form a satellite population around nearby massive galaxies rather than an isolated field sample. Assuming that the candidates over the whole survey are satellites of the nearby galaxies, we demonstrate that the MATLAS dwarfs follow the same scaling relations as dwarfs in the Local Group as well as the Virgo and Fornax clusters. We also find that the nucleated fraction increases with $M_g$, and find evidence of a morphology-density relation for dwarfs around isolated massive galaxies.

The occurrence of planetary nebulae (PNe) in globular clusters (GCs) provides an excellent chance to study low-mass stellar evolution in a special (low-metallicity, high stellar density) environment. We report a systematic spectroscopic survey for the [O{\sc iii}] 5007 emission line of PNe in 1469 Virgo GCs and 121 Virgo ultra-compact dwarfs (UCDs), mainly hosted in the giant elliptical galaxies M87, M49, M86, and M84. We detected zero PNe in our UCD sample and discovered one PN ($M_{5007} = -4.1$ mag) associated with an M87 GC. We used the [O{\sc iii}] detection limit for each GC to estimate the luminosity-specific frequency of PNe, $\alpha$, and measured $\alpha$ in the Virgo cluster GCs to be $\alpha \sim 3.9_{-0.7}^{+5.2}\times 10^{-8}\mathrm{PN}/L_\odot$. $\alpha$ in Virgo GCs is among the lowest values reported in any environment, due in part to the large sample size, and is 5--6 times lower than that for the Galactic GCs. We suggest that $\alpha$ decreases towards brighter and more massive clusters, sharing a similar trend as the binary fraction, and the discrepancy between the Virgo and Galactic GCs can be explained by the observational bias in extragalactic surveys toward brighter GCs. This low but non-zero efficiency in forming PNe may highlight the important role played by binary interactions in forming PNe in GCs. We argue that a future survey of less massive Virgo GCs will be able to determine whether PN production in Virgo GCs is governed by internal process (mass, density, binary fraction), or is largely regulated by external environment.

The globular cluster (GC) specific frequency ($S_N$), defined as the number of GCs per unit galactic luminosity, represents the efficiency of GC formation (and survival) compared to field stars. Despite the naive expectation that star cluster formation should scale directly with star formation, this efficiency varies widely across galaxies. To explore this variation we measure the z-band GC specific frequency ($S_{N,z}$) for 43 early-type galaxies (ETGs) from the Hubble Space Telescope (HST)/Advanced Camera for Surveys (ACS) Fornax Cluster Survey. Combined with the homogenous measurements of $S_{N,z}$ in 100 ETGs from the HST/ACS Virgo Cluster Survey from Peng et al. (2008), we investigate the dependence of $S_{N,z}$ on mass and environment over a range of galaxy properties. We find that $S_{N,z}$ behaves similarly in the two galaxy clusters, despite the clusters' order-of-magnitude difference in mass density. The $S_{N,z}$ is low in intermediate-mass ETGs ($-20<M_z<-23$), and increases with galaxy luminosity. It is elevated at low masses, on average, but with a large scatter driven by galaxies in dense environments. The densest environments with the strongest tidal forces appear to strip the GC systems of low-mass galaxies. However, in low-mass galaxies that are not in strong tidal fields, denser environments correlate with enhanced GC formation efficiencies. Normalizing by inferred halo masses, the GC mass fraction, $\eta=(3.36\pm0.2)\times10^{-5}$, is constant for ETGs with stellar masses $\mathcal{M}_\star \lesssim 3\times10^{10}M_\odot$, in agreement with previous studies. The lack of correlation between the fraction of GCs and the nuclear light implies only a weak link between the infall of GCs and the formation of nuclei.

(Abridged) The Maunakea Spectroscopic Explorer (MSE) is an end-to-end science platform for the design, execution and scientific exploitation of spectroscopic surveys. It will unveil the composition and dynamics of the faint Universe and impact nearly every field of astrophysics across all spatial scales, from individual stars to the largest scale structures in the Universe. Major pillars in the science program for MSE include (i) the ultimate Gaia follow-up facility for understanding the chemistry and dynamics of the distant Milky Way, including the outer disk and faint stellar halo at high spectral resolution (ii) galaxy formation and evolution at cosmic noon, via the type of revolutionary surveys that have occurred in the nearby Universe, but now conducted at the peak of the star formation history of the Universe (iii) derivation of the mass of the neutrino and insights into inflationary physics through a cosmological redshift survey that probes a large volume of the Universe with a high galaxy density. MSE is positioned to become a critical hub in the emerging international network of front-line astronomical facilities, with scientific capabilities that naturally complement and extend the scientific power of Gaia, the Large Synoptic Survey Telescope, the Square Kilometer Array, Euclid, WFIRST, the 30m telescopes and many more.

The so-called ultra-diffuse galaxy NGC~1052-DF2 was announced to be a galaxy lacking dark matter based on a spectroscopic study of its constituent globular clusters. Here we present the first spectroscopic analysis of the stellar body of this galaxy using the MUSE integral-field spectrograph at the (ESO) Very Large Telescope. The MUSE datacube simultaneously provides DF2's stellar velocity field and systemic velocities for seven globular clusters (GCs). We further discovered three planetary nebulae (PNe) that are likely part of this galaxy. While five of the clusters had velocities measured in the literature, we were able to confirm the membership of two more candidates through precise radial velocity measurements, which increases the measured specific frequency of GCs in DF2. The mean velocity of the diffuse stellar body, 1792.9$^{-1.8}_{+1.4}$~\kms, is consistent with the mean globular cluster velocity. We detect a weak but significant velocity gradient within the stellar body, with a kinematic axis close to the photometric major-axis, making it a prolate-like rotator. We estimate a velocity dispersion from the clusters and PNe of $\sigma_{\mathrm{int}}=10.6^{+3.9}_{-2.3}$~\kms. The velocity dispersion $\sigma_{\rm{DF2}\star}$(\re) for the stellar body within one effective radius is $10.8^{-4.0}_{+3.2}$~\kms. Considering various sources of systemic uncertainties this central value varies between 5 and 13~\kms, and we conservatively report a 95\% confidence upper limit to the dispersion within one \re\ of 21~\kms. We provide updated mass estimates based on these dispersions corresponding to the different distances to NGC~1052-DF2 that have been reported in the recent literature.

We develop a comprehensive approach to simulate the deformation of mirrors and lenses due to thermal and mechanical stresses that couples efficiently to photon-based optics simulations. This expands upon previous work where we demonstrated a comprehensive ab initio approach to simulate astronomical images using a photon Monte Carlo method. We apply elasticity theory and estimate thermal effects by adapting a three-dimensional numerical method. We also consider the effect of active optics control systems and active cooling systems in further correcting distortions in the optics. We validate the approach by showing convergence to analytic estimates, and then apply the methodology to the WIYN 3.5m telescope primary mirror. We demonstrate that changes in the soak temperature result in second order point spread function (PSF) defocusing, the gravitational sag and positioning errors result in highly structured PSF distortions, and large-scale thermal gradients result in an elliptical PSF distortion patterns. All three aspects of the environment are larger than the intrinsic optical aberrations of the design, and further exploration with a variety of telescopes should lead to detailed PSF size and shape, astrometric distortion, and field variation predictions. The simulation capabilities developed in this work is publicly available with the Photon Simulation (PhoSim) package.

NGC 1052-DF2, an ultra diffuse galaxy (UDG), has been the subject of intense debate. Its alleged absence of dark matter, and the brightness and number excess of its globular clusters (GCs) at an initially assumed distance of 20Mpc, suggested a new formation channel for UDGs. We present the first systematic spectroscopic analysis of both the stellar body and the GCs (six of which were previously known, and one newly confirmed member) of this galaxy using MUSE@VLT. Even though NGC 1052-DF2 does not show any spatially extended emission lines we report the discovery of three planetary nebulae (PNe). We conduct full spectral fitting on the UDG and the stacked spectra of all GCs. The UDG's stellar population is old, 8.9$\pm$1.5 Gyr, metal-poor, with [M/H] = $-$1.07$\pm$0.12 with little or no $\alpha$-enrichment. The stacked spectrum of all GCs indicates a similar age of 8.9$\pm$1.8 Gyr, but lower metallicity, with [M/H] = $-$1.63$\pm$0.09, and similarly low $\alpha$-enrichment. There is no evidence for a variation of age and metallicity in the GC population with the available spectra. The significantly more metal-rich stellar body with respect to its associated GCs, the age of the population, its metallicity and alpha enrichment, are all in line with other dwarf galaxies. NGC 1052-DF2 thus falls on the same empirical mass-metallicity relation as other dwarfs, for the full distance range assumed in the literature. We find that both debated distance estimates (13 and 20 Mpc) are similarly likely, given the three discovered PNe.

We study the intrinsic shapes of a sample of over 400 quiescent galaxies in the cores of the Virgo and Fornax clusters with luminosities $10^{6} \leq L_{g}/L_{\odot} \leq 10^{8}$. Similar to satellites of the Local Group and Centaurus A, these faint, low surface brightness cluster galaxies are best described as a family of thick ($C/A > 0.5$), oblate-triaxial spheroids. However, the large sample size allows us to show that the flattening of their stellar distributions depends both on luminosity and on the presence of a nuclear star cluster. Nucleated satellites are thicker at all luminosities compared to their non-nucleated counterparts, and fainter galaxies are systematically thicker as well, regardless of nucleation. Once nucleation is accounted for, we find no evidence that the environment the satellites live in plays a relevant role in setting their three-dimensional structure. We interpret both the presence of stellar nuclei and the associated thicker shapes as the result of preferential early and rapid formation, effectively making these faint nucleated galaxies the first generation of cluster satellites.

We have observed the giant elliptical galaxy M87 during the Virgo Environmental Survey Tracing Galaxy Evolution (VESTIGE), a blind narrow-band Halpha+[NII] imaging survey of the Virgo cluster carried out with MegaCam at the Canada French Hawaii Telescope (CFHT). The deep narrow-band image confirmed the presence of a filament of ionised gas extending up to ~ 3 kpc in the north-western direction and ~ 8 kpc to the south-east, with a couple of plumes of ionised gas, the weakest of which, at ~ 18 kpc from the nucleus, was previously unknown. The analysis of deep optical images taken from the NGVS survey confirms that this gas filament is associated with dust seen in absorption which is now detected up to ~ 2.4 kpc from the nucleus. We have also analysed the physical and kinematical properties of the ionised gas filament using deep IFU MUSE data covering the central 4.8 x 4.8 kpc^2 of the galaxy. The spectroscopic data confirms a perturbed kinematics of the ionised gas, with differences in velocity of ~ 700-800 km s^-1 on scales of < 1 kpc. The analysis of 2D diagnostic diagrams and the observed relationship between the shock-sensitive [OI]/Halpha line ratio and the velocity dispersion of the gas suggest that the gas is shock-ionised.

Using deep, high resolution optical imaging from the Next Generation Virgo Cluster Survey we study the properties of nuclear star clusters (NSCs) in a sample of nearly 400 quiescent galaxies in the core of Virgo with stellar masses $10^{5}\lesssim M_{*}/M_{\odot} \lesssim10^{12}$. The nucleation fraction reaches a peak value $f_{n}\approx90\%$ for $M_{*} \approx 10^{9} M_{\odot}$ galaxies and declines for both higher and lower masses, but nuclei populate galaxies as small as $M_{*} \approx5\times10^{5} M_{\odot}$. Comparison with literature data for nearby groups and clusters shows that at the low-mass end nucleation is more frequent in denser environments. The NSC mass function peaks at $M_{NSC}\approx7\times10^{5} M_{\odot}$, a factor 3-4 times larger than the turnover mass for globular clusters (GCs). We find a nonlinear relation between the stellar masses of NSCs and of their host galaxies, with a mean nucleus-to-galaxy mass ratio that drops to $M_{NSC}/M_{*}\approx3.6\times10^{-3}$ for $M_{*} \approx 5\times10^{9} M_{\odot}$ galaxies. Nuclei in both more and less massive galaxies are much more prominent: $M_{NSC}\propto M_{*}^{0.46}$ at the low-mass end, where nuclei are nearly 50% as massive as their hosts. We measure an intrinsic scatter in NSC masses at fixed galaxy stellar mass of 0.4 dex, which we interpret as evidence that the process of NSC growth is significantly stochastic. At low galaxy masses we find a close connection between NSCs and GC systems, including a very similar occupation distribution and comparable total masses. We discuss these results in the context of current dissipative and dissipationless models of NSC formation.

We explore a scenario where metal poor globular clusters (GCs) form at the centres of their own dark matter halos in the early universe before reionization. This hypothesis leads to predictions about the abundance, distribution and kinematics of GCs today that we explore using cosmological N-body simulations and analytical modelling. We find that selecting the massive tail of collapsed objects at $z\gtrsim 9$ as GC formation sites leads to four main predictions: i) a highly clustered population of GCs around galaxies today, ii) a natural scaling between number of GCs and halo virial mass that follows roughly the observed trend, iii) a very low number of free floating GCs outside massive halos and iv) GCs should be embedded within massive and extended dark matter (sub)halos. We find that the strongest constraint to the model is given by the combination of (i) and (ii): a mass cut to tagged GCs halos which accounts for the number density of metal poor GCs today predicts a radial distribution that is too extended compared to recent observations. On the other hand, a mass cut sufficient to match the observed half number radius could only explain 60% of the metal poor population. In all cases, observations favour early redshifts for GC formation ($z\geq 15$).

Intra-cluster (IC) populations are expected to be a natural result of the hierarchical assembly of clusters, yet their low space densities make them difficult to detect and study. We present the first definitive kinematic detection of an IC population of globular clusters (GCs) in the Virgo cluster, around the central galaxy, M87. This study focuses on the Virgo core for which the combination of NGVS photometry and follow-up spectroscopy allows us to reject foreground star contamination and explore GC kinematics over the full Virgo dynamical range. The GC kinematics changes gradually with galactocentric distance, decreasing in mean velocity and increasing in velocity dispersion, eventually becoming indistinguishable from the kinematics of Virgo dwarf galaxies at $\mathrm{R>320\, kpc}$. By kinematically tagging M87 halo and intra-cluster GCs we find that 1) the M87 halo has a smaller fraction ($52\pm3\%$) of blue clusters with respect to the IC counterpart ($77\pm10\%$), 2) the $(g'-r')_{0}$ vs $(i'-z')_{0}$ color-color diagrams reveal a galaxy population that is redder than the IC population that may be due to a different composition in chemical abundance and progenitor mass, and 3) the ICGC distribution is shallower and more extended than the M87 GCs, yet still centrally concentrated. The ICGC specific frequency, $S_{N,\mathrm{ICL}}=10.2\pm4.8$, is consistent with what is observed for the population of quenched, low-mass galaxies within 1~Mpc from the cluster's center. The IC population at Virgo's center is thus consistent with being an accreted component from low-mass galaxies tidally stripped or disrupted through interactions, with a total mass of $\mathrm{M_{ICL,tot}=10.8\pm0.1\times10^{11}M_{\odot}}$.

Ultra-diffuse galaxies (UDGs) are unusual galaxies with low luminosities, similar to classical dwarf galaxies, but sizes up to $\sim\!5$ larger than expected for their mass. Some UDGs have large populations of globular clusters (GCs), something unexpected in galaxies with such low stellar density and mass. We have carried out a comprehensive study of GCs in both UDGs and classical dwarf galaxies at comparable stellar masses using HST observations of the Coma cluster. We present new imaging for 33 Dragonfly UDGs with the largest effective radii ($>2$ kpc), and additionally include 15 UDGs and 54 classical dwarf galaxies from the HST/ACS Coma Treasury Survey and the literature. Out of a total of 48 UDGs, 27 have statistically significant GC systems, and 11 have candidate nuclear star clusters. The GC specific frequency ($S_N$) varies dramatically, with the mean $S_N$ being higher for UDGs than for classical dwarfs. At constant stellar mass, galaxies with larger sizes (or lower surface brightnesses) have higher $S_N$, with the trend being stronger at higher stellar mass. At lower stellar masses, UDGs tend to have higher $S_N$ when closer to the center of the cluster, i.e., in denser environments. The fraction of UDGs with a nuclear star cluster also depends on environment, varying from $\sim\!40$\% in the cluster core, where it is slightly lower than the nucleation fraction of classical dwarfs, to $\lesssim20\%$ in the outskirts. Collectively, we observe an unmistakable diversity in the abundance of GCs, and this may point to multiple formation routes.

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.

We describe a program to measure surface brightness fluctuation (SBF) distances to galaxies observed in the Next Generation Virgo Cluster Survey (NGVS), a photometric imaging survey covering $104~deg^2$ of the Virgo cluster in the ${u}^*,g,i,z$ bandpasses with the Canada-France Hawaii Telescope. We describe the selection of the sample galaxies, the procedures for measuring the apparent $i$-band SBF magnitude $\bar{i}$, and the calibration of the absolute $\bar{M}_i$ as a function of observed stellar population properties. The multi-band NGVS data set provides multiple options for calibrating the SBF distances, and we explore various calibrations involving individual color indices as well as combinations of two different colors. Within the color range of the present sample, the two-color calibrations do not significantly improve the scatter with respect to wide-baseline, single-color calibrations involving $u^{*}$. We adopt the ${u}^*{-}z$ calibration as reference for the present galaxy sample, with an observed scatter of 0.11 mag. For a few cases that lack good ${u}^*$ photometry, we use an alternative relation based on a combination of $g{-}i$ and $g{-}z$ colors, with only a slightly larger observed scatter of 0.12 mag. The agreement of our measurements with the best existing distance estimates provides confidence that our measurements are accurate. We present a preliminary catalog of distances for 89 galaxies brighter than $B_T\approx13.0$ mag within the survey footprint, including members of the background M and W Clouds at roughly twice the distance of the main body of the Virgo cluster. The extension of the present work to fainter and bluer galaxies is in progress.

The Virgo Environmental Survey Tracing Ionised Gas Emission (VESTIGE) is a blind narrow-band H$\alpha$+[NII] imaging survey carried out with MegaCam at the CFHT. During pilot observations we have observed NGC 4330, an intermediate mass, edge-on star forming spiral currently falling into the core of the Virgo cluster. New deep observations revealed a low surface brightness 10 kpc tail exhibiting a peculiar filamentary structure. The filaments are remarkably parallel one another and clearly indicate the direction of motion of the galaxy in the Virgo potential. Motivated by the detection of these features, indicating ongoing gas stripping, we collected literature photometry in 15 bands from the far-UV to the far-IR and deep optical long slit spectroscopy using the FORS2 instrument at the ESO Very Large Telescope. Using a newly developed Monte Carlo code that jointly fits spectroscopy and photometry, we reconstructed the star formation histories in apertures along the major axis of the galaxy. Our results have been validated against the output of CIGALE, a fitting code which has been previously used for similar studies. We found a clear outside-in gradient with radius of the time when the quenching event started: the outermost radii have been stripped 500 Myr ago, while the stripping has reached the inner 5 kpc from the center in the last 100 Myr. Regions at even smaller radii are currently still forming stars fueled by the presence of HI and H2 gas. When compared to statistical studies of the quenching timescales in the local Universe we find that ram pressure stripping of the cold gas is an effective mechanism to reduce the transformation times for galaxies falling into massive clusters. Future systematic studies of all the active galaxies observed by VESTIGE in the Virgo cluster will extend these results to a robust statistical framework.

The Virgo Environmental Survey Tracing Ionised Gas Emission (VESTIGE) is a blind narrow-band Halpha+[NII] imaging survey carried out with MegaCam at the Canada-France-Hawaii Telescope. The survey covers the whole Virgo cluster region from its core to one virial radius (104 deg^2). The sensitivity of the survey is of f(Halpha) ~ 4 x 10^-17 erg sec-1 cm^-2 (5 sigma detection limit) for point sources and Sigma (Halpha) ~ 2 x 10^-18 erg sec^-1 cm^-2 arcsec^-2 (1 sigma detection limit at 3 arcsec resolution) for extended sources, making VESTIGE the deepest and largest blind narrow-band survey of a nearby cluster. This paper presents the survey in all its technical aspects, including the survey design, the observing strategy, the achieved sensitivity in both the narrow-band Halpha+[NII] and in the broad-band r filter used for the stellar continuum subtraction, the data reduction, calibration, and products, as well as its status after the first observing semester. We briefly describe the Halpha properties of galaxies located in a 4x1 deg^2 strip in the core of the cluster north of M87, where several extended tails of ionised gas are detected. This paper also lists the main scientific motivations of VESTIGE, which include the study of the effects of the environment on galaxy evolution, the fate of the stripped gas in cluster objects, the star formation process in nearby galaxies of different type and stellar mass, the determination of the Halpha luminosity function and of the Halpha scaling relations down to ~ 10^6 Mo stellar mass objects, and the reconstruction of the dynamical structure of the Virgo cluster. This unique set of data will also be used to study the HII luminosity function in hundreds of galaxies, the diffuse Halpha+[NII] emission of the Milky Way at high Galactic latitude, and the properties of emission line galaxies at high redshift.

Future weak lensing surveys potentially hold the highest statistical power for constraining cosmological parameters compared to other cosmological probes. The statistical power of a weak lensing survey is determined by the sky coverage, the inverse of the noise in shear measurements, and the galaxy number density. The combination of the latter two factors is often expressed in terms of $n_{\rm eff}$ -- the "effective number density of galaxies used for weak lensing measurements". In this work, we estimate $n_{\rm eff}$ for the Large Synoptic Survey Telescope (LSST) project, the most powerful ground-based lensing survey planned for the next two decades. We investigate how the following factors affect the resulting $n_{\rm eff}$ of the survey with detailed simulations: (1) survey time, (2) shear measurement algorithm, (3) algorithm for combining multiple exposures, (4) inclusion of data from multiple filter bands, (5) redshift distribution of the galaxies, and (6) masking and blending. For the first time, we quantify in a general weak lensing analysis pipeline the sensitivity of $n_{\rm eff}$ to the above factors. We find that with current weak lensing algorithms, expected distributions of observing parameters, and all lensing data ($r$- and $i$-band, covering 18,000 degree$^{2}$ of sky) for LSST, $n_{\rm eff} \approx37$ arcmin$^{-2}$ before considering blending and masking, $n_{\rm eff} \approx31$ arcmin$^{-2}$ when rejecting seriously blended galaxies and $n_{\rm eff} \approx26$ arcmin$^{-2}$ when considering an additional 15% loss of galaxies due to masking. With future improvements in weak lensing algorithms, these values could be expected to increase by up to 20%. Throughout the paper, we also stress the ways in which $n_{\rm eff}$ depends on our ability to understand and control systematic effects in the measurements.

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.

We present a survey of X-ray point sources in the nearest and dynamically young galaxy cluster, Virgo, using archival Chandra observations that sample the vicinity of 80 early-type member galaxies. The X-ray source populations at the outskirt of these galaxies are of particular interest. We detect a total of 1046 point sources (excluding galactic nuclei) out to a projected galactocentric radius of $\sim$40 kpc and down to a limiting 0.5-8 keV luminosity of $\sim$$2\times10^{38}{\rm~erg~s^{-1}}$. Based on the cumulative spatial and flux distributions of these sources, we statistically identify $\sim$120 excess sources that are not associated with the main stellar content of the individual galaxies, nor with the cosmic X-ray background. This excess is significant at a 3.5 $\sigma$ level, when Poisson error and cosmic variance are taken into account. On the other hand, no significant excess sources are found at the outskirt of a control sample of field galaxies, suggesting that at least some fraction of the excess sources around the Virgo galaxies are truly intra-cluster X-ray sources. Assisted with ground-based and HST optical imaging of Virgo, we discuss the origins of these intra-cluster X-ray sources, in terms of supernova-kicked low-mass X-ray binaries (LMXBs), globular clusters, LMXBs associated with the diffuse intra-cluster light, stripped nucleated dwarf galaxies and free-floating massive black holes.

Small distortions in the images of Einstein rings or giant arcs offer the exciting prospect of detecting dark matter haloes or subhaloes of mass below $10^9$M$_{\odot}$, most of which are too small to have made a visible galaxy. A very large number of such haloes are predicted to exist in the cold dark matter model of cosmogony; in contrast other models, such as warm dark matter, predict no haloes below a mass of this order which depends on the properties of the warm dark matter particle. Attempting to detect these small perturbers could therefore discriminate between different kinds of dark matter particles, and even rule out the cold dark matter model altogether. Globular clusters in the lens galaxy also induce distortions in the image which could, in principle, contaminate the test. Here, we investigate the population of globular clusters in six early type galaxies in the Virgo cluster. We find that the number density of globular clusters of mass $\sim10^6$M$_{\odot}$ is comparable to that of the dark matter perturbers (including subhaloes in the lens and haloes along the line-of-sight). We show that the very different degrees of mass concentration in globular clusters and dark matter haloes result in different lensing distortions. These are detectable with milli-arcsecond resolution imaging which can distinguish between globular cluster and dark matter halo signals.

Using a cosmological $N$-body simulation, we investigate the origin and distribution of stars in the intracluster light (ICL) of a Fornax-like cluster. In a dark matter only simulation we identify a halo which, at $z=0$, has $M_200 \simeq 4.1 \times 10^{13}M_{sun}$ and $r_{200} = 700kpc$, and replace infalling subhalos with models that include spheroid and disc components. As they fall into the cluster, the stars in some of these galaxies are stripped from their hosts, and form the ICL. We consider the separate contributions to the ICL from stars which originate in the haloes and the discs of the galaxies. We find that disc ICL stars are more centrally concentrated than halo ICL stars. The majority of the disc ICL stars are associated with one initially disc-dominated galaxy that falls to the centre of the cluster and is heavily disrupted, producing part of the cD galaxy. At radial distances greater than 200kpc, well beyond the stellar envelope of the cD galaxy, stars formerly from the stellar haloes of galaxies dominate the ICL. Therefore at large distances, the ICL population is dominated by older stars.

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.

We present here parameter distributions and fundamental scaling relations for 190 galaxies as part of the Spectroscopy and H-bang Imaging of Virgo cluster galaxies (SHIVir) survey. We find the distribution of galaxy velocities to be bimodal about $V_{\rm circ} \sim 125$ km ${\rm s^{-1}}$, hinting at the existence of dynamically unstable modes in the inner regions of galaxies. An analysis of the Tully-Fisher relation (TFR) of late-type galaxies (LTGs) and fundamental plane (FP) of early-type galaxies (ETGs) is also presented, yielding a compendium of galaxy scaling relations. The slope and zero-point of the Virgo TFR match those of field galaxies, while scatter differences likely reflect distinct evolutionary histories. The velocities minimizing scatter for the TFR and FP are measured at large apertures where the baryonic fraction becomes subdominant. While TFR residuals remain independent of any galaxy parameters, FP residuals (i.e. the FP "tilt") correlate strongly with the dynamical-to-stellar mass ratio, yielding stringent galaxy formation constraints. Furthermore, we construct a stellar-to-total mass relation (STMR) for ETGs and LTGs and find linear but distinct trends over the range $M_{*} = 10^{8-11} M_{\odot}$. Stellar-to-halo mass relations (SHMRs), which probe the extended dark matter halo, can be scaled down to masses estimated within the optical radius, showing a tight match with the Virgo STMR at low masses; however, possibly inadequate halo abundance matching prescriptions and broad radial scalings complicate this comparison at all masses. While ETGs appear to be more compact than LTGs of the same stellar mass in projected space, their mass-size relations in physical space are identical. The trends reported here call for validation through well-resolved numerical simulations.

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.

The Next Generation Virgo Cluster Survey is a deep (with a $2\sigma$ detection limit $\mu_g$ = 29~mag~arcsec$^{-2}$ in the $g-$band) optical panchromatic survey targeting the Virgo cluster from its core to virial radius, for a total areal coverage of 104 square degrees. As such, the survey is well suited for the study of galaxies' outskirts, haloes and low surface brightness features that arise from dynamical interactions within the cluster environment. We report the discovery of extremely faint ($\mu_g$ $>$ 25 mag arcsec$^{-2}$) shells in three Virgo cluster early-type dwarf galaxies, VCC~1361, VCC~1447 and VCC~1668. Among them, VCC~1447 has an absolute magnitude M$_{g}$ = -11.71 mag and is {\it the least massive galaxy with a shell system discovered to date}. We present a detailed study of these low surface brightness features. We detect between three and four shells in each of our galaxies. Within the uncertainties, we find no evidence of a color difference between the galaxy main body and shell features. The observed arcs of the shells are located upto several effective radii of the galaxies. We further explore the origin of these low surface brightness features with the help of idealized numerical simulations. We find that a near equal mass merger is best able to reproduce the main properties of the shells, including their quite symmetric appearance and their alignment along the major axis of the galaxy. The simulations provide support for a formation scenario in which a recent merger, between two near-equal mass, gas-free dwarf galaxies forms the observed shell systems.

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.

The mean alpha-to-iron abundance ratio ([$\alpha$/Fe]) of galaxies is sensitive to the chemical evolution processes at early time, and it is an indicator of star formation timescale ($\tau_{{\rm SF}}$). Although the physical reason remains ambiguous, there is a tight relation between [$\alpha$/Fe] and stellar velocity dispersion ($\sigma$) among massive early-type galaxies (ETGs). However, no work has shown convincing results as to how this relation behaves at low masses. We assemble 15 data sets from the literature and build a large sample that includes 192 nearby low-mass ($18<\sigma<80$~\kms) ETGs. We find that the [$\alpha$/Fe]-$\sigma$ relation generally holds for low-mass ETGs, except in extreme environments. Specifically, in normal galaxy cluster environments, the [$\alpha$/Fe]-$\sigma$ relation and its intrinsic scatter are, within uncertainties, similar for low-mass and high-mass ETGs. However, in the most massive relaxed galaxy cluster in our sample, the zero point of the relation is higher and the intrinsic scatter is significantly larger. By contrast, in galaxy groups the zero point of the relation offsets in the opposite direction, again with substantial intrinsic scatter. The elevated [$\alpha$/Fe] of low-mass ETGs in the densest environments suggests that their star formation was quenched earlier than in high-mass ETGs. For the low-mass ETGs in the lowest density environments, we suggest that their more extended star formation histories suppressed their average [$\alpha$/Fe]. The large scatter in [$\alpha$/Fe] may reflect stochasticity in the chemical evolution of low-mass galaxies.

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.

Diffuse star clusters (DSCs) are old and dynamically hot stellar systems that have lower surface brightness and more extended morphology than globular clusters (GCs). Using the images from HST/ACS Fornax Cluster Survey, we find that 12 out of 43 early-type galaxies (ETGs) in the Fornax cluster host significant numbers of DSCs. Together with literature data from the HST/ACS Virgo Cluster Survey, where 18 out of 100 ETGs were found to host DSCs, we systematically study the relationship of DSCs with GCs, and their host galaxy environment. Two DSC hosts are post-merger galaxies, with most of the other hosts either having low mass or showing clear disk components. We find that while the number ratio of DSCs to GCs is nearly constant in massive galaxies, the DSC-to-GC ratio becomes systematically higher in lower mass hosts. This suggests that DSCs may be more efficient at forming (or surviving) in low density environments. DSC hosts are not special either in their position in the cluster, or in the galactic color-magnitude diagram. Why some disk and low-mass galaxies host DSCs while others do not is still a puzzle, however. The mean ages of DSC hosts and non-hosts are similar at similar masses, implying that formation efficiency, rather than survival, is the reason behind different DSC number fractions in early-type galaxies.

MSE is an 11.25m aperture observatory with a 1.5 square degree field of view that will be fully dedicated to multi-object spectroscopy. More than 3200 fibres will feed spectrographs operating at low (R ~ 2000 - 3500) and moderate (R ~ 6000) spectral resolution, and approximately 1000 fibers will feed spectrographs operating at high (R ~ 40000) resolution. MSE is designed to enable transformational science in areas as diverse as tomographic mapping of the interstellar and intergalactic media; the in-situ chemical tagging of thick disk and halo stars; connecting galaxies to their large scale structure; measuring the mass functions of cold dark matter sub-halos in galaxy and cluster-scale hosts; reverberation mapping of supermassive black holes in quasars; next generation cosmological surveys using redshift space distortions and peculiar velocities. MSE is an essential follow-up facility to current and next generations of multi-wavelength imaging surveys, including LSST, Gaia, Euclid, WFIRST, PLATO, and the SKA, and is designed to complement and go beyond the science goals of other planned and current spectroscopic capabilities like VISTA/4MOST, WHT/WEAVE, AAT/HERMES and Subaru/PFS. It is an ideal feeder facility for E-ELT, TMT and GMT, and provides the missing link between wide field imaging and small field precision astronomy. MSE is optimized for high throughput, high signal-to-noise observations of the faintest sources in the Universe with high quality calibration and stability being ensured through the dedicated operational mode of the observatory. (abridged)

Observations of nearby galaxy clusters at low surface brightness have identified galaxies with low luminosities, but sizes as large as L* galaxies, leading them to be dubbed "ultra-diffuse galaxies" (UDGs). The survival of UDGs in dense environments like the Coma cluster suggests that UDGs could reside in much more massive dark halos. We report the detection of a substantial population of globular clusters (GCs) around a Coma UDG, Dragonfly 17 (DF17). We find that DF17 has a high GC specific frequency of S_N=26+/-13. The GC system is extended, with an effective radius of 12"+/-2", or 5.6+/-0.9 kpc at Coma distance, 70% larger than the galaxy itself. We also estimate the mean of the GC luminosity function to infer a distance of 97 (+17/-14) Mpc, providing redshift-independent confirmation that one of these UDGs is in the Coma cluster. The presence of a rich GC system in DF17 indicates that, despite its low stellar density, star formation was intense enough to form many massive star clusters. If DF17's ratio of total GC mass to total halo mass is similar to those in other galaxies, then DF17 has an inferred total mass of ~10^11 solar masses, only ~10% the mass of the Milky Way, but extremely dominated by dark matter, with M/L_V~1000. We suggest that UDGs like DF17 may be "pure stellar halos", i.e., galaxies that formed their stellar halo components, but then suffered an early cessation in star formation that prevented the formation of any substantial central disk or bulge.

We present new HST optical and near-infrared (NIR) photometry of the rich globular cluster (GC) system of NGC 4874, the cD galaxy in the core of the Coma cluster (Abell 1656). NGC 4874 was observed with the HST Advanced Camera for Surveys in the F475W (g) and F814W (I) passbands and the Wide Field Camera 3 IR Channel in F160W (H). The GCs in this field exhibit a bimodal optical color distribution with more than half of the GCs falling on the red side at g-I > 1. Bimodality is also present, though less conspicuously, in the optical-NIR I-H color. Consistent with past work, we find evidence for nonlinearity in the g-I versus I-H color-color relation. Our results thus underscore the need for understanding the detailed form of the color-metallicity relations in interpreting observational data on GC bimodality. We also find a very strong color-magnitude trend, or "blue tilt," for the blue component of the optical color distribution of the NGC 4874 GC system. A similarly strong trend is present for the overall mean I-H color as a function of magnitude; for M_814 < -10 mag, these trends imply a steep mass-metallicity scaling with $Z\propto M_{\rm GC}^{1.4\pm0.4}$, but the scaling is not a simple power law and becomes much weaker at lower masses. As in other similar systems, the spatial distribution of the blue GCs is more extended than that of the red GCs, partly because of blue GCs associated with surrounding cluster galaxies. In addition, the center of the GC system is displaced by 4+/-1 kpc towards the southwest from the luminosity center of NGC 4874, in the direction of NGC 4872. Finally, we remark on a dwarf elliptical galaxy with a noticeably asymmetrical GC distribution. Interestingly, this dwarf has a velocity of nearly -3000 km/s with respect to NGC 4874; we suggest it is on its first infall into the cluster core and is undergoing stripping of its GC system by the cluster potential.

We analyze the kinematics of six Virgo cluster dwarf early-type galaxies (dEs) from their globular cluster (GC) systems. We present new Keck/DEIMOS spectroscopy for three of them and reanalyze the data found in the literature for the remaining three. We use two independent methods to estimate the rotation amplitude (Vmax) and velocity dispersion (sigma_GC) of the GC systems and evaluate their statistical significance by simulating non-rotating GC systems with the same number of GC satellites and velocity uncertainties. Our measured kinematics agree with the published values for the three galaxies from the literature and, in all cases, some rotation is measured. However, our simulations show that the null hypothesis of being non-rotating GC systems cannot be ruled out. In the case of VCC1861, the measured Vmax and the simulations indicate that it is not rotating. In the case of VCC1528, the null hypothesis can be marginally ruled out, thus, it might be rotating although further confirmation is needed. In our analysis, we find that, in general, the measured Vmax tends to be overestimated and the measured sigma_GC tends to be underestimated by amounts that depend on the intrinsic Vmax/sigma_GC, the number of observed GCs (N_GC), and the velocity uncertainties. The bias is negligible when N_GC>~20. In those cases where a large N_GC is not available, it is imperative to obtain data with small velocity uncertainties. For instance, errors of <2km/s lead to Vmax<10km/s for a system that is intrinsically not rotating.

(Abridged) We investigate the intrinsic shapes of low-luminosity galaxies in the central 300 kpc of the Virgo cluster using deep imaging obtained as part of the NGVS. We build a sample of nearly 300 red-sequence cluster members in the yet unexplored $-14 < M_{g} < -8$ magnitude range. The observed distribution of apparent axis ratios is then fit by families of triaxial models with normally-distributed intrinsic ellipticities and triaxialities. We develop a Bayesian framework to explore the posterior distribution of the model parameters, which allows us to work directly on discrete data, and to account for individual, surface brightness-dependent axis ratio uncertainties. For this population we infer a mean intrinsic ellipticity E=0.43, and a mean triaxiality T=0.16. This implies that faint Virgo galaxies are best described as a family of thick, nearly oblate spheroids with mean intrinsic axis ratios 1:0.94:0.57. We additionally attempt a study of the intrinsic shapes of Local Group satellites of similar luminosities. For the LG population we infer a slightly larger mean intrinsic ellipticity E=0.51, and the paucity of objects with round apparent shapes translates into more triaxial mean shapes, 1:0.76:0.49. We finally compare the intrinsic shapes of NGVS low-mass galaxies with samples of more massive quiescent systems, and with field, star-forming galaxies of similar luminosities. We find that the intrinsic flattening in this low-luminosity regime is almost independent of the environment in which the galaxy resides--but there is a hint that objects may be slightly rounder in denser environments. The comparable flattening distributions of low-luminosity galaxies that have experienced very different degrees of environmental effects suggests that internal processes are the main drivers of galaxy structure at low masses--with external mechanisms playing a secondary role.

The Next Generation Virgo Cluster Survey (NGVS) is a deep $u^*giz$ survey targeting the Virgo cluster of galaxies at 16.5~Mpc. This survey provides high-quality photometry over an $\sim$ 100 deg$^2$ region straddling the constellations of Virgo and Coma Berenices. This sightline through the Milky Way is noteworthy in that it intersects two of the most prominent substructures in the Galactic halo: the Virgo Over-Density (VOD) and Sagittarius stellar stream (close to its bifurcation point). In this paper, we use deep $u^*gi$ imaging from the NGVS to perform tomography of the VOD and Sagittarius stream using main-sequence turnoff (MSTO) stars as a halo tracer population. The VOD, whose centroid is known to lie at somewhat lower declinations ($\alpha \sim 190^\circ$, $\delta \sim -5^\circ$) than is covered by the NGVS, is nevertheless clearly detected in the NGVS footprint at distances between $\sim$ 8 and 25~kpc. By contrast, the Sagittarius stream is found to slice directly across the NGVS field at distances between 25 and 40~kpc, with a density maximum at $\simeq$~35~kpc. No evidence is found for new substructures beyond the Sagittarius stream, at least out to a distance of $\sim$~90 kpc --- the largest distance to which we can reliably trace the halo using MSTO stars. We find clear evidence for a distance gradient in the Sagittarius stream across the $\sim 30$~deg of sky covered by the NGVS and its flanking fields. We compare our distance measurements along the stream to those predicted by leading stream models.

We explore the environmental dependence of star formation timescales in low mass galaxies using the [$\alpha$/Fe] abundance ratio as an evolutionary clock. We present integrated [$\alpha$/Fe] measurements for 11 low mass ($M_\star \sim 10^9~M_\odot$) early-type galaxies (ETGs) with a large range of cluster-centric distance in the Virgo Cluster. We find a gradient in [$\alpha$/Fe], where the galaxies closest to the cluster center (the cD galaxy, M87) have the highest values. This trend is driven by galaxies within a projected radius of 0.4~Mpc (0.26 times the virial radius of Virgo~A), all of which have super-solar [$\alpha$/Fe]. Galaxies in this mass range exhibit a large scatter in the [$\alpha$/Fe]--$\sigma$ diagram, and do not obviously lie on an extension of the relation defined by massive ETGs. In addition, we find a correlation between [$\alpha$/Fe] and globular cluster specific frequency ($S_N$), suggesting that low-mass ETGs that formed their stars over a short period of time, were also efficient at forming massive star clusters. The innermost low-mass ETGs in our sample have [$\alpha$/Fe] values comparable to that of M87, implying that environment is the controlling factor for star formation timescales in dense regions. These low-mass galaxies could be the surviving counterparts of the objects that have already been accreted into the halo of M87, and may be the link between present-day low-mass galaxies and the old, metal-poor, high-[$\alpha$/Fe], high-$S_N$ stellar populations seen in the outer halos of massive ETGs.

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.

We report on the properties of the most massive ultra-compact dwarf galaxy (UCD) in the nearby Virgo Cluster of galaxies using imaging from the Next Generation Virgo Cluster Survey (NGVS) and spectroscopy from Keck/DEIMOS. This object (M59-UCD3) appears to be associated with the massive Virgo galaxy M59 (NGC 4621), has an integrated velocity dispersion of 78 km/s, a dynamical mass of $3.7\times10^8 M_\odot$, and an effective radius ($R_e$) of 25 pc. With an effective surface mass density of $9.4\times10^{10} M_\odot/kpc^2$, it is the densest galaxy in the local Universe discovered to date, surpassing the density of the luminous Virgo UCD, M60-UCD1. M59-UCD3 has a total luminosity of $M_{g'}=-14.2$ mag, and a spectral energy distribution consistent with an old (14 Gyr) stellar population with [Fe/H]=0.0 and [$\alpha$/Fe]=+0.2. We also examine deep imaging around M59 and find a broad low surface brightness stream pointing towards M59-UCD3, which may represent a tidal remnant of the UCD progenitor. This UCD, along with similar objects like M60-UCD1 and M59cO, likely represents an extreme population of tidally stripped galaxies more akin to larger and more massive compact early-type galaxies than to nuclear star clusters in present-day dwarf galaxies.

We investigate a scenario in which feedback from black-hole X-ray binaries (BHXBs) sometimes begins inside young star clusters before strong supernova feedback. Those BHXBs could reduce the gas fraction inside embedded young clusters whilst maintaining virial equilibrium, which may help globular clusters (GCs) to stay bound when supernova-driven gas ejection subsequently occurs. Adopting a simple toy model with parameters guided by BHXB population models, we produce GC formation efficiencies consistent with empirically-inferred values. The metallicity dependence of BHXB formation could naturally explain why GC formation efficiency is higher at lower metallicity. For reasonable assumptions about that metallicity dependence, our toy model can produce a GC metallicity bimodality in some galaxies without a bimodality in the field-star metallicity distribution.

We report the discovery of three large (R29 >~ 1 arcminute) extremely low surface brightness (mu_(V,0) ~ 27.0) galaxies identified using our deep, wide-field imaging of the Virgo Cluster from the Burrell Schmidt telescope. Complementary data from the Next Generation Virgo Cluster Survey do not resolve red giant branch stars in these objects down to i=24, yielding a lower distance limit of 2.5 Mpc. At the Virgo distance, these objects have half-light radii 3-10 kpc and luminosities L_V=2-9x10^7 Lsun. These galaxies are comparable in size but lower in surface brightness than the large ultradiffuse LSB galaxies recently identified in the Coma cluster, and are located well within Virgo's virial radius; two are projected directly on the cluster core. One object appears to be a nucleated LSB in the process of being tidally stripped to form a new Virgo ultracompact dwarf galaxy. The others show no sign of tidal disruption, despite the fact that such objects should be most vulnerable to tidal destruction in the cluster environment. The relative proximity of Virgo makes these objects amenable to detailed studies of their structural properties and stellar populations. They thus provide an important new window onto the connection between cluster environment and galaxy evolution at the extremes.

The Next Generation Virgo Cluster Survey has recently determined the luminosity function of galaxies in the core of the Virgo cluster down to unprecedented magnitude and surface brightness limits. Comparing simulations of cluster formation to the derived central stellar mass function, we attempt to estimate the stellar-to-halo-mass ratio (SHMR) for dwarf galaxies, as it would have been before they fell into the cluster. This approach ignores several details and complications, e.g., the contribution of ongoing star formation to the present-day stellar mass of cluster members, and the effects of adiabatic contraction and/or violent feedback on the subhalo and cluster potentials. The final results are startlingly simple, however; we find that the trends in the SHMR determined previously for bright galaxies appear to extend down in a scale-invariant way to the faintest objects detected in the survey. These results extend measurements of the formation efficiency of field galaxies by two decades in halo mass, or five decades in stellar mass, down to some of the least massive dwarf galaxies known, with stellar masses of $\sim 10^5 M_\odot$.

We present a kinematic analysis of the globular cluster systems and diffuse stellar light of four intermediate luminosity (sub-$L^{\ast}$) early-type galaxies in the Virgo cluster based on Gemini/GMOS data. Our galaxy sample is fainter ($-23.8<M_K<-22.7$) than most previous studies, nearly doubling the number of galaxies in this magnitude range that now have GC kinematics. The data for the diffuse light extends to $4R_e$, and the data for the globular clusters reaches 8--$12R_e$. We find that the kinematics in these outer regions are all different despite the fact that these four galaxies have similar photometric properties, and are uniformly classified as "fast rotators" from their stellar kinematics within $1R_e$. The globular cluster systems exhibit a wide range of kinematic morphology. The rotation axis and amplitude can change between the inner and outer regions, including a case of counter-rotation. This difference shows the importance of wide-field kinematic studies, and shows that stellar and GC kinematics can change significantly as one moves beyond the inner regions of galaxies. Moreover, the kinematics of the globular cluster systems can differ from that of the stars, suggesting that the formation of the two populations are also distinct.

We present a comprehensive methodology for the simulation of astronomical images from optical survey telescopes. We use a photon Monte Carlo approach to construct images by sampling photons from models of astronomical source populations, and then simulating those photons through the system as they interact with the atmosphere, telescope, and camera. We demonstrate that all physical effects for optical light that determine the shapes, locations, and brightnesses of individual stars and galaxies can be accurately represented in this formalism. By using large scale grid computing, modern processors, and an efficient implementation that can produce 400,000 photons/second, we demonstrate that even very large optical surveys can be now be simulated. We demonstrate that we are able to: 1) construct kilometer scale phase screens necessary for wide-field telescopes, 2) reproduce atmospheric point-spread-function moments using a fast novel hybrid geometric/Fourier technique for non-diffraction limited telescopes, 3) accurately reproduce the expected spot diagrams for complex aspheric optical designs, and 4) recover system effective area predicted from analytic photometry integrals. This new code, the photon simulator (PhoSim), is publicly available. We have implemented the Large Synoptic Survey Telescope (LSST) design, and it can be extended to other telescopes. We expect that because of the comprehensive physics implemented in PhoSim, it will be used by the community to plan future observations, interpret detailed existing observations, and quantify systematics related to various astronomical measurements. Future development and validation by comparisons with real data will continue to improve the fidelity and usability of the code.

We present Gemini GMOS-IFU data of eight compact low-mass early-type galaxies (ETGs) in the Virgo cluster. We analyse their stellar kinematics, stellar population, and present two-dimensional maps of these properties covering the central 5"x 7" region. We find a large variety of kinematics: from non- to highly-rotating objects, often associated with underlying disky isophotes revealed by deep images from the Next Generation Virgo Cluster Survey. In half of our objects, we find a centrally-concentrated younger and more metal-rich stellar population. We analyze the specific stellar angular momentum through the lambdaR parameter and find six fast-rotators and two slow-rotators, one having a thin counter-rotating disk. We compare the local galaxy density and stellar populations of our objects with those of 39 more extended low-mass Virgo ETGs from the SMAKCED survey and 260 massive ($M>10^{10}$\Msun) ETGs from the A3D sample. The compact low-mass ETGs in our sample are located in high density regions, often close to a massive galaxy and have, on average, older and more metal-rich stellar populations than less compact low-mass galaxies. We find that the stellar population parameters follow lines of constant velocity dispersion in the mass-size plane, smoothly extending the comparable trends found for massive ETGs. Our study supports a scenario where low-mass compact ETGs have experienced long-lived interactions with their environment, including ram-pressure stripping and gravitational tidal forces, that may be responsible for their compact nature.

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 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.

The Next Generation Virgo Cluster Survey is an optical imaging survey covering 104 deg^2 centered on the Virgo cluster. Currently, the complete survey area has been observed in the u*giz-bands and one third in the r-band. We present the photometric redshift estimation for the NGVS background sources. After a dedicated data reduction, we perform accurate photometry, with special attention to precise color measurements through point spread function-homogenization. We then estimate the photometric redshifts with the Le Phare and BPZ codes. We add a new prior which extends to iAB = 12.5 mag. When using the u*griz-bands, our photometric redshifts for 15.5 \le i \lesssim 23 mag or zphot \lesssim 1 galaxies have a bias |\Delta z| < 0.02, less than 5% outliers, and a scatter \sigma_{outl.rej.} and an individual error on zphot that increase with magnitude (from 0.02 to 0.05 and from 0.03 to 0.10, respectively). When using the u*giz-bands over the same magnitude and redshift range, the lack of the r-band increases the uncertainties in the 0.3 \lesssim zphot \lesssim 0.8 range (-0.05 < \Delta z < -0.02, \sigma_{outl.rej} ~ 0.06, 10-15% outliers, and zphot.err. ~ 0.15). We also present a joint analysis of the photometric redshift accuracy as a function of redshift and magnitude. We assess the quality of our photometric redshifts by comparison to spectroscopic samples and by verifying that the angular auto- and cross-correlation function w(\theta) of the entire NGVS photometric redshift sample across redshift bins is in agreement with the expectations.

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.

In this paper we introduce the deepest and most extensive ultraviolet extragalactic source catalogs of the Virgo Cluster area to date. Archival and targeted GALEX imaging is compiled and combined to provide the deepest possible coverage over ~120 deg^2 in the NUV (lambda_eff=2316 angstroms) and ~40 deg^2 in the FUV (lambda_eff=1539 angstroms) between 180 deg <= R.A. <= 195 deg and 0 deg <= Decl. <= 20 deg. We measure the integrated photometry of 1770 extended UV sources of all galaxy types and use GALEX pipeline photometry for 1,230,855 point-like sources in the foreground, within, and behind the cluster. Extended source magnitudes are reliable to m_UV ~22, showing ~0.01 sigma difference from their asymptotic magnitudes. Point-like source magnitudes have a 1 sigma standard deviation within ~0.2 mag down to m_uv ~23. The point-like source catalog is cross-matched with large optical databases and surveys including the SDSS DR9 (> 1 million Virgo Cluster sources), the Next Generation Virgo Cluster Survey (NGVS; >13 million Virgo Cluster sources), and the NED (~30,000 sources in the Virgo Cluster). We find 69% of the entire UV point-like source catalog has a unique optical counterpart, 11% of which are stars and 129 are Virgo cluster members neither in the VCC nor part of the bright CGCG galaxy catalog (i.e., m_pg < 14.5). These data are collected in four catalogs containing the UV extended sources, the UV point-like sources, and two catalogs each containing the most relevant optical parameters of UV-optically matched point-like sources for further studies from SDSS and NGVS. The GUViCS catalogs provide a unique set of data for future works on UV and multiwavelength studies in the cluster and background environments.

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.

We use our integrated SDSS photometry for 96 globular clusters in $g$ and $z$, as well as $r$ and $i$ photometry for a subset of 56 clusters, to derive the integrated colour-metallicity relation (CMR) for Galactic globular clusters. We compare this relation to previous work, including extragalactic clusters, and examine the influence of age, present-day mass function variations, structural parameters and the morphology of the horizontal branch on the relation. Moreover, we scrutinise the scatter introduced by foreground extinction (including differential reddening) and show that the scatter in the colour-metallicity relation can be significantly reduced combining two reddening laws from the literature. In all CMRs we find some low-reddening young GCs that are offset to the CMR. Most of these outliers are associated with the Sagittarius system. Simulations show that this is due less to age than to a different enrichment history. Finally, we introduce colour-metallicity relations based on the infrared Calcium triplet, which are clearly non-linear when compared to $(g^\prime-i^\prime)$ and $(g^\prime-z^\prime)$ colours.

We present $g^\prime$ and $z^\prime$ aperture photometry for 96 Galactic Globular Clusters, making this the largest homogeneous catalog of photometry for these objects in the SDSS filter system. For a subset of 56 clusters we also provide photometry in $r^\prime$ and $i^\prime$. We carry out comparisons with previous photometry as well as with the SDSS dataset. The data will be useful for a series of applications in Galactic and extragalactic astrophysics. Future papers will analyse the colour-metallicity relation, colour-magnitude diagrams, and structural parameters. The compilation of results based on this dataset will be collected in the Galactic Globular Cluster Catalog (G2C2).

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.

We study the azimuthal distribution of globular clusters (GCs) in early-type galaxies and compare them to their host galaxies using data from the ACS Virgo Cluster Survey. We find that in host galaxies with visible elongation (epsilon > 0.2) and intermediate to high luminosities (M_z<-19), the GCs are preferentially aligned along the major axis of the stellar light. The red (metal-rich) GC subpopulations show strong alignment with the major axis of the host galaxy, which supports the notion that these GCs are associated with metal-rich field stars. The metal-rich GCs in lenticular galaxies show signs of being more strongly associated with disks rather than bulges. Surprisingly, we find that the blue (metal-poor) GCs can also show the same correlation. If the metal-poor GCs are part of the early formation of the halo and built up through mergers, then our results support a picture where halo formation and merging occur anisotropically, and where the present day major axis is an indicator of the preferred merging axis.

Recent observations have shown that compact nuclear star clusters (NSCs) are present in up to 80% of galaxies. However, detailed studies of their dynamical and chemical properties are confined mainly to spiral galaxy hosts, where they are more easily observed. In this paper we present our study of the NSC in FCC 277, a nucleated elliptical galaxy in the Fornax cluster. We use a combination of adaptive optics assisted near-infrared integral field spectroscopy, Hubble Space Telescope imaging, and literature long slit data. We show that while the NSC does not appear to rotate within our detection limit of ~6 km/s, rotation is detected at larger radii, where the isophotes appear to be disky, suggesting the presence of a nuclear disk. We also observe a distinct central velocity dispersion drop that is indicative of a dynamically cold rotating sub-system. Following the results of orbit-based dynamical modelling, co-rotating as well as counter-rotating stellar orbits are simultaneously needed to reproduce the observed kinematics. We find evidence for varying stellar populations, with the NSC and nuclear disk hosting younger and more metal rich stars than the main body of the galaxy. We argue that gas dissipation and some level of merging have likely played an important role in the formation of the nucleus of this intermediate-mass galaxy. This is in contrast to NSCs in low-mass early- type galaxies, which may have been formed primarily through the infall of star clusters.

The ability to measure metallicities and {\alpha}-element abundances in individual red giant branch (RGB) stars using medium-resolution spectra ($R \approx 6000$) is a valuable tool for deciphering the nature of Milky Way dwarf satellites and the history of the Galactic halo. Extending such studies to more distant systems like Andromeda is beyond the ability of the current generation of telescopes, but by co-adding the spectra of similar stars, we can attain the necessary signal-to-noise ratio to make detailed abundance measurements. In this paper, we present a method to determine metallicities and {\alpha}-element abundances using the co-addition of medium resolution spectra. We test the method of spectral co-addition using high-S/N spectra of more than 1300 RGB stars from Milky Way globular clusters and dwarf spheroidal galaxies obtained with the Keck II telescope/DEIMOS spectrograph. We group similar stars using photometric criteria and compare the weighted ensemble average abundances ([Fe/H], [Mg/Fe], [Si/Fe], [Ca/Fe] and [Ti/Fe]) of individual stars in each group with the measurements made on the corresponding co-added spectrum. We find a high level of agreement between the two methods, which permits us to apply this co-added spectra technique to more distant RGB stars, like stars in the M31 satellite galaxies. This paper outlines our spectral co-addition and abundance measurement methodology and describes the potential biases in making these measurements.

The statistics of peak counts in reconstructed shear maps contain information beyond the power spectrum, and can improve cosmological constraints from measurements of the power spectrum alone if systematic errors can be controlled. We study the effect of galaxy shape measurement errors on predicted cosmological constraints from the statistics of shear peak counts with the Large Synoptic Survey Telescope (LSST). We use the LSST image simulator in combination with cosmological N-body simulations to model realistic shear maps for different cosmological models. We include both galaxy shape noise and, for the first time, measurement errors on galaxy shapes. We find that the measurement errors considered have relatively little impact on the constraining power of shear peak counts for LSST.

A main science goal for the Large Synoptic Survey Telescope (LSST) is to measure the cosmic shear signal from weak lensing to extreme accuracy. One difficulty, however, is that with the short exposure time ($\simeq$15 seconds) proposed, the spatial variation of the Point Spread Function (PSF) shapes may be dominated by the atmosphere, in addition to optics errors. While optics errors mainly cause the PSF to vary on angular scales similar or larger than a single CCD sensor, the atmosphere generates stochastic structures on a wide range of angular scales. It thus becomes a challenge to infer the multi-scale, complex atmospheric PSF patterns by interpolating the sparsely sampled stars in the field. In this paper we present a new method, PSFent, for interpolating the PSF shape parameters, based on reconstructing underlying shape parameter maps with a multi-scale maximum entropy algorithm. We demonstrate, using images from the LSST Photon Simulator, the performance of our approach relative to a 5th-order polynomial fit (representing the current standard) and a simple boxcar smoothing technique. Quantitatively, PSFent predicts more accurate PSF models in all scenarios and the residual PSF errors are spatially less correlated. This improvement in PSF interpolation leads to a factor of 3.5 lower systematic errors in the shear power spectrum on scales smaller than $\sim13'$, compared to polynomial fitting. We estimate that with PSFent and for stellar densities greater than $\simeq1/{\rm arcmin}^{2}$, the spurious shear correlation from PSF interpolation, after combining a complete 10-year dataset from LSST, is lower than the corresponding statistical uncertainties on the cosmic shear power spectrum, even under a conservative scenario.

The complete 10-year survey from the Large Synoptic Survey Telescope (LSST) will image $\sim$ 20,000 square degrees of sky in six filter bands every few nights, bringing the final survey depth to $r\sim27.5$, with over 4 billion well measured galaxies. To take full advantage of this unprecedented statistical power, the systematic errors associated with weak lensing measurements need to be controlled to a level similar to the statistical errors. This work is the first attempt to quantitatively estimate the absolute level and statistical properties of the systematic errors on weak lensing shear measurements due to the most important physical effects in the LSST system via high fidelity ray-tracing simulations. We identify and isolate the different sources of algorithm-independent, \textit{additive} systematic errors on shear measurements for LSST and predict their impact on the final cosmic shear measurements using conventional weak lensing analysis techniques. We find that the main source of the errors comes from an inability to adequately characterise the atmospheric point spread function (PSF) due to its high frequency spatial variation on angular scales smaller than $\sim10'$ in the single short exposures, which propagates into a spurious shear correlation function at the $10^{-4}$--$10^{-3}$ level on these scales. With the large multi-epoch dataset that will be acquired by LSST, the stochastic errors average out, bringing the final spurious shear correlation function to a level very close to the statistical errors. Our results imply that the cosmological constraints from LSST will not be severely limited by these algorithm-independent, additive systematic effects.

The Advanced Camera for Surveys (ACS) Fornax Cluster Survey is a Hubble Space Telescope program to image 43 early-type galaxies in the Fornax cluster, using the F475W and F850LP bandpasses of the ACS. We employ both 1D and 2D techniques to characterize the properties of the stellar nuclei in these galaxies, defined as the central "luminosity excesses" relative to a Sersic model fitted to the underlying host. We find 72+/-13% of our sample (31 galaxies) to be nucleated, with only three of the nuclei offset by more than 0.5" from their galaxy photocenter, and with the majority of nuclei having colors bluer than their hosts. The nuclei are observed to be larger, and brighter, than typical Fornax globular clusters, and to follow different structural scaling relations. A comparison of our results to those from the ACS Virgo Cluster Survey reveals striking similarities in the properties of the nuclei belonging to these different environments. We briefly review a variety of proposed formation models and conclude that, for the low-mass galaxies in our sample, the most important mechanism for nucleus growth is probably infall of star clusters through dynamical friction, while for higher mass galaxies, gas accretion triggered by mergers, accretions and tidal torques is likely to dominate, with the relative importance of these two processes varying smoothly as a function of galaxy mass. Some intermediate-mass galaxies in our sample show a complexity in their inner structure that may be the signature of "hybrid nuclei" that arose through parallel formation channels.

Context: We study the peculiar interacting galaxy system of VCC1249/M49 located in the core of the Virgo B subcluster. Owing to a recent interaction between the dwarf galaxy VCC1249 and the halo gas of the gE M49, neutral hydrogen has been displaced from the interstellar medium of this dwarf into the Virgo ICM. Observations also reveal multiple compact star-forming regions that are embedded in this HI cloud, with a projected separation up to 13 kpc from VCC1249 in the northwest direction. Aims: Motivated by recent NUV imaging from GUViCS of the VCC1249/M49 system that shows significant ongoing/recent star formation in the compact regions, we aim to constrain the origin of these outlying HII regions with a multi-wavelength approach. Methods: Using deep optical (u, g, i, z) imaging from NGVS and new Halpha imaging obtained at the San Pedro Martir observatory together with Keck long-slit spectroscopy, we characterize the SFR, ages, and metallicity of VCC1249 and its outlying compact regions. Moreover, we analyze the color and luminosity profile of the galaxy to investigate its recent interaction with M49. Results: Our new observations indicate that VCC1249 underwent a recent interaction with M49 in which both ram-pressure stripping and tidal interaction occured. The joint action of the two mechanisms led to the removal of the HI gas from the ISM of VCC1249, while the gravitational tides triggered the stellar tail and counter-tail of VCC1249. Our SED analysis reveals that the star formation in this galaxy was truncated around 200 Myr ago and that the outlying HII regions were born in situ about 10 Myr ago out of pre-enriched gas removed from the dwarf galaxy. These observations also reveal that interactions between central and satellite galaxies similar to the one between VCC1249/M49 may be an effective way of dispersing metals into the halos of massive galaxies.

(ABRIDGED) We use high resolution (~0.1") F814W ACS images from the HST ACS Treasury survey of the Coma cluster at z~0.02 to study bars in massive disk galaxies (S0s), and in dwarf galaxies in the Coma core. Our study helps constrain the evolution of bars and disks in dense environments and provides a comparison point for studies in lower density environments and at higher redshifts. (1) We characterize the fraction and properties of bars in a sample of 32 bright (M_V <= -18, M_* > 10^9.5 M_sun) S0 galaxies, which dominate the population of massive disk galaxies in the Coma core. Measuring the S0 bar fraction must be handled carefully, as the results depend on the method used: the bar fraction for bright S0s in the Coma core is 50%+/-11%, 65%+/-11%, and 60%+/-11% for three methods of bar detection: strict ellipse fit criteria, relaxed ellipse fit criteria, and visual classification. (2) We compare the S0 bar fraction across different environments (Coma core, A901/902, Virgo). We find that the bar fraction among bright S0 galaxies does not show a statistically significant variation across environments spanning two orders of magnitude in galaxy number density (n~300-10,000 gal/Mpc^3). We speculate that the S0 bar fraction is not significantly enhanced in rich clusters because S0s in rich clusters are less prone to bar instabilities as they are dynamically hot and gas poor due to ram pressure stripping and accelerated star formation. In addition, high-speed encounters in rich clusters may be less effective than slow, strong encounters in inducing bars. (3) We analyze a sample of 333 faint (M_V > -18) dwarf galaxies in the Coma core. Using unsharp-masking, we find only 13 galaxies with bar and/or spiral structure. The paucity of disk structures in Coma dwarfs suggests that either disks are not common in these galaxies, or that any disks present are too hot to develop instabilities.

We combine new Wide Field Camera~3 IR Channel (WFC3/IR) F160W (H) imaging data for NGC1399, the central galaxy in the Fornax cluster, with archival F475W (g), F606W (V), F814W (I), and F850LP (z) optical data from the Advanced Camera for Surveys (ACS). The purely optical g-I, V-I, and g-z colors of NGC1399's rich globular cluster (GC) system exhibit clear bimodality, at least for magnitudes $I_814 > 21.5$. The optical-IR I-H color distribution appears unimodal, and this impression is confirmed by mixture modeling analysis. The V-H colors show marginal evidence for bimodality, consistent with bimodality in V-I and unimodality in I-H. If bimodality is imposed for I-H with a double Gaussian model, the preferred blue/red split differs from that for optical colors; these "differing bimodalities" mean that the optical and optical-IR colors cannot both be linearly proportional to metallicity. Consistent with the differing color distributions, the dependence of I-H on g-I for the matched GC sample is significantly nonlinear, with an inflection point near the trough in the g-I color distribution; the result is similar for the I-H dependence on g-z colors taken from the ACS Fornax Cluster Survey. These g-z colors have been calibrated empirically against metallicity; applying this calibration yields a continuous, skewed, but single-peaked metallicity distribution. Taken together, these results indicate that nonlinear color-metallicity relations play an important role in shaping the observed bimodal distributions of optical colors in extragalactic GC systems.

One of the conundrums in extragalactic astronomy is the discrepancy in observed metallicity distribution functions (MDFs) between the two prime stellar components of early-type galaxies-globular clusters (GCs) and halo field stars. This is generally taken as evidence of highly decoupled evolutionary histories between GC systems and their parent galaxies. Here we show, however, that new developments in linking the observed GC colors to their intrinsic metallicities suggest nonlinear color-to-metallicity conversions, which translate observed color distributions into strongly-peaked, unimodal MDFs with broad metal-poor tails. Remarkably, the inferred GC MDFs are similar to the MDFs of resolved field stars in nearby elliptical galaxies and those produced by chemical evolution models of galaxies. The GC MDF shape, characterized by a sharp peak with a metal-poor tail, indicates a virtually continuous chemical enrichment with a relatively short timescale. The characteristic shape emerges across three orders of magnitude in the host galaxy mass, suggesting a universal process of chemical enrichment among various GC systems. Given that GCs are bluer than field stars within the same galaxy, it is plausible that the chemical enrichment processes of GCs ceased somewhat earlier than that of field stellar population, and if so, GCs preferentially trace the major, vigorous mode of star formation events in galactic formation. We further suggest a possible systematic age difference among GC systems, in that the GC systems in more luminous galaxies are older. This is consistent with the downsizing paradigm of galaxies and supports additionally the similar nature shared by GCs and field stars. Our findings suggest that GC systems and their parent galaxies have shared a more common origin than previously thought, and hence greatly simplify theories of galaxy formation.

We have undertaken a spectroscopic search for ultra compact dwarf galaxies (UCDs) in the dense core of the dynamically evolved, massive Coma cluster as part of the HST/ACS Coma Cluster Treasury Survey. UCD candidates were initially chosen based on color, magnitude, degree of resolution within the ACS images, and the known properties of Fornax and Virgo UCDs. Follow-up spectroscopy with Keck/LRIS confirmed 27 candidates as members of the Coma Cluster, a success rate > 60% for targeted objects brighter than M_R = -12. Another 14 candidates may also prove to be Coma members, but low signal-to-noise spectra prevent definitive conclusions. An investigation of the properties and distribution of the Coma UCDs finds these objects to be very similar to UCDs discovered in other environments. The Coma UCDs tend to be clustered around giant galaxies in the cluster core and have colors/metallicity that correlate with the host galaxy. With properties and a distribution similar to that of the Coma cluster globular cluster population, we find strong support for a star cluster origin for the majority of the Coma UCDs. However, a few UCDs appear to have stellar population or structural properties which differentiate them from the old star cluster populations found in the Coma cluster, perhaps indicating that UCDs may form through multiple formation channels.

Intracluster stellar populations are a natural result of tidal interactions in galaxy clusters. Measuring these populations is difficult, but important for understanding the assembly of the most massive galaxies. The Coma cluster is one of the nearest truly massive galaxy clusters, and is host to a correspondingly large system of globular clusters (GCs). We use imaging from the HST/ACS Coma Cluster Survey to present the first definitive detection of a large population of intracluster GCs (IGCs) that fills the Coma cluster core and is not associated with individual galaxies. The GC surface density profile around the central massive elliptical galaxy, NGC 4874, is dominated at large radii by a population of IGCs that extend to the limit of our data (R<520 kpc). We estimate that there are 47000+/-1600 (random) +4000/-5000 (systematic) IGCs out to this radius, and that they make up ~70% of the central GC system, making this the largest GC system in the nearby Universe. Even including the GC systems of other cluster galaxies, IGCs still make up ~30-45% of the GCs in the cluster core. Observational limits from previous studies of the intracluster light (ICL) suggest that the IGC population has a high specific frequency. If the IGC population has a specific frequency similar to high-S_N dwarf galaxies, then the ICL has a total stellar mass of ~10^12 M_sun within the cluster core. The ICL makes up approximately half of the stellar luminosity and one-third of the stellar mass of the central (NGC4874+ICL) system. The color distribution of the IGC population is bimodal, with blue, metal-poor GCs outnumbering red, metal-rich GCs by a ratio of 4:1. The fraction of red IGCs (20%), and the red color of those GCs, implies that IGCs can originate from the halos of relatively massive, L* galaxies, and not solely from the disruption of dwarf galaxies. (Abridged)

We use the largest homogeneous sample of globular clusters (GCs), drawn from the ACSVCS and ACSFCS, to investigate the color gradients of GC systems in 76 early-type galaxies. We find that most GC systems possess an obvious negative gradient in g-z color (bluer outwards). For GC systems displaying color bimodality, both metal-rich and metal-poor GC subpopulations present shallower but significant color gradients on average, and the mean gradients of these two subpopulations are of roughly equal strength. The FOV of ACS mainly restricts us to measuring the inner gradients of GC systems. These gradients, however, can introduce an aperture bias when measuring the mean colors of GC subpopulations from relatively narrow central pointings. Inferred corrections to previous work imply a reduced significance for the relation between the mean color of metal-poor GCs and their host galaxy luminosity. The GC color gradients also show a dependence with host galaxy mass where the gradiens are weakest at the ends of the mass spectrum--in massive galaxies and dwarf galaxies--and strongest in galaxies of intermediate mass, around a stellar mass of M_stellar~10^10M_sun. We also measure color gradients for field stars in the host galaxies. We find that GC color gradients are systematically steeper than field star color gradients, but the shape of the gradient-mass relation is the same for both. If gradients are caused by rapid dissipational collapse and weakened by merging, these color gradients support a picture where the inner GC systems of most intermediate-mass and massive galaxies formed early and rapidly with the most massive galaxies having experienced greater merging. The lack of strong gradients in the GC systems of dwarfs, which probably have not experienced many recent major mergers, suggests that low mass halos were inefficient at retaining and mixing metals during the epoch of GC formation.

Although early observations with the Hubble Space Telescope (HST) pointed to a sharp dichotomy among early-type galaxies in terms of the logarithmic slope gamma' of their central surface brightness profiles, several studies in the past few years have called this finding into question. In particular, recent imaging surveys of 143 early-type galaxies belonging to the Virgo and Fornax Clusters using the Advanced Camera for Surveys (ACS) on board HST have not found a dichotomy in gamma', but instead a systematic progression from central luminosity deficit to excess relative to the inward extrapolation of the best-fitting global Sersic model. Given that earlier studies also found that the dichotomy persisted when analyzing the deprojected density profile slopes, we investigate the distribution of the three-dimensional luminosity density profiles of the ACS Virgo and Fornax Cluster Survey galaxies. Having fitted the surface brightness profiles with modified Sersic models, we then deproject the galaxies using an Abel integral and measure the inner slopes gamma_3D of the resulting luminosity density profiles at various fractions of the effective radius R_e. We find no evidence of a dichotomy, but rather, a continuous variation in the central luminosity profiles as a function of galaxy magnitude. We introduce a parameter, Delta_3D, that measures the central deviation of the deprojected luminosity profiles from the global Sersic fit, showing that this parameter varies smoothly and systematically along the luminosity function.

The Coma cluster was the target of a HST-ACS Treasury program designed for deep imaging in the F475W and F814W passbands. Although our survey was interrupted by the ACS instrument failure in 2007, the partially completed survey still covers ~50% of the core high-density region in Coma. Observations were performed for 25 fields that extend over a wide range of cluster-centric radii (~1.75 Mpc) with a total coverage area of 274 arcmin^2. The majority of the fields are located near the core region of Coma (19/25 pointings) with six additional fields in the south-west region of the cluster. In this paper we present reprocessed images and SExtractor source catalogs for our survey fields, including a detailed description of the methodology used for object detection and photometry, the subtraction of bright galaxies to measure faint underlying objects, and the use of simulations to assess the photometric accuracy and completeness of our catalogs. We also use simulations to perform aperture corrections for the SExtractor Kron magnitudes based only on the measured source flux and half-light radius. We have performed photometry for ~73,000 unique objects; one-half of our detections are brighter than the 10-sigma point-source detection limit at F814W=25.8 mag (AB). The slight majority of objects (60%) are unresolved or only marginally resolved by ACS. We estimate that Coma members are 5-10% of all source detections, which consist of a large population of unresolved objects (primarily GCs but also UCDs) and a wide variety of extended galaxies from a cD galaxy to dwarf LSB galaxies. The red sequence of Coma member galaxies has a constant slope and dispersion across 9 magnitudes (-21<M_F814W<-13). The initial data release for the HST-ACS Coma Treasury program was made available to the public in 2008 August. The images and catalogs described in this study relate to our second data release.

As part of the HST/ACS Coma Cluster Treasury Survey, we have undertaken a Keck/LRIS spectroscopic campaign to determine membership for faint dwarf galaxies. In the process, we discovered a population of Ultra Compact Dwarf galaxies (UCDs) in the core region of the Coma cluster. At the distance of Coma, UCDs are expected to have angular sizes 0.01 < R_e < 0.2 arcsec. With ACS imaging, we can resolve all but the smallest ones with careful fitting. Candidate UCDs were chosen based on magnitude, color, and degree of resolution. We spectroscopically confirm 27 objects as bona fide UCD members of the Coma cluster, a 60% success rate for objects targeted with M_R < -12. We attribute the high success rate in part to the high resolution of HST data and to an apparent large population of UCDs in Coma. We find that the UCDs tend to be strongly clustered around giant galaxies, at least in the core region of the cluster, and have a distribution and colors that are similar to globular clusters. These findings suggest that UCDs are not independent galaxies, but rather have a star cluster origin. This current study provides the dense environment datapoint necessary for understanding the UCD population.

We measure surface brightness fluctuation (SBF) magnitudes in the F814W filter and g-I colors for nine bright early-type Fornax cluster galaxies imaged with the Hubble Space Telescope Advanced Camera for Surveys (ACS). The goal is to achieve the first systematic SBF calibration for the ACS/F814W bandpass. Because of its much higher throughput, F814W is more efficient for SBF studies of distant galaxies than the ACS/F850LP bandpass that has been used to study nearby systems. Over the color range spanned by the sample galaxies, 1.06<g-I<1.32 (AB mag), the dependence of SBF magnitude mbar_I on g-I is linear to a good approximation, with slope $\sim2$. When the F850LP SBF distance measurements from the ACS Fornax Cluster Survey are used to derive absolute Mbar_I magnitudes, the dependence on g-I becomes extremely tight, with a slope of $1.8\pm0.2$ and scatter of 0.03 mag. The small observed scatter indicates both that the estimated random errors are correct, and that the intrinsic deviations from the SBF-color relation are strongly correlated between the F814W and F850LP bandpasses, as expected. The agreement with predictions from stellar population models is good, both in slope and zero point, indicating that our mean Fornax distance of 20 Mpc is accurate. The models predict curvature in the relation beyond the color limits of our sample; thus, the linear calibration should not be extrapolated naively. In the Appendices, we reconsider the Tonry ground-based and Jensen NICMOS SBF distance catalogues; we provide a correction formula to ameliorate the small apparent bias in the former, and the offset needed to make the latter consistent with other SBF studies. We also tabulate two new SBF distances to galaxies observed in the ACS Virgo Cluster Survey.

We use a highly homogeneous set of data from 132 early-type galaxies in the Virgo and Fornax clusters in order to study the properties of the globular cluster luminosity function (GCLF). The globular cluster system of each galaxy was studied using a maximum likelihood approach to model the intrinsic GCLF after accounting for contamination and completeness effects. The results presented here update our Virgo measurements and confirm our previous results showing a tight correlation between the dispersion of the GCLF and the absolute magnitude of the parent galaxy. Regarding the use of the GCLF as a standard candle, we have found that the relative distance modulus between the Virgo and Fornax clusters is systematically lower than the one derived by other distance estimators, and in particular it is 0.22mag lower than the value derived from surface brightness fluctuation measurements performed on the same data. From numerical simulations aimed at reproducing the observed dispersion of the value of the turnover magnitude in each galaxy cluster we estimate an intrinsic dispersion on this parameter of 0.21mag and 0.15mag for Virgo and Fornax respectively. All in all, our study shows that the GCLF properties vary systematically with galaxy mass showing no evidence for a dichotomy between giant and dwarf early-type galaxies. These properties may be influenced by the cluster environment as suggested by cosmological simulations.

We measure the half-light radii of globular clusters (GCs) in 43 galaxies from the ACS Fornax Cluster Survey (ACSFCS). We use these data to extend previous work in which the environmental dependencies of the half-light radii of GCs in early type galaxies in the ACS Virgo Cluster Survey (ACSVCS) were studied, and a corrected mean half-light radius (corrected for the observed environmental trends) was suggested as a reliable distance indicator. This work both increases the sample size for the study of the environmental dependencies, and adds leverage to the study of the corrected half-light radius as a possible distance indicator (since Fornax lies at a larger distance than the Virgo cluster). We study the environmental dependencies of the size of GCs using both a Principal Component Analysis as well as 2D scaling relations. We largely confirm the environmental dependencies shown in Jordan et al. (2005), but find evidence that there is a residual correlation in the mean half-light radius of GC systems with galaxy magnitude, and subtle differences in the other correlations - so there may not be a universal correction for the half-light radii of lower luminosity galaxy GC systems. The main factor determining the size of a GC in an early type galaxy is the GC color. Red GCs have <r_h> = 2.8+/-0.3 pc, while blue GCs have <r_h> = 3.4+/-0.3 pc. We show that for bright early-type galaxies (M_B < -19 mag), the uncorrected mean half-light radius of the GC system is by itself an excellent distance indicator (with error ~11%), having the potential to reach cosmologically interesting distances in the era of high angular resolution adaptive optics on large optical telescopes.

We use ACS data from the HST Treasury survey of the Coma cluster (z~0.02) to study the properties of barred galaxies in the Coma core, the densest environment in the nearby Universe. This study provides a complementary data point for studies of barred galaxies as a function of redshift and environment. From ~470 cluster members brighter than M_I = -11 mag, we select a sample of 46 disk galaxies (S0--Im) based on visual classification. The sample is dominated by S0s for which we find an optical bar fraction of 47+/-11% through ellipse fitting and visual inspection. Among the bars in the core of the Coma cluster, we do not find any very large (a_bar > 2 kpc) bars. Comparison to other studies reveals that while the optical bar fraction for S0s shows only a modest variation across low-to-intermediate density environments (field to intermediate-density clusters), it can be higher by up to a factor of ~2 in the very high-density environment of the rich Coma cluster core.

We investigate the color-magnitude relation for globular clusters (GCs) -- the so-called "blue tilt" -- detected in the ACS Fornax Cluster Survey and using the combined sample of GCs from the ACS Fornax and Virgo Cluster Surveys. We find a tilt of gamma_z=d(g-z)/dz=-0.0257 +- 0.0050 for the full GC sample of the Fornax Cluster Survey (~5800 GCs). This is slightly shallower than the value gamma_z=-0.0459 +- 0.0048 found for the Virgo Cluster Survey GC sample (~11100 GCs). The slope for the merged Fornax and Virgo datasets (~16900 GCs) is gamma_z=-0.0293 +- 0.0085, corresponding to a mass-metallicity relation of Z ~ M^0.43. We find that the blue tilt sets in at GC masses in excess of M ~ 2*10^5 M_sun. The tilt is stronger for GCs belonging to high-mass galaxies (M_* > 5 * 10^10 M_sun) than for those in low-mass galaxies (M_* < 5 * 10^10 M_sun). It is also more pronounced for GCs with smaller galactocentric distances. Our findings suggest a range of mass-metallicity relations Z_GC ~ M_GC^(0.3-0.7) which vary as a function of host galaxy mass/luminosity. We compare our observations to a recent model of star cluster self-enrichment with generally favorable results. We suggest that, within the context of this model, the proto-cluster clouds out of which the GCs formed may have had density profiles slightly steeper than isothermal and/or star formation efficiencies somewhat below 0.3. We caution, however, that the significantly different appearance of the CMDs defined by the GC systems associated with galaxies of similar mass and morphological type pose a challenge to any single mechanism that seeks to explain the blue tilt. We therefore suggest that the merger/accretion histories of individual galaxies have played a non-negligible role determining the distribution of GCs in the CMDs of individual GC systems.

Two recent empirical developments in the study of extragalactic globular cluster (GC) populations are the color-magnitude relation of the blue GCs (the "blue tilt") and the nonlinearity of the dependence of optical GC colors on metallicity. The color-magnitude relation, interpreted as a mass-metallicity relation, is thought to be a consequence of self-enrichment. Nonlinear color-metallicity relations have been shown to produce bimodal color distributions from unimodal metallicity distributions. We simulate GC populations including both a mass-metallicity scaling relation and nonlinear color-metallicity relations motivated by theory and observations. Depending on the assumed range of metallicities and the width of the GC luminosity function (GCLF), we find that the simulated populations can have bimodal color distributions with a "blue tilt" similar to observations, even though the metallicity distribution appears unimodal. The models that produce these features have the relatively high mean GC metallicities and nearly equal blue and red peaks characteristic of giant elliptical galaxies. The blue tilt is less apparent in the models with metallicities typical of dwarf ellipticals; the narrower GCLF in these galaxies has an even bigger effect in reducing the significance of their color-magnitude slopes. We critically examine the evidence for nonlinearity versus bimodal metallicities as explanations for the characteristic double-peaked color histograms of giant ellipticals and conclude that the question remains open. We discuss the prospects for further theoretical and observational progress in constraining the models presented here and for uncovering the true metallicity distributions of extragalactic globular cluster systems.

Metal-poor globular clusters (GCs) are our local link to the earliest epochs of star formation and galaxy building. Studies of extragalactic GC systems using deep, high-quality imaging have revealed a small but significant slope to the color-magnitude relation for metal-poor GCs in a number of galaxies. We present a study of the M87 GC system using deep, archival HST/ACS imaging with the F606W and F814W filters, in which we find a significant color-magnitude relation for the metal-poor GCs. The slope of this relation in the I vs. V-I color-magnitude diagram ($\gamma_I=-0.024\pm0.006$) is perfectly consistent with expectations based on previously published results using data from the ACS Virgo Cluster Survey. The relation is driven by the most luminous GCs, those with $M_I<-10$, and its significance is largest when fitting metal-poor GCs brighter than $M_I=-7.8$, a luminosity which is ~1 mag fainter than our fitted Gaussian mean for the luminosity function (LF) of blue, metal-poor GCs (~0.8 mag fainter than the mean for all GCs). These results indicate that there is a mass scale at which the correlation begins, and is consistent with a scenario where self-enrichment drives a mass-metallicity relationship. We show that previously measured half-light radii of M87 GCs from best-fit PSF-convolved King models are consistent with the more accurate measurements in this study, and we also explain how the color-magnitude relation for metal-poor GCs is real and cannot be an artifact of the photometry. We fit Gaussian and evolved Schechter functions to the luminosity distribution of GCs across all colors, as well as divided into blue and red subpopulations, finding that the blue GCs have a brighter mean luminosity and a narrower distribution than the red GCs. Finally, we present a catalog of astrometry and photometry for 2250 M87 GCs.

We present a new mass estimate of a well-studied gravitational lensing cluster, Abell 1689, from deep Chandra observations with a total exposure of 200 ks. Within r=200 h-1 kpc, the X-ray mass estimate is systematically lower than that of lensing by 30-50%. At r>200 h-1 kpc, the mass density profiles from X-ray and weak lensing methods give consistent results. The most recent weak lensing work suggest a steeper profile than what is found from the X-ray analysis, while still in agreement with the mass at large radii. Previous studies have suggested that cooler small-scale structures can bias X-ray temperature measurements or that the northern part of the cluster is disturbed. We find these scenarios unlikely to resolve the central mass discrepancy since the former requires 70-90% of the space to be occupied by these cool structures and excluding the northern substructure does not significantly affect the total mass profiles. A more plausible explanation is a projection effect. We also find that the previously reported high hard-band to broad-band temperature ratio in A1689, and many other clusters observed with Chandra, may be resulting from the instrumental absorption that decreases 10-15% of the effective area at ~1.75 keV.

We present g-z color and z-band surface brightness fluctuations (SBF) measurements for 43 early-type galaxies in the Fornax cluster imaged with HST/ACS. These are combined with measurements for Virgo cluster galaxies to derive a revised, nonlinear calibration of the z-band SBF absolute magnitude as a function of g-z. In all, we tabulate recalibrated SBF distances for 134 galaxies in Virgo, Fornax, the Virgo W' group, and NGC4697 in the Virgo Southern Extension. The calibration procedure yields a precise relative distance modulus for Fornax with respect to Virgo of 0.42+/-0.03 mag, or a distance ratio 1.214 +/- 0.017. The resulting Fornax distance modulus is 31.51 +/-0.03 +/-0.15 mag, 20.0 +/-0.3 +/-1.4 Mpc, where the second set of error bars reflects systematic uncertainty from our assumed Virgo distance of 16.5 Mpc. The rms distance scatter for the early-type Fornax galaxies is 0.49^{+0.11}_{-0.15} Mpc, consistent with its compact appearance on the sky. This translates to a depth scatter smaller than the intrinsic, or "cosmic", scatter sigma_cosmic in the SBF calibration, unlike in Virgo. As a result, we are able to place the first tight constraints on the value of sigma_cosmic, finding 0.06 +/- 0.01 mag, with a firm upper limit of sigma_cosmic<0.08 mag, at least for galaxies with g-z>1.02. We also present an alternative SBF calibration based on the `fluctuation count' parameter Nbar, a proxy for galaxy mass, which gives a consistent relative distance. Finally, we find no evidence for systematic trends of the galaxy distances with position or velocity (e.g., no current infall); the Fornax cluster appears both compact and well virialized.[abridged]

We have obtained high-resolution spectra of 23 ultra-compact dwarf galaxies (UCDs) in the Fornax cluster with -10.4>M_V>-13.5 mag (10^6<M/M_*<10^8), using FLAMES/Giraffe at the VLT. This is the largest homogeneous data set of UCD internal dynamics assembled to date. We derive dynamical M/L ratios for 15 UCDs covered by HST imaging. In the M_V-sigma plane, UCDs with M_V<-12 mag are consistent with the extrapolated Faber-Jackson relation for luminous ellipticals, while fainter UCDs are closer to the extrapolated globular cluster (GC) relation. At a given metallicity, Fornax UCDs have on average 30-40% lower M/L ratios than Virgo UCDs, suggesting possible differences in age or dark matter content between Fornax and Virgo UCDs. For our sample of Fornax UCDs we find no significant correlation between M/L ratio and mass. We combine our data with available M/L ratio measurements of compact stellar systems with 10^4<M/M_*<10^8, and normalise all M/L estimates to solar metallicity. We find that UCDs (M > 2*10^6 M_*) have M/L ratios twice as large as GCs (M < 2*10^6 M_*). We show that stellar population models tend to under-predict dynamical M/L ratios of UCDs and over-predict those of GCs. Considering the scaling relations of stellar spheroids, UCDs align well along the 'Fundamental Manifold', constituting the small-scale end of the galaxy sequence. The alignment for UCDs is especially clear for r_e >~ 7 pc, which corresponds to dynamical relaxation times that exceed a Hubble time. In contrast, GCs exhibit a broader scatter and do not appear to align along the manifold. We argue that UCDs are the smallest dynamically un-relaxed stellar systems, with M > 2*10^6 M_* and 7<r_e<100 pc. Future studies should aim at explaining the elevated M/L ratios of UCDs and the environmental dependence of their properties.

We have performed a spectroscopic study of globular clusters (GCs) in the giant elliptical NGC 5128 using the 2dF facility at the Anglo-Australian telescope. We obtained integrated optical spectra for a total of 254 GCs, 79 of which are newly confirmed on the basis of their radial velocities and spectra. In addition, we obtained an integrated spectrum of the galaxy starlight along the southern major axis. We derive an empirical metallicity distribution function (MDF) for 207 GCs (~14 of the estimated total GC system) based upon Milky Way GCs. This MDF is multimodal at high statistical significance with peaks at [Z/H]~-1.3 and -0.5. A comparison between the GC MDF and that of the stellar halo at 20 kpc (~4 Reff) reveals close coincidence at the metal-rich ends of the distributions. However, an inner 8 kpc stellar MDF shows a clear excess of metal-rich stars when compared to the GCs. We compare a higher S/N subsample (147 GCs) with two stellar population models which include non-solar abundance ratio corrections. The vast majority of our sample (~90%) appears old, with ages similar to the Milky Way GC system. There is evidence for a population of intermediate-age (~4-8 Gy) GCs (<15% of the sample) which are on average more metal-rich than the old GCs. We also identify at least one younger cluster (~1-2 Gy) in the central regions of the galaxy. Our observations are consistent with a picture where NGC 5128 has undergone at least two mergers and/or interactions involving star formation and limited GC formation since z=1, however the effect of non-canonical hot stellar populations on the integrated spectra of GCs remains an outstanding uncertainty in our GC age estimates.

The fraction of stellar mass contained in globular clusters (GCs), also measured by number as the specific frequency, is a fundamental quantity that reflects both a galaxy's early star formation and its entire merging history. We present specific frequencies, luminosities, and mass fractions for the globular cluster systems of 100 early-type galaxies in the ACS Virgo Cluster Survey, the largest homogeneous catalog of its kind. We find that 1) GC mass fractions can be high in both giants and dwarfs, but are universally low in galaxies with intermediate luminosities. 2) The behavior of specific frequency across galaxy mass is dominated by the blue GCs. 3) The GC fractions of low-mass galaxies exhibit a dependence on environment. Nearly all dwarf galaxies with high GC fractions are within 1 Mpc of the cD galaxy M87, presenting the first strong evidence that GC formation in dwarfs is biased toward dense environments. 4) GC formation in central dwarfs is biased because their stars form earliest and most intensely. Comparisons to the Millennium Simulation show that central dwarfs have older stellar populations and form more stars at higher star formation rates (SFRs) and SFR surface densities. The SFR surface density in simulated dwarfs peaks before the total SFR, naturally producing GC populations that are older and more metal-poor than the field stars. 5) Dwarfs within ~40 kpc of the giant ellipticals M87 and M49 are red and have few or no GCs, suggesting that they have been tidally stripped and have contributed their GCs to the halos of their giant neighbors. The central dwarfs with high GC mass fractions are thus likely to be the survivors most similar to the protogalaxies that assembled the rich M87 globular cluster system.(Abridged)

We present results of Hubble Space Telescope and Chandra X-ray Observatory observations of globular clusters (GCs) and low-mass X-ray binaries (LMXBs) in the central regions of Centaurus A. Out of 440 GC candidates we find that 41 host X-ray point sources that are most likely LMXBs. We fit King models to our GC candidates in order to measure their structural parameters. We find that GCs that host LMXBs are denser and more compact, and have higher encounter rates and concentrations than the GC population as a whole. We show that the higher concentrations and masses are a consequence of the dependence of LMXB incidence on central density and size plus the general trend for denser GCs to have higher masses and concentrations. We conclude that neither concentration nor mass are fundamental variables in determining the presence of LMXBs in GCs, and that the more fundamental parameters relate to central density and size.

We analyse HST surface brightness profiles for 143 early-type galaxies in the Virgo and Fornax Clusters. Sersic models provide accurate descriptions of the global profiles with a notable exception: the observed profiles deviate systematically inside a characteristic "break" radius of R_b ~ 0.02R_e where R_e is the effective radius of the galaxy. The sense of the deviation is such that bright galaxies (M_B < -20) typically show central light deficits with respect to the inward extrapolation of the Sersic model, while the great majority of low- and intermediate-luminosity galaxies (-19.5 < M_B < -15) show central light excesses; galaxies occupying a narrow range of intermediate luminosities (-20 < M_B < -19.5) are usually well fitted by Sersic models over all radii. The slopes of the central surface brightness profiles, when measured at fixed fractions of R_e, vary smoothly as a function of galaxy luminosity in a manner that depends sensitively on the choice of measurement radius. We show that a recent claim of strong bimodality in slope is likely an artifact of the galaxy selection function used in that study. To provide a more robust characterization of the inner regions of galaxies, we introduce a parameter that describes the central luminosity deficit or excess relative to the inward extrapolation of the outer Sersic model. We find that this parameter varies smoothly over the range of ~ 720 in blue luminosity spanned by the Virgo and Fornax sample galaxies, with no evidence for a dichotomy. We argue that the central light excesses (nuclei) in M_B > -19 galaxies may be the analogs of the dense central cores that are predicted by some numerical simulations to form via gas inflows. (ABRIDGED)

M85 optical transient 2006-1 (M85 OT 2006-1) is the most luminous member of the small family of V838 Mon-like objects, whose nature is still a mystery. This event took place in the Virgo cluster of galaxies and peaked at an absolute magnitude of I~-13. Here we present Hubble Space Telescope images of M85 OT 2006-1 and its environment, taken before and after the eruption, along with a spectrum of the host galaxy at the transient location. We find that the progenitor of M85 OT 2006-1 was not associated with any star forming region. The g and z-band absolute magnitudes of the progenitor were fainter than about -4 and -6 mag, respectively. Therefore, we can set a lower limit of ~50 Myr on the age of the youngest stars at the location of the progenitor that corresponds to a mass of <7 solar mass. Previously published line indices suggest that M85 has a mean stellar age of 1.6+/-0.3 Gyr. If this mean age is representative of the progenitor of M85 OT 2006-1, then we can further constrain its mass to be less than 2 solar mass. We compare the energetics and mass limit derived for the M85 OT 2006-1 progenitor with those expected from a simple model of violent stellar mergers. Combined with further modeling, these new clues may ultimately reveal the true nature of these puzzling events.

A new kinematic and dynamic study of the halo of the giant elliptical galaxy, NGC 5128, is presented. From a spectroscopically confirmed sample of 340 globular clusters and 780 planetary nebulae, the rotation amplitude, rotation axis, velocity dispersion, and the total dynamical mass are determined for the halo of NGC 5128. The globular cluster kinematics were searched for both radial dependence and metallicity dependence by subdividing the globular cluster sample into 158 metal-rich ([Fe/H] > -1.0) and 178 metal-poor ([Fe/H] < -1.0) globular clusters. Our results show the kinematics of the metal-rich and metal-poor subpopulations are quite similar. The kinematics are compared to the planetary nebula population where differences are apparent in the outer regions of the halo. The total mass of NGC 5128 is found using the Tracer Mass estimator (Evans et al. 2003), to determine the mass supported by internal random motions, and the spherical component of the Jeans equation to determine the mass supported by rotation. We find a total mass of (1.0+/-0.2) x 10^(12) Msun from the planetary nebulae data out to a projected radius of 90 kpc and (1.3+/-0.5) x 10^(12) Msun from the globular clusters out to a projected radius of 50 kpc. Lastly, we present a new and homogeneous catalog of known globular clusters in NGC 5128. This catalog combines all previous definitive cluster identifications from radial velocity studies and HST imaging studies, as well as 80 new globular clusters from a study of M.A. Beasley et al. (2007, in preparation).

The ACS Virgo Cluster Survey consists of HST ACS imaging for 100 early-type galaxies in the Virgo Cluster, observed in the F475W and F850LP filters. We derive distances for 84 of these galaxies using the method of surface brightness fluctuations (SBFs), present the SBF distance catalog, and use this database to examine the three-dimensional distribution of early-type galaxies in the Virgo Cluster. The SBF distance moduli have a mean (random) measurement error of 0.07 mag (0.5 Mpc), or roughly 3 times better than previous SBF measurements for Virgo Cluster galaxies. Five galaxies lie at a distance of ~23 Mpc and are members of the W' cloud. The remaining 79 galaxies have a narrow distribution around our adopted mean distance of 16.5+/-0.1 (random mean error) +/-1.1 Mpc (systematic). The rms distance scatter of this sample is 0.6+/-0.1 Mpc, with little dependence on morphological type or luminosity class (i.e., 0.7+/-0.1 and 0.5+/-0.1 Mpc for the giants and dwarfs, respectively). The back-to-front depth of the cluster measured from our sample of galaxies is 2.4+/-0.4 Mpc (i.e., +/-2sigma of the intrinsic distance distribution). The M87 (cluster A) and M49 (cluster B) subclusters are found to lie at distances of 16.7+/-0.2 and 16.4+/-0.2 Mpc, respectively. There may be a third subcluster associated with M86. A weak correlation between velocity and line-of-sight distance may be a faint echo of the cluster velocity distribution not having yet completely virialized. In three dimensions, Virgo's early-type galaxies appear to define a slightly triaxial distribution, with axis ratios of (1:0.7:0.5). The principal axis of the best-fit ellipsoid is inclined ~20-40 deg. from the line of sight, while the galaxies belonging to the W' cloud lie on an axis inclined by ~10-15 deg.

(Abridged) Chandra and Hubble ACS observations of 11 early-type galaxies probe the low-mass X-ray binary (LMXB) - globular cluster (GC) connection. We explore the optical properties of 270 GCs with LMXBs and 6,488 GCs without LMXBs. More massive, redder, and more compact GCs are more likely to contain LMXBs. Unlike Galactic GCs, a large number of GCs with LMXBs have half-mass relaxation times > 2.5 Gyr. We fit the dependence of the expected number of LMXBs per GC, \lambda_t, on the GC mass M, color (g-z), and half-mass radius r_{h,cor}, and find that \lambda_t \propto M^{1.24\pm0.08} 10^{0.9^{+0.2}_{-0.1} (g-z)} r_{h,cor}^{-2.2^{+0.3}_{-0.4}}. Our fit rules out that the number of LMXBs per GC is linearly proportional to GC mass and that most GCs with high X-ray luminosities contain a single LMXB. The detailed dependence of \lambda_t on GC properties appears essentially equivalent to a dependence on the encounter rate \Gamma_h and the metallicity Z, \lambda_t \propto \Gamma_h^{0.82\pm0.05} Z^{0.39\pm0.07}. Our analysis provides strong evidence that dynamical formation and metallicity play the primary roles in determining the presence of an LMXB in extragalactic GCs. The shallower than linear dependence for our sample requires an explanation by theories of dynamical binary formation. The abundance dependence is consistent with a metallicity-dependent variation in the number o neutron stars and black holes per unit mass GC, effects from irradiation induced winds or suppression of magnetic braking in metal-poor stars.

We analyze the luminosity function of the globular clusters (GCs) belonging to the early-type galaxies observed in the ACS Virgo Cluster Survey. We have obtained estimates for a Gaussian representation of the GC luminosity function (GCLF) for 89 galaxies. We have also fit the GCLFs with an "evolved Schechter function", which is meant to reflect the preferential depletion of low-mass GCs, primarily by evaporation due to two-body relaxation, from an initial Schechter mass function similar to that of young massive clusters. We find a significant trend of the GCLF dispersion with galaxy luminosity, in the sense that smaller galaxies have narrower GCLFs. We show that this narrowing of the GCLF in a Gaussian description is driven by a steepening of the GC mass function above the turnover mass, as one moves to smaller host galaxies. We argue that this behavior at the high-mass end of the GC mass function is most likely a consequence of systematic variations of the initial cluster mass function. The GCLF turnover mass M_TO is roughly constant, at ~ 2.2 x 10^5 M_sun in bright galaxies, but it decreases slightly in dwarfs with M_B >~ -18. We show that part of the variation could arise from the shorter dynamical friction timescales in smaller galaxies. We probe the variation of the GCLF to projected galactocentric radii of 20-35 kpc in the Virgo giants M49 and M87, finding that M_TO is essentially constant over these spatial scales. Our fits of evolved Schechter functions imply average dynamical mass losses (Delta) over a Hubble time that fall in the range 2 x 10^5 <~ (Delta/M_sun) < 10^6 per GC. We agree with previous suggestions that if the full GCLF is to be understood in more detail GCLF models will have to include self-consistent treatments of dynamical evolution inside time-dependent galaxy potentials. (Abridged)

The Fornax Cluster is a conspicuous cluster of galaxies in the southern hemisphere and the second largest collection of early-type galaxies within <~ 20 Mpc after the Virgo Cluster. In this paper,we present a brief introduction to the ACS Fornax Cluster Survey - a program to image, in the F475W (g_475) and F850LP (z_850) bandpasses, 43 early-type galaxies in Fornax using the Advanced Camera for Surveys (ACS) on the Hubble Space Telescope. Combined with a companion survey of Virgo, the ACS Virgo Cluster Survey, this represents the most comprehensive imaging survey to date of early-type galaxies in cluster environments in terms of depth, spatial resolution, sample size and homogeneity. We describe the selection of the program galaxies, their basic properties, and the main science objectives of the survey which include the measurement of luminosities, colors and structural parameters for globular clusters associated with these galaxies, an analysis of their isophotal properties and surface brightness profiles, and an accurate calibration of the surface brightness fluctuation distance indicator. Finally, we discuss the data reduction procedures adopted for the survey.

We present results from a study of the globular cluster luminosity function (GCLF) in a sample of 89 early-type galaxies observed as part of the ACS Virgo Cluster Survey. Using a Gaussian parametrization of the GCLF, we find a highly significant correlation between the GCLF dispersion, sigma, and the galaxy luminosity, M_B, in the sense that the GC systems in fainter galaxies have narrower luminosity functions. The GCLF dispersions in the Milky Way and M31 are fully consistent with this trend, implying that the correlation between sigma and galaxy luminosity is more fundamental than older suggestions that GCLF shape is a function of galaxy Hubble type. We show that the sigma - M_B relation results from a bonafide narrowing of the distribution of (logarithmic) cluster masses in fainter galaxies. We further show that this behavior is mirrored by a steepening of the GC mass function for relatively high masses, M >~ 3 x 10^5 M_sun, a mass regime in which the shape of the GCLF is not strongly affected by dynamical evolution over a Hubble time. We argue that this trend arises from variations in initial conditions and requires explanation by theories of cluster formation. Finally, we confirm that in bright galaxies, the GCLF "turns over" at the canonical mass scale of M_TO ~ 2 x 10^5 M_sun. However, we find that M_TO scatters to lower values (~1-2 x 10^5 M_sun) in galaxies fainter than M_B >~ -18.5, an important consideration if the GCLF is to be used as a distance indicator for dwarf ellipticals.

We examine the correlation between globular cluster (GC) color and magnitude using HST/ACS imaging for a sample of 79 early-type galaxies (-21.7<M_B<-15.2 mag) with accurate SBF distances from the ACS Virgo Cluster Survey. Using the KMM mixture modeling algorithm, we find a highly significant correlation, d(g-z)/dz = -0.037 +- 0.004, between color and magnitude for the subpopulation of blue GCs in the co-added GC color-magnitude diagram of the three brightest Virgo galaxies (M49, M87 and M60): brighter GCs are redder than their fainter counterparts. For the single GC systems of M87 and M60, we find similar correlations; M49 does not appear to show a significant trend. There is no correlation between (g-z) and M_z for GCs of the red subpopulation. The correlation d(g-z)/dg for the blue subpopulation is much weaker than d(g-z)/dz. Using Monte Carlo simulations, we attribute this to the fact that the blue subpopulation in M_g extends to higher luminosities than the red subpopulation, which biases the KMM fits. The correlation between color and M_z thus is a real effect. This conclusion is supported by biweight fits to the same color distributions. We identify two environmental dependencies of the color-magnitude relation: (1) the slope decreases in significance with decreasing galaxy luminosity; and (2) the slope is stronger for GCs at smaller galactocentric distances. We examine several mechanisms that might give rise to the observed color-magnitude relation: (1) presence of contaminators; (2) accretion of GCs from low-mass galaxies; (3) stochastic effects; (4) capture of field stars by individual GCs; and (5) GC self-enrichment. We conclude that self-enrichment and field-star capture, or a combination of these processes, offer the most promising means of explaining our observations.

Recent spectroscopic observations of planetary nebulae (PNe) in several elliptical galaxies have revealed structural and kinematical properties of the outer stellar halo regions. In order to elucidate the origin of the properties of these planetary nebula systems (PNSs), we consider the merger scenario in which an elliptical galaxy is formed by merging of spiral galaxies. Using numerical simulations, we particularly investigate radial profiles of projected PNe number densities, rotational velocities, and velocity dispersions of PNSs extending to the outer halo regions of elliptical galaxies formed from major and unequal-mass merging. We find that the radial profiles of the project number densities can be fitted to the power-law and the mean number density in the outer halos of the ellipticals can be more than an order of magnitude higher than that of the original spiral's halo. The PNSs are found to show a significant amount of rotation (V/sigma >0.5) in the outer halo regions ($R$ $>$ $5R_{\rm e}$) of the ellipticals. Two-dimensional velocity fields of PNSs are derived from the simulations and their dependences on model parameters of galaxy merging are discussed in detail. We compare the simulated kinematics of PNSs with that of the PNS observed in NGC 5128 and thereby discuss advantages and disadvantages of the merger model in explaining the observed kinematics of the PNS. We also find that the kinematics of PNSs in elliptical galaxies are quite diverse depending on the orbital configurations of galaxy merging, the mass ratio of merger progenitor spirals, and the viewing angle of the galaxies. This variation translates directly into possible biases by a factor of more than two in observational mass estimation (abridged). ~ ~

Imaging surveys with the Hubble Space Telescope (HST) have shown that roughly 50-80% of low- and intermediate-luminosity galaxies contain a compact stellar nucleus at their center, regardless of host galaxy morphological type. We combine HST imaging for early-type galaxies from the ACS Virgo Cluster Survey with ground-based long-slit spectra from KPNO to show that the masses of compact stellar nuclei in Virgo Cluster galaxies obey a tight correlation with the masses of the host galaxies. The same correlation is obeyed by the supermassive black holes (SBHs) found in predominantly massive galaxies. The compact stellar nuclei in the Local Group galaxies M33 and NGC205 are also found to fall along this same scaling relation. These results indicate that a generic by-product of galaxy formation is the creation of a Central Massive Object (CMO) -- either a SBH or a compact stellar nucleus -- that contains a mean fraction, ~ 0.2%, of the total galactic mass. In galaxies with masses greater than a few times 10^10 solar masses, SBHs appear to be the dominant mode of CMO formation.

(Abridged) The ACS Virgo Cluster Survey is an HST program to obtain high-resolution, g and z-band images for 100 early-type members of the Virgo Cluster, spanning a range of ~460 in blue luminosity. Based on this large, homogeneous dataset, we present a sharp upward revision in the frequency of nucleation in early-type galaxies brighter than M_B ~ -15 (66 < f_n < 82%), and find no evidence for nucleated dwarfs to be more concentrated to the center of Virgo than their non-nucleated counterparts. Resolved stellar nuclei are not present in galaxies brighter than M_B ~ -20.5, however, there is no clear evidence from the properties of the nuclei, or from the overall incidence of nucleation, for a change at M_B ~ -17.6, the traditional dividing point between dwarf and giant galaxies. On average, nuclei are ~3.5 mag brighter than a typical globular cluster and have a median half-light radius ~4.2 pc. Nuclear luminosities correlate with nuclear sizes and, in galaxies fainter than M_B ~ -17.6, nuclear colors. Comparing the nuclei to the "nuclear clusters" found in late-type spiral galaxies reveals a close match in terms of size, luminosity and overall frequency, pointing to a formation mechanism that is rather insensitive to the detailed properties of the host galaxy. The mean nuclear-to-galaxy luminosity ratio is indistinguishable from the mean SBH-to-bulge mass ratio, calculated in early-type galaxies with detected supermassive black holes (SBHs). We argue that compact stellar nuclei might be the low-mass counterparts of the SBHs detected in the bright galaxies, and that one should think in terms of "Central Massive Objects" -- either SBHs or compact stellar nuclei -- that accompany the formation of almost all early-type galaxies and contain a mean fraction ~0.3% of the total bulge mass.

(Abridged) We present a detailed analysis of the morphology, isophotal parameters and surface brightness profiles for 100 early-type members of the Virgo Cluster, from dwarfs (M_B = -15.1 mag) to giants (M_B = -21.8 mag). Each galaxy has been imaged in two filters, closely resembling the Sloan g and z passbands, using the Advanced Camera for Surveys on board the Hubble Space Telescope. Dust and complex morphological structures are common, with kiloparsec-scale stellar disks, bars, and nuclear stellar disks seen in 60% of galaxies with intermediate luminosity (-20 < M_B < -17), and dust seen in 42% of galaxies brighter than M_B = -18.9 mag. Dust morphologies range from faint wisps and patches on tens of parsec scales, to regular, highly organized kpc-scale dust disks, often showing evidence of recent star formation. Surface brightness profiles and isophotal parameters are derived typically within 8 kpc from the center for the brightest galaxies, and 1.5 kpc for the faintest systems, with a resolution (FWHM) of 7 pc. Based on a parametrization of the surface brightness profiles in terms of a Sersic or core-Sersic model, we find that 1) there is no evidence of a bimodal behavior of the slope, gamma, of the profile in the innermost regions; 2) although the brightest galaxies have shallow inner profiles, the shallowest profiles (lowest gamma values) are found in faint dwarf systems; 3) the widely adopted separation of early-type galaxies between "core" and "power-law" types, which had originally been prompted by the claim of a clear bimodal distribution of gamma values, is untenable; and 4) there is no evidence of a structural dichothomy between dwarf and regular ellipticals.

We use HST/ACS imaging of 100 early-type galaxies in the ACS Virgo Cluster Survey to investigate the nature of diffuse star clusters (DSCs). Compared to globular clusters (GCs), these star clusters have moderately low luminosities (M_V > -8) and a broad distribution of sizes (3 < r_h < 30 pc), but they are principally characterized by their low mean surface brightnesses which can be more than three magnitudes fainter than a typical GC (mu_g > 20 mag arcsec^-2). The median colors of diffuse star cluster systems are red, 1.1 < g-z < 1.6, which is redder than metal-rich GCs and often as red as the galaxy itself. Most DSC systems thus have mean ages older than 5 Gyr or else have super-solar metallicities. We find that 12 galaxies in our sample contain a significant excess of diffuse star cluster candidates -- nine are lenticulars (S0s), and five visibly contain dust. We also find a substantial population of DSCs in the halo of the giant elliptical M49, associated with the companion galaxy VCC 1199. Most DSC systems appear to be both aligned with the galaxy light and associated with galactic disks, but at the same time many lenticulars do not host substantial DSC populations. Diffuse star clusters in our sample share similar characteristics to those identified in other nearby lenticular, spiral, and dwarf galaxies, and we suggest that DSCs preferentially form, survive, and coevolve with galactic disks. Their properties are broadly consistent with those of merged star cluster complexes, and we note that despite being 3 - 5 magnitudes brighter than DSCs, ultra-compact dwarfs have similar surface brightnesses. The closest Galactic analogs to the DSCs are the old open clusters. We suggest that if a diffuse star cluster population did exist in the disk of the Milky Way, it would be very difficult to find. (Abridged)

We present the color distributions of globular cluster (GC) systems for 100 Virgo cluster early-type galaxies observed in the ACS Virgo Cluster Survey. The color distributions of individual GC systems are consistent with continuous trends across galaxy luminosity, color, and stellar mass. On average, almost all galaxies possess a component of metal-poor GCs, with the average fraction of metal-rich GCs ranging from 15 to 60%. The colors of both subpopulations correlate with host galaxy luminosity and color, with the red GCs having a steeper slope. To convert color to metallicity, we also introduce a preliminary (g-z)-[Fe/H] relation calibrated to Galactic, M49 and M87 GCs. This relation is nonlinear with a steeper slope for [Fe/H] < -0.8. As a result, the metallicities of the metal-poor and metal-rich GCs vary similarly with respect to galaxy luminosity and stellar mass, with relations of [Fe/H]_MP ~ L^0.16 ~ M_star^0.17 and [Fe/H]_MR ~ L^0.26 ~ M_star^0.22, respectively. Although these relations are shallower than the mass-metallicity relation predicted by wind models and observed for dwarf galaxies, they are very similar to the mass-metallicity relation for star forming galaxies in the same mass range. The offset between the two GC populations varies slowly (~ M_star^0.05) and is approximately 1 dex across three orders of magnitude in mass, suggesting a nearly universal amount of enrichment between the formation of the two populations of GCs. We also find that although the metal-rich GCs show a larger dispersion in color, it is the *metal-poor GCs* that have an equal or larger dispersion in metallicity. Like the color-magnitude relation, these relations derived from globular clusters present stringent constraints on the formation and evolution of early-type galaxies. (Abridged)

We have measured half-light radii, r_h, for globular clusters (GCs) belonging to the 100 early-type galaxies observed in the ACS Virgo Cluster Survey and the elliptical galaxy NGC 4697. An analysis of the dependencies of the measured r_h on both the properties of the GCs themselves and their host galaxies reveals that the average r_h increases with increasing galactocentric distance or, alternatively, with decreasing galaxy surface brightness. For the first time, we find that the average r_h decreases with the host galaxy color. We also show that there is no evidence for a variation of r_h with the luminosity of the GCs. Finally, we find in agreement with previous observations that the average r_h depends on the color of GCs, with red GCs being ~17% smaller than their blue counterparts. We show that this difference is probably a consequence of an intrinsic mechanism, rather than projection effects, and that it is in good agreement with the mechanism proposed in Jordan (2004). We discuss these findings in light of two simple pictures for the origin of the r_h of GCs and show that both lead to a behavior in rough agreement with the observations. After accounting for the dependencies found we show that the average GC half-light radii <r_h> can be successfully used as a standard ruler for distance estimation. We outline the methodology, and provide a calibration for its use. We find <r_h> = 2.7 +- 0.35 pc for GCs with (g-z)=1.2 mag in a galaxy with color (g-z)_{gal}=1.5 mag and at an underlying surface z-band brightness of mu_z = 21 mag arcsec^{-2}. Using this technique, we place an upper limit of 3.4 Mpc on the 1-sigma line-of-sight depth of the Virgo Cluster. Finally, we examine the form of the r_h distribution for our sample galaxies and provide an analytic expression which successfully describes this distribution. (Abridged)

As part of the Advanced Camera for Survey (ACS) Virgo Cluster Survey, we have measured Surface Brightness Fluctuations (SBF) in a sample of 100 early-type Virgo galaxies. Distances derived from these measurements are needed to explore the three-dimensional structure of the Virgo Cluster, study the intrinsic parameters of globular clusters associated with the program galaxies, and compare with the galaxy distances derived from globular cluster luminosity functions. Our SBF measurements have been performed in the F850LP bandpass of the Wide Field Channel of the ACS on the Hubble Space Telescope. These are the first measurements of this kind, and we present the first SBF calibration for this bandpass. The measured fluctuations depend on galaxy stellar population variations, which we quantify by galaxy color (g_475-z_850)_0, where g_475 and z_850 are the galaxy magnitudes, respectively, in the F475W and F850LP ACS filters. We derive the following calibration for the absolute SBF magnitude M_850: M_850 = -2.06 +/- 0.04 + (2.0 +/ 0.2) [ (g_475-z_850)_0-1.3 ] in the range 1.3 <(g_475-z_850)_0 <= 1.6, and M_850 = -2.06 +/- 0.04 + (0.9 +/- 0.2) [ (g_475-z_850)_0-1.3 ] in the range 1.0 <=(g_475-z_850)_0 <= 1.3. Physically, the two different color regimes correspond to different galaxy types: giant ellipticals and S0s at the red end, and early-type dwarfs at the blue end. For the first time in SBF studies, we are able to provide a firm empirical calibration of SBF in early-type dwarf galaxies. Our results agree with stellar population model predictions from Bruzual & Charlot (2003).

Using the HST Wide Field Camera of the Advanced Camera for Surveys (ACS), we have obtained deep (V,I) photometry of an outer-halo field in NGC 5128, to a limiting magnitude of I~29. Our photometry directly reveals the core helium-burning stellar population (the ``red clump'' or horizontal branch) in a giant E/S0 galaxy for the first time. The color-magnitude diagram displays a very wide red giant branch (RGB), an asymptotic giant branch (AGB) bump, and the red clump; no noticeable population of blue HB stars is present, confirming previous suggestions that old, very metal-poor population is not ubiquitous in the halo of this galaxy. From the upper RGB we derive the metallicity distribution, which we find to be very broad and moderately metal-rich, with average [M/H]=-0.64 and dispersion 0.49 dex. The MDF is virtually identical to that found in other halo fields observed previously with the HST, but with an enhanced metal-rich population which was partially missed in the previous surveys due to V-band incompleteness for these very red stars. Combining the metallicity sensitive colors of the RGB stars with the metallicity and age sensitive features of the AGB bump and the red clump, we infer the average age of the halo stars to be 8^{+3}_{-3.5} Gy. As part of our study, we present an empirical calibration of the ACS F606W and F814W filters to the standard V and I bands, achieved with ground-based observations of the same field made from the EMMI camera of the New Technology Telescope of the ESO La Silla Observatory.

We investigate the connection between globular clusters and ultra-compact dwarf galaxies (UCDs) by examining the properties of several compact objects associated with M87, all of which were previously classified as globular clusters. Combining imaging from the Hubble Space Telescope with ground-based Keck spectroscopy, we find two objects to have half-light radii, velocity dispersions and mass-to-light ratios that are consistent with the predictions of population synthesis models for old, metal-rich, luminous globular clusters. Three other objects are much larger, with half-light radii of approximately 20pc, and have V-band mass-to-light ratios in the range 6-9. These objects, which we consider to be UCDs, resemble the nuclei of nucleated dwarf elliptical galaxies in Virgo, having similar mass-to-light ratios, luminosities and colors. These UCDs are found to obey the extrapolated scaling relations of galaxies more closely than those of Galactic globular clusters. There appears to be a transition between the two types of stellar systems at a mass of about two million solar masses. If the UCDs are gravitationally bound, then we suggest that the presence of dark matter is the fundamental property distinguishing globular clusters from UCDs. More than half of the UCD candidates uncovered in the ACS Virgo Cluster Survey are associated with a single galaxy -- M87 -- which suggests that proximity to the Virgo center may be of critical importance for the formation of these objects. These results show that distinguishing bonafide UCDs from bright globular clusters requires a careful analysis of their detailed structural and dynamical properties, particularly their mass-to-light ratios. (ABRIDGED)

The Advanced Camera for Surveys (ACS) Virgo Cluster Survey is a large program to image 100 early-type Virgo galaxies using the F475W and F850LP bandpasses of the Wide Field Channel of the ACS instrument on the Hubble Space Telescope (HST). The scientific goals of this survey include an exploration of the three-dimensional structure of the Virgo Cluster and a critical examination of the usefulness of the globular cluster luminosity function as a distance indicator. Both of these issues require accurate distances for the full sample of 100 program galaxies. In this paper, we describe our data reduction procedures and examine the feasibility of accurate distance measurements using the method of surface brightness fluctuations (SBF) applied to the ACS Virgo Cluster Survey F850LP imaging. The ACS exhibits significant geometrical distortions due to its off-axis location in the HST focal plane; correcting for these distortions by resampling the pixel values onto an undistorted frame results in pixel correlations that depend on the nature of the interpolation kernel used for the resampling. This poses a major challenge for the SBF technique, which normally assumes a flat power spectrum for the noise. We investigate a number of different interpolation kernels and show through an analysis of simulated galaxy images having realistic noise properties that it is possible, depending on the kernel, to measure SBF distances using distortion-corrected ACS images without introducing significant additional error from the resampling. We conclude by showing examples of real image power spectra from our survey.

The study of planetary nebulae (PNe) in the nearby post-merger elliptical galaxy NGC 5128 (Cen A) now has a history of nearly twenty years. As the nearest giant elliptical, it is a prime target for extragalactic PN studies. These studies have addressed many issues including the galaxy's distance, dark matter content, halo structure, merger history, and stellar populations. We review the main PN studies that have been conducted in NGC 5128, and introduce a new study where we measure the [NII]/H-alpha ratio for 134 PNe. We find that there are no PNe in our sample that are obviously of Type I, supporting the idea that the last major star formation event in the galaxy halo occurred over 1-2 Gyr ago.

The ACS Virgo Cluster Survey is a large program to carry out multi-color imaging of 100 early-type members of the Virgo Cluster using the Advanced Camera for Surveys (ACS) on the Hubble Space Telescope. Deep F475W and F850LP images (~ SDSS g and z) are being used to study the central regions of the program galaxies, their globular cluster systems, and the three-dimensional structure of Virgo itself. In this paper, we describe in detail the data reduction procedures used for the survey, including image registration, drizzling strategies, the computation of weight images, object detection, the identification of globular cluster candidates, and the measurement of their photometric and structural parameters.

The ACIS instrument on board the Chandra X-ray Observatory has been used to carry out the first systematic study of low-mass X-ray binaries (LMXBs) in M87. We identify 174 X-ray point-sources, of which ~150 are likely LMXBs. This LMXB catalog is combined with deep F475W and F850LP images taken with ACS on HST to examine the connection between LMXBs and globular clusters in M87. Of the 1688 globular clusters in our catalog, f_X = 3.6 +- 0.5% contain a LMXB and we find that the metal-rich clusters are 3 +- 1 times more likely to harbor a LMXB than their metal-poor counterparts. In agreement with previous findings for other galaxies, we find that brighter, more metal-rich clusters are more likely to contain a LMXB. For the first time, however, we are able to demonstrate that the probability, p_X, that a given cluster will contain a LMXB depends sensitively on the dynamical properties of the host cluster. Specifically, we use the HST images to measure the half-light radius, concentration index and central density, \rho_0, for each globular, and define a parameter, \Gamma, which is related to the tidal capture and binary-neutron star exchange rate. Our preferred form for p_X is then p_X \propto \Gamma \rho_0^{-0.42\pm0.11} (Z/Z_{\odot})^{0.33\pm0.1}. We argue that if the form of p_X is determined by dynamical processes, then the observed metallicity dependence is a consequence of an increased number of neutron stars per unit mass in metal-rich globular clusters. Finally, we find no compelling evidence for a break in the luminosity distribution of resolved X-ray point sources. Instead, the LMXB luminosity function is well described by a power law with an upper cutoff at L_X ~ 10^39 erg/s. (abridged)

We utilize Sloan Digital Sky Survey imaging and spectroscopy of ~53,000 star-forming galaxies at z~0.1 to study the relation between stellar mass and gas-phase metallicity. We derive gas-phase oxygen abundances and stellar masses using new techniques which make use of the latest stellar evolutionary synthesis and photoionization models. We find a tight (+/-0.1 dex) correlation between stellar mass and metallicity spanning over 3 orders of magnitude in stellar mass and a factor of 10 in metallicity. The relation is relatively steep from 10^{8.5} - 10^{10.5} M_sun, in good accord with known trends between luminosity and metallicity, but flattens above 10^{10.5} M_sun. We use indirect estimates of the gas mass based on the H-alpha luminosity to compare our data to predictions from simple closed box chemical evolution models. We show that metal loss is strongly anti-correlated with baryonic mass, with low mass dwarf galaxies being 5 times more metal-depleted than L* galaxies at z~0.1. Evidence for metal depletion is not confined to dwarf galaxies, but is found in galaxies with masses as high as 10^{10} M_sun. We interpret this as strong evidence both of the ubiquity of galactic winds and of their effectiveness in removing metals from galaxy potential wells.

The Virgo Cluster is the dominant mass concentration in the Local Supercluster and the largest collection of elliptical and lenticular galaxies in the nearby universe. In this paper, we present an introduction to the ACS Virgo Cluster Survey: a program to image, in the F475W and F850LP bandpasses (Sloan g and z), 100 early-type galaxies in the Virgo Cluster using the Advanced Camera for Surveys on the Hubble Space Telescope. We describe the selection of the program galaxies and their ensemble properties, the choice of filters, the field placement and orientation, the limiting magnitudes of the survey, coordinated parallel observations of 100 "intergalactic" fields with WFPC2, and supporting ground-based spectroscopic observations of the program galaxies. In terms of depth, spatial resolution, sample size and homogeneity, this represents the most comprehensive imaging survey to date of early-type galaxies in a cluster environment. We briefly describe the main scientific goals of the survey which include the measurement of luminosities, metallicities, ages, and structural parameters for the many thousands of globular clusters associated with these galaxies, a high-resolution isophotal analysis of galaxies spanning a factor of ~ 450 in luminosity and sharing a common environment, the measurement of accurate distances for the full sample of galaxies using the method of surface brightness fluctuations, and a determination of the three-dimensional structure of Virgo itself.

We present the results of our analysis of the photometric data of globular clusters in the elliptical galaxy NGC 5128 (Cen A). We show that the integrated colour U-B can be an effective metallicity indicator for simple stellar populations. This is because it is sensitive to metallicity via the opacity effect but relatively insensitive to the effective main sequence turn-off temperature of the population (and thus to age) when Teff ~ 7000-12000K, that is, when the Balmer Jump is strong. This flat U-B vs Teff relation is a result of the fact that the blueing effect of the optical continuum with increasing temperature is temporarily stopped by the Balmer Jump which becomes greater with temperature in this range. In this study we use U-B and B-V as metallicity and age indicators, respectively. We first show that the use of the U-B vs B-V two-colour diagram roughly yields the metallicities and ages of the Milky Way globular clusters independently determined, and then apply the technique to the clusters in NGC 5128. There is a large range in U-B, which corresponds to [Fe/H] of -2.0 through over +0.3. The large uncertainties from the data and the models currently prevent us from pinning down their ages and metallicities. Although a constant age for all these clusters cannot be ruled out, there is a hint of the metal-rich clusters being younger. Significance of these results and caveats of the analysis are discussed.

We present the results of a photometric and spectroscopic survey of the globular cluster system of NGC 5128 (Centaurus A), a galaxy whose proximity makes it an important target for early-type galaxy studies. We imaged three fields in UBVRI that extend 50 and 30 kpc along the major and minor axes, respectively. We used both color and size information to develop efficient selection criteria for differentiating between star clusters and foreground stars. In total, we obtained new velocities for 138 globular clusters, nearly tripling the number of known clusters, and bringing the confirmed total in NGC 5128 to 215. We present a full catalog of all known GCs, with their positions, photometry, and velocities. In addition, we present catalogs of other objects observed, such as foreground stars, background galaxies, three Galactic white dwarfs, seven background QSOs, and 52 optical counterparts to known X-ray point sources. We also report an observation of the cluster G169, in which we confirm the existence of a bright emission line object. This object, however, is unlikely to be a planetary nebula, but may be a supernova remnant.

We present a study of the nearby post-merger giant elliptical galaxy, NGC 5128 (Centaurus A), in which we use the properties of its globular cluster (GC) and planetary nebula (PN) systems to constrain its evolution. Using photometric and spectroscopic data for 215 GCs presented in Paper I, we study trends in age, metallicity, and kinematics for the GC system. We confirm that the GC metallicity distribution is bimodal, and show that these two sub-populations have different properties. Using spectral line index measurements of the brightest clusters, the metal-poor GCs have old ages like the Milky Way globular clusters, while the metal-rich GCs have H-beta line-strengths that could be interpreted as a mean age of ~5 (+3/-2) Gyr. Both populations appear to have [Mg/Fe] ratios consistent with that of the Galactic GC system, although this quantity is not very well-constrained. The kinematics of the metal-rich GCs are similar to those of the planetary nebulae, exhibiting significant rotation about a misaligned axis, while the metal-poor GCs have a higher velocity dispersion and show a weaker kinematic correlation with the field stars. The total gravitating mass of NGC 5128 derived from the GCs is in excellent agreement with the value derived from stellar (PN) kinematics. We suggest that these and other data support a picture in which the main body of NGC 5128 was formed 3-8 Gyr ago by the dissipational merger of two unequal-mass disk galaxies supplemented by the continual accretion of both gas-rich and gas-poor satellites.

The halos of elliptical galaxies are faint and difficult to explore, but they contain vital clues to both structure and formation. We present the results of an imaging and spectroscopic survey for planetary nebulae (PNe) in the nearby elliptical NGC 5128. We extend the work of Hui et al.(1995) well into the halo of the galaxy--out to distances of 100 and 50 kpc along the major and minor axes. We now know of 1141 PNe in NGC 5128, 780 of which are confirmed. Of these 780 PNe, 349 are new from this survey, and 148 are at radii beyond 20 kpc. PNe exist at distances up to 80 kpc (~15 r_e), showing that the stellar halo extends to the limit of our data. This study represents by far the largest kinematic study of an elliptical galaxy to date, both in the number of velocity tracers and in radial extent. We confirm the large rotation of the PNe along the major axis, and show that it extends in a disk-like feature into the halo. The rotation curve of the stars flattens at ~100 km/s with V/sigma between 1 and 1.5, and with the velocity dispersion of the PNe falling gradually at larger radii. The two-dimensional velocity field exhibits a zero-velocity contour with a pronounced twist, showing that the galaxy potential is likely triaxial in shape, tending toward prolate. The total dynamical mass of the galaxy within 80 kpc is ~5 x 10^{11} M_sun, with M/L_B ~ 13. This mass-to-light ratio is much lower than what is typically expected for elliptical galaxies.

We present a determination of the `Cosmic Optical Spectrum' of the Universe, i.e. the ensemble emission from galaxies, as determined from the red-selected Sloan Digital Sky Survey main galaxy sample and compare with previous results of the blue-selected 2dF Galaxy Redshift Survey. Broadly we find good agreement in both the spectrum and the derived star-formation histories. If we use a power-law star-formation history model where star-formation rate $\propto (1+z)^\beta$ out to z=1, then we find that $\beta$ of 2 to 3 is still the most likely model and there is no evidence for current surveys missing large amounts of star formation at high redshift. In particular `Fossil Cosmology' of the local universe gives measures of star-formation history which are consistent with direct observations at high redshift. Using the photometry of SDSS we are able to derive the cosmic spectrum in absolute units (i.e.$ W \AA$^{-1}$ Mpc$^{-3}$) at 2--5\AA resolution and find good agreement with published broad-band luminosity densities. For a Salpeter IMF the best fit stellar mass/light ratio is 3.7--7.5 $\Msun/\Lsun$ in the r-band (corresponding to $\omstars h = 0.0025$--0.0055) and from both the stellar emission history and the H$\alpha$ luminosity density independently we find a cosmological star-formation rate of 0.03--0.04 h $\Msun$ yr$^{-1}$ Mpc$^{-3}$ today.

We study the relations between stellar mass, star formation history, size and internal structure for a complete sample of 122,808 galaxies drawn from the Sloan Digital Sky Survey. We show that low-redshift galaxies divide into two distinct families at a stellar mass of 3 \times 10^10 M_sol. Lower mass galaxies have young stellar populations, low surface mass densities and the low concentrations typical of disks. A significant fraction of the lowest mass galaxies in our sample have experienced recent starbursts. At given stellar mass, the sizes of low mass galaxies are log- normally distributed with dispersion sigma(ln R_50) \sim 0.5, in excellent agreement with the idea that they form with little angular momentum loss through cooling and condensation in a gravitationally dominant dark matter halo. Their median stellar surface mass density scales with stellar mass as mu* propto M_*^0.54, suggesting that the stellar mass of a disk galaxy is proprtional to the three halves power of its halo mass. This suggests that the efficiency of the conversion of baryons into stars in low mass galaxies increases in propor- tion to halo mass, perhaps as a result of supernova feedback processes. At stellar masses above 3 \times 10^10 M_sol, there is a rapidly increasing frac- tion of galaxies with old stellar populations, high surface mass densities and high concentrations typical of bulges. In this regime, the size distribution is log-normal, but its dispersion decreases rapidly with increasing stellar mass and the median mass surface density is approximately constant. This suggests that the star formation efficiency decreases in the highest mass halos, and that little star formation occurs in massive galaxies once they have assembled.

We develop a new method to constrain the star formation histories, dust attenuation and stellar masses of galaxies. It is based on two stellar absorption line indices, the 4000 Angstrom break strength and the Balmer absorption line index Hdelta_A. Together, these indices allow us to constrain the mean stellar ages of galaxies and the fractional stellar mass formed in bursts over the past few Gyr. A comparison with broad band photometry then yields estimates of dust attenuation and of stellar mass. We generate a large library of Monte Carlo realizations of different star formation histories, including starbursts of varying strength and a range of metallicities. We use this library to generate median likelihood estimates of burst mass fractions, dust attenuation strengths, stellar masses and stellar mass-to-light ratios for a sample of 122,208 galaxies drawn from the Sloan Digital Sky Survey. The typical 95% confidence range in our estimated stellar masses is +-40%. We study how the stellar mass-to-light ratios of galaxies vary as a function of absolute magnitude, concentration index and photometric pass-band and how dust attenuation varies as a function of absolute magnitude and 4000 Angstrom break strength. We also calculate how the total stellar mass of the present Universe is distributed over galaxies as a function of their mass, size, concentration, colour, burst mass fraction and surface mass density. We find that most of the stellar mass in the local Universe resides in galaxies that have stellar masses \~5\times 10^10 M_sol, half light radii ~3 kpc, and half- light surface mass densities ~10^9 M_sol/kpc^2. The distribution of D(4000) is strongly bimodal, showing a clear division between galaxies dominated by old stellar populations and galaxies with more recent star formation.

Stellar streams in galaxy halos are the natural consequence of a history of merging and accretion. We present evidence for a blue tidal stream of young stars in the nearest giant elliptical galaxy, NGC 5128 (Centaurus A). Using optical UBVR color maps, unsharp masking, and adaptive histogram equalization, we detect a blue arc in the northwest portion of the galaxy that traces a partial ellipse with an apocenter of 8 kpc. We also report the discovery of numerous young star clusters that are associated with the arc. The brightest of these clusters is spectroscopically confirmed, has an age of 350 Myr, and may be a proto-globular cluster. It is likely that this arc, which is distinct from the surrounding shell system and the young jet-related stars in the northeast, is a tidally disrupted stellar stream orbiting the galaxy. Both the age derived from the integrated optical colors of the stream and its dynamical disruption timescale have values of 200-400 Myr. We propose that this stream of young stars was formed when a dwarf irregular galaxy, or similar sized gas fragment, underwent a tidally triggered burst of star formation as it fell into NGC 5128 and was disrupted 300 Myr ago. The stars and star clusters in this stream will eventually disperse and become part of the main body of NGC 5128, suggesting that the infall of gas-rich dwarfs plays a role in the building of stellar halos and globular cluster systems.

Planetary nebulae (PNe) can be used to trace intermediate and old stellar populations in galaxies and the intracluster medium out to approximately 20 Mpc. PNe can be easily identified with narrow-band surveys, and their chemical abundances and radial velocities can be measured with multi-slit or multi-fiber spectrographs. We briefly review what has been accomplished to date, and include tables (with references) that summarize the observations of PNe in 42 galaxies and in the intracluster media of the Virgo and Fornax clusters. The surveys for intracluster PNe suggest that between 20% and 50% of the clusters' total stellar mass is in the intracluster medium! We illustrate the utility of PNe for dynamical studies by presenting results from a recently completed survey of Cen A (NGC 5128) with the NOAO mosaic camera on the CTIO 4-m telescope. The PNe and globular cluster surveys extend respectively to distances of 80 Kpc and 50 Kpc. We compare the kinematics of the stars in Cen A (736 PN velocities) to the kinematics of the red and blue globular clusters (188 clusters with measured velocities). The latter include 125 newly confirmed globular clusters from our survey.

We report the discovery of five quasars with redshifts of 4.67 - 5.27 and z'-band magnitudes of 19.5-20.7 M_B ~ -27. All were originally selected as distant quasar candidates in optical/near-infrared photometry from the Sloan Digital Sky Survey (SDSS), and most were confirmed as probable high-redshift quasars by supplementing the SDSS data with J and K measurements. The quasars possess strong, broad Lyman-alpha emission lines, with the characteristic sharp cutoff on the blue side produced by Lyman-alpha forest absorption. Three quasars contain strong, broad absorption features, and one of them exhibits very strong N V emission. The amount of absorption produced by the Lyman-alpha forest increases toward higher redshift, and that in the z=5.27 object (D_A ~ 0.7) is consistent with a smooth extrapolation of the absorption seen in lower redshift quasars. The high luminosity of these objects relative to most other known objects at z >~ 5 makes them potentially valuable as probes of early quasar properties and of the intervening intergalactic medium.

Recent discussion of the effects of finite source size on high magnification microlensing events due to MACHOs motivates a study into the feasibility of observing such effects and extracting the source radius. Simulated observations are generated by adding Gaussian error to points sampled on theoretical microlensing light curves for a limb darkened, extended source. These simulated data sets are fitted in an attempt to see how well the fits extract the radius of the source. The source size can be fitted with reasonable accuracy only if the impact parameter of the event, $p$, is less than the stellar radius, $R_{\star}$. It is possible to distinguish ``crossing'' events, ones where $p<R_{\star}$, from ``non--crossing'' events if the light curve is well sampled around the peak and photometric error is small --- i.e.~$\geq 3$ observations while the lens transits the disk of the source, and $\sigma_{phot}<0.08\;{\rm mag}$. These requirements are just within the reach of current observational programs; the use of an early-warning system and multiple observing sites should increase the likelihood that $R_{\star}$ can be fitted. The programs used to simulate and fit data can be obtained via anonymous {\tt ftp}.