We present new Hubble Space Telescope (HST)-NIC3, near-infrared H-band photometry of globular clusters (GC) around NGC 4365 and NGC 1399 in combination with archival HST-WCPC2 and ACS optical data. We find that NGC 4365 has a number of globular clusters with bluer optical colors than expected for their red optical to near-infrared colors and an old age. The only known way to explain these colors is with a significant population of intermediate-age (2-8 Gyr) clusters in this elliptical galaxy. In contrast, NGC 1399 reveals no such population. Our result for NGC 1399 is in agreement with previous spectroscopic work that suggests that its clusters have a large metallicity spread and are nearly all old. In the literature, there are various results from spectroscopic studies of modest samples of NGC 4365 globular clusters. The spectroscopic data allow for either the presence or absence of a significant population of intermediate-age clusters, given the index uncertainties indicated by comparing objects in common between these studies and the few spectroscopic candidates with optical to near-IR colors indicative of intermediate ages. Our new near-IR data of the NGC 4365 GC system with much higher signal-to-noise agrees well with earlier published photometry and both give strong evidence of a significant intermediate-age component. The agreement between the photometric and spectroscopic results for NGC 1399 and other systems lends further confidence to this conclusion, and to the effectiveness of the near-IR technique.

We study the dynamical evolution of the M87 globular cluster system (GCS) with a number of numerical simulations. We explore a range of different initial conditions for the GCS mass function (GCMF), for the GCS spatial distribution and for the GCS velocity distribution. We confirm that an initial power-law GCMF like that observed in young cluster systems can be readily transformed through dynamical processes into a bell-shaped GCMF. However,only models with initial velocity distributions characterized by a strong radial anisotropy increasing with the galactocentric distance are able to reproduce the observed constancy of the GCMF at all radii.We show that such strongly radial orbital distributions are inconsistent with the observed kinematics of the M87 GCS. The evolution of models with a bell-shaped GCMF with a turnover similar to that currently observed in old GCS is also investigated. We show that models with this initial GCMF can satisfy all the observational constraints currently available on the GCS spatial distribution,the GCS velocity distribution and on the GCMF properties.In particular these models successfully reproduce both the lack of a radial gradient of the GCS mean mass recently found in an analysis of HST images of M87 at multiple locations, and the observed kinematics of the M87 GCS.Our simulations also show that evolutionary processes significantly affect the initial GCS properties by leading to the disruption of many clusters and changing the masses of those which survive.The preferential disruption of inner clusters flattens the initial GCS number density profile and it can explain the rising specific frequency with radius; we show that the inner flattening observed in the M87 GCS spatial distribution can be the result of the effects of dynamical evolution on an initially steep density profile. (abridged)

The discovery of young globular clusters in merging galaxies and other environments provides an opportunity to study directly the process of globular cluster formation. Empirically it appears that globular cluster formation occurs preferentially in regions in which star formation occurs at a high rate and efficiency. Further, the interstellar medium in such regions is likely to be at a higher pressure than less active star-forming environments. An additional observational clue to the globular cluster formation process is that young globular clusters have little or no mass-radius relationship. In this paper I argue that high pressure and high star-formation efficiency are responsible for current globular cluster formation. I suggest that the precursors to globular clusters are molecular clouds and that the mass-radius relationship exhibited by such clouds is wiped out by a variable star formation efficiency.

We present a preliminary analysis of HST-WFPC2 observations of globular cluster systems in the two brightest galaxies, UGC 9799 (cD) and NGC 1129 (non-cD), located in the center of rich clusters.

We explore the constraints on globular cluster formation provided by the observed conditions in starbursts where globulars are currently forming, and by the observed properties of young and old globular clusters. We note that the pressure in the ISM of starbursts and mergers implies that molecular clouds in these environments have radii similar to those of globular clusters. Such molecular clouds are therefore viable precursors to globular clusters if the star formation efficiency in the clouds is high. A high star formation efficiency may be a consequence of the high density and associated high binding energy and short dynamical timescale of molecular clouds in such environments. We also note that the apparent lack of a mass-radius relationship in young and old globular cluster systems places important constraints on globular cluster formation models. This is because molecular clouds are observed to follow a virial scaling relation between mass and radius. We suggest that a variable star formation efficiency may weaken or eliminate the mass-radius relation of molecular clouds as they fragment to form globular clusters. We attribute the absence of young globular clusters in the disks of ordinary galaxies such as the Milky Way to the relatively low ambient pressures in such systems.

