Source author record

Chengze Liu

Chengze Liu appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

34works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

34 published item(s)

preprint2022arXiv

An Extended Halo-based Group/Cluster finder: application to the DESI legacy imaging surveys DR8

We extend the halo-based group finder developed by \citet[][]{Yang2005a} to use data {\it simultaneously} with either photometric or spectroscopic redshifts. A mock galaxy redshift survey constructed from a high-resolution N-body simulation is used to evaluate the performance of this extended group finder. For galaxies with magnitude ${\rm z\le 21}$ and redshift $0<z\le 1.0$ in the DESI legacy imaging surveys (the Legacy Surveys), our group finder successfully identifies more than 60\% of the members in about $90\%$ of halos with mass $\ga 10^{12.5}\msunh$. Detected groups with mass $\ga 10^{12.0}\msunh$ have a purity (the fraction of true groups) greater than 90\%. The halo mass assigned to each group has an uncertainty of about 0.2 dex at the high mass end $\ga 10^{13.5}\msunh$ and 0.40 dex at the low mass end. Groups with more than 10 members have a redshift accuracy of $\sim 0.008$. We apply this group finder to the Legacy Surveys DR8 and find 5.2 Million groups with at least 3 members. About 387,000 of these groups have at least 10 members. The resulting catalog containing 3D coordinates, richness, halo masses, and total group luminosities, is made publicly available.

preprint2022arXiv

Groups and protocluster candidates in the CLAUDS and HSC-SSP joint deep surveys

Using the extended halo-based group finder developed by Yang et al. (2021), which is able to deal with galaxies via spectroscopic and photometric redshifts simultaneously, we construct galaxy group and candidate protocluster catalogs in a wide redshift range ($0 < z < 6$) from the joint CFHT Large Area $U$-band Deep Survey (CLAUDS) and Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) deep data set. Based on a selection of 5,607,052 galaxies with $i$-band magnitude $m_{i} < 26$ and a sky coverage of $34.41\ {\rm deg}^2$, we identify a total of 2,232,134 groups, within which 402,947 groups have at least three member galaxies. We have visually checked and discussed the general properties of those richest groups at redshift $z>2.0$. By checking the galaxy number distributions within a $5-7\ h^{-1}\mathrm{Mpc}$ projected separation and a redshift difference $Δz \le 0.1$ around those richest groups at redshift $z>2$, we identified a list of 761, 343 and 43 protocluster candidates in the redshift bins $2\leq z<3$, $3\leq z<4$ and $z \geq 4$, respectively. In general, these catalogs of galaxy groups and protocluster candidates will provide useful environmental information in probing galaxy evolution along the cosmic time.

preprint2022arXiv

The Next Generation Virgo Cluster Survey. XXXIII. Stellar Population Gradients in the Virgo Cluster Core Globular Cluster System

We present a study of the stellar populations of globular clusters (GCs) in the Virgo Cluster core with a homogeneous spectroscopic catalog of 692 GCs within a major axis distance $R_{\rm maj} = $ 840 kpc from M87. We investigate radial and azimuthal variations in the mean age, total metallicity, [Fe/H], and $α$-element abundance, of blue (metal-poor) and red (metal-rich) GCs using their co-added spectra. We find that the blue GCs have a steep radial gradient in [Z/H] within $R_{\rm maj} =$ 165 kpc, with roughly equal contributions from [Fe/H] and [$α$/Fe], and flat gradients beyond. By contrast, the red GCs show a much shallower gradient in [Z/H], which is entirely driven by [Fe/H]. We use GC-tagged Illustris simulations to demonstrate an accretion scenario where more massive satellites (with more metal- and $α$-rich GCs) sink further into the central galaxy than less massive ones, and where the gradient flattening occurs because of the low GC occupation fraction of low-mass dwarfs disrupted at larger distances. The dense environment around M87 may also cause the steep [$α$/Fe] gradient of the blue GCs, mirroring what is seen in the dwarf galaxy population. The progenitors of red GCs have a narrower mass range than those of blue GCs, which makes their gradients shallower. We also explore spatial inhomogeneity in GC abundances, finding that the red GCs to the northwest of M87 are slightly more metal-rich. Future observations of GC stellar population gradients will be useful diagnostics of halo merger histories.

