Source author record

William E. Harris

William E. Harris 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

47works
5topics
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

47 published item(s)

preprint2022arXiv

Light from the Darkness: Detecting Ultra-Diffuse Galaxies in the Perseus Cluster through Over-densities of Globular Clusters with a Log-Gaussian Cox Process

We introduce a new method for detecting ultra-diffuse galaxies by searching for over-densities in intergalactic globular cluster populations. Our approach is based on an application of the log-Gaussian Cox process, which is a commonly used model in the spatial statistics literature but rarely used in astronomy. This method is applied to the globular cluster data obtained from the PIPER survey, a \textit{Hubble Space Telescope} imaging program targeting the Perseus cluster. We successfully detect all confirmed ultra-diffuse galaxies with known globular cluster populations in the survey. We also identify a potential galaxy that has no detected diffuse stellar content. Preliminary analysis shows that it is unlikely to be merely an accidental clump of globular clusters or other objects. If confirmed, this system would be the first of its kind. Simulations are used to assess how the physical parameters of the globular cluster systems within ultra-diffuse galaxies affect their detectability using our method. We quantify the correlation of the detection probability with the total number of globular clusters in the galaxy and the anti-correlation with increasing half-number radius of the globular cluster system. The Sérsic index of the globular cluster distribution has little impact on detectability.

preprint2021arXiv

Clearing the hurdle: The mass of globular cluster systems as a function of host galaxy mass

Current observational evidence suggests that all large galaxies contain globular clusters (GCs), while the smallest galaxies do not. Over what galaxy mass range does the transition from GCs to no GCs occur? We investigate this question using galaxies in the Local Group, nearby dwarf galaxies, and galaxies in the Virgo Cluster Survey. We consider four types of statistical models: (1) logistic regression to model the probability that a galaxy of stellar mass $M_{\star}$ has any number of GCs; (2) Poisson regression to model the number of GCs versus $M_{\star}$, (3) linear regression to model the relation between GC system mass ($\log{M_{gcs}}$) and host galaxy mass ($\log{M_{\star}}$), and (4) a Bayesian lognormal hurdle model of the GC system mass as a function of galaxy stellar mass for the entire data sample. From the logistic regression, we find that the 50% probability point for a galaxy to contain GCs is $M_{\star}=10^{6.8}M_{\odot}$. From post-fit diagnostics, we find that Poisson regression is an inappropriate description of the data. Ultimately, we find that the Bayesian lognormal hurdle model, which is able to describe how the mass of the GC system varies with $M_{\star}$ even in the presence of many galaxies with no GCs, is the most appropriate model over the range of our data. In an Appendix, we also present photometry for the little-known GC in the Local Group dwarf Ursa Major II.

preprint2020arXiv

Implementing Primordial Binaries in Simulations of Star Cluster Formation with a Hybrid MHD and Direct N-Body Method

The fraction of stars in binary systems within star clusters is important for their evolution, but what proportion of binaries form by dynamical processes after initial stellar accretion remains unknown. In previous work, we showed that dynamical interactions alone produced too few low-mass binaries compared to observations. We therefore implement an initial population of binaries in the coupled MHD and direct N-body star cluster formation code Torch. We compare simulations with, and without, initial binary populations and follow the dynamical evolution of the binary population in both sets of simulations, finding that both dynamical formation and destruction of binaries take place. Even in the first few million years of star formation, we find that an initial population of binaries is needed at all masses to reproduce observed binary fractions for binaries with mass ratios above the $q \geq 0.1$ detection limit. Our simulations also indicate that dynamical interactions in the presence of gas during cluster formation modify the initial distributions towards binaries with smaller primary masses, larger mass ratios, smaller semi-major axes and larger eccentricities. Systems formed dynamically do not have the same properties as the initial systems, and systems formed dynamically in the presence of an initial population of binaries differ from those formed in simulations with single stars only. Dynamical interactions during the earliest stages of star cluster formation are important for determining the properties of binary star systems.

preprint2020arXiv

The PIPER Survey: I. An Initial Look at the Intergalactic Globular Cluster Population in the Perseus Cluster

We describe the goals and first results of a Program for Imaging of the PERseus cluster of galaxies (PIPER). The first phase of the program builds on imaging of fields obtained with the Hubble Space Telescope (HST) ACS/WFC and WFC3/UVIS cameras. Our PIPER target fields with HST include major early-type galaxies including the active central giant NGC 1275; known Ultra-Diffuse Galaxies; and the Intracluster Medium. The resulting two-color photometry in F475W and F814W reaches deep enough to resolve and measure the globular cluster (GC) populations in the Perseus member galaxies. Here we present initial results for eight pairs of outer fields that confirm the presence of Intergalactic GCs (IGCs) in fields as distant as 740 kpc from the Perseus center (40\% of the virial radius of the cluster). Roughly 90% of these IGCs are identifiably blue (metal-poor) but there is a clear trace of a red (metal-rich) component as well, even at these very remote distances.

preprint2020arXiv

The Strikingly Metal-Rich Halo of the Sombrero Galaxy

The nature of the Sombrero galaxy (M 104 = NGC 4594) has remained elusive despite many observational studies at a variety of wavelengths. Here we present Hubble Space Telescope imaging of two fields at $\sim$16 and 33 kpc along the minor axis to examine stellar metallicity gradients in the extended spheroid. We use this imaging, extending more than 2 mag below the tip of the red giant branch (TRGB), in combination with artificial star tests to forward model observed color-magnitude diagrams (CMDs), measuring metallicity distribution functions (MDFs) at different radii along the minor axis. An important and unexpected result is that the halo of the Sombrero is strikingly metal-rich: even the outer field, located at $\sim$17 effective radii of the bulge, has a median metallicity [Z/H]$\sim$-0.15 and the fraction of stars with [Z/H]<-1.0 is negligible. This is unprecedented among massive galaxy halos studied to date, even among giant ellipticals. We find significant radial metallicity gradients, characterized by an increase in the fraction of metal-poor stars with radius and a gradient in median metallicity of $\sim$-0.01 dex/kpc. The density profile is well fit by power laws with slopes that exhibit a dependence on metallicity, with flatter slopes for more metal-poor stars. We discuss our results in the context of recent stellar MDF studies of other nearby galaxies and potential formation scenarios for the Sombrero galaxy.

preprint2019arXiv

Globular Cluster Systems and X-Ray Atmospheres in Galaxies