We present radial velocities for 144 globular clusters around the Virgo gE NGC 4472 (M49), and ages and metallicities from co-added spectra of 131 clusters. We confirm our earlier finding that the metal-poor clusters have a significantly higher velocity dispersion than the metal-rich clusters, and we find little or no rotation in the metal-rich clusters. The velocity dispersion profile is consistent with isotropic orbits for the clusters and the NGC 4472 mass distribution inferred from X-ray data. Our sample of clusters spans a metallicity range of -1.6 < [Fe/H] < 0 dex. The metal-poor and metal-rich cluster populations are coeval within the errors, and all clusters are older than 6 Gyr at 95% confidence. Our findings are consistent with a merger origin for NGC 4472, but other elliptical formation models cannot be ruled out.

We present new radial velocities for 87 globular clusters around the elliptical galaxy NGC 4472, and combine these with our previously published data to create a data set of velocities for 144 globular clusters around NGC 4472. We utilize this data set to analyze the kinematics of the NGC 4472 globular cluster system. The new data confirms our previous discovery that the metal-poor clusters have significantly higher velocity dispersion than the metal-rich clusters in NGC 4472. We also find very little or no rotation in the more spatially concentrated metal-rich population, with a resulting upper limit for this population of (v/sigma)_{proj} < 0.34 at a 99% confidence level. The very small angular momentum in the metal-rich population requires efficient angular momentum transport during the formation of this population which is spatially concentrated and chemically enriched. Such angular momentum transfer can be provided by galaxy mergers, but has not been achieved in other extant models of elliptical galaxy formation that include dark matter halos. We also calculate the velocity dispersion as a function of radius, and show that it is consistent with roughly isotropic orbits for the clusters and the mass distribution of NGC 4472 inferred from X-ray observations of the hot gas around the galaxy.

We have derived ages and metallicities from co-added spectra of 131 globular clusters associated with the giant elliptical galaxy NGC 4472. Based upon a calibration with Galactic globular clusters, we find that our sample of globular clusters in NGC 4472 span a metallicity range of approximately -1.6 <= [Fe/H] <= 0 dex. There is evidence of a radial metallicity gradient in the globular cluster system, which is steeper than that seen in the underlying starlight. Determination of the absolute ages of the globular clusters is uncertain, but formally, the metal-poor population of globular clusters has an age of 14.5 +/- 4 Gyr and the metal-rich population is 13.8 +/- 6 Gyr old. Monte Carlo simulations indicate that the globular cluster populations present in these data are older than 6 Gyr at the 95% confidence level. We find that within the uncertainties, the globular clusters are old and coeval, implying that the bimodality seen in the broadband colours primarily reflects metallicity and not age differences.

We present images of the recent galaxy merger NGC 3256 obtained with the WFPC2 of the Hubble Space Telescope in B and I filters. We show that there is a large population of more than 1,000 compact, bright blue objects in this galaxy within the 7 kpc x 7 kpc region studied. These objects have sizes, colors and luminosities like those expected for young Galactic globular clusters with ages ranging from a few to several hundred Myr. On this basis, we identify at least some fraction of the compact, bright blue objects in NGC 3256 as young globular clusters. The young cluster system comprises 15-20% of the luminosity of the galaxy in B and half that in I within the region studied, indicating a high efficiency of cluster formation on a galaxy-wide scale. We find that the intrinsic color distribution is broad, and there is no significant trend of color with magnitude. The combination of the broad range of observed colors and the lack of a trend of redder colors at fainter magnitudes can not be fit solely by a broad age distribution and/or differential reddening, although the latter is clearly present. The observations can be accounted for by either the preferential depletion/destruction of lower mass clusters as they age, or a very young age (<= 20 Myr) for the cluster population, comparable to or less than the dynamical time of the region in which the clusters are observed. We also find that the luminosity function of the young cluster system can be roughly fit by a power-law with an exponent of -1.8, with tentative evidence that it flattens at faint magnitudes. The clusters are compact in size, with typical estimated half-light radii of 5-10 parsecs, but there is no obvious cut-off for larger radii, and only a shallow trend of size with luminosity. (slightly abridged).