preprint2021arXiv

A giant central red disk galaxy at redshift $z=0.76$: challenge to theories of galaxy formation

We report a giant red central disk galaxy in the XMM-LSS north region. The region is covered with a rich variety of multiband photometric and spectroscopic observations. Using the photometric data of the Canada-France-Hawaii Telescope Legacy Survey (CFHTLS) and spectroscopic observation of the Baryon Oscillation Spectroscopic Survey (BOSS), we find that the galaxy has a stellar mass of $\sim 10^{11.6}$ times of the solar mass $M_\odot$. The galaxy has a red color and has an old stellar population, and thus its star formation has stopped. With the photometric image data of Hyper Suprime-Cam (HSC) Subaru Strategic Program, we demonstrate that its luminosity profile is perfectly described by a Sérsic form with $n=1.22$ indicating disk morphology. We also analyze its environment based on the VIMOS Public Extragalactic Redshift Survey (VIPERS) photometric catalog, and find that its close neighbors are all less massive, indicating that our observed galaxy is sitting at the center of its host halo. Existence of the giant red central disk galaxy seriously challenges the current standard paradigm of galaxy formation, as there is no known physical mechanism to explain the quenching of its star formation. This conclusion is supported by state-of-the-art hydrodynamical simulations of galaxy formation.

preprint2021arXiv

GOLDRUSH. IV. Luminosity Functions and Clustering Revealed with ~4,000,000 Galaxies at z~2-7: Galaxy-AGN Transition, Star Formation Efficiency, and Implication for Evolution at z>10

We present new measurements of rest-UV luminosity functions and angular correlation functions from 4,100,221 galaxies at z~2-7 identified in the Subaru/Hyper Suprime-Cam survey and CFHT Large-Area U-band Survey. The obtained luminosity functions at z~4-7 cover a very wide UV luminosity range of ~0.002-2000L*uv combined with previous studies, revealing that the dropout luminosity function is a superposition of the AGN luminosity function dominant at Muv<-24 mag and the galaxy luminosity function dominant at Muv>-22 mag, consistent with galaxy fractions based on 1037 spectroscopically-identified sources. Galaxy luminosity functions estimated from the spectroscopic galaxy fractions show the bright end excess beyond the Schechter function at >2sigma levels, which is possibly made by inefficient mass quenching, low dust obscuration, and/or hidden AGN activity. By analyzing the correlation functions at z~2-6 with halo occupation distribution models, we find a weak redshift evolution (within 0.3 dex) of the ratio of the star formation rate (SFR) to the dark matter accretion rate, SFR/(dMh/dt), indicating the almost constant star formation efficiency at z~2-6, as suggested by our earlier work at z~4-7. Meanwhile, the ratio gradually increases with decreasing redshift at z<5 within 0.3 dex, which quantitatively reproduces the redshift evolution of the cosmic SFR density, suggesting that the evolution is primarily driven by the increase of the halo number density due to the structure formation, and the decrease of the accretion rate due to the cosmic expansion. Extrapolating this calculation to higher redshifts assuming the constant efficiency suggests a rapid decrease of the SFR density at z>10 with $\propto10^{-0.5(1+z)}$, which will be directly tested with JWST.

preprint2021arXiv

The Color Gradients of the Globular Cluster Systems in M87 and M49