We compare the empirical relationships between the mass of a galaxy's globular system M_GCS, the gas mass in the hot X-ray atmosphere M_X within a fiducial radius of 5 r_e, the total gravitational mass M_grav within 5 r_e, and lastly the total halo mass M_h calibrated from weak lensing. We use a sample of 45 early-type galaxies (ETGs) for which both GCS and X-ray data are available; all the galaxies in our sample are relatively high-mass ones with M_h > 10^12 M_sun. We find that M_X ~ M_h^1.0, similar to the previously known scaling relation M_GCS ~ M_h^1.0. Both components scale much more steeply than the more well known dependence of total stellar mass M_star ~ M_h^0.3 for luminous galaxies. These results strengthen previous suggestions that feedback had little effect on formation of the globular cluster system. The current data are also used to measure the relative mass fractions of baryonic matter and dark matter (DM) within 5 r_e. We find a strikingly uniform mean of <f_DM> = 0.83 with few outliers and an rms scatter of +-0.07. This result is in good agreement with two recent suites of hydrodynamic galaxy formation models.

preprint2016arXiv

Globular Cluster Scale Sizes in Giant Galaxies: Orbital Anisotropy and Tidally Under-filling Clusters in M87, NGC 1399, and NGC 5128

We investigate the shallow increase in globular cluster half-light radii with projected galactocentric distance $R_{gc}$ observed in the giant galaxies M87, NGC 1399, and NGC 5128. To model the trend in each galaxy, we explore the effects of orbital anisotropy and tidally under-filling clusters. While a strong degeneracy exists between the two parameters, we use kinematic studies to help constrain the distance $R_β$ beyond which cluster orbits become anisotropic, as well as the distance $R_{fα}$ beyond which clusters are tidally under-filling. For M87 we find $R_β> 27$ kpc and $20 < R_{fα} < 40$ kpc and for NGC 1399 $R_β> 13$ kpc and $10 < R_{fα} < 30$ kpc. The connection of $R_{fα}$ with each galaxy's mass profile indicates the relationship between size and $R_{gc}$ may be imposed at formation, with only inner clusters being tidally affected. The best fitted models suggest the dynamical histories of brightest cluster galaxies yield similar present-day distributions of cluster properties. For NGC 5128, the central giant in a small galaxy group, we find $R_β> 5$ kpc and $R_{fα} > 30$ kpc. While we cannot rule out a dependence on $R_{gc}$, NGC 5128 is well fitted by a tidally filling cluster population with an isotropic distribution of orbits, suggesting it may have formed via an initial fast accretion phase. Perturbations from the surrounding environment may also affect a galaxy's orbital anisotropy profile, as outer clusters in M87 and NGC 1399 have primarily radial orbits while outer NGC 5128 clusters remain isotropic.

preprint2016arXiv

Where Are Most of the Globular Clusters in Today's Universe?

The total number of globular clusters (GCs) in a galaxy rises continuously with the galaxy luminosity L, while the relative number of galaxies decreases with L following the Schechter function. The product of these two very nonlinear functions gives the relative number of GCs contained by all galaxies at a given L. It is shown that GCs, in this universal sense, are most commonly found in galaxies within a narrow range around $L_{\star}$. In addition, blue (metal-poor) GCs outnumber the red (metal-richer) ones globally by 4 to 1 when all galaxies are added, pointing to the conclusion that the earliest stages of galaxy formation were especially favorable to forming massive, dense star clusters.

preprint2015arXiv

Estimating the Galactic Mass Profile in the Presence of Incomplete Data

A powerful method to measure the mass profile of a galaxy is through the velocities of tracer particles distributed through its halo. Transforming this kind of data accurately to a mass profile M(r), however, is not a trivial problem. In particular, limited or incomplete data may substantially affect the analysis. In this paper we develop a Bayesian method to deal with incomplete data effectively; we have a hybrid-Gibbs sampler that treats the unknown velocity components of tracers as parameters in the model. We explore the effectiveness of our model using simulated data, and then apply our method to the Milky Way using velocity and position data from globular clusters and dwarf galaxies. We find that in general, missing velocity components have little effect on the total mass estimate. However, the results are quite sensitive to the outer globular cluster Pal 3. Using a basic Hernquist model with an isotropic velocity dispersion, we obtain credible regions for the cumulative mass profile M(r) of the Milky Way, and provide estimates for the model parameters with 95 percent Bayesian credible intervals. The mass contained within 260 kpc is 1.37x10^12 solar masses, with a 95 percent credible interval of (1.27,1.51)x10^12 solar masses. The Hernquist parameters for the total mass and scale radius are 1.55 (+0.18/-0.13)x10^12 solar masses and 16.9 (+4.8/-4.1) kpc, where the uncertainties span the 95 percent credible intervals. The code we developed for this work, Galactic Mass Estimator (GME), will be available as an open source package in the R Project for Statistical Computing.

preprint2015arXiv

Globular Cluster Systems in Brightest Cluster Galaxies. II: NGC 6166

We present new deep photometry of the globular cluster system (GCS) around NGC 6166, the central supergiant galaxy in Abell 2199. HST data from the ACS and WFC3 cameras in F475W, F814W are used to determine the spatial distribution of the GCS, its metallicity distribution function (MDF), and the dependence of the MDF on galactocentric radius and on GC luminosity. The MDF is extremely broad, with the classic red and blue subpopulations heavily overlapped, but a double-Gaussian model can still formally match the MDF closely. The spatial distribution follows a Sersic-like profile detectably to a projected radius of at least $R_{gc} = 250$ kpc. To that radius, the total number of clusters in the system is N_{GC} = 39000 +- 2000, the global specific frequency is S_N = 11.2 +- 0.6, and 57\% of the total are blue, metal-poor clusters. The GCS may fade smoothly into the Intra-Cluster Medium of A2199; we see no clear transition from the core of the galaxy to the cD halo or the ICM. The radial distribution, projected ellipticity, and mean metallicity of the red (metal-richer) clusters match the halo light extremely well for R > 15 kpc, both of them varying as σ_{MRGC} ~ σ_{light} ~ R^-1.8. By comparison, the blue (metal-poor) GC component has a much shallower falloff σ_{MPGC} ~ R^-1.0 and a more nearly spherical distribution. This strong difference in their density distributions produces a net metallicity gradient in the GCS as a whole that is primarily generated by the population gradient. With NGC 6166 we appear to be penetrating into a regime of high enough galaxy mass and rich enough environment that the bimodal two-phase description of GC formation is no longer as clear or effective as it has been in smaller galaxies.