We present a study of the globular cluster systems of two edge-on spiral galaxies, NGC4565 and NGC5907, from WFPC2 images in the F450W and F814W filters. The globular cluster systems of both galaxies appear to be similar to the Galactic globular cluster system. In particular, we derive total numbers of globular clusters of N_{GC}(4565)= 204+/-38 {+87}{-53} and N_{GC}(5907)=170+/-41 {+47}{-72} (where the first are statistical, the second potential systematic errors) for NGC4565 and NGC5907, respectively. This determination is based on a comparison to the Milky Way system, for which we adopt a total number of globular clusters of 180+/-20. The specific frequency of both galaxies is S_N~0.6: indistinguishable from the value for the Milky Way. The similarity in the globular cluster systems of the two galaxies is noteworthy since they have significantly different thick disks and bulge-to-disk ratios. This would suggest that these two components do not play a major role in the building up of a globular cluster system around late-type galaxies.

1057 globular cluster candidates have been identified in a WFPC2 image of the inner region of M87. The Globular Cluster Luminosity Function (GCLF) can be well fit by a Gaussian profile with a mean value of m_V^0=23.67 +/- 0.07 mag and sigma=1.39 +/- 0.06 mag (compared to m_V^0=23.74 mag and sigma=1.44 mag from an earlier study using the same data by Whitmore it et al. 1995). The GCLF in five radial bins is found to be statistically the same at all points, showing no clear evidence of dynamical destruction processes based on the luminosity function (LF), in contradiction to the claim by Gnedin (1997). Similarly, there is no obvious correlation between the half light radius of the clusters and the galactocentric distance. The core radius of the globular cluster density distribution is R_c=56'', considerably larger than the core of the stellar component (R_c=6.8''). The mean color of the cluster candidates is V-I=1.09 mag which corresponds to an average metallicity of Fe/H = -0.74 dex. The color distribution is bimodal everywhere, with a blue peak at V-I=0.95 mag and a red peak at V-I=1.20 mag. The red population is only 0.1 magnitude bluer than the underlying galaxy, indicating that these clusters formed late in the metal enrichment history of the galaxy and were possibly created in a burst of star/cluster formation 3-6 Gyr after the blue population. We also find that both the red and the blue cluster distributions have a more elliptical shape (Hubble type E3.5) than the nearly spherical galaxy. The average half light radius of the clusters is ~2.5 pc which is comparable to the 3 pc average effective radius of the Milky Way clusters, though the red candidates are ~20% smaller than the blue ones.

Optical multi-slit spectra have been obtained for 47 globular clusters surrounding the brightest Virgo elliptical NGC 4472 (M49). Including data from the literature, we analyze velocities for a total sample of 57 clusters and present the first tentative evidence for kinematic differences between the red and blue cluster populations which make up the bimodal colour distribution of this galaxy. The redder clusters are more centrally concentrated and have a velocity dispersion of 240 km/sec compared with 320 km/sec for the blue clusters. The origin of this difference appears to be a larger component of systematic rotation in the blue cluster system. The larger rotation in the more extended blue cluster system is indicative of efficient angular momentum transport, as provided by galaxy mergers. Masses estimated from the globular cluster velocities are consistent with the mass distribution estimated from X-ray data, and indicate that the M/L rises to 50 times solar at 2.5 R_e.

We review the observed properties of globular cluster systems and their implications for models of galaxy formation. Observations show that globular clusters form in gas-rich mergers, and that bimodal metallicity distributions are common in the globular cluster systems of ellipticals, with the metal-poor population more extended than the metal-rich one. These are three of the four predictions of the simple merger model of Ashman & Zepf (1992). The fourth prediction concerns the properties of the globular cluster systems of spirals, and is still to be tested by observation. Adopting Occam's razor, the confirmation of the fundamental predictions of the merger model from both young and old globular cluster systems is strong evidence that typical elliptical galaxies formed from the mergers of spiral galaxies. However, the simplifying assumptions of the Ashman-Zepf merger model limit its applicability to certain complex situations such as the formation of cD galaxies. We conclude this review by introducing new observational and theoretical programs that will further the understanding of the physical mechanisms of globular cluster and galaxy formation.