Combining data from the ACS Virgo Cluster Survey (ACSVCS) and the Next Generation Virgo cluster Survey (NGVS), we extend previous studies of color gradients of the globular cluster (GC) systems of the two most massive galaxies in the Virgo cluster, M87 and M49, to radii of $\sim 15~R_e$ ($\sim 200$ kpc for M87 and $\sim 250$ kpc for M49). We find significant negative color gradients, i.e., becoming bluer with increasing distance, out to these large radii. The gradients are driven mainly by the outwards decrease of the ratio of red to blue GC numbers. The color gradients are also detected out to $\sim 15~R_e$ in the red and blue sub-populations of GCs taken separately. In addition, we find a negative color gradient when we consider the satellite low-mass elliptical galaxies as a system, i.e., the satellite galaxies closer to the center of the host galaxy usually have redder color indices, both for their stars and GCs. According to the "two phase" formation scenario of massive early-type galaxies, the host galaxy accretes stars and GCs from low-mass satellite galaxies in the second phase. So the accreted GC system naturally inherits the negative color gradient present in the satellite population. This can explain why the color gradient of the GC system can still be observed at large radii after multiple minor mergers.

preprint2020arXiv

A discrete chemo-dynamical model of M87's globular clusters: Kinematics extending to ~ 400 kpc

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}/σ\sim$ 0.4, and blue GCs have weak rotation with $V_{\rm max}/σ\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.

preprint2020arXiv

Populating HI gas in dark matter halos: I. method

We combine data from the Sloan Digital Sky Survey (SDSS) and the Arecibo Legacy Fast ALFA Survey (ALFALFA) to establish an empirical model for the HI gas content within dark matter halos. A cross-match between our SDSS DR7 galaxy group sample and the ALFALFA HI sources provides a catalog of 16,520 HI-galaxy pairs within 14,270 galaxy groups (halos). Using these matched pairs, we model the HI gas mass distributions within halos using two components: 1) {\it in situ} galaxy relations that involve the HI masses, colors $({\rm g-r})$ and stellar masses 2) an {\it ex situ} dependence of the HI mass on the halo mass/environment. We find that if we solely use galaxy associated scaling relations to predict the HI gas distribution (solely component 1), the number of HI detections is significantly over-predicted with respect the ALFALFA observations. We introduce a concept for the survival of the HI masses/members within halos of different masses labelled as the `efficiency' factor, in order to describe the probability that a halo has in retaining its HI detections. Taking the above consideration into account we construct a `halo based HI mass model' which does not only predict the HI masses of galaxies, but also yields similar number, stellar, halo mass and satellite fraction distributions to the HI detections retrieved from observational data.

preprint2020arXiv

Revealing the formation mechanism of the shell galaxy NGC474 with MUSE

Stellar shells around galaxies could provide precious insights into their assembly history. However, their formation mechanism remains poorly empirically constrained, in particular the type of galaxy collisions at their origin. We present MUSE@VLT data of the most prominent outer shell of NGC 474, to constrain its formation history. The stellar shell spectrum is clearly detected, with a signal-to-noise ratio of around 65 pix$^{-1}$. We use a full spectral fitting method to determine the line-of-sight velocity and the age and metallicity of the shell and associated point-like sources within the MUSE field of view. We detect six GC candidates and eight PN candidates which are all kinematically associated to the stellar shell. We show that the shell has an intermediate metallicity, [M/H] = $-0.83^{+0.12}_{-0.12}$ and a possible $α$-enrichment, [$α$/Fe] ~ 0.3. Assuming the material of the shell comes from a lower mass companion, and that the latter had no initial metallicity gradient, such a stellar metallicity would constrain the mass of the progenitor to be around 7.4 x 10^8 M$_\odot$, implying a merger mass ratio of about 1:100. However our census of PNs and earlier photometry of the shell would suggest a much higher ratio, around 1:20. Given the uncertainties, this difference is significant only at the ~1sigma level. We discuss the characteristics of the progenitor, in particular whether the progenitor could also be composed of stars from the low metallicity outskirts from a more massive galaxy. Ultimately, the presented data does not allow us to put a firm constraint on the progenitor mass. We show that at least two globular cluster candidates possibly associated with the shell are quite young, with ages below 1.5~Gyr. We also note the presence of a young (~1Gyr) stellar population in the center of NGC 474. The two may have resulted from the same event.

preprint2020arXiv

Star Formation in Massive Galaxies at Redshift $z \sim 0.5$

