International Journal of Astronomy and Astrophysics, 2012, 2, 167-173
http://dx.doi.org/10.4236/ijaa.2012.23020 Published Online September 2012 (http://www.SciRP.org/journal/ijaa)
Examining the M67 Classification as an Open Cluster
Shim’on Naim1, Evgeny Griv2
1Ilan Ramon Physics Center, Ben-Gurion University of the Negev, Beer-Sheva, Israel
2Department of Physics, Ben-Gurion University of the Negev, Beer-Sheva, Israel
Email: shimonnaim@yahoo.com, griv@bgu.ac.il
Received May 22, 2012; revised June 28, 2012; accepted July 15, 2012
ABSTRACT
The cluster M67 (= NGC 2682) in Cancer is a rich stellar cluster, usually classified as an open cluster. Using our own
observations with the 0.4 m telescope, we show that M67 is a tight group of about 1200 stars. The actual radius of the
cluster is about 3.1 pc and the average mass of a star in the system is about 1.33
M
. We also show that the ratio of the
mean kinetic energy of the cluster to its mean gravitational potential energy clust clust0.3KU , while the value pre-
dicted by the virial theorem is equal to . So the system is a gravitationally bound. This value of
0.5clust clust
KU is
considered as an evidence of quasi-stability of the cluster and allows us to use the Chandrasekhar-Spitzer relaxation
time for M67 chand 30
Myr as a characteristic dynamical relaxation time of the system under study. As the cluster is
almost twice older its half-life time hal chand
tf70
0
, it is argued that M67 was in the past (about 4 Gyr ago, close to its
formation) a relatively small ( stars) globular cluster, but got “open cluster” shape due to the dynamical
evaporation of the majority of its stars.
~ 10,00
Keywords: Photometers; Techniques; Photometric; Open Clusters and Associations; M67; Galaxies; Star Clusters
1. Introduction
Usually classified as an open cluster, the dense stellar
cluster M67 in Cancer is one of the oldest known open
clusters in our Galaxy with an age clust of
Gyr, that is, close to the age of the sun [1-3]. In the
present paper we adopt the value clust Gyr [4,5]. It
is important to note that open clusters usually get de-
structed much faster, over yr timescales. It is
some pc distant the Sun. Following Ref. [6],
we adopt the distance clust pc. The total number
of stars in M67 is estimated at over 700 [7]. M67 is the
nearest old open cluster. The cluster is also anomalous
with respect to the distance of clust pc to the
Galactic equatorial plane, where most open clu-
sters are found. This distance clust is even larger than
scale height pc of the local thick disk of
the Galaxy (see, e.g. Ref. [8]). Other open clusters lie
much closer to the Galactic plane. M67 current mass is
clust
t
=4
60
3.5 -4.5
t
8
z
~10
= 770
z
750- 850
=
h
1400
r
=0
400
0
z
-300
M
M and its initial mass was to be appro-
ximately
p
astclust [9]. Mass segregation has
been suggested as the cause for the observed structural
properties of M67 in Refs. [10,11]. This process has long
been suspected for globular star clusters (and has been
directly seen in globular cluster 47 Tucanae [12]), but
has never been suggested for open clusters. We conclude
that M67 has properties of globular clusters.
=10MM
Here, we study M67 with the use of the Meade LX-
200R robotic 0.4 m telescope of the Ilan Ramon Physics
Center Observatory in Israel (http://ilanramon.bgu.ac.il)
(Figure 1). The telescope operates from a small dome
and is equipped with Meade DSI III pro 1360 1024
CCD camera (Figure 2). In the last four years this
telescope was intensively used in a variety of monitoring
projects, including the study of stellar clusters, in par-
ticular, M67 and M35 (=NGC 2168).
2. Observations
The cluster M67 was photographed in two filters, V and
B, on 7th April 2010. An (f/6.3) local reducer was used,
and therefore the scale is 0.52 arcsec/pixel. Two sets of
exposure times were used in the observations, namely,
long (10 s in both B and V filters) and short (3 s in both B
and V filters).
Figure 3 shows a part of M67. This photograph is a
combination of 13 frames from each exposure, 10 s and 3
s, summed after a registration. The picture was taken
with a blue Johnson filter. First, we estimated the angular
size of the cluster = 27.63.6
= 3.1
R
6

arcsec. Then, adopt-
ing the revised distance to the cluster pc, its
actual radius was found clust pc. In order
to cast more light on the structure of M67, in Figure 4
we show a wider field 81 photograph of M67,
clust = 770r
0.26
C
opyright © 2012 SciRes. IJAA
S. NAIM, E. GRIV
168
Figure 1. The Meade LX200R robotic 16" telescope of the
Ilan Ramon Physics Center Observatory in Beer-Sheva
(Israel).
Figure 2. The dome of the Meade LX200R robotic 16" tele-
scope.
Figure 3. This photograph shows 173 stars of M67 cluster in
our observations.
Figure 4. A wider field 8' × 16' photograph of M67, taken
from the Atlas Image website.
taken from the Atlas Image website. Using these photo-
graphs (Figures 3 and 4) and this Rclust, finally we es-
timated the total number of stars in the system
clust 1200N
.
Figure 5 shows the observed color-magnitude,
vis . -mvsBV, diagram (or the HR diagram) for M67,
based on our observations. Each dot denotes a star. The
purple color indicates the positions of the main sequence
stars. The red color indicates the positions of the hori-
zontal branch stars. The yellow color indicates the posi-
tions of the giant stars. The positions of blue stragglers
are indicated in the diagram by the blue color. As is seen,
the color-magnitude diagram for M67 shows a short
main sequence extending from , where the
turn-off point is located, up to . After

-0BV

-BV
.4
.70
-0.7BV it starts merging in the field region stars
and is not seen clearly. The giant branch is also clearly
seen. Figure 5 evidently illustrates that M67 is quite old
cluster. There are stars that are brighter and bluer than
the main sequence stars and giant branch stars. These
stars are a part of an exotic class of blue stragglers. A
spectroscopic study of the blue stragglers in M67 has
been done in Ref. [13]. Notice that the HR diagram
shown in Figure 5 is more typical for globular clusters
than for open clusters.
Figure 6 presents our photometric color-magnitude
diagram for the “ordinary” open cluster M35. (The clu-
ster age is about 130 Myr [14].) In sharp contrast to M67,
the color-magnitude diagram for M35 clearly shows a
well-defined long main sequence stretching from
-BV 0.1
to
-1.BV4, but does not show
horizontal branch stars, the giant branch stars, and blue
stragglers. Such an HR diagram is typical for open clu-
sters.
Thus, the following parameters were derived from our
observations: =27.63.6
clust = 770r arcsec, Rclust = 3.1 ± 0.26
pc (adopting pc), and .
clust 1200N
Copyright © 2012 SciRes. IJAA
S. NAIM, E. GRIV 169
Figure 5. The color-magnitude diagram for M67, showing
the apparent magnitude mvis of stars versus their color in-
dex (B - V).
Figure 6. The color-magnitude diagram for M35 (cf. Ref.
[14]).
3. Total Energy of M67
One of the principal features of stellar clusters is their
quasi-stability: Hundreds of millions of years may elapse
before the cluster will begin to disintegrate by some stars
acquiring velocities exceeding the escape velocity esc .
This is because the total energy of a cluster in a quasi-
steady state is negative. Let us show that the assumption
of a quasi-steady state for M67 is reasonable.
v
The average mass of a main-sequence star in the
cluster star was evaluated by using the standard mass-
luminosity relation . From the definition
abs , one gets the following mass-
absolute-magnitude (
M
=2.5l
m
3.5
LM
onst
L
abs
og c
M
m) relation:


8.75
abs, abs
abs =10 .
(1)
mm
Mm M
From the HR diagram (Figure 5) and Equation (1),
one can estimate the average mass of a star in the cluster
clust
star abs,
=1
clust
1
=.
N
i
i
MNM
N
i
m
(2)
As is seen, however, this diagram is too noisy for the
dimmest stars (and some of the dimmest stars are not
shown in Figure 5 due to a high noise). Therefore, we
prefer to use a much more simple model. Namely,
examining our and other (e.g. Ref. [1]) HR diagrams, one
can see that the main sequence of this cluster is occupied
mainly by stars with absolute magnitude ranging from
to abs , with more stars closer to
abs . Also, from Figure 5 we obtained that the
number of stars
abs =3m
=3m=6m

=3 2=6Nm
=3m
abs
Nm abs . Therefore,
in the lowest approximation the cluster consists only of
stars with absolute magnitude abs . In the next
approximation, the average mass of a star is obviously
 
abs abs
star
2=3 =
.
3
Mm Mm
M
6
Taking abs, 4.85m
, we get star 1.33
M
M. Fi-
nally, using g, we have
33
10= 1.989M
star 2.650.910 g.M
33 (3)
This average mass of a star is in good agreement with
most recent studies of M67.
In order to estimate the kinetic energy, an average
speed of a star was calculated as
22
star =
r
vvv.
(4)
The typical radial velocity is
s –1 and the transverse velocity is

3
=0.48 0.15 10
r
v
3
=0.81 0.1 10v
 m·s–1 [7]. This value of v
is
based on the distance clust pc [7]. We corrected = 870r
v
to a new one based on the distance clust = 770100r
pc by taking
70= 0.8v
17708 km·s–1. From Equ-
ation (4) one easily obtains
3
star=0.860.210 ms.v
 
1
(5)
The kinetic energy of the cluster would then be
clust
K

2
39
star star
clust clust
==1.10.610 J.
2
Mv
KN  (6)
The average gravitational potential energy clust was
calculated using the well-known relation for energy of a
self-attracting spherically symmetric mass M. The gra-
vitational potential energy of the mass M of the radius R
is given by
U
2
3.
5
GM
UR
 (7)
For the cluster we use , clust
and get
clustclust star
=MNM =RR
22
clustclust clust
UG . This imme-
diately implies that the potential energy of the cluster is
=35 MNR
Copyright © 2012 SciRes. IJAA
S. NAIM, E. GRIV
170
39
clust =4.3310 JU . (8)
Thus, the total energy of M67,
39
clustclust clust
=3.2EKU10 J,
is negative but small in its absolute value in a sense that
the kinetic energy and the potential energy are of the
same order of magnitude. Moreover, the value of
clust clust0.3KU , while the value predicted by the virial
theorem for a stable system consisting of N particles,
bound by potential forces, is clust clust
KU . From
the above, one concludes that the assumption of a
quasi-steady state for M67 is indeed reasonable. This is
an important fact from the point of view of stellar
dynamics. One understands, however, that even though
the cluster is dynamically bound, clust, escaping
of stars from the system due to the Maxwellian-like
velocity spread is possible (see Section 4 below for an
explanation).
=0.5
<0E
4. Evaporation of Stars from M67
Many factors cause stars to dissipate from dynamically
bound clusters. In Ref. [15], it has pointed out that
because of the scattering of stars in clusters on the stars
in the galactic field the total internal energy of a cluster
should rise, ultimately leading to a complete breakup of
the cluster. It turned out that the characteristic time for
this process of destruction of dense stellar clusters is at
most 1010 yr. Spitzer [16,17] has calculated the increase
of energy leading to desintegration of a cluster due to
encounters with interstellar clouds. Spitzer has argued
that the characteristic time for a process of this sort is
also 1010 yr. The importance of the galactic tidal force on
the evolution and stability of clusters has been also
discussed (e.g. Ref. [18]). The characteristic decay time
for this process is more than 109 yr. We conclude that the
time yr might be regarded as an upper limit
on the age of clusters in the Galaxy.
10
age =10t
Other investigations have studied the most important
mechanism for the destruction of dense stellar clusters,
namely, the dissipation which is caused by the interaction
of stars as they approach one another [19-22]. See also
Refs. [23-25] for a discussion of the problem. Evidently,
in this process energy exchange among cluster members
will cause individual stars to acquire supercritical velo-
cities (escape velocities) of the order of
2
esc 2vv, where 2
v is the rms velocity) in ex-
cess of the parabolic velocity and actually to leave the
cluster, whose total mass therefore decreases. The escap-
ing stars carry away a positive energy and a state of
statistical equilibrium is impossible for a cluster contain-
ing a finite number of stars if the interactions between
stars is strictly taken into account [26,27]. See also Ref.
[28] for a discussion.
By the time of relaxation of star systems, generally
speaking we mean the characteristic time of approach of
the distribution function
,,
f
trv of the stars with re-
spect to peculiar (random) velocities to a Maxwellian
distribution [22]. (The distribution function
v
,,
f
trv
ddvr
,d
is
defined by the condition that is the
number of stars at time t in volume element

,,ftrv
rr
M
r
and in velocity element .) The time of estab-
lishment of a Maxwellian distribution of random velo-
cities as a result of encounters between test stars test
and field stars field (the “collisional” relaxation time)
has the order of magnitude [22,23]
,dvvv
M
3
chand 22
field
=,
10 log
v
GM n
(9)
where is the relative velocity, which is approximately
the velocity dispersion of the lighter stars, n is the
number density of the field stars, and field are
the masses of test and field stars, is the so-called
Coulomb (Newton) logarithm, by means of which the long-
range nature of the gravitational force is taken into ac-
count, clust
v
test
M
log M
0.4N
, and field is the total number of
field stars in the system. If we assume that the stellar
distribution in the system tends to a Maxwell-Boltzmann
distribution as a result of stellar encounters, during re-
laxation time the system loses approximately 0.0074 of
its members (e.g. Ref. [28]). Thus we can properly regard
half chand
N
70t
as a measure of the half-life time of a
cluster. (Although the general conclusions on the relax-
ation time reached by Chandrasekhar, Spitzer, and others
are correct, the discussions were based on a simplified
“molecular-kinetic” theory and not on an explicit solu-
tion of the time-dependent problem. In Appendix 8, we
follow numerically the relaxation of an isotropic system
of like stars through small-angle Coulomb encounters by
using the more accurate approach.)
In M67, assuming fieldteststar =1.33
M
MM M
=0.86
(Equation (3)), star
vv
km·s–1 (Equation (5)),
and 3
4πR
9
~4 10
~10,0
clustclust , one obtains yr
for star-star gravitational encounters, So that the half-life
time of the cluster is relatively short, of the order of
yr only. (We estimate that the dynamical
Chandrasekhar-Spitzer relaxation time for M67 is
yr, which implies that the cluster age is 133
times its relaxation time. In contrast for the “ordinary”
open cluster M35, the cluster age yr is only
1.6 times its relaxation time [14].) We argue therefore
that in the past ( yr ago) the cluster M67 was a
relatively small ( stars) globular cluster, but
got “open cluster” shape due to the dynamical eva-
poration of a substantial portion or even the majority of
its stars. Interestingly, the direct N-body model with
36,000 stars for M67, evolved from zero age to 4 Gyr,
=3nN 7
chand 310