preprint2014arXiv

Dark Matter Halos in Galaxies and Globular Cluster Populations

We combine a new, comprehensive database for globular cluster populations in all types of galaxies with a new calibration of galaxy halo masses based entirely on weak lensing. Correlating these two sets of data, we find that the mass ratio $η\equiv M_{GCS}/M_{h}$ (total mass in globular clusters, divided by halo mass) is essentially constant at $\langle η\rangle \sim 4 \times 10^{-5}$, strongly confirming earlier suggestions in the literature. Globular clusters are the only known stellar population that formed in essentially direct proportion to host galaxy halo mass. The intrinsic scatter in $η$ appears to be at most 0.2 dex; we argue that some of this scatter is due to differing degrees of tidal stripping of the globular cluster systems between central and satellite galaxies. We suggest that this correlation can be understood if most globular clusters form at very early stages in galaxy evolution, largely avoiding the feedback processes that inhibited the bulk of field-star formation in their host galaxies. The actual mean value of $η$ also suggests that about $1/4$ of the \emph{initial} gas mass present in protogalaxies collected into GMCs large enough to form massive, dense star clusters. Finally, our calibration of $\langle η\rangle$ indicates that the halo masses of the Milky Way and M31 are $(1.2\pm0.5)\times 10^{12} M_{\odot}$ and $(3.9\pm1.8)\times 10^{12} M_{\odot}$ respectively.

preprint2014arXiv

Globular Cluster Systems in Brightest Cluster Galaxies: A Near-Universal Luminosity Function?

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.

preprint2014arXiv

Red giants in the outer halo of the elliptical galaxy NGC 5128 / Centaurus A

We used VIMOS on VLT to perform $V$ and $I$ band imaging of the outermost halo of NGC 5128 / Centaurus A ($(m-M)_0=27.91\pm0.08$), 65 kpc from the galaxy's center and along the major axis. The stellar population has been resolved to $I_0 \approx 27$ with a $50\%$ completeness limit of $I_0 = 24.7$, well below the tip of the red-giant branch (TRGB), which is seen at $I_0 \approx 23.9$. The surface density of NGC 5128 halo stars in our fields was sufficiently low that dim, unresolved background galaxies were a major contaminant in the source counts. We isolated a clean sample of red-giant-branch (RGB) stars extending to $\approx 0.8$ mag below the TRGB through conservative magnitude and color cuts, to remove the (predominantly blue) unresolved background galaxies. We derived stellar metallicities from colors of the stars via isochrones and measured the density falloff of the halo as a function of metallicity by combining our observations with HST imaging taken of NGC 5128 halo fields closer to the galaxy center. We found both metal-rich and metal-poor stellar populations and found that the falloff of the two follows the same de Vaucouleurs' law profiles from $\approx 8$ kpc out to $\approx$ 70 kpc. The metallicity distribution function (MDF) and the density falloff agree with the results of two recent studies of similar outermost halo fields in NGC 5128. We found no evidence of a "transition" in the radial profile of the halo, in which the metal-rich halo density would drop rapidly, leaving the underlying metal-poor halo to dominate by default out to greater radial extent, as has been seen in the outer halo of two other large galaxies. If NGC 5128 has such a transition, it must lie at larger galactocentric distances.

preprint2014arXiv

The Effect of Orbital Eccentricity on the Dynamical Evolution of Star Clusters

We use N-body simulations to explore the influence of orbital eccentricity on the dynamical evolution of star clusters. Specifically we compare the mass loss rate, velocity dispersion, relaxation time, and the mass function of star clusters on circular and eccentric orbits. For a given perigalactic distance, increasing orbital eccentricity slows the dynamical evolution of a cluster due to a weaker mean tidal field. However, we find that perigalactic passes and tidal heating due to an eccentric orbit can partially compensate for the decreased mean tidal field by energizing stars to higher velocities and stripping additional stars from the cluster, accelerating the relaxation process. We find that the corresponding circular orbit which best describes the evolution of a cluster on an eccentric orbit is much less than its semi-major axis or time averaged galactocentric distance. Since clusters spend the majority of their lifetimes near apogalacticon, the properties of clusters which appear very dynamically evolved for a given galactocentric distance can be explained by an eccentric orbit. Additionally we find that the evolution of the slope of the mass function within the core radius is roughly orbit-independent, so it could place additional constraints on the initial mass and initial size of globular clusters with solved orbits. We use our results to demonstrate how the orbit of Milky Way globular clusters can be constrained given standard observable parameters like galactocentric distance and the slope of the mass function. We then place constraints on the unsolved orbits of NGC 1261,NGC 6352, NGC 6496, and NGC 6304 based on their positions and mass functions.

preprint2014arXiv

The Effects of Orbital Inclination on the Scale Size and Evolution of Tidally Filling Star Clusters

We have performed N-body simulations of tidally filling star clusters with a range of orbits in a Milky Way-like potential to study the effects of orbital inclination and eccentricity on their structure and evolution. At small galactocentric distances Rgc, a non-zero inclination results in increased mass loss rates. Tidal heating and disk shocking, the latter sometimes consisting of two shocking events as the cluster moves towards and away from the disk, help remove stars from the cluster. Clusters with inclined orbits at large Rgc have decreased mass loss rates than the non-inclined case, since the strength the disk potential decreases with Rgc. Clusters with inclined and eccentric orbits experience increased tidal heating due to a constantly changing potential, weaker disk shocks since passages occur at higher Rgc, and an additional tidal shock at perigalacticon. The effects of orbital inclination decrease with orbital eccentricity, as a highly eccentric cluster spends the majority of its lifetime at a large Rgc. The limiting radii of clusters with inclined orbits are best represented by the rt of the cluster when at its maximum height above the disk, where the cluster spends the majority of its lifetime and the rate of change in rt is a minimum. Conversely, the effective radius is independent of inclination in all cases.

preprint2014arXiv

The pre-processing of subhaloes in SDSS groups and clusters

We investigate pre-processing using the observed quenched fraction of group and cluster galaxies in the Yang et al. (2007) SDSS-DR7 group catalogue in the redshift range of 0.01 < z < 0.045. We categorize group galaxies as virialized, infall or backsplash and we apply a combination of the Dressler-Shectman statistic and group member velocities to identify subhaloes. On average the fraction of galaxies that reside in subhaloes is a function of host halo mass, where more massive systems have a higher fraction of subhalo galaxies both in the overall galaxy and infall populations. Additionally, we find that between 2 < r_200 < 3 the quiescent fraction is higher in the subhalo population with respect to both the field and non-subhalo populations. At these large radii (2 < r_200 < 3), the majority of galaxies (~ 80 %) belong to the infall population and therefore, we attribute the enhanced quenching to infalling subhalo galaxies, indicating that pre-processing has occurred in the subhalo population. We conclude that pre-processing plays a significant role in the observed quiescent fraction, but only for the most massive (M_halo > 10^14.5 M_sun) systems in our sample.