It is believed that massive galaxies have quenched their star formation because of active galactic nucleus feedback. However, recent studies have shown that some massive galaxies are still forming stars. We analyze the morphology of star formation regions for galaxies of stellar mass larger than 10$^{11.3}$ M$_{\odot}$ at around redshift $z_r=0.5$ using $u-z$ color images. We find that about $20\%$ of the massive galaxies are star-forming (SF) galaxies, and most of them ($\sim 85\%$) have asymmetric structures induced by recent mergers. Moreover, for these asymmetric galaxies, we find that the asymmetry of the SF regions becomes larger for bluer galaxies. Using the Illustris simulation, we can qualitatively reproduce the observed relation between asymmetry parameter and color. Furthermore, using the merger trees in the simulation, we find a correlation between the color of the main branch galaxies at $z_r=0.5$ and the sum of the Star Formation Rates (SFRs) of the recently accreted galaxies, which implies that star formation of the accreted galaxies has contributed to the observed star formation of the massive (host) galaxies (ex situ star formation). Furthermore, we find two blue and symmetric galaxies, candidates for massive blue disks, in our observed sample, which indicates that about $\sim 10\%$ of massive SF galaxies are forming stars in the normal mode of disk star formation (in situ star formation). With the simulation, we find that the disk galaxies at $z_r \approx 0.5$ should have experienced few major mergers during the last 4.3 Gyrs.

preprint2020arXiv

The Next Generation Virgo Cluster Survey (NGVS). XXX. Ultra-Diffuse Galaxies and their Globular Cluster Systems

We present a study of ultra-diffuse galaxies (UDGs) in the Virgo Cluster based on deep imaging from the Next Generation Virgo Cluster Survey (NGVS). Applying a new definition for the UDG class based on galaxy scaling relations, we define samples of 44 and 26 UDGs using expansive and restrictive selection criteria, respectively. Our UDG sample includes objects that are significantly fainter than previously known UDGs: i.e., more than half are fainter than $\langleμ\rangle_e \sim27.5$ mag arcsec$^{-2}$. The UDGs in Virgo's core region show some evidence for being structurally distinct from "normal" dwarf galaxies, but this separation disappears when considering the full sample of galaxies throughout the cluster. UDGs are more centrally concentrated in their spatial distribution than other Virgo galaxies of similar luminosity, while their morphologies demonstrate that at least some UDGs owe their diffuse nature to physical processes---such as tidal interactions or low-mass mergers---that are at play within the cluster environment. The globular cluster (GC) systems of Virgo UDGs have a wide range in specific frequency ($S_N$), with a higher mean $S_N$ than "normal" Virgo dwarfs, but a lower mean $S_N$ than Coma UDGs at fixed luminosity. Their GCs are predominantly blue, with a small contribution from red clusters in the more massive UDGs. The combined GC luminosity function is consistent with those observed in dwarf galaxies, showing no evidence of being anomalously luminous. The diversity in their morphologies and their GC properties suggests no single process has given rise to all objects within the UDG class. Based on the available evidence, we conclude that UDGs are simply those systems that occupy the extended tails of the galaxy size and surface brightness distributions.

preprint2019arXiv

The Next Generation Virgo Cluster Survey. XVII. A Search for Planetary Nebulae in Virgo Cluster Globular Clusters

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, $α$, and measured $α$ in the Virgo cluster GCs to be $α\sim 3.9_{-0.7}^{+5.2}\times 10^{-8}\mathrm{PN}/L_\odot$. $α$ 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 $α$ 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.

preprint2016arXiv

A Chandra Study of the Image Power Spectra of 41 Cool Core and Non-Cool Core Galaxy Clusters