8
1.3 10t
00
9
half =210t
7
310
Copyright © 2012 SciRes. IJAA
S. NAIM, E. GRIV 171
has been presented in Ref. [9]. It has been shown that at
4 Gyr the total mass has reduced to 2000
M
1200N
as a
result of mass loss and stellar escapes.
5. Concluding Remark
Summarizing, from our own observations, we built the
HR diagram for the cluster M67. We showed that this
diagram is more typical for relatively old globular
clusters than for open clusters. We also estimated the
total number of stars in the cluster clust , its
actual radius pc, and the average mass of a
star star
clust 3.1R
1.33
M
M
. We studied the dynamical features
of M67 and argued that
The total energy of the cluster clust clustust
U
cl
=EK
is negative, 39
clust . Moreover
3.210 JE
clust clust0.3KU. The latter may be considered as
evidence of the quasi-stability of the system under
consideration.
The cluster is a tight group of about 1200 stars which
are gravitationally bound, <0E.
clust
As the cluster is quasi-stable, its characteristic rela-
xation time is the Chandrasekhar-Spitzer time
7 yr and the half-life time is equal to
9
2 10yr, respectively.
chand310

70t
half chand
As the cluster is almost twice older its half-life time,
we conclude that many stars have left the cluster
since its formation. In the past (9
~4 10 yr ago) the
cluster was a relatively small (~ 10,000 stars) glo-
bular cluster, but got “open cluster” shape due to the
dynamical evaporation of the majority of its stars.
6. Acknowledgements
This work was began while the first author (S.N.) was an
undergraduate student at the Ben-Gurion University. S.N.
thanks the Head of Astronomical Programs of the Ilan
Ramon Physics Center, Netzach Farbiash, and the Scien-
tific Director of the Ilan Ramon Physics Center, Moshe
Schechter, for making possible observations with Meade
LX200R telescope. We thank Amir Bernat, Michael Ge-
dalin, and Nimrod Nissim for useful discussions in rela-
tion to this work and both anonymous referees for nu-
merous comments that improved the presentation of the
paper. Atlas Image obtained as part of the Two Micron
All Sky Survey, a joint project of the University of Mas-
sachusetts and the Infrared Processing and Analysis
Center/California Institute of Technology, funded by the
National Aeronautics and Space Administration and the
National Science Foundation. This work was sponsored
in part by the Israel Science Foundation, the Israeli Mi-
nistry of Immigrant Absorption in the framework of the
program “KAMEA,” and the Binational US-Israel Sci-
ence Foundation.
REFERENCES
[1] G. Meynet, J.-C. Mermilliod and A. Maeder, “New Dat-
ing of Galactic Open Clusters,” Astronomy and Astrophy-
sics Supplement Series, Vol. 98, No. 3, 1993, pp. 477-
504.
[2] S. Randich, F. Primas, L. Pasquini, P. Sestito and R. Pal-
lavicini, “Tracing Mixing in Stars: New Beryllium Ob-
servations of the Open Clusters NGC 2516, Hyades, and
M67,” Astronomy & Astrophysics, Vol. 469, No. 1, 2007,
pp. 163-172. doi:10.1051/0004-6361:20066218
[3] A. Sarajedini, A. Dotter and A. Kirkpatrick, “Deep 2MASS
Photometry of M67,” The Astrophysical Journal, Vol.
698, No. 2, 2009, pp. 1872-1878.
doi:10.1088/0004-637X/698/2/1872
[4] X. Fan, D. Burstein, J.-S. Chen, J. Zhu, et al., “Deep Wide-
Field Spectrophotometry of the Open Cluster M67,” The
Astrophysical Journal, Vol. 112, No. 2, 1996, pp. 628-
648. doi:10.1086/118039
[5] D. A. VandenBerg and P. B. Stetson, “On the Old Open
Clusters M67 and NGC 188,” Publications of the Astro-
nomical Society of Japan, Vol. 116, No. 825, 2004, pp.
997-1011. doi:10.1086/426340
[6] V. Laugalys, A. Kazlauskas, R. P. Boyle, F. J. Vrba, P. A.
Davis and V. Straižys, “CCD Photometry of the M67
Cluster in the Vilnius System. II. New Photometry of High
Accuracy,” Baltic Astron., Vol. 13, 2004, pp. 1-33.
[7] T. M. Girard, W. M. Grundy, C. E. Lopez and W. F. van
Altena, “Relative Proper Motions and the Stellar Velocity
Dispersion of the Open Cluster M67,” Astronomical
Journal, Vol. 98, No. 1, 1989, pp. 227-243.
doi:10.1086/115139
[8] E. K. Kharadze, R. A. Bartaya, O. B. Dluzhnevskaya, A.
E. Piskunov and E. D. Pavlovskaya, “Population of the
Galactic Disc in the Solar Neighbourhood. I—Parameters
of Spatial Distribution for Stellar Groups of A-K Spectral
and III-V Luminosity Classes,” Astrophysics and Space
Science, Vol. 151, No. 2, 1989, pp. 319-334.
doi:10.1007/BF00648388
[9] J. R. Hurley, O. R. Pols, S. J. Aarseth and C. A. Tout, “A
Complete N-body Model of the Old Open Cluster M67,”
Monthly Notices of the Royal Astronomical Society, Vol.
363, No. 1, 2005, pp. 293-314.
doi:10.1111/j.1365-2966.2005.09448.x
[10] C. Bonatto and E. Bica, “Mass Segregation in M67 with
2MASS,” Astronomy & Astrophysics, Vol. 405, No. 2,
2003, pp. 525-530. doi:10.1051/0004-6361:20030205
[11] J. R. A. Davenport and E. L. Sandquist, “Death of a
Cluster: The Destruction of M67 as Seen by the Sloan
Digital Sky Survey,” The Astrophysical Journal, Vol. 711,
No. 2, 2010, pp. 559-572.
doi:10.1088/0004-637X/711/2/559
[12] D. E. McLaughlin, J. Anderson, G. Meylan, K. Gebhardt,
C. Pryor, D. Minniti and S. Phinney, “Hubble Space
Telescope Proper Motions and Stellar Dynamics in the
Core of the Globular Cluster 47 Tucanae,” The Astro-
physical Journal Supplement Series, Vol. 166, No. 1,
2006, pp. 249-297. doi:10.1086/505692
[13] G. Q. Liu, L. Deng, M. Chavez, E. Bertone, A. H. Davo
Copyright © 2012 SciRes. IJAA
S. NAIM, E. GRIV
Copyright © 2012 SciRes. IJAA
172
and M. D. Mata-Chavez, “A Spectroscopic Study of the
Blue Stragglers in M67,” Monthly Notices of the Royal
Astronomical Society, Vol. 390, No. 2, 2008, pp. 665-
674.
[14] J. S. Kalirai, G. G. Fahlman, H. B. Richer and P. Ventura,
“The CFHT Open Star Cluster Survey. IV. Two Rich,
Young Open Star Clusters: NGC 2168 (M35) and NGC
2323 (M50),” The Astrophysical Journal, Vol. 126, No. 2,
2003, pp. 1402-1414. doi:10.1086/377320
[15] S. Rosseland, “Scattering of Stars in Clusters,” Zeitschrift
fur Astrophysik, Vol. 4, No. 1, 1932, pp. 235-254.
[16] L. Spitzer, “Distribution of Galactic Clusters,” The As-
trophysical Journal, Vol. 127, No. 1, 1958, pp. 17-27.
doi:10.1086/146435
[17] L. Spitzer and R. Harm, “Evaporation of Stars from Iso-
lated Clusters,” The Astrophysical Journal, Vol. 127, No.
1, 1958, pp. 544-550. doi:10.1086/146486
[18] E. Terlevich, “Evolution of N-body Open Clusters,” Month-
ly Notices of the Royal Astronomical Society, Vol. 224,
No. 1, 1987, pp. 193-225.
[19] V. A. Ambartsumyan, “Stability of Stellar Clusters,”
Uchenye Zapiski LGU, Vol. 22, No. 1, 1938, pp. 19-29.
[20] L. Spitzer, “The Stability of Isolated Clusters,” Monthly
Notices of the Royal Astronomical Society, Vol. 100, No.
1, 1940, pp. 396-413.
[21] T. A. Agekyan, “The Velocity Distribution Function and
the Rate of Dissipation in Systems of Gravitating Bod-
ies,” Soviet Astronomy, Vol. 3, No. 1, 1959, pp. 280-290.
[22] S. Chandrasekhar, “Principles of Stellar Dynamics,” Do-
ver, New York, 1960.
[23] L. Spitzer and M. H. Hart, “Random Gravitational En-
counters and the Evolution of Spherical Systems. I. Me-
thod,” The Astrophysical Journal, Vol. 164, No. 1, 1971,
pp. 399-409. doi:10.1086/150855
[24] R. de Grijs, “The Long-Term Survival Chances of Young
Massive Star Clusters,” Astrophysics and Space Science,
Vol. 324, No. 2-4, 2009, pp. 283-291.
doi:10.1007/s10509-009-0100-0
[25] J. Schneider, P. Amaro-Seoane and R. Spurzem, “Higher-
Order Moment Models of Dense Stellar Systems: Appli-
cations to the Modelling of the Stellar Velocity Distribu-
tion Function,” Monthly Notices of the Royal Astronomi-
cal Society, Vol. 410, No. 1, 2011, pp. 432-454.
doi:10.1111/j.1365-2966.2010.17454.x
[26] I. King, “The Escape of Stars from Clusters 11. A Simple
Theory of the Evolution of an Isolated Cluster,” The As-
trophysical Journal, Vol. 63, No. 1258, 1958, pp. 114-
117. doi:10.1086/107702
[27] S. von Hoerner, “Die Auflsungszeit Offener Sternhau-
fen,” Zeitschrift fur Astrophysik , Vol. 44, No. 1, 1958, pp.
221-242.
[28] L. P. Ossipkov, “On the Jubilee of Academician V. A.
Ambartsumyan Statistical Mechanics of Stellar Systems:
From Ambartsumyan Onward,” Astrophysics, Vol. 51,
No. 4, 2008, pp. 428-442.
doi:10.1007/s10511-008-9032-6
[29] A. F. Alexandrov, L. S. Bogdankevich and A. A. Rukhadze,
“Principles of Plasma Electrodynamics,” Springer, Berlin,
1984. doi:10.1007/978-3-642-69247-5
[30] G. G. Kuzmin, “The Effect of Star Encounters and the
Evolution of Star Clusters,” Publications of Tartu Ob-
servatory, Vol. 33, No. 2, 1957, pp. 75-102.
[31] M. N. Rosenbluth, W. M. MacDonald and D. L. Judd,
“Fokker-Planck Equation for an Inverse-Square Force,”
Physical Review, Vol. 107, No. 1, 1957, pp. 1-6.
doi:10.1103/PhysRev.107.1
[32] D. C. Montgomery and D. A. Tidman, “Plasma Kinetic
Theory,” McGraw-Hill, New York, 1964.
[33] B. A. Trubnikov, “Particle Interactions in a Fully Ionized
Plasma,” In: M. A. Leontovich, Ed., Reviews of Plasma
Physics, Consultants Bureau, New York, 1965, pp. 105-
140.
[34] E. Griv, “Numerical Integration of the Landau Kinetic
Equation,” In: G. Byrd, Eds., “Order and Chaos in Stellar
and Planetary Systems,” ASP, San Francisco, 2004, pp.
352-356.
[35] W. M. MacDonald, M. N. Rosenbluth and W. Chuck,
“Relaxation of a System of Particles with Coulomb In-
teractions,” Physical Review, Vol. 107, No. 2, 1957, pp.
350-353. doi:10.1103/PhysRev.107.350
[36] J. Killeen and A. H. Futch, “Numerical Solution of the
Fokker-Planck Equations for a Hydrogen Plasma Formed
by Neutral Injection,” Journal of Computational Physics,
Vol. 2, No. 3, 1968, pp. 236-254.
S. NAIM, E. GRIV 173
Appendix
Relaxation of a System of Stars:
Numerical Solution of the Fokker-Planck
Equation
Let us consider a system of stars interacting through in-
verse-square-law Coulombian (Newtonian) forces. As is
known, the evolution of an isotropic distribution of like
particles (stars) through small-angle Coulomb-like col-
lisions is described by the Fokker-Planck equation for an
inverse-square force:
2
22 4
23
00
2
00
2
21
=2πlnd d
3
4dd11
3
2
v
v
ff
Gmufu ufu
tvv
fu
ufuufu
vvv v
f