A set of 76 open clusters with abundances based upon DDO photometry and/or moderate dispersion spectroscopy has been transformed to a common [Fe/H] scale and used to study the local structure and evolution of the galactic disk. The metallicity distribution of clusters with R_GC is best described by two distinct zones. Between R_GC = 6.5 and 10 kpc, the distribution has a mean [Fe/H] = 0.0 and a dispersion of 0.1 dex; there is only weak evidence for a shallow abundance gradient over this distance range. Beyond R_GC = 10 kpc, the metallicity distribution has a dispersion between 0.10 and 0.15 dex, but with a mean [Fe/H] = -0.3, implying a sharp discontinuity at R_GC = 10 kpc. After correcting for the discontinuity, no evidence is found for a gradient perpendicular to the plane. Adopting the clusters interior to 10 kpc as a representative sample of the galactic disk over the last 7 Gyr, the cluster metallicity range is found to be about half that of the field stars. When coupled with the discontinuity in the galactocentric gradient, the discrepancy in the metallicity distribution is interpreted as an indication of significant diffusion of field stars into the solar neighborhood from beyond 10 kpc. These results imply that the sun is NOT atypical of the stars formed in the solar circle 4.6 Gyr ago. It is suggested that the discontinuity reflects the edge of the initial galactic disk as defined by the disk globular cluster system and the so-called thick disk; the initial offset in [Fe/H] created by the differences in the chemical history on either side of the discontinuity has carried through to the current stage of galactic evolution. If correct, diffusion coupled with the absence of an abundance gradient could make the separation of field stars on the basis of galactocentric origin difficult.

We have obtained spectra for globular cluster candidates in M104 with LDSS-2 on the William Herschel Telescope, confirming 34 objects as M104 globular clusters. We find a cluster velocity dispersion of $\sim$ 260 km/sec, and the Projected Mass Estimator gives a mass of 5.0 (3.5,6.7) $\times 10^{11}$ M$_\odot$ for M104 within a projected radius of $\sim 330^{''}$ (14 kpc for D=8.55 Mpc). Our best estimate for the mass-to-light ratio is M/L$_{V_T}$= 16$^{+5.5}_{-5.0}$ within the same radius. Considering all of the possible sources of uncertainty, we find a lower limit of M/L_V= 5.3, which is larger than the M/L$_V$ found from rotation curve analyses inside 180$^{''}$. We thus conclude that the mass-to-light ratio increases with radius, or in other words that M104 possesses a dark matter halo. There is a marginal detection of rotation in the M104 cluster system at the 92.5% confidence level; larger samples will be needed to investigate this possibility. Interestingly, the M104 globular cluster and planetary nebulae (PNe) kinematics are roughly consistent inside $\sim 100^{''}$. Finally, we find a mean cluster metallicity of [Fe/H] = $-$0.70 $\pm$ 0.3, which is more typical of clusters in gE/cD galaxies than it is of clusters in other spirals.

We investigate the universality of the globular cluster luminosity function (GCLF) and the use of this function as an extragalactic distance indicator. Previous studies have found an offset between GCLF distances and those obtained with other techniques. We show that this offset can be understood in terms of a metallicity effect. Specifically, the globular cluster systems used in distance scale studies have traditionally been those around elliptical galaxies. These systems have higher mean metallicities than the Milky Way globular cluster system. Consequently, the peak of the GCLF in the systems around ellipticals is significantly fainter in $B$ and $V$ than the GCLF peak in the Milky Way. We calculate the shift in the peak of the GCLF relative to the Milky Way globulars in $B$, $V$, $R$, $I$ and $J$ for a range of globular cluster metallicities. Applying these corrections, we find good agreement between GCLF distances and those obtained using the surface brightness fluctuations method. The similarity between metallicity-corrected GCLFs suggests that the underlying mass function of globular cluster systems is remarkably constant from one galaxy to another. Our results allow the GCLF to be employed as an improved distance indicator.