In this work we propose a new diagnostic to segregate cool core (CC) clusters from non-cool core (NCC) clusters by studying the two-dimensional power spectra of the X-ray images observed with the Chandra X-ray observatory. Our sample contains 41 members ($z=0.01\sim 0.54$), which are selected from the Chandra archive when a high photon count, an adequate angular resolution, a relatively complete detector coverage, and coincident CC-NCC classifications derived with three traditional diagnostics are simultaneously guaranteed. We find that in the log-log space the derived image power spectra can be well represented by a constant model component at large wavenumbers, while at small wavenumbers a power excess beyond the constant component appears in all clusters, with a clear tendency that the excess is stronger in CC clusters. By introducing a new CC diagnostic parameter, i.e., the power excess index (PEI), we classify the clusters in our sample and compare the results with those obtained with three traditional CC diagnostics. We find that the results agree with each other very well. By calculating the PEI values of the simulated clusters, we find that the new diagnostic works well at redshifts up to 0.5 for intermediately sized and massive clusters with a typical Chandra or XMM pointing observation. The new CC diagnostic has several advantages over its counterparts, e.g., it is free of the effects of the commonly seen centroid shift of the X-ray halo caused by merger event, and the corresponding calculation is straightforward, almost irrelevant to the complicated spectral analysis.

preprint2016arXiv

An empirical model to form and evolve galaxies in dark matter halos

Based on the star formation histories (SFH) of galaxies in halos of different masses, we develop an empirical model to grow galaxies in dark mattet halos. This model has very few ingredients, any of which can be associated to observational data and thus be efficiently assessed. By applying this model to a very high resolution cosmological $N$-body simulation, we predict a number of galaxy properties that are a very good match to relevant observational data. Namely, for both centrals and satellites, the galaxy stellar mass function (SMF) up to redshift $z\simeq4$ and the conditional stellar mass functions (CSMF) in the local universe are in good agreement with observations. In addition, the 2-point correlation is well predicted in the different stellar mass ranges explored by our model. Furthermore, after applying stellar population synthesis models to our stellar composition as a function of redshift, we find that the luminosity functions in $^{0.1}u$, $^{0.1}g$, $^{0.1}r$, $^{0.1}i$ and $^{0.1}z$ bands agree quite well with the SDSS observational results down to an absolute magnitude at about -17.0. The SDSS conditional luminosity functions (CLF) itself is predicted well. Finally, the cold gas is derived from the star formation rate (SFR) to predict the HI gas mass within each mock galaxy. We find a remarkably good match to observed HI-to-stellar mass ratios. These features ensure that such galaxy/gas catalogs can be used to generate reliable mock redshift surveys.

preprint2016arXiv

ELUCID - Exploring the Local Universe with reConstructed Initial Density field III: Constrained Simulation in the SDSS Volume

A method we developed recently for the reconstruction of the initial density field in the nearby Universe is applied to the Sloan Digital Sky Survey Data Release 7. A high-resolution N-body constrained simulation (CS) of the reconstructed initial condition, with $3072^3$ particles evolved in a 500 Mpc/h box, is carried out and analyzed in terms of the statistical properties of the final density field and its relation with the distribution of SDSS galaxies. We find that the statistical properties of the cosmic web and the halo populations are accurately reproduced in the CS. The galaxy density field is strongly correlated with the CS density field, with a bias that depend on both galaxy luminosity and color. Our further investigations show that the CS provides robust quantities describing the environments within which the observed galaxies and galaxy systems reside. Cosmic variance is greatly reduced in the CS so that the statistical uncertainties can be controlled effectively even for samples of small volumes.

preprint2016arXiv

Mapping the real space distributions of galaxies in SDSS DR7: I. Two Point Correlation Functions

