The Threebody Problem Valtonen M Karttunen H PDF Download
The Threebody Problem Valtonen M Karttunen H PDF Download
download
https://ebookbell.com/product/the-threebody-problem-valtonen-m-
karttunen-h-2048332
https://ebookbell.com/product/the-threebody-problem-cixin-liu-56481884
https://ebookbell.com/product/the-threebody-problem-cixin-liu-56888538
https://ebookbell.com/product/the-threebody-problem-and-the-equations-
of-dynamics-poincars-foundational-work-on-dynamical-systems-
theory-1st-edition-henri-poincar-auth-5882734
https://ebookbell.com/product/the-threebody-problem-from-pythagoras-
to-hawking-1st-ed-mauri-valtonen-6745112
The Threebody Problem Omnibus Cixin Liu
https://ebookbell.com/product/the-threebody-problem-omnibus-cixin-
liu-35112578
https://ebookbell.com/product/the-threebody-problem-trilogy-
remembrance-of-earths-past-omnibus-cixin-liu-232240994
https://ebookbell.com/product/the-threebody-problem-remembrance-of-
earths-past-cixin-liu-35470160
https://ebookbell.com/product/the-threebody-problem-liu-cixin-11003028
https://ebookbell.com/product/the-dark-forest-the-threebody-problem-
liu-cixin-35996694
THE THREE-BODY PROBLEM
How do three celestial bodies move under their mutual gravitational attraction? This
problem was studied by Isaac Newton and other leading mathematicians over the
last two centuries. Poincaré’s conclusion, that the problem represents an example
of chaos in nature, opens the new possibility of using a statistical approach. For the
first time this book presents these methods in a systematic way, surveying statistical
as well as more traditional methods.
Cambridge University Press has no responsibility for the persistence or accuracy of s
for external or third-party internet websites referred to in this publication, and does not
guarantee that any content on such websites is, or will remain, accurate or appropriate.
Contents
Preface page ix
1 Astrophysics and the three-body problem 1
1.1 About the three-body problem 1
1.2 The three-body problem in astrophysics 5
1.3 Short period comets 8
1.4 Binary stars 12
1.5 Groups of galaxies 15
1.6 Binary black holes 17
2 Newtonian mechanics 20
2.1 Newton’s laws 20
2.2 Inertial coordinate system 21
2.3 Equations of motion for N bodies 22
2.4 Gravitational potential 24
2.5 Constants of motion 25
2.6 The virial theorem 27
2.7 The Lagrange and Jacobi forms of the equations of motion 29
2.8 Constants of motion in the three-body problem 31
2.9 Moment of inertia 32
2.10 Scaling of the three-body problem 34
2.11 Integration of orbits 34
2.12 Dimensions and units of the three-body problem 38
2.13 Chaos in the three-body problem 39
2.14 Rotating coordinate system 43
Problems 45
3 The two-body problem 47
3.1 Equations of motion 47
3.2 Centre of mass coordinate system 48
3.3 Integrals of the equation of motion 49
v
vi Contents
ix
x Preface
has generously provided research tools for the calculations in this book. Other
members of the Turku research group have also helped us with illustrations. Most
of all, we would like to thank Sirpa Reinikainen for typing much of the final text,
including the great many mathematical formulae.
Financial support for this project has been provided by Finland’s Society for
Sciences and Letters and the Academy of Finland (project ‘Calculation of Orbits’),
which gave the opportunity for one of us (M. V.) to concentrate on writing the book
for a period of two years. The generous support by the Department of Computer
Science, Mathematics and Physics (in Barbados) and the Department of Physics
(in Trinidad) of the University of the West Indies made it possible to carry out
the writing in optimal surroundings. Parts of the text were originally published in
Finnish: H. Karttunen, Johdatus Taivaanmekaniikkaan, Helsinki: Ursa, 2001.
Finally, M. V. would like to express his appreciation to Sverre Aarseth who taught
him how to calculate orbits (and much else), and to his wife Kathleen, the Caribbean
link, whose encouragement was vital for the accomplishment of the book.
1
Astrophysics and the three-body problem
1
2 Astrophysics and the three-body problem
solution (Lagrange 1778). The theory of the restricted three-body problem was
further developed by Jacobi (1836), and it was used for the purpose of identifying
comets by Tisserand (1889, 1896) and reached its peak in the later nineteenth
century with the work of Hill (1878) and Delaunay (1860). The ‘classical’ period
reached its final phase with Poincaré (1892–1899).
In spite of these successes in special cases, the solution of the general three-body
problem remained elusive even after two hundred years following the publication
of Principia. In the general three-body problem all three masses are non-zero and
their initial positions and velocities are not arranged in any particular way. The
difficulty of the general three-body problem derives from the fact that there are no
coordinate transformations which would simplify the problem greatly. This is in
contrast to the two-body problem where the solutions are found most easily in the
centre of mass coordinate system. The mutual force between the two bodies points
towards the centre of mass, a stationary point in this coordinate system. Thus the
solution is derived from the motion in the inverse square force field. Similarly, in
the restricted three-body problem one may transfer to a coordinate system which
rotates at the same rate around the centre of mass as the two primary bodies. Then
the problem is reduced to the study of motions in two stationary inverse square
force fields. In the general problem, the lines of mutual forces do not pass through
the centre of mass of the system. The motion of each body has to be considered
in conjunction with the motions of the other two bodies, which made the problem
rather intractable analytically before the age of powerful computers.
At the suggestion of leading scientists, the King of Sweden Oscar II established
a prize for the solution of the general three-body problem. The solution was to
be in the form of a series expansion which describes the positions of the three
bodies at all future moments of time following an arbitrary starting configuration.
Nobody was able to claim the prize for many years and finally it was awarded in
1889 to Poincaré who was thought to have made the most progress in the subject
even though he had not solved the specific problem. It took more than twenty years
before Sundman completed the given task (Sundman 1912). Unfortunately, the
extremely poor convergence of the series expansion discovered by Sundman makes
this method useless for the purpose of calculating the orbits of the three bodies.
Now that the orbits can be calculated quickly by computer, it is quite obvious why
this line of research could not lead to a real solution of the three-body problem: the
orbits are good examples of chaos in nature, and deterministic series expansions
are utterly unsuitable for their description. Poincaré was on the right track in this
regard and with the current knowledge was thus a most reasonable recipient of the
prize.
At about the same time, a new approach began which has been so successful in
recent years: the integration of orbits step by step. In orbit integration, each body,
1.1 About the three-body problem 3
P1 (1, 3)
m1 = 3
in turn, is moved forward by small steps. In the most basic scheme, the step is
calculated on the basis of the accelerations caused by the two other bodies during
that step while they are considered to remain fixed. There is an error involved when
only one body is moved at a time, and others move later, but this can be minimised
by taking short steps and by other less obvious means.
Burrau (1913) considered a well defined, but in no way special, initial configu-
ration of three bodies which has since become known as the Pythagorean problem
since the three bodies are initially at the corners of a Pythagorean right trian-
gle. The masses of the three bodies are 3, 4 and 5 units, and they are placed
at the corners which face the sides of the triangle of the corresponding length
(Fig. 1.1). In the beginning the bodies are at rest. Burrau’s calculation revealed
the typical behaviour of a three-body system: two bodies approach each other,
have a close encounter, and then recede again. Subsequently, other two-body en-
counters were calculated by Burrau until he came to the end of his calculating
capacity. Only after the introduction of modern computers and new orbit integra-
tion methods was the celestial dance in the Pythagorean problem followed to its
conclusion.
Later work has shown that the solution of the Pythagorean problem is quite typical
of initially bound three-body systems. After many close two-body approaches, a
configuration arises which leads to an escape of one body and the formation of
a binary by the other two bodies (see Fig. 1.2, Szebehely and Peters 1967). A
theoretical treatment of a three-body system of this kind is given in Chapter 7. In
the following chapter, situations are discussed where a third body comes from a
large distance, meets a binary, and perhaps takes the place of one of the binary
members which escapes. Such orbits were calculated already in 1920 (Becker
1920). Sometimes the third body is always well separated from the binary; then the
situation is best described by perturbations on the binary caused by the third body.
Some examples of these systems are discussed in Chapter 10.
4 Astrophysics and the three-body problem
6 6
4 4
2 2
0 0
-2 -2
-4 -4
-6 -6
-6 -4 -2 0 2 4 6 -6 -4 -2 0 2 4 6
t = 40 t = 60
6 6
4 4
2 2
0 0
-2 -2
-4 -4
-6 -6
-6 -4 -2 0 2 4 6 -6 -4 -2 0 2 4 6
t = 20 t = 70
6 6
4 4
2 2
0 0
-2 -2
-4 -4
-6 -6
-6 -4 -2 0 2 4 6 -6 -4 -2 0 2 4 6
t = 10 t = 80
Figure 1.2 Trajectories of the Pythagorean problem. The orbits up to time 10, 20,
40, 60, 70 and 80 units are shown. The last two panels are identical since the escape
has already happened and the bodies are outside the frame.
In recent years there has been increasing demand for solutions of the general
three-body problem in various astrophysical situations. For example, binaries and
their interactions with single stars play a major role in the evolution of star clusters
(Aarseth 1973). Triple stellar systems are another obvious astrophysical three-body
1.2 The three-body problem in astrophysics 5
problem. Many other astrophysical bodies, ranging from compact bodies to galax-
ies, occur in triple systems. These will be discussed in Chapter 11.
Before then, we have to do some preparatory work which is partly found in
standard textbooks on mechanics and celestial mechanics. The proper discussion
of the three-body problem is thus delayed until Chapter 5 where special cases are
first introduced. This leads to consideration of the general three-body problem in
the later chapters, the main objective of this book.
may want to skip much of the first chapters. Also the Hamiltonian techniques of
Chapter 4 are commonly studied in courses of advanced mechanics. They form
such an essential part of the treatment of weak perturbations of binaries that it is
necessary to introduce the Hamiltonian concepts also in this book.
A planet going around the Sun is an example of a binary. The third body could then
be another planet, a moon, an asteroid or a comet. Because of the large differences
between the masses of the bodies, from the dominant Sun down to asteroids of
kilometre size or even smaller, Solar System dynamics is special in many ways. A
very up to date treatment exists in this area (Murray and Dermott 1999). We discuss
only a small class of Solar System problems which are related to stellar dynamics
and therefore form a suitable introduction to later studies of three-body systems
with more equal masses. Thus Chapter 5 contains many topics which readers may
have encountered earlier.
At the present time, three-body astrophysics is primarily motivated by the need
to understand the role of binaries in the evolution of stellar systems. For most of the
time, a binary acts just like a single star in a stellar system. The distances between
stars are large compared with the sizes of the stars and even compared with the
sizes of close binary orbits. For a relatively brief moment a binary and a third star
interact strongly, a ‘new’ binary forms, and a ‘new’ third star leaves the scene.
The importance of the process lies in its ability to redistribute energy and angular
momentum efficiently; the population of binaries may become more and more
tightly bound as time and three-body scattering go on, while the population of single
stars may gain speed and become ‘heated’. This will have profound consequences on
the structure and evolution of a star cluster; for example single stars and sometimes
also binaries escape from the cluster; binary orbits may shrink to form contact
binaries, and also triple stars may form where the third body remains bound to
the binary. The end products of the three-body process may appear as sources of
radio jets, X-rays, gamma rays or as other kinds of ‘exotic’ objects (Hut et al.
2003).
These various scenarios can be reproduced by numerical orbit calculations. The
orbits of thousands of stars can be calculated in a simulation of a star cluster.
Even though these simulations have now reached a great level of complication
and trustworthiness (e.g. Heggie and Hut 2003, Aarseth 2003a), it is still useful to
examine the three-body process to see how much is understood from elementary
principles. Together with the simulations of large numbers of bodies one may attain
a deeper understanding of the evolutionary process.
Chapter 7 starts with the discussion of initially strongly interacting three-body
systems. We will learn that such systems have a limited lifetime, and we do not
expect to find very many of them in nature. But they are important in the descrip-
tion of the intermediate state between the impact of a third body on a binary and
1.2 The three-body problem in astrophysics 7
the departure of a ‘new’ third body from a ‘new’ binary. The techniques used in
Chapter 7, assuming a complete reshuffle of the positions and velocities of the three
bodies in the phase space, work surprisingly well.
By ‘work well’ we mean that a large body of numerical orbit calculations can
be described in a statistical sense by simple physical principles. That this should
be so is not at all obvious at the outset. Therefore we devote considerable space
to comparing numerical calculations with the theory. The beauty of it is that once
the theory is established for certain parts of the parameter space, we have good
reasons to expect that it applies more generally. Also we do not need to go back to
calculating millions of orbits when a slightly different astrophysical problem arises
but we can use the theory directly with a fair amount of confidence.
In addition to these practical considerations, it gives the reader a certain satisfac-
tion to learn that simple analytical solutions of the three-body problem exist, even
though only in a statistical sense and for a limited part of the parameter space. It
will also become clear that these are the only solutions of any significance in large
parts of the phase space due to the chaotic nature of the problem.
Chapters 9 and 10 try to cover the remaining parameter space, i.e. when the
binary and the third body are so well separated throughout the interaction that the
principle of complete chaos is not productive. At one extreme there is the very slow
and gentle perturbation of the binary known as the Kozai mechanism. Then only
the binary eccentricity and inclination change periodically while the orbital sizes
are unaffected. At the other extreme we have a high speed intruder which gives the
binary a ‘shock treatment’ during its brief encounter with the binary. In between, a
binary is strongly perturbed at close encounters but is able to maintain its identity
and not break up or exchange members. A stability boundary is derived which tells
us where the perturbation treatment ends and the chaos theory begins.
In this way we can give a rather complete coverage of the astrophysical
three-body problem. As in the case of the chaos theory of Chapters 7 and 8, also
in the perturbation theory of Chapters 9 and 10, a great deal of space is dedicated
to comparing numerical results with the theory. This is necessary since it is not
always obvious, in the absence of exact theory, what approximations should lead to
the best understanding of the experiments. Often we even find that in the final steps
we just have to accept the guidance of the numerical experiments without clear
justification of the theory. This is not because the theories could not be pushed any
further but more because we like to keep the theory at a rather simple level (and it
may appear quite complicated to some readers already as it is). But also we have to
remember that in the general three-body problem with strong interactions no exact
theory exists, and we should not spend too much effort towards this elusive goal.
Throughout the development of the theory we will look at some small astrophys-
ical problems which are easily solved at this stage. In the final chapter a couple of
8 Astrophysics and the three-body problem
larger issues are discussed which require solutions from different parts of the text
as well as other astrophysical information. To give the reader an idea of what sorts
of problems we are dealing with we outline a couple of astrophysical examples in
the next sections.
Figure 1.3 Halley’s comet, the most famous periodic comet, photographed in May
13, 1910. The big round object is Venus and the stripes are city lights of Flagstaff.
Image Lowell Observatory.
to change the orbit either to a hyperbolic escape orbit or to a more strongly bound
short period elliptical orbit. Basically we need to solve the three-body problem
consisting of the Sun, a planet and a comet. Since the comet has much smaller
mass than the other two bodies, the problem is restricted only to the question of
the motion of the comet. Generally we may assume that the planet goes around the
Sun in a circular orbit, and remains in this orbit independent of what happens to
the comet. This is an example of the restricted circular three-body problem which
will be discussed in Chapter 5.
It is a straightforward procedure to make use of an orbit integrator, a computer
code which calculates orbits, and to calculate the orbit of an Oort Cloud comet
through the Solar System, past various planets, perhaps through millions of orbital
cycles, until the comet escapes from the Solar System, or until it collides with the
Sun or one of the planets, or until the comet disintegrates. The calculation can be
long but it is possible with modern computers. However, questions remain: how
representative is this orbit and how accurate is the solution? Indeed, we do not know
the exact starting conditions for any Oort Cloud comets and neither can we carry
out the calculation over the orbital time of millions of years without significant loss
of accuracy.
To some extent these problems can be avoided by using the method of sampling.
We take a large sample of orbits with different initial conditions and calculate them
through the necessary period of time. We may then look at observational samples
10 Astrophysics and the three-body problem
and compare them with samples obtained by orbit calculation. In practice we need
millions of orbits in order to get satisfactory statistics of the captured short period
comets. It becomes a major computational challenge.
We can learn quite a lot by studying the interaction of a comet with only one
of the planets. A typical Oort Cloud comet has a high inclination relative to the
plane of the Solar System. It comes from one side of the plane (say, above), dives
through the plane, turns around below the plane, and crosses the plane again. A
strong three-body interaction takes place only if a planet happens to be close to
one of the crossing points at the right time. Most likely no planet is there at all, but
when a close encounter happens, almost always it is only with one planet only at
one of the two crossing points. Therefore we have a three-body problem.
But even then the exact nature of the three-body encounter is unclear. The best
we can do is to develop a statistical theory of how the comet is likely to react to
the presence of the planet close to the crossing point. Öpik (1951) developed such
a theory where the comet is assumed to follow an exact two-body orbit around the
Sun until it comes to the sphere of gravitational influence of the planet. At this
point the comet starts to follow an exact two-body orbit relative to the planet. After
leaving the sphere of influence of the planet, the orbit is again a (different) two-body
orbit around the Sun. A theory along these lines will be discussed in Chapter 6,
together with some more recent work. It is the most basic form of a solution to
the three-body problem. Notice that we will be discussing probabilities; this is the
recurring theme of the solutions of astrophysical three-body problems.
The statistical properties of comets which we will have to confront are primarily
the distribution of their orbital sizes (semi-major axes), their perihelion distances
and the distribution of orbital inclinations relative to the Solar System plane. In
terms of orbital sizes, the comets can be classified as being either Jupiter family
comets (orbital period below 20 years), Halley type comets (period above 20 but
below 200 years) or long period comets (longer periods than 200 years). Oort Cloud
comets have orbital periods in excess of a million years. Within these groups, the
inclination distributions vary (Fig. 1.4). One of the aims of a successful theory of
short period comets is to explain their observed numbers (28 Halley type, 183 Jupiter
family comets; Marsden and Williams 1999) in relation to the rate of comets coming
from the Oort Cloud. Also it should explain how the differences in the inclination
distributions originate.
If it turns out that a decrease in the semi-major axis of a comet is associated
with a decrease in its inclination (and this will be shown in Chapter 6) then we
might propose that comets evolve from Oort Cloud comets to Halley type comets
to Jupiter family comets by successively decreasing the semi-major axis and the
inclination of the orbit. How this works out in detail will be discussed in Chapter 11.
1.3 Short period comets 11
Inclinations in comet families
45
40
Jupiter family
35
30
number
25 Oort Cloud
20
15
10
5 Halley types
0
0 20 40 60 80 100 120 140 160 180
inclination
Figure 1.4 The inclination distributions in three classes of comets: Oort Cloud
comets, the orbits of which agree with random orientations (number proportional
to sin i), Halley type comets which are more often in direct orbits than in retrograde
orbits, and Jupiter family comets which are in direct orbits close to the Solar System
plane.
Another idea which is as old as the Oort Cloud theory is the concept of a disk
of comets left over from the formation of the Solar System. This is usually called
the Kuiper Belt or Edgeworth–Kuiper Belt (Edgeworth 1949, Kuiper 1951). These
comets reside mostly by the orbit of Neptune or beyond (Transneptuneans). Comets
do not develop a tail at this distance from the Sun and they can be observed only if
the nucleus of the comet (the solid body from which the tail originates) is quite large,
greater than about 100 km in diameter. Therefore our knowledge of the comets at the
Kuiper Belt is limited to the very largest bodies among them. By extrapolation it has
been deduced that there may exist as many as 107 ‘ordinary’ comets in the Kuiper
Belt and in its vicinity (Levison and Duncan 1997), even though only somewhat
over 300 very bright ones are known. Some of the Transneptunean comets certainly
approach the Sun at some point in their orbital evolution and become visible by
their bright tails; however, it is difficult to know exactly which comets have this
origin.
The Oort Cloud comets are very loosely bound to the Solar System and therefore
the Oort Cloud may require replenishment. A more tightly bound and dynamically
more robust cloud of comets has been suggested to lie inside the Oort Cloud, the
so-called Inner Oort Cloud (Hills 1981). Numerical simulations of the origin of
the Oort Cloud by ejection of cometary bodies from the region of the outer planets
also suggest the existence of the Inner Oort Cloud (Duncan et al. 1987). Passing
stars and molecular clouds perturb the orbits of the Inner Oort Cloud comets and
12 Astrophysics and the three-body problem
Figure 1.5 In the course of years the components of the binary star Krüger 60 are
seen to orbit their common centre of mass. Image Yerkes Observatory
some of them end up in the Oort Cloud proper. The Inner Oort Cloud is yet another
potential source region for short period comets.
log[dN/d(logP)]
1
0.8
0.6
0.4
0.2
0
-2 -1 0 1 2 3 4 5 6
log(P/yr)
Figure 1.6 The distribution of binary star periods P, plotted per logarithmic inter-
val d log P (Abt 1977). The horizontal dashed line represents the so-called Öpik
law.
Because of the wide range of physical sizes of binary systems, studies of their
distributions require different methods on different scales. Various selection effects
also complicate these studies. When all the studies on different scales and with
corrections for selection effects have been combined, the following results have
been obtained. The period distribution is flat over five orders of magnitude when
plotted per logarithmic interval of the period P (Fig. 1.6; Abt 1977). A distribution of
this type is sometimes called ‘scale free’. Written in the usual way the distribution is
f (P) dP ∝ P −1 dP (1.1)
which corresponds to
in terms of the binary orbital energy E B . This form of binary energy distribution
is known as Öpik’s law (Öpik 1924). One of our aims is to explain the origin of
this law; we will come back to it in Chapter 11.
The distribution of orbital eccentricities e depends on the period: short period
orbits have low eccentricities while long period orbits are more eccentric with the
weight of the distribution at the high end of the range of e (Fig. 1.7; Duquennoy
and Mayor 1991). In the latter case the data are consistent with
4
P < 3 yr
0
0 0.3 0.6 0.9
e
Figure 1.7 The observed distribution of eccentricities in a sample of binary stars.
The data for binaries of longer orbital periods (over 3 years) are shown by points
with error bars while the data for binaries of short orbital periods (less than 3 years)
are indicated by the broken line (Duquennoy and Mayor 1991). The straight dotted
line represents a possible interpretation of the long orbital period data.
operates most strongly during the pericentre of the orbit, when the two stars are
closest to each other. The resulting loss of orbital speed at the pericentre leads to
a rather circular orbit. But the low eccentricity orbits are dominant even among
binaries whose orbital size is as large as one astronomical unit (AU), the distance
between the Earth and the Sun. Tidal friction does not play any role at these sepa-
rations since the radius of the Sun (and a typical main sequence star) is only about
0.005 AU. We have to look for other explanations for the different eccentricity
distributions.
The mass ratio distribution depends on the type of the primary star (Abt 1983).
The most massive stars (primaries of spectral type O) tend to be in rather equal
mass pairs, m 2 /m 1 >
∼ 0.3 while the reverse is true for less massive upper main
sequence stars. For the low mass, solar type primaries the mass ratio peaks at about
m 2 /m 1 = 0.25 (Duquennoy and Mayor 1991).
To some extent the differences can be understood by the form of the mass dis-
tribution function for single stars:
f (m) ∝ m −α (1.4)
where α = 2.35 over most of the range (Salpeter 1955). If the more massive com-
ponent is taken to be m 1 = 1, then the distribution for m 2 should follow Eq. (1.4),
and the same should also apply to the mass ratio m = m 2 /m 1 . However, obser-
vations give a less steep function. For example, Kuiper (1935a, b) found that
f (m) = 2(1 + m)−2 gives a good fit to data on close binaries while Trimble (1990)
prefers f (m) ∝ m −1 . Thus the masses of the binary components are more equal
1.5 Groups of galaxies 15
Mass ratios in binary stars
1.8
1.6
1.4
log(number) 1.2
1
0.8
0.6
0.4
0.2
0
0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
m_2/m_1
Figure 1.8 The distribution of the mass ratio m 2 /m 1 (m 2 < m 1 ) in a B spectral
type binary star sample (points with error bars; Evans 1995). In comparison, the
dashed line gives the expected distribution if the companion of the B-type star has
been picked at random among ordinary stars lighter than the primary.
than we would predict by picking a companion for the primary at random from the
single star mass function (Fig. 1.8; Evans 1995).
On the other hand, the mass function for low mass stars (below the mass of the
Sun) has a lower power-law index α 1.25. Therefore the mass ratio distribution is
also expected to be less steep, f (m) ∝ m −1.25 . This agrees with the observed flatter
distribution of m for solar type primaries. Quantitative comparisons are difficult
because the single star mass function below about 0.5 solar mass is poorly known.
At the other end, the mass function for O-type stars is steeper than Eq. (1.4)
with α = 2.35. A power-law index α = 3.2 may be used at the high mass end,
above m 10M (Mihalas and Binney 1981). This leads to an even steeper mass
ratio distribution, in contrast with observations which go in the opposite direction.
Thus there is much more to the mass ratio distribution than simply picking a pair at
random. Even if the pairs might have formed like this initially, there has obviously
been strong evolution which has modified the mass ratios. We will discuss the
three-body process in Chapter 11 as a possible cause of this evolution.
Figure 1.9 Coma Cluster of galaxies. Picture taken with the KPNO 4-m Mayall
Telescope. Image NOAO/AURA/NSF.
requires an understanding of the multiple systems. First and foremost, the total
masses of galaxies, including their dark halos, are best derived from the study of
multiple systems.
The method which is commonly used is to assume that the multiple system
is in a bound state, i.e. that none of the members is currently escaping from the
grouping. In addition, some assumptions are made about the nature of the orbits
(e.g. circular, radial or intermediate) as well as about projections onto the plane of
the sky. With these assumptions it is possible to connect observed properties, i.e.
projected separations and radial velocities with inferred quantities such as the total
mass of the system.
The total amount of light emitted by all the member galaxies is also measurable.
Then the calculated mass (in solar mass) may be divided by the measured light
output (in solar luminosity), and the resulting number is called the mass-to-light
ratio. It is usually assumed that it is rather narrowly confined to a value around
30–50, at least for spiral galaxies. It is still an unsolved problem in astronomy as to
how frequently galaxies in a complete volume sample fall into this range of mass-to-
light ratio, and to what extent there is variation from one class of galaxies to another.
1.6 Binary black holes 17
Figure 1.10 Radio galaxy Cygnus A, as seen in radio waves. The two bright
emission regions straddle the galaxy which is centred on the spot of radio emission
at the centre of the picture. Image courtesy NRAO/AUI.
scale is very long, of the order of the current age of the Universe (Hubble time).
Occasionally the binary merger may happen in a shorter time scale through an
increase of eccentricity in the binary, and a collision of black holes at the pericentre
of the binary orbit (Aarseth 2003b). The problem of calculating this time scale will
be discussed in the final chapter using three-body theory.
Multiple mergers of galaxies also take place, and these may produce triple or
quadruple supermassive black hole systems. The evolution of the triple systems
is also discussed in the last chapter. Then we will find that supermassive black
holes should occasionally fly out of their parent galaxies, or sometimes remain as
oscillators in the galaxy. These are processes not well understood at present since
no unique identification of a supermassive black hole outside a galaxy has yet been
made. However, the currently popular CDM cosmology leads us to expect that
the black hole escape process was extremely common in the early universe when
large galaxies were assembled from their small progenitors.
Supermassive black holes are most likely associated with strong radio sources,
radio galaxies and quasars. Radio galaxies typically have a double lobe structure
where the regions of strongest radio emission lie on diametrically opposite sides of
the galaxy (Fig. 1.10). Two categories of theories are usually advanced to explain
the double radio source phenomenon: beam theory (Blandford and Rees 1974) and
slingshot theory (Saslaw et al. 1974). According to the former theory, supermassive
black holes remain stationary in the centres of galaxies, or at most move around
each other in binary orbits, and send out beams of relativistic plasma in two oppo-
sitely directed beams. These beams supposedly reach outside the galaxy and cause
radio emission there. In the slingshot model supermassive black holes themselves
1.6 Binary black holes 19
are thrown out of the centre of the galaxy via the three-body process. Due to con-
servation of linear momentum, a single black hole escapes in one direction and the
binary in the opposite direction from the galaxy, and each produces radio emission
along its escape path.
As long as no supermassive black holes have been directly observed or their
existence in the radio lobes has not been disproven, it is difficult to verify either
theory. However, there are indirect ways to do so. The most clear-cut test is related
to the escape direction of radio lobes. In the beam theory the direction is along the
rotational axis of the disk of gas (the accretion disk) which surrounds the black hole.
Even though the disk usually cannot be imaged directly, it is a fair assumption that
it arises from a more extended disk of gas in the galaxy via accretion. Therefore the
observed orientation of the extended gas in the galaxy should also give information
about the accretion disk of the central black hole. The projection of a circular disk
in the sky is an ellipse; the direction of the minor axis of the ellipse shows the
projection of the axis of symmetry of the disk.
Another way to find the orientation of the accretion disk is to observe small
scale jets which emanate at right angles to the disk. Again the jet line shows the
projection of the axis of symmetry in the sky. In the beam theory the radio lobes
should lie further along the jet line, and the projected angle in the sky between the
jet line and the double lobe radio axis should be zero. Thus the radio axis should
be in the direction of the minor axis of the image of the extended gas.
It is also possible that the distribution of stars in the galaxy is spheroidal and
flattened along the same axis as the gas disk. This is certainly true of spiral galaxies,
but may also pertain to elliptical galaxies which are hosts of double radio sources.
The alignment of the gas disk and the stellar disk could have its origin in mergers of
galaxies: the angular momentum of both components should be strongly influenced
by the orbital angular momentum of the binary galaxy which has merged.
In contrast to the beam theory, the slingshot theory predicts the escape of radio
lobes along or close to the central plane of the galaxy, and perpendicular to the
minor axis of the gas distribution and perpendicular also to the jet directions. In
Chapter 7 we will discuss the escape directions in the three-body problem in more
detail, and also compare them with observations.
2
Newtonian mechanics
This chapter introduces the basic concepts of Newtonian mechanics. We will empha-
size the areas which are most useful in the three-body problem, and also familiarise
ourselves with a system of units and scaling laws. The calculation of orbits using
Newton’s laws is a central theme of this book, and therefore a brief introduction to
the methods follows. It is not the purpose of this work to teach the latest orbit cal-
culation techniques; therefore only a brief introduction is given. Finally, we discuss
the connection of Newtonian mechanics to chaos. It may come as a surprise that
the introduction of just one more body to the well behaved two-body system brings
about a chaotic, unpredictable dynamical system. This was realised by Poincaré
well before the concept of deterministic chaos became a popular topic.
First law If there are no external forces, an object will maintain its state of motion,
i.e. it will stay at rest or continue rectilinear motion at constant velocity.
20
2.2 Inertial coordinate system 21
In more modern language this means that in the presence of several forces F 1 , F 2 ,
. . . , F n the body will behave as if there were only a single force F, which is the
vector sum of the individual forces: F = F 1 + F 2 + · · · + Fn .
These laws determine how objects will move in the presence of forces. Another
law was discovered by Newton to provide the force:
Law of gravity If the bodies A and B have masses m A and m B , respectively, and
if their mutual distance is r , A will act on B with a force that is directed towards A
and has a magnitude
Gm A m B /r 2 , (2.2)
where G is a constant, the constant of gravity, the value of which depends on the
units chosen.
When we calculate the gravitational force or acceleration due to an object, or
the effect of gravity on a body, we assume that all the mass is concentrated in the
centre of mass, i.e. the objects are pointlike. This is valid also for real bodies if they
are spherically symmetric, which is true for many celestial objects.
This is almost all one needs to know about physics; from here on, mainly math-
ematical methods are required.
where
rij = ri − r j,
(2.8)
ri j = |r i − r j |.
r i = xi ê x + yi ê y + z i ê z . (2.9)
and
1 n n
1
∇i V = − G m j m k ∇i .
2 j=1 k=1,k= j r jk
Since r jk depends only on the coordinates of the bodies j and k, only the terms
with j = i or k = i are non-zero:
n n
1 1 1 1
∇i V = − G m i m k ∇i − G m j m i ∇i
2 k=1,k=i rik 2 j=1, j=i r ji
n
1
=−G m j m i ∇i ,
j=1, j=i
ri j
where
1
∇i = ∇i ((x i − x j )2 + (yi − y j )2 + (z i − z j )2 )−1/2
ri j
(xi − x j )ê x + (yi − y j )ê y + (zi − z j )êz
=−
ri3j
rij
=− 3.
ri j
Then immediately
n
rij
∇i V = G mi m j . (2.12)
j=1, j=i ri3j
Comparing this with the equation of motion (2.7) one may write
where
1 n n
m j mk
V =− G . (2.14)
2 j=1 k=1,k= j r jk
24 Newtonian mechanics
n
mj
U (r ) = −G (2.15)
r − r j
j=1
has the dimension of energy per unit mass, and the potential energy of a mass m at
r is U (r )m. The equation of motion of this object is
r̈ = −∇U. (2.16)
where dV is the volume element and the integration is extended over the whole
mass distribution. Evaluation of the potential of an arbitrary object often leads to
complicated integrals.
Example 2.1 Find the potential of a homogeneous disk with radius R at a distance
z along the axis of the disk.
Let the surface density of the disk be ρ (kg/m2 ). Then its total mass is ρπ R 2 . In
cylindrical coordinates (r, φ) the potential at a distance z is
R 2π
ρr dφ
U (z) = −G dr √
0 0 z2 + r 2
R
r dr
= −2π Gρ √
0 z2 + r 2
R
= −2π Gρ z2 + r 2
0
√
= −2π Gρ z2 + R2 − z2
R2
= −2π Gρ √
z2 + R2 + z
GM 2
=− .
z 1 + (R/z)2 + 1
2.5 Constants of motion 25
If the distance z is much greater than the diameter of the disk, we can replace the
square root with the two first terms of its Taylor expansion:
GM 1
U (z) ≈ −
z 1 + 4 (R/z)2
1
GM 1 R 2
≈− 1− .
z 4 z
The potential of a point mass would be −G M/z. It is the shape of the mass
distribution that gives rise to the second term. Thus we cannot replace such an
object with a point mass.
The gradient of the potential gives the force (per unit mass). In this case only the
component parallel to the z axis remains. Using the previous series approximation,
2
dU GM 3 R
Fz = − ≈ 2 1− .
dz z 4 z
n
n
n
rij
m i r¨i = −G mi m j .
i=1 i=1 j=1, j=i ri3j
26 Newtonian mechanics
The sum on the right hand side consists of pairs (r i − r j )/ri3j and (r j − r i )/r 3ji .
These terms cancel each other, and hence
n
m i r¨i = 0.
i=1
where a and b are constant vectors. Comparison with the definition of the centre
of mass (2.4) shows that the sum on the left hand side is M R. Thus we have
at + b
R= . (2.19)
M
Thus the centre of mass or barycentre moves at constant velocity along a straight
line. The obvious consequence of this is that the barycentric coordinate frame is an
inertial frame.
The total angular momentum of the system is
n
L= m i r i × r˙i . (2.20)
i=1
Again the terms on the right hand side cancel each other pairwise, because the
vector product is anticommutative. Thus
L̇ = 0, (2.21)
Finally the total energy of the system is derived. Taking the scalar product of the
equation of motion of the ith body and r˙i , and adding all equations, one obtains:
n
n
m i r¨i · r˙i = − r˙i · ∇i V
i=1 i=1
n
∂V ∂V ∂V
=− x˙i + ẏi + z˙i
i=1
∂ xi ∂ yi ∂z i
dV
=− .
dt
Integration with respect to time gives
1 n
m i r˙i · r˙i = −V + E, (2.22)
2 i=1
where E is a constant of integration. The left hand side of the equation gives the
kinetic energy of the system
1 n
T = m i vi2 , (2.23)
2 i=1
E = T + V, (2.24)
Consider the time average in the interval t ∈ [0, τ ]. Denoting this average value
by ,
1 τ n
Ȧ = Ȧ dt = 2T + Fi · r i . (2.26)
τ 0 i=1
If the system remains bounded, its virial cannot increase without limit. Thus the
integral in (2.26) remains finite. When the length of the interval increases, Ȧ
approaches zero, and thus
n
2 T + Fi · r i = 0. (2.27)
i=1
This is the virial theorem in its general form. If the forces Fi are due to mutual
gravitational forces only, the sum becomes
n
n
n
ri − r j
Fi · r i = −G mi m j · ri.
i=1 i=1 j=1, j=i
ri3j
It is possible to write the double sum using only the vector differences r i − r j . We
start by writing the sum twice but by changing the order of the summation indices
i and j. We add the two identical sums and divide by two:
n
Fi · r i
i=1
1 n n
ri − r j r j − ri
=− G mi m j · r i + m j mi ·rj
2 i=1 j=1, j=i ri3j r 3ji
1 n n
ri − r j
=− G mi m j · (r i − r j )
2 i=1 j=1, j=i ri3j
1 n n
mi m j
=− G .
2 i=1 j=1, j=i ri j
n n
The double sum i=1 j=1 , where i = j, gives every term twice, due to the
symmetry relative to interchange between i and j. In order to have every term only
2.7 The Lagrange and Jacobi forms of the equations of motion 29
once, the latter sum should start from j = i + 1. We make this change and multiply
the above expression by two:
n n n
mi m j
F i · r i = −G = V.
i=1 i=1 j=i+1
ri j
Thus
1
T = − V. (2.28)
2
The virial theorem (2.28) can be used to study, for example, the stability of
clusters of stars. The kinetic energy can be determined from the velocity dispersion
and the potential energy from the mass distribution derived from the brightness of
the objects. If the kinetic energy is much higher than the average value given by the
virial theorem, the system is unstable. The virial theorem indicates that clusters of
galaxies should be unstable, which contradicts observational evidence about their
stability. This is one reason to believe that such systems must contain large amounts
of dark matter.
r r31
1 r3
rB r1
Figure 2.1 Coordinates of the Jacobi system: the bodies 1 and 2 form a binary
while the third body (3) is more distant.
The corresponding equations for the acceleration between bodies 2 and 3 as well
between bodies 3 and 1 are obtained from (2.31) by changing the subscripts ac-
cordingly.
In the Jacobi system of coordinates, the three-body system is considered to be
composed of two parts: a clearly defined binary, and a somewhat distant third body
(Fig. 2.1). The centre of mass of the binary is used as a reference point whose
position is r B . The relative position vector of the binary is r 21 = r . By definition
of the centre of mass (Eq. (2.4))
1
rB = (m 1r 1 + m 2r 2 )
m1 + m2
m1 m2 m2 m2
= r1 + r2 + r1 − r1 (2.33)
m1 + m2 m1 + m2 m1 + m2 m1 + m2
m2
= r1 + r
m1 + m2
or, alternatively,
m1 m2 m1 m1
rB = r1 + r2 + r2 − r2
m1 + m2 m1 + m2 m1 + m2 m1 + m2
m1 (2.34)
= r2 − r.
m1 + m2
Let us call the position vector of body 3 relative to the centre of mass R3 = r 3 − r B .
Then
m2
r 31 = r 3 − r 1 = R3 + r B − r 1 = R3 + r,
m1 + m2
m1 (2.35)
r 32 = r 3 − r 2 = R3 + r B − r 2 = R3 − r.
m1 + m2
2.8 Constants of motion in the three-body problem 31
Also
m1 m2
r 31 + r 32
m1 + m2 m1 + m2
m1 m2 m 1m 2 m 1m2
= + R3 + − r
m1 + m2 m1 + m2 (m 1 + m 2 )2 (m 1 + m 2 )2
= R3 . (2.36)
Using Eqs. (2.31) and (2.36), the equations of motion are:
r r 31 r 32
r̈ = −G (m 1 + m 2 ) 3 + m 3 3
− 3 (2.37)
r r31 r32
and
m1 m2
R̈3 = r̈ 31 + r̈ 32
m1 + m2 m1 + m2
m1 m 1 r 31
= G m2 W−M
m1 + m2 m 1 + m 2 r313
(2.38)
m2 m 2 r 32
+ m1 (−W) − M
m1 + m2 m 1 + m 2 r32
3
m 1 r 31 m 2 r 32
= − GM + .
m 1 + m 2 r31
3 m 1 + m 2 r323
L = m 1r 1 × ṙ 1 + m 2r 2 × ṙ 2 + m 3r 3 × ṙ 3 . (2.41)
32 Newtonian mechanics
In the Lagrangian form the angular momentum in the centre of mass coordinate
system reads
m 1 m 2 m 3 r 12 × ṙ 12 r 23 × ṙ 23 r 31 × ṙ 31
L= + + (2.42)
M m3 m1 m2
and in the Jacobi form
n
ṙ i j r i j (r i j · ṙ i j )
Ḟ = −G mj 3
−3 ,
j=1, j=i ri j ri5j
n
r̈ i j ṙ i j (r i j · ṙ i j )
F̈ = −G mj 3
−6 (2.54)
j=1, j=i
ri j ri5j
r i j (r i j · ṙ i j )2
+3 5 5 − ṙ i j · ṙ i j − r i j · r̈ i j .
ri j ri2j
2.11 Integration of orbits 35
1 1
ṙ 1 = ṙ + F t + Ḟt 2 + F̈t 3 + · · ·. (2.55)
2 6
By repeating the short time step for all of the bodies we may hope to calculate the
time evolution of the whole system.
There are obvious problems with this scheme. Complicated calculations are car-
ried out at every time step which consume computer time and which cause round-off
errors. Thus, in order to finish the calculation in a reasonable amount of computer
time, the time step should not be very small. But then the convergence of the Taylor
series becomes poorer. We realise that during a relatively long time step the acceler-
ation can change considerably from what it was at the beginning of the step, because
in reality all bodies move simultaneously. Some form of prediction of motions of
all bodies has to be introduced prior to taking the time step. The prediction can be
based on a low order Taylor series which does not require complicated calculations:
1 1
r 1 = r + ṙ t + F t 2 + F̂ t 3 . (2.56)
2 6
where t−1 is the length of the preceding time step, and F −1 is the value of accel-
eration at the beginning of that step. In general, the computer algorithm becomes
faster when the derivatives (2.54) are constructed from the corresponding back-
ward differences. The values of F only need to be stored for the several previous
steps, and the derivatives of F are then calculated by simple operations from these
quantities. In particular, the prediction makes use of the values of the acceleration
at two successive points where the acceleration has to be calculated in any case,
and thus the prediction requires only modest extra calculation effort.
Other features which save computer time are the use of variable time step and
individual time step. At some points in the orbit a body may be advanced by
long time steps, at other more critical times a small step is necessary to preserve
accuracy. Every body has a different situation in this regard; therefore it is useful to
assign a different, individual, time step to each of them. The prediction of the orbits
of all bodies up to a common moment of time then becomes crucial. When the
acceleration and its derivatives change rapidly, it is necessary to go to a short time
step. Therefore one may calculate the length of the step by requiring that the ratio
of successive terms in the Taylor series is of the order of a small number η 1.
36 Newtonian mechanics
For example,
1
| Ḟ|t 3
6
= η,
1
2
|F|t 2
1 r 2 µ
h= − ,
2 r2 r (2.62)
r 1
e = − − 2 r × r × r .
r µr
2.11 Integration of orbits 37
where
r · r r r r
= 2· r −µ +r f
2
r2 r r r
r r 2 µ
= 3 − 2 r + r · f .
r r
Thus
h = r · f . (2.64)
After some calculations we find also
1
e = − r (r · f ) + f (r · r ) − 2r h (2.65)
µ
(Problem 2.11). Therefore a regularized equation of motion has been derived which
has no singular terms when r → 0 as long as f → 0. Quantities h and e arise again
in the next chapter where their physical significance is clarified (Eqs. (3.11) and
(3.13)).
The above regularisation scheme (Burdet 1967, Heggie 1973) is not the only
possibility. Even though the Burdet–Heggie regularisation is known to function
perfectly well in numerical solutions of the N -body problem (Saslaw et al. 1974,
Heggie 1975, Valtonen and Heggie 1979), actually another transformation called
the K-S regularisation is more widely used (Kustaanheimo 1964, Kustaanheimo
and Stiefel 1965). However, the important point to us is that there are methods to
handle close two-body encounters without slowing down the orbit calculation and
without loss of numerical accuracy.
How does one then measure the accuracy of the orbit calculation? One may
test the method with a problem where the analytic solution is known, such as
the two-body problem (Chapter 3, Problem 3.8). Alternatively, one may fol-
low the conservation of quantities which are known to be strictly constant in
the N -body problem, i.e. the centre of mass and momentum, the total energy
and angular momentum. Perhaps the most stringent test is based on the time-
reversability of the orbits. If at any point along the orbit the velocities of all bodies
are reversed, the bodies should retrace their trajectories exactly. Computational
38 Newtonian mechanics
errors make this unlikely to happen in numerical orbit integration. The ability
to retrace the orbit backward is a very demanding test for the accuracy of the
method.
L 20 |E 0 |/G 2 m 50 (2.66)
is a dimensionless quantity. Here m 0 signifies the average mass of the three bodies,
to be defined more specifically in Chapter 7 (Eq. (7.28)). If this quantity is given
the value 6.25,
G = 1 ties the units such that corresponding to pc and km/s the mass unit is
1/(4.3 × 10−3 M ) = 233M . Since
1 km/s ≈ 1 pc/Myr (2.69)
the corresponding unit of time is very nearly one megayear (1 Myr = 106 yr).
Obviously the calculation can be carried out in other units but this does not bring
a new dimension to the three-body problem. In practice, the time is frequently
given in the natural time unit of the three-body system which is the crossing time
Tcr = G M 5/2 /(2|E 0 |)3/2 (see Eq. (8.26) and the discussion there). This is roughly
the time that it takes a body to travel through the three-body system.
In principle, the three-body problem may be solved by going through the 11-
dimensional phase space with a fine tooth comb, and by calculating the orbital
evolution at each phase space point. Let us briefly consider how practical this
proposal is. Construct a grid in the phase space composed of 100 values in each
coordinate and calculate the orbits at each grid point. It is in total 10011 = 1022
orbits. If every orbit takes 1 second to calculate in a fast machine, the project would
last 1022 seconds, i.e. about 3 × 1014 yr. Even so, the solution would be rather rough
since a mesh of 100 grid points is far too coarse to cover the phase space adequately.
Therefore it is not possible to find a complete solution of the three-body problem by
‘brute force’, even if we have perfect integration methods and powerful computers
in use.
The alternative, an analytical solution for example in the form of a series has
been proposed and searched for. It has been estimated that 108 000 000 terms are
required in the Sundman series to represent a three-body solution with the accuracy
commonly used in the ephemeris calculations. The computer time requirement
would be impossible, and certainly the representation could not be said to be simple.
It is still possible that, in the future, computers that use mathematical manipulators
(e.g. Maple, Mathematica) may be used to identify important terms in long series
expansions.
In recent years it has become clear that solutions of the three-body problem can
be found and represented in simple form by using a combination of analytical and
computational methods. Analytical methods give the basic functional forms which
can be tested and improved by computer experiments. They also provide useful
physical insights into the problem. We will describe this avenue of research mainly
in Chapters 7 to 10. Before that, some computer experiments will be described
which are useful to develop the analytical theory.
“I have just been down,” answered he, “to the herb-market. I had
made my bargain, and was coming away, when I met one of old
Capito’s men, driving an ass laden with articles from the country. So
I asked if he was carrying a present to his master’s brother. He said
he had brought nothing for Lucius but a letter; and that he believed
its purport was to invite the two young ladies, to come out to-day
and enjoy the beauty of the season. I no sooner got this
information, than I ran hither as swiftly as my legs would carry me.
You can easily go out, as if by chance, to pay your respects to the
Patrician.”
“Ah, Sempronia!” sighed Sextus, “shall I approach you at last?—
What will she think when she sees me there?—Oh, how will she
speak to me?”
Here Dromo, perceiving that his young master was a good deal
confused, came forward and said, “From observing your looks last
night, when I was making a handle of yon barbarian to torture our
friend of the porch, I think you are a good-natured person, who
would not willingly bring any of us into trouble. The truth is, that
Licinius wishes my young master here to marry a certain lady, who
has already had wet eyes over the ashes of a first husband; but who
is of noble birth, and very rich. Now Sextus, being only eighteen,
does not like this great lady so well as she likes him; and has, in
fact, lost his heart elsewhere.”—“Dromo,” answered I, taking young
Sextus by the hand as I spoke, “this is a pretty common sort of
story; but I shall take no side till I have seen both of the ladies; and
the sooner your ingenuity can bring that about, the more shall I be
beholden to you.”—“We shall try,” replied the slave, observing that I
had overcome the reluctance of the lover; “but in the meantime I
observe that the clients are beginning to assemble in the porch, to
await the forthcoming of Licinius. Go, therefore, and get some
breakfast, for, by and by, you will both be expected to accompany
the Senator to the Forum, to hear him plead; which, between
ourselves, will be a six hours’ job for you, unless you manage
matters dexterously.”
With this, arranging his gown, and putting himself into a dignified
attitude, he ordered the porter, who stood chained by the door, to
throw wide its massy valves; which being done, the litigants and
consulters, who were without, received the orator with acclamations,
and surrounded him on all sides. Some of the poorer ones, I
observed kissing the hem of his garment, and dodging wistfully at
his elbows, without ever attracting a word or look from him; while
those of a higher class came forward more familiarly, seeking to
impress particular circumstances upon his memory, and [pg
25]paying him compliments on the appearance he had made the day
before in the Centumviral Court. Encircled by this motley group, he
walked towards the great Forum, followed at a little distance by
Sextus, the preceptor, myself, and some freedmen of his household.
In moving on, we passed, by accident, the door of another great
pleader, by name Bruttianus, who stood there attended in a similar
manner. When he perceived Licinius, this man took from his door-
post a green palm-branch, and waved it towards us in a vaunting
manner; but our friend, saluting him courteously, cried out, with his
sharp and cutting voice, “We shall try it again.” Whereon
Xerophrastes, immediately stepping up to his patron, began thus,
“How this vain-glorious person exposes himself!—he is certainly a
weak man; and his tones, by Hermes, are more detestable than
those of an African fowl.”—At which words, Sextus tipped me the
wink; but I did not observe that Licinius was at all displeased with
them. Yet, soon after, Bruttianus having overtaken us, the
processions were joined, and the two pleaders walked the rest of the
way together in a loving manner, exchanging complimentary
speeches; to which Xerophrastes listened with edifying gravity of
visage.
The keen and lively gestures of the fervid Licinius, whose soul
seemed to speak out of every finger he moved, and who appeared
to be altogether immersed in the cause he pleaded, were succeeded
by the solemn and somewhat pompous stateliness of Bruttianus,
who made a brief pause between every two sentences, as if he were
apprehensive that the mind of the judge could not keep pace with
the stream of his illustrations, and looked round ever and anon upon
the spectators with a placid and assured smile, rather, as it seemed
to me, to signify his approbation of their taste in applauding him,
than his own pleasure in their applauses. Nevertheless, he also was
a splendid speaker, and his affectation displeased the more, because
it was evidently unworthy of his understanding. While he was
speaking, I observed that the Stoic preceptor was frequently shifting
his place among the crowd, and muttering every where expressions
of high contempt. But this did not disgust me so much as the fixed
attitude of ecstasy in which he listened to the discourse of his own
patron, and the pretended involuntary exclamations of his delight.
“Oh, admirable cadence!” he would say, “I feel as if I were draining a
honey-comb. Oh, harmonious man, where have I, or any other
person here, sucked in such sweetness!” These absurd phrases,
however, were caught up forthwith, and repeated by the numerous
young men who hung upon the skirts of [pg 28]the orator, and
seemed, indeed, to be drinking in nectar from the speech, if one
might judge from their countenances. From their taking notes in
their tablets from time to time, and from the knowing looks they
assumed at the commencement of every new chain of argument, I
guessed that these might be embryo jurisconsults, preparing
themselves by their attendance for future exertions of the same
species; and, indeed, when I listened to their conversation at the
close of every speech, I thought I could perceive in their tones and
accents, studied mimicry of the natural peculiarities of Licinius,
Bruttianus, and the other orators. Altogether, the scene was as full of
amusement as of novelty, and I could willingly have remained to the
end of the discussion. But my eyes chanced to fall upon young
Sextus, and I could not but see that his mind was occupied in
matters remote from the business of the Forum. He stood with his
arms folded in his gown, and his eyes fixed upon the ground, only
lifting them up from time to time with an impatient air towards a
side entrance, or to observe by the shadows on the porticos what
progress the sun was making.
[pg 30]
CHAPTER IV.
A lofty wall protected the fields of this retirement from the intrusive
eyes of passengers. We entered by a small side-door, and found
ourselves, as if by some magical delusion, transported from the glare
of a Roman highway, into the depth and silence of some primeval
forest. No nicely trimmed path conducted our feet. Every thing had
at least the appearance of being left as nature had formed it. The
fern rustled beneath us as we moved; the ivy was seen spreading its
careless tresses from tree to tree; the fawn bounded from the
thicket. By degrees, however, the gloom lessened, till at length, over
an open space of lawn, we perceived the porch of entrance, and a
long line of colonnade. We passed under the porch, and across a
paved court where a fountain was playing, into the great hall, which
commanded all the other side of the place—a noble prospect of
elaborate gardens gradually rising into shady hills, and lost in a
distance of impenetrable wood.
As we passed on, our host from time to time directed the attention
of his visiters, more particularly of the two Greeks, to the statues of
bronze and marble, which were placed at convenient intervals along
the terraces of his garden. The symmetry of these figures, and the
graceful simplicity of their attitudes, inspired me with I know not
what of calm and soothing pleasure such as I had never before
tasted, so that I thought I could have lingered for ever amidst these
haunts of philosophic luxury. The images were, for the most part,
portraits of illustrious men—Greeks, Romans—sages and heroes; but
beautiful female forms were not wanting, nor majestic
representations of gods and demi-gods, and all the ethereal
imaginations of the Grecian poets. Seeing the [pg 34]name of
Jupiter inscribed upon one of the pedestals, I paused for a moment
to contemplate the glorious personification of might and wisdom,
depositing, at the same time, a handful of roses at the feet of the
statue—upon which I could observe that my behaviour furnished
some mirth to the Epicurean Demochares; while, on the contrary,
Euphranor, the disciple of the Porch, approved of what I did, and
rebuked his companion for saying any thing that might even by
possibility disturb the natural piety of an innocent youth. But the
Roman Stoic stood by with a smile of stately scorn; and utter
indifference was painted on the countenance of Capito. At another
time, Sextus having staid behind to examine the beauties of a
certain statue of Diana, which represented the goddess stretched
out in careless slumber on the turf, with a slender grayhound at her
feet, the Epicurean began to rally me on having a taste inferior to
that of my friend, whose devotion, he said, could not be blamed,
being paid to an exquisite imitation of what the great Nature of
things had decreed should ever be the most agreeable of all objects
in the eyes of a person of his age.—“Whereas you,” continued he,
“appear to be more occupied with deep-hung eye-brows, ambrosial
beards, and fantastic thunderbolts, and the other exuberances of
Homeric imagination.”
The other, Athanasia, was not a dazzling beauty. Taller than her
cousin, and darker-haired, but with eyes rather light than otherwise,
of a clear, somewhat melancholy gray—with a complexion paler than
is usual in Italy, [pg 36]a demeanour hovering between cheerfulness
and innocent gravity, and attired with a vestal simplicity in the old
Roman tunic, and cloak of white cloth—it is possible that most men
might have regarded her less than the other. A single star of
diamonds, planted high up among her black hair, was the only
ornament she wore.
“In me,” he continued, “you behold one that has gone through the
experience necessary to produce an entire acquiescence in these
doctrines. I am one of those, Valerius, who have resolved to
concentrate, after this fashion, the whole of my dreams upon the
hour that is. There are not wanting, indeed, here and elsewhere,
persons who profess the same theories, only in the view of finding
excuse and shelter for the practice of vice. But till it be proved that
the practice of vice [pg 39]is the best means of enjoyment, in vain
shall it be asserted that our doctrine is essentially adverse to virtue.
The mistakes or the misdeeds of individuals must be estimated for
nothing; for where is the doctrine that may not be shewn to have
been defended by impure livers? The founder of our sect is
acknowledged, by its most virulent enemies, to have been the most
blameless of men, and they, I must take leave to believe, can never
be sincere friends of virtue, who doubt, that he who is a true
worshipper of pleasure, may also be the worshipper of virtue.”
They were both at a little distance behind the rest of the party, when
Euphranor addressed himself to me, saying, “Is not this youth, your
companion, the same that is under the guidance of a certain
Xerophrastes?”—“The same,” said I, “and a wary, sage-looking
Athenian is his tutor. I believe he also is of the Porch.”—“No doubt,”
interrupted Demochares; “he has a beard that Zeno might have
been proud of, and walks as if he conceived himself to be the chief
pillar of the Porch, if not the Porch itself.”—“Who shall prevent
Demochares from having his jest?” replied Euphranor. “The man is
by birth a Thessalian, and his gutturals still remind one strongly of
his native hills.”—I would gladly have heard more of it, but he was
interrupted by the nearer approach of the rest.
[pg 41]
CHAPTER V.
Before the hour of taking the bath, we exercised ourselves for some
time in the tennis-court, where I could not but admire, especially
after having heard Capito philosophise, the vigour and agility
displayed by him as well as his companions. I was then conducted
into the baths, where, after being washed and perfumed in the most
luxurious manner, I was arrayed in an elegant supper-garment by
one of the slaves of our host. At table we were joined again by the
ladies, who both reclined on the same couch with their uncle. Three
comely youths attended us, in short tunics, and girt with napkins of
fine linen; but, during the repast, an ancient female slave stood in
silence behind the couch of the young ladies. A small fountain of
alabaster played between two tall candelabra of the same material,
at the farther end of the apartment; and a young damsel stood
beside them, swinging slowly from time to time a silver censer, from
which clouds of delicate odour rolled up to the mirrored roof.
In all things the feast was splendid; but there was no appearance of
useless or vain ostentation. Every thing was conducted in a style of
great calmness and order, without the least formality. The repast
inter[pg 42]rupted not the conversation, which went on in a manner
to me equally instructive as entertaining; although I must confess
the presence of Athanasia sometimes rendered me inattentive to
what was spoken. I could not divest myself of the idea, that some
unknown circumstance was pressing on the mind of the fair
creature, and that when she smiled upon those who addressed her,
it was sometimes to conceal her ignorance of that which had been
said.
Being asked by Capito, I endeavoured, among other things, to
inform him and his friends, as far as I could, concerning the then
condition of this island, which, more particularly after the exploits of
Agricola, had come to be a subject of some interest. In return, the
chief topics which then occupied the capital were discussed by them,
as I perceived, in a great measure on my account; and I listened
with delight to the praises, which they all agreed on bestowing on
the new Emperor. Many anecdotes were narrated, which tended to
strengthen the feelings of admiration, with which I had already been
accustomed to contemplate his character. But others were told, as
the conversation went on, which I could not so easily reconcile with
the idea I had formed of him.
[pg 44]
“I hope,” replied Capito, “I shall never regard with indifference any
question, in which the interest of the empire and the honour of
Trajan are concerned. But if you mean only that I am indifferent
about the nature of this Syrian superstition, you are in the right. I
have no knowledge of its dogmas, nor desire to have. I presume
they have their share of that old eastern barbarity, in the shady
places of which the elder Greeks used to think they could discover
the outlines of something really grand and majestic.”
Capito said, he could not think of our going into the city that
evening, and that we must all make up our minds to remain in the
villa. The countenance of [pg 46]Sextus brightened up, and he
looked to me as if to ask my assent. I was easily persuaded, and our
host despatched a messenger to inform Licinius of the cause of our
absence. The old man then led us into another apartment, which
was richly furnished with books and paintings. Here he read for
some time out of one of the poets, to a party, none of whom, I am
afraid, were very attentive in listening to him, till, the hour of rest
being come, we were conducted to our several apartments, Sextus
and myself, indeed, being lodged in the same chamber.
After we had both retired to our beds, and the lights were
extinguished, we still continued for some time to talk over the
incidents of our visit, and the future prospects of Sextus and his
love; until at length sleep overpowered us in easy bonds, and
agreeable dreams followed, I doubt not, in the hearts of us both, the
thoughts and sights of a delightful day. Mine surely were delightful,
for they were all of Athanasia. Yet, even in the visions of the night, I
could never gaze [pg 47]on her face without some strange
impression of mystery. I saw her placid smile—I heard the sweet low
cadence of her voice—but I felt, and I could not feel it without a
certain indescribable anxiety, that her deep thoughts were far away.
[pg 48]
CHAPTER VI.
I awoke early, and drew near to the bed of Sextus; but seeing that
he was fast asleep, and that a quiet smile was on his lips, I could not
think of awakening him. The sun shone bright into the apartment,
and I resolved to walk forth and breathe the balmy air of the
garden.
The moisture was still heavy on the green paths, and the birds were
singing among glittering leaves; the god-like statues stood
unscathed in their silent beauty. I walked to and fro, enjoying the
enchantment of the scene;—a new feeling of the beauty of all things
seemed to have been breathed into my soul; and the pensive grace
of Athanasia hovered over my imagination, like some presiding
genius of the groves.
I found myself near the favourite grotto, and had stood over against
its entrance for some space, contemplating the augmented stream
as it fell from the superincumbent rock, and regretting the ravage
which the nightly tempest had made among the delicate flowers
round its basin. Twice I thought I heard the murmurs of a voice, and
twice I persuaded myself that it was only the rippling of the waters;
but the third time I was satisfied that some person must be near. I
passed [pg 49]between the water and the rock, and beheld the fair
creature that had been occupying so many of my thoughts, kneeling
far within the grotto, as it seemed, in supplication. To disturb her by
advancing farther, would have been impious; to retire, without the
risk of disturbing her, almost impossible; but I remained there fixed
to the spot, without perhaps considering all these things as I should
have done. The virgin modesty of her attitude was holy in my eyes,
and the thought never occurred to me, that I might be doing wrong
in permitting myself to witness the simple devotions of Athanasia.
“Great God, listen to my prayers,” was all I understood of what she
said; but she whispered for some moments in a lowly and fervent
tone, and I saw that she kissed something with her lips ere she
arose from her knees. She then plunged her hands into the well, by
whose brink she had knelt, and turned round to the light.
“Athanasia, forgive me,” was already on my lips; but on seeing me,
she uttered a faint cry and fell prostrate upon the marble. I rushed
forward, lifted up her head, and laved water from the fountain, till I
saw her lips tremble. At last she opened her eyes, and after gazing
on me wildly for a moment, she gathered her strength, and stood
quite upright, supporting herself against the wall of the grotto.
“Great heavens!” cried I, “in what have I offended, that I should be
rendered the cause of affliction to Athanasia? Speak, lady, and say
that you forgive me.”
I sat down beside the dark well, wherein she had dipt her hands,
and mused in a most disturbed mood on all the particulars of this
strange and unexpected interview. Every motion of her features—
every modulation of her voice, was present with me; I had gathered
them all into my heart, and I felt that I must cherish them there for
ever. From the first moment I saw her, my eyes had been
constrained to gaze upon her with an interest quite novel to me; but
now I knew that she could not smile, without making my heart faint
within me, and that the least whisper of her voice was able to bring
tears into mine eyes. Now I thought of my own unworthiness, and
could not help saying to myself, “Why should a poor ignorant
provincial, such as I am, be torturing myself with the thoughts of
such a creature as this?” Then, again, some benign glance of hers
would return before me, and I could not help having some faint
hopes, that her [pg 52]innocent heart might be won to me by
faithful unwearied love. But what always threw me back into despair,
was the recollection of the mystery that I knew hung over her mind,
although what it was I could not know. That she had been saying
something in her prayers which could not be overheard without
betraying her, she had herself confessed. What could be this secret,
so cherished in dread, and in darkness?—A crime?—No crime could
sully the clear bosom of her innocence. No consciousness of guilt
could be concealed beneath that heavenly visage. But perhaps she
had been made the confidante of some erring,—some unhappy
friend. Perhaps, in her prayer, she had made mention of another’s
name, and implored the pardon of another’s guilt. Last of all, why
might it not be so, that the maiden loved, and was beloved again;
that she might have some reason to regard any casual betrayal of
her affection as a calamity; and that, having uttered the name of her
lover in her secret supplications, her terrors might all have been
occasioned by her apprehensions of my having overheard it? And yet
there was something in the demeanour of Athanasia, that I could
not bring myself to reconcile entirely with any one of these
suppositions. Had she feared that I had overheard any confession of
guilt,—even of the guilt of another,—surely some semblance of
shame would have been mingled with her looks of terror. Had she
apprehended only the discovery of an innocent love, surely her
blushes would have been deeper, and her boldness less. Yet the last
solution of the difficulty was that which haunted me the most
powerfully.
When I came forth into the open air, I perceived [pg 53]that the sun
was already high in heaven, and proceeded in haste towards the
villa, not doubting that Sextus and Capito would be astonished by
the length of my absence. I found them and the ladies walking
under the northern colonnade, having returned, as they told me,
from a fruitless search after me through almost the whole of the
garden. I looked to Athanasia, as if to signify that she well knew
where I might have been found; but, although I saw that she
understood my meaning, she said nothing in explanation. Sextus
drew me aside shortly after, and told me, that his father had sent to
inform him, that our presence was necessary in the city before
supper-time, to attend a great entertainment which was to be given
that evening by the lady whose cause he had successfully pleaded in
the Forum on the preceding day; which lady, I now for the first time
learned, was no other than the same Marcia Rubellia, to whom his
father was very anxious the youth should be married. The success of
this pleading had increased very much the wealth of the lady, and,
of course, as Sextus well knew, the anxiety of Licinius for the
proposed union; and to remain at the villa any longer, was, he said,
entirely impossible, since he already suspected his father had not
been quite pleased with him for leaving the Forum the day before,
without staying to hear out a cause in which his duty, if not his
inclination, ought to have made him feel so greatly interested.
We bade adieu, therefore, to our kind host and the young ladies, not
without more reluctance than either of us durst express, and ready
promises to return soon again to the villa. We found Dromo and
Boto waiting [pg 54]for us at the gate, the former of whom looked
unutterable things, while the latter appeared to be as joyful in
seeing me again, as if we had been parted for a twelvemonth. The
two slaves were mounted on asses, but they led horses for our
conveyance; so we mounted with all speed, and were soon beyond
the beautiful enclosures of Capito. As soon as we were fairly out of
sight of the house, Dromo began to ply Sextus with innumerable
questions about the result of the visit, all of them in bad Greek; that,
as he said, there might be no chance of what passed being
understood by the Druid; for by that venerable designation, he
informed us, the primitive Boto had already come to be best known
in the vestibule of Licinius. “Ah!” quoth he, “there is no need for
many words; I am sure my young master has not been behindhand
with himself. If he has, it is no fault of mine, however. I put
Opportunity into his hands, and she, you know, as the poets say, has
only one lock of hair, and that is in front.”
Welcome to our website – the perfect destination for book lovers and
knowledge seekers. We believe that every book holds a new world,
offering opportunities for learning, discovery, and personal growth.
That’s why we are dedicated to bringing you a diverse collection of
books, ranging from classic literature and specialized publications to
self-development guides and children's books.
ebookbell.com