preprint2013arXiv

A Catalog of Globular Cluster Systems: What Determines the Size of a Galaxy's Globular Cluster Population?

We present a catalog of 422 galaxies with published measurements of their globular cluster (GC) populations. Of these, 248 are E galaxies, 93 are S0 galaxies, and 81 are spirals or irregulars. Among various correlations of the total number of GCs with other global galaxy properties, we find that N_GC correlates well though nonlinearly with the dynamical mass of the galaxy bulge M_dyn = 4 σ_e^2 R_e /G, where σ_e is the central velocity dispersion and R_e the effective radius of the galaxy light profile. We also present updated versions of the GC specific frequency S_N and specific mass S_M versus host galaxy luminosity and baryonic mass. These graphs exhibit the previously known U-shape: highest S_N or S_M values occur for either dwarfs or supergiants, but in the midrange of galaxy size (10^9 - 10^10 L_Sun) the GC numbers fall along a well defined baseline value of S_N ~ 1 or S_M ~ 0.1, similar among all galaxy types. Along with other recent discussions, we suggest that this trend may represent the effects of feedback, which systematically inhibited early star formation at either very low or very high galaxy mass, but which had its minimum effect for intermediate masses. Our results strongly reinforce recent proposals that GC formation efficiency appears to be most nearly proportional to the galaxy halo mass M_halo. The mean "absolute" efficiency ratio for GC formation that we derive from the catalog data is M_GCS/M_halo = 6 \times 10^-5. We suggest that the galaxy-to-galaxy scatter around this mean value may arise in part because of differences in the relative timing of GC formation versus field-star formation. Finally, we find that an excellent empirical predictor of total GC population for galaxies of all luminosities is N_GC \sim (R_e σ_e)^1.3$, a result consistent with Fundamental Plane scaling relations.

preprint2013arXiv

Cluster formation in molecular clouds: I. stellar populations, star formation rates, and ionizing radiation

We present a model for the radiative output of star clusters in the process of star formation suitable for use in hydrodynamical simulations of radiative feedback. Gas in a clump, defined as a region whose density exceeds 10^4 cm^-3, is converted to stars via the random sampling of the Chabrier IMF. A star formation efficiency controls the rate of star formation. We have completed a suite of simulations which follow the evolution of accretion-fed clumps with initial masses ranging from 0 to 10^5 M_sol and accretion rates ranging from 10^-5 to 10^-1 M_sol yr^-1. The stellar content is tracked over time which allows the aggregate luminosity, ionizing photon rate, number of stars, and star formation rate (SFR) to be determined. For a fiducial clump of 10^4 M_sol, the luminosity is ~4x10^6 L_sol with a SFR of roughly 3x10^-3 M_sol yr^-1. We identify two regimes in our model. The accretion-dominated regime obtains the majority of its gas through accretion and is characterized by an increasing SFR while the reservoir-dominated regime has the majority of its mass present in the initial clump with a decreasing SFR. We show that our model can reproduce the expected number of O stars, which dominate the radiative output of the cluster. We find a nearly linear relationship between SFR and mass as seen in observations. We conclude that our model is an accurate and straightforward way to represent the output of clusters in hydrodynamical simulations with radiative feedback.

preprint2013arXiv

Do group dynamics play a role in the evolution of member galaxies?