Using a method to correct redshift space distortion (RSD) for individual galaxies, we mapped the real space distributions of galaxies in the Sloan Digital Sky Survey (SDSS) Data Release 7 (DR7). We use an ensemble of mock catalogs to demonstrate the reliability of our method. Here as the first paper in a series, we mainly focus on the two point correlation function (2PCF) of galaxies. Overall the 2PCF measured in the reconstructed real space for galaxies brighter than $^{0.1}{\rm M}_r-5\log h=-19.0$ agrees with the direct measurement to an accuracy better than the measurement error due to cosmic variance, if the reconstruction uses the correct cosmology. Applying the method to the SDSS DR7, we construct a real space version of the main galaxy catalog, which contains 396,068 galaxies in the North Galactic Cap with redshifts in the range $0.01 \leq z \leq 0.12$. The Sloan Great Wall, the largest known structure in the nearby Universe, is not as dominant an over-dense structure as appears to be in redshift space. We measure the 2PCFs in reconstructed real space for galaxies of different luminosities and colors. All of them show clear deviations from single power-law forms, and reveal clear transitions from 1-halo to 2-halo terms. A comparison with the corresponding 2PCFs in redshift space nicely demonstrates how RSDs boost the clustering power on large scales (by about $40-50\%$ at scales $\sim 10 h^{-1}{\rm {Mpc}}$) and suppress it on small scales (by about $70-80\%$ at a scale of $0.3 h^{-1}{\rm {Mpc}}$).

preprint2016arXiv

New Constraints on a complex relation between globular cluster colors and environment

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.

preprint2016arXiv

The Next Generation Virgo Cluster Survey (NGVS). XIII. The Luminosity and Mass Function of Galaxies in the Core of the Virgo Cluster and the Contribution from Disrupted Satellites

We present measurements of the galaxy luminosity and stellar mass function in a 3.71 deg$^2$ (0.3 Mpc$^2$) area in the core of the Virgo cluster, based on $ugriz$ data from the Next Generation Virgo Cluster Survey (NGVS). The galaxy sample consists of 352 objects brighter than $M_g=-9.13$ mag, the 50% completeness limit of the survey. Using a Bayesian analysis, we find a best-fit faint end slope of $α=-1.33 \pm 0.02$ for the g-band luminosity function; consistent results are found for the stellar mass function as well as the luminosity function in the other four NGVS bandpasses. We discuss the implications for the faint-end slope of adding 92 ultra compact dwarfs galaxies (UCDs) -- previously compiled by the NGVS in this region -- to the galaxy sample, assuming that UCDs are the stripped remnants of nucleated dwarf galaxies. Under this assumption, the slope of the luminosity function (down to the UCD faint magnitude limit, $M_g = -9.6$ mag) increases dramatically, up to $α= -1.60 \pm 0.06$ when correcting for the expected number of disrupted non-nucleated galaxies. We also calculate the total number of UCDs and globular clusters that may have been deposited in the core of Virgo due to the disruption of satellites, both nucleated and non-nucleated. We estimate that ~150 objects with $M_g\lesssim-9.6$ mag and that are currently classified as globular clusters, might, in fact, be the nuclei of disrupted galaxies. We further estimate that as many as 40% of the (mostly blue) globular clusters in the core of Virgo might once have belonged to such satellites; these same disrupted satellites might have contributed ~40% of the total luminosity in galaxies observed in the core region today. Finally, we use an updated Local Group galaxy catalog to provide a new measurement of the luminosity function of Local Group satellites, $α=-1.21\pm0.05$.

preprint2016arXiv

The Next Generation Virgo cluster Survey (NGVS). XXV. Fiducial panchromatic colors of Virgo core globular clusters and their comparison to model predictions

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.

preprint2016arXiv

The Next Generation Virgo Cluster Survey XVI. The Angular Momentum of Dwarf Early-Type Galaxies from Globular Cluster Satellites

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.

preprint2016arXiv