2
u


 

 
 


(10)
where

,
f
vt is the distribution function of stars,
which is defined by the condition that
2
0
=dnvf
v
is the number density of stars. In this approximation the
distribution function depends only on v, the magnitude of
the velocity, and t, the time. Equation (10) has been
derived by L. D. Landau (e.g. Ref. [29]) and then fully
independent in Ref. [30] for gravitational systems and in
Refs. [31,32] for an ionized gas. See Refs. [33,34] for a
discussion. It can be shown that, in the absence of sources
and losses, the number of particles (stars), as well as total
momentum and energy, are conserved within the for-
malism of the Fokker-Planck equation. The solution of
Equation (10), under the same conditions, for a steady
state (=0ft ) is the Maxwellian distribution,
2
2
() exp.
2
v
fv v





Following Refs. [35,36], Equation (10) may be put
into dimensionless form
2
2
=
fff
,
A
BCf
x
x
 


(11)
where 0
=
x
vv , v0 is a constant and is a characteristic
velocity,
is the dimensionless time, and quantities A,
B, and C are defined in Ref. [36]. We integrate this
equation directly using the difference equation given in
Ref. [36] (see also Ref. [35]).
Figure 7. The relaxation of a system of like stars to the
equilibrium Maxwellian distribution. The dashed line shows
an initial distribution and the dash-dot line shows a Max-
wellian of the same density and energy. The time is nor-
malized so that τ = 1 corresponds to a Chandrasekhar-
Spitzer relaxation time tchand (Equation (9)).
We agree with the second referee of this paper that
Equation (11) “is neither matching an actual multi-com-
ponent system nor a spatially inhomogeneous stellar
system.” This equation can serve only as a convenient
starting point for more realistic computations.
We solve Equation (11) numerically, subject to the
initial condition that at time =0
, the initial distri-
bution is
 
2
,0exp0.3 0.3fx x
. The result-
ing solution for
,
f
x
is plotted in Figure 7 for
various values of time
. The final Maxwellian distri-
bution function at
 is also shown. In fair agree-
ment with the simple Chandrasekhar-Spitzer theory [22,
23], at 1
the distribution function of stars becomes
Maxwellian-like. As calculations show, near 0.25x
and 0.45x
the distribution function departs slightly
from the equilibrium Maxwellian. However, at =1
the lower energy portion of the spectrum is somewhat
underpopulated, while the higher energy part is over-
populated. Thus, the higher- and lower-energy parts of
the distribution are filled in at a much later time (at times
10
, Figure 7) as one may expect from the Chan-
drasekhar-Spitzer theory. We conclude that in quanti-
tative terms the collisional relaxation of the high- and
low-energy parts of the distribution is only poorly des-
cribed by the Chandrasekhar-Spitzer simple “molecular-
kinetic” theory.
Copyright © 2012 SciRes. IJAA