A&A 455, L17-L20 (2006)
DOI: 10.1051/0004-6361:20065567
LETTER TO THE EDITOR
H. J. G. L. M. Lamers1,2 - M. Gieles1
1 - Astronomical Institute, Utrecht University, Princetonplein 5, 3584CC Utrecht, The Netherlands
2 - SRON Laboratory for Space Research, Sorbonnelaan 2, 3584CC, Utrecht, The Netherlands
Received 8 May 2006 / Accepted 26 June 2006
Abstract
We predict the survival time of initially bound star clusters in the solar
neighbourhood taking into account: (1) stellar evolution, (2) tidal stripping, (3) shocking by spiral arms and (4) encounters with giant molecular clouds.
We find that the predicted dissolution time is
Gyr for clusters in the mass range of
.
The resulting predicted shape of the logarithmic age distribution agrees very well with the empirical one, derived from
a complete sample of clusters in the solar neighbourhood within 600 pc.
The required scaling factor implies a star formation rate of 4
102
Myr-1 within 600 pc from the Sun or a surface formation rate of 3.5
10-10
yr-1 pc-2 for stars in bound clusters
with an initial mass in the range of 102 to 3
104
.
Key words: Galaxy: open clusters and associations: general - Galaxy: solar neighbourhood - Galaxy: disk - Galaxy: general - Galaxies: star clusters - Galaxy: kinematics and dynamics
The first empirical determination of the lifetime of clusters in the
solar neighbourhood is by Oort (1958), who noticed the
lack of clusters older than a few Gyr in the solar neighbourhood.
Later, Wielen (1971) derived a mean dissolution time of 0.2 Gyr from the age distribution of clusters. Since most of the observed
clusters within about 1 kpc from the sun have a mass in the range of 102 to a few 103
the value derived by Wielen is for clusters in that mass range. Theory
predicts that the dissolution time of clusters depends on their
initial mass in that massive clusters survive longer than low mass
clusters (e.g. Spitzer 1958;
Wielen 1985; Chernoff & Weinberg 1990;
Gnedin & Ostriker 1997, and references
therein). Baumgardt & Makino (2003) (hereafter BM03) showed from
N-body simulations that the dissolution time of clusters in the tidal
field of the galaxy depends on their initial mass, Mi, as
.
Independently, this same power-law dependence was
also derived empirically in a study of cluster samples in four galaxies by Boutloukos & Lamers (2003).
The dissolution time of clusters in the solar neighbourhood was
recently redetermined by Lamers et al. (2005) (hereafter L05),
based on a new cluster sample of Kharchenko et al. (2005).
They found a dissolution time of
with
0.5 Gyr
for clusters with
102 < M < 104.
This is a factor 5 shorter than the
Gyr
that follows from the
relation derived
by the N-body simulations of BM03 for clusters more massive than 4500
at a Galactocentric distance of 8.5 kpc. The simulations of BM03 include a realistic stellar mass function, stellar evolution, two-body relaxation, a detailed treatment of binary evolution and close encounters of stars. Part of this difference may be due to the fact that in low mass clusters,
with lifetimes shorter than about 1 Gyr, the dynamical evolution is affected
by the presence of massive stars during most of their lifetime (see Fig. 5
of BM03). It is doubtfull that this effect alone can fully explain the
difference between the results of L05 and BM03. In fact, the large discrepancy suggests
that other, probably external, disruptive effects must play an important role in
destroying star clusters in the solar neighbourhood.
In this paper we explain the lifetime of clusters in the solar neighbourhood, by taking into account the combined effects of stellar evolution, tidal stripping, encounters with giant molecular clouds (GMCs) and spiral arm shocks. We use stellar population models to describe the stellar evolution and the results of BM03 for tidal stripping. For the effects of GMCs and spiral arms we adopt the new estimates from the recent studies by Gieles et al. (2006c) (hereafter GPZB06) and Gieles et al. (2006a) (hereafter GAPZ06), which are based on N-body simulations.
The structure of the paper is as follows. In Sect. 2 we discuss the predicted mass loss from star clusters by stellar evolution, tidal stripping, encounters with GMCs and spiral arm shocks. We calculate the mass evolution of clusters due to these four effects. In Sect. 3 we compare our results with the observed age distribution of clusters in the solar neighbourhood. The discussion and conclusions are given in Sect. 4.
The mass loss from clusters due to stellar evolution has been
calculated for cluster evolution models by several groups. We adopt
the GALEV models for single stellar populations with a Salpeter type
mass function in the range of
(Anders & Fritze-v. Alvensleben 2003; Schulz et al. 2002). These models are
based on stellar evolution tracks from the Padova group, which include
overshooting, mass loss due to stellar winds and supernovae.
(Bertelli et al. 1994; Girardi et al. 2000). L05
have shown that the fraction of the initial cluster mass that is lost
by stellar evolution,
,
can be approximated
accurately by
BM03 have calculated a grid of N-body simulations of clusters in
circular and elliptical orbits in the tidal field of
a galaxy for different initial cluster masses, galactocentric distances R, and
different cluster density profiles. The stars follow a Kroupa initial mass function
and stellar evolution is taken into account during the evolution.
Gieles et al. (2004) have shown that for all models with
clusters with
of BM03 the dissolution time can be expressed as a function of the initial cluster mass as
GAPZ06 studied the dissolution of star clusters by spiral arms by
means of N-body simulations. They used and adjusted the analytical
expression of Ostriker et al. (1972) for the dissolution time of
star clusters due to disk shock, to derive an expression for the
dissolution time of star clusters by spiral arms (
).
Mass loss by spiral arm shocks will occur just at the moment the cluster
crosses the spiral arm. Assuming that spiral arms move with a constant pattern speed (
)
and that the matter in the disk has a constant circular
velocity (
), the relative velocity between the two (
)
depends on the location in the galaxy (R). Density waves that pass with a low velocity have a large effect on the star clusters (e.g. Ostriker et al. 1972). Therefore, the disruptive
effect of spiral arm shocks is most important close to the corotation
radius (
), i.e. the point where the disk and the spiral
arms have the same rotational velocity. We adopt the "average''
spiral arm model of GAPZ06, which is based on the study
of Elmegreen et al. (1989) of the spiral galaxies M81 and M100, to
derive the density contrast of the spiral arm.
Dias & Lépine (2005) found
for the spiral arms in the Galaxy, from a study of the nearby star
clusters, and a corotation radius (
)
almost coinciding with the
solar radius
0.08. Based on the adopted values of R0 =
8.5 kpc,
and the assumption that our Galaxy has 4 spiral arms (Vallée 2005), GAPZ06 used
and
Gyr.
Taking into account the ratio
between
the energy gain and the mass loss, predicted by GAPZ06 we find for
the solar neighbourhood that
GPZB06 studied the encounters between GMCs and clusters with N-body
simulations. They derived an expression for the energy gain and the
resulting mass loss for the full range of encounter distances, from
head-on to distant encounters. Adopting
a mean GMC density in the galactic plane near the sun of
pc-3, a surface density of GMCs
pc-2 (Solomon et al. 1987) and a mean
velocity dispersion of clusters and GMCs of
,
they derived a dissolution time (
)
for clusters by GMC encounters in the solar neighbourhood of
The decrease of mass due to the combined effects of stellar evolution,
tidal stripping, spiral arm shocks and GMC encounters can then be described as
![]() |
Figure 1:
The mass evolution of a cluster with an initial mass of
![]() |
Open with DEXTER |
Figure 2 shows the ages of clusters when their remaining
mass is 0 and 100
as a function of the initial mass.
The almost linear part from log (Mi/
) = 3.5 to 5 has a slope of about 0.67.
The figure also shows the dissolution times by the Galactic
tidal field, predicted by BM03 for clusters with an initial concentration factor W0=5in a circular orbit at R0=8.5 kpc. Our predicted timescales are about a factor 5 smaller, which agrees with the empirically determined dissolution time (L05).
![]() |
Figure 2:
The predicted dissolution times of clusters in the solar
neighbourhood due to the combined effects of stellar evolution,
tidal field, spiral arm
shocks and encounters with GMCs, as a function of the initial mass.
Full line: total dissolution time. Dashed line: time when the
remaining mass is 100 ![]() |
Open with DEXTER |
Given the initial mass distribution of the clusters, their formation
rate, CFR(t), and the time it takes for a dissolving cluster to fade
below the detection limit, we can predict
the distribution of observable clusters as a function of age
or mass. L05 have derived an expression for the general
case of a cluster sample that is set by a magnitude limit. Here we are
interested in the prediction for a cluster sample that is complete
down to a mass limit of 100 ,
because this is the mass limit of
the unbiased sample of clusters within 600 pc of Kharchenko et al. (2005) (see L05).
![]() |
Figure 3:
The observed age distribution of an unbiased sample
of clusters with
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
For a constant CFR and a power law cluster IMF with
a slope of
(Lada & Lada 2003) the number of clusters with
as a function of age is
Figure 3 shows a comparison between the observed age distribution of clusters
with
within 600 pc (from L05) with the predicted distribution for
.
This value of
is adopted because the
observed distribution shows a steep drop at
(with
only one cluster in the last bin) and Fig. 2 shows that this corresponds to
.
The predicted relation for twice higher or lower values of
agree worse with the observed relation. However, see discussion in Sect. 4. We have also calculated the expected age distribution in case there was no mass-radius relation for the clusters, i.e. for
.
The downward slope of the resulting distribution (not shown here) is significantly less steep than the one for
and does not fit the observed distribution.
The flattening of the predicted distribution at the low age end is due to the
fact that clusters with an initial mass in the range of
about 100 to 300
quickly reach 100
(see Fig. 2). The bump in the observed distribution around
is due to a local starburst (see L05 and
Piskunov et al. 2006). Notice the good agreement in the shapes of the predicted and observed
distributions!
The vertical shift that is applied to the predicted curve
to match the observed one gives a value of
C=10-4.15 in
Eq. (10), which corresponds to a SFR of 4
for bound clusters in the range of
104 within 600 pc from the sun.
We studied the dissolution of star clusters in the solar neighbourhood due to four effects: stellar evolution, tidal stripping, spiral arm shocks and encounters with GMCs. For this study we adopted the descriptions of GAPZ06 and GPZB06 for the dissolution of star clusters by spiral arms and encounters with GMCs. We found that the last effect plays a dominant role in the solar neighbourhood.
The cluster dissolution time due to spiral arms and GMCs depends
on the density of the clusters, i.e. on
.
This implies that
the dissolution time is
with
if the radius of a cluster depends on its mass
as
.
We adopted
as found by
Larsen (2004) for clusters in spiral galaxies and so
for dissolution by both spiral arms and GMCs. This
value is very similar to
predicted for dissolution by
the tidal field only (BM03) and empirically derived for cluster
samples in four galaxies by Boutloukos & Lamers (2003).
If there was no mass radius dependence for clusters in the solar neighbourhood,
e.g.
,
then the predicted age distribution would have
a shallower slope than shown in Fig. 3, since
more old (massive) clusters would have survived.
Our calculated dissolution times of clusters in the solar neighbourhood are about a factor five smaller than predicted by BM03 for clusters in the tidal field of our Galaxy, with stellar evolution, binaries and two-body relaxation taken into account. This is reminiscent of the short dissolution time of clusters in the central region of the interacting galaxy M51, where the empirical dissolution time is even ten times shorter than can be explained by stellar evolution and tidal fields (Gieles et al. 2005). GMCs severely limit the lifetime of clusters in that galaxy also (see the discussion in GPZB06).
The steep drop in the observed age distribution at
Gyr
can be explained by an upper mass limit for the initial cluster
mass in the solar neighbourhood of about 3
104
.
However, this
value is uncertain because it depends crucially on the completeness of
the used sample at ages above 1 Gyr. (The sample contains only six clusters older than 1 Gyr.) The mass versus age distribution of
our adopted sample, shown in Fig. 8 of L05, suggests that the lower
mass limit of the observed clusters increases steeply for clusters
older than 1 Gyr. Since the predicted value of
at any age
depends on
and
as given in
Eq. (10), and
increases when the lower mass limit
increases, an increase in this limit implies an increase in
derived from the observed age distribution. Based on this argument and the small number of clusters older than about 1 Gyr
in the observed sample, we conclude that the derived value of
should be considered as a lower limit of
the real maximum initial mass.
The vertical shift applied to the predicted age distributions to
match the observed one indicates a star formation rate of 4
102
Myr-1 in bound clusters of
within a distance of 600 pc, corresponding to a surface formation rate of 3.5
10-10
yr-1 pc-2. This is a factor 2 to 3 smaller than the SFR
derived from the study of embedded stars by Lada & Lada (2003)
because many of the stars are born in unbound clusters that dissolve
within 10 Myr.
The very good agreement between the predicted and observed age distribution of clusters shows that dissolution of clusters in the solar neighbourhood is dominated by encounters with GMCs, as was already suggested by Oort (1958). In fact, the good agreement may be slightly fortuitous because we have underestimated the dissolution by two-body relaxation (see Sect. 1) and slightly overestimated the dissolution by encounters with GMCs, because we adopted the midplane density of GMCs whereas clusters may spend a faction of their lifetime above or under the galactic disk. Both effects are expected to be smaller than a factor two and may partially cancel out.
Acknowledgements
This work is supported by a grant from the Netherlands Research School for Astronomy (NOVA).