A&A 448, 881-891 (2006)
DOI: 10.1051/0004-6361:20054177
N. Bastian1,2 - R. P. Saglia3 - P. Goudfrooij4 - M. Kissler-Patig2 - C. Maraston5 - F. Schweizer6 - M. Zoccali7
1 - Department of Physics and Astronomy, University College London, Gower Street, London, WC1E 6BT, UK
2 - European Southern Observatory, Karl-Schwarzschild-Strasse 2,
85748 Garching b. München, Germany
3 -
Max-Planck-Institut für Extraterrestrische
Physik, Giessenbachstrasse, 85748 Garching, Germany
4 -
Space Telescope Science Institute, 3700 San Martin
Drive, Baltimore, MD 21218, USA
5 -
University of Oxford, Denys Wilkinson Building,
Keble Road, Oxford, OX13RH, UK
6 -
Carnegie Observatories, 813 Santa Barbara Str.,
Pasadena, CA 91101-1292, USA
7 -
Pontificia Universidad Católica de Chile,
Departamento de Astronomía y Astrofísica, Av. Vicuña
Mackenna 4860, 782-0436 Macul, Santiago, Chile
Received 8 September 2005 / Accepted 28 October 2005
Abstract
We present high-dispersion spectra of two extremely
massive star clusters in galactic merger remnants,
obtained using the UVES spectrograph mounted on the ESO Very Large
Telescope. One cluster, W30, is located in the 500 Myr
old merger remnant NGC 7252 and has a velocity dispersion and
effective radius of
km s-1 and
pc, respectively.
The other cluster, G114, located in the
3 Gyr old merger
remnant NGC 1316, is much more compact,
pc, and has
a velocity dispersion of
km s-1. These measurements
allow an estimate of the virial mass of the two clusters, yielding
and
.
Both clusters are
extremely massive, being more than three times heavier than the
most massive globular clusters in the Galaxy. For both
clusters we measure light-to-mass ratios, which when compared to
simple stellar population (SSP) models of the appropriate age, are
consistent
with a Kroupa-type stellar mass function. Using measurements from
the literature we find a strong age dependence on how well SSP
models (with underlying Kroupa or Salpeter-type stellar mass
functions) fit the light-to-mass ratio of clusters. Based on this
result we suggest that the large scatter in the light-to-mass ratio
of the youngest clusters is not due to variations in the underlying
stellar mass function, but instead to the rapidly changing
internal dynamics
of young clusters. Based on sampling statistics we argue that while W30 and G114 are
extremely massive, they are consistent with being the most massive
clusters formed in a continuous power-law cluster mass
distribution. Finally, based on the positions of old globular
clusters, young massive
clusters (YMCs), ultra-compact dwarf galaxies (UCDs) and
dwarf-globular transition objects (DGTOs) in
-space we
conclude that 1) UCDs and DGTOs are consistent with the high mass end of star
clusters and 2) YMCs occupy a much larger parameter space than old
globular clusters, consistent with the idea of preferential
disruption of star clusters.
Key words: galaxies: star clusters - galaxies: interactions - galaxies: individual: NGC 1316 - galaxies: individual: NGC 7252
Our concept of star clusters has changed rapidly during the past two decades. The first resolved young clusters with masses comparable to those of the traditional globular clusters (taken with the Hubble Space Telescope Holtzman et al. 1992) confirmed the suggestions of Schweizer (1987) that mergers of galaxies may produce "young'' globular cluster sized objects. These results were rapidly followed by the discovery of additional young massive clusters (YMCs) in other galaxy mergers, as well as in dwarf, starburst, and normal galaxies (see reviews by Whitmore 2003 and Larsen 2004). Even our own galaxy is producing YMCs with comparable masses and sizes to those observed in merging galaxies, e.g. Westerlund 1 (Clark et al. 2005). The apparent ubiquity of these objects has raised the question of how "universal'' their detailed properties are, in particular concerning their formation and subsequent evolution.
In order to address this and other questions regarding YMCs, we have
begun a programme to obtain kinematic and structural properties of
star clusters which lie at the extreme high-end of the distribution of
observed (luminous) masses. Our first result, presented in Maraston et al. (2004) was for the extremely luminous star cluster, W3, in the
galactic merger remnant NGC 7252. Combining the velocity dispersion
measured with UVES on the VLT ( km s-1) with the
size determined from HST images (
)
led to the dynamical mass estimate of
.
This mass and that estimated from photometric
methods (see Maraston et al. 2004, for details) were in excellent
agreement, arguing that the stellar mass function within W3 was
Salpeter-like.
Using similar techniques as were employed in the above work, some studies have suggested that the stellar mass function in star clusters can vary substantially (e.g. Smith & Gallagher 2001), while others have reported standard Kroupa (2002) or Salpeter (1955) type stellar mass functions (e.g. Larsen et al. 2004). These discrepant results have left the question of the variance of the stellar mass function of massive YMCs open to debate.
Additionally, detailed knowledge of the internal dynamics and structural parameters of YMCs has allowed a comparison between them and other gravitationally bound systems. In Maraston et al. (2004) we showed that W3 is too diffuse for its mass when compared with old globular clusters, whereas it is too compact relative to dwarf galaxies. However, we also showed that W3's properties were extremely similar to those of massive point-like objects discovered in Fornax (Hilker et al. 1999). Based on this similarity, Bastian et al. (2005b) have suggested a mechanism which may allow massive star clusters to exist far from the main body of the host galaxy, namely the formation of massive clusters in the tidal debris of galactic interactions/mergers.
Following up on these results, we have obtained high resolution UVES optical spectra of two additional highly luminous star cluster
candidates in galactic merger remnants. The first is the cluster W30,
the second brightest cluster in NGC 7252. W30 has an estimated age of
300-500 Myr and a
metallicity between half solar and solar, estimated from
optical spectra (Schweizer & Seitzer 1998) as well as from optical
and near-infrared photometry (Maraston et al. 2001). The observed
magnitude of W30 is mV(W30)=19.46 mag and it has a (V-I) colour
of 0.63 mag
(Miller et al. 1997). The second cluster in this study is G114 in
NGC 1316. Based on optical and near-infrared photometry along with
near-infrared spectroscopy, Goudfrooij et al. (2001a,b) estimate an
age of
Gyr for G114. This cluster, the brightest one
in NGC 1316, has an observed
magnitude mB(G114)=19.63 mag and (B-I) colour 1.87 mag. Throughout
this work we adopt the distances to NGC 7252 and NGC 1316 which were
used in Maraston et al. (2004) and Goudfrooij et al. (2001a,b)
respectively, namely 64.4 and 22.9 Mpc. Figures 1 and 2 show the HST/WFPC2 planetary camera chip images of NGC 7252 and NGC 1316, respectively, along with the slit sizes and positions used. In
Fig. 1 we also mark the massive star cluster W3.
The luminosities of G114 and W30 (assuming a Salpeter or
Kroupa-type stellar IMF) imply that they have extremely high
masses, more
than 3 times that of the most massive globular cluster in the
Galaxy, Cen which has a mass of
(Meylan & Mayor 1986).
In this work we investigate the structural and kinematic properties of these two massive star clusters in order to determine how they relate to "normal'' young and old globular clusters. In Sect. 2 we measure the effective radii of the two clusters using high-resolution HST imaging, and in Sect. 3 we determine their velocity dispersions. We combine these results in Sect. 4 to estimate the dynamical mass of the clusters. In Sect. 5 we discuss the implications of our results in terms of the underlying stellar mass function of the clusters, and compare their properties with other bound stellar systems. Finally, we summarise the results and present our conclusions in Sect. 6.
![]() |
Figure 1:
F555W image of the centre of NGC 7252 showing the slit
size and position, as well as cluster W3 for reference. The
image is 740 pixels on a side which corresponds to a linear distance of ![]() |
Open with DEXTER |
![]() |
Figure 2:
F450W image of the centre of NGC 1316 showing the slit
position. The image is 800 pixels on a side which corresponds to a linear distance of ![]() |
Open with DEXTER |
Structural parameters for clusters G114 and W30 were measured on Wide Field Planetary Camera 2 (WFPC2) and on Advanced Camera for Surveys (ACS) HST images. The images for W30 are presented in detail in Miller et al. (1997). Here we note that W30 is located on the Planetary Camera chip. For G114 we measured the size on both the WFPC2 (data presented in Goudfrooij et al. 2001b) and ACS (presented in Goudfrooij et al. 2004) images.
Sizes were found using the ISHAPE routine of Larsen (1999). This routine convolves the PSF with a specified model profile of varying sizes and fits it directly to the images. The outputs of this routine for the best fitting model are the FWHM of the major axis, the minor to major axis ratio, and the goodness of fit. For the WFPC2 images we used a PSF generated by TinyTim (Krist & Hook 1997) at the exact location of the cluster on the chip. We used the drizzled ACS images to eliminate the geometric distortion of the camera. The PSF for each filter of the ACS images was constructed using sources from ACS observations of the globular cluster 47-Tuc.
For each cluster we fit multiple profiles to the images in each of the observed filters, including King, Moffat, and Gaussian profiles.
Using the F450W WFPC2 observations, we fit eight profiles to the image of the cluster. The results are shown in Table 1. In addition to the set profiles, we also fit a Moffat profile with a variable index. The best fitting model is one with an index 1.76. Using this profile we measure an effective radius of 4.08 pc, which is remarkably close to the mean of the size determined using all the other profiles (4.10 pc). We estimate the error on the size as the standard deviation of the size measured for the different profiles, 0.25 pc.
We have also measured the size of G114 on F555W and F814WHST-ACS images. In Table 1 we show the
results for the fits using the best fitting profile from the WFPC2
observations. We find that the size determined on the F555WACS images is 15% larger than that found on the
WFPC2 images. However the size measured on the F814WACSimages is
15% smaller than on the WFPC2 images.
As the three images all give approximately the same result we conclude
that G114 is resolved, and we assign the size of
(G114
pc.
Table 1:
Effective radius of NGC 1316:G114 for various models as
determined from the F450W WFPC2 image. See text for details. The
best fitting model (lowest )
is shown in bold.
Table 2:
Effective radius of NGC 7252:W30 for various models as
determined from the F555W WFPC2 image (PC chip). See text for details. The
best fitting model (lowest )
is shown in bold.
We observed NGC 7252:W30 and NGC 1316:G114 with the UltraViolet
Echelle Spectrograph (UVES) mounted on the ESO/VLT on the nights of
Sept. 13-16th, 2004. We used the red arm, CD#3 grating centred on 5200 Å. This resulted in a wavelength coverage from 4200 Å to 6200 Å and a resolution of 5 km s-1 at 5200 Å. The data were reduced and extracted using the on-line UVES
pipeline with the relevant bias and flat-frames. Cosmic-rays were
also removed using the pipeline. After extraction,
each spectrum was corrected to the helio-centric velocity frame, and
summed to create the total spectrum for each cluster. The total
exposure times were 8.67 h on G114 and 25.1 h on W30.
Figures 1 and 2 show the positions and lengths of
the UVES slits for clusters W30 and G114 respectively, superimposed on
HST planetary camera images. We note that the background near
the positions of each of these clusters is devoid of spurious sources, which allowed
a clear background subtraction (determined using a spline function) at
the position of the clusters.
Additionally, we observed several template stars which were used to complement our existing template catalogue (Maraston et al. 2004), namely HR 35 (F4 V), HR 8709 (A4 V), HD 203638 (K0 III), and HD 212574 (A6 V) where the designation in the brackets refers to the spectral type of the star. The stars were reduced in the same way as described above.
The determination of the velocity dispersion of each of the clusters was carried out in exactly the same way as was done for W3, which is described in detail in Maraston et al. (2004). In summary we used an adapted version of the Fourier Correlation Quotient (FCQ, Bender 1990) method as implemented by Bender et al. (1994), using templates chosen to match the stellar populations within each cluster. The templates were chosen to have temperatures and gravities (i.e. luminosity classes) appropriate to stars at the main-sequence turn-off point and giants stars in synthetic stellar populations of the same age as each cluster. Their contributions to the composite template were weighted using the weights predicted by the SSPs at the appropriate age (Maraston 2005, e.g. her Fig. 13).
Table 3: Velocity dispersion measurements of NGC 1316:G114.
The results of the determination of the velocity dispersion for G114
is shown in Table 3. Due to the high S/N ratio of
the data and the large number of metal lines in the optical part of
the spectrum (due to the dominance of cool stars at the cluster's age
of 3 Gyr) we were able
to use the full spectral range to determine the velocity dispersion.
We note that the results do not depend crucially on the assumed stellar
template. The adopted one-dimensional velocity dispersion for G114
is
km s-1, which is the average over the
full wavelength range of Templates 2 and 3, which should be the
closest match to the actual stellar population based on the models of
Maraston (2005). Figure 3 shows the blue section of the
observed spectrum of G114
(black), the best fitting broadened stellar
template (red) and the difference between the two (green). All
spectra shown in this work have been divided by the continuum and had
a value of one subtracted from them, to place the average value at zero.
We have also measured the heliocentric line-of-sight velocity of G114,
![]() |
Figure 3:
The observed (observed divided by the continuum minus 1
and smoothed by 15 pixels) spectrum of G114
is plotted in
black. The red shows the best fitting template, broadened by 42.1 km s-1. The green is the difference between the two, shifted
downward by a constant for clarity. This is just one
region that we fitted, Table 3 shows all the
regions that were used. The large scale undulations
(![]() |
Open with DEXTER |
Based on comparisons between the optical/near-infrared colours (Miller
et al. 1997; Maraston et al. 2001) and optical spectroscopy (Schweizer
& Seitzer 1998), W30 and W3 appear to have very similar ages and
metallicities. Because of this, we have used the same stellar
template to determine the velocity dispersion of W30 as we used for W3
(Maraston et al. 2004). Similarly as was done for W3, we have limited our analysis to
the region red-ward of H,
in particular concentrating on the
region around the Mg triplet (
Å)
and Fe lines (
Å). Table 4 shows the results of the fitting
on specific regions of the spectrum. We adopt the value
km s-1 which is the average between the value found for
fitting solely on the Mg lines and fitting on the full region of
interest (between 5150 Å and 5350 Å).
Figure 4 shows the spectrum of cluster W30 in the fitting region (black), the best fitting broadened template (red), and the difference between the two (green).
The measured line of sight velocity of W30 is
![]() |
Figure 4: The observed (smoothed by 5 pixels) spectrum of W30 is plotted in black. The red shows the best fitting template, broadened by 28.0 km s-1. The green in the difference between the two, shifted downward by a constant for clarity. The top two panels show the entire wavelength range which was used in the fitting, while the bottom figure shows a blow up of the region including the Mg triplet, which are the strongest lines in this region. |
Open with DEXTER |
Assuming that the clusters are in virial equilibrium, we can estimate
their virial masses using the relation
Table 4: Derived velocity dispersion for different portions of the spectra of W30 in NGC 7252. All values are given in km s-1.
We note, however, that the value of the constant in
Eq. (1) known as ,
may
change as a function of age of the cluster (Boily et al. 2005). This
effect will be most dramatic in the first 30 Myr of a cluster's
lifetime, and the size of the effect depends on the surface mass density of
the cluster in the sense that the clusters with the highest densities
will be the most affected. We will return to this point in
Sect. 5.1.1.
These results confirm the extremely large masses of G114 and W30.
However, it is important to put these clusters in the context of the
mass functions of the full cluster systems of their respective galaxies
in order to test to what extent they really are outliers.
Table 5: The properties of the massive star clusters.
To do this we have performed a series of monte-carlo tests of the
mass function of cluster populations. We assume an initial
power-law mass distribution of the clusters within each galaxy of the form
,
with
(e.g. Miller et al. 1997) and also that this distribution gets filled
randomly. In the case of the NGC 7252 system, W3 (the most massive
cluster in the NGC 7252) is
5 times
more massive than W30 (the second most massive cluster in this
system). Under these conditions, we expect that
20% of the
realizations of the
cluster populations will have a factor of five or greater between the
most massive and the second most massive clusters within the
system
.
The third brightest cluster in NGC 7252 system (W6) is
0.2 mag
fainter than W30, corresponding to a mass difference of only
17%
(assuming a common age). Hence, we see that while W3 and W30 are
extremely massive clusters, they are compatible with being the most massive
clusters of a continuous power-law distribution, which we note
continues to the detection limit.
The intermediate aged cluster population of NGC 1316 can also be
readily explained by the same argument as above, as the second and
third most massive clusters in this system (G114 is the most massive)
are only 0.47 and 0.54 mag fainter respectively. This
corresponds to less than a factor of two in the luminosity (and mass
assuming that the clusters have similar ages and metallicities). A
difference of two or greater in the ratio of the most massive and second most
massive clusters was found in 50% of the realizations of a
cluster population. Goudfrooij et al. (2004) have reported that the bright end of the luminosity
function (which we assume to represent the mass function) is well
approximated by a power-law of the type used above.
Since NGC 7252:W3, NGC 7252:W30, and NGC 1316:G114 can be readily understood through sampling statistics we will assume that they are simply the most massive clusters of a continuous distribution. In Sect. 5.2 we will use this assumption and the detailed properties of these clusters to understand more enigmatic objects, namely the dwarf galaxy transition objects (DGTOs) and ultra-compact dwarf galaxies (UCDs).
A comparison between the mass of a cluster derived using photometric methods and the mass determined through kinematic arguments, allows an independent test of the assumptions that went into each estimate. The assumption that has garnered the most attention in recent years, is that of the underlying stellar mass function (MF) of the cluster. In order to estimate the photometric mass of a cluster, mass-to-light ratios from simple stellar population (SSP) models are required. These, in turn, are heavily dependent on the assumed stellar mass function. Therefore, assuming that all the other assumptions are valid (such as the state of equilibrium, correct extinction determination, and stellar evolutionary tracks) any difference between the mass of a cluster derived in these two ways is caused by the difference between the input stellar mass function and that of the cluster.
Studies that have used this technique have reported strongly divergent
results. For example, Smith & Gallagher (2001) and McCrady et al. (2005) have
reported that the YMC M82F is deficient in low-mass stars, relative to
the standard Salpeter (1955) (a single power-law mass distribution from the lower to
the upper mass limit) or Kroupa (2002) type MFs (a single power-law above
1
and significantly flatter below this limit). However, Maraston et al. (2004), Larsen et al. (2004), and Larsen & Richtler (2004) have
shown that clusters in a wide variety of galactic environments are
consistent with a Salpeter or Kroupa-type MF.
We therefore carry out this experiment for the two massive clusters W30 and G114. In Fig. 5 we show the light-to-mass ratio from the Maraston (2005) models for solar metallicity and a Salpeter (dashed line) & Kroupa (solid line) stellar mass functions. Over-plotted in red are the two clusters in the present study as well as W3 from our previous study. These three clusters all lie impressively close to the value using a Kroupa MF. Hence they are most likely not deficient in low mass stars.
![]() |
Figure 5:
The derived light-to-mass (V-band) ratios as a
function of the age of the clusters. Over-plotted are the
![]() ![]() ![]() |
Open with DEXTER |
In order to compare our results to other young clusters, we have taken a sample of clusters with velocity dispersion and radii measurements from the literature. The clusters, their parameters, and the corresponding references are listed in Table 6. We have taken the fundamental parameters (age, extinction, brightness, velocity dispersion, radius, and distance modulus) directly from the given reference. In some cases the V-band magnitude was not given, although we note that all clusters with ages greater than 20 Myr have observed V-band magnitudes. These older clusters will constitute the main part of our comparison. In those cases where V-band magnitudes were not available we transformed the given magnitude to the V-band using the colours in the Maraston (2005) SSP models at the appropriate age (which assume a Salpeter IMF).
We have estimated the mass of each cluster using Eq. (1)
and used their V-band magnitudes (and given ages) to place them in
Fig. 5 (blue points). From this figure it is clear that
the amount of deviation from standard stellar mass functions (Kroupa
or Salpeter-type) is heavily age dependent, with the older clusters
(with the exception of M82F)
all consistent with a Kroupa or Salpeter-IMF and the youngest clusters
showing a large amount of scatter. Below we outline three possible
explanations for this.
A first possibility for the age-dependent scatter in
Fig. 5 is that
(the parameter
in the numerator of Eq. (1)) changes as a
function of time (Boily et al. 2005). This is caused by
mass-segregation in young clusters and further internal dynamical
evolution of the cluster. The variation of
is expected to also
be heavily dependent on the surface density of the star cluster, with
higher surface densities leading to larger variations of
.
Figure 6 shows the mean surface density within the
half-light radius (using the
estimated virial mass of the cluster) of star clusters vs. their
measured velocity dispersions. The small blue filled triangles represent old globular
clusters in our galaxy (McLaughlin & van der Marel 2005), the small
green filled circles
are globular clusters in M 31 (a collation of data from McLaughlin &
van der Marel in prep.), the small upside-down
magenta triangles are old GCs in NGC 5128 (Martini & Ho 2004), and the
large red circles are young massive star
clusters in a variety of galaxies (listed in
Table 6). All of the YMCs in
Fig. 6 have surface densities above
pc-2 and most are above
pc-2,
which is the regime where
is expected to vary strongly (Boily
et al. 2005).
Table 6: The properties of young massive clusters taken from the literature.
A second possibility is that many of the youngest clusters
are not in dynamical equilibrium. This could be due to external
gravitational effects (e.g. close passages to massive GMCs). As
clusters are born in gas-rich environments this is a likely possibility.
A lack of equilibrium could also be caused by the rapid expulsion of
the gas left over from the star formation process (assuming a
non-100% star formation efficiency). This can have a severe
influence on a young cluster (e.g. Boily & Kroupa 2003). Bastian et al. (2005a) have suggested that rapid gas removal may be responsible for the
dissolution of 70-90% of clusters within the first 10 Myr of their
lives, independent of cluster mass. The lack of dynamical
equilibrium is also supported by the work of de Grijs et al. (2005) who showed that the clusters which deviate the most from
the old globular cluster MV vs. log
relation are found
in the highest density environments, and hence are the most likely to
be affected by interactions with the external environment.
Finally, a third possibility for the observed age-dependent scatter in
Fig. 5 is that only clusters with Kroupa
or Salpeter-type stellar mass functions survive for more than 100 Myr. Smith & Gallagher (2001) suggest that if the
cluster M82 F has a significantly top-heavy stellar IMF
(i.e. truncated below 2-3
)
it will lack the gravitational
potential to remain bound due to stellar evolutionary mass loss after
2-3 Gyr. However, in Fig. 5 the clusters which deviate
the most from Salpter or Kroupa-type stellar IMFs have light-to-mass
ratios below the expected value. This implies that they are
over-abundant in low mass stars relative to Salpeter or Kroupa-type
IMFs. Since low mass, long-lived stars provide the gravitational
potential to keep a cluster bound, we would expect these clusters to
be long-lived and hence to see old clusters
which also have light-to-mass ratios below that expected for standard
IMFs. Such clusters are clearly lacking in Fig. 5
arguing against this possibility.
![]() |
Figure 6:
The measured velocity dispersion of star clusters
versus their mean surface density within the half light radius.
The blue triangles, small green circles, magenta upside-down
triangles and the large red
circles represent old globular clusters in the Galaxy, M 31, NGC
5128, and young massive star clusters in a variety of galaxies,
respectively. Note that the masses have been estimated assuming
virial equilibrium.
However, many of the young clusters have surface
densities greater than
![]() ![]() |
Open with DEXTER |
Following the analysis by Maraston et al. (2004) we attempt to place
W30 and G114 in the broader context of gravitationally bound stellar
systems. For this we exploit the re-definition of the fundamental
plane known as -space (Bender et al. 1992) which combines the
three fundamental observable parameters (radius, surface brightness,
and velocity dispersion) into physically motivated values.
Figure 7 shows the position of W3 (asterisk), W30
(upward triangle), and G114 (downward triangle) in the
plane. In this space
traces the mass
of the system, while
measures the compactness of a system
for a given mass (the product of
/L and surface brightness). For
comparison we also show the mean positions of bulges and
ellipticals (B+E), dwarf ellipticals (dE), M 32 (all taken from Burstein et al. 1997, and assuming H0=75 km s-1 Mpc-1). Next, we add the
point-like objects in Fornax (the ultra-compact dwarf galaxies or UCDs) (the
average value of the four objects presented in Drinkwater et al. 2003). We also plot the positions of old globular
clusters in the Milky Way (small blue triangles) (McLaughlin & van der Marel 2005), M 31 (small green dots) (a collation of data from McLaughlin &
van der Marel in prep.), and NGC 5128
(magenta up-side triangles) (Martini & Ho 2004), for which we assumed a
constant (B-V) colour for old metal poor GCs, namely 0.7 mag.
We also add the young clusters taken from the literature (see
Table 6), which are shown as red points. The
large filled green squares represent the Nuclear Star Clusters (NCs)
in bulge-less disk galaxies (Böker et al. 2004; Walcher et al. 2005, see Table 6 for details). Finally, we
add the dwarf galaxy transition objects
(DGTOs) found in the Virgo galaxy cluster (Hasegan et al. 2005) as
filled blue squares. For the DGTOs we assumed that (B-V)=1, typical
of a 10 Gyr solar metallicity SSP.
The arrows which begin at W3, W30, and G114 represent the evolution
of the cluster in this space when the clusters are "aged'' to a common
age of 10 Gyr using the Maraston (2005) SSP models. Note,
however, that the SSP model tracks do not take mass loss (and hence
fading) by evaporation or due to external perturbations into
account.
![]() |
Figure 7:
![]() ![]() ![]() ![]() |
Open with DEXTER |
As was found for the massive cluster NGC 7252:W3 by Maraston et al. (2004), we find that NGC 7252:W30 and NGC 1316:G114 will evolve
into the region of -space occupied by the UCDs and DGTOs.
This shows a strong similarity between the most massive star clusters and
these enigmatic objects, and may suggest that they formed through
similar mechanisms.
Additionally, we can estimate the amount of mass loss which is
expected to occur within W30 and G114. From the SSP models of
Maraston (2005) we see that a cluster (assuming solar metallicity and
a Kroupa stellar IMF) is expected to lose 18% of its mass
between the ages of 400 Myr and 10 Gyr (i.e. between the present age
of W30 and the age of globular clusters). G114, with an age of
3 Gyr is only expected to lose
5% of its current mass to
stellar evolution. Using the analytic expressions for mass loss in a
tidal field of Lamers et al. (2005, Eq. (6)), we note that due to the strong
dependence on cluster mass, neither of the clusters studied here are
expected to lose a significant amount of mass due to disruption. W30
is expected to lose
8% of its mass over the next 10 Gyr while G114 is expected to lose just
5% of its mass over the next 7 Gyr. For this calculation we have assumed t4 (the
average time for a
cluster to disrupt) to be 1 Gyr,
based on the galactic GC population (Boutloukos & Lamers 2003).
However, we note that the conclusions reached are not significantly
affected by the choice of t4. As these small changes
would barely be visible in Fig. 7 and would only add
confusion, we choose not to show them.
Figure 8 again shows the
projection of
-space, except with all of the YMCs (red points)
aged to 10 Gyr using SSP models, again assuming only passive stellar
evolution of the cluster. Here we see that many of the YMCs have
evolved "past'' the globular cluster region and into the space occupied
by W3, W30, the UCDs and the UGTOs. Burstein et al. (1997) suggest
that the tightness of the GC relation in
-space may be due to
the preferential destruction of star clusters outside a rather narrow
region of parameter space (e.g. mass and radius, see also Fall & Rees 1977). Along the same lines,
Gnedin & Ostriker (1997) show that only a narrow region of the
mass-radius plane of GCs is stable against disruption, and suggest
that the initial parameter distribution may have been much larger than
what is observed today.
This may be what we are seeing in Fig. 8 where the
young clusters occupy a much larger region of
-space (in terms
of mass, radius, and compactness) than their older globular cluster
counterparts.
We note that many young star clusters are not expected to
survive for more than 100 Myr (e.g. Bastian et al. 2005a; Fall et al. 2005) due to internal and external disruption mechanisms. Thus,
it can be expected that many of the youngest YMCs in our sample will
never survive to comparable ages as the Galactic GCs. We do not wish
to imply in Fig. 8 that all the YMCs will survive
to comparable ages, only that once differences in their stellar
populations are taken into account, the YMCs occupy a much larger
parameter space than their old GC counterparts. In particular we note
that YMCs tend to display extended envelopes in contrast to the
tidally truncated older GCs (e.g. Schweizer 2004). The loss of such
extended envelopes (i.e. due to
interaction with their environment) may significantly change the
position of the young clusters in Fig. 8. As
NGC 1316:G114 is
3 Gyr old it is likely that it has already lost
its extended envelope, which may explain why it falls on the
relation for old GCs in Fig. 8.
The similarities between young massive clusters and old globular
clusters have been shown in a number of recent works
(e.g. Kissler-Patig 2004; de Grijs et al. 2005). In particular,
Kissler-Patig (2004) showed that YMCs will follow a very similar
relation (which is one
projection of
-space) as old GCs once the fading of the
young clusters (due to stellar evolution) is taken into account.
However, the dangers of using such a projection
can be seen when comparing Fig. 1 in Kissler-Patig (2004) with
Figs. 7 and 8 in the present
work. In the
projection, the UCDs appear to be quite
consistent with the distribution of young and old clusters. However,
in
-space (which also includes information on the size) we see
that the UCDs are disjoint from the globular clusters. Additionally,
as stated above, in
-space the young clusters clearly occupy a
much larger parameter space than the old GCs, contrary to what is seen
in the
projection. Finally, we note that the
UCDs, DGTOs, and YMCs with masses above a few
all appear
very similar with respect to their scaling relations (i.e. mass,
velocity dispersion, size, and mean mass density), suggesting a common
formation mechanism (Kissler-Patig et al. 2005).
![]() |
Figure 8:
![]() ![]() |
Open with DEXTER |
We have presented velocity dispersion, effective radius, and hence dynamical mass measurements of two extremely massive clusters in galactic merger remnants. These results confirm that galactic mergers can produce star clusters with masses well in excess of the most massive globular clusters in the Milky Way. However, we have also shown that while these clusters are extremely massive, they are consistent with being the high-mass end of a continuous power-law distribution of star clusters, suggesting that cluster formation is (mass) scale independent.
Comparing the light-to-mass ratios of W30 (in NGC 7252) and G114 (in NGC 1316) to those predicted by simple stellar population models (at the ages of the clusters), shows that both of these clusters are consistent with having Kroupa-type stellar mass functions. Applying the same analysis to other young clusters taken from the literature shows that there is a significant age dependence on how well SSP models fit the light-to-mass ratios of young clusters. Therefore, it is possible that the deviation from the light-to-mass ratio of young clusters from that predicted by SSP models is not due to a varying stellar mass function, but instead reflects the state of equilibrium of the youngest clusters.
We have shown that W30 and G114 currently reside at the high-mass tip
of the old globular cluster distribution in -space (a
re-definition of the fundamental plane). Both clusters, along with
many young clusters taken from the literature, are likely to
evolve into the space occupied by the so-called ultra-compact dwarf
galaxies (UCDs) and the dwarf galaxy transition objects (DGTOs). This
shows that young massive clusters and UCDs/DGTOs share many similar
properties and suggests that the enigmatic UCDs/DGTOs may have formed
in a similar manner as the most massive globular clusters in mergers,
i.e., under rather violent circumstances.
Additionally, we showed that young massive star clusters will occupy a
much larger region of -space than presently occupied by old
globular clusters. This is consistent with the interpretation that
star clusters are born with a larger range of parameters (e.g. radius,
mass, and compactness) than displayed by globular clusters, and
destructive processes whittle away at the initial full distribution
with only clusters which have parameters within a small range surviving
to old ages.
Acknowledgements
We thank Marcelo Mora for his help in generating the ACS WFC PSF. Dean McLaughlin is gratefully acknowledged for providing a uniform electronic table of the parameters of the globular clusters in the Milky Way, M 31, and NGC 5128. We also thank Mark Gieles for insightful discussions, as well as Linda Smith for her help in understanding the properties of cluster M82F. F.S. gratefully acknowledges partial support from the National Science Foundation through grant AST-0205994.