The Next Generation Virgo Cluster Survey. XIX. Tomography of Milky Way Substructures in the NGVS Footprint

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 ($α\sim 190^\circ$, $δ\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.

preprint2015arXiv

A Chandra Study of Radial Temperature Profiles of the Intra-Cluster Medium in 50 Galaxy Clusters

In order to investigate the spatial distribution of the ICM temperature in galaxy clusters in a quantitative way and probe the physics behind, we analyze the X-ray spectra of a sample of 50 galaxy clusters, which were observed with the Chandra ACIS instrument in the past 15 years, and measure the radial temperature profiles out to $0.45r_{500}$. We construct a physical model that takes into account the effects of gravitational heating, thermal history (such as radiative cooling, AGN feedback, and thermal conduction) and work done via gas compression, and use it to fit the observed temperature profiles by running Bayesian regressions. The results show that in all cases our model provides an acceptable fit at the 68% confidence level. To further validate this model we select nine clusters that have been observed with both Chandra (out to $\gtrsim 0.3r_{500}$) and Suzaku (out to $\gtrsim 1.5r_{500}$), fit their Chandra spectra with our model, and compare the extrapolation of the best-fits with the Suzaku measurements. We find that the model profiles agree with the Suzaku results very well in seven clusters. In the rest two clusters the difference between the model and observation is possibly caused by local thermal substructures. Our study also implies that for most of the clusters the assumption of hydrostatic equilibrium is safe out to at least $0.5r_{500}$, and the non-gravitational interactions between dark matter and its luminous counterpart is consistent with zero.

preprint2015arXiv

A Gemini/GMOS Study of Intermediate Luminosity Early-Type Virgo Cluster Galaxies. I. Globular Cluster and Stellar Kinematics

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.

preprint2015arXiv

Galaxies at the extremes: Ultra-diffuse galaxies in the Virgo Cluster

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.

preprint2015arXiv

The Most Massive Ultra-Compact Dwarf Galaxy in the Virgo Cluster

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 [$α$/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.

preprint2015arXiv

The Next Generation Virgo Cluster Survey. VI. The Kinematics of Ultra-compact Dwarfs and Globular Clusters in M87

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.

preprint2015arXiv

The Next Generation Virgo Cluster Survey. X. Properties of Ultra-Compact Dwarfs in the M87, M49 and M60 Regions

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.

preprint2014arXiv

The GALEX Ultraviolet Virgo Cluster Survey (GUViCS) III. The Ultraviolet Source Catalogs

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.

preprint2014arXiv

The Next Generation Virgo cluster Survey. V. modelling the dynamics of M87 with the Made-to-Measure method

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.

preprint2014arXiv

The Next Generation Virgo Cluster Survey. VIII. The Spatial Distribution of Globular Clusters in the Virgo Cluster

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)

preprint2013arXiv

The Next Generation Virgo Cluster Survey - Infrared (NGVS-IR): I. A new Near-UV/Optical/Near-IR Globular Cluster selection tool

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σ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 μ_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.

preprint2012arXiv

Exploring the Cosmic Reionization Epoch in Frequency Space: An Improved Approach to Remove the Foreground in 21 cm Tomography

Aiming to correctly restore the redshifted 21 cm signals emitted by the neutral hydrogen during the cosmic reionization processes, we re-examine the separation approaches based on the quadratic polynomial fitting technique in frequency space to investigate whether they works satisfactorily with complex foreground, by quantitatively evaluate the quality of restored 21 cm signals in terms of sample statistics. We construct the foreground model to characterize both spatial and spectral substructures of the real sky, and use it to simulate the observed radio spectra. By comparing between different separation approaches through statistical analysis of restored 21 cm spectra and corresponding power spectra, as well as their constraints on the mean halo bias $b$ and average ionization fraction $x_e$ of the reionization processes, at $z=8$ and the noise level of 60 mK we find that, although the complex foreground can be well approximated with quadratic polynomial expansion, a significant part of Mpc-scale components of the 21 cm signals (75% for $\gtrsim 6h^{-1}$ Mpc scales and 34% for $\gtrsim 1h^{-1}$ Mpc scales) is lost because it tends to be mis-identified as part of the foreground when single-narrow-segment separation approach is applied. The best restoration of the 21 cm signals and the tightest determination of $b$ and $x_e$ can be obtained with the three-narrow-segment fitting technique as proposed in this paper. Similar results can be obtained at other redshifts.

preprint2012arXiv

Stripped gas as fuel for newly formed HII regions in the encounter between VCC1249 and M49: a unified picture from NGVS and GUViCS

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.

preprint2010arXiv

The ACS Fornax Cluster Survey. X. Color Gradients of Globular Cluster Systems in Early-Type Galaxies

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.