CAAO Tutorial
CAAO Tutorial
Armin Hodaei
2023
Preface
This guide is intended for students who wish to participate in the Canadian
Astronomy and Astrophysics Olympiad (CAAO). It serves as an introductory-
level introduction to the Olympiad for interested students across Canada. Each
year, the highest-achieving CAAO students are selected to represent Team
Canada in the International Astronomy and Astrophysics Olympiad (IOAA),
and are provided with a training program to prepare them for the international
competition.
While this guide contains a great deal of information, we recommend that
students supplement their learning with other resources listed in the reference
section. The guide includes numerous practice problems designed to complement
students’ learning path. Additionally, we highly encourage students to solve
past CAAO problems available in a separate file.
To participate in the Canadian astronomy Olympiad, students should have
a solid foundation in high school-level physics and mathematics. However,
for the international Olympiads, students will need to develop an advanced
understanding of physics and mathematics beyond what is typically taught in
high school.
The International Astronomy and Astrophysics Olympiad (IOAA) is a
prestigious international competition for high school students. Each year, the
brightest students from around the world compete in this event, and Canada
has been participating in the IOAA since 2013.
i
Resources
iii
Contents
Preface i
Resources iii
Contents v
List of Tables ix
I Introduction to CAAO 1
1 Basic Concepts 3
1.1 Parallax . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3
1.2 Stellar Luminosity . . . . . . . . . . . . . . . . . . . . . . . . . 4
1.3 Brightness (radiant flux) . . . . . . . . . . . . . . . . . . . . . 4
1.4 Magnitude system . . . . . . . . . . . . . . . . . . . . . . . . . 4
1.5 Limiting magnitude . . . . . . . . . . . . . . . . . . . . . . . . 5
2 Telescopes 7
2.1 Optical Telescopes . . . . . . . . . . . . . . . . . . . . . . . . . 7
2.2 Refracting Telescopes . . . . . . . . . . . . . . . . . . . . . . . 7
2.3 Reflecting Telescopes . . . . . . . . . . . . . . . . . . . . . . . 8
2.4 F-number . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
2.5 Images . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
2.6 Lens and the focal length . . . . . . . . . . . . . . . . . . . . . 10
2.7 Telescope resolving power . . . . . . . . . . . . . . . . . . . . 12
4 Spherical Astronomy 21
v
Contents
5 Celestial Mechanics 33
5.1 Newton’s Law of Gravitation . . . . . . . . . . . . . . . . . . . 33
5.2 Linear Momentum . . . . . . . . . . . . . . . . . . . . . . . . . 34
5.3 Angular Speed . . . . . . . . . . . . . . . . . . . . . . . . . . . 34
5.4 Angular Momentum . . . . . . . . . . . . . . . . . . . . . . . . 35
5.5 Conservation of Angular Momentum . . . . . . . . . . . . . . 35
5.6 Kepler’s laws . . . . . . . . . . . . . . . . . . . . . . . . . . . . 36
5.7 Velocities in Different Orbits in Celestial Mechanics . . . . . . 39
6 Practice problems 43
Appendices 59
A Math Appendix 61
A.1 Circles . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 61
A.2 Parabolas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 62
A.3 Ellipses . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 64
A.4 Hyperbolas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 64
A.5 Angles . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 66
A.6 Trigonometric Identities . . . . . . . . . . . . . . . . . . . . . 68
A.7 Polar Coordinates . . . . . . . . . . . . . . . . . . . . . . . . . 68
References 71
vi
List of Figures
2.1 (a) Objective lens gathers the light and forms an image. (b)
Eyepiece magnifies the image formed by the objective. (c) The
focal length of the objective lens. . . . . . . . . . . . . . . . . . . . 8
2.2 Newtonian reflecting telescope . . . . . . . . . . . . . . . . . . . . 9
2.3 Cassegrain reflecting telescope. . . . . . . . . . . . . . . . . . . . . 9
2.4 Focal plane. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
2.5 Telescope Magnification . . . . . . . . . . . . . . . . . . . . . . . . 11
2.6 Field of view . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
vii
List of Figures
5.2 (a) The law of areas is demonstrated by the equal shaded areas,
which are traversed by a line connecting a planet to the Sun in equal
time intervals. (b) During a time interval t, the line connecting a
planet to the Sun sweeps through an angle theta (θ) while covering
an area A. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 38
viii
List of Tables
ix
PART I
Introduction to CAAO
CHAPTER 1
Basic Concepts
1.1 Parallax
Measuring the intrinsic brightness of stars is linked with determining their
distances. On Earth, the distance to the peak of a remote mountain can be
determined by measuring that peak’s angular position from two observation
points separated by a known baseline distance. Simple trigonometry then
supplies the distance to the peak. Finding the distance even to the nearest
stars requires a longer baseline. As Earth orbits the Sun, two observations of
the same star made 6 months apart employ a baseline equal to the diameter
of Earth’s orbit. These measurements reveal that a nearby star exhibits an
annual back-and-forth change in its position against the stationary background
of much more distant stars. a measurement of the parallax angle p (one-half of
the maximum change in angular position) allows the calculation of the distance
d to the star.
3
1. Basic Concepts
4
1.5. Limiting magnitude
the faintest visible were ranked as 6th magnitude. In other words, the brightest
stars were assigned the smallest number, the faintest the largest number. And 6
divisions were used because of the mysticism about 6, which is the first perfect
number. The brightness ratio of rank first and sixth is 100:
√
K 5 = 100 → K = 5 100 = 2.5118 → bb12 = 2.5118m1 −m2
b1
m1 − m2 = −2.5 log , (1.4)
b2
Where m is the magnitude of the stars. The magnitude system is based on
the comparison; This means that you need to know the magnitude of a certain
star and by comparing its brightness with other stars you can determine the
magnitude of the star.
The absolute magnitude, M , is defined to be the apparent magnitude a
star would have if it were located at 10 pc. Recall that a difference of 5
magnitudes between the apparent magnitudes of two stars corresponds to the
smaller magnitude star being 100 times brighter than the larger-magnitude
star. This allows us to find an equation for the absolute magnitude just like
the apparent magnitude:
L1
M1 − M2 = −2.5 log , (1.5)
L2
Where M is the absolute magnitude and L is the luminosity of the star. To
use this equation, we need to know a specific star’s absolute magnitude and
luminosity to be able to compare it with other stars. We have defined two
different magnitudes, absolute and apparent. We defined absolute magnitude
as the apparent magnitude at a certain distance (10 parsecs), therefore, there
should be a connection between the star distances and their magnitudes. This
relation is called Distance Modulus:
m − M = 5 log d − 5, (1.6)
5
1. Basic Concepts
it is able to gather more light than your eyes can. Therefore, we have the
equation below to determine the limiting magnitude of a telescope:
Dt
me − mt = −5 log , (1.7)
De
Where me is the limiting magnitude of the naked eye (me ≈ 6.5), mt is the
telescope’s limiting magnitude, De is the pupil’s diameter (De ≈ 6 mm) and
Dt is the telescope’s diameter.
6
CHAPTER 2
Telescopes
7
2. Telescopes
Figure 2.1: (a) Objective lens gathers the light and forms an image. (b)
Eyepiece magnifies the image formed by the objective. (c) The focal length of
the objective lens.
8
2.4. F-number
2.4 F-number
Telescopes are often described by both their aperture size and f-number. The f-
number is the ratio of the focal length of the main lens or mirror to the aperture.
These specifications are important because the brightness, size, and clarity of
the image produced by a telescope depend on the aperture and focal length
of its main lens or mirror. For example, a “150-mm (6-inch), f /8 reflector”
means the primary mirror is 150 mm (6-inches) in diameter and has a focal
length of 1200 mm (8 × 150), or 48 inches (8 × 6).
2.5 Images
All stars except our Sun are so far away that they appear as dots of light in
a telescope. The Moon and planets appear as small disks. Image size is
proportional to the focal length of the telescope’s main lens or mirror.
For example, a mirror with a focal length of 2.5 m (100 inches) produces
an image of the Moon that measures about 2.5 cm (1 inch) across. You know
that the 5 m (200-inch), f /3.3 mirror has a focal length of 16.5 m (660-inches),
which is over six times as long. Hence, it produces an image of the Moon that
is about six times as big or 15 cm (6-inches) across.
Lenses and mirrors form real images that are upside down. (A real image
is formed by the actual convergence of light rays.) Since inverted images do
not matter in astronomical work and righting them would require additional
light-absorbing optics, nothing is done to turn images upright in telescopes.
9
2. Telescopes
In the f igure 2.4, two stars separated by a small angle θ on the sky have
images that are separated by a physical distance d on the focal plane. Another
useful parameter, in addition to the focal length, is the scale of the image on
the focal plane, known for historical reasons as the plate scale. Specifically, an
angular distance θ on the celestial sphere is related to a physical distance d on
the image plane by the plate scale s:
10
2.6. Lens and the focal length
where f is the focal ratio, and D is the diameter of the telescope’s aperture.
As an example, the famous “forty-inch” Yerkes Telescope (at Williams Bay,
Wisconsin) has an aperture D = 1.02 m and a focal ratio f = 19. The plate
scale of the Yerkes Telescope is thus:
206.265
s= arcsec/mm = 10.6 arcsec/mm, (2.6)
19 × (1.02)
therefore an image of the full Moon produced by the Yerkes Telescope is 170
mm across, about the size of a salad plate.
The major optical component of a refracting telescope is the primary or
objective lens of focal length fobj . The purpose of the objective lens is to collect
as much light as possible and with the greatest possible resolution, bringing the
light to a focus at the focal plane. A photographic plate or other detector may
be placed at the focal plane to record the image, or the image may be viewed
with an eyepiece, which serves as a magnifying glass. The eyepiece would be
placed at a distance from the focal plane equal to its focal length, feye , causing
the light rays to be refocused at infinity. The figure below shows the path of
rays coming from a point source lying off the optical axis at an angle θ. The
rays ultimately emerge from the eyepiece at an angle ϕ from the optical axis.
The angular magnification produced by this arrangement of lenses can be shown
to be:
fobj
m= . (2.7)
feye
In astronomy, the field of view is the amount of sky you can see, whether
with your unaided vision, binoculars, or a telescope. If you had eyes on all sides
of your head, you would have a 360◦ field of view. (Some insects actually do!)
If you include peripheral vision, your naked eye field of view is nearly 180◦ , but
with varying quality across this field. A telescope will have a much smaller field
of view, but it has significant advantages, such as greater magnification and
light-gathering power.
Field of view (FOV) is the diameter of a region of the sky that you can see
using a specific instrument. The FOV would change with the magnification of
the telescope you are using.
11
2. Telescopes
12
CHAPTER 3
Figure 3.1: Stellar evolution for stars with different initial masses
13
3. Observing the Universe
Formation of Stars
Stars are formed from large clouds of gas and dust, known as nebulae. These
clouds are primarily composed of hydrogen and helium, with smaller amounts
of heavier elements. The process of star formation begins when a region of a
nebula becomes dense enough for gravity to take over. As the gas and dust
collapse under their own weight, the temperature and pressure at the center of
the cloud begin to rise. Eventually, the temperature and pressure become high
enough to trigger nuclear fusion, and a new star is born.
Red Dwarfs
Red dwarfs are the most common type of star in the universe. They are small,
cool stars with masses less than about 0.5 times that of the Sun. Because they
are so small and cool, red dwarfs can burn hydrogen for a very long time, with
some estimated to live up to trillions of years.
White Dwarfs
White dwarfs are the remnants of low-mass stars, such as red dwarfs or main
sequence stars with masses less than about 8 times that of the Sun. When these
stars run out of fuel, they no longer have the outward pressure from nuclear
fusion to balance the inward pull of gravity, and they collapse inward. This
collapse causes the outer layers of the star to be expelled in a planetary nebula,
leaving behind a hot, dense core known as a white dwarf.
Supernovae
Supernovae are some of the most violent events in the universe. They occur
when a massive star runs out of fuel and can no longer generate enough pressure
to resist the inward pull of gravity. The core of the star collapses, creating a
14
3.1. Stellar Evolution
shock wave that causes the outer layers of the star to explode outward in a
brilliant display of light and energy. This explosion can briefly outshine an
entire galaxy and release more energy than the Sun will produce in its entire
lifetime.
There are two types of supernovae: Type I and Type II. Type I supernovae
occur when a white dwarf in a binary star system accretes enough matter from
its companion to reach a critical mass and undergo a runaway fusion reaction.
Type II supernovae occur when a massive star runs out of fuel and collapses,
triggering a shock wave that causes the outer layers of the star to explode.
Figure 3.2: The turbulent disk of gas churning around a black hole takes on a
crazy double-humped appearance
15
3. Observing the Universe
3.2 Galaxies
Galaxies are vast systems of stars, gas, and dust held together by gravity.
They come in a variety of shapes and sizes, and play a crucial role in our
understanding of the universe. We will explore the different types of galaxies
and their properties, as well as some fascinating facts about these cosmic
structures.
Elliptical Galaxies
Elliptical galaxies are the most common type of galaxy in the universe. They
are shaped like ellipsoids and contain mostly old, red stars. Elliptical galaxies
can range in size from small dwarf galaxies to giant ellipticals that are up to 20
times larger than the Milky Way.
The largest known elliptical galaxy, IC 1101, has a diameter of over 6 million
light years and contains trillions of stars.
Spiral Galaxies
Spiral galaxies are characterized by their prominent spiral arms, which are made
up of young, blue stars, gas, and dust. The Milky Way is a spiral galaxy, and
our Sun is located in one of its spiral arms. Spiral galaxies come in a variety of
shapes, from tight and compact to loose and open.
The Andromeda Galaxy, our closest neighboring galaxy, is a spiral galaxy
that is expected to collide with the Milky Way in about 4 billion years.
Irregular Galaxies
Irregular galaxies have no well-defined shape and contain a mix of old and
young stars, gas, and dust. They are often small and found in the vicinity of
larger galaxies, which can disrupt their structure through tidal forces.
16
3.2. Galaxies
Composition
Galaxies are composed of a variety of elements, including hydrogen, helium,
and heavier elements such as carbon and oxygen. The relative abundances of
these elements can provide clues about the formation and evolution of a galaxy.
The stars in the Andromeda Galaxy contain more heavy elements than
those in the Milky Way, suggesting that it has had a more active history of star
formation.
Black Holes
Many galaxies, including the Milky Way, contain supermassive black holes
at their centers. These black holes can have masses ranging from millions to
billions of times that of the Sun and can profoundly influence the evolution of
their host galaxies.
The largest known black hole, TON 618, has a mass of about 66 billion
times that of the Sun and is located in a quasar, an extremely luminous object
powered by material falling onto the black hole.
Galaxies are incredibly diverse and complex systems that continue to
fascinate astronomers and the public alike. By studying their properties and
evolution, we can gain insights into the history and structure of the universe.
Galaxy Evolution
Galaxies are not static objects, but rather are constantly changing over time.
The study of galaxy evolution seeks to understand how galaxies form, grow,
and change over time.
Galaxies are thought to have formed from the gravitational collapse of
primordial gas clouds in the early universe. These clouds were primarily
composed of hydrogen and helium, with small amounts of heavier elements.
As the gas clouds collapsed, the temperature and pressure at the center of the
cloud increased, triggering the formation of the first stars and galaxies.
Over time, galaxies grow by merging with other galaxies. When galaxies
merge, their stars and gas clouds interact gravitationally, causing them to lose
energy and fall towards the center of the newly-formed galaxy. As more and
more galaxies merge together, the resulting galaxy becomes larger and more
massive.
17
3. Observing the Universe
Galaxies also evolve through the process of star formation. As stars form and
die within a galaxy, they release heavy elements into the surrounding gas. Over
time, the concentration of heavy elements within a galaxy increases, leading to
changes in the way the gas clouds behave. This can lead to changes in the rate
of star formation within the galaxy.
Observations of distant galaxies suggest that they were more active in the
past than they are today. Galaxies in the early universe were more likely to be
forming stars at a rapid rate and were often much more irregular in shape than
galaxies we see today. This suggests that galaxies have undergone significant
evolution over time.
Galaxy evolution is a complex and ongoing area of research in astronomy.
By studying the properties of galaxies at different points in cosmic history,
astronomers hope to gain a better understanding of how galaxies form and
evolve over time, and how they contribute to the overall structure and evolution
of the universe.
18
3.4. Hubble’s Law
Where c is the speed of light. Its non-relativistic formula can be written as:
vr
z= (3.4)
c
v = H0 D, (3.5)
The fact that Mercury and Venus are never seen more than 28◦ and 47◦ ,
respectively, east or west of the Sun clearly shows that their orbits are located
19
3. Observing the Universe
inside the orbit of Earth. These planets are referred to as inferior planets,
and their maximum angular separations east or west of the Sun are known as
the greatest eastern elongation and greatest western. elongation, respectively
(see F igure 3.4).
Opposition occurs when the Earth lies between the Sun and the superior
planet. That is, the Sun and planet are 180° apart on the celestial sphere as
seen from the Earth. A conjunction occurs when the Sun lies between the
Earth and the superior planet. That is, the Sun and planet are 0° apart as seen
from the Earth. Quadrature occurs when the Sun and the superior planet
are 90° apart as seen from the Earth. The quadrature can be either eastern,
when the planet appears 90° east of the Sun in the sky, or western when the
planet appears 90° west of the Sun. Inferior conjunction occurs when the
inferior planet lies between the Earth and the Sun. Superior conjunction
occurs when the Sun lies between the Earth and the inferior planet.
The relative orbital motions of Earth and the other planets mean that
the time interval between successive oppositions or conjunctions can differ
significantly from the amount of time necessary to make one complete orbit
relative to the background stars. The former time interval (between oppositions)
is known as the synodic period (S), and the latter time interval (measured
relative to the background stars) is referred to as the sidereal period (P ). The
relationship between the two periods is given by:
1 1 1
= − , (3.6)
S TP TE
where TP is the period of rotation of a random planet, TE is the period
of Earth (1 year), and S is the synodic period (or the period between two
successive same situations of the two planets).
20
CHAPTER 4
Spherical Astronomy
We have seen that the observer who views the heavens at night gets the
impression that they are at the centre of a great hemisphere onto which the
heavenly bodies are projected. The moon, planets, and stars seem to lie on this
celestial hemisphere, their directions defined by the positions they have on its
surface. For many astronomical purposes, the distances are irrelevant so the
radius of the sphere can be chosen at will. The description of the positions of
bodies on it, considering positional changes with time, necessarily involves the
use of special coordinate and timekeeping systems. The relationship between
the positions of bodies requires a knowledge of the geometry of the sphere.
This branch of astronomy, called spherical astronomy, is in one sense the
oldest branch of the subject, its foundations dating back at least 4000 years.
Its subject matter is still essential and never more so than today when the
problem arises of observing or calculating the position of an artificial satellite
or interplanetary probe. We, therefore, begin by considering the geometry of
the sphere.
21
4. Spherical Astronomy
We can draw infinite circles on a sphere, some may have a radius of the
sphere (great circles) and others will have a smaller radius (small circles). In
the figure on the right ANBM, CNDM, and APBQ are all great circles, while
EFG is a small circle.
where all angles must be written in terms of Radians and R is the radius
of the sphere. Just as the formulas of plane trigonometry can be used to
perform calculations in plane geometry, special trigonometrical formulas for use
in spherical geometry can be established. There are many such formulas but
four are more often used than any of the others. They are the relations between
the sides and angles of a spherical triangle and are invaluable in solving the
problems that arise in spherical astronomy.
22
4.2. Position on the Earth’s surface
ABC is a spherical triangle with sides AB, BC, and CA of lengths c, a, and
b, respectively, and with angles ∠CAB,∠ABC, and ∠BCA hereafter referred
to as angles A, B and C respectively. The four formulas are:
Sine formula:
sin a sin b sin c
= = (4.2)
sin A sin B sin C
Cosine formula:
cos a = cos b cos c + sin b sin c cos A
cos b = cos a cos c + sin a sin c cos B (4.3)
cos c = cos a cos b + sin a sin b cos C
Four-parts formula:
23
4. Spherical Astronomy
The latitude of a point on the Earth’s surface is its angular distance from
the equator, measured along a great circle perpendicular to the Earth’s equator.
Latitude is measured in degrees, arc-minutes, and arc-seconds, as is longitude.
Thus, the use of latitude and longitude does not require knowing the size of
the Earth in kilometers or any other unit of length. The longitude may be
expressed in angular measure or in time units related to each other by the table
on the right.
24
4.3. The horizontal (alt-azimuth) system
altitude of 90◦ . Points on the horizon circle are at an altitude of 0◦ . The nadir
is at an altitude of −90◦ , but in practice, negative altitudes are seldom used,
since they represent objects that are hidden by the Earth. The longitude-like
coordinate in the horizon coordinate system is called the azimuth
The azimuth of this star is the red angle shown in the figure: A = 360◦ − N\OM .
Azimuth is usually expressed from North to East. But if the star is located in
the western hemisphere (like the star in the figure), we can express the azimuth
from North to West: A = N \ OM W .
For any point on the celestial sphere, half a great circle can be drawn from
the zenith, through the point in question, to the nadir. The half-circle that
runs through the north point on the horizon circle acts as the ıprimemeridianȷ
in the horizon coordinate system. The azimuth is measured in degrees running
from north to east. An object due north of an observer has an azimuth of 0°,
an object due east has an azimuth of 90°, and so forth. If you know the altitude
and azimuth of any object in your horizon coordinate system, you know where
to point your telescope to see it.
25
4. Spherical Astronomy
26
4.4. The equatorial system
sphere is called the celestial equator. The celestial equator passes through
the zenith for an observer on the Earth’s equator.
On the Earth’s surface, a point’s latitude is its angular distance north or
south of the equator. Similarly, on the celestial sphere, a point’s declination (δ)
is its angular distance north or south of the celestial equator. For points north
of the celestial equator, the declination is positive (0◦ < δ ≤ 90◦ ), and for points
south of the celestial equator, the declination is negative (−90◦ ≤ δ < 0◦ ).
Right ascension α is analogous to longitude and is measured eastward
along the celestial equator from the vernal equinox (γ) to its intersection
with the object’s hour circle (the great circle passing through the object being
considered and through the north celestial pole). Right ascension is traditionally
measured in hours, minutes, and seconds. The coordinates of the right ascension
and declination are also indicated in the figure below. Since the equatorial
coordinate system is based on the celestial equator and the vernal equinox,
changes in the latitude and longitude of the observer do not affect the values of
right ascension and declination. Values of and are similarly unaffected by the
annual motion of Earth around the Sun.
27
4. Spherical Astronomy
argument, E is the east point. Any great semicircle through P and Q is called
a meridian. The meridian through the celestial object X is the great semicircle
P XBQ cutting the celestial equator in B.
28
4.5. The ecliptic system
through the star from 0h to 24h or from 0◦ to 360◦ . Consequently, the hour
angle increases by 24h each sidereal day for a star. Having both coordinates on
the sphere, using Zenith, the North celestial pole, and the star (three points) we
are able to create a spherical triangle (figure below). We need to use spherical
trigonometry to solve any spherical triangle.
29
4. Spherical Astronomy
Figure 4.13: The Celestial sphere used for coordinate system conversion
30
4.5. The ecliptic system
Let’s assume the equatorial coordinates of a star are known, and we want
to determine its ecliptic coordinates. This means α and δ are given. Using the
spherical triangle above, we can use cosine law:
cos (90 − β) = cos ϵ cos (90 − δ) + sin ϵ sin (90 − δ) cos (90 + α)
(4.8)
−→ sin β = cos ϵ sin δ − sin ϵ cos δ cos α
cos (90 − δ) = cos ϵ cos (90 − β) + sin ϵ sin (90 − β) cos (90 − λ)
−→ sin δ = cos ϵ sin β − sin ϵ cos β cos λ (4.9)
−→ sin λ = sin sin
δ−cos ϵ sin β
ϵ cos β
31
CHAPTER 5
Celestial Mechanics
Gravitational force is one of the four fundamental forces of nature, along with
electromagnetic force, weak nuclear force, and strong nuclear force. It is the
force that causes objects with mass to be attracted to one another. In this
article, we will explore the concept of gravitational force and its mathematical
description.
The English physicist Sir Isaac Newton was the first to describe the nature
of gravitational force. He formulated his law of gravitation in 1687, which
states that the force of attraction between two objects with masses m1 and m2 ,
separated by a distance r, is given by:
m1 m2
FG = G , (5.1)
r2
33
5. Celestial Mechanics
Gravitational Field
The gravitational force can also be described in terms of a gravitational field.
A gravitational field is a region of space where an object with mass experiences
a force due to the presence of another object with mass. The gravitational field
strength at a point in space is defined as the force per unit mass experienced
by a small test mass placed at that point.
The gravitational field strength at a distance r from a point mass M is given
by:
GM
g= , (5.2)
r2
where G is the gravitational constant. The gravitational field strength is a
vector quantity, pointing towards the point mass M .
GM m
U =− , (5.3)
r
The negative sign indicates that the gravitational force is attractive, and the
potential energy is lower at closer distances.
p = mv, (5.4)
where p is the linear momentum, m is the mass of the object, and v is its
velocity.
Linear momentum is conserved in an isolated system, meaning that the total
linear momentum of a system remains constant if no external forces act upon
it. This principle is known as the law of conservation of linear momentum.
34
5.4. Angular Momentum
θ
ω= , (5.5)
t
where ω is the angular speed, θ is the angular displacement of the object, and t
is the time taken for the object to complete the rotation.
Angular speed is measured in radians per second (rad/s). It is important to
note that angular speed is not the same as linear speed, which is the distance
traveled per unit time.
L = Iω, (5.6)
where L is the angular momentum, I is the moment of inertia of the object,
and ω is its angular velocity.
Angular momentum can also be expressed as the product of the mass of the
object, its tangential velocity, and the distance from the axis of rotation:
L = mrv, (5.7)
where m is the mass of the object, v is its tangential velocity, and r is the
distance from the axis of rotation.
Angular momentum is also a vector quantity, meaning it has both magnitude
and direction. Its direction is perpendicular to the plane of rotation. The
moment of inertia is a measure of an object’s resistance to rotational motion
and depends on both the mass and the distribution of mass relative to the axis
of rotation.
Like linear momentum, angular momentum is conserved in an isolated system.
This principle is known as the law of conservation of angular momentum. The
law states that if no external torques act upon an isolated system, the total
angular momentum of the system remains constant. Mathematically, this can
be expressed as:
dL
= τnet , (5.8)
dt
where dL/dt is the rate of change of angular momentum and τnet is the net
external torque acting on the system. If there is no net external torque, then
dL/dt is zero and the angular momentum of the system is conserved.
35
5. Celestial Mechanics
36
5.6. Kepler’s laws
Figure 5.1: A planet of mass m moving in an elliptical orbit around the Sun
with mass M
2. The Law of Areas: dictates that, during equal intervals of time, the
imaginary line that connects a planet to its central star will cover equal areas.
Figure 4.2 serves to visually demonstrate this concept, and implies that an
orbiting object will move with greater velocity when it is nearer to the central
body than when it is further away. It can be proven that the Law of Areas is
in fact equivalent to the Law of Conservation of Angular Momentum.
If we examine the small area increment A that is traversed during a time
interval t, as illustrated in Figure 4.2, we can see that the area of the triangular
wedge is roughly equal to half of its base, r∆θ multiplied by its height r. We
can then calculate the rate at which this area is swept out:
dA ∆A 1 ∆θ 1
= lim = lim r2 = r2 ω (5.10)
dt ∆t→0 ∆t ∆t→0 2 ∆t 2
If we make the assumption that the more massive body M can be considered
stationary, then the angular momentum of the orbiting body m can be described
relative to the origin at the central body as:
Lz = Iω = mr2 ω (5.11)
Thus:
dA Lz
= (5.12)
dt 2m
If the M and m system is isolated and there is no external torque acting
on it, then the angular momentum Lz remains constant. This means that the
derivative of the area A with respect to time t is also constant, as stated in the
equation. Consequently, during each interval of time dt, the line connecting m
and M sweeps out an equal area dA, which confirms Kepler’s second law. The
increase in speed of a comet as it passes close to the Sun is an example of this
effect and is directly related to the law of conservation of angular momentum.
37
5. Celestial Mechanics
Figure 5.2: (a) The law of areas is demonstrated by the equal shaded areas,
which are traversed by a line connecting a planet to the Sun in equal time
intervals. (b) During a time interval t, the line connecting a planet to the Sun
sweeps through an angle theta (θ) while covering an area A.
GM m v2
= m . (5.13)
r2 r
Then replacing the speed v with 4πr/T , where T is the rotational period (the
time for a full orbit), we obtain:
4π 2 3
T2 = ( )r . (5.14)
GM
The same outcome can be achieved for orbits that are elliptical, where the
radius r is substituted with the semi-major axis a. The constant ratio between
38
5.7. Velocities in Different Orbits in Celestial Mechanics
Circular Orbit
A circular orbit is a special case of an elliptical orbit where the semi-major
axis a is equal to the radius r. Kepler’s third law states that the square of the
orbital period T is proportional to the cube of the semi-major axis a:
4π 2 3
T2 = a . (5.15)
GM
The velocity of a circular orbit can be derived by equating the centripetal force
Fc with the gravitational force Fg :
Fc = Fg , (5.16)
which can be written as:
mv 2 GM m
= . (5.17)
r r2
Simplifying this expression, we obtain the velocity of a circular orbit:
r
GM
v= . (5.18)
r
39
5. Celestial Mechanics
Elliptical Orbit
An elliptical orbit is described by the semi-major axis a and the eccentricity e,
where e is the ratio of the distance between the foci of the ellipse to the length
of the major axis. Kepler’s second law states that the area swept out by the
radius vector in a given time is constant, which implies that the speed of the
orbiting body varies along the orbit.
We can derive the velocity of an elliptical orbit by using the conservation of
angular momentum, which states that the product of the mass m, the velocity
v, and the distance r from the center of mass to the orbiting body is constant:
mvr = h, (5.19)
where h is the specific angular momentum. We can express h in terms of the
semi-major axis a and the eccentricity e:
h= GM a(1 − e2 ). (5.20)
p
Using Kepler’s second law, we can express the speed v at any point in the orbit
as:
h 1
v= , (5.21)
r 1 + e cos θ
where θ is the true anomaly, which is the angle between the position vector of
the orbiting body and the pericenter of the orbit.
We can simplify this expression by expressing the distance r in terms of the
semi-major axis a and the eccentricity e:
a(1 − e2 )
r= . (5.22)
1 + e cos θ
Substituting this expression into the equation for the velocity, we obtain the
velocity of an elliptical orbit:
s r
2 1 2GM
GM
v = GM − = − . (5.23)
r a r a
Parabolic Orbit
A parabolic orbit is an orbit in which the distance between the two bodies
approaches infinity. The velocity of a parabolic orbit can be derived using the
concept of specific energy, which is the sum of the kinetic and potential energy
per unit mass of the orbiting body. For a parabolic orbit, the specific energy is
zero, which means that the kinetic energy is equal in magnitude to the potential
energy. Thus, the total energy is also zero.
Using the conservation of energy, we can equate the kinetic energy to the
negative potential energy:
1 GM m
mv 2 = − . (5.24)
2 r
Solving for the velocity, we obtain the velocity of a parabolic orbit:
40
5.7. Velocities in Different Orbits in Celestial Mechanics
r
2GM
v= . (5.25)
r
In conclusion, the velocity of conic section orbits can be derived using the
principles of classical mechanics and the laws of gravity. Circular orbits have
a constant velocity, while elliptical orbits have varying speeds along the orbit.
Parabolic and hyperbolic orbits have specific energies that result in unique
velocities. By understanding the velocity of conic section orbits, we can better
understand the motion of celestial objects and their interactions with each
other.
41
CHAPTER 6
Practice problems
1. The luminosity of the Sun is L = 3.85 × 1026 W . If the distance of the Earth
from Sun is equal to 1 AU . Determine what is the flux that Earth receives from
Sun.(1AU = 1.496 × 1011 m)
Solution The flux (or intensity) that the Earth receives from the Sun can
be calculated using the inverse square law of radiation, which states that the
intensity of radiation decreases with the square of the distance from the source.
The flux at a distance r from a source with luminosity L is given by:
L
F = ,
4πr2
where F is the flux received per unit area.
Substituting the given values, we get:
3.85 × 1026 W
F =
4π(1.496 × 1011 m)2
Simplifying this expression, we get:
2
F = 1361 W/m
Therefore, the flux that the Earth receives from the Sun is 1361 W/m2 . This
quantity is known as the solar constant and is an important parameter in
climate science and solar energy applications.
2. Determine the brightness of a star with the magnitude of 2. You can use
Sun as the known star to compare. (The apparent magnitude of our Sun is
-26.8)
Solution The magnitude system used to measure the brightness of stars is
logarithmic, meaning that a difference of 1 magnitude corresponds to a factor
of 2.512 in brightness.
the apparent magnitude of the Sun (m⊙ = −26.8) as a reference, we can
calculate the ratio of the brightness of the star to the brightness of the Sun:
L⋆
= 2.512m⊙ −m⋆ ,
L⊙
where L⋆ and L⊙ are the luminosities of the star and the Sun, respectively, and
m⋆ is the apparent magnitude of the star.
Plugging in the values, we get:
43
6. Practice problems
L⋆
= 2.512−26.8−2 ≈ 3.98 × 10−10 .
L⊙
Therefore, the star is about 3.98 × 10−10 times as bright as the Sun.
To convert this ratio to a measure of brightness (flux), we can use the formula:
F⋆ L⋆
=
F⊙ 4πd2⋆ L⊙
where F⋆ and F⊙ are the fluxes (energy per unit area per unit time) received
from the star and the Sun, respectively, and d⋆ is the distance to the star.
that the star is at the same distance from us as the Sun (d⋆ = 1 AU
= 1.496 × 1011 m), we can plug in the values and solve for F⋆ :
F⋆ 3.98 × 10−10
= ≈ 4.38 × 10−5
F⊙ 4π(1.496 × 1011 )2
Therefore, the star has a flux (brightness) of about 4.38 × 10−5 times that of
the Sun at Earth’s distance.
8 × 25
mt = 6.5 + 5 log = 14.1
6
Therefore, the limiting magnitude of an 8-inch telescope is approximately
mlimit = 14.1. This means that the telescope can detect stars with an apparent
magnitude of 14.1 or brighter.
44
For the Ca II H line, we have:
v = H0 d,
where v is the recessional velocity, d is the distance, and H0 is the Hubble
constant. The recessional velocity can be calculated from the redshift using the
formula:
v = cz,
where c is the speed of light. Therefore, we have:
v cz
d= =
H0 H0
The current value of the Hubble constant is a matter of debate and has been
measured to be around 73; km/s/Mpc by some recent studies. Using this value,
we can calculate the distance to the galaxy as:
45
6. Practice problems
conjunction of Mars will occur after half of the synodic period, which is 390
days from the date of the observation.
(b) The synodic period of Jupiter is the time it takes for Jupiter to return to the
same relative position with respect to Earth and the Sun. Since we know that
Jupiter is currently in western quadrature with Earth, we can use the synodic
period to determine when the next opposition will occur.
synodic period of Jupiter is approximately 398.9 Earth days. To find the
time between the current western quadrature and the next opposition, we need
to find the difference between the current phase angle (i.e., the angle between
the Sun-Jupiter line and the Sun-Earth line) and 180 degrees, which is the angle
between the Sun-Jupiter line and the Earth-Jupiter line at opposition.
western quadrature, the elongation angle is 90 degrees. Therefore, the time
between western quadrature and opposition is:
cos−1 ( 5.2
1 AU
AU )
∆t = 360 × 398.9 days = 87.44 days
6. The time interval between two successive oppositions of Mars is 779.9 days.
Calculate the semi-major axis of Mars’ orbit.
Solution First, we can use the given time interval between two successive
oppositions of Mars to find the synodic period Ps :
1
Ps = 1 1
PM ars − PEarth
where PM ars is the orbital period of Mars and PEarth is the orbital period of
Earth.
The orbital period of Earth is approximately 365.25 days, and we are given
that the time interval between two successive oppositions of Mars is 779.9 days.
Therefore, we can solve for the synodic period:
1
Ps = 1 1 = 779.9 days
PM ars − 365.25
1/3 !1/3
2
GP 2 6.674 × 10−11 m3 /kg/s × (687 days × 24 hours/day × 3600 s/hour)2
a= =
4π 2 4π 2
a ≈ 2.28 × 1011 m
Therefore, the semi-major axis of Mars’ orbit is approximately 2.28 × 1011
meters, which is equal to 1.52 AU
46
7. The quasar SDSS 1030 + 0524 produces a hydrogen emission line of
wavelength λrest = 121.6 nm. On Earth, this emission line is observed to have
a wavelength of λobs = 885.2 nm:
(a) What is the redshift of this quasar?
(b) Determine the radial velocity of the quasar.
(c) Determine the distance of this quasar.
Solution (a) Using the formula z = λobsλ−λ rest
rest
, we have:
z= 885.2 nm−121.6;nm
121.6 nm ≈ 6.266
Therefore, the redshift of this quasar is z = 6.266.
(b) The radial velocity of the quasar can be calculated using the formula v = cz,
where c is the speed of light. Substituting the values, we get:
v = (2.998 × 108 m/s) × 6.266 ≈ 1.88 × 109 m/s
Therefore, the radial velocity of the quasar is v ≈ 1.88 × 109 m/s.
(c) We can use the Hubble’s law to estimate the distance of the quasar:
H0 = 73.3 (km/s)/M pc
v = H0 × d
where d is the distance to the quasar. Converting v to meters per second and
H0 to meters per second per meter, we get:
v = 1.88 × 109 m/s
H0 = 73.3 (km/s)/M pc ≈ 2.39 × 10−18 (m/s)/m
Substituting the values, we get:
1.88×109 m/s
d = Hv0 ≈ 2.39×10 −18 (m/s)/m ≈ 7.87 × 10
27
m
Therefore, the distance to the quasar is d ≈ 7.87 × 1027 m.
8. Solve completely the spherical triangle ABC, and find the area of the
triangles: (assume R=1)
(a) a = 34◦ 46′ , b = 27◦ 22′ , C = 72◦ 31′
(b) b = 98◦ 18′ , C = 24◦ 49′ , A = 68◦ 36′
(c) a = 14◦ 03′ , b = 53◦ 32′ , c = 124◦ 14′
(d) A = 23◦ 32′ , B = 102◦ 38′ , C = 34◦ 44′
Solution To solve a spherical triangle completely, we need to find all three
sides and all three angles. We will use the following formulas for the spherical
law of cosines and sines:
Law of cosines:
cos a = cos b cos c + sin b sin c cos A
cos b = cos a cos c + sin a sin c cos B
cos c = cos a cos b + sin a sin b cos C
Law of sines:
sin a sin b sin c
= =
sin A sin B sin C
The area of a spherical triangle can be found using the following formula:
Area = R2 · (A + B + C − π),
47
6. Practice problems
Area = R2 · (A + B + C − π) ≈ 0.0042
(b) b = 98◦ 18′ ; , C = 24◦ 49′ ; , A = 68◦ 36′ We have one side and two angles
opposite to it. Using the law of sines, we can find the other two sides, and then
use the law of cosines to find the remaining angle.
Let a be the side opposite to angle A, and c be the side opposite to angle C.
Then, from the law of sines:
sin a sin c sin A sin C
= ⇒ sin a = ≈ 0.9907
sin A sin C sin b
sin c sin b sin C sin b
= ⇒ sin c = ≈ 0.6478
sin C sin B sin B
Now, using the law of cosines to find the third side:
Area = R2 · (A + B + C − π) ≈ 0.0806
(c) a = 14◦ 03′ , b = 53◦ 32′ , c = 124◦ 14′ We have all three sides. Using the
law of cosines, we can find all three angles.
48
cos c − cos a cos b
cos C = ≈ −0.8295
sin a sin b
Since all three cosines are negative, we have an obtuse triangle. Therefore, we
have to use the supplementary formula for each angle:
Area = R2 · (A + B + C − π) ≈ 0.0245
(d) A = 23◦ 32′ , B = 102◦ 38′ , C = 34◦ 44′ We have all three angles. Using
the law of sines, we can find all three sides.
Let a be the side opposite to angle A, b be the side opposite to angle B, and c
be the side opposite to angle C. Then, from the law of sines:
a b c
= =
sin A sin B sin C
Thus, we have:
sin A
a= c ≈ 0.4501
sin C
sin B
b= c ≈ 1.2467
sin C
sin C
c= c=R≈1
sin C
To find the area, we can use the formula:
Area = R2 · (α + β + γ − π) = R2 · (A + B + C − π) ≈ 0.2796
9. Two cities A and B on the same parallel of latitude ϕ = 43◦ 39′ N are
127◦ 22′ apart in longitude. Calculate in kilometers:
(a)their distance apart along the parallel. (On the small circle between the
cities with same latitude)
(b)the great circle distance AB.
(c)Determine the highest latitude of the great circle between two cities.
Solution We can use the following formulas for the calculations:
For a sphere with radius R, the length L of a small circle with radius r and
central angle θ is L = Rθ. The great circle distance d between two points on a
sphere with radius R and colatitudes ϕ1 and ϕ2 and difference in longitudes
∆λ is d = R cos−1 (sin ϕ1 sin ϕ2 + cos ϕ1 cos ϕ2 cos ∆λ).
(a) To calculate the distance between two points along a parallel of latitude,
we can assume that the Earth is a sphere with radius R = 6, 371 km and use
the formula L = Rθ, where θ is the central angle between the two points
on the sphere. The distance between the two cities along the parallel is
49
6. Practice problems
simply the arc length of the parallel between the two longitudes, which is
◦
22′
θ = 127
360◦ · 2πR cos ϕ ≈ 14, 105 km.
(b) To calculate the great circle distance between the two cities, we can use
the formula d = R cos−1 (sin ϕ1 sin ϕ2 + cos ϕ1 cos ϕ2 cos ∆λ), where ϕ1 = ϕ2 =
43◦ 39′ N and ∆λ = 127◦ 22′ . Then, we have:
d = 6, 371 km cos−1 (sin 43◦ 39′ N sin 43◦ 39′ N+cos 43◦ 39′ N cos 43◦ 39′ N cos 127◦ 22′ )
d ≈ 9, 325 km
Therefore, the great circle distance between the two cities is approximately
9, 325 km.
12. Show that the point on the horizon at which a star rises is:
Solution To solve this problem, we will need to use some basic trigonometric
identities and the geometry of the celestial sphere. Let’s begin by defining some
terms:
ϕ is the observer’s latitude, measured in degrees or radians. δ is the star’s
declination, measured in degrees or radians. H is the star’s hour angle, measured
in hours or radians. h is the altitude of the star, measured in degrees or radians.
Using these definitions, we can write the following equations:
sin ϕ sin δ
cos H = −
cos ϕ cos δ
the identity sec θ = cos θ ,
1
we can write this as:
50
sin ϕ sin δ sin ϕ sin δ
cos H = − =− · = − tan ϕ tan δ
cos ϕ cos δ cos ϕ cos δ
the inverse sine of both sides, we get:
sin ϕ sin δ
sin−1 (− tan ϕ tan δ) = sin−1 − = sin−1 (sec ϕ sin δ)
cos ϕ cos δ
is the desired result, which shows that the point on the horizon at which a
star rises is given by sin−1 (sec ϕ sin δ).
13. We have the coordinates of Vega δV = 38◦ 47′ , αV = 18h 36m . A person in
Toronto (ϕT oronto = 43.65◦ N ) is observing this star:
(a) Determine the hour angle of Vega when it rise/set.
(b) What is the Azimuth of rise and set of Vega in Toronto’s horizon?
(c) Determine its maximum altitude in Toronto.
(d) Determine the total time Vega is above horizon.
(e) On which date does Vega rise at the same time as the Sun in Toronto?
14. The parallax angle of a star is measured to be 0.4 arcseconds. What is the
distance to the star in parsecs? Assume that the distance to the star is much
greater than the radius of the Earth’s orbit.
Solution The parallax angle θ is related to the distance d to the star by the
formula:
θ = 1 AU
d
where 1 AU is the mean distance between the Earth and the Sun, which is
approximately 1.496 × 1011 meters.
Converting the parallax angle to radians:
θ = 0.4 arcseconds × 180×3600
π
radians per arcsecond = 1.184 × 10−6 radians
Substituting into the formula, we get:
11
d = 1 AU
θ = 1.496×10
1.184×10−6 = 1.262 parsecs
m
Note that the assumption that the distance to the star is much greater than
the radius of the Earth’s orbit is necessary to ensure that the parallax angle is
small enough to measure accurately.
51
6. Practice problems
15. Star A has an apparent magnitude of 3.5 and a parallax of 0.05 arcseconds.
Star B has an apparent magnitude of 2.0 and a parallax of 0.02 arcseconds.
Which star is closer to Earth, and by how much in parsecs?
Solution We can use the formula for converting parallax to distance:
d = p1 ,
where d is the distance in parsecs and p is the parallax in arcseconds.
For Star A, p = 0.05 arcseconds, so dA = 0.05 1
= 20 parsecs.
For Star B, p = 0.02 arcseconds, so dB = 0.02 = 50 parsecs.
1
Since Star A is farther away from Earth than Star B, it must be dimmer due
to the inverse square law of
light.
We can use the formula:
m2 − m1 = −2.5 log10 F1 ,
F2
where m1 and m2 are the apparent magnitudes of Star A and Star B, respectively,
and F1 and F2 are their corresponding fluxes. We can assume that their intrinsic
brightnesses are equal, so their flux ratio is simply the ratio of their distances
squared:
2
F2
F1 = d1
d2 .
Plugging in the values, weget:
d2
m2 − m1 = −2.5 log10 d12 ,
22
m2 − m1 = −2.5 log10 20 502 = −1.95.
Therefore, Star A is dimmer than Star B by 1.95 magnitudes.
(c) To find the radius R of the star in units of the Sun’s radius, we can use the
formula for luminosity:
L = 4πR2 σT 4
Rearranging
q for R, we have:
R= L
4πσT 4
52
Substituting
q in 2the values we have:
10 L⊙
R = 4πσ(5778 K)4 = 6.48R⊙
Therefore, the star has a radius of 6.48 times the radius of the Sun.
18. A telescope has a focal length of 1000mm and an eyepiece with a focal length
of 20mm. What is the magnification of the telescope? If the telescope is used
to observe an object at a distance of 2000m, what is the angular magnification?
Solution The magnification of the telescope is given by:
fobj
M = feyepiece
where fobj is the focal length of the objective lens/mirror, and feyepiece is the
focal length of the eyepiece.
Substituting the given values, we get:
M = 100020 = 50
Therefore, the magnification of the telescope is 50.
The angular magnification is given by:
θobj
Mang = M × θeyepiece
where θobj is the angular size of the object as seen through the objective
lens/mirror, and θeyepiece is the angular size of the image formed by the eyepiece.
The angular size of an object is given by:
θ= D d
where d is the actual size of the object, and D is the distance from the observer
to the object.
Substituting the given values, we get:
θobj = 2000
0.01
= 5 × 10−6 radians
The angular size of the image formed by the eyepiece is given by:
deyepiece
θeyepiece = feyepiece
where deyepiece is the distance between the eyepiece and the image formed by
the objective lens/mirror.
Assuming that the telescope is in normal adjustment, we have:
deyepiece ≈ fobj
Substituting the given values, we get:
53
6. Practice problems
−3
θeyepiece = 25×10
20 = 1.25 × 10−3 radians
Substituting these values into the formula for angular magnification, we get:
5×10−6
Mang = 50 × 1.25×10 −3 ≈ 0.2
19. A telescope has a focal length of 2000 mm and a plate scale of 1 arcsecond
per pixel. If we want to image the planet Jupiter, which has an angular size
of approximately 50 arcseconds, what should be the size of the imaging sensor
to capture the entire planet? What is the approximate pixel resolution of the
image? Assume that the telescope is diffraction limited.
Solution The plate scale of a telescope is defined as the angular size of one
pixel in arcseconds per millimeter. In other words:
Plate scale = focal length×pixel
206,265
size
where the focal length is in millimeters and the pixel size is in millimeters.
We can rearrange this equation to solve for the pixel size:
Pixel size = focal length×plate
206,265
scale
Plugging in the given values, we get:
Pixel size = 2000 mm×1 arcsecond/pixel ≈ 0.103 mm/pixel
206,265
To capture the entire planet Jupiter, we need an imaging sensor that is at least
50 arcseconds across. Since we know the plate scale and pixel size, we can
calculate the required number of pixels:
Number of pixels = PlateAngular size
scale×Pixel size
Plugging in the values, we get:
Number of pixels = 1 arcsecond/pixel×0.103
50 arcseconds
mm/pixel ≈ 485 pixels
Note that this is only an approximate value, and the actual number of pixels
needed may be higher due to other factors such as image cropping and
interpolation.
Finally, we can calculate the approximate pixel resolution of the image by
dividing the angular size of one pixel by the distance to Jupiter:
Pixel resolution = Distance
Plate scale
to Jupiter
Since Jupiter is about 5.2 astronomical units (AU) from Earth on average, we
have:
Pixel resolution = 1 arcsecond/pixel
5.2 AU ≈ 0.01 arcseconds/pixel
This means that features on Jupiter that are at least 0.01 arcseconds in size
can be resolved by the telescope.
20. A telescope with a focal length of 1000 mm is used to observe Jupiter. The
angular diameter of Jupiter is 40 arcseconds.
a) What is the plate scale of the telescope in arcseconds/mm?
b) If the telescope is equipped with an eyepiece with a magnification of 100x,
what is the apparent diameter of Jupiter in the eyepiece?
c) If Jupiter has an average distance from the Sun of 5.2 astronomical units
(AU), what is its actual diameter in kilometers?
Solution a) The plate scale of a telescope is the angular size of an object in
the sky per unit distance on the detector or the focal plane. It is given by the
formula:
54
Angular Diameter of Object
Plate Scale =
Focal Length of Telescope
Substituting the given values, we get:
40′′
Plate Scale = = 0.04′′ /mm
1000 mm
b) The apparent diameter of Jupiter in the eyepiece is given by the formula:
4000′′
Actual Diameter = ≈ 139, 822 km
0.04′′ /mm × 1000 mm/km × 778.72 million km
Therefore, the actual diameter of Jupiter is approximately 139, 822 km.
21. A telescope with a focal length of 1200 mm and a plate scale of 0.5
arcseconds/mm is used to observe the full moon. The field of view of the
telescope is 30 arcminutes. The magnitude of the moon is -12.7.
a) What is the apparent size of the full moon in the telescope’s field of view?
b) How many pixels will the full moon occupy in an image taken with a
camera that has a pixel size of 5 microns?
55
6. Practice problems
22. A small asteroid is in an elliptical orbit around the Sun. The asteroid has
a semimajor axis of 2.5 AU and an eccentricity of 0.3.
a) What is the period of the asteroid’s orbit?
b) What is the asteroid’s speed when it is closest to the Sun (at perihelion)?
c) What is the asteroid’s speed when it is farthest from the Sun (at aphelion)?
d) What is the asteroid’s orbital energy?
Solution a) The period T of an object in an elliptical orbit can be calculated
using Kepler’s third law:
4π 2 3
T2 = a
GM
where G is the gravitational constant, M is the mass of the central object (in
this case, the Sun), and a is the semimajor axis of the orbit. Substituting the
given values, we get:
4π 2
T2 = (2.5 AU)3
(6.674 × 10−11 m3 /kg s2 )(1.989 × 1030 kg)
Converting AU to meters, and taking the square root, we get:
T = 4.33 years
Therefore, the period of the asteroid’s orbit is approximately 4.33 years.
) At perihelion, the asteroid is closest to the Sun, and therefore its speed is
at a maximum. The speed of an object in an elliptical orbit can be calculated
using the following equation:
56
2 1
v 2 = GM −
r a
where r is the distance between the asteroid and the Sun, and a is the semimajor
axis of the orbit. At perihelion, the distance r is equal to (1 − e)a, where e is
the eccentricity of the orbit. Substituting the given values, we get:
s
2 1
vperi = (6.674 × 10−11 m3 /kg s2 )(1.989 × 1030 kg) −
(1 − 0.3)(2.5 AU) 2.5 AU
s
2 1
vapo = (6.674 × 10−11 m3 /kg s2 )(1.989 × 1030 kg) −
(1 + 0.3)(2.5 AU) 2.5 AU
Simplifying, we get:
E = −6.35 × 1022 J
Therefore, the asteroid’s orbital energy is approximately −6.35 × 1022 J.
57
Appendices
APPENDIX A
Math Appendix
x2 y2
x2 + y 2 = 1 y = x2 + 1 + = 1 x2 − y 2 = 1
9 4
which represents a circle, a parabola, an ellipse, and a hyperbola, respectively.
The graph of such an equation in x and y is the set of all points (x, y) that
satisfy the equation; it gives a visual representation of the equation. Conversely,
given a curve in the xy-plane, we may have to find an equation that represents
it, that is, an equation satisfied by the coordinates of the points on the curve
and by no other point. This is the other half of the basic principle of analytic
geometry formulated by Descartes and Fermat. The idea is that if an algebraic
equation can represent a geometric curve, then the rules of algebra can be used
to analyze the geometric problem.
A.1 Circles
As an example of this type of problem, let’s find an equation of the circle with
radius r and center (h, k). By definition, the circle is the set of all points P (x, y)
whose distance from the center C(h, k) is r. (See Figure A1.)
Thus P are on the circle if and only if |P C| = r. From the distance formula,
we have:
(x − h)2 + (y − k)2 = r,
p
61
A. Math Appendix
(x − h)2 + (y − k)2 = r2 ,
This is the desired equation.
Equation of a circle: An equation of the circle with center (h, k) and radius r is
(x − h)2 + (y − k)2 = r2 ,
In particular, if the center is the origin (0, 0), the equation is
x2 + y 2 = r 2 .
A.2 Parabolas
We regard a parabola as a graph of an equation of the form y = ax2 + bx + c.
Let’s draw the graph of the parabola y = x2 . We set up a table of values, plot
points, and join them by a smooth curve to obtain the graph in Figure A2.
x y = x2
0 0
± 21 1
4
±1 1
±2 4
±3 9
Figure A.3 shows the graphs of several parabolas with equations of the form
y = ax2 for various values of the number a. In each case the vertex, the point
where the parabola changes direction, is the origin. We see that the parabola
y = ax2 opens upward if a > 0 and downward if a < 0 (as in Figure A.3).
Notice that if (x, y) satisfies y = ax2 , then so does (−x, y). This corresponds
to the geometric fact that if the right half of the graph is reflected about the
y-axis, then the left half of the graph is obtained. We say that the graph is
symmetric with respect to the y-axis.
The graph of an equation is symmetric with respect to the y-axis if the equation
is unchanged when x is replaced by −x.
62
A.2. Parabolas
63
A. Math Appendix
A.3 Ellipses
The curve with equation:
x2 y2
+ =1
a2 b2
where a and b are positive numbers, is called an ellipse in standard position.
Equation above is unchanged if x is replaced by −x or y is replaced by −y, so
the ellipse is symmetric with respect to both axes. As a further aid to sketching
the ellipse, we find its intercepts.
The x-intercepts of a graph are the x-coordinates of the points where the
graph intersects the x-axis. They are found by setting y = 0 in the equation of
the graph.
The y-intercepts are the y-coordinates of the points where the graph intersects
the y-axis. They are found by setting x = 0 in its equation.
If we set y = 0 in the equation of ellipse, we get x2 = a2 and so the x-intercepts
are ±a. Setting x = 0, we get y 2 = b2 , so the y-intercepts are ±b. Using this
information, together with symmetry, we sketch the ellipse in Figure A.6 . If
a = b, the ellipse is a circle with radius a.
x2 y2
Figure A.6: Ellipse with equation a2 + b2 =1
A.4 Hyperbolas
The curve with equation
x2 y2
− =1
a2 b2
is called a hyperbola in standard position. Again, Equation above is unchanged
when x is replaced by −x or y is replaced by −y, so the hyperbola is symmetric
with respect to both axes. To find the x-intercepts we set y = 0 and obtain
x2 = a2 and x = ±a. However, if we put x = 0 in above equation, we get
y 2 = −b2 , which is impossible, so there is no y-intercept. In fact, we obtain:
x2 y2
2
=1+ 2 ⩾1
a b
64
A.4. Hyperbolas
√
which shows that x2 ⩾ a2 and so |x| = x2 ⩾ a. Therefore we have x ⩾ a
or x ⩽ −a. This means that the hyperbola consists of two parts, called its
branches. It is sketched in Figure A.7.
x2 y2
Figure A.7: Hyperbola with equation a2 − b2 =1
In drawing a hyperbola it is useful to draw first its asymptotes, which are the
lines y = (b/a)x and y = −(b/a)x shown in Figure A.7. Both branches of the
hyperbola approach the asymptotes; that is, they come arbitrarily close to the
asymptotes.
By interchanging the roles of x and y we get an equation of the form
y2 x2
− =1
a2 b2
which also represents a hyperbola and is sketched in Figure A.8:
x2 y2
Figure A.8: Hyperbola with equation a2 − b2 =1
65
A. Math Appendix
A.5 Angles
Angles can be measured in degrees or in radians (abbreviated as rad). The
angle given by a complete revolution contains 360◦ , which is the same as 2π
rad. Therefore:
π rad = 180◦
and
◦
180
π
1rad = ≈ 57.3◦ 1◦ = rad ≈ 0.017rad
π 180
Example 1
(a) Find the radian measure of 60◦ . (b) Express 5π/4 rad in degrees.
Solution
(a) From Equation 1 or 2 we see that to convert from degrees to radians we
multiply by π/180. Therefore
π π
60◦ = 60 = rad
180 3
(b) To convert from radians to degrees we multiply by 180/π. Thus
5π 5π 180
rad = = 225◦
4 4 π
In calculus we use radians to measure angles except when otherwise indicated.
The following table gives the correspondence between degree and radian
measures of some common angles.
Degrees 0◦ 30◦ 45◦ 60◦ 90◦ 120◦ 135◦ 150◦ 180◦ 270◦ 360◦
Radians 0 π
6
π
4
π
3
π
2
2π
3
3π
4
5π
6 π 3π
2 2π
Figure A.9 shows a sector of a circle with central angle θ and radius r subtending
an arc with length a. Since the length of the arc is proportional to the size of
the angle, and since the entire circle has circumference 2πr and central angle
2π, we have:
θ a
=
2π 2πr
66
A.5. Angles
Remember that above equations are valid only when θ is measured in radians.
In particular, putting a = r in above equation, we see that an angle of 1 rad is
the angle subtended at the center of a circle by an arc equal in length to the
radius of the circle (see Figure A.10).
Figure A.10: Sector of a circle with its radius equal to the arc
The standard position of an angle occurs when we place its vertex at the origin
of a coordinate system and its initial side on the positive x-axis as in Figure
A.11 . A positive angle is obtained by rotating the initial side counterclockwise
until it coincides with the terminal side. Likewise, negative angles are obtained
by clockwise rotation as in Figure A.11.
3π 5π 3π 11π
− 2π = − + 2π =
4 4 4 4
and 2πrad represents a complete revolution.
67
A. Math Appendix
sin2 x + cos2 x = 1
tan2 x + 1 = sec2 x
1 + cot2 x = csc2 x
Double-Angle Formulas
Half-Angle Formulas
1 − cos x
r
x
sin = ±
2 2
1 + cos x
r
x
cos = ±
2 2
x sin x
tan =
2 1 + cos x
68
A.7. Polar Coordinates
p
r= x2 + y 2
y
θ = tan−1
x
Note that the angle θ must be adjusted to lie in the appropriate quadrant.
x = r cos θ
y = r sin θ
Position Vector
The position vector in polar coordinates is given by
r = rr̂
r̂ = cos θ, sin θ
Velocity Vector
To derive the velocity vector in polar coordinates, we differentiate the position
vector with respect to time:
dr d dr̂
v= = (rr̂) = ṙr̂ + r
dt dt dt
The time derivative of the unit vector r̂ can be found using the chain rule:
dr̂ d dθ
= cos θ, sin θ = − sin θ, cos θ = − sin θ, cos θ θ̇
dt dt dt
v = ṙr̂ + rθ̇θ̂
θ̂ = − sin θ, cos θ
69
A. Math Appendix
Acceleration Vector
To derive the acceleration vector in polar coordinates, we differentiate the
velocity vector with respect to time:
dv
a= = (r̈ − rθ̇2 )r̂ + (rθ̈ + 2ṙθ̇)θ̂
dt
d2 r d2 θ
where r̈ = dt2 and θ̈ = dt2 are the second derivatives of r and θ, respectively.
70
References
8) Roy, Archie E., and David Clarke. Astronomy: Principles and Practice.
CRC Press, 2018.
10) Schneider, Stephen E., and Thomas Arny. Pathways to Astronomy. McGraw-
Hill Education, 2021.
71