We present spectra of the brightest member of the population of compact blue objects discovered in the peculiar galaxy NGC 1275 by Holtzman et al. using Hubble Space Telescope images. These spectra show strong Balmer absorption lines like those observed in A-type stars, as expected if the object is a young globular cluster. The age estimated from the strength of the Balmer lines is about 0.5 Gyr, although ages ranging from 0.1 Gyr to 0.9 Gyr can not be confidently excluded given current models of stellar populations. If these estimated ages are adopted for the young cluster population of NGC 1275 as a whole, the fading predicted by stellar populations models gives a luminosity function which is consistent with that of the Galactic globular cluster system convolved with the observational selection function for the NGC 1275 system. We also use the equivalent widths of the Mg $b$ and Fe 5270 features to constrain the metallicity of the young cluster. Combining these absorption-line widths with the age estimates from the Balmer lines and stellar population models, we find a metallicity of roughly solar based on the Mg $b$ index and somewhat higher for the Fe 5270 index. The radial velocity of the absorption lines of the cluster spectrum is offset from the emission lines of the galaxy spectrum at the same position by $-130$ \kms, providing further evidence for the identification of the object as a globular cluster, and opening up the future possibility of studying the kinematics of young cluster systems. The discovery of objects with the characteristics of young globular clusters in NGC 1275, which shows evidence of a recent interaction or merger, supports the hypothesis that galaxy interactions and mergers are favorable sites for the formation of globular clusters.

We present a study of the colors of globular clusters associated with the elliptical galaxy NGC 3923. Our final sample consists of Washington system C and T_1 photometry for 143 globular cluster candidates with an expected contamination of no more than 10%. We find that the color distribution of the NGC 3923 globular cluster system (GCS) is broad and appears to have at least two peaks. A mixture modeling analysis of the color distribution indicates that a two component model is favored over a single component one at a high level of confidence (>99%). This evidence for more than one population in the GCS of NGC 3923 is similar to that previously noted for the four other elliptical galaxies for which similar data have been published. Furthermore, we find that the NGC 3923 GCS is redder than the GCSs of previously studied elliptical galaxies of similar luminosity. The median metallicity inferred from our (C-T_1)_0 colors is [Fe/H]_{med} = -0.56, with an uncertainty of 0.14 dex arising from all sources of uncertainty in the mean color. This is more metal rich than the median metallicity found for the GCS of M87 using the same method, [Fe/H]_{med} = -0.94. Since M87 is more luminous than NGC 3923, this result points to significant scatter about any trend of higher GCS metallicity with increasing galaxy luminosity. We also show that there is a color gradient in the NGC 3923 GCS corresponding to about -0.5 dex in Delta[Fe/H]\Delta(log r). We conclude that the shape of the color distribution of individual GCSs and the variation in mean color among the GCSs of ellipticals are difficult to understand if elliptical galaxies are formed in a single protogalactic

We present new data on the globular cluster system of the elliptical galaxy NGC 3923 which show that it has the most globular clusters per unit luminosity of any non-cluster elliptical yet observed, with $S_N=6.4 \pm1.4$. NGC 3923 is also among the brightest ellipticals outside of a galaxy cluster for which the number of globular clusters has been determined. Our observation of a large number of clusters per unit luminosity (high $S_N$ value) for a bright elliptical in a sparse environment is consistent with the suggestion of Djorgovski \& Santiago that the number of globular clusters is a power law function of the luminosity with an exponent greater than one. We relate this higher specific frequency of globular clusters in more luminous galaxies to other observations which indicate that the physical conditions within elliptical galaxies at the time of their formation were dependent on galaxy mass.

We discuss statistical techniques for detecting and quantifying bimodality in astronomical datasets. We concentrate on the KMM algorithm, which estimates the statistical significance of bimodality in such datasets and objectively partitions data into sub-populations. By simulating bimodal distributions with a range of properties we investigate the sensitivity of KMM to datasets with varying characteristics. Our results facilitate the planning of optimal observing strategies for systems where bimodality is suspected. Mixture-modeling algorithms similar to the KMM algorithm have been used in previous studies to partition the stellar population of the Milky Way into subsystems. We illustrate the broad applicability of KMM by analysing published data on globular cluster metallicity distributions, velocity distributions of galaxies in clusters, and burst durations of gamma-ray sources. PostScript versions of the tables and figures, as well as FORTRAN code for KMM and instructions for its use, are available by anonymous ftp from kula.phsx.ukans.edu.