We examine galaxy groups from the present epoch to z = 1 to explore the impact of group dynamics on galaxy evolution. We use group catalagues from the Sloan Digital Sky Survey (SDSS), the Group Environment and Evolution Collaboration (GEEC) and the high redshift GEEC2 sample to study how the observed member properties depend on galaxy stellar mass, group dynamical mass and dynamical state of the host group. We find a strong correlation between the fraction of non-star-forming (quiescent) galaxies and galaxy stellar mass, but do not detect a significant difference in the quiescent fraction with group dynamical mass, within our sample halo mass range of 10^13-10^14.5 M_sun, or with dynamical sate. However, at a redshift of approximately 0.4 we do see some evidence that the quiescent fraction in low mass galaxies (log(M_star/M_sun) < 10.5) is lower in groups with substructure. Additionally, our results show that the fraction of groups with non-Gaussian velocity distributions increases with redshift to roughly z = 0.4, while the amount of detected substructure remains constant to z = 1. Based on these results, we conclude that for massive galaxies (log(M_star/M_sun_ > 10.5), evolution is most strongly correlated to the stellar mass of a galaxy with little or no additional effect related to either the group dynamical mass or dynamical state. For low mass galaxies, we do see some evidence of a correlation between the quiescent fraction and the amount of detected substructure, highlighting the need to probe further down the stellar mass function to elucidate the role of the environment in galaxy evolution.

preprint2013arXiv

Globular Cluster Scale Sizes in Giant Galaxies: The Case of M87 and the Role of Orbital Anisotropy and Tidal Filling

We present new Hubble Space Telescope imaging of the outer regions of M87 in order to study its globular cluster (GC) population out to large galactocentric distances. We discuss particularly the relationship between GC effective radii $r_h$ and projected galactocentric distance $R_{gc}$. The observations suggest a shallow trend $r_h \propto R_{gc}^{0.14}$ out to $R_{gc} \sim 100$ kpc, in agreement with studies of other giant elliptical galaxies. To theoretically reproduce this relationship we simulate GC populations with various distributions of orbits. For an isotropic distribution of cluster orbits we find a steeper trend of $r_h \propto R_{gc}^{0.4}$. Instead we suggest that (a) if the cluster system has an orbital anisotropy profile, where orbits become preferentially radial with increasing galactocentric distance, AND (b) if clusters become more tidally under-filling with galactocentric distance, the observed relationship can be recovered. We also apply this approach to the red and blue GC populations separately and predict that red clusters are preferentially under-filling at large $R_{gc}$ and have a more isotropic distribution of orbits than blue clusters.

preprint2013arXiv

Globular clusters and supermassive black holes in galaxies: further analysis and a larger sample

We explore several correlations between various large-scale galaxy properties, particularly total globular cluster population (N_GCS), the central black hole mass (M_BH), velocity dispersion (nominally sigma_e), and bulge mass (M_dyn). Our data sample of 49 galaxies, for which both N_GC and M_BH are known, is larger than used in previous discussions of these two parameters and we employ the same sample to explore all pairs of correlations. Further, within this galaxy sample we investigate the scatter in each quantity, with emphasis on the range of published values for sigma_e and effective radius (R_e). We find that these two quantities in particular are difficult to measure consistently and caution that precise intercomparison of galaxy properties involving R_e and sigma_e is particularly difficult. Using both chi^2 and Monte Carlo Markov Chain (MCMC) fitting techniques, we show that quoted observational uncertainties for all parameters are too small to represent the true scatter in the data. We find that the correlation between M_dyn and N_GC is stronger than either the M_BH-sigma_e or M_BH-N_GC relations. We suggest that this is because both the galaxy bulge population ans N_GC were fundamentally established at an early epoch during the same series of star-forming events. By contrast, although the seed for M_BH was likely formed at a similar epoch, its growth over time is less similar from galaxy to galaxy and thus less predictable.

preprint2013arXiv

The Influence of Orbital Eccentricity on Tidal Radii of Star Clusters

We have performed N-body simulations of star clusters orbiting in a spherically symmetric smooth galactic potential. The model clusters cover a range of initial half-mass radii and orbital eccentricities in order to test the historical assumption that the tidal radius of a cluster is imposed at perigalacticon. The traditional assumption for globular clusters is that since the internal relaxation time is larger than its orbital period, the cluster is tidally stripped at perigalacticon. Instead, our simulations show that a cluster with an eccentric orbit does not need to fully relax in order to expand. After a perigalactic pass, a cluster re-captures previously unbound stars, and the tidal shock at perigalacticon has the effect of energizing inner region stars to larger orbits. Therefore, instead of the limiting radius being imposed at perigalacticon, it more nearly traces the instantaneous tidal radius of the cluster at any point in the orbit. We present a numerical correction factor to theoretical tidal radii calculated at perigalacticon which takes into consideration both the orbital eccentricity and current orbital phase of the cluster.

preprint2012arXiv

N-body models of globular clusters: metallicity, half-light radii and mass-to-light ratios

Size differences of approx. 20% between red (metal-rich) and blue (metal-poor) sub-populations of globular clusters have been observed, generating an ongoing debate as to weather these originate from projection effects or the difference in metallicity. We present direct N-body simulations of metal-rich and metal-poor stellar populations evolved to study the effects of metallicity on cluster evolution. The models start with N = 100000 stars and include primordial binaries. We also take metallicity dependent stellar evolution and an external tidal field into account. We find no significant difference for the half-mass radii of those models, indicating that the clusters are structurally similar. However, utilizing observational tools to fit half-light (or effective) radii confirms that metallicity effects related to stellar evolution combined with dynamical effects such as mass segregation produce an apparent size difference of 17% on average. The metallicity effect on the overall cluster luminosity also leads to higher mass-to-light ratios for metal-rich clusters.

preprint2012arXiv

Substructure in the Most Massive GEEC Groups: Field-like Populations in Dynamically Active Groups

The presence of substructure in galaxy groups and clusters is believed to be a sign of recent galaxy accretion and can be used not only to probe the assembly history of these structures, but also the evolution of their member galaxies. Using the Dressler-Shectman (DS) Test, we study substructure in a sample of intermediate redshift (z ~ 0.4) galaxy groups from the Group Environment and Evolution Collaboration (GEEC) group catalog. We find that 4 of the 15 rich GEEC groups, with an average velocity dispersion of ~525 km s-1, are identified as having significant substructure. The identified regions of localized substructure lie on the group outskirts and in some cases appear to be infalling. In a comparison of galaxy properties for the members of groups with and without substructure, we find that the groups with substructure have a significantly higher fraction of blue and star-forming galaxies and a parent colour distribution that resembles that of the field population rather than the overall group population. In addition, we observe correlations between the detection of substructure and other dynamical measures, such as velocity distributions and velocity dispersion profiles. Based on this analysis, we conclude that some galaxy groups contain significant substructure and that these groups have properties and galaxy populations that differ from groups with no detected substructure. These results indicate that the substructure galaxies, which lie preferentially on the group outskirts and could be infalling, do not exhibit signs of environmental effects, since little or no star-formation quenching is observed in these systems.

preprint2012arXiv

The Observational and Theoretical Tidal Radii of Globular Clusters in M87

Globular clusters have linear sizes (tidal radii) which theory tells us are determined by their masses and by the gravitational potential of their host galaxy. To explore the relationship between observed and expected radii, we utilize the globular cluster population of the Virgo giant M87. Unusually deep, high signal-to-noise images of M87 are used to measure the effective and limiting radii of approximately 2000 globular clusters. To compare with these observations, we simulate a globular cluster population that has the same characteristics as the observed M87 cluster population. Placing these simulated clusters in the well-studied tidal field of M87, the orbit of each cluster is solved and the theoretical tidal radius of each cluster is determined. We compare the predicted relationship between cluster size and projected galactocentric distance to observations. We find that for an isotropic distribution of cluster velocities, theoretical tidal radii are approximately equal to observed limiting radii for Rgc < 10 kpc. However, the isotropic simulation predicts a steep increase in cluster size at larger radii, which is not observed in large galaxies beyond the Milky Way. To minimize the discrepancy between theory and observations, we explore the effects of orbital anisotropy on cluster sizes, and suggest a possible orbital anisotropy profile for M87 which yields a better match between theory and observations. Finally, we suggest future studies which will establish a stronger link between theoretical tidal radii and observed radii.

preprint2012arXiv

The Size Difference Between Red And Blue Globular Clusters Is NOT Due To Projection Effects

Metal-rich (red) globular clusters in massive galaxies are, on average, smaller than metal-poor (blue) globular clusters. One of the possible explanations for this phenomenon is that the two populations of clusters have different spatial distributions. We test this idea by comparing clusters observed in unusually deep, high signal-to-noise images of M87 with a simulated globular cluster population in which the red and blue clusters have different spatial distributions, matching the observations. We compare the overall distribution of cluster effective radii as well as the relationship between effective radius and galactocentric distance for both the observed and simulated red and blue subpopulations. We find that the different spatial distributions does not produce a significant size difference between the red and blue subpopulations as a whole, or at a given galactocentric distance. These results suggest that the size difference between red and blue globular clusters is likely due to differences during formation or later evolution

preprint2012arXiv

Unearthing Foundations of a Cosmic Cathedral: Searching the Stars for M33's Halo

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

preprint2011arXiv

The M33 Globular Cluster System with PAndAS Data: The Last Outer Halo Cluster?

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

preprint2010arXiv

A New Catalog of Globular Clusters in the Milky Way

A new revision of the McMaster catalog of Milky Way globular clusters is available. This is the first update since 2003 and the biggest single revision since the original version of the catalog published in 1996. The list now contains a total of 157 objects classified as globular clusters. Major upgrades have been made especially to the cluster coordinates, metallicities, and structural profile parameters, and the list of parameters now also includes central velocity dispersion. NB: This paper is a stand-alone publication available only on the astro-ph archive; it will not be published separately in a journal.

preprint2010arXiv

Chemodynamics of Compact Stellar Systems in NGC 5128: How similar are Globular Clusters, Ultra-Compact Dwarfs, and Dwarf Galaxies?

Velocity dispersion measurements are presented for luminous GCs in NGC 5128 derived from high-res. UVES spectra. The measurements are made with the pPXF code that parametrically recovers line-of-sight velocity dispersions. Combining the measured velocity dispersions with surface photometry and structural parameter data from HST enables both dynamical masses and M/L ratios to be derived. The fundamental plane relations of these clusters are investigated in order to fill the apparent gap between the relations of Local Group GCs and more massive early-type galaxies. It is found that the properties of these massive stellar systems match those of nuclear clusters in dwarf elliptical galaxies and UCDs better than those of Local Group GCs, and that all objects share similarly old (>8 Gyr) ages, suggesting a possible link between the formation and evolution of dE,Ns, UCDs and massive GCs. We find a very steep correlation between dynamical (M/L) ratio and dynamical mass of the form (M/L)_dyn ~ M_dyn^(0.24+/-0.02) above M_dyn = 2x10^6 Msol. Formation scenarios are investigated with a chemical abundance analysis using absorption line strengths calibrated to the Lick/IDS index system. The results lend support to two scenarios contained within a single general formation scheme. Old, massive, super-solar [alpha/Fe] systems are formed on short (<100 Myr) timescales through the merging of single-collapse GCs which themselves are formed within single, giant molecular clouds. More intermediate- and old-aged (~3-10 Gyr), solar- to sub-solar [alpha/Fe] systems are formed on much longer (~Gyr) timescales through the stripping of dE,Ns in the 10^13-10^15 Msol potential wells of massive galaxies and galaxy clusters.

preprint2010arXiv

Consequences of cosmic microwave background-regulated star formation

It has been hypothesized that the cosmic microwave background (CMB) provides a temperature floor for collapsing protostars that can regulate the process of star formation and result in a top-heavy initial mass function at high metallicity and high redshift. We examine whether this hypothesis has any testable observational consequences. First we determine, using a set of hydrodynamic galaxy formation simulations, that the CMB temperature floor would have influenced the majority of stars formed at redshifts between z=3 and 6, and probably even to higher redshift. Five signatures of CMB-regulated star formation are: (1) a higher supernova rate than currently predicted at high redshift; (2) a systematic discrepancy between direct and indirect measurements of the high redshift star formation rate; (3) a lack of surviving globular clusters that formed at high metallicity and high redshift; (4) a more rapid rise in the metallicity of cosmic gas than is predicted by current simulations; and (5) an enhancement in the abundances of alpha elements such as O and Mg at metallicities -2 < [Fe/H] < -0.5. Observations are not presently able to either confirm or rule out the presence of these signatures. However, if correct, the top-heavy IMF of high-redshift high-metallicity globular clusters could provide an explanation for the observed bimodality of their metallicity distribution.

preprint2010arXiv

Possible Subgroups of Globular Clusters and Planetary Nebulae in NGC 5128

We use recently compiled position and velocity data for the globular cluster and planetary nebula subsystems in NGC 5128, the nearby giant elliptical, to search for evidence of past dwarf-satellite accretion events. Beyond a 10 arcmin (~11 kpc) radius in galactocentric distance, we find tentative evidence for 4 subgroups of globular clusters and 4 subgroups of planetary nebulae. These each have > 4 members within a search radius of 2 arcmin and internal velocity dispersion of < 40 km/s, typical parameters for a dwarf galaxy. In addition, 2 of the globular cluster groupings overlap with 2 of the planetary nebulae groupings, and 2 subgroupings also appear to overlap with previously known arc and shell features in the halo light. Simulation tests of our procedure indicate that the probability of finding false groups due to chance is < 1%.

preprint2010arXiv

The Globular Cluster/Central Black Hole Connection in Galaxies

We explore the relation between the total globular cluster population in a galaxy (N_GC) and the the mass of its central black hole (M_BH). Using a sample of 33 galaxies, twice as large as the original sample discussed by Burkert & Tremaine (2010), we find that N_GC for elliptical and spiral galaxies increases in almost precisely direct proportion to M_BH. The S0-type galaxies by contrast do not follow a clear trend, showing large scatter in M_BH at a given N_GC. After accounting for observational measurement uncertainty, we find that the mean relation defined by the E and S galaxies must also have an intrinsic or "cosmic" scatter of +-0.2 in either logN_GC or logM_BH. The residuals from this correlation show no trend with globular cluster specific frequency. We suggest that these two types of galaxy subsystems (central black hole and globular cluster system) may be closely correlated because they both originated at high redshift during the main epoch of hierarchical merging, and both require extremely high-density conditions for formation. Lastly, we note that roughly 10% of the galaxies in our sample (one E, one S, and two S0) deviate strongly from the main trend, all in the sense that their M_BH is at least 10x smaller than would be predicted by the mean relation.

preprint2010arXiv

The Inner Halo of M87: A First Direct View of the Red-Giant Population

An unusually deep (V,I) imaging dataset for the Virgo supergiant M87 with the Hubble Space Telescope ACS successfully resolves its brightest red-giant stars, reaching M_I(lim) = -2.5. After assessing the photometric completeness and biasses, we use this material to estimate the metallicity distribution for the inner halo of M87, finding that the distribution is very broad and likely to peak near [m/H] ~ -0.4 and perhaps higher. The shape of the MDF strongly resembles that of the inner halo for the nearby giant E galaxy NGC 5128. As a byproduct of our study, we also obtain a preliminary measurement of the distance to M87 with the TRGB (red-giant branch tip) method; the result is (m-M)_0 = 31.12 +- 0.14 (d = 16.7 +- 0.9 Mpc). Averaging this result with three other recent techniques give a weighted mean d(M87) = (16.4 +- 0.5) Mpc.

preprint2010arXiv

The Kinematics of the Globular Cluster System of NGC 5128 with a New, Large Sample of Radial Velocity Measurements

New radial velocity measurements for previously known and newly confirmed globular clusters (GCs) in the nearby massive galaxy NGC 5128 are presented. We have obtained spectroscopy from LDSS-2/Magellan, VIMOS/VLT, and Hydra/CTIO from which we have measured the radial velocities of 218 known, and identified 155 new, GCs. The current sample of confirmed GCs in NGC 5128 is now 605 with 564 of these having radial velocity measurements. We have performed a new kinematic analysis of the GC system that extends out to 45 arcmin in galactocentric radius. We have examined the systemic velocity, projected rotation amplitude and axis, and the projected velocity dispersion of the GCs as functions of galactocentric distance and metallicity. Our results indicate that the metal-poor GCs have a very mild rotation signature of (26 pm 15) km/s. The metal-rich GCs are rotating with a higher, though still small signature of (43 pm 15) km/s around the isophotal major axis of NGC 5128 within 15 arcmin. Their velocity dispersions are consistent within the uncertainties and the profiles appear flat or declining within 20 arcmin. We note the small sample of metal-rich GCs with ages less than 5 Gyr in the literature appear to have different kinematic properties than the old, metal-rich GC subpopulation. The mass and mass-to-light ratios have also been estimated using the GCs as tracer particles for NGC 5128. Out to a distance of 20 arcmin, we have obtained a mass of (5.9 pm 2.0) x 10^(11) M_(sun) and a mass-to-light ratio in the B-band of 16 M_(sun)/L_(B,sun). Combined with previous work on the ages and metallicities of its GCs, as well as properties of its stellar halo, our findings suggest NGC 5128 formed via hierarchical merging over other methods of formation, such as major merging at late times.

preprint2010arXiv

Ultra-Compact Dwarfs in the Core of the Coma Cluster

We have discovered both a red and a blue subpopulation of Ultra-Compact Dwarf (UCD) galaxy candidates in the Coma galaxy cluster. We analyzed deep F475W (Sloan g) and F814W (I) Hubble Space Telescope images obtained with the Advanced Camera for Surveys Wide Field Channel as part of the Coma Cluster Treasury Survey and have fitted the light profiles of ~5000 point-like sources in the vicinity of NGC 4874, one of the two central dominant galaxies of the Coma cluster. Although almost all of these sources are globular clusters that remain unresolved, we found that 52 objects have effective radii between ~10 and 66 pc, in the range spanned by Dwarf Globular Transition Objects (DGTO) and UCDs. Of these 52 compact objects, 25 are brighter than M_V ~-11 mag, a magnitude conventionally thought to separate UCDs and globular clusters. The UCD/DGTO candidates have the same color and luminosity distribution as the most luminous globular clusters within the red and blue subpopulations of the immensely rich NGC 4874 globular cluster system. Unlike standard globular clusters, blue and red UCD/DGTO subpopulations have the same median effective radius. The spatial distribution of UCD/DGTO candidates reveal that they congregate towards NGC 4874, and are not uniformly distributed. We find a relative deficit of UCD/DGTOs compared with globular clusters in the inner 15 kpc around NGC 4874, however at larger radii UCD/DGTO and globular clusters follow the same spatial distribution.

preprint2009arXiv

Further Definition of the Mass-Metallicity Relation in Globular Cluster Systems Around Brightest Cluster Galaxies

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

preprint2009arXiv

Massive star clusters in galaxies

The ensemble of all star clusters in a galaxy constitutes its star cluster system. In this review, the focus of the discussion is on the ability of star clusters, particularly the systems of old massive globular clusters (GCSs), to mark the early evolutionary history of galaxies. I review current themes and key findings in GCS research, and highlight some of the outstanding questions that are emerging from recent work.

preprint2009arXiv

Statistical Tools for Classifying Galaxy Group Dynamics

The dynamical state of galaxy groups at intermediate redshifts can provide information about the growth of structure in the universe. We examine three goodness-of-fit tests, the Anderson--Darling (A-D), Kolmogorov and chi-squared tests, in order to determine which statistical tool is best able to distinguish between groups that are relaxed and those that are dynamically complex. We perform Monte Carlo simulations of these three tests and show that the chi-squared test is profoundly unreliable for groups with fewer than 30 members. Power studies of the Kolmogorov and A-D tests are conducted to test their robustness for various sample sizes. We then apply these tests to a sample of the second Canadian Network for Observational Cosmology Redshift Survey (CNOC2) galaxy groups and find that the A-D test is far more reliable and powerful at detecting real departures from an underlying Gaussian distribution than the more commonly used chi-squared and Kolmogorov tests. We use this statistic to classify a sample of the CNOC2 groups and find that 34 of 106 groups are inconsistent with an underlying Gaussian velocity distribution, and thus do not appear relaxed. In addition, we compute velocity dispersion profiles (VDPs) for all groups with more than 20 members and compare the overall features of the Gaussian and non-Gaussian groups, finding that the VDPs of the non-Gaussian groups are distinct from those classified as Gaussian.

preprint2009arXiv

Stochastic self-enrichment, pre-enrichment, and the formation of globular clusters

We develop a model for stochastic pre-enrichment and self-enrichment in globular clusters (GCs) during their formation process. GCs beginning their formation have an initial metallicity determined by the pre-enrichment of their surrounding protocloud, but can also undergo internal self-enrichment during formation. Stochastic variations in metallicity arise because of the finite numbers of supernova. We construct an analytic formulation of the combined effects of pre-enrichment and self-enrichment and use Monte Carlo models to verify that the model accurately encapsulates the mean metallicity and metallicity spread among real GCs. The predicted metallicity spread due to self-enrichment alone, a robust prediction of the model, is much smaller than the observed spread among real GCs. This result rules out self-enrichment as a significant contributor to the metal content in most GCs, leaving pre-enrichment as the viable alternative. Self-enrichment can, however, be important for clusters with masses well above 10^6 Msun, which are massive enough to hold in a significant fraction of their SN ejecta even without any external pressure confinement. This transition point corresponds well to the mass at which a mass-metallicity relationship ("blue tilt") appears in the metal-poor cluster sequence in many large galaxies. We therefore suggest that self-enrichment is the primary driver for the mass-metallicity relation. Other predictions from our model are that the cluster-to-cluster metallicity spread decreases amongst the highest mass clusters; and that the red GC sequence should also display a more modest mass-metallicity trend if it can be traced to similarly high mass.

preprint2009arXiv

Structural Parameters of the M87 Globular Clusters

We derive structural parameters for ~2000 globular clusters in the giant Virgo elliptical M87 using extremely deep Hubble Space Telescope images in F606W (V) and F814W (I) taken with the ACS/WFC. The cluster scale sizes (half-light radii r_h) and ellipticities are determined from PSF-convolved King-model profile fitting. We find that the r_h distribution closely resembles the inner Milky Way clusters, peaking at r_h~2.5 pc and with virtually no clusters more compact than r_h ~ 1 pc. The metal-poor clusters have on average an r_h 24% larger than the metal-rich ones. The cluster scale size shows a gradual and noticeable increase with galactocentric distance. Clusters are very slightly larger in the bluer waveband V a possible hint that we may be beginning to see the effects of mass segregation within the clusters. We also derived a color magnitude diagram for the M87 globular cluster system which show a striking bimodal distribution.

preprint2009arXiv

The Ages, Metallicities and Alpha Element Enhancements of Globular Clusters in the Elliptical NGC 5128: A Homogeneous Spectroscopic Study with Gemini/GMOS

We present new integrated light spectroscopy of globular clusters (GCs) in NGC 5128 in order to measure radial velocities and derive ages, metallicities, and alpha-element abundance ratios. Using Gemini-S 8-m/GMOS, we obtained spectroscopy in the range of ~3400-5700 AA for 72 GCs with S/N > 30 /AA and we have also discovered 35 new GCs within NGC 5128 from our radial velocity measurements. We measured and compared the Lick indices from HdeltaA through Fe5406 with the single stellar population (SSP) models of Thomas et al.(2003,2004). We also measure Lick indices for 41 Milky Way GCs from Puzia et al. (2002) and Schiavon et al. (2005) with the same methodology for direct comparison. Our results show that 68% of the NGC 5128 GCs have old ages (> 8 Gyr), 14% have intermediate ages (5-8 Gyr), and 18% have young ages (< 5 Gyr). However, when we look at the metallicity of the GCs as a function of age, we find 92% of metal-poor GCs and 56% of metal-rich GCs in NGC 5128 have ages > 8 Gyr, indicating that the majority of both metallicity subpopulations of GCs formed early, with a significant population of young and metal-rich GCs forming later. Our metallicity distribution function generated directly from spectroscopic Lick indices is clearly bimodal, as is the color distribution of the same set of GCs. Thus the metallicity bimodality is real and not an artifact of the color to metallicity conversion. The [alpha/Fe] values are supersolar with a mean value of 0.14pm0.04, indicating a fast formation timescale. However, the GCs in NGC 5128 are not as [alpha/Fe] enhanced as the Milky Way GCs also examined in this study. Our results support a rapid, early formation of the GC system in NGC 5128, with subsequent major accretion and/or GC and star forming events in more recent times (abridged).

preprint2009arXiv

The Color-Magnitude Relation for Metal-Poor Globular Clusters in M87: Confirmation From Deep HST/ACS Imaging

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 ($γ_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.

preprint2009arXiv

The Distance to NGC 5128 (Centaurus A)

In this paper we review the various high precision methods that are now available to determine the distance to NGC 5128. These methods include: Cepheids, TRGB (tip of the red giant branch), PNLF (planetary nebula luminosity function), SBF (surface brightness fluctuations) and Long Period Variable (LPV) Mira stars. From an evaluation of these methods and their uncertainties, we derive a best-estimate distance of 3.8 +- 0.1 Mpc to NGC 5128 and find that this mean is now well supported by the current data. We also discuss the role of NGC 5128 more generally for the extragalactic distance scale as a testbed for the most direct possible comparison among these key methods.

preprint2009arXiv

The Globular Cluster System in M87: A Wide-Field Study with CFHT/Megacam

CFHT Megacam data in (g',r',i') are used to obtain deep, wide-field photometry of the globular cluster system (GCS) around M87. A total of 6200 GCs brighter than M_I = -8.5 are included in the study, essentially containing almost the entire bright half of the total GC population in the galaxy. The classic bimodal metal-poor and metal-rich sequences of GCs show up clearly; while the spatial distribution of the GCs can be traced detectably outward to more than 100 kpc, the blue, metal-poor subpopulation is very much more spatially extended than the red subpopulation. Both the red and blue GC subsystems have radial metallicity gradients, where mean heavy-element abundance scales with projected galactocentric distance as Z ~ R^{-0.12} (blue) and R^{-0.17} (red). The blue sequence exhibits a strongly significant mass/metallicity relation (MMR) in which the mean metallicity gradually increases with cluster luminosity as Z ~ L^{0.25 \pm 0.05} for the luminosity range M_I < -10. However, this relation is also clearly nonlinear, becoming vertical at fainter levels. The MMR can be understood as the result of self-enrichment within the most massive metal-poor GCs during their formation. In giant elliptical galaxies, the red GCs are likely to represent a broad composite population formed during several major starbursts. If so, the red sequence might display a population-based MMR that could in principle be either positive or negative.

preprint2009arXiv

The remnants of galaxy formation from a panoramic survey of the region around M31

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

preprint2007arXiv

Wide-Field Survey of Globular Clusters in M31. I. A Catalog of New Clusters

We present the result of a wide-field survey of globular clusters (GCs) in M31 covering a 3deg x 3deg field c. We have searched for GCs on CCD images taken with Washington CMT1 filters at the KPNO 0.9 m telescope using steps: (1) inspection of morphological parameters given by the SExtractor package such as stellarity, full maximum, and ellipticity; (2) consulting the spectral types and radial velocities obtained from spectra takena spectrograph at the WIYN 3.5 m telescope; and (3) visual inspection of the images of each object. We have and GC candidates, of which 605 are newly found GCs and GC candidates and 559 are previously known GCs. Amoects there are 113 genuine GCs, 258 probable GCs, and 234 possible GCs, according to our classification critee known objects there are 383 genuine GCs, 109 probable GCs, and 67 possible GCs. In total there are 496 genprobable GCs and 301 possible GCs. Most of these newly found GCs have T1 magnitudes of 17.5 - 19.5 mag, [17.9 < V < 19.9 mag assuming (C-T1) ~ 1.5], and (C-T1) colors in the range 1 - 2.