Astronomy & Astrophysics
Astronomy & Astrophysics
Astronomy & Astrophysics
AND
ASTROPHYSICS
ii TABLE OF CONTENTS
The CD-ROM that accompanies this book may only be used on a single PC. This
license does not permit its use on the Internet or on a network (of any kind). By
purchasing or using this book/CD-ROM package(the “Work”), you agree that this
license grants permission to use the products contained herein, but does not give you
the right of ownership to any of the textual content in the book or ownership to any of
the information or products contained on the CD-ROM. Use of third party software
contained herein is limited to and subject to licensing terms for the respective products,
and permission must be obtained from the publisher or the owner of the software in order
to reproduce or network any portion of the textual material or software (in any media)
that is contained in the Work.
INFINITY SCIENCE PRESS LLC (“ISP” or “the Publisher”) and anyone involved in the
creation, writing or production of the accompanying algorithms, code, or computer
programs (“the software”) or any of the third party software contained on the CD-ROM
or any of the textual material in the book, cannot and do not warrant the performance
or results that might be obtained by using the software or contents of the book. The
authors, developers, and the publisher have used their best efforts to insure the accuracy
and functionality of the textual material and programs contained in this package; we,
however, make no warranty of any kind, express or implied, regarding the performance
of these contents or programs. The Work is sold “as is” without warranty (except for
defective materials used in manufacturing the disc or due to faulty workmanship);
The authors, developers, and the publisher of any third party software, and anyone
involved in the composition, production, and manufacturing of this work will not be
liable for damages of any kind arising out of the use of (or the inability to use) the
algorithms, source code, computer programs, or textual material contained in this
publication. This includes, but is not limited to, loss of revenue or profit, or other
incidental, physical, or consequential damages arising out of the use of this Work.
The sole remedy in the event of a claim of any kind is expressly limited to replacement of
the book and/or the CD-ROM, and only at the discretion of the Publisher.
The use of “implied warranty” and certain “exclusions” vary from state to state, and
might not apply to the purchaser of this product.
TABLE OF CONTENTS iii
ASTRONOMY
AND
ASTROPHYSICS
A. B. BHATTACHARYA
S. JOARDAR
R. BHATTACHARYA
Revision & Reprint ©Copyright 2008 by INFINITY SCIENCE PRESS LLC. All rights reserved.
Original ©Copyright 2008 by Overseas Press India Pvt. Ltd. All rights reserved.
This publication, portions of it, or any accompanying software may not be reproduced in any
way, stored in a retrieval system of any type, or transmitted by any means or media, electronic or
mechanical, including, but not limited to, photocopy, recording, Internet postings or scanning,
without prior permission in writing from the publisher.
TABLE OF CONTENTS
Preface xv
v
vi TABLE OF CONTENTS
Chapter 4. Origin of the Sun, the Earth, and the Solar System 33
4.1 Introduction 33
4.1.1 The Solar Nebula Theory 34
4.1.2 The Planetsimal Theory 34
4.2 Age of the Earth 35
4.2.1 Radioactive Dating 35
4.3 Exposure Ages of Meteoroids 37
4.4 Age of the Radioactive Elements 38
4.5 Motions of the Planets 39
4.6 Emergence of Life on Earth 40
4.6.1 The Evolution of the Earth’s Atmosphere 41
4.6.2 The Formation of the Ozone Layer 42
4.7 The Role of Life in Changing the Earth’s Atmosphere 43
Exercises 45
Chapter 16. Some Unsolved Problems of Solar Meterwave Radio Bursts and
Some International Solar Radio Instruments 217
16.1 Introduction 217
16.2 Influence of Solar Radiations on the Earth’s Atmosphere 217
16.3 Global Electrical Circuit 218
16.4 Types of Solar Radiations 218
16.4.1 Solar-ionizing Radiation 219
16.4.2 Solar Radiation and D-Region of the Ionosphere 219
16.5 Solar Modulation of Atmospheric Electrification 220
16.6 Diffuse Solar Radiation and Associated Parameters 222
16.6.1 Instrumentation 222
16.6.2 Factors Controlling Net Radiation 223
16.6.3 Distribution of Sunshine hours, Normalized D and IFIA 224
16.7 Effects of Geomagnetic Disturbances of Solar Origin 226
16.8 Sprite Flashes 226
16.8.1 Visual Reports 226
xii TABLE OF CONTENTS
Chapter 18. The Post Main Sequence Period of the Stars 253
18.1 Introduction 253
18.2 The Red Giant Phase 254
18.2.1 The Red Giant Branch (RGB) 254
18.2.2 The Triple-alpha Process 256
18.2.3 The Horizontal Branch (HB) 256
18.2.4 The Asymptotic Giant Branch (AGB) 257
18.2.5 Thermal Pulses in AGB Stars 257
18.2.6 Mass Loss from AGB Stars 258
18.3 All Stars Pass Through the Red Giant Phase 259
18.3.1 Small Stars 260
18.3.2 Medium Mass Stars 261
18.3.3 Massive Stars 261
Exercises 262
TABLE OF CONTENTS xiii
Chapter 19. Brilliant Phenomena: Black Holes, Pulsars, Quasars, Dark Matters,
Galaxies, and Clusters 265
19.1 Introduction 265
19.2 Supernova 266
19.2.1 Type-I Supernova (Nova) 266
19.2.2 Type-II Supernova 268
19.3 Compact Objects of the Sky 269
19.3.1 White-Dwarfs 269
19.3.2 Neutron Stars 271
19.3.3 Pulsars (A Special Case of Neutron Stars) 273
19.3.4 Black Holes 274
19.3.5 Quasars 278
19.4 Evidences of Dark Matter 280
19.5 Galaxies 282
19.5.1 Classification of Galaxies Based on Hubble Sequence 283
19.5.2 General Classification and Properties of Galaxies 284
19.5.3 Our Milkyway Galaxy 289
19.5.4 Rotation of Galaxies 293
19.6 Galaxy Groups, Clusters, and Super Clusters 294
19.7 Gravitational Lensing (Einstein Effect) 296
19.8 Red-Shifts 299
19.8.1 Mechanisms for Red-Shifts and Blue-Shifts 299
19.9 Hubble’s Law in Relation to the Expanding Universe 301
Exercises 301
References 345
Index 349
TABLE OF CONTENTS xv
PREFACE
T
he subjects of astronomy and astrophysics have grown tremendously in
scope during the last three decades. In fact, over the past few years there has
been an explosive growth of these subjects, in particular, about new ideas
and concepts as well as a large increase in available information and data in all
of the solar technologies. This demands for a comprehensive textbook describing
the present state of knowledge, offering the different techniques and instruments
developed in this relatively young science. It is hoped that the present book will
fulfill this requirement to a large extent.
The subject may be treated traditionally under a single banner; however, at
the present level of our understanding it can be well classified into two distinct
parts: (i) different techniques and instruments employed and the results obtained
thereof and (ii) purely theoretical approaches. The present volume deals in a
comprehensive manner both the categories where suitable texts are highly lacking.
This volume has been prepared primarily for post-graduate students of different
universities, though many of the sections will prove useful for practicing engineers
and scientists actively engaged in advanced research.
To ensure the highest degree of reliability, besides a large number of reference
books, we have consulted freely with many specialists who have provided an
up-to-date overview of their areas and have given very valuable information
and ideas. The authors are grateful for their indirect contribution to the text. We
would like to acknowledge Sri Dibyendu Joardar, for drawing many complicated
diagrams. Finally, we would like to take this opportunity to express our deepest
gratitude to many of the astronomers and other scientists whose works are directly
or indirectly included in the text.
— Authors
xv
iv TABLE OF CONTENTS
Chapter 1
ASTRONOMICAL UNIVERSE AND BASIC
TECHNIQUES OF RADIO ASTRONOMY
1.1 INTRODUCTION
T
he recorded history of astronomy reveals that in the beginning there was
almost no difference between astronomy and astrology. The first references
are found in the Indian astronomy-astrology, in the Rig Veda (2000 B.C1).
Around 1800 BC, Yajnavalkya advanced a 95-year cycle to synchronize the
motions of the Sun and the Moon. Lagadha wrote a text on astronomy in 1350 BC.
Horoscopes contained the positions of Sun, planets and nodes of the Moon (Rahu
2
and Ketu) . The stars were recognized as the Sun, with the Sun being the center
of the Universe. Aryabhatta in the 5th century propounded the heliocentric theory
of gravitation, later translated as “Aryabhatta’s Magnum Opusthus.” He wrote
1,582,237,500 rotations of the Earth equals 57,753,336 lunar orbits. It is perhaps
the oldest accurate fundamental astronomical ratio. The 7th century BC book by
1
The Vedas existed much before 2000 BC but not in the form of written text. In fact they were
handed down by teaching and memorizing from generation to generation.
2
The two points formed by the intersection of the ecliptic’s great circle with the same of the
Moon’s orbit as seen from the Earth (eclipse points).
1
2 ASTRONOMY AND ASTROPHYSICS
3
1 Yojan = 7.2 miles.
ASTRONOMICAL UNIVERSE AND BASIC TECHNIQUES OF RADIO ASTRONOMY 3
Mercury, Venus, Earth, Mars, Jupiter, and Saturn are only six of the nine planets
that revolve around the Sun; the other three are Uranus, Neptune, and Pluto, though
Pluto is controversial as a planet. However, planets and their satellites are not the
major constituents of the Universe. The Sun is a star, the closest one to Earth, and
it contains about 99.9 percent of the mass of our solar system. Our Galaxy contains
11
more than 10 stars and it is only one of myriads in the observable Universe. The
observable Universe contains about 1010 galaxies. Thus, in the observable Universe
there are about 1010 galaxies × 1011 stars/galaxy = 1021 stars.
Today, we know that the Universe is expanding and the major constituents of the
matter are hydrogen and helium with a smidgen of heavier elements like oxygen,
carbon, and iron. If we extrapolate back in time by 10 to 20 million years, we find
that the density of the Universe was almost infinite. The entire space was filled
with thermal radiation at 3 K. In 1965, Penzias and Wilson discovered the cosmic
microwave background radiation. These facts agree with the possibility that
the Universe started with a sudden decompression, or a Big Bang. It happened
simultaneously everywhere in space. The standard model of the Big Bang Theory
is explained in the following eight steps.
1. Before the Universe: Before space and time existed, the entity of the Universe
was compressed into the confines of an atomic nucleus also known as a
singularity. The temperature was extremely high of the order of trillions of
degrees or more and the density was almost infinite.
2. The Big Bang: The creation process of the Universe began with an expansion
and not an explosion as the name suggests.
3. The Inflationary Epoch (immediately after the Big Bang): There is no
compelling physical model for the first 10–35 seconds. However, we give a picture
for some possible happening during this period. The four forces of nature:
strong nuclear, weak nuclear, electromagnetic, and gravity were combined
as a single super force called Wald. Elementary particles known as quarks
begin to bond in trios forming photons, leptons, antileptons, gluons, photons,
neutrinos, and antineutrinos. There were smaller amounts of neutrons and
protons, approximately 1 in every one billion photons, electrons, or neutrinos.
At this stage, the density of the Universe was assumed to be nearly 1094 g/cm3
with the majority of this being radiation. The Universe was expanding much
faster than the speed of light during this creation and annihilation of these
particles. The rate of inflation of the Universe was assumed to be doubling
every 10–5 seconds.
4 ASTRONOMY AND ASTROPHYSICS
4. Nanoseconds after the Big Bang: At this point the Universe was in a ionized
plasma state with radiation and matter inseparable. There were equal amounts
14
of particles and antiparticles. As the temperature fell below 10 K, the quarks,
antiquarks, and gluons annihilated and transformed into less massive particles.
The quarks were slightly more in number than the antiquarks and as a result
some quarks were left behind to form protons and neutrons. The leptons and
antileptons were annihilated with the fall in temperature.
5. Between a millisecond to a second after the Big Bang: As the temperature
fell around 1010, the density of the Universe drastically reduced to neutrinos
to effectively interact with the matter. These non-interacting decoupled
neutrinos formed a universal gas within the Universe.4 Immediately after
this, the annihilation of electron-positron pairs removed all positrons and
most electrons.
6. 100 seconds after the Big Bang: As the temperature fell to 109 K, protons and
neutrons slowed down enough to allow nucleosynthesis to take place. Atomic
nuclei of helium were produced as two protons and neutrons each bonded. For
every helium nucleus formed, there were about ten protons left over allowing
for 25 percent of the Universe to be comprised of helium. The light nuclei were
formed by the combination of neutrons and protons.
7. Thirty minutes after the Big Bang: At this stage an important phase of the
expansion occurred when the creation of photons increased through the
annihilation of electron-positron pairs. It is fortunate that the Universe began
with slightly more electrons than positrons which insured that our Universe
was able to form the way it has.
8. Three-hundred thousand years after the Big Bang: The Universe continued
to expand and cool until it reached a temperature of 10,000 K. It allowed for
helium nuclei to absorb free-floating electrons and to form helium atoms.
Also, hydrogen atoms were bonding together and forming lithium. As the
temperature reduced to 4000 K, hydrogen and helium nuclei combined with
electrons to form neutral hydrogen and helium atoms. It was here that the
density of the Universe had reduced so much that light could be perceived.
Until this point photons continued to be trapped within matter. Thus, the
Universe became transparent to electromagnetic radiation.5
4
Because of the expansion of space this gas has cooled to a temperature of 2 K as of today.
However, until today it has not been possible to detect this universal background of neutrinos.
5
The radiation freed from interacting with matter has now cooled to a temperature of 3 K due to
the expansion of space, which is the cosmic microwave background radiation.
ASTRONOMICAL UNIVERSE AND BASIC TECHNIQUES OF RADIO ASTRONOMY 5
We now consider the scales of some astronomical objects. Listed in Table 1.1 are
the very rough sizes, masses, and ages of the Sun, the Galaxy, and the observable
Universe. For comparison, the quantities are listed here with that of a small child.
17
Table 1.1 gives the age of the Universe as somewhat greater than 10 sec.
Astronomers believe that the time since creation is between 10 and 20 billion
years. Table 1.1 notes that these estimates apply to the observable Universe. This
distinction between the Universe and the observable Universe is required because
astronomers have not yet settled the question of whether the Universe is open
or closed. If the Universe is open, space-time (and the number of stars in the
Universe) is infinite. As time had a beginning, but no end, the amount of space
and the number of stars in it are actually infinite.
The inhabitants of the Universe are mainly stars of which the Sun is only one
example. Besides suns of all colors there are the tremendously distended red giants
and the mysteriously tiny white dwarfs. The Sun itself is an average star, neither
very massive nor very light.
Varieties of stars live in our Galaxy either singly or in pairs. Sometimes,
hundreds or thousands of stars can be found in loose groups called open
clusters. The Pleiades is a fairly young open cluster, and the reflection nebulosity
surrounds the seven sisters of this cluster. The oldest stars in our Galaxy are found
in tighter groups called globular clusters, which may contain more than a million
members.
The space between the stars is not completely empty. The diffuse matter
between the stars is called the interstellar medium. Clouds of gas and dust, as well
6 ASTRONOMY AND ASTROPHYSICS
Undoubtedly, the human eye is a poor photon detection and storage instrument.
All types of electronic equipment help the astronomer (i) measure the brightness of
astronomical objects (photometry), (ii) calculate the wavelength distribution and
spectral features (spectrophotometry and spectroscopy), (iii) assess the amount
of polarization (polarimetry), and (iv) fix the position of the object in the sky
(positional astronomy). The one basic tool of modern astronomy is definitely the
computer. This can be used to run many instruments more or less automatically,
guide the telescope, reduce data, and help construct theories. The basic astronomical
measurements are listed in Table 1.2.
Proper motion The real angular shift of nearby objects against the positions
of distant objects. This shift yields the component of velocity
perpendicular to the line of sight provided the distance to
the object is known.
ASTRONOMICAL UNIVERSE AND BASIC TECHNIQUES OF RADIO ASTRONOMY 7
The basic tool for exploring the radio Universe is a radio telescope. Essentially,
it consists of a reflector antenna, a feed, and a receiver system. The reflector
focuses the radio waves from a celestial object on to the feed. The feed is usually
a dipole or a horn that collects the radio signals and feeds these weak signals to
a sensitive low-noise amplifier where it is further processed by a sophisticated
electronic receiver system. The two fundamental characteristics of a telescope are
(i) high sensitivity to detect the weak radio waves from distant celestial objects and
(ii) high angular resolution to distinguish features close together in the sky. The
larger the collecting area or size of a telescope the greater its sensitivity as well as
its resolving power.
Owing to the very long wavelengths of radio waves, ranging between about 1 cm
and 10 m, the resolution of even a large, single radio telescope is relatively much
poorer. At a wavelength of 1 m, for instance, a parabolic dish with a diameter of
about 1 km would be required to obtain even the modest resolution of the human
eye. Radio astronomers have overcome this limitation by development of an
interferometry technique in which the voltage outputs of a pair of modestly sized
radio telescopes are brought together to a common point and multiplied together.
By using an array of telescopes spread over tens of kilometers such interferometric
techniques can be utilized to synthesize much larger telescopes with effective
diameters equivalent to the maximum separations or baselines between the most
distant antennas in array. This process of aperture synthesis is greatly aided by the
Earth’s rotation, which constantly changes the orientation and effective lengths of
different interferometric baselines as viewed from the radio source.
Radio astronomers can now obtain resolutions as high as a thousandth of an arc
second through telescopes by the technique of Very Long Baseline Interferometry
(VLBI). In VLBI radio signals are independently recorded by telescopes separated
by thousands of kilometers and subsequently correlated in a computer. The
resolution is, of course, much higher than possible with optical telescopes.
A major part of the pioneering work in Radio Astronomy after World War II was
carried out at meter wavelengths. However, the development of interferometric
10 ASTRONOMY AND ASTROPHYSICS
techniques and low-noise receivers found a steady shift in emphasis toward shorter
wavelengths, where the Galactic background noise is smaller and higher angular
resolution can be attained for a given antenna size or interferometer baseline. Most
powerful radio-astronomical facilities in the world today operate in the centimeter
and decimeter part of the radio spectrum even though there are many exciting
astrophysical problems that are best studied at meter wavelengths. The importance
of this part of the radio spectrum has been highlighted by a variety of astronomical
research carried out using instruments built for meter-wave operation. In general,
these instruments are limited by sensitivity and resolution and limited frequency
coverage.
The largest radio telescopes constructed over the past two decades have been
designed to operate mainly at centimeter and decimeter wavelengths. There are
several reasons for this: (i) a radio telescope of a given aperture provides higher
resolution at shorter wavelengths; (ii) a larger collecting area is necessary to
reach similar sensitivities at larger wavelengths due to a steep increase of the
galactic background temperature with wavelength; (iii) ionospheric irregularities
introduce increasingly serious distortion in the wave-front from a distant radio
source at wavelengths longer than about 20 cm; (iv) there is a great deal of
man-made radio interference in the allocated radio astronomy bands at longer
wavelengths that at shorter wavelengths, particularly in the western countries.
Meter-wave astronomy, on the other hand, offers certain distinct advantages. It is
sufficient to achieve only a moderate resolution at meter wavelengths in order to
study the diffuse emission arising from the oldest relativistic electrons in galactic
and extra-galactic radio sources. The radiation from these becomes dominant at
meter wavelengths and is a useful probe of the evolution of the parent objects.
Moreover, large antennas required for meter-wave operation can be constructed
quite economically as their reflecting surface does not require high precision.
Harmful effects of ionospheric irregularities at longer wavelengths can be largely
overcome through the application of methods like self-calibration that have been
developed for use with synthesis radio telescopes over the last decade.
To locate stars and other objects in the sky we consider all the objects to lie on
the surface of a sphere with ourselves at the center of it. This imaginary sphere
is known as the celestial sphere. For simplicity, we fix the radius of the sphere as
unity. The direction of an object in the sky is then given in terms of its position on
the celestial sphere. The point and the direction of the object can then be assigned
two angle coordinates. As the principle underlying the astronomical coordinate
systems is the same as the terrestrial system, we consider it first.
ASTRONOMICAL UNIVERSE AND BASIC TECHNIQUES OF RADIO ASTRONOMY 11
Figure 1.2: The coordinate system used for astronomical observations from Earth.
In Fig. 1.2, let NQSP be the globe where N represents the North Pole and S
the South Pole. The plane of the circle QLKP is perpendicular to the line joining
N and S and passes through the center of the sphere. This plane is the equator.
The circles whose planes pass through the center of the sphere are known as great
circles while if the plane of the circle does not pass through the center, it is called
a small circle. Any half-circle through N and S is a meridian and the meridian
through Greenwich is the principal or prime meridian. The point of intersection of
the prime meridian with the equator is the reference point to measure arc lengths
along the equator. The meridian through the fixed point X cuts the equator at L.
The longitude (λ) of X is then called the spherical angle GNX or the angle KOL in
the plane of the equator or the arc length KL. The other coordinate of X is called
the latitude (φ). It is the arc length LX or the angle LOX. Latitude ranges from 0 to
90 degrees North or South, accordingly, as the plane is in the North or the South
hemisphere. Likewise, the longitudes range from 0 to 180 degrees East or West
depending on whether the plane is to the East or to the West of Greenwich. The 360
degrees of longitude is equal to a complete rotation of the Earth on its axis, which
takes 24 hours. Hence, 15 degrees of longitude becomes equal to one hour of time
(i.e., 1 degree of longitude becomes equal to 4 minutes of time).
The longitude-latitude system provides the essential requirements of an
astronomical coordinate system, which are (a) a fundamental great circle (e.g., the
12 ASTRONOMY AND ASTROPHYSICS
equator), and (b) a reference point or origin (e.g., the point K). Depending on the
fundamental circle and reference point, the astronomical systems of coordinates
are: (i) the Horizon System, (ii) the Equatorial System, (iii) the Ecliptic System, and
(iv) the Galactic System.
Fig. 1.3 (a), (b), (c), and (d) represent, respectively, the above four coordinate
systems used in astronomy. Fig. 1.3 (a) illustrates the Horizontal System of the
coordinates for astronomical objects. In it, SZNZ is the celestial sphere while point
L is the intersection of the great circle through X and the horizon. The azimuth A is
to be measured from the North point N. Fig. 1.3 (b) illustrates the equatorial system
of coordinates, where Y is the first point of Aries. In it, PXQ is a great circle through
X. The horizon system and the spherical triangle ZXP are shown. This gives the
relation between the Equatorial and Horizontal systems of coordinates. Fig. 1.3 (c)
explains the Ecliptic System of coordinates. The spherical triangle MXP establishes
the relation between the Ecliptic and the Equatorial systems of coordinates.
Fig. 1.3 (d), on the other hand, shows the Galactic system of coordinates. In it, point
A indicates the direction of the center of the Galaxy.
Different units are used for measuring the distances between astronomical objects
depending on their positions. For measurements within our solar system, or other
solar systems, the common unit used is the astronomical unit. 1 astronomical unit
= the average distance between the Sun and the Earth (≈ 150,000,000 km). This is
very small when measurements are to be made for galactic objects like the stars
where another unit called parsec (pc) is used. 1 pc = 3.26156378 light years (≈ 206265
AU), which is nearly equal to the distance of the nearest star. The term parsec is the
short form of parallax of one arc-second. Note that if an object has a parallax of 0.5
arc-second, it will be twice the distance of an object with 1 arc-second of parallax,
or about 2 parsecs away (6.52 light years). Conversely, an object with 2 arc-seconds
of parallax will be twice as close as an object with 1 arc-second of parallax, or
about half a parsec away (1.63 light years). For measuring distances within our
galaxy or beyond (other galaxies) kilo-parsec (kpc) is used. 1 kpc = 1000 pc. For
distances between galaxies and the cosmology Mega-parsecs (Mpc) are used.
1 Mpc = 1,000,000 pc. It may be noted that though a light-year is a large distance,
it is generally not used in present day astronomy.
ASTRONOMICAL UNIVERSE AND BASIC TECHNIQUES OF RADIO ASTRONOMY 13
Figure 1.3: The various coordinate systems. (a) Horizon System, (b) Equatorial System, (c) Ecliptic System,
and (d) Galactic System of coordinates for astronomical objects.
EXERCISES
1. What was the explanation of Copernicus for the complicated motion of the planets
in terms of the stationary Sun? How did Galileo and Kepler adopt the Copernicus
system?
2. How did Newton establish the laws of planetary motion in a Sun-centered system?
3. How many stars are there in our Galaxy and in the observable universe?
4. What is Big Bang?
5. Explain the standard model of the Big Bang theory in eight different steps.
6. Who first discovered the cosmic microwave background?
7. Provide in a tabular form: very rough sizes, masses and ages of the Sun, the Galaxy
and the observable universe. You may compare the quantities with respect to a small
child.
8. Which are the oldest stars in our Galaxy?
14 ASTRONOMY AND ASTROPHYSICS
2.1 INTRODUCTION
A
s we approach the solar system the first bodies we encounter are frozen
balls of gas and dust that exist in the billions. At least one of those frozen
balls has wandered too close to the Sun. It has taken the familiar form of
a comet. Pluto6 is the most distant planet from the Sun, but its eccentric orbit has
actually carried it inside the orbit of Neptune. Also, inside Neptune’s orbit lies
Uranus, the seventh planet of the solar system. Pluto, Neptune, and Uranus are
all poorly known.
The sixth planet, Saturn, is bright enough to be visible. The next planet, Jupiter,
is more massive. Its face is covered by an intricate system of zones and belts. Jupiter
is surrounded by a retinue of satellites including 15 known moons.
Inside the orbit of Jupiter lies a vast wasteland, the asteroids. Closer to the Sun
lies the fourth planet, Mars. The third planet is the Earth, the second planet Venus,
and closest to the Sun is Mercury.
6
Currently, Pluto is under controversy in the astrophysical world. It is not accepted as a planet
anymore, though we have given its old heritage in various places in this book.
15
16 ASTRONOMY AND ASTROPHYSICS
At the center of the solar system lies the mighty Sun. Thus, as we swing around
the Sun, we pass Mercury and Venus, and approach the third planet, the Earth,
which is nothing but an insignificant speck of dust in this immense Universe.
Table 2.1 gives the properties of the nine planets that circle around the Sun. It
shows that the planets divide into two categories: (i) the inner or terrestrial
planets, which are small and have a mean density of 4 to 5 gm/cm3; and (ii) the
outer or Jovian planets, which are large (except Pluto) and have a mean density of
1 to 2 gm/cm3.
Each of the Jovian planets, i.e., Jupiter, Saturn, Uranus, and Neptune, has more
than one moon. In fact, Jupiter has at least 15 moons. Saturn has at least 22 moons.
Uranus has five known moons; and Neptune, two. In addition to satellite systems,
the Jovian planets may also all possess rings (e.g., the rings of Uranus are shown
in Fig. 2.1).
Figure 2.1: The rings of Uranus were revealed by studies of the occultation of a distant star. The lower
light curve shows schematically the observed light intensity of the star as its position relative
to a ringed planet changes. Such light curves can reveal not only the presence of occulting
ring material but also much information about the composition and structure of the planets
atmospheric layers. A central flash occurs when the star is directly behind the center of the
planet and is caused by refraction and extinction in the planets atmosphere acting like a lens in
this configuration.
The terrestrial planets, i.e., Mercury, Venus, Earth, and Mars, have no rings and
very few moons. Mercury and Venus have no moon, Earth has one (Luna or the
Moon), and Mars has two (Phobos and Deimos). Furthermore, in comparison with
the large moons of the Jovian planets, the moons of the terrestrial planets, except
for Luna, are quite small (Table 2.2). As the two small moons of Mars are captured
asteroids, the Earth-Moon system is unique among the terrestrial planets.
18 ASTRONOMY AND ASTROPHYSICS
Mercury none - - - -
Venus none - - - -
⊕
A gravitational mass estimate for Phobos is available from Viking Orbiter and suggests that
Phobos has chondritic material.
⊗
Voyager-1 measurements showed that Titan is slightly smaller than Ganymede contrary to
previous belief. It remains for spacecrafts to settle whether Triton or Ganymede is the solar
system’s biggest moon.
THE SOLAR SYSTEM 19
heat up from the friction generated in the passage. It will then appear as a fiery
shooting star called a meteor. If the mass of the meteor is less than 10–10 gm, it may
slow down so fast that it survives the flight; if so, it is called a micrometeorite.
On the other hand, if the mass of the meteor is greater than 103 gm, the object is
called a meteorite. Meteorites come in three types depending on their chemical
composition: stones, stony irons, and irons. Stones resemble rocks; stony irons
have some metal-rich inclusions, while irons contain mostly metals like iron and
nickel.
Figure 2.2: The collision of two asteroids shatters them into many pieces. Some of these may suffer resonant
interactions with Jupiter and within 104 to 105 years, enter Earth-crossing meteoroids.
2.4 COMETS
The word comet is derived from the Latin word cometa, which means long-haired.
Comets are astronomy’s most spectacular displays. Fig. 2.3 shows schematically
some features found in comets; a head consisting of two parts: a nucleus and a
coma. In addition it has one or two tails. The nucleus of the comet is the essential
part, as it is the ultimate source of all the mass and is generally less than 50 km
across. The coma could be larger than the Sun and the ion tails have been observed
to extend 150 million km (1 astronomical unit) or more. Whipple proposed that the
20 ASTRONOMY AND ASTROPHYSICS
nucleus is composed of chunks of dust and frozen ices of compounds like methane
(CH4), ammonia (NH3), water (H2O), and carbon dioxide (CO2). Comets are often
called dirty snow balls based on their composition.
Before our solar system came into existence, the interstellar cloud out of
which the comets were formed had uncountable trillions of these bodies mixed
in with a cloud of various gases of carbon monoxide, hydrogen, formaldehyde,
ammonia, and so on. The process by which this diffused gas clumped into comets
still unknown. These bodies probably measured from centimeters to tens of
kilometers. Today, all the remnants of the original interstellar cloud around the
solar system are these comets. Sometimes comets visit the inner solar system while
moving in their orbit. We shall continue to see these bodies even after most of the
asteroids have vanished. These objects can also be a threat to our existence as they
could collide with the Earth like the comet Shoemaker-Levy, which collided with
Jupiter.
Figure 2.3: The structure of a comet. (a) Details of a comet. (b) The orbit of a comet near the Sun.
The orbits of comets are constantly changing since their origin is outside the
solar system. Some comets move into Sun-grazing orbits and get destroyed when
they near the Sun, while others get thrown out of the solar system permanently.
Comets are identified from their periods, i.e., short-period comets or long-period
comets. The former have orbits of less than 200 years, while the latter have longer
orbits (thousands of years) but remain gravitationally bound to the Sun. These are
THE SOLAR SYSTEM 21
also called periodic comets since their orbits are elliptical. Due to their low masses
and elliptical orbit, the comets are pulled close to the giant planets and their orbits
are often perturbed. Non-periodic or single apparition comets have parabolic or
hyperbolic orbits, which make them permanently exit the solar system after a
single pass by the Sun.
As the comet approaches near the Sun, both the coma and tail gets illuminated
as they reflect sunlight and could be visible from the Earth. In general, most comets
are extremely faint and can be seen only through telescopes. It is surprising that
comet nuclei are among the darkest objects known to exist in the solar system
since they reflect only 4 percent of the light falling on them. In 1996 it was also
found that comets emit x-rays. These are thought to be generated by the interaction
between comets and the solar wind when highly charged ions penetrate the comet’s
atmosphere and collide with the atoms and molecules. The collision process results
in the capture of one or more electrons by the ions leading to emission of x-rays
and far ultraviolet photons.
Some well-known comets are Comet Borelly, Comet Enche, Comet Hale-Bopp,
Comet Halley (also known as Halley’s comet), Comet Hymason, Comet Ikeya-
Seki, Comet Mrkos, and Comet Shoemaker-Levy 9.
EXERCISES
3.1 INTRODUCTION
W
e know that the atmosphere of a terrestrial planet is the mass of gas that
lies above the surface of the planet. But a Jovian planet does not have
a solid surface. Also, a star like the Sun has no solid surface. But we
have no difficulty distinguishing between its interior and its atmosphere. By the
atmosphere of an astronomical body we mean the layer where the photons make
a transition from walking to flying. This definition is fairly unambiguous for a star
as most of the ambient photons have similar wavelengths. For a planet, sunlight
provides a source of incoming radiation, which has photons in the visual and
near-ultraviolet regime. This radiation field affects the atmosphere qualitatively.
Radio, infrared, optical, and ultraviolet wavelengths also affect the atmosphere of
a planet.
Mars has only a very thin atmosphere while Mercury and the Moon have no
atmosphere to speak of. The temperature structure of the Martian atmosphere is
23
24 ASTRONOMY AND ASTROPHYSICS
affected by the direct absorption of sunlight by dust particles swept into the air
(Fig. 3.1). The terrestrial planets with gaseous atmospheres are Venus and Earth.
Fig. 3.2 shows a schematic plot of the temperature of the Earth’s atmosphere
as a function of height above the ground. The thermal structure of the Earth’s
atmosphere has four distinct regions:
(i) The troposphere is the region that adjoins the ground. It is the region with
which we associate weather. Vertical heat transfer in the layer consists
mostly of the upward diffusion of infrared photons and fluid transport
by convection. The temperature is highest near the ground and steadily
decreases to the tropopause, the top of the troposphere.
(ii) The stratosphere is heated by the ultraviolet component of sunlight in
forming and destroying the ozone O3. This direct absorption leads to an
inversion of the temperature above the tropopause, so that the stratosphere
is actually warmer than the troposphere. This increase in temperature with
height continues until the stratopause is reached.
Figure 3.1: Thermal structure of the atmosphere of Mars is affected by the amount of dust swept into
the air by the winds of that planet. The shaded region between the curves demarcates the
theoretical range of temperatures observed in calculations using typical daytime conditions for
a dust-laden and clear Martian atmosphere. The dashed curves indicate the limits found by
aircraft observations of Mars.
of this heat input leads to a rise of the temperature with increasing height
in the thermosphere until temperatures of about 1,000 K are reached.
The primary energy transfer in the lower thermosphere is a downward
conduction of heat into the mesosphere.
As the atmosphere is transparent to optical photons, sunlight strikes the
ground, with little interaction with the troposphere except for scattering by
molecules, dust, and reflection from cloud tops. The troposphere can be regarded
as a gaseous atmosphere. It is heated from below by the emanation of infrared
radiation from the ground. The troposphere is partially opaque to this emergent
infrared radiation. In planetary astronomy, this is known as the greenhouse effect,
but the basic principle is the same as that which governs the radiative transfer
inside a star.
Figure 3.2: The thermal structure of the Earth’s atmosphere. The various regions are presumably duplicated
qualitatively in the other planets. The detailed mechanisms of ultraviolet heating in the upper
atmospheres of the other planets may, however, differ from those of the Earth due to the
different chemical compositions involved.
As the Sun does not shine equally on the Earth’s equator and poles, we
expect some variation of this temperature with latitude. Measurements show that
the tropopause at the tropics and polar caps have different temperatures. The
tropopause is warmer at the poles (about 225 K) than at the tropics (about 195 K).
26 ASTRONOMY AND ASTROPHYSICS
Figure 3.3: The greenhouse effect in planetary atmospheres involves the additional warming of the surface
layers due to the partial trapping of the remitted infra-red radiation. Once it was assumed that
the name for this effect might be a misnomer as greenhouses may not be kept warm by the
same mechanism. However, a reanalysis by Silverstein reveals that the greenhouse effect may
be nontrivial even for greenhouses. In any case, the effect is the same as that which keeps the
interiors of stars warmer in their interiors than at their surfaces.
Venus is 1.38 times closer to the Sun, but it reflects about 72 percent of the
incident sunlight while the Earth reflects 39 percent. Thus, the average flux of
sunlight absorbed by Venus is (1.38)2(0.28/0.61) = 0.87 of that absorbed by Earth
1/4
(i.e., the effective temperature of Venus should only be (0.87) that of Earth,
THE ATMOSPHERES OF THE PL ANETS 27
or about 238 K). Venus has whopping greenhouse effects and has a whopping
thick atmosphere. Surface pressure is about 90 times that at sea level on Earth.
This atmosphere is composed of 96 percent of carbon dioxide (CO2), which is an
efficient absorber in the infra-red. Venus has a much thicker atmosphere than the
Earth with a composition heavily weighted toward CO2.
The crucial factor to be accounted in this connection is the existence of liquid water
on the surface of Earth and not on Venus. Sedimentary rocks have a capacity to
absorb carbon dioxide if there is liquid water present. The Earth has an estimated
equivalent of 70 atmospheres of CO2 locked up in rocks like lime-stone (CaCO3).
The corresponding amount of CO2 in Venus is actually in the atmosphere.
Simpson showed that as long as there is liquid water on the surface of a planet,
the top of the convection zone in the troposphere will hold as much water vapor
as it physically can. This is because more water can be held in warm air than in
cold air. The layers above the top of the convection zone are convectively stable.
So there is no descent of dry air into the top of such a convection zone. This means
that the air at the top tends to be exactly saturated, with the excess water going
into the formation of clouds of small droplets of liquid water or particles of ice.
Water vapor has a saturation temperature T, which is a function of the partial
pressure PH O. This can be expressed as,
2
T = TS(PH O). (3.1)
2
The temperature above the top of the convection zone, on the contrary, follows
the radiative equilibrium law
1/4
⎡3 2 ⎤
T = Te ⎢ ⎛⎜ τ + ⎞⎟ ⎥ . (3.2)
⎣ ⎝
4 3 ⎠⎦
Here, Te is the planetary effective temperature. If the radiative equilibrium
keeps the top of the convection zone too warm to allow water vapor to saturate
at the top of the convection zone, more water will be drawn from the surface by
the ascending convection currents. The dwindling supply of liquid water to form
sedimentary rocks would release more CO2 into the atmosphere. The increase
of both CO2 and H2O concentrations in the atmosphere would increase the
greenhouse effect. The runaway would not be halted until all the available surface
water is released into the atmosphere and the surface of the planet is left in an arid
and sweltering state.
Scientists have speculated that the opposite of a runaway greenhouse effect
may have occurred on Mars. Spectroscopic and photometric studies suggest a
28 ASTRONOMY AND ASTROPHYSICS
surface pressure of about 0.6 to 0.7 percent that of the Earth. A surface temperature
varies greatly with latitude and between day and night, but it never exceeds 300 K
and may fall below 150 K (i.e., the freezing point of carbon dioxide when it forms
the dry ice at the surface pressure of Mars). Mars’ atmosphere consists of 95 percent
carbon dioxide (CO2), 2.7 percent nitrogen molecules (N2), 1.6 percent argon (Ar),
and small amounts of free oxygen (O2) and water vapor (H2O). The average
surface pressure and temperature on Mars are too low to allow the presence of
liquid water. However, it is likely that Mars once had a much denser atmosphere,
and much higher surface temperatures. Most of the water and carbon dioxide now
seems to be locked up either in the polar caps or in carbonate rocks.
An earlier and denser atmosphere on Mars did allow the presence of liquid
water on the planet. However, the gradual locking up of the atmospheric carbon
dioxide into sedimentary rocks led to a thinning of the atmosphere, with a
consequent lowering of the temperature. This caused more water vapor to precipitate
out of the atmosphere resulting in even more carbon dioxide being locked up.
Fig. 3.4 shows a schematic diagram of the greenhouse effect for Venus, Earth,
and Mars. The solid curves represent the surface temperature of the planet,
assuming an opacity associated with various values of the vapor pressure of water
at ground level. The dashed curves show the loci of thermodynamic equilibrium
possible for the three phases of water. Water opacity cannot raise the surface
temperature on Mars higher than a point where its vapor becomes saturated with
respect to ice, and on Earth, with respect to liquid water. For Venus, a runaway
greenhouse effect occurs because the solid curve is incompatible with the existence
of liquid water on the surface of the planet.
Figure 3.4: A schematic diagram of the greenhouse effect for Venus, Earth, and Mars.
THE ATMOSPHERES OF THE PL ANETS 29
Molecules cannot escape from the lower levels of a planetary atmosphere any
more than a photon in the deep interior of the Sun can fly out in a straight line
(Fig. 3.5). A fast molecule would collide with another molecule before it had made
any substantial progress toward escape. The path followed by any molecule in
these lower regions is a tortuous random walk. Only in the rarefied outer portions
of a planetary atmosphere can a molecule fly freely and have a chance to escape
the gravitational clutch of the planet. It is in these regions where a newly formed
hydrogen atom can escape before recombining with another hydrogen atom.
The region where fast molecules change from walking to flying is known as the
exosphere. Exospheres are to molecules on planets what photospheres are to
photons on stars. However, all photons can escape, as they travel at the speed
of light. With molecules, only that fraction in the exosphere can actually escape
whose velocity exceeds the escape speed the planet.
Figure 3.5: Deep in a planetary atmosphere, the mean free path for molecule collisions is so small that even
a very fast molecule has no chance to escape from the planet. Above the base of the exosphere,
the density of molecules has dropped low enough that the probability for escape becomes high
for fast molecules.
1/2
Here, vT = ⎛⎜
kT ⎞
⎟ is the thermal speed. The dimensionless function f(xe) is
expressed as ⎝m⎠
GMm
f(xe) = xe(1 + xe) exp (–xe) with xe = ⋅ (3.4)
re kT
The function f(xe) is plotted in Fig. 3.6 for the relevant regime xe appreciably
greater than unity. In this regime, f(xe) is a steeply declining function of increasing
xe. With everything being equal, the Jeans escape mechanisms drastically favor the
escape of light molecules, since vT and f(xe) are both larger for smaller values of m
in Eq. (3.4). This gives the number of molecules lost per second. Hydrogen atoms
are particularly susceptible to evaporation from a planet.
Figure 3.6: The function f (xe) is a rapidly declining function of its argument for xe considerably greater than
unity.
The substitution of realistic numbers into Eq. (3.4) reveals that the Jeans
mechanism is inadequate and to see this, we consider the Earth. Let us adopt for
8
T the temperature 4,000 K, which results in re = 7.0 × 10 cm (about 600 km above
sea level). The collisional cross-section of hydrogen atoms by hydrogen atoms is
about 1015 cm2. To take into account collisions with other molecules, we adopt a
collisional cross-section per hydrogen atom equal to A = 6 × 10–14 cm2. With these
values we obtain from Eq. (3.3) a total loss rate of hydrogen atoms by Jeans escape
from the Earth of Nm = 1 × 1027 sec–1.
This value is comparable to observations of Lyman alpha emission from
the Earth’s upper atmosphere. During the 4.6 billion years of its existence, the
Earth would have lost at the above rate, about 2 × 1020 gm of hydrogen, roughly
THE ATMOSPHERES OF THE PL ANETS 31
three orders of magnitude less hydrogen than is contained in the form of H2O
in the oceans. The hydrogen loss rate measured for Venus is comparable to the
above rate. The Jeans mechanism thus seems incapable of giving a much larger
rate in the past. If Venus ever had as much liquid water as the Earth, it probably
could have lost only 10–3 of it by the Jeans mechanism. Hunten has shown that
any increase of the exospheric evaporation rate ultimately is limited by the ability
of the lower layers to re-supply the lost particles by diffusion. From the regions
where the material is well-mixed, the hydrogen or hydrogen-carrying compounds
has a limited capacity to rise diffusively relative to the other constituents in the
atmosphere.
EXERCISES
4.1 INTRODUCTION
T
he issue of the origin of the solar system is not the same as that of the origin
of the Universe. The Universe is sufficiently larger than the solar system. For
the issue of origins, geology supports that the solar system is appreciably
younger than the oldest objects in the Universe. The solar system thus has formed
against the backdrop of an already mature Universe. There are two popular
theories about formation of the solar system. The first is popularly known as the
solar nebula theory, originally proposed by Immanuel Kant during 1775 and later
in 1796 by Pierre-Simon Laplace. The second was originally proposed by Buffon,
taken up by Chamberlian and Moulton and developed by Jeans and Jeffreys and
by 1930 it became the theory of the origin of the solar system. This is known as the
planetsimal theory. However, these theories do not give proper evidence about the
age of the solar system or Earth. More reliable methods are to find the ages of the
solar system rocks by the study of the long-lived radioactive elements which they
contain.
33
34 ASTRONOMY AND ASTROPHYSICS
theory of the origin of the solar system, although serious objections against the
theory still remain today.
The geologists of the eighteenth and nineteenth centuries showed how the sequence
of rock strata and fossils could be arranged in a qualitative and semi-quantitative
sequence. However, absolute estimates of ages were stymied by the uncertain and
variable rates of weathering, volcanism, erosion, and sedimentation. Estimates of
the age of the Earth from theoretical considerations of its rate of cooling suffered
from lack of knowledge about the composition of the Earth’s interior, the origin
of its interior heat, and the role of solar heating for keeping the surface warm. By
the nineteenth century, it was clear that if the rate of geological processes in the
past was not different from that in the present, then the Earth must be much older
than roughly 6,000 years. Geologists insisted that the Earth had to be even older
than the 20 or 30 million years that Kelvin and Helmholtz deduced from their
contraction theory for the power of the Sun. Rutherford applied the definitive
method for determining the age of the Earth by radioactive dating.
Figure 4.1: An example of radioactive dating of uranium ores. The number of uranium-238 atoms N(238U)
has declined; its present value is given by the solid circle, by a known law of exponential decay.
Its ultimate decay product lead-206 has a number, N(206Pb), whose present value is given by the
open circle, can also be measured. The sum N(238U) + N(206Pb) must remain a constant when
extrapolated back in time. Therefore, the age of the ore cannot be older than a time when this
sum consists of a contribution only from 238U and nothing from 206Pb.
ORIGIN OF THE SUN, THE EARTH, AND THE SOLAR SYSTEM 37
By techniques described above, geochemists have found the ages of the oldest
rocks on Earth to be about 3.8 billion years. It is, however, believed that the Earth
itself is older than this. Rocks older than 3.8 billion years have apparently not
survived to this day, or else there were none that solidified in an early phase of the
evolution of the Earth. The radioactive dating of lunar rocks shows that the oldest
lunar rocks are about 4.2 billion years old. Scientists speculate that still older rocks
may lie on the other side of the Moon, where less destruction of lunar rocks by
meteoroid bombardment may have occurred. In contrast with the spread of ages
found for Earth and Moon rocks, the radioactive dating of meteorites show them
to be all the same age, to within the errors of measurement, between 4.5 and 4.7
billion years old. The interpretation has grown that the first objects solidified in the
solar system about 4.6 billion years ago. This is the probable age of the entire solar
system.
Figure 4.2: A sophisticated method for dating rocks that contain uranium. As in Fig. 4.1 we begin by
measuring the present amounts of 238U and 206Pb, given by the solid and open circles, respectively.
We also measure the amount of lead-204 (204Pb), which is not a product of radioactive decay,
and multiply that amount by R, the ratio [206Pb]/[204Pb] of lead-206 to lead-204 found in
non-radioactive rocks. The product R.N. (204Pb) is then assumed to be the amount of primitive
lead-206 incorporated into the original rock at its solidification.
The value of 4.6 billion years refers to the time since the asteroids solidified. However,
meteoroids are fragments of asteroids. There is another age for meteoroids; the
38 ASTRONOMY AND ASTROPHYSICS
elapsed time since the shattering that produced the meteoroid. This latter age is
possible to obtain as the solar system suffers a continuous bombardment by high-
energy, cosmic-ray particles. These high-energy particles can smash the heavy
atomic nuclei that exist in an iron meteorite, and produce lighter nuclei, some
of which may be radioactive. Before the asteroid was shattered, the particular
piece of meteorite, which eventually fell to Earth, may have been deep inside the
main body, and therefore shielded from the external cosmic-ray bombardment.
Sophisticated techniques exist to measure the exposure ages of meteoroids, based
on the idea that longer exposures lead to a greater concentration of long-lived
isotopes relative to short-lived ones, as the latter attains an equilibrium between
production and decay. Such measurements yield exposure ages of hundreds of
millions of years for the iron meteorites and tens of millions of years for the stones.
Discrepancies arise from the use of radioactive species with different half-lives.
Such apparent discrepancies originate in a variable flux of galactic cosmic rays
encountered by the solar system.
conclusion that we can draw is that the age of the radioactive elements is probably
somewhere between 7 and 15 billion years. In order of magnitude, this delineation
is in good agreement with astronomers’ estimates of the age of globular cluster
stars and the age of the Universe. The three independent methods of establishing
the ages of ancient objects produce the same number within a factor of two or
so, which is very reassuring (Fig. 4.3). We can state with considerable confidence
that the Sun and the solar system are probably less than half the age of the oldest
objects in the Universe.
Fig. 4.3 shows estimates of the ages of three ancient objects. The horizontal
bars indicate the spread of estimates of ages of the Universe, globular clusters, and
radioactive elements that most astronomers would consider still to be permissible
given the uncertainties of the measurements. Notice that there is significant overlap
in the 9 to 15 billion year range. Thus, the age of the solar system, 4.6 billion years,
is likely to be between one-half and one-quarter of the age of the oldest objects in
the Universe.
Since Kepler’s work in the early part of the seventeenth century, it has been
known that the planets follow his three general laws of motion, which are listed in
Table 4.1.
40 ASTRONOMY AND ASTROPHYSICS
1. The orbit of a planet forms an ellipse with the Sun at one focus.
2. The Sun-planet radius vector sweeps out equal areas in equal times.
The geometric meaning of Kepler’s three laws can be explained easily. The
first law states that the orbits of all planets with long orbital periods and short
orbital periods take the form of ellipses with the Sun at one focus. The second law
states that equal areas as indicated by the shaded regions in the bottom diagram
(Fig. 4.4) are swept out in equal times. The third law states that the square of the
2 3
period is proportional to the cube of the semi-major axis: τ ∝ a with a constant
of proportionality, which is the same for all planets. After Kepler, it was Newton
who coupled his laws of mechanics with his theory of gravity and explained the
Kepler’s empirical discoveries in terms of dynamics.
Figure 4.4: Kepler’s laws for planetary motions around the Sun.
This topic is, an unfinished story, since humans have a great potential to change
the fact of the Earth. One may continue here to the relatively early history of the
Earth, where the changes that occurred were in the atmosphere of the Earth or,
more generally, what is called the biosphere.
ORIGIN OF THE SUN, THE EARTH, AND THE SOLAR SYSTEM 41
Figure 4.5: World geography as it may be taught 50 million years from now.
It is now believed that the change of the Earth’s atmosphere — from a reducing
to an oxidizing one — was a two-stage process. The first step is inorganic and
depends on the production of ozone, while the second step is organic and depends
on the appearance of life. Water plays a vital role in both steps.
42 ASTRONOMY AND ASTROPHYSICS
Figure 4.6: Schematic formation of the ozone layer. UV light, which strikes water molecules in either the
primitive oceans or the atmosphere dissociates them into hydrogen and oxygen.
ORIGIN OF THE SUN, THE EARTH, AND THE SOLAR SYSTEM 43
The role of the ozone layer in filtering out the UV rays of the Sun cannot be
overemphasized. UV rays are very harmful to many lifeforms, inducing skin
cancer in humans. This explains the concern over our release into the atmosphere
of fluorocarbons, which have an enormous potential for destroying ozone. In
geological times, the dissociation of water by ultraviolet light was halted by the
appearance of the ozone shield. The full conversion of the Earth’s atmosphere
from a reducing one to an oxidizing one had to wait for the development of life.
Life played a major role in the evolution of the Earth’s atmosphere. Walcott was
the first to draw attention to thinly layered deposits of limestone rock called
stromatolites and interpreted these structures as fossilized reefs of angle growth.
The angles are now gone, but they have left behind ‘beds of stone.’ Paleontologists
remained skeptical until 1954, when Tyler and Barghoorn reported the finding of a
large number of fossil microorganisms in similar structures. Some of the stromatolite
deposits in Africa and Australia are known to date back to 3.5 billion years ago.
This suggests that life developed on Earth soon after the crust re-solidified after
the melting of the Earth. In 1980, Awramik and Schopf discovered microfossils of
long filamental chains of bacterial cells in rocks from the region of the Australian
stromatolites. Around 3.5 billion years ago, unicellular organisms existed on Earth
that could perform photosynthesis, a process by which organisms that contain
chlorophyll convert carbon dioxide and water into carbohydrates. Carbohydrates
are the basic source of chemical energy for all organisms. Schematically, for
example, the production of glucose, C6H12O6, by photosynthesis can be expressed
as
Chlorophyll
6CO2 + 6H2O + light energy ⎯⎯⎯⎯⎯⎯
→ C6H12O6 + 6O2.
The carbon dioxide on the left-hand side could have come from geological
activity when carbonate rocks are heated or from biological activity via respiration.
Respiration involves the burning of carbohydrates with oxygen to release chemical
energy, e.g.,
free oxygen was not substantial until two billion years ago. In rocks more than
two billion years old, geologists find reduced minerals, which implies that these
rocks were not exposed to much free oxygen when they formed, while in rocks
less than two billion years old, geologists find oxidized minerals. This means that
in the intervening period, within a few hundred million years, there was a major
shift of the Earth’s atmosphere from reducing to oxidizing. The shift is responsible
for all the iron in the oceans to rust and sink (i.e., to precipitate to the ocean floor).
When the iron buffer in the oceans had been swept clean, the oxygen released
could begin to rise to modern atmospheric concentrations (Fig. 4.7). The shift from
a reducing to an oxidizing atmosphere allowed a wild proliferation in varieties
of microorganisms. The large number of microfossils turned out to be amazingly
diverse. The organisms have shapes that include long filaments, umbrellas,
spheroids, and stars. Fossil stromatolites became quite common after about 2.3
billion years.
Figure 4.7: Schematic plot of free oxygen in the Earth’s atmosphere in the past as a fraction of the present
supply.
versatility than cells without nuclei. For one thing, two cells with the same nuclei
structure can serve entirely different specialized biological functions.
Figure 4.8: Examples of eukaryotic cells in rocks aged less than 1.5 billion years.
EXERCISES
1. What are the two popular theories about the formation of the solar system?
2. Explain clearly the ‘Solar Nebula Theory’ and the ‘Planetsimal Theory’.
3. How was the age of the Earth determined by geologists of the eighteenth and nineteenth
centuries?
4. How was the age of the Earth determined (by Rutherford) using radioactive dating?
5. Consider the radioactive dating of Uranium ores to its ultimate decay product of
lead, and then explain the radioactive dating technique for determining the age of the
Earth.
6. What is the age of the oldest rocks found on Earth by geochemists by applying a
sophisticated method of dating rocks?
7. What are meteoroids? What is the idea on which the exposure ages of meteoroids are
based?
8. Discuss how radioactive dating permits astrophysicists to measure the ages of
radioactive elements.
9. The age of the solar system is likely to be between one-half and one-quarter of the age
of the oldest objects in the universe. Please explain.
46 ASTRONOMY AND ASTROPHYSICS
10. Write down the three laws of Kepler for planetary motion and then explain them
briefly.
92U and 82Pb how many total α and β decays need to occur for a transform-
11. For 238 206
12. Does fluorescence arise from the excitation of the nucleus of an atom?
13. Explain why radium is used in watch dials.
14. Given the ratio of 235U to 238U in uranium ore is 0.007. Assume the formation by a single
event and calculate the elapsed time since the creation of uranium. The half lives of
235
U and 238U are ts = 7.1 × 108 year and ts = 4.6 × 109 respectively.
Chapter 5
THE STARS IN THE SKY
5.1 INTRODUCTION
A
star can be considered a huge sphere made of very hot gas mainly hydrogen
and helium, which generates its energy through nuclear fusion. The self-
gravitational pull of the star is so high that it keeps its atmosphere of hot
gases intact. A fraction of the energy generated in the star is radiated out. The stars
are categorized depending on size or mass (i.e., white dwarfs, main sequence,
giants, and super-giants).
It is believed that the Universe was created about 10 to 20 billion years ago with
the Big Bang (see Chapter 1). The cosmic expansion resulted in the spreading of
helium and hydrogen elements, which are the major constituents of a star. The
distribution of the gas varied with time and the expansion of space. In the regions
where the gas concentration was large, the mutual gravitational attractions of these
gas molecules led to the growth of the first generation of stars.9 These swirling
9
A prominent star-forming region is within the Orion Nebula and is called Orion KL. Most of the
star systems are separated by many light years.
47
48 ASTRONOMY AND ASTROPHYSICS
gaseous clouds gradually became massive, and began to collapse under self-
gravity, thereby forming a dense central region called a ‘protostar,’ which warmed
up due to compression. As more and more gases are pulled in the process, the
internal pressure mounts up, thereby igniting nuclear reactions in its core. Thus, a
star is born. The protostars generally eject gaseous streams in opposite directions.
Some of them generate highly focused or collimated outflows, while the rest
produce more oval-shaped winds. It is not clear what produces these outflows,
but astronomers believe that the magnetic fields and the growing accretion disk
around the newborn star play crucial roles. As more and more material is attracted,
the protostar spins faster and faster. The increasing rotational energy might direct
the stellar winds along the axis of rotation. The powerful magnetic fields might
further constrain this outflow thus resulting in bipolar jets extending over several
light years. The bipolar outflows disrupt the surrounding gases to a large extent.
The motion of these gases may be studied under radio spectrum to understand
their effect on the protostar. The heat released by the protostars in the form of
ultraviolet radiation breaks apart the surrounding molecular gases into atoms,
which are ejected in bipolar outflows. These outflowing gases cool down and
recombine into molecules, which emit millimeter waves in radio. These emissions
could be mapped and more information on the activity around the protostellar
region could be obtained.
It must be noted that the clouds of interstellar gas are very large and their
collapse does not originate just one star, but a set of stars called a star cluster after
having fragmented into smaller clouds. Each fragment may give rise to a single
star or a set of stars that orbit around a common barycenter. For example, in our
Milky Way Galaxy, single stars are nearly 50% of the total. The majority among the
remaining stars are gathered into binary systems, and in some cases, into multiple
systems consisting of even six stars.
Until the nuclear furnace in a star remains alive (i.e., until the hydrogen fuel
exists in it), the star keeps shining and is said to be a main sequence star. After
the fuel ends, the outer layers of the star can no longer resist the central force of
gravity, which is pulling the star’s outer matter inward toward its core. The length
of a star’s life also depends on its mass. While some of the stars may explode
violently through a supernova and rapidly create even heavier elements and spew
much of the stellar material into space, in other cases, the process could be slower,
and instead of an explosion the elements from the star’s interior zones rises to
the surface and get lost in space when the outer layers blow off. In either case the
end results are similar (i.e., the space between the stars gets enriched with heavy
elements, many of which condense to form small solid grains). These processes
of the birth and death of stars occur over and over again, with each successive
generation of stars starting off with a greater quantity of heavy elements than the
previous generation.
THE STARS IN THE SKY 49
⎛ ⎞
2 h ν3 ⎜ 1 ⎟
I(ν, T) = ⎜ hν ⎟ (5.1.1)
c 2 ⎜ kT ⎟
⎝ e − 1⎠
⎛ ⎞
2 h c2 ⎜ 1 ⎟
I(ν, T) = 5 ⎜ ν ⎟ (5.1.2)
λ ⎜ λhkT ⎟
⎝ e − 1 ⎠
The quantity I is called the spectral radiance or spectral intensity, which is the
energy per unit time per unit surface area per unit solid angle per unit frequency
and is measured in W m–2 sr–2 Hz–1, ν is the frequency of radiation measured in hertz,
T is the temperature of the black body in K, h is the Planck’s constant (6.6260693
× 10–34 Joules/Hz), c is the velocity of light (2.998 × 108 m/sec), e is the base of
natural logarithm (2.718282), and k is the Boltzmann constant (1.3806505 × 10–23
joule/kelvin). If the black body is in thermal equilibrium with the surroundings,
then Kirchoff’s law states that the emissivity of the body is equal to its absorptivity.
The radiated energy increases rapidly with increases in temperature and can be
expressed by Steffan-Boltzmann law as
F = σT 4 W/m2. (5.2)
Here, F is called the energy flux or power per unit radiated from an object and
σ = 5.7 × 10–8 W m–2 K–4 is the Steffan-Boltzmann constant for all wavelengths.
The total energy radiated from an object per unit time is defined as the luminosity.
The surface area of a star can be calculated from its radius R as 4πR2 and thus the
luminosity L of the star can be expressed as
L = 4 π R2 σ T 4 W. (5.3)
Eq. (5.3) clearly indicates that the luminosity of a star is directly proportional to
its surface area and to the square of the square of its temperature. The wavelength at
which the emission spectra of black body is maximum, i.e., λpeak is also dependent
on its temperature and is governed by Wien’s displacement law as expressed in
Eq. (5.4).
2.9 × 10 −3
λpeak = m. (5.4)
T
50 ASTRONOMY AND ASTROPHYSICS
When viewed from a distance like from Earth, the observed flux of an object,
also known as its apparent brightness, is the power that we actually receive from
it. This quantity is generally measured in Wm–2 and is dependent on the distance
of the object. If r is the distance between the object and point of measurement, then
apparent brightness or observed flux f can be expressed as
L
f = ⋅ (5.5)
4 π r2
The brightness of the star is sometimes expressed logarithmically and called
magnitude. These could be (i) apparent or (ii) absolute, respectively, denoted by m
and M. The apparent magnitude m of a star is a measure of its apparent brightness
as seen by an observer on Earth. Since the amount of light received actually depends
on the thickness of the Earth’s atmosphere coming in the line of sight to the object,
the apparent magnitudes are normalized to the value it would have outside the
atmosphere. The apparent magnitude mx in the spectral band x is defined as
mx = –2.5 log10 (fx) + C (5.6)
Here, fx is the flux at the spectral band x and C is a constant, which depends
on the units of the flux and the band. Thus, the dimmer an object appears, the
higher its apparent magnitude. On the other hand, absolute magnitude refers to
a measurement at a fixed distance of 10 parsecs. That is, if the luminous object
is placed at a fixed distance of 10 parsecs from Earth, the magnitude measured
would be its absolute magnitude. If the actual distance of the object is known, then
the absolute magnitude M is expressed as
Mx = mx – 5(log10 DL – 1). (5.7)
Here, DL is the luminosity distance in parsecs. For the nearby objects like
the stars in our galaxy, the luminosity distance DL is more or less identical to the
real distance, since the space-time is almost Euclidean. For much more distant
objects, this approximation is invalid and the general relativity must be applied for
calculating. Table 5.1 lists the apparent and absolute magnitudes of some stars.
Although in a true sense a star is not strictly a black body, it may be approxi-
mated with it within some range of frequency. The reason being is that the star
contains various elements, which produce absorption lines at different wavelengths
depending on the elemental composition of its photosphere and chromospheres.
Thus, if the luminosity of a star is plotted as a function of wavelengths, absorbing
type glitches would be found on the curve. Moreover, the star also produces sudden
active emissions such as radio bursts, etc., whose intensity are much higher than
when the star is quiet. If the glitches in the spectrum produced out of absorption
lines in a star are neglected, then within a small range of spectrum the quiet stars
may be approximated to a black body as shown in Fig. 5.1.
THE STARS IN THE SKY 51
Figure 5.1: Approximate black body models of stars at different temperatures using Planck’s law. Note that
the wavelength at peak emission is highly temperature dependent and is governed by Wien’s
displacement law.
A star like our Sun may be approximated to a black body when the glitches
due to spectral absorption are neglected and when it is quiet for within a range
of 10–7 to 10–4 m wavelengths as shown in Fig. 5.2. The Sun produces an extra
amount of electromagnetic radiations when certain activities in it are alive. These
Figure 5.2: The Sun is approximated with a black body at 5800 K. This approximation is reasonable within
the range of the wavelength between 10–7 to 10–4 meters.
52 ASTRONOMY AND ASTROPHYSICS
confer to the large bursts in the gamma-rays and flares in upper ultraviolet and
x-rays, microwaves, and radio waves. The shaded regions in the diagram show
the most probable ranges of intensity fluctuation of these special radiations, which
vary with time and are unpredictable. However, at certain times, the amplitude
range could exceed.
Table 5.2: Morgan-Keenan classification of stars based upon temperature and Spectra.
Before proceeding further, we shall discuss a bit more about the features of
stars based on their class. The strength of the hydrogen line is related to the surface
temperature of the star. The spectrum of stars are highly size dependent. The super-
giant stars often swing between O or B (blue) and K or M (red). But they do not
stay in the G classification which is an extremely unstable place for a super-giant.
The classes are further subdivided between 0 and 9. For example, G0 is the hottest
and G9 is the coolest within the class-G. The Sun is a G2 star.
(1) Class-O: The O stars are bluish and very hot and luminous with most of
their emission in the ultraviolet range. They constitute only 1 out of 32,000
main sequence stars and shine a million times brighter than our Sun. They
produce prominent helium lines (ionized and neutral) and weak hydrogen
lines. Due to their huge size, their fuel gets consumed quickly thereby
leaving the main sequence stars. Examples of O stars are Zeta Puppis and
Epsilon Orionis.
(2) Class-B: The B stars are blue and extremely luminous consisting of neutral
helium and moderate hydrogen lines in their spectra. The B stars also live
as a main sequence only for a short time like the O stars due to their huge
size and large fuel consumption rate. They constitute about 1 out of 770
stars (0.13%) in the main sequence. Owing to their short life span they
do not have sufficient time to stray far from the region of their births.
10
Thus, they tend to form clusters known as OB1 associations. These are
associated with giant molecular clouds. Our Milky Way Galaxy has one of
10
These are gravitationally bound groups of stars. Two types of clusters known as globular and
open. The globular clusters constitute hundreds of thousands of very old stars. The open clusters
generally constitute less than a few hundred stars and are generally very young.
THE STARS IN THE SKY 55
its entire spiral arm as OB1 association known as Orion OB1 association.
Examples of B stars are Rigel and Spica.
(3) Class-A: These are white or bluish-white stars commonly seen with the
naked eye. Their spectra consists of strong hydrogen lines and ionized
metals. They constitute 1 out of 158 stars (0.63%) in the main sequence.
Examples of A stars are Vega and Sirius.
(4) Class-F: These are white stars with a slight tinge of yellow. Their spectra
consists of weak hydrogen lines and ionized metals. These constitute 1 out
of 32 stars (3.1%) in the main sequence. Examples of F stars are Canopus
and Procyon.
(5) Class-G: These are yellow or light-yellow stars. Their spectra contain weak
hydrogen lines compared to the F stars along with metals (ionized and
neutral). These constitute 1 out of 12 or 13 stars (8%) in the main sequence.
Examples of G stars are Sun and Capella.
(6) Class-K: These are orangish or light orangish stars. Their temperatures
are slightly less than our Sun. Their spectra consists of extremely weak
hydrogen lines or nil along with metals (mostly neutral). These constitute
of roughly 1 out of 8 stars (13%) in the main sequence. Examples of K stars
are Arcturus and Aldebaran. Some of these could be giants or super-giants
and may not be a main sequence star.
(7) Class-M: These are reddish or reddish orange stars. Their spectra consists
of lines belonging to molecules and neutral metals, but the hydrogen is
usually absent. More than 78% of the stars in the main sequence are red
dwarfs. Examples of M stars are Betelgeuse and Barnard’s star. Some of
the stars could be super-giants such as Antares and Betelgeuse.
Additionally, a number of new classes have been identified (also called
extended classes) from the newly discovered ones which are not described here.
Figure 5.3: The Hertzsprung-Russel diagram showing the spectal classification of the stars. Only a few stars
are highlighted by their names including our Sun. Stars falling on the angled lines take the radii
indicated on the line. In the main sequence the lifespan of the stars are marked.
An examination of the HRD reveals that the stars tend to fall within certain
regions of the diagram. The most predominant is the diagonal belt called the main
sequence. The upper left regions hold the stars with very high temperature in the
blue region. Note that the masses are very high. In the central region lies our Sun.
The bottom right region of the main sequence consist of smaller stars with mass
less than our Sun and are relatively much cooler.
The luminosity of a star is dependent on its mass and its current state of existence.
The star passes through different states of existence like protostar, main sequence,
red-dwarfs, super-giants, white-dwarf etc. These are very much dependent on
their mass. The luminosity also varies accordingly. For example, a star like our
Sun while it was forming, probably it was 100 to 200 times brighter than now. In
its first billion years in the main sequence it was probably 2/3 times brighter than
present. At the end of the main sequence its brightness would probably reduce to
50% of today following which it would probably swell up to become a Red Giant
with its diameter exceeding the orbit of Earth and its brightness increasing about
1000 times of today’s brightness. After this it would probably shrink to become a
THE STARS IN THE SKY 57
white-dwarf 1000 to 10,000 times fainter than today. In general stars have a huge
variation of luminosities, from red-dwarfs which are 50,000 times fainter than the
Sun to super-giants which are 1,000,000 times brighter than the Sun. But if we
restrict our study of brightness for the main sequence alone, we would probably
be able to quantify the mass luminosity relationship.
Figure 5.4: The mass luminosity relationship data for various stars in the main sequence.
Fig. 5.4 shows the luminosities of main sequence stars possessing different
masses. It may be observed that luminosity steadily increases with mass in a non-
uniform manner. In other words, the luminosity is proportional to some power
of mass. Further the variations are different for mass above and below 0.43 MΘ.
Based on these the luminosities are related to the masses as follows
3.5
⎛ M ⎞
L = LΘ ⎜ ⎟ for M ≥ 0.43 MΘ (5.8)
⎝ MΘ ⎠
2.3
⎛ M ⎞
L = LΘ ⎜ ⎟ for M < 0.43 MΘ (5.9)
⎝ MΘ ⎠
58 ASTRONOMY AND ASTROPHYSICS
In all the above we have seen that the luminosity of the star is dependent on its
temperature which in turn is dependent on the mass of the star. Thus there should
be a minimum mass below which the temperature of the star’s core11 will not be at
sufficient temperature to conduct the nuclear reactions. This is known to be about
0.08 MΘ or 80 times than the mass of Jupiter. Similarly the limits on the maximum
mass of a star based on theoretical modeling is thought to be somewhere between
150 and 200 MΘ. These very much short lived and rarely found.
EXERCISES
11
It is the innermost region of a star where the energy is generated through nuclear fusion.
THE STARS IN THE SKY 59
17. Calculate the numerical value of luminosity of the Sun from the following data: r = 1.5
× 1013 cm and f = 1.36 × 106 erg sec–1 cm–2.
18. Assume the luminosity of the Sun is 3.9 × 1033 erg sec–1; find the linear size of the Sun
and its effective temperature. Given: r = 1.5 × 1013 cm and the angular diameter of
the Sun θ = 32 minutes of arc.
19. Prove that the total energy of the Earth’s orbital motion is equal to one-half of its
potential energy.
20. Establish that the numerical value for the Sun’s central temperature is equal to
1.5 × 107 K.
21. Calculate the ratio of plasma-energy density to the radiation energy density at the
center of the Sun.
60 ASTRONOMY AND ASTROPHYSICS
Chapter 6
THE SUN AS A STAR
6.1 INTRODUCTION
I
t is the high quality of sunshine, which maintains life on Earth. The quantity
or sunshine is equally important for maintaining the comfortable average
temperature of about 290 K on the surface of the Earth. A happy combination
of quality and quantity of sunshine makes Earth a thriving garden.
To get an idea about how high-grade sunshine is, let us start with directly
measurable quantities and deduce the luminosity and effective temperature of
the Sun. From direct measurements, the energy which falls per unit time per unit
area perpendicularly outside the Earth’s atmosphere (i.e., the solar constant) is
obtained as
61
62 ASTRONOMY AND ASTROPHYSICS
33 –1
LΘ = 3.90 × 10 erg sec . (6.3)
The value of LΘ represents an enormous total power. The Earth intercepts less
than one billionth of this power, and even this is mostly wasted by immediate re-
radiation in the infra-red back into space.
large burst on the Sun. The total energy contained at these wavelengths is small,
in the quiet Sun. Near the peak of the Sun’s emission a 5,800 K, black body gives
a representation of most of the output of the Sun. Even at optical wavelengths,
Fig. 6.1 shows detectable departures in the spectral-energy distribution of the
continuum radiation, at higher dispersions dark absorption lines appear in the
solar spectrum. These clues suggest that the Sun is not at a uniform temperature,
and that its surface layers are not in perfect thermo-dynamic equilibrium with its
radiation field. The temperature 5,800 K characterizes only part of the surface layers.
Since the zeroth law of thermodynamics states that heat flows from hot to cold,
we may infer that the Sun is hotter on the inside. This can be checked observationally
by the phenomenon known as limb darkening. The edge of the Sun, called the limb
in astronomical nomenclature, appears darker than its center. The explanation is
given in Fig. 6.2, which shows two lines of sight for an observer; one directed
toward the limb of the Sun and the other toward its center. For the same optical
path lengths traversed through the outer layers of the Sun, the eye looks closer
to point C along a line of sight that passes through the Sun’s center than along
the Sun’s limb. As the deeper layers are hotter, and intrinsically brighter, the limb
seems dark relative to the center. This is the origin of limb darkening. The evidence of
absorption lines can be interpreted to give the same conclusion. We assume that the
atmosphere of the Sun is characterized by cooler rarefied regions, which overlie
hotter dense regions (Fig. 6.3). Thermal photons flow from the deeper layers with a
continuous distribution of wavelengths. Some of them with the proper wavelengths
will be absorbed by atoms situated in the upper layers (Fig. 6.3 b). The excited atom
may then de-excite in two different ways as shown in (Fig. 6.3 c).
Figure 6.2: The explanation for limb darkening in the Sun. For the same optical path lengths traversed
through the outer layers of the Sun (the short solid lines), the eye sees deeper and hotter layers
along a line of sight aimed through the Sun’s center than along one aimed through the Sun’s
limb.
(i) It may de-excite collisionally with the excess energy carried off as extra
kinetic energy of the two colliding partners. Since the temperature of the
gas in the upper layers is lower than the temperature of the radiation
64 ASTRONOMY AND ASTROPHYSICS
field emerging from the lower layers, statistically there are more
radiative excitations followed by collisional de-excitations than there are
collisional excitations followed by radiative de-excitations. So, there is
more absorption than emission. The net effect is the removal of some line
photons from the underlying continuum radiation field. This mechanism
for producing the intensity at the wavelength of the line which is darker
than the intensity at the wavelengths of the neighboring continuum is
known as true absorption.
Figure 6.3: Formation of a dark line in a stellar photosphere. (a) A continuum of photons with different
wavelengths flow from high temperatures to low temperatures toward the surface. (b) The
absorption of a line photon causes an atom to enter an excited state. (c) The excited atom
subsequently de-excites either collisionally (true absorption) or radiatively by the emission of a
photon in a different direction (scattering).
THE SUN AS A STAR 65
(ii) The radiatively excited atom may also de-excite radiatively, with the
emission of a line photon traveling in a different direction than the original
photon. This process of resonant scattering does not constitute true absorption.
It can also lead to the appearance of a dark line. The statistical condition
of detailed balance is equivalent to the condition of true thermodynamic
equilibrium. However, in a stellar atmosphere, the presence of a boundary
prevents detailed balance. There exists a basic asymmetry at the surface
of a star. A hot, bright star lies to one side while cold, dark space lies to
the other. The net loss of line photons from our line of sight compared to
the neighboring continuum photons leads to a dark scattering line. Dark
lines are produced generally by a combination of true absorption and pure
scattering. The quantitative study of the formation of the continuum
radiation and the pattern of absorption that lies in the outer layers of a
star is part of the discipline called stellar atmosphere. The layers that can
be analyzed in this way only constitute the outer 0.1 percent of the radius
of the Sun. These layers, where optical photons can fly to us more-or-less
directly, are collectively called the photosphere of the Sun.
In order to discover the conditions in the deep interior of the Sun astronomers
have to resort to theory as there are no direct observations of the Sun’s interior.
This theory is complex, though the main qualitative ideas are simple.
For getting the mass of the Sun we know that the Earth orbits the Sun
approximately in a circle with radius r = 1.5 × 1013 cm and period P = 1 year =
3.16 × 107 sec. Therefore, the speed of the Earth’s motion about the Sun is
circumference 2 π r
v = = = 2.98 × 106 arc sec–1. (6.4)
period P
If we ignore the motion of the Sun, we have F = ma in the form where a is the
2
centripetal acceleration (magnitude = v /r) and F is the gravitational force (F =
2
G MΘ /r ) with m = mass of the Earth. Solving we get the mass of the Sun MΘ as
rv 2
MΘ = = 2.0 × 1033. (6.5)
G
The Sun is thus the most massive body in the solar system. It is about 300,000
times more massive than the Earth. The Sun is about 100 times larger than the
Earth. So, the mean density of the Sun is not very different from the mean density
of the Earth. The mean density of the Sun is
66 ASTRONOMY AND ASTROPHYSICS
mass MΘ
mean desnity of the sun = = = 1.4 gm/cm3. (6.6)
volume 4 π RΘ3
3
This is to be compared with the average density of the Earth, 5.5 gm/cm .
To estimate the pressure, we consider the observation that the Sun is not
shrinking or expanding rapidly. For the Sun to remain mechanically static, all
forces must be very nearly in balance, the condition of hydrostatic equilibrium.
This hydrostatic equilibrium involves the balance of gravitational force and
pressure force at each level in the star (Fig. 6.4). In hydrostatic equilibrium, the
pressure at radius r equals the weight of the vertical column of matter of unit area
above it. In this way, when the forces are balanced at all values of r, no net force
acts on any element of gas, and the star remains mechanically steady.
Figure 6.4: Hydrostatic equilibrium involves the balance of gravitational force and pressure force at each
level in the star.
The high mass of the Sun indicates that the pressure needed at the center to
support the overlying layers against its own weight must be enormous. In such
high-pressure regions, neutral atoms and molecules would be pushed so hard
against one another that their electronic shells could not possibly survive intact.
Thus, the matter in the interior of the Sun must be composed of bare nuclei and
free electrons, a state of matter called plasma. Although such a mixture of ions
THE SUN AS A STAR 67
and electrons is electrically neutral on average, the sizes of the particles are much
smaller than atoms or molecules. Consequently, even at the relatively high mean
density of the interior of the Sun, there exists plenty of room between the individual
particles (Fig. 6.5) so we may say that the plasma in the Sun behaves as a perfect
gas.
The pressure of a perfect gas is associated with the random thermal motions
of the particles. If the gas has a number density n and a temperature T, the perfect
gas law states that the total pressure
P = nKT. (6.7)
Figure 6.5: The solar plasma contains the bare nuclei of atoms (open circles) plus their stripped orbital
electrons (black dots). The mean separation between particles is much larger than their
individual sizes and the medium is electrically neutral as a whole. The electric energy associated
with neighbors is small in comparison to the thermal energy. Thus, the solar plasma can be
considered a perfect gas.
The central mass density ρc in the Sun is about 110 times its average value 1.4
gm/cm3 (i.e., ρc = 150 gm/cm3). The mean particle mass mc at the center is less than
the mass mp of the proton: mc =1.5 × 10–24 gm. Thus, the central number density
nc = ρc/mc = 1.0 × 1026 cm–3 with the pressure P = 2.1 × 1017 erg cm–3. Eq. (6.7)
implies that the Sun’s central temperature must equal 15 million degrees Kelvin.
average before they crash into a material particle and are either absorbed or
scattered into a completely different direction.
years. So there must be some other, much larger supply of energy that explains the
Sun’s power.
Although thermonuclear energy is a very powerful source, it does not have
infinite capacity to offset radiation from the surface of a star. Thus, sooner or later,
every star gets hotter and hotter while losing more and more energy to a cold
surface. This confrontation provides the central plot to the life stores of the stars.
Figure 6.6: Radiative transfer in a star in the stellar interior. Photons suffer many interactions with matter
and slowly diffuse outward by a tortuous series of ‘random walks,’ which eventually bring them
to the photosphere. In the process, one x-ray photon is dragged into many optical photons.
One source of energy is gravitational energy. The basic idea is based on the
virial theorem. The Sun is a statistically steady object. The total energy E of the Sun
is numerically equal to half its gravitational potential energy. In the Sun, it is the
thermal energy of the plasma rather than the kinetic energy of orbital motion that
produces the resistance to gravitational contraction. By the virial theorem the total
energy E, which is the sum of the thermal energy and the gravitational potential
energy, is equal to half the gravitational potential energy, i.e.,
Figure 6.7: A self-gravitating star differs from a cooling ember in one very important respect.
The Sun has a perfectly sharp outer edge, coinciding approximately with the
location of its photosphere. The total thickness of the photosphere occupies only
0.1 percent of the radius of the Sun. Therefore, the Sun looks as if it has a sharp
edge, which explains the conventional definition of the photosphere of the Sun as
its surface. There is gas belonging to the Sun, which lies above the photosphere,
in the chromosphere and the corona. The chromosphere itself is a transition zone
between the properties of the gas in the photosphere and those of the gas in the
corona (Fig. 6.9). The gas temperature as we go outward from the photosphere
initially decreases, and later it begins to climb dramatically, until temperatures
around 2 million degrees Kelvin are reached in the corona. This temperature is
much higher than the photospheric temperature of around 5,800 K. However, the
light output from the photosphere vastly dominates that from the corona. This is
because the corona gas is so rarefied that it has little emissive power. The gas in the
chromosphere and the corona is optically thin. For this reason, it is difficult to see
the corona in optical light unless the light from the photosphere is blocked, either
naturally by an eclipse of the Sun by the Moon, or artificially by a coronagraph.
When the light from the solar photosphere is blocked by a solar eclipse, light
arises from the chromosphere wherein hot gas lies above cool gas. This situation
gives rise to emission lines instead of absorption lines. This is the origin of the
flash spectrum associated with emission lines in the chromosphere observed during
solar eclipses. That element is helium. Helium is quite rare on Earth, but it is now
known to be the second most abundant element in the Universe. Helium is the
second simplest chemical element. This is a second important clue to the creation
of the Universe. The 2 million degree temperature of the solar corona means that
this gas can emit x-rays. This explains the presence of x-rays in the solar spectrum.
X-ray photographs of the Sun look quite different from the usual optical pictures.
More information about the structure and dynamics of the solar corona is obtained
from satellite x-ray observations of the solar corona.
Astronomers believe that solar corona is heated by acoustic waves and other
kinds of mechanical waves generated in the solar convection zone lying below the
solar photosphere. These waves propagate upward into the more rarefied corona
where they steepen into shock waves and dissipate their energy into heat. This
heats the corona, which has very little heat capacity, to such high temperatures
that the gravitational field of the Sun finds it hard to hold onto the gas in the outer
parts of the solar corona. Theoretically, it was shown that these portions blow
continuously from the Sun to form a solar wind, which pervades the solar system
(Fig. 6.10). It is now possible to collect the particles in this wind by satellites, and
to verify that the outer parts of the Sun are made primarily of hydrogen.
Figure 6.10: Comparison of the magnetic fields of a bar magnet and that produced by an electric current
circulating through the coil of a wire.
THE SUN AS A STAR 73
In addition to the solar wind, the solar atmosphere also manifests other violent
activity in the form of flares, loops, arches, etc. This activity may provide a
substantial input of energy and particles into the expanding solar corona. The
active regions are associated with the presence of magnetic fields in the surface
layers of the Sun. Magnetic field activity can be studied qualitatively by observing
sunspots, or by measuring the Zeeman effect splitting in atomic lines.
Fig. 6.10 makes comparison of the magnetic fields of a bar magnet and that
produced by an electric current circulating through the coil of a wire. The magnetic
field lines never end on an object. Magnetic field lines could end on or originate
from an object only if particles called magnetic monopoles existed. Such particles
have not been found in nature.
Fig. 6.11 shows a schematic energy versus a magnetization curve for a bar
of ferromagnetic material. At high temperatures, the state of minimum total
energy corresponds to zero magnetization. At low temperatures the state of zero
magnetization is a local maximum of total energy. There are two minima of total
energy, which correspond to spontaneous magnetization with the north pole of the
magnetic dipole at the left or right end of the bar. The actual bar magnet will settle
on one of these two states as denoted by dots. The phenomenon of spontaneous
magnetization is another example of our general rule for the thermodynamic
behavior of matter. At high temperatures, the state of least energy corresponds
to more freedom for the atomic magnets to twist and turn. At low temperatures,
ferromagnetic materials prefer to gain more magnetic binding energy by giving up
some freedom for the atomic magnets to twist and turn.
Figure 6.11: Schematic energy versus magnetization curve for a bar of ferromagnetic material.
74 ASTRONOMY AND ASTROPHYSICS
Fig. 6.12 shows the magnetic-field structure of the Sun, as deduced from
ground-based observations. As seen it is usually quite chaotic. The average
strength of the magnetic field on the surface of the Sun is comparable to that on
the surface of the Earth. In localized spots, however, the solar magnetic field can
rise to thousands of times higher than the average value.
Figure 6.12: Magnetic field structure of the Sun, as deduced from ground-based observations.
The present radius of the Sun has the right value to maintain a central temperature
to provide a rate of nuclear energy generation that exactly balances the leakage
rate due to the random walk of the photons from center to surface. A summary of
the important concepts of stellar structure is given in Table 6.1.
2. Energy Transfer. Photons in the interior of the star carry heat outward by random
walking from regions of high temperature to regions of low temperature. At the
photosphere of the star, the photons make a transition from walking to flying. If the
luminosity required to be carried out is too large for radiative transfer to handle
stably, convection results.
Fig. 6.13 shows schematically the structure of the Sun from core to photosphere
as understood from detailed modeling. The solar core is defined as that region
where the temperatures are high enough to yield appreciable rates of hydrogen
fusion into helium. The solar photosphere is defined as the outer layers where the
photons stop walking (i.e., optically thick) and begin flying (i.e., optically thin). In the
sub-photospheric layers lies the solar convection zone. The Sun generates nuclear
energy in its core at a rate equal to LΘ = 3.9 × 1033 ergs/sec. This luminosity is
ultimately radiated from the photosphere according to the formula
L = 4 π R2 σ Te .
4
(6.11)
Here, L is LΘ and R is RΘ for the Sun.
EXERCISES
1. Define the solar constant. What is the value of the solar constant obtained from direct
measurement?
2. How can you determine the effective temperature of the Sun?
3. What is the origin of “limb darkening” of the Sun? Explain with a suitable diagram.
4. Explain the formation of a dark line in a stellar photosphere.
5. How were the conditions in the deep interior of the Sun discovered?
6. How are the mass, mean density and pressure inside the Sun estimated?
7. What do you mean by solar plasma? How can the solar plasma be considered as a
perfect gas?
8. Explain how the Sun’s safety valve for thermonuclear fusion works.
9. “There is a net slow effect performed by photons inside stars like the Sun.” Considering the
dance of the photons as “random walk,” explain the statement.
10. How does a self-gravitating star differ from a cooling ember?
11. Discuss elaborately the two outer layers- ‘Chromosphere’ and ‘Corona’ of the Sun.
12. Critically discuss the magnetic field structure of the sun as deduced from ground-based
observations.
13. Compare the magnetic fields of a bar magnet and that produced by an electric current
circulating through the coil of a wire.
14. Give the summary of the important concepts of “stellar.”
15. Show schematically the structure of the Sun from core to photosphere, as understood
from detailed modeling and discuss.
16. A typical amateur’s telescope is given to you. It has a 1 meter objective focal length
and a 40 mm eyepiece. What would be its magnification? Does the magnification apply
only to observation of double stars?
17. Find the diffraction limit in sec of an arc of a 5-meter telescope corresponding to a
wavelength of 5000 Å. Compare the calculated value to the seeing limit.
Chapter 7
ASTRONOMICAL TELESCOPES
7.1 INTRODUCTION
A
stronomical telescopes can be constructed using either reflection or
refraction to gather light over a large collecting surface and thereby to
focus it to a much smaller area. This property of gathering more light
and focusing it, not magnification, is the primary purpose of most astronomical
telescopes.
7.2 REFLECTION
If a thin layer of metal is deposited onto a glass surface and is polished, a mirror
is formed. The reflecting property of the mirror is a result of the microscopic fact
that metals contain many free electrons, which are not bound to the atoms in
the lattice of the metal. These free electrons respond readily to the presence of
electromagnetic waves. Also, they allow almost no radiation of sufficiently long
wavelengths to pass through the metal. A mirror, therefore, reflects almost all the
optical light that falls on its surface. The geometric properties of a flat mirror can
be deduced from the knowledge of simple wave properties of light. Light reflects
from a flat surface in such a way that the angle of reflectance is equal to the angle
77
78 ASTRONOMY AND ASTROPHYSICS
Figure 7.2: An elliptical mirror brings the light from a point source at one focus to the other focus.
Let us now assume that the point source of light is at infinity. If the two foci of
an ellipse become infinitely separated the ellipse becomes a parabola. Therefore,
to bring the light rays from a very distant point source like a star to a single focus,
the objective of a reflecting telescope should be a paraboloidal mirror. To collect
as much light as possible, we know that the aperture of the paraboloid should be
made large.
Figure 7.3: Part of a complete ellipsoidal mirror will capture only part of the light that leaves a point
source.
80 ASTRONOMY AND ASTROPHYSICS
Figure 7.4: An x-ray telescope, which works by grazing incidence reflections. The focus is brought to the
back of the mirror arrangement.
One disadvantage of a reflector is that its prime focus is in front of the mirror.
Consequently, for an astronomer to observe at the prime focus, an appreciable
fraction of the incoming light may be blocked (Fig. 7.5 a). This problem can,
however, be partially solved in several ways. The first solution was invented
by Newton. In that arrangement (Fig. 7.5 b) a secondary mirror is used to focus
the light at the side of the telescope tube. This arrangement tends to make the
system lopsided. A more satisfactory solution was devised by Cassegrain. In
the Cassegrain arrangement (Fig. 7.5 c) a secondary mirror is used to focus the
light behind the telescope tube through a hole in the primary mirror. However,
none of these arrangements are very suitable for making precise spectroscopic
observations. This is because the bulky and heavy spectrographs would strain on
the telescope tube if they were suspended in front, to the side, or in the back of
the primary mirror. For such studies the Coude arrangement was devised where
a series of secondary mirrors bring the light down the polar axis of the telescope
drive to a fixed focus in a room below the telescope (Fig. 7.5 d). In the Coude
room, large pieces of equipment can be mounted properly in a stationary
configuration.
ASTRONOMICAL TELESCOPES 81
Figure 7.5: Various arrangements for reflecting telescopes: (a) The prime focus in front of the primary
mirror. (b) The Newtonian focus at the side of the telescope. (c) The Cassegrain focus at the
back of the telescope. (d) The Coude focus underneath the telescope.
7.4 REFRACTION
We know that the speed of light in air va is greater than the speed of light in
glass vg. A whole collection of photons are incident as a plane wave on an air-
glass interface at an oblique angle (Fig. 7.6). The photons are imagined to come
from a distant source so that they arrive in waves whose fronts are connected by
the dashed lines. The initial direction of the propagation of the photons and the
wave-front is represented by the heavy top arrow. As the photons pass through the
glass interior, they are effectively slowed down and travel less distance in a given
interval of time. The net effect is to turn the direction of the wave-front, so that the
propagation of the wave-front in the glass is in a direction closer to the normal.
This turning is a property of the entire collection of photons. An individual photon
cannot be suddenly turned in direction upon entering the glass.
In a refracting telescope the main optical element is the objective lens. It collects
light over the entire area of the aperture and forms an image of the observed object
at the focus. In order to discuss the focusing ability of a lens we imagine a source
like a star, which is a single point of light. Since the star is at infinity, the rays from
the star arrive at Earth nearly on parallel lines. A perfect lens would bring these
parallel rays to a single point at the focus. This point is the correct image for our
idealized star. In practice, starlight is made up of light of different colors. Since
blue light is bent more by a glass lens than red light, chromatic aberration results
(Fig. 7.8). A photographic plate placed at the blue focus will exhibit a blurred red
image, and vice versa at the red focus. We ignore this chromatic aberration to discuss
light-gathering power. A lens can be used in reverse of the above description; i.e.,
if an object (or the image of a real object) were placed at the focus of a lens, the
rays from it, in passing through the lens, would leave in parallel lines (Fig. 7.9).
A human eyeball placed behind this lens would perceive the object or the image
as if it were at infinity. Thus, the combination of a large objective lens plus an
eyepiece (Fig. 7.10) serves as an astronomical telescope (i.e., refractors) and this
enables an astronomer to gather starlight over a much larger area than his or
her own eyeball could. Electronic devices (i.e., image intensifiers) with or without
photographic film allow the astronomer to observe to much fainter limits than is
possible with naked-eye observations, even on telescopes with very large apertures.
Different apparent sizes of the final image, for example, the apparent separation
of a double star result from using different eyepieces. The light-gathering power
is found by the square of the diameter of the objective, not by the eyepiece. For an
extended object, spreading the same amount of light over a larger area decreases
the apparent surface brightness of the observed object. Thus, high magnification
may be useful for solar-system observations, but it is not required for other types
of astronomical investigations.
We assume that there is another star lying at an angle θo with respect to the first.
The rays from both the stars will pass through the center of the objective lens on a
straight line. These rays will be brought into focus at a distance equal to its focal
length fo behind the objective. If the eyepiece is placed at its focal length fe further
back, the double-star image will pass through the eyepiece and be collimated into
two parallel beams, at an angle θe, to each other. Behind the eyepiece an astronomer
will perceive an inverted image of two stars separated by an angle θe, whereas the
actual double star is only separated by the angle θo. The angular separation thus
has been magnified by the ratio θe/θo. To calculate this magnification in terms of the
properties of the telescope, as indicated in Fig. 7.11, we may consider the relation
Radio waves do not suffer from this atmospheric blurring. The angular
resolution of modern radio telescopes is limited by the diffraction produced by the
wave nature of light. We know that the ray-tracing method to determine the focusing
of electromagnetic waves assumes that the wavelength is too small compared with
the aperture of the telescope. However, on angular scales comparable to the ratio
of the wavelength λ of light to the diameter D of the telescopes collecting surface,
diffraction effects produce an unavoidable blurring even of point sources. This
is due to the fact that the secondary waves produced by the interaction of light
and matter at the edges of the objective interfere to give a pattern of light and
dark fringes which, on small angular scales, is different from that predicted by
ray tracing. For a circular aperture, the blurring produced by diffraction limits the
angular resolution by an amount
86 ASTRONOMY AND ASTROPHYSICS
Figure 7.13: Principle of aperture synthesis illustrated by an imaginary line of radio telescopes at the North
Pole. During one 12-hour period, the rotation of the Earth causes this line to sweep out in a
full circle.
ASTRONOMICAL TELESCOPES 87
EXERCISES
8.1 INTRODUCTION
T
here are some primary reasons for observing the Sun (e.g., the Sun is a star
closest to us, it is a variable star, and we are more or less dependent on its
behavior). This chapter provides the important facts about the behavior of
the Sun and describes how interesting observations may be made. This is covered
in two sub-sections: (i) observation of the Sun in integrated (white) light and
(ii) observation in monochromatic light (primarily Hα). The discussion includes
certain general features such as sunspots, groups of spots, and photospheric
activity, flares, prominences, and chromospheric activity. Specific observational
techniques suitable for individual or group programmers are described, which
may be undertaken for the observer’s personal pleasure, or with scientific purposes
in mind.
The Sun is a star that varies over the course of time. It is the closest star, and
the light we receive from it was emitted only 8 minutes earlier. The largest modern
telescopes can detect light from distant celestial objects, which was emitted
thousands of million of years ago. The fact that the Sun is so near provides us with
a considerable amount of light, enabling us to make detailed measurements with
short exposures without having to use large apertures. It is particularly favorable
for amateur observation. In fact, it is useful to have a long focal length because
89
90 ASTRONOMY AND ASTROPHYSICS
the diameter of the primary solar image is about 1/100th of an objective’s focal
length.
For the Sun, it is the quality of the observing site that ultimately determines
the largest aperture to be used. We have learned that solar observations may be
carried out with small diameter telescopes and that 150 to 200 mm is adequate.
The main problem is to reject the excess light, and thereby to reduce the heating of
the interior of the instrument and to protect the eyes of the observer.
Even a small amateur telescope may be used to examine the surface of our daytime
star and to study some of its features. But one should know the dangers posed by
direct observation of the Sun. Let us summarize the precautions to be taken to
avoid any possibility of damage to the eyes.
Permanent changes to the retina due to the unprotected observation of the
Sun, particularly during eclipses, have been analyzed by Quintel for amateur
astronomers. When an observation is taken using any magnifying optical system,
the amount of light entering the eye is proportional to the surface area of the
instrument’s objective. For a luminous body like the Sun, the flux received on the
retina is considerably increased by the optical system, but the eye can protect itself
particularly by reducing the diameter of the iris to a minimum, which is about
1.5 mm for most subjects. If the light is not too strong, one is only dazzled, and no
lesion occurs. As soon as the bright light is removed, the eye regains its sense of
vision after a few minutes. On the other hand, if the light is intense, the amount
of heat absorbed becomes more than the cooling capacity of the blood supply to
the retina can handle. There will be an increase of temperature of the cells, leading
to definite damage, but the person will not feel any sensation of pain that might
serve as a warning. The severity of the burns to the retina depends on the flux
received and the length of the exposure. It is probable that a weak exposure
without sufficient protection may cause the death of some cells in the center of the
retina. It is therefore essential to ensure that the eyes are properly protected even if
it seems unlikely that one will be dazzled. The fact of being dazzled, which leads
to a decrease in visual response, may, paradoxically, give the impression that the
observation is comfortable.
Therefore, one should never put his or her eye to the eyepiece before having
checked to ensure that a suitable filter or other device has been fitted to the
instrument. The second is to cover the objective of the finder to avoid the risk
of burns. In addition one should be wary of the solar filters provided with most
commercial telescopes, which are located at the eyepiece end, and absorb the light
HOW TO OBSERVE THE SUN 91
at the point where it is closest to the focus of the objective. Such a filter should
always be used behind an objective filter. Otherwise, it may shatter in a few seconds
through the intense heating to which it is subject. The whole of the luminous flux
may suddenly reach the retina and cause serious damage before the observer has
time to move away from the eyepiece. Moreover, it is to be ensured that the filter
absorbs over a sufficiently wide spectral band. Some filters pass a considerable
amount of infra-red or near-ultraviolet radiation, which is invisible but may still
cause burning of the retina.
Here, we describe three methods of observing the Sun: direct visual observation,
observation by projection, and photographic observation, which are explained
below.
the solar image formed on a projection screen. The speed of this process, when the
setting is carried out in advance, means that drift may be eliminated and a drive
is no longer indispensable. This method demands considerable care in aligning
and setting the telescope properly to get a truly circular image. The method is
extremely sensitive to atmospheric turbulence.
Figure 8.1: An example of a narrow bandwidth obtained by combining the transmission curves of the filter
(a Schott RG 610) and the spectral sensitivity of the film (Kodak Panatomic-X).
94 ASTRONOMY AND ASTROPHYSICS
Figure 8.3: Defining the events to find the orientation of the observation of an observation made with an
Altazimuth instrument.
The path followed by the Sun across the sky during the day is the circular arc
R.M.S. If the angle is between the geographical North/South direction (i.e., the
meridian NPASP) and the vertical then it can be shown that
Fig. 8.4 shows how the North/South line rotates with respect to the vertical
during a day. The plane of the diagram is the plane tangent to the celestial sphere
at the center of the solar disk. The rotation of the North/South line takes place in
the same direction as the hands of a clock.
Figure 8.4: Variation in the geographical North/South direction as a function of time of day for an
Altazimuth instrument.
Figure 8.5: Determination of the direction of the Sun’s axis of rotation with respect to geographical
North/South line.
HOW TO OBSERVE THE SUN 97
the visible region and the radio window at longer wavelengths. These two spectral
regions include the wavelengths where most of the energy is radiated from the
photosphere into the lower corona. The chromosphere is almost transparent and
is visible against the disk. After the light has passed through the photosphere and
chromosphere, its spectrum exhibits a large number of absorption lines called
the Fraunhofer spectrum having highly variable widths, relative intensities, and
profiles. Measurement and analysis of their precise wavelengths throughout the
spectrum enable us to find the Sun’s chemical composition. Solar absorption lines
contain more than 25,000 different lines in the region accessible to photography
alone. Identification is made by comparison with laboratory spectra of atoms and
molecules and about 80 different elements have been recognized in this way, the
most predominant being hydrogen, calcium, and helium. The strongest lines are
the H and K lines of ionized calcium.
8.7 SPECTRO-HELIOGRAPHS
In the design of a spectrograph an objective forms an image of the object and a slit
selects a narrow region of this image. The light from this strip is dispersed by a
prism or grating preceded by a collimating lens and followed by an appropriate
camera lens. The dispersion varies proportionately to the wavelength and it takes
place in a direction perpendicular to the axis of the slit. This produces a spectrum
at a specific wavelength of the area covered by the entry slit. If one spectral line is
isolated by a second slit, the image falling on the latter represents the area of the
Sun selected by the first slit at the wavelength selected by the second slit. Thus, the
result obtained is a monochromatic image. The image of the Sun may be trailed
slowly across the entry slit and at a result it sweeps the whole of the solar image at
constant speed. If the film behind the second slit is moved simultaneously in the
same direction and at the same speed, it captures an image of the solar surface at
the desired wavelength. Such a device is called spectro-heliograph (Fig. 8.7).
Solar spectro-heliographs are heavy instruments and cannot be mounted on a
refractor. The spectro-heliograph is generally a fixed instrument and light is fed to
it via a coelostat. The following are the sources of error with spectro-heliographs:
100 ASTRONOMY AND ASTROPHYSICS
(i) Permanent. If the displacements of the slits are not precisely parallel, or
if the first and second slits have unequal residual curvature the errors are
permanent.
(ii) Random. If there is a variation in the rate at which the slits move relative
to the fixed parts of the equipment causing wear in any part of the variable
speed drive train, periodic error arises from the lead-screws and the clutch
for the worm and nut, or any other periodic error arises from inaccuracies
in any part of the whole drive train.
8.7.1 Coelostat
The coelostat has two mirrors. The first one called the primary rotates at half the
Sun’s diurnal rate around an axis parallel to the Earth’s axis and located in the same
plane as the reflecting surface; consequently, it reflects light into a fixed direction.
A secondary mirror, on the other hand, subsequently reflects the beam toward the
optical axis of the objective at the entrance to the spectroscope. In order to allow
for the Sun’s changes in declination, the secondary mirror may be shifted in a
North/South direction so that it intercepts the beam from the primary mirror. In
the winter season when the Sun is low, the secondary is very close to the primary
and may even cast a shadow on it. When this occurs, the primary is shifted toward
HOW TO OBSERVE THE SUN 101
the East in the morning and toward the West in the afternoon. The field rotation
that results needs a correction to the Sun’s position angle to be made.
If undertaken on an amateur scale, a portable coelostat may permit the
less robust part of the instrument to be permanently mounted under cover. But
alignment should be done with great care since the coelostat may introduce errors,
particularly if observation includes long exposures. These are as follows:
(i) Permanent. Errors owing to the setting of the hour angle if observations are
not always carried out at the same time of day and errors due to the poor
mechanical state of the tangent screw in the drive strain are considered
permanent. Another error is that of variations in the speed of rotation
of the hour axis due to the variations in the frequency of the electricity
supply.
(ii) Random. An error in alignment is not highly significant in summer when
the primary mirror coincides with the optical axis of the instrument. The
same does not apply in the winter, when the coelostat is East or West of
the optical axis. Any error is responsible for the field rotation and thus
changes the position angle.
8.8 HELIOGRAPHS
EXERCISES
1. What are the dangers posed by direct observations of the Sun? What are the precautions
to be taken to avoid any possibility of damage to the eyes?
2. Describe the method of direct visual observation of the Sun. What are two reasons for
which this method is not recommended?
3. How is the Sun observed by projection? Describe the method.
4. Why is the direct photography of the solar image (at the primary focus of the instrument)
not wise? What is a filter and where it is to be placed?
104 ASTRONOMY AND ASTROPHYSICS
5. Why is high contrast not the most important factor in solar photography? How is the
geographical North-South determined using the photographic method?
6. Discuss the mounting and orientation of the image.
7. Explain with a suitable diagram the variation in the geographical North-South direction
as a function of time of day for an Altazimuth instrument.
8. How are sunspots watched by amateur astronomers?
9. How is the Wolf Number related to the number of individual spots? How is the active
area calculated?
10. How do you measure the proper motion of sunspots?
11. What are the instruments developed to get monochromatic images of the Sun? Give the
underlying principle.
12. Describe the construction of spectro-heliographs.
Chapter 9
OBSERVING THE SUN
WITH A CORONAGRAPH
9.1 INTRODUCTION
T
he corona is only about one millionth of the brightness of the Sun, so
its observation with an ordinary telescope, outside of eclipses, is quite
impossible. It was realized about a century ago that the corona is the
Sun’s outer atmosphere. Many astronomers tried to study it outside eclipses.
Total eclipses are rare and visible from very restricted areas of the Earth. Initially,
astronomers succeeded only to observe prominences, at the very base of the
corona. It observed at wavelengths corresponding to their emission line, that they
look much brighter. The instrument that allows observation of the corona outside
eclipses is the coronagraph. Invented by Bernard Lyot in 1930, it succeeded in
eliminating the source of scattered, stray light, which is stronger than the corona if
no precautions are taken.
Observing white-light corona using a coronagraph is a difficult operation. It
requires both a perfect instrument and an ideal site, but it is still possible to observe
monochromatic emission from the corona under slightly less perfect observing
conditions with a suitable monochromatic filter. This requirement became less
105
106 ASTRONOMY AND ASTROPHYSICS
costly with the development of the interference filter around 1950. Subsequently,
a small-sized coronagraph suitable for amateurs was devised at the beginning of
the 1960s. Twenty-five years later the study of the corona has taken on new life
with the development of radio methods and through ultraviolet observations from
artificial satellites.
Figure 9.1: An arrangement often initially considered to observe the Sun’s outer atmosphere.
OBSERVING THE SUN WITH A CORONAGRAPH 107
Figure 9.2: Using this laboratory arrangement, the various forms of stray light that occur within the lens L1
become visible.
We now turn to the brilliant spot of light in the center of the image of the lens
L1 that is observed on the screen E. Fig. 9.3 shows a cross-section of the objective
L1, which is a simple lens, normally chosen to be plano-convex. The major part of
the light incident on the lens will form an images S′ of the light source, and this
108 ASTRONOMY AND ASTROPHYSICS
will be hidden behind the disk D. But transmission of a beam of light through a
lens is not perfect, and partial reflections arise at the entry and exit surfaces. In the
objective, a small fraction of light will be reflected by the outermost surface and
will be sent back toward the Sun. This will, of course, not cause any problems. An
approximately equal amount of light will be reflected at the glass-air exit surface
of Fig. 9.3. This portion will suffer a second partial reflection at the entry surface,
and this will send a small fraction of the beam back to form a secondary image of
the Sun, somewhere near S′. This image is very faint as it is formed by a double
reflection, and its intensity is approximately 5/100 × 5/100 = 25/10000 of that of
the principal beam. But it is particularly inconvenient as it should be formed at
a point where it cannot be hidden by the opaque disk D. The bright spot seen in
the center of the image of L1 (Fig. 9.2) is caused by this light that has undergone a
double reflection from the surfaces of the objective.
Figure 9.3: The geometrical paths of rays doubly reflected by the surfaces of a simple lens.
given by the objective will be a slightly widened, bright line, bordered by weaker
fringes with minima where the intensity does not completely drop to zero. Finally,
we consider the image of a luminous area that is sharply bounded at its right-
hand edge (Fig. 9.4, bottom). We may consider the luminous area as consisting of
numerous narrow, luminous lines side by side. The image given by the objective
cannot be properly represented in a drawing, but it will gradually diminish. This
means that a certain amount of light will reach areas of the image that ought to be
completely black.
The disk of the Sun resembles a luminous spot, the brightness of which falls
to zero as soon as we pass the limb. Although the Sun may be like that, because
of the diffraction effects, the Sun’s image will not be sharply defined but will fall
off gradually. This stray light will not be blocked by the opaque disk inside the
coronagraph; rather, it will be mingled with the weak light of the corona we want
to investigate.
Figure 9.4: Due to diffraction phenomena, even a perfect lens does not produce an infinitely small
luminous point, but instead produces the well known pattern (top) of a disk surrounded by
rings. Similarly, the image of an infinitely fine line is enlarged and bordered by figures (center).
Further, the image of a bright area with a sharp edge will appear as a more gradual decline
(bottom). This causes considerable degradation of observations at the solar limb.
distinctly smaller diaphragm (about 2/3 or 3/4 of the size of L1’s image) is used. It
is wise to lose a little light rather than to run the risk that any slight misalignment
will allow part of the bright fringe.
Fig. 9.5 shows the design of the coronagraph as devised by Lyot. The objective L1
is a simple plano-convex lens, perfectly polished, with a focal ratio of about 1/20.
Due to this long ratio, geometrical aberrations are negligible. In addition, problems
of heating at the primary focus may be controlled without much complication.
The metal disk D is calculated to be the same size as, or slightly larger than the
image of the solar disk that it is to hide. The apparent diameter of the Sun varies
during the course of the year as the distance of the Earth from the Sun varies. So a
different disk must be provided to correspond to the size at different times, and to
cater to the various wavelengths being studied. This second point arises because
the objective L1, being a simple lens, suffers from chromatic aberration.
Fig. 9.6 reveals a method of mounting the reflecting cone and occulting the
disk D that has proven successful in practice. The two are held by a central rod
several centimeters in front of the lens L2, which has been drilled through the center.
An insulating sleeve protects the glass from the heat transmitted by the cone. To
avoid thermal problems, the cone should be highly polished, or even aluminized.
Lens L2 forms an image of L1 at the diaphragm E. Fig. 9.7 shows schematically the
image that should be visible at E when the coronagraph is properly adjusted. The
diaphragm should have a diameter somewhat less than the full aperture of L1.
This means that some light is lost, but more stray light will be lost than light from
the corona, which is fundamental. The resulting arrangement is also reasonably
insensitive to small errors in alignment. The central disk should not be in contact
with any optical component.
Behind the diaphragm E in Fig. 9.5, all the instrumental sources of scattered
light have been eliminated. It only forms an image of the eclipsed Sun onto a
recording medium, e.g., onto a film F, by a photographic objective L3. To avoid
OBSERVING THE SUN WITH A CORONAGRAPH 111
Figure 9.6: Method of mounting the cone and occulting disk D in front of lens L2 in a coronagraph.
The bold line in Fig. 9.8 shows the variation in intensity from the center of the solar
disk to the white-light corona. The scale here is logarithmic, where the lower of the
fine lines shows the brightness of the sky close to the Sun during a total eclipse.
This is about 1/1000th of the intensity of the inner corona. On the right-hand side of
the figure, the brightness of a blue sky, several lens of degrees away from the Sun,
is indicated. This light characteristic of a clear blue sky arises from the scattering
of sunlight by the molecules forming the Earth’s atmosphere. It is seen from the
figure that its brightness is more or less equal to that of the white-light corona. The
112 ASTRONOMY AND ASTROPHYSICS
sky close to the Sun is not blue; the disk of the Sun behind a distant object is a very
brilliant, white halo around the object. This is indeed an atmospheric effect and
its contribution is shown in Fig. 9.8 (upper line). This is, in fact, an atmospheric
corona, but the term halo is used to differentiate it from the solar corona.
Figure 9.7: Appearance at the plane of the Lyot diaphragm in a properly adjusted coronagraph.
Figure 9.8: Ordinate shows the luminous intensity on a logarithmic scale with respect to the intensity at the
center of the solar disk. While the abscissa reveals distances from the center of the Sun, expressed
in solar radii, the curves show the variation in the intensity of the Sun’s outer atmosphere
(the white-light corona) and that of the sky through which it is observed both during and
outside an eclipse.
The halo, due to scattering by aerosols, is not the only detrimental effect
arising in the Earth’s atmosphere. Fig. 9.9 shows (solid line) the brightness of
the atmospheric halo as observed at a high-altitude site with a good quality
coronagraph that does not add any scattered light. The ordinate gives the
luminosity expressed on a logarithmic scale with respect to the center of the solar
OBSERVING THE SUN WITH A CORONAGRAPH 113
disk. The abscissa represents the distance from the Sun’s limb in seconds and
minutes of arc, also on a logarithmic scale. The solid line shows the variation in the
scattered corona at the final focus in a properly adjusted coronagraph. The dashed
line shows how this scattering changes between good seeing conditions and those
with a sky that causes a slight degree of scattering. The objective of a coronagraph
is generally a simple lens (Fig. 9.10), so it is impossible to block out the whole solar
disk at all wavelengths. Here, the occulting disk is placed at the focal plane for red
light (R), and it can be seen how part of the blue image (B) is not obstructed.
Figure 9.9: Ordinate shows the luminosity expressed on a logarithmic scale with respect to the center of
the solar disk, while the abscissa represents the distance from the Sun’s limb in seconds and
minutes of arc, also on a logarithmic scale.
red-light image of the Sun, there will be a flux of non-occulted blue light within
the instrument. In the best case (Fig. 9.9) the scattered halo at 1 arc-minute from the
solar limb has an intensity of 1/100000th of the solar disk. Because the white light
corona is about 1/1000000th of the brightness of the Sun, it is difficult to detect
without employing special methods.
Figure 9.10: The objective of a coronagraph is generally a simple lens. Here, the occulting disk is kept at the
focal plane for red light (R). Notice how part of the blue image (B) is not obstructed.
Fig. 9.11 shows a tracing of a spectrum of the corona. Apart from the continuum,
there are narrow but strong emission lines in the spectrum. The principal one is in
the green at 530.3 nm; another is visible in the red at 637.5 nm. At the wavelength of
an emission line, the corona is as much as 100 times stronger than the continuum.
The diffraction halo arising in the Earth’s atmosphere typically shows a continuum
with narrow absorption lines. Although the corona is weaker than the scattered
halo, it predominates if it is possible to make observations at the exact wavelength
of emission lines.
Figure 9.11: The solid line is a schematic representation of part of the coronal spectrum. In the presence of
scattered light, it is easier to see the corona if it is possible to isolate wavelengths close to an
emission line using a filter.
are far from being sufficiently selective (curve G in Fig. 9.12). On the other hand,
interference filters give maximum reduction of the scattered halo (curve 1 in
Fig. 9.12). It should be noted that it is difficult to use filters of bandwidths less than
1 nm because narrowband filters only give the desired bandwidth at a specific
temperature and for a truly parallel beam of light. To detect coronal emission it
is essential to reduce the scattered light as much as possible, so a filter with a
bandwidth of 0.2 nm must be used. This will require the use of a temperature-
controlled enclosure. The filter should be installed in a beam that converges or
diverges as little as possible. A solution (Fig. 9.13) is to use two doublets instead of
the lens L3 and to mount the interference filter between them.
Figure 9.12: General shape of the bandwidth of a filter made from colored gelatin (G) and of an interference
filter (I).
Figure 9.13: Interference filter should be located where the beam of light is essentially parallel. This is
achieved by splitting lens L3 as indicated.
shows a schematic representation of the most intense emission lines in the visible
region. Observations may be made easily if the Hα hydrogen-line can be isolated
by using a filter centered on 656.3 nm. A bandwidth of 5 to 10 nm suffices to show
prominences clearly, although a narrower filter allows observation when there is
greater scattering from the sky.
Figure 9.14: In the visible region the spectra from prominences mainly consists of hydrogen lines (Hα, Hβ,
Hγ, etc.) and a helium line (D3). The red Hα line is the strongest, and so it is easiest to see.
In all solar investigations, observers must not forget that precautions are to be taken
so that the intense flux may not damage the instrument, or worse, the observer’s
eye. These warnings are valid for observations with a coronagraph. Because the
blinding solar image is normally hidden by the occulting disk, it is tempting to
look at the prominences without any additional attenuation of the beam of light.
But any accidental decent ring of the coronagraph by an involuntary movement, a
gust of wind, or failure of the drive motor, etc., would uncover the Sun and could
immediately cause a burn, which might be severe. This only serves to emphasize
that special attention should be taken to the equatorial drive of a coronagraph.
A good drive is, of course, essential if good photographs are to be taken. Here
the potential problem is more serious. If someone is observing a planet through
a telescope that is tracking badly, there will be a slow shift of the object in the
eyepiece and this is not too inconvenient for visual observation or short exposures.
With a coronagraph, any fault in tracking will cause the Sun to appear in the form
of a dazzling crescent, preventing any useful observation, and possibly putting
the observer’s eyesight in danger. So the golden rule for the observer is not to put
an eye to the eyepiece unless the reduction in the beam of light is such that the
non-occulted image of the Sun can be born without any danger. If an interference
filter is used, then the reduction in the amount of sunlight is achieved by virtue
OBSERVING THE SUN WITH A CORONAGRAPH 117
The curve shown in Fig. 9.9 shows the variation in the brightness of the halo of
stray light observed at the final focus of a good coronagraph. This is repeated in
Fig. 9.15, which shows the variation in brightness of prominences observed in Hα,
Hβ, and D3. In Fig. 9.15, the relative intensities of the stray light and the emissions
correspond to those observed through narrow-pass interference filters. There is
no doubt that observation of the corona in monochromatic light is not simple
and should only be attempted by an experienced observer having the benefit of
Figure 9.15: The intensity of scattered light compared with the intensities of prominence and coronal
emission lines observed through an interference filter.
118 ASTRONOMY AND ASTROPHYSICS
an excellent site. We know that observation of the white-light corona when the
atmospheric corona cannot be reduced by proper filtration is scarcely feasible.
The images seen in monochromatic emission only show coronal material at a very
specific temperature, depending on the wavelength studied, so the visible features
entirely depend on the filter employed. The white-light corona exhibits material
over a wide range of temperatures simultaneously and thus gives a broader-based
picture.
The best way to become familiar with a corona graph is to observe prominences
in Hα light. It may be noted that on some days the Sun is practically featureless.
Prominences are constantly varying features in the solar atmosphere. Some of their
characteristics are described below.
Figure 9.16: Due to the rotation of the Sun, prominences may be seen with a coronagraph at the eastern
limb, then observed against the disk with spectro-heliograph, and finally become visible at the
western limb before disappearing behind the Sun.
OBSERVING THE SUN WITH A CORONAGRAPH 119
Figure 9.17: Bands of latitudes where prominences most frequently appear (shown by the shaded zones in
the figure) vary over the course of the 11-year cycle.
very high prominence does not remain visible for more than about 2 to 3 days each
time it crosses the limb. The lifetimes of prominences may exceed a month. Hence,
it is possible to see a prominence disappear over the western limb and then recover
it at the same latitude at the eastern limb 14 days later after the Sun has completed
half a rotation. Despite this, the appearance is very different after a period of two
weeks. The third time-scale of changes in prominences arises from the 11-year
solar cycle. Prominences are less frequent around the maximum, when the disk
has most spots. Moreover, the latitudes at which prominences appear depend on
the phase in the cycle. If the latitudes are plotted throughout the cycle, we may get
a diagram shown in Fig. 9.17. At a minimum there are two bands of prominences
on the Sun, at about 45° latitude in each hemisphere. Prominences may therefore
be seen crossing the limb in the Northeast, Northwest, Southeast, and Southwest.
As the cycle progress, some of the prominences shift toward higher and higher
latitudes forming what is known as the polar cap and ultimately more or less
reaching the poles at the time of maximum. It is then possible to see prominences
rotating in one spot over a considerable period of time. A second zone corresponds
to the latitudes at which spots are seen (i.e., starting at about 40° at the beginning
of the cycle and ending at about 10° just before the minimum). The pattern shown
in Fig. 9.17 permits one to predict reasonably well which regions of the Sun are
likely to show prominences.
Figure 9.18: The model illustrates that the relatively cold, dense material in prominences is held in a sort of
magnetic trap by the lines of force.
Anyone viewing the Sun through a coronagraph should not expect to get similar
events. Two main categories of events deserve to be monitored. When an active
OBSERVING THE SUN WITH A CORONAGRAPH 121
area crosses the solar limb it shows small, active prominences, which exhibit
rapid changes, often in the form of mass ejections of greater or lesser frequency.
In specialized works these are usually referred to as surges, sprays, etc. When
an area’s activity is most intense, one may see a chromospheric flare. Through a
coronagraph one may observe, very close to the Sun, arches and large and small
clumps of gas that appear very bright. Some hours later, these may have turned
into spectacular loops, which look like fountains of light, with material appearing
in the corona and raining down onto the surface of the Sun. Another category of
event is the sudden disappearance of a quiescent prominence. This may happen to
prominences at any latitude, but most frequently affects high prominences that are
close to active areas. Through the coronagraph, the prominence may be seen rising
to considerable heights and progressively dissipating while some of the material
appears to fall back on each side.
Figure 9.19: Optical diagram of a coronagraph used at the Pic du Midi to photograph the solar corona.
Details of one of the coronagraphs used at the Pic du Midi for photographing
monochromatic emission from the solar corona are presented here. Fig. 9.19 shows
a schematic diagram of the various elements. Table 9.1 shows the corresponding
various parameter settings of the coronagraph.
Final image of the Sun. This is formed about 360 mm behind L3″ on the film, which
is carried in an Exakta 24 × 36 mm body. The diameter of the (occulted) solar disk is
18 mm, with a longitudinal adjustment of L3″. Focusing is carried out by longitudinal
adjustment of L3″.
EXERCISES
5. Using a diagram, show the design of the coronagraph as devised by Lyot. Next
describe the different sections pointing out a method of mounting the reflecting cone,
and occulting the disk in front of the lens in a coronagraph.
6. Show the variation in intensity from the center of the solar disk to the white-light
corona, and explain.
7. Discuss how the stray light is removed owing to chromatic aberration.
8. Draw a coronal spectrum and explain the role of a filter.
9. Explain how one can detect the coronal emission successfully.
10. How is the interference filter located where the beam of light is essentially parallel?
11. What are the precautions used to safeguard the observer’s eyes when studying with a
coronagraph?
12. Draw the curves showing the variation in the brightness of the halo, if stray light is
observed at the final focus of a good coronagraph, and discuss the results.
13. What is the best way of becoming familiar with a coronagraph?
14. Discuss some interesting characteristics of the telescopic appearance and
prominences.
15. Briefly discuss the changes in the appearance of prominences over three typical time-
scales.
16. Explain the modes to illustrate the relatively cold, dense material in prominences.
17. What are the active and quiescent prominences?
18. Present the details of the coronagraph used at the Pic du Midi for photographing
monochromatic emission from the solar corona.
124 ASTRONOMY AND ASTROPHYSICS
Chapter 10
SOLAR ECLIPSES AND
STUDYING THE CORONA
10.1 INTRODUCTION
T
otal solar eclipses are impressive astronomical events that suddenly reveal
the solar corona. During totality the black disk of the Moon is ringed by
bright pink light and surrounded by a silvery halo. The phenomenon enables
extremely valuable astrophysical observations to be made. The solar corona is the
source of the solar wind and strongly scattered light within our solar system. Total
eclipses of the Sun allow thorough investigations of the corona.
Fig. 10.1 shows the Earth and the Moon in space illuminated by the rays of the
Sun. The Moon completes one rotation around the Earth in 27.3 days and has
an orbit inclined at an appreciable angle (5°09’) to the plane of the Earth’s orbit
around the Sun. The Earth-Moon distance varies around the lunar orbit and the
apparent size of the lunar disk is similar to that of the Sun (about 32 minutes of
125
126 ASTRONOMY AND ASTROPHYSICS
arc). This happy coincidence arises from the fact that the ratio of the distance from
the Earth to the actual diameter has a similar value for each of the two bodies.
When the apparent size of the Moon is larger than that of the Sun, there is a total
eclipse. If the apparent size is smaller, the eclipse will be annular. The duration of
totality is governed by the ratio of the apparent diameters of the two bodies, which
determine the magnitude of the eclipse. When the Earth’s rotation is taken into
account this duration never exceeds 8 minutes. These few minutes are extremely
precious for observers, as they offer the finest view of the corona. Anyone within
th
the shadow cone observes the brightness of the daytime sky fall to about 1/10000
of its normal value. Prominences may sometimes be seen around the lunar disk.
Surrounding everything is the magnificent halo of light scattered by the plasma
in the solar corona and by particles in the interplanetary medium. The conditions
required for any eclipse to be total are: (i) it must be a New Moon; (ii) the Moon
must be near the line of nodes; and (iii) the apparent size of the lunar disk must be
larger than the apparent size of the Sun.
Figure 10.1: A schematic diagram showing the conditions under which the Moon’s shadow may be
intercepted by the Earth, producing a total eclipse of the Sun.
Fig. 10.2 shows the various phases of a total eclipse, assuming the apparent disks of
the Sun (filled) and of the Moon (open) are perfectly circular. The disk of the Moon
crosses in front of the disk of the Sun, so the latter is completely covered between
the 2nd and 3rd contacts. This lasts for several seconds because the lunar limb is
very irregular (Fig. 10.3). Variations in the profile of the Maria and Highlands may
permit a certain amount of the solar photosphere to be seen, giving rise to Baily’s
SOLAR ECLIPSES AND STUDYING THE CORONA 127
Beads. The difference between the apparent motions of the two bodies across
the sky gives a typical relative velocity of 1 second of arc in 2 seconds of time.
When variations in the profile are considered, it is seen that the duration of the
visibility of Baily’s Beads may be taken into account. This is one of the difficulties
of measuring the duration of eclipses.
Published ephemeredes for the Sun and the Moon permit the average
apparent diameters of the two disks as seen from the Earth to be determined
very accurately from the topo-centric Earth-Moon and Earth-Sun distances. By
making the assumption that the diameter of the Moon remains constant over a
time-scale of a few centuries and by measuring the precise interval between the
2nd and 3rd contacts, it is possible to determine if the diameter of the Sun is truly
constant at least from one eclipse to the next. There are some doubts about the
constancy of the solar diameter. Provided allowance can be made for the profile of
the lunar limb, a precise measurement is possible. The intervals to be measured are
generally of the order of a few minutes and accuracy of the order of 0.02 seconds
of time is expected. This corresponds to measuring the Sun’s diameter to 0.01 arc-
seconds. If the observer is not on the central line, it may or may not be necessary to
introduce a correction for the relative positions of the disks. Using the calculated
ephemeredes, the interval of time may be converted into a value for the solar
diameter. The methods of carrying out this experiment are similar to those used
for studying the form of the Moon and for verifying ephemeredes.
The times of contact for a specific place may be calculated from ephemeredes, which
have a limited accuracy. This produces errors in the times or contact that are of the
order of a few seconds of time. One of the sources of error lies in our insufficient
knowledge of the form of the Moon, and also of the diameter of the Sun. Several
other factors also affect the calculations. As a result, the times of contact calculated
from the two sets of ephemeredes may differ by two seconds of time. Contact
times may be measured to a high degree of accuracy in a simple experiment. A
photometer was required in the past to measure the solar luminosity during the
128 ASTRONOMY AND ASTROPHYSICS
few seconds before and after totality. It is also possible to use a microcomputer
to provide an adequate time-base. In recent years it has become possible to do
far better, either by using a cine camera synchronized with the clock, or a video
camera. This method has an advantage over the photometer method, where the
total flux from the solar crescent is measured. It allows the problem posed by the
irregular edge of the Moon (Fig. 10.3) to be tackled by providing a visual record
of the Baily’s Beads produced by gaps in the Moon’s profile. To get a good-quality
Figure 10.3: Typical profile of the lunar limb as calculated for the total eclipse of November 23, 1984. Note
that the radial scale is greatly exaggerated (about 60 times).
The observation of the corona during totality should be made with the greatest
care so that experiments function without fail. It is preferable to undertake reliable
experiments than to attempt some special feat. Even if the resulting photograph
SOLAR ECLIPSES AND STUDYING THE CORONA 129
Figure 10.4: A schematic diagram of an experiment designed to measure the 2nd and 3rd contacts of an
eclipse.
is unique it cannot be used unless properly calibrated and all the instrumental
factors are well known. Fig. 10.5 shows the origin of the sky illumination during an
eclipse. A represents primary scattering by the ground and the lower atmosphere
of sunlight, which then illuminates B, the line of sight in the Earth’s atmosphere.
This further scatters some of this radiation toward the observer at C.
Fig. 10.6 shows why the solar corona may be observed during totality. The
sky becomes sufficiently dark for the contrast between the various features to be
perceptible. Even the faint Earthshine on the face of the Moon may be measured at
that time.
Figure 10.6: The radial distribution of intensity close to the Sun of the daytime sky and of various objects
observed during an eclipse.
super-spicules, etc., and very hot ones like coronal arches and flares. These extremes
are separated by only a few hundred kilometers. Recent observations show highly
variable features, which may change rapidly with waves of ionization, excitation,
condensation, etc. This is particularly noticeable when they are cool and low in
density. Others are violent ejections of material at velocities of several hundred
km/s forming very narrow, long jets. Above active regions, coronal enhancements
show that changes are continually taking place especially in the shape of temporary
loops that form and dissipate repeatedly and that may even be the site of violent
Figure 10.7: Diagram showing the main structures in the inner corona and the chromospheric extensions
over various quiet (A and B), or active areas (C).
SOLAR ECLIPSES AND STUDYING THE CORONA 131
We come to the field involving study of the outermost region of the solar
atmosphere; the K or plasma corona, the region that consists of completely ionized
material. It contains mainly protons and electrons trapped by the Sun’s overall
magnetic field lines. This also entails investigation of mass loss from the Sun and
of its counterpart in the interplanetary medium, the solar wind. Fig. 10.8 shows
the principal features involved; very narrow features, like the exceptionally long
coronal steamers, which fade at several solar radii into another, tenuous corona.
Light arises from scattering of solar radiation by free electrons in the ionized gas
(Thomson scattering). This radiation is so polarized and measurements of the
degree of polarization enable the location of the features in space to be determined
and enable them to be seen out to greater distances. The shape of the corona is of
great interest as the shape depends on the distribution of the magnetic field lines.
In fact, this shape is difficult to determine during an eclipse, as all that can be seen
during an eclipse is the projection of the various features around the Sun onto the
plane of the sky. Different features appear superimposed on one another along the
line of sight. The features may, in fact, be intertwined. Most features do have one
predominant direction, and this is generally radial, but they are also nearly always
curved to a greater or lesser extent and many even form arches. These arches are
difficult to detect. They demonstrate the existence of greatly extended magnetic
field lines, the bases of which lie in the deeper layers of the Sun. It is possible to
SOLAR ECLIPSES AND STUDYING THE CORONA 133
Figure 10.8: Schematic representation of the solar corona and its extensions into the interplanetary
medium.
Interesting studies may be carried out with modest equipment, if the film is
calibrated correctly. By making certain assumptions we may deduce the total mass
of the corona at the time of observation. Also, the records may be used to create a
series of isophotes from which the degree of flattening of the plasma corona may
134 ASTRONOMY AND ASTROPHYSICS
Figure 10.9: A schematic diagram for an eclipse experiment designed to photograph the features in the
solar corona with a high spatial resolution.
SOLAR ECLIPSES AND STUDYING THE CORONA 135
in the corona gives rise to a strong radial intensity gradient. To avoid gross over-
exposure of the inner corona, a radial mask or a radial neutral-density filter is
required. The structures within the corona exhibit a very low contrast.
Figure 10.10: An example of the arrangement of the various elements in an eclipse experiment designed to
carry out spectroscopic analysis of the solar corona.
136 ASTRONOMY AND ASTROPHYSICS
EXERCISES
1. Draw a diagram showing the Earth and the Moon in space illuminated by the rays of
the Sun, and then explain the circumstances of eclipses.
2. Show and explain the various phases of a total eclipse (for the measurement of solar
diameter).
3. What are the required conditions for any eclipse to be total?
4. Draw a schematic diagram of an experimental design to measure the second and third
contacts of an eclipse, and then interpret briefly.
5. What cautions should be taken for the observation of corona during the totality of an
eclipse?
6. Explain why the solar corona may be observed during totality.
7. Show the main structures in the inner corona and the chromospheric extensions over
the various quiet or active areas, and discuss.
8. To study the photosphere why is it preferable to work in the red and near infra-red?
9. Discuss briefly the method of observing the lunar corona. Explain how the photographing
of the inner corona is done.
10. Explain with a suitable diagram how the solar corona may be extended to a interplanetary
medium.
11. Describe the method and equipment for photographing the plasma corona, pointing
out the calibration to be made.
12. Why is the photograph of a plasma corona complicated? What type of special
photographic camera is recommended for the purpose?
13. Discuss the arrangement of the various components in an experiment to observe part
of the coronal spectrum.
Chapter 11
MEASUREMENT OF SOLAR
FLUX DENSITY
11.1 INTRODUCTION
W
hen solar radio waves arrive at the Earth, one can measure the following
properties: (i) flux density (i.e., intensity of the solar radiation) at various
single frequencies, (ii) frequency spectrum, (iii) polarization, and (iv) source
brightness distribution at one or more frequencies. Due to the transient nature of
many solar phenomena it is necessary to observe and record these parameters
simultaneously and continuously. The CSIRO observatory at Culgoora has
comprehensive facilities for ground-based meter-wavelength solar observations.
Its main instruments are:
(i) Radiometers for measuring flux density at 43, 80, 160, and 327 MHz.
(ii) Two-antenna short-baseline interferometers for radiometer calibration and for
quiet-Sun flux density measurements.
(iii) Scanning radio-spectrograph, recording from 8 to 8000 MHz the dynamic
spectrum. The recording is on 70-mm film and on a facsimile paper chart.
(iv) Acousto-optical radio-spectrograph, recording the dynamic spectrum of very
weak bursts in selected frequency bands between 30 and 335 MHz.
137
138 ASTRONOMY AND ASTROPHYSICS
In this section first we will describe the basic principles of antennas as applicable
to radio astronomy. Next, we will discuss the antenna requirements from the point
of solar observations.
GW
F = ⋅ (11.1)
4 π r2
The constant of proportionality G is the gain of the antenna. G is a function of
the direction of transmission and has a maximum value G0 along the axis of the
main antenna beam.
In receiving antennas, including radio telescopes, it is more appropriate to
specify its properties in terms of the effective area. If w is the power received by an
antenna of effective area A (in square meters) when exposed to a plane wave with
flux density F, then
w = A F. (11.2)
The effective area is a function of direction and is proportional to the
corresponding gain
4π A
G = ⋅ (11.3)
λ2
Here, λ is the wavelength.
In solar radio astronomy we are concerned with sources of wideband noise-
like radiation. The received power w can be specified by equating it to the thermal
noise power obtained by replacing the antenna with a resistor at the antenna
temperature Ta and is expressed as,
MEASUREMENT OF SOLAR FLUX DENSITY 139
w = K Ta Δf. (11.4)
(K = Boltzmann’s constant and Δf = the receiving bandwidth.)
For a noise source, we may define the flux density per unit bandwidth at the
receiving antenna as S (in W m–2 Hz–1), so that
F = S Δf. (11.5)
It is seen from Eq. (11.2) through (11.5) that, if the receiving aperture area is A
perpendicular to the direction of the source and the source is small compared to
the antenna beam width, then
S A = K T. (11.6)
The received power is obtained by the angular distribution of brightness in the
source. For each element of the source, the brightness temperature TB is defined as
the temperature of a black body for which the brightness of the thermal radiation
would equal the brightness observed actually. TB may change with frequency, but,
like S, it is assumed to be constant over the receiving bandwidth Δf. The power
received from a source element subtending a solid angle dΩ is
KTB A
dw = 2
Δf dΩ. (11.7)
λ
Hence, the total received power is given by
K Δf
w =
λ2 ∫∫ T
4π
B A dΩ. (11.8)
These equations are used to find expressions for the signal received from a
uniformly bright source of angular size considerably less than that of the antenna
beam. The source is assumed to be located on the axis of the beam, where the
effective area has its maximum value A0. Then
K TB Ω
S = (where Ω = angular size) (11.9)
λ2
and
A0TB Ω
Ta = ⋅ (11.10)
λ2
For the special case of a circular source of angular diameter d,
K TB π d 2
S = (11.11)
λ2 4
140 ASTRONOMY AND ASTROPHYSICS
A0TB π d 2
Ta = ⋅ (11.12)
λ2 4
Ta = TB. (11.13)
Hence, we can write,
∫∫ A d Ω
2
=λ. (11.14)
4π
Eq. (11.14) shows a relation between the antenna beam-width to the effective
area (or to the gain). Therefore, for beams of circular cross-section we have,
4 λ2
A0 = k (11.15)
π α2
and
16 k
G0 = 2 ⋅ (11.16)
α
Here, α is the half-power beam-width, k is a constant for a given beam shape,
with 1 > k > 0.
The need for uniform response sets a minimum beam-width and thus limits the
receiving area. The absolute sensitivity of the receiving system is obtained by
the relation between the wanted signal from the Sun and the unwanted statistical
fluctuations owing to noise from the bright sky within the antenna beam and from
the receiver itself. This noise component is expressed as an equivalent antenna
temperature TN. The receiver contribution to TN at meter wavelengths is typically
of the order of 300 K, since the sky brightness temperature TSK is usually uniform
over a solid angle greater than the antenna beam. The effect of the sky radiation
is to increase the antenna temperature by an amount TSK. Sky temperatures are
shown in Fig. 11.1. The figure shows the following: (a) upper limit and (b) lower
limit of sky noise component; (c) typical receiver noise, expressed as an equivalent
antenna temperature; (d) and antenna temperature corresponding to quiet Sun at
sunspot minimum, for an antenna with an effective area of 20 m.
waves. At elevation angles less than about 40° the vertically polarized ground
reflection is from 3 to 20 dB weaker than the horizontally polarized component.
With a suitable choice of dimensions this array operates over a band slightly
less than the range between the resonant frequencies of the longest and shortest
dipoles. Bandwidths of from 3:1 to 10:1 are readily obtained. At any one frequency
the antenna current is concentrated in the dipoles, which are close to resonance.
The LPA exhibits almost constant gain (typically 5 to 10 dB) and beam-width over
2
its operating band; the effective area is directly proportional to λ .
Figure 11.2: A log periodic dipole array. The spacing between adjacent dipoles is proportional to the dipole
length. The feeder (twin-wire transmission line) is transposed as indicated between successive
dipoles to provide the required phase relations. The feeder is terminated by a short-circuit
about λmax/8 beyond the longest dipole, where λmax is the longest operating wavelength.
The conical helix is a wire helix wound to fit a conical surface. It is a broad-
band antenna somewhat analogous to the LPA. The receiver is connected to the
apex of the helix. The response is circularly polarized corresponding to that of
the winding. Antennas of this type have been used as the elements of the Clark
Lake radio-heliograph and the operating wavelength range is 2.5 to 15 m. At
wavelengths up to a few meters a paraboloid reflector with an LPA as its feed
is a good broadband receiving antenna. The effective area is constant but the
beam-width is proportional to wavelength. The position of the phase center of
the LPA feed varies with wavelength causing defocusing. However, the resulting
deterioration in the antenna pattern can, in general, be made very small.
144 ASTRONOMY AND ASTROPHYSICS
11.6 RADIOMETERS
Radiometers are designed for measuring accurately the total solar flux density
at a single frequency. At meter wavelengths this measurement presents problems
not encountered at shorter wavelengths. The flux from the quiet Sun is smaller,
the background sky brightness is greater, and the receiving antennas required are
larger. The high fluxes encountered during solar bursts also introduce problems of
calibration and dynamic range.
The parameters involved are shown in Fig. 11.3. This shows typical flux
densities for quiet Sun and for large solar bursts. Maximum and minimum galactic
background flux densities are also revealed. The variation of the background
radiation around the ecliptic at 80 MHz is shown in Fig. 11.4. In the presence of
the Sun this background is difficult to measure as it varies with a scale size similar
to the beam-widths of typical antennas. Other characteristics of solar bursts
affecting the design of radiometer receivers are the rapid variation of flux with
MEASUREMENT OF SOLAR FLUX DENSITY 145
4
both time and frequency. Rise rate as high as 10 solar flux units per second are
sometimes noticed. Time constants and sampling intervals between 1 second and
1 millisecond are usually chosen. Good rejection of signals is also important as
instantaneous variations in flux up to 30 dB in a few megahertz span may also
occur. Individual bursts have bandwidths in the range 0.1 to 100 MHz or more.
The receiver bandwidth is thus a compromise between sensitivity and the need to
measure the true flux of narrow-band features. Bandwidths from 0.1 to 10 MHz are
used in general.
Figure 11.3: The flux densities of the quiet Sun (at sunspot minimum) and of large solar bursts of
several types. The range of the galactic background flux is also shown (20° beam-width is
assumed).
In order to accommodate a wide range of fluxes encountered between quiet Sun and
intense bursts, radiometer receivers need dynamic ranges of ≈ 40 dB. The recording
and presentation of data require special attention. If the range of fluxes is 40 dB
and the lowest flux is recorded with a precision of 1%, then a total recording range
6
of 10 :1 is required if a proportional response is used. Even with digital recording
146 ASTRONOMY AND ASTROPHYSICS
an elaborate system using 20-bit binary numbers would be needed to give the
precision over the whole range. One method of compressing the output range is
to use a form of automatic gain control. In it the receiver gain is automatically
reduced in one or more steps when the input signal exceeds predetermined levels.
This method is unsuitable if direct analog recording, e.g., chart or film, is used. The
automatic gain changes make such records almost indecipherable when rapidly
varying bursts are recorded. Alternatively, two separate receiving and recording
channels can be used. One of these records the weaker signals accurately while the
latter are recorded on the second channel of much lower gain.
Figure 11.4: Galactic background flux and brightness temperature measured around the ecliptic at 80
MHz. Both measurements are averaged over a 20° beam-width. The position is given as right
ascension on the ecliptic; the upper horizontal scale indicates the right ascension of the Sun
month by month.
Very high precision is not required for the measurement of bursts. For these
the recording scale can be greatly compressed by using a receiver with logarithmic
response (i.e., one with the output voltage proportional to the logarithm of the
input power). Since high precision is required for quiet-Sun signal levels, a suitable
compromise is to use a receiver response that follows a steep logarithmic law at
low signal levels and changes progressively to a flatter one as the signal increases.
Such a response is ideal for display purposes. It has a minor disadvantage of
requiring a complicated step-by-step calibration procedure, as shown in Fig. 11.5.
MEASUREMENT OF SOLAR FLUX DENSITY 147
Figure 11.5: A typical calibration for the Culgoora radiometers, which employ dual-slope logarithmic
responses. The output of the 6 × 106 K Zener diode noise generator is attenuated in steps of
1.5 dB.
11.7.1 Calibration
To convert radiometer output signals into absolute values of flux density a
calibration of the instrument is required. This can be done by comparing the output
signal produced by the Sun with that obtained when the antenna input is replaced
by the signal from an absolutely calibrated noise generator. This method requires
accurate calibration of the noise source and also the measurement or calculation
of the effective area of the antenna. Both these intermediate steps must be
repeated periodically to guard against errors due to any changes in the equipment
parameters so that this absolute calibration method is laborious and subject to error.
It is more reliable to use as a calibrator one of the strong cosmic sources of small
angular size for which the flux density has been accurately determined. Usually it
is not possible to observe the cosmic sources during the solar observing period; the
Sun is much stronger than the calibration sources, so that solar radiation picked
up in the antenna side-lobes prevents an accurate daytime measurement of the
calibrator. It is therefore better to use an artificial noise generator as a secondary
standard, comparing its output first with the calibrator and then with the Sun.
When this is done, it is necessary for the noise generator to be stable in output over
periods of hours or days; it need not be calibrated absolutely.
If an artificial noise generator is used as the primary calibration standard for
a radiometer with a linearly polarized receiving antenna, the calculated solar flux
148 ASTRONOMY AND ASTROPHYSICS
must be doubled to allow for the rejected polarization component. This correction
is, of course, not to be applied when the primary standard is a cosmic radio source.
The widely used and by far the most intense calibration source is Cassiopeia-A.
This is a recent supernova remnant. Its intensity is decreasing at about 2% per
year. The flux density and its secular variation have been accurately determined.
Accepted values are:
The spectra of both Hercules-A and Virgo-A maintain fairly constant slope.
The flux density at 178 MHz has been used, together with the slope between 38
and 750 MHz:
The problems of accurately correcting for sky noise are more serious with weak
cosmic sources than with the Sun. This means that total power observations of
cosmic sources are simply not possible with solar radiometers. An alternative is to
determine the effective area of the antenna by observing cosmic sources using two
identical antennas as an interferometer. This can be done even when the source flux
densities are very much less than the background sky flux. Provided the antenna
system is stable the measured effective area can then be used to determine solar
flux densities when one of the antennas is used alone as a total power radiometer.
MEASUREMENT OF SOLAR FLUX DENSITY 149
Two types of artificial noise source are frequently used at meter wavelengths:
(i) Thermionic diode and (ii) Zener diode noise sources.
11.8 SPECTROMETRY
(24 to 666 MHz) the dynamic range required is very large. The spectrum in these
ranges is recorded at two sensitivity levels differing by 10 dB.
12
In each frequency span, the spectrum is sampled at 200 points.
MEASUREMENT OF SOLAR FLUX DENSITY 151
Figure 11.6: Block diagram of a scanning spectrograph showing components of a single range unit, A-L (in
the dashed box) and the display and recording system, M-R used at Culgoora.
Figure 11.7: The Culgoora spectrograph showing in schematic form, which shows the various control
waveforms applied to a single range unit (inside the dashed box) and the x and y waveforms
applied to the CRO display. Other range units are multiplexed to the CRO display via the
integrator and channel samplers shown. The separate spectra can be displayed on the screen
in parallel; by use of a system of mirrors these spectra are arranged in a single line and in the
correct frequency sequence to be photographed.
For all practical purposes the available recording media for modern solar
radio-spectrographs, with their high data rates, are photographic film, facsimile
records, and magnetic tape.
A suitable antenna for solar radio telescopes should have sufficient directivity to
discriminate against non-solar radiation, while at the same time its beam must be
wide enough to give uniform coverage to all parts of the radio Sun at the operating
wavelength. In conjunction with the receiver to be used the antenna should have
sufficient receiving area to detect the quiet Sun.
The threshold antenna gain Gt required to detect the quiet Sun in the presence
of other noise sources is
4π(Tg + Tr )
Gt = ⋅ (11.22)
ΩTd τ Δf
Here, Tg is the sky noise temperature, Tr receiver noise temperature, Td
apparent disk temperature of the quiet Sun, Ω solid angle subtended by the Sun, τ
integration time, and Δf receiver bandwidth.
Fig. 11.10 (a) shows the characteristics of reflectors of various diameters, and
of several arrays of log-periodic antennas. Fig. 11.10 (b) is a performance diagram
for the antennas actually used in the Culgoora spectrograph. Other antenna
configurations can also be used. For example, the Harvard spectrograph used a
MEASUREMENT OF SOLAR FLUX DENSITY 155
single 8.4-m paraboloid for the frequency ranges 100 to 580 MHz and 2000 to 4000
MHz, and a simple array for 25 to 100 MHz.
Figure 11.10: (a) Gain-versus-frequency characteristic for paraboloid dishes and various combinations
of log-periodic antennas. (b) The antennas used for the 8 to 8000 MHz Culgoora radio-
spectrograph compared with the threshold level of the quiet Sun.
EXERCISES
1. What are the properties that can be measured from solar radio waves arriving at the
earth?
2. Mention the comprehensive facilities of the CSIRO observatory at Culgoora for ground-
based meter wavelength solar observation.
3. Describe the basic principles of antennas as applicable to the radio astronomy.
4. What is ‘noise-like’ radiation?
5. Establish a relationship between the antenna beam-width to the effective area of the
antenna.
6. Summarize the antenna requirements for solar observations.
7. The ground-reflection problem is common to all mechanically steerable radio telescopes.
Explain how this can be minimized through a careful design.
8. Discuss why a steerable paraboloidal reflector is a convenient form for a solar radio
telescope antenna at wavelengths from 1 to 5 meters.
9. What type of antennas is preferable for spectrographs and multi-frequency
radiometers?
156 ASTRONOMY AND ASTROPHYSICS
10. Give the principle of broadband antenna, clearly explaining the geometry. What is the
utility of dipole arrays?
11. Discuss how radiometers are suitably used for measuring accurately the total solar flux
density at a single frequency.
12. In order to accommodate a wide range of fluxes, what are the necessary modifications
to be made in radiometer receivers?
13. To convert radiometer output signals into absolute values of flux density, a calibration
of the instrument is required. Discuss how this can be done.
14. Describe a radio-spectrograph. Classify the three different types of radio-
spectrographs.
15. Which was the first solar radio-spectrograph?
16. Mention the characteristics of antenna at different frequency bands of a radio-
spectrograph.
17. Give the block diagrams of a scanning spectrograph showing the different
components.
18. Explain how multi-channel spectrographs can be used to observe weak solar bursts.
19. What is an acousto-optical spectrograph? Draw a diagram to explain the principle.
20. What are the qualities to be maintained in an antenna for a radio telescope?
21. Show the gain vs. frequency characteristic for paraboloidal dishes and various
combinations of log periodic antennas, and explain.
Chapter 12
MEASUREMENT OF SOLAR
BRIGHTNESS DISTRIBUTION
12.1 INTRODUCTION
M
ost solar features, particularly at meter wavelengths, develop and
change rapidly so that high time resolution is necessary. The elegant
but slow methods of Earth-rotation aperture synthesis are usually not
applicable. In this chapter, we discuss radio-heliographs (i.e., radio telescopes that
produce images of the Sun in real time and in rapid succession). To get useful
information on the position and structure of localized radio sources on the Sun
it is required to achieve an antenna beam-width no greater than a few minutes of
arc. We need apertures of at least 1000 λ. At radio wavelengths longer than about
10 cm it is impracticable to build single filled-aperture antennas of sufficient size.
At meter wavelengths, spaced arrays of many small antennas have been used to
give the requisite overall dimensions. The elements of such arrays are directional,
and for hour-angle coverage they must be steered during observations so that
their beams are always pointed at the Sun. In many cases facilities are included for
measuring the polarization of the radio sources by providing the antennas with
crossed-dipole feeds and quadrature phase switches.
157
158 ASTRONOMY AND ASTROPHYSICS
Figure 12.1: Coordinate system for the array geometry. The array is in the x-y plane. PO is the direction of
the incident ray from the source; l, m are direction cosines of PO.
MEASUREMENT OF SOLAR BRIGHTNESS DISTRIBUTION 159
For a horizontal array, θ is the zenith distance and φ is the azimuth. It follows
from the reciprocity theorem that the response pattern in reception is of the same
form as the radiated field pattern in transmission. An example of a linear co-phased
array of small antennas is shown in Fig. 12.2. The field pattern consists of multiple
fan-shaped lobes. With a uniform grading the pattern exhibits large side-lobes,
which can be reduced by altering the grading (i.e., changing the weighting of the
fourier components). In this example the required reweighting can be produced by
tapering the current distribution from the center to the ends of the array (Fig. 12.2
(d) and (e)).
Figure 12.2: Grading and field response patterns for an array of eight small elements uniformly spaced
along the x-axis: (a) uniform grading; (b) gives the pattern; (c) a suitably tapered grading;
(d) gives a response; and (e) with low side-lobes. For nearly normal incidence the direction
cosine z angle of incidence.
where c(u, v) is the transfer function of the antenna and is equal to the autocorrelation
function of the grading and g* is the complex conjugate of g. For an array of small
spaced elements the grading g is zero everywhere in the x-y plane except at the
location of each element. It is seen from Eq. (12.3) that the transfer function c(u, v)
is non-zero only at a set of discrete points in the u-v plane, each at a position where
the u and v coordinates are respectively equal to the x- and y-components of the
line joining a pair of the array elements. Fig. 12.3 shows the transfer function of an
annular array with uniform grading.
this spectrum and can be computed by one of the standard methods of numerical
Fourier transformation.
The above method of image formation has some attractive features. As the
Fourier components are separately accessible for processing, various operations
such as the reweighting of components and the application of phase and amplitude
correction can be carried out directly before performing the Fourier transform.
All the correlations are simultaneously measured and thus no loss of sensitivity
due to time-sharing occurs. In a typical receiving system each correlator has a
circuit for forming the instantaneous product of the two signals at a sampling
rate of the order of twice the radio frequency bandwidth. For solar observations
the value is usually a few megahertz. For practical purposes two correlators are
used for each pair of channels with a 90° difference in the relative phases of the
input signals. From the two outputs the phase and magnitude of the complex
correlation product can be achieved. Successive values of each correlation product
are summed to improve the signal-to-noise value. At intervals equal to the time
between successive images the integrated values are transferred to the device to
perform the Fourier transform. Usually it is a general-purpose digital computer.
The correlator input signals are quantized by representing any positive value
by 1 and any negative value by 0. This one-bit correlation procedure is valid for
noise-like signals. After correction the final time-averaged result is the same as if
exact values were used for the correlator inputs.
Figure 12.4: Power polar diagram of a uniformly fed grating array, with N = 16, d = 57.2 λ. The broken line
is the polar diagram of the array elements, which are directed toward the zenith.
For solar applications the beam spacing should be greater than the diameter of
the radio Sun. Generally, beam spacing of the order of 1° has generally been used
at wavelengths up to 1 m, so that the element spacing d is ≈ 60 λ. At longer meter
wavelengths, where the diameter of the radio Sun is over 1°, spacings of 20 to 30 λ
are more appropriate. Since angular resolution of a few minutes of arc is required
for solar studies, the aperture D requires the order of 1,000 λ. These considerations
determine the required number of elements which has been from 8 to 32 in grating
array radio-heliographs.
An early meter-wavelength grating radio-heliograph was developed at Nancay,
France. This East/West array was operated at 169 MHz and used 32 paraboloid
antennas on a 1550-m (870 λ) baseline, giving 3′.8 fan beams, 2°.01 apart.
2πa⎞
A(θ) = J 02 ⎛⎜ ⎟. (12.5)
⎝ λ ⎠
Here, a is the array radius, θ is the angular distance from zenith, and J0(x) is
the zero-order Bessel function of the first kind. As the array uses spaced elements,
there are further grating responses. These are complicated ring patterns of radius
θ1, 2θ1, 3θ1, …, centered on the main pencil beam:
164 ASTRONOMY AND ASTROPHYSICS
θ1 = λ . (12.6)
d
Here, d represents the spacing between adjacent array elements.
The dimensions of the array were chosen to give a field of view of approximately
2 clear of grating lobes, at the original operating frequency of 80 MHz. The fields
and beam-widths now obtained at the four observing frequencies are shown in
Fig. 12.7. At any one of these frequencies two images per second are built up,
each with 60 × 48 picture points. A combination of multiple beams and phase
scanning is used. A branching network forms a North/South row of 48 pencil-
beam responses and these are shifted from East to West in 60 steps to form the
picture. The two pictures revealing LH and RH circular polarization are produced
quasi-simultaneously. The signals for each column of picture points are integrated
for about 8 ms, then polarization switches and the antennas are reversed before
the next 8-ms integrations. Thus, columns of points are recorded for each picture
in turn. The operating frequency is changed from second to second in a sequence
selected by the operator.
It is clear from Fig. 12.6 that the J02 pencil beam has unacceptably high side-
lobes. Therefore, some methods of correction equivalent to reweighting the spatial
MEASUREMENT OF SOLAR BRIGHTNESS DISTRIBUTION 165
Figure 12.6: The power polar diagram of the Culgoora 96-element circular array with co-phased elements:
(a) shows the radial cross-section of the pattern, with the central J02-beam and first-order
grating response. The pattern in two dimensions is produced in (b), which is a photograph of
the diffraction pattern of an optical model of the array.
and determines which picture points are brighter than a certain threshold. The
intensities and positions of bright points are recorded on the magnetic tape. When
the recorded picture is reconstituted from the record the intensity of each of these
discarded points is set at zero. The computer finds the threshold level for each
picture as a prescribed fraction of the intensity of the brightest point in the picture.
Fainter points are below the level of side-lobes from the strongest source and so
unreliable.
Figure 12.7: The field of view and picture format for the Culgoora radio-heliograph. For each frequency the
recorded field of view is shown with the outline of the optical Sun. The black filled circle shows
the half-power contour of the array beam on the same scale.
outputs of the two arrays are connected to a common receiver via a switch, which
reverses their relative phase at some convenient audio frequency. The switching
frequency component of the output is isolated by a synchronous demodulator
circuit. If V1 and V2 are the instantaneous signal voltages from the two arrays it is
easy to show that the amplitude of this modulated component is proportional to
the time-average of V1V2. Thus, the modulated receiver output corresponds to a
power polar diagram, which is the product of the voltage polar diagrams of the
two arrays. An example of a cross-correlation system is the Mills Cross, which
uses two orthogonal linear arrays with a common phase center. This is shown in
Fig. 12.8 (a) while typical voltage polar diagram is shown in Fig. 12.8 (b). The
product of the two fan-beam patterns is essentially zero except where they intersect
in the central region, so that the cross-correlation power polar diagram is a pencil
beam corresponding to the hatched area of Fig. 12.8 (b).
Figure 12.8: (a) A pair of linear antennas arranged as a Mills cross. (b) A cross-section of the fan-beam
field response patterns of the two arrays. The pencil-beam power polar diagram of the cross-
correlation response is revealed by the hatched area.
Figure 12.9: The crossed grating interferometer. (a) Geometry of a crass array. Each baseline occurs twice,
e.g., A1B1 and A2B2. (b) The equivalent T array. (c) Central part of power polar diagram.
Circles show positions and relative responses of pencil beams. Alternate pencil beams give
positive and negative responses—filled circles represent negative response.
required separation between adjacent fan beams. The second array CD consists
of two widely spaced elements. The cross-correlations between AB and CD give
Fourier components corresponding to all the multiples of the unit spacing b. The
length of the array is thus increased effectively from (n–1)a to 2na, and the angular
resolution is approximately doubled by the addition of only two extra elements.
Several compound arrays have been used as radio-heliographs.
EXERCISES
13.1 INTRODUCTION
T
o determine the polarization state of a radio signal, it is required to measure
four independent quantities. In the case of solar radio emission we know that
with the receiver bandwidths that are normally used there is no possibility
of detecting any linearly polarized component and we are concerned solely with
the analysis of circular polarization.
13.2 POLARIMETERS
The quantities measured are the received signal intensities IR and IL using right-
hand (RH or R) and left-hand (LH or L) circularly polarized antennas, respectively.
In radio astronomy, circular polarization of the incident wave is said to be RH if
the field vector rotates clockwise when viewed in the direction of propagation.
The degree of circular polarization P is defined as
(I R − I L )
P = . (13.1)
(I R + I L )
173
174 ASTRONOMY AND ASTROPHYSICS
the color film of the spectro-polarimeter because color is used to represent the
polarization. However, the latitude of color film is very narrow, and changes in
density tend to cause color variations also. Thus, it is wise to record polarization
on the color film and to retain the separate black-and-white intensity record.
EXERCISES
1. What are the four independent quantities required to determine the polarization state
of radio signals.
2. Define mathematically the degree of circular polarization. What is percentage
polarization?
3. Draw a nomogram for a quick calculation of the degree of polarization, and discuss.
4. Discuss single-frequency polarimeters.
5. “The greatness problem associated with the observation in solar emission is the
presence of background components.” Explain how the problems can be minimized in
a radio-heliograph study.
6. What is a swept-frequency polarimeter?
7. Give a simplified schematic diagram of a spectro-polarimeter, and discuss its
operation.
180 ASTRONOMY AND ASTROPHYSICS
Chapter 14
RECEIVERS AND ANTENNAS FOR
R ADIO ASTRONOMY
14.1 INTRODUCTION
R
adio astronomy is the study of the radio waves naturally emitted by
astronomical objects like the Sun, Moon, planets, our own galaxy and other
galaxies, and clusters of galaxies far distant from our own. The amounts
of radio energy received are very small. The science has grown mainly because
experimenters have built antennas of very large collecting area and have observed
with receivers of very high sensitivity and stability. This chapter describes the
principle behind the design of such antennas and receivers and shows how such
instruments have been built and used in recent years.
A radio astronomy receiver (or radiometer) is the device for measuring the radio
frequency power available from the antenna system. The receiver amplifies the
power over a range of frequencies centered on the frequency being used. It thereby
181
182 ASTRONOMY AND ASTROPHYSICS
detects the power, and investigates or filters the result before recording it. Fig. 14.1
shows the general elements of such a receiver. We define some of the important
parameters of the receiver below:
(i) Central Frequency (f MHz). The output of the antenna (Fig. 14.1) is fed to the
preamplifier to provide gain over a range of radio frequencies. The central
frequency is the frequency at which observations are made and is usually at
the center of the preamplifier frequency response.
(ii) Gain (G). The receiver is intended to convert the very small change in power
of about 10–17 W due to a radio astronomy source in the antenna beam into a
power capable of being recorded. The gain may be defined as,
(iv) Time Constant (τ sec). The signal leaving the detector is still fluctuating with
time. It is required to smooth out the fluctuations before recording the signal.
In its simplest form, this smoothing is done by a resistor and capacitor network
forming a low-pass filter or by making a time integral of the detector output
over a fixed integration time. Whatever means of smoothing is adopted, the
result is to provide a number of statistically independent measures of the
detector output per second of time. The reciprocal of this number is referred to
as the time constant (τ) of the integrator.
(v) Effective Input Noise Temperature (Te° K). Practical receivers generate noise
within themselves and also the antenna may generate or receive unwanted
noise. A good radiometer is one which allows radio astronomical signals
detected most easily and measured most accurately in the presence of such
noise. To describe this noise we consider the behavior of a simple receiver
(Fig. 14.2 (a)). Assuming that the receiver generates no noise, and that it is
connected not to an antenna but to a matched resistor in an enclosure at
T° K, the available noise power, which such a resistor can deliver to the receiver
within its bandwidth B is
P = kT B watts. (14.2)
In the above equation, k is Boltzmann’s constant and the output meter will
show a deflection due to this input power. Eq. (14.2) is derived from Nyquist’s
formula for the thermal noise power generated in a resistor. If the temperature of
the resistor were reduced to 0° K no noise power would be delivered to the noiseless
receiver (Fig. 14.2 (b)) and the output meter deflection would be zero. Under these
conditions, assume the receiver will become noisy (Fig. 14. 2 (c)), generating its own
noise power and deflecting the recorder by δ. Then we assume that the receiver
returned to its noise-free state and the temperature of the resistor increased to
a value Te so that the same deflection δ occurs (Fig. 14.2 (d)) as with the noisy
receiver. The effective input noise temperature of the noisy receiver is defined as Te.
Te = (F – 1) × 290° K (14.3)
It can be said that all radio astronomy signals are white noise, which has the
following important characteristics:
(a) Its power spectral density function is constant over all frequencies.
(b) Its autocorrelation function is a constant times the unit impulse function.
If the radio astronomy signal were truly white noise as it impinges on the
antenna, its characteristics would be modified by the antenna and receiving
system. Before the noise signal goes the detector, it is passed through various parts
of the receiver. It is now nearly true to describe the noise power as bandwidth-
emitted white noise. The degree to which this approximation approaches the truth
depends on:
(a) The power spectral density function of the noise before it reaches the
antenna, and
(b) The width of the pre-detector bandwidth.
Evidence for (a) is derived by measuring the power received from a radio
source within a finite bandwidth whose center is placed at various frequencies
throughout the usable radio spectrum. The type of spectra measured is shown in
Fig. 14.3, which illustrates the fact that though no radio source is a white-noise
emitter, the noise power received from such sources can be described as bandwidth-
limited white noise so long as the pre-detector bandwidth is chosen to be so small
that the spectral density function of the source is nearly constant over the chosen
bandwidth. So for various measurements, different bandwidths are chosen. The
condition is met in cases (a) to (c) of Fig. 14.3 if bandwidths of a few MHz are
used. For the study of neutral hydrogen (Fig. 14.3 (d)) bandwidths of 10 to 20 kHz
are usually used. It follows that the characteristics of the radio astronomy noise
powers to be measured are those of bandwidth-limited white noise. Such signals
are measured in the presence of noise originating in the atmosphere, the antenna,
and the whole receiver system. It is adequate to describe all these unwanted noise
powers also as bandwidth-limited white noise.
When a radio astronomy signal is received, a method of measuring the
available power from the antenna terminals is that of antenna temperature. The
power available from the antenna is related to an antenna temperature TA as
RECEIVERS AND ANTENNAS FOR RADIO ASTRONOMY 185
To detect the smallest signals the radiometer must be sensitive and to measure
these signals accurately it must be stable. In fact, some stability is needed even to
detect a source, but it of course depends on the techniques that must be used.
(i) Sensitivity. The sensitivity of a radiometer is measured by the smallest change
in TA, which an instrument can detect. Even if the stability were perfect, the
total noise (i.e., noise from the radiometer Te together with antenna noise and
atmospheric noise) limits the sensitivity of the instrument. If the total noise in
the radiometer system be represented by Top ° K, where Top is the operating noise
temperature, then the root mean square (rms) fluctuations of the radiometer
output signal (ΔTA), measured in terms of changes of antenna temperature is
expressed as,
186 ASTRONOMY AND ASTROPHYSICS
Top
TA = c . (14.5)
τB
In Eq. (14.5), c is a constant of the order of unity and τ is the time constant
or integration time used. For high sensitivity T should be reduced, while B
and τ should be increased. The reduction of Top by the development of low-
noise receiving techniques is one of the main advances in radio astronomical
techniques.
(ii) Stability. Changes of the gain, the bandwidth, or the noise temperature of
a radiometer are examples of possible unstable behavior. The techniques of
improving stability by good design and construction of circuits, the use of
switched radiometer systems, and other devices forms the second area of
technical advances.
The problem of stabilizing radiometers can be discussed in the light of
Eq. (14.5), which shows that ΔTA, the rms noise fluctuations should have, for a
given Top, bandwidth and time constant. It is desirable that observations should
be limited only by ΔTA, but in practice the changes in gain, noise temperature,
or bandwidth of the radiometer may be the limiting factor in radio astronomy
observations. There are many ways in which the effect of these changes has been
minimized, but it should still be recorded that in practice very few radiometers
operated on radio telescopes give fluctuations in their outputs as low as those
indicated by Eq. (14.5). It is evident that the overall gain of the system must be
stable to a high degree. If no means were adopted to remove the effects of the
operating noise temperature Top, any gain change ΔG would cause a change at
the output corresponding to a ΔT at the input of Top ΔG/G. With earlier receivers,
where Top was about 1000° K, this required gain stability of 0.1% to keep ΔT down
to 0.1° K, while even with these receivers the rms noise fluctuations (B = 10 MHz,
T = 10 sec) could be only 0.1° K. When Top is reduced with the newest receivers
there is low-noise fluctuations, and to take advantage of the reduced fluctuation
noise still requires an equivalent fractional gain stability. The means by which this
stability is achieved are:
(a) Total power radiometers.
(b) Dicke radiometers.
(c) D-C comparison radiometers.
(d) Frequency-switched radiometers.
(e) Phase-switched radiometers.
(f ) Correlation radiometers.
RECEIVERS AND ANTENNAS FOR RADIO ASTRONOMY 187
Te 2 T
Te = Te1 e 3 ... . (14.6)
G1 G1G2
Here, Te is the overall noise temperature due to amplifier stages with noise
temperatures Te1, Te2, etc., and gains G1, G2, etc. Eq. (14.6) shows that the noise in
the first amplifier appears directly with the signal, but noise in later stages is of less
importance since the signal and first-stage noise have been amplified by the earlier
stages. Since power gains of about 100 (20 dB) are usually achieved in the first
stage of a receiver, the contributions to Te from later stages become small and can
be neglected. Average values for the contribution to the temperature of an antenna
due to radiation from the galactic background and from the atmosphere are shown
in Fig. 14.4.
Figure 14.4: Average value for the contribution to the temperature of an antenna due to radiation from the
galactic background and from the atmosphere is shown.
188 ASTRONOMY AND ASTROPHYSICS
In this radiometer, as shown in Fig. 14.5, we must be careful to maintain the gain,
Te, and bandwidth constant by suitable circuit design. The rectified signal at the
detector output due to the noise of the receiving system is balanced against a
constant voltage and the difference signal amplified at DC before recording. Stable
supplies to both heaters and anodes of vacuum tubes, the inclusion of large cathode
resistors to a negative supply, temperature control, and the use of temperature-
insensitive tuned circuits have all been employed. The results are satisfactory, with
gain stabilities of 0.1% maintained over several hours, and reasonable periods
existing when the stability is perhaps even five times better. The total-power
radiometer observes the source continuously, which is not true for most switched
radiometers, so that for a given duration of observing time the noise fluctuation is
lower than for the switched system.
This type of radiometer is shown in Fig. 14.6. It is the most widely used radiometer.
The input of the receiver is switched between the antenna and a constant comparison
source of noise power, usually a resistance at a constant temperature. To a first
approximation the receiver noise is the same in both switched positions, so that
the output of the detector is a signal corresponding to the constant system noise
carrying on it a square wave modulation of peak-to-peak amplitude corresponding
to the difference (TA – Tc) between the antenna temperature and the comparison-
noise-source temperature. This signal is applied to the phase-sensitive detector,
fed also with a square wave at the switch frequency, and the output of this detector
smoothed, amplified, and recorded, is a measure of the difference between the
antenna and noise-source temperatures.
RECEIVERS AND ANTENNAS FOR RADIO ASTRONOMY 189
In Eq. (14.7) M is the gain modulator ratio and Te the noise temperature of the
receiver after the switch. Thus, without a gain modulator TA and Tc should
therefore be the same but with a gain modulator they can still be different and
yet the radiometer can have good stability against small gain changes.
(b) Changes in Te. A change in receiving noise temperature will show a change
in receiver output in all practical conditions unless TA = Tc and thus M = 1.
However, this condition cannot be achieved easily in a normal Dicke receiver.
The need for correct measurements of the 1420 MHz radiation from neutral
hydrogen in our own galaxy and in external galaxies has led to the development of
special radiometers for the purpose. One of the best may be taken as the technique
of frequency switching. Fig. 14.7 shows a block diagram of the radiometer in its
simplest form. Later developments have added more channels and reduced the
bandwidth of the channels.
The first stage of the receiver is a crystal mixer. The frequency switching is
accomplished by supplying a local oscillator signal, alternately, 400 times per
second through an electronic switch (ES) from the two crystal-controlled sources
O1 and O2. The frequencies of O1 and O2 supplied to the mixer differ by 1.08 MHz.
The output of the mixer, after IF amplification at 28 to 33 MHz passes to the second
mixer, the local oscillator of which O3 can be swept in frequency from 34 to 39 MHz.
The output after amplification over 4 to 9 MHz is fed to two separate channels;
one centered on 4.96 MHz and the other on 6.04 MHz. The difference frequency
between these is chosen to be 1.08 MHz—the same as that between O1 and O2.
These two channels give outputs through a third mixer and 540 kHz amplifiers
over a bandwidth of 35 kHz to the detectors and synchronous detector (SD), which
RECEIVERS AND ANTENNAS FOR RADIO ASTRONOMY 191
is also supplied with the 400 Hz switching wave. If a 1420 MHz signal is traced
through the receiver it will be seen to emerge through the 4.96 MHz channel when
O1 is the LO, and through the 6.04 MHz channel when O2 is the LO; O3 is assumed
to be set at 35.764 MHz for this to occur.
It is seen that if 1420 MHz is the required frequency it is always available at
one of the two channels while the other channel gives a 35 kHz wide band 1.08
MHz above or below 1420 MHz. In this way the synchronous detector output is
a measure of the difference between the noise power at 1420 MHz and the mean
of the noise power at frequencies equally spaced above and below 1420 MHz.
Fig. 14.8 shows how the bands lie with respect to a hydrogen profile to be observed.
It is clear that by tuning O3 across the band, the sensitive band of the receiver may
be swept across the hydrogen-line profile.
channels each of 95 kHz wide are used. For galactic hydrogen more channels of
narrower bandwidth are necessary.
The use of phase switching is well known and is illustrated in Fig. 14.10 for a
two-element interferometer. The phase switch in one of its two switch positions
introduces an extra half-wavelength of cable into the line from one antenna to
the receiver and thus changing the phase of the signal from that antenna by 180°.
The phase-sensitive detector is driven in synchronism with the phase switch. The
phase-switching results in a movement on the sky at the switching-frequency of
the lobes of the interference pattern as shown in Fig. 14.11. Thus, if a point source
is being observed against a uniform sky background only the power from the
point source contributes to the recorded fringe pattern. It is also known that for an
extended source, the amplitude of the recorded trace is proportional to one term of
the Fourier transform of the angular distribution of the intensity across the source.
For these two reasons the practice of phase switching has been widely used in
RECEIVERS AND ANTENNAS FOR RADIO ASTRONOMY 193
radio astronomy. The operation of phase switching removes from the output any
contributions from the receiver noise as well as from a uniform sky background,
so that the system has the additional stability associated with this characteristic.
V1 (t) V2 (t)
ρ = 2
. (14.8)
V
2
Here, V is the mean value of [V1(t)]2 and [V2(t)]2. When voltages V1(t) and
V2(t) are completely correlated, ρ = 1. When they are completely uncorrelated,
ρ = 0. A multiplier as shown in Fig. 14.12 with appropriate smoothing will give the
average value and normalizing will provide a measure of ρ.
The most valuable and versatile antenna in radio astronomy is the parabolic
reflector. Such instruments of diameters up to 300 feet have been built on a variety
of mountings. In this section we describe these antennas by outlining the choices
that lead to the design of such an instrument. Also, we describe some typical
antennas, which have already been built and tested.
The polar-mounted 85-ft dish with its simple drive and control system is
an excellent telescope for radio astronomy. The larger dishes are cheaper if
they are on Az/El mounts and cheapest on a simple meridian transit mount.
Coordinate converters, which may be mechanical or optical or which may use
either analogue or digital computer techniques, have all been used in existing
radio telescopes to convert from Az/El to polar coordinates and to permit
the telescopes to be moved and pointed in the celestial coordinate system.
The drives for such Az/El telescopes must be servo mechanisms of high
performance.
5. The Drive and Position-indicator System. The most important characteristics
required of the drive and indicator system are that it can direct the antenna
beam to the chosen direction in space and move the beam either to track
the apparent motion of a radio source or to scan a selected area of the sky.
In radio astronomy the accuracy with which the telescope can be moved
has been defined in terms of the telescope and pointing accuracies of one-
twentieth of the half-power beam-width are considered to be satisfactory. For
an instrument that will be used over a range of frequencies, the smallest beam-
width associated with the highest usable frequency sets this pointing and
guiding accuracy. Radio astronomers have accepted the fact that the telescopes
will deflect under their structural dead load, and so large telescopes have to
be calibrated to find the effect of these deflections on the drive and positional
accuracy.
6. Environment. All parabolic dishes used for radio astronomy have been built
to withstand the local weather environment. Haystack was the first operated
within a radome, although many antennas used for radar and communications
are so housed. Some savings of cost have been made by requiring precise
operation of the radio telescope in winds only up to about 25 mph and thus
accepting the fact that observations may have to be suspended in high winds.
Survival conditions are chosen to suit the particular telescope location.
The greatest gain for a perfect parabolic dish would be achieved with a uniform
illumination by the feed horn across the dish aperture and no illumination outside
the aperture. Such a feed design cannot be achieved, and most paraboloids have
used feed-horn patterns, which diverge from this maximum-gain criterion.
Compromises have to be made between the need for gain and the desire to reduce
the side-lobes and spill-over in the antenna. The main departure from perfection
in parabolic antennas is due to the irregularities of the reflector surface, partial
transmission through the surface, and blocking of the total antenna aperture by
the feed support legs and the equipment near the focal point of the telescope.
All these imperfections lead to a loss of gain of the instrument. The effect due to
transmission through the reflector is kept small by making use of closely spaced
mesh or solid-surface reflectors. The criterion accepted for transmission is that
less than 1% of the energy incident on the reflector surface should be transmitted.
This, in turn, has led to choices of solid-surface reflector materials for antenna
reflectors designed to work at 10 GHz and of expanded aluminum material for
lower frequencies.
The effects of aperture blocking have not been studied, in general, for parabolic
radio telescopes. Diffraction theory can be applied to this problem if the obstacles
to the wave front are large compared with the wavelength used and if the blocking
is assumed to be a geometrical obstruction both of the plane wavefront entering the
telescope aperture and of the ray paths from the surface of the dish to the feed. The
aperture blocking may then be represented as a reduction to zero of the field over
the obstructed part of the aperture and the gain and pattern of the antenna may
be calculated for this aperture illumination. A calculation of this type made for the
case of an aperture blocked by a tetrapod feed support and a circular electronics
canister at the focus showed a loss of gain of 0.5 dB. The antenna pattern was also
modified by the blocking in that the first side-lobes were slightly increased in size
and the next slightly decreased. It is important to know the effect on the gain of
a parabolic antenna of irregularities in the reflector surface. A theoretical analysis
of the problem can be made in the case where only small random phase errors are
produced in the wave front due to irregularities of the reflecting surface. These
irregularities are correlated over finite distances across the surface of the reflector.
The theory gives simple expressions for the cases where the correlation interval
(c), the distance at which on the average the errors become independent of one
another, is much greater or much less than the wavelength used (λ). In large radio
astronomy dishes c is always greater than λ and the gain (G) compared to that for
a perfect surface (G0) is for small phase errors:
G/Go = 1 δ 2 . (14.9)
2
R ≥ 2D (14.10)
λ
for this criterion to be satisfied. Phase path errors larger than λ/16 can be tolerated,
but deformed radiation patterns result. R is impractically large for man-made
sources to be used to measure patterns of large radio telescopes. The Australian
210 ft and the NRAO 300 ft only move from the zenith to an elevation angle of 30°.
The values of R for these instruments near their shortwave limits are both about
50 miles, and at 30° elevation an artificial source would have to be at a height of
25 miles. The use of the stronger radio astronomical sources is possible, provided
those of small enough angular diameter are chosen. The Sun is both too large and
RECEIVERS AND ANTENNAS FOR RADIO ASTRONOMY 199
too variable in intensity for reliable pattern measurements of large dishes. The most
intense of the other radio sources give adequate power for good measurements of
main beam shape.
EXERCISES
15.1 INTRODUCTION
A
t metric and decimetric wavelengths partly a suitable characteristic is the
so-called dynamic spectrum of the radio emission (i.e., the time dependence
of the frequency spectrum). The form of the frequency spectrum and its
behavior in time provide important information on the generation mechanism of
the radio emission.
Wild and McCready (1950) divided the bursts in the metric band into three
types and started the classification of the sporadic radio emission by the nature
of the dynamic spectra. Two further components, Types IV and V radio emission,
have now been added to the Types I, II, and III bursts. Together with the enhanced
radio emission owing to the sunspots these five components are the basic group of
phenomena observed in the metric band. Enhanced radio emission accompanied
by Type I bursts is known as a noise storm while Type II bursts are called outbursts.
The microwave radio emission of the disturbed Sun is not divided into its individual
components in accordance with the nature of their dynamic spectrum largely. This
is due to the lack of spectral observations in this range. On the contrary, it is quite
201
202 ASTRONOMY AND ASTROPHYSICS
possible that at centimetric wavelengths and also in the short-wave part of the
decimetric band the dynamic spectrum as a whole is not a decisive characteristic
for identifying the observed radio emission as it is at metric wavelengths. This is
due to the fact that the dynamic spectrum of the microwave radio emission does
not have the variety of details that distinguish it at the longer wavelengths.
This chapter discusses the features of the sporadic solar radio emission and
the connection between these components and different kinds of geophysical
phenomena. We have thus presented a general picture of the phenomenon of
sporadic radio emission when a center of activity appears on the Sun.
The slowly varying component on a recording of the flux of the Sun’s radio emission
looks like smooth rise in the signal received above the emission level of the quiet
Sun with a characteristic duration of the order of tens of days. This component
is observed mainly in the range of wavelengths from 3 to 50 cm and occupies
a broad spectrum of frequencies of the order of several gigacycles. The relative
contribution from the S-component outside the range λ ≈ 3 to 50 cm is very small.
Observation of the S-component is complicated because of the masking action
of the more intense enhanced radio emission connected with sunspots. The time
dependence of the S-component’s flux reveals a periodicity with a characteristic
time of 27 days, equal to the Sun’s period of rotation. This indicates that the source
of this component is not only the Sun but individual local regions on its disk.
The variations in the enhanced level at centimeter and decimeter wavelengths
correlate well with the relative number of sunspots W known as the Wolf number.
The corresponding correlation coefficient goes up to 0.8. In relative units the radio
emission flux Sw at λ = 10 cm is connected with w by the approximate relation
–3
Sw = 1 + 1.10 W. (15.1)
The intensity of the S-component also correlates with the total area occupied
by spots on the Sun’s disk. In general, the intensity S of the S-component is directly
proportional to the visible area σ of the spots. The coefficient of proportionality
depends on wavelength. In the decimetric band the best index of the radio emission
when compared with the magnitude of σ is the so-called complex area of the spots
σcomp. The correlation between Sw and σcomp implies that radio emission continues
from a center of activity even after a group of spots has decayed.
THE SUN’S SPORADIC RADIO EMISSION 203
Figure 15.1: One-dimensional distribution of radio brightness over the Sun’s disk at λ = 3.2 cm. The figure
also shows the positions of sunspots on the disk during the observations.
Figure 15.2: Total radio emission flux from the bright part of local sources as a function of wavelength:
(a) observations of April 19, 1958; observations of February 15, 1961. (b) The line of dots and
dashes indicates the spectrum of a local source located above a group of flocculi.
THE SUN’S SPORADIC RADIO EMISSION 205
band is the same as the area of the spots and groups of spots and does not change
significantly with wavelengths, then the following conclusions can be drawn. In
the decimetric band Teff for large sources reaches values of about 10° K, and they
are of the same order of magnitude when the wavelength is reduced right down
to 5 to 10 cm. In the centimetric band at wavelengths λ < 5 to 10 cm the effective
2
temperature decreases together with the wavelength more rapidly than λ . This
2
can be confirmed by considering Teff ∝ I λ . For shorter wavelengths the intensity
decreases, and Teff drops more rapidly than λ2. It is possible that the spectrum of
Teff (λ) is not the same for the bright part of a source and the halo. For the latter it
is perhaps close to the spectrum of local sources above flocculi, which are free of
spots.
In maps of radio isophots it can be seen that the sources of the S-component often
extend beyond the optical disk of the Sun. This means that they are situated at a
considerable altitude above the photosphere. At λ ≈ 3 cm the angular size of a radio
emission source is close to the size of the sunspots. This causes the radio emission
to appear before the group of spots comes from behind the Sun’s disk and to
disappear after the group is hidden behind the disk because of the Sun’s rotation.
At λ ≈ 3 cm the interval between the appearance of a local source and a group
9
of spots is about a day. This corresponds to an altitude of the order of 2 × 10 cm
above the level of the photosphere. This interval increases with the wavelength,
reaching two or three days at λ = 10 cm, which can be explained by the increase
in the altitude of the local sources at longer wavelengths. It is possible that the
effect observed is caused by the extent of the sources increasing in heliographic
longitude as the wavelength rises. In the decimetric band the center of a source of
emission located near the central meridian coincides with the corresponding plage
in the decimetric band and with the nuclei of spots at centimetric wavelengths.
The difference in the altitudes at which the local sources and the optical centers of
activity are located leads to a difference in the rates of motion over the Sun’s disk.
The local sources lag behind the optical features connected with them right up to
the moment they intersect the central meridian while after passing through the
meridian the local sources are ahead of the spots and plages.
Despite the fact that the local sources of the S-component are connected with
spots and flocculi, measurements of source altitudes indicate that the generation
regions, being located in the higher layers belonging to the inner corona, are not
identical with these photospheric and chromospheric formations. From optical
observations in continuous light and in the coronal lines it is seen that in the corona
206 ASTRONOMY AND ASTROPHYSICS
above the spots and flocculi there are regions of enhanced electron concentration,
which have a slightly greater kinetic temperature than the ambient medium. The
connection between the local sources of radio emission and the features of the
solar corona in the λ = 5303 Å line and in continuous light has been investigated
by various authors.
Figure 15.3: Comparison of the Sun’s radio isophots at λ = 21 cm with the polar diagram characterizing
the intensity of the K-corona at a distance 1.06 R0 from the center of the Sun’s disk. The
aperture of the coronagraph is shown by the circle in the left-hand figure.
Some results that allow comparison of the features of the coronal emission in
the optical region and the radio frequency band are shown in Figs. 15.3 and 15.4.
Fig. 15.3 shows that the radio brightness distribution at 21 cm above the Sun’s limb
follows the variations in intensity of the K-corona’s glow along the line, which is
shown by dashes. Fig. 15.4 illustrates the presence of a general correspondence
between the radio brightness distribution at λ = 21 cm and the coronal isophots in
the 5303 Å line.
Figure 15.4: Radio isophots of the Sun at 21 cm (a) and isophots of the corona on the 5303 Å line (b). The
radio brightness contours are drawn at intervals of 105° K. The black areas in the right-hand
figure are groups spots; wavy lines are boundaries of plages.
THE SUN’S SPORADIC RADIO EMISSION 207
In recent years the importance of microwave bursts has been recognized due
to their connection with sudden ionospheric disturbances, geomagnetic storms
and cosmic ray variations, bursts of gamma radiation, aurora, etc. The limited
information about this component of the solar radio emission is felt in the absence
of spectrographic observations, which is the usual method of studying the
sporadic radio emission at metric wavelengths. It should be pointed out that the
term microwave (centimetric) bursts only stresses the most characteristic range in
which these events occur. They are also observed at millimetric and decimetric
wavelengths and closely connected with sporadic phenomena in the metric band.
Profiles of typical microwave bursts are shown in Fig. 15.5. Although the bursts
frequently appear in groups their component elements can apparently be discussed
individually so that the grouping of the majority of centimetric bursts is not their
characteristic sign. At rare periods of irregular activity numerous relatively weak
bursts of comparable intensity (Fig. 15.5 a) appear for several hours. The bursts,
which have only one maximum, make up the greater majority of all events. Only
a tenth have a complex structure with two or three or more maxima. In the clearly
defined twin bursts the individual components, which are separated by intervals
of the order of a minute, often show very similar characteristics of duration,
intensity, and size of the source of emission. The microwave bursts consist chiefly
of elementary events of three types: A, B and C.
Type A bursts (simple bursts as shown in Fig. 15.5 b) have the nature of a pulse,
being distinguished by a rapid rise in intensity to a maximum value and a slower
falling away to the pre-burst level. The ratio of the time elapsing from the beginning
of the burst to the maximum to the overall duration of the burst is 0.1 to 0.5 with
a mean value of 0.25 while the duration of a burst is 1 to 5 min. Measurements of
interference have shown that regions with a diameter of about 1 minute are the source
of weak type A bursts with a radio emission flux less than 0.1 from the level of the
quiet Sun. In the case of strong bursts the size of the generation increases slightly to
1.6 minute. The corresponding effective radio emission temperature is of the order of
106° K for weak bursts and of the order of 107 to 109° K for strong bursts. The
angular diameter of the source changes in the process of the development of the
burst and reaches its minimum value at the time of the burst’s maximum phase.
Type B bursts (post-bursts as shown in Fig. 15.5 c) follow an individual simple
burst or a group busts. Their duration varies from a few minutes to a few hours.
The generation region is large in size and the effective temperature of the source
of emission is from 105 to 107° K. Type C bursts (gradual rise and fall as shown in
Fig. 15.5 d) are distinguished by a slow rise in intensity to a maximum value and
a slow fall to the pre-burst level. The ratio of the duration of the first period to the
life of the burst averages 0.4. The duration of these bursts varies from 10 minutes to
several hours. The sources of the bursts are localized in regions of small diameter
with an effective temperature below 106° K. Type A, A–B, and C bursts may be
superimposed on each other or exist separately.
Fig. 15.6 reveals the first dynamic spectrum of an outburst obtained on a radio
spectrograph in the 7.5 to 15 cm band (Kundu and Haddock, 1961). It shows the
recording of the burst’s intensity while observing at a wavelength of 10.7 cm
and the dynamic spectrum of Type III bursts occurred at the same time as metric
wavelengths. The microwave outburst consisted of an intense Type A pulse on
which another peak of low intensity was superimposed. It was followed by a Type
B burst, during which one more Type A simple burst appeared. The spectrum of
the Type B radio emission cannot be seen in the figure due to its low intensity. The
spectrum of the Type A burst clearly indicates that this emission occupies a broad
spectrum of frequencies. Type III bursts at metric wavelengths appeared at the
time of the maximum phase of the microwave outburst. However, the nature of
the dynamic spectra and the duration of these bursts are different. Type III bursts
are thus not a simple continuation of the spectrum of the microwave outbursts into
the metric waveband. In certain cases the impression is created that on the radio
THE SUN’S SPORADIC RADIO EMISSION 209
Figure 15.6: Dynamic spectrum of the outburst of radio emission of July 29, 1959 in the 2000 to 4000
MHz band. At the top is the profile of the outburst at a frequency of 2800 MHz.
emission with a spectrum of the type shown in Fig. 15.7 there is superimposed
broadband emission with a well-defined maximum of spectral intensity located at
Figure 15.7: Profiles (a) and frequency spectrum (b) of the burst of radio emission of April 18, 1959; at
the time of the maximum (11h 40m); at the time of minimum (11h 41.6m); at the time of the
maximum (11h 43.3m); during the post-bursts decay (11h 45m).
210 ASTRONOMY AND ASTROPHYSICS
frequencies where the rise in emission alternates with a region of constant energy
(Fig. 15.8). The presence of the maximum of this broadband emission sometimes
affects the frequency spectrum in the interval where the intensity increases with
the frequency.
Figure 15.8: The same as in Fig. 14.7 (a) for the burst of June 9, 1959; at the time of the maximum
(17h 40m); at the time of attenuation (17h 50m).
It was reported that at λ = 3 cm the emission of two-thirds of all bursts are partially
circular or elliptically polarized. The degree of polarization varies from a few
percent to 30% and in exceptional cases is as much as 50%. Sometimes quasi-linear
polarization was also noted. At λ = 8 cm with a polarimeter whose band was
10 MHz linear polarization was not found in a single burst among more than a
hundred bursts studied. The nature of the polarization of Type B bursts is similar
to the nature of the polarization of the preceeding simple bursts. In the majority of
cases Type C bursts are partially circularly polarized similar to the S-component.
The polarization of outbursts is typical for type A and B events. A study of the
radio brightness distribution over the sources of outbursts has shown that the
generation area consists of two regions: a comparatively small sector with a low
intensity and a long life responsible for the creation of polarized emission, and an
extensive region with a high, strongly varying intensity and a short life.
THE SUN’S SPORADIC RADIO EMISSION 211
Figure 15.9: Radio emission flux and degree of polarization of burst of July 15, 1957 at various
frequencies.
The sources of the microwave bursts are localized in the regions of the Sun’s disk
from which the S-component of the radio emission originates. A connection is
also known between the intensity of this component and the probability of the
appearance of centimetric bursts.
212 ASTRONOMY AND ASTROPHYSICS
Figure 15.10: Radio emission intensity at λ = 9.4 cm and brightness in the Hα line for two bursts and the
flares connected with them (in percentages of the level of the radio emission and brightness
existing before these events). The strong fluctuations in the Hα line in the second figure were
probably caused by the effect of the Earth’s atmosphere.
THE SUN’S SPORADIC RADIO EMISSION 213
is a very clear connection between the radio emission flux in a microwave burst
and the width of the Hα emission line in the flare. A close analogy between the
profile of microwave bursts and the nature of the variation in the width exists for
the great majority of events studied.
The term noise storm combines two closely interconnected events: (i) slow rises in
radio emission intensity of the order of hours and days appearing at frequencies
of the order of hundreds of megacycles, and (ii) brief bursts of the order of seconds
and fractions of a second following rapidly upon one another. It should be noted
that there are bursts of two different types (I and III) among surges of this kind.
Noise storms containing Type III bursts are more rare than Type I storms.
Type I bursts are the most frequent sporadic phenomenon in the radio band.
Each burst appears and stays for the whole of its life in a narrow frequency range
whose width may be much less than the radio emission frequency. The duration of
individual bursts can be obtained from recordings of the radio emission intensity
with high-speed devices having a small time constant. An example of one of these
recordings is shown in Fig. 15.11. Type I surges exhibit a noticeable tendency to
collect in groups of two to three bursts. The interval between two successive bursts
in a group is about 0.5 seconds on the average.
Figure 15.11: High-speed recording of solar radio emission flux during a noise storm.
When a noise storm starts, the fluxes from the solar radio-emission within the
metric band increase by a factor of 102 to 103 or more.
Type II bursts are rarest events in the Sun’s radio emission. Even at a period of
maximum solar activity in every 50 to 100 hours of observations on the average
only one burst of this type is noticed. On the contrary, the frequency of appearance
of Type I and III bursts is one burst per minute. The number of Type II bursts
observed shows considerable variations over a solar cycle and decreasing near its
minimum phase. Type II bursts, in general, appear at the time of chromospheric
flares near the maximum of their radiation in the Hα line or slightly later, and
generally stop earlier than the flare itself. In some cases, the appearance of
bursts is not accompanied by any noticeable effects in the optical region. Flares
of importance Type I or over have been observed in approximately 60% to 80%
of the cases when Type II bursts have appeared. However, the probability of the
appearance of a Type II burst when a flare arises on the Sun’s disk rises rapidly
with the class of flare—from 2% for a Type I flare to 30% for a strong Type III
flare. The eruptive activity of the flares (i.e., ejection of material into the corona)
is not apparently a sufficient condition for the generation of Type II bursts since
the ejection of prominences from the flares is a more frequent phenomenon than
bursts with slow frequency drift. During Type II bursts a rise in the radio emission
–19 –2 –1
flux at metric wavelengths to 10 W m Hz may not be considered out of the
way. Cases are known when the radio emission intensity was even higher. The
data on the nature of the polarization of the Type II bursts are very incomplete but
some say that Type II bursts are either not polarized or are polarized very weakly.
The most outstanding feature of Type II bursts (and also Type III) found when
investigating the dynamic spectra of bursts is the so-called second harmonic. In
the dynamic spectra of the Type II bursts the radio emission occupies two bands
several tens of megacycles wide whose contours repeat each other in general
outline and in many details. Since the structure of one band at the frequency is
repeated in the other band approximately at a frequency 2ω, this served as a basis
for calling the first band of the radio emission the fundamental frequency and the
second band the second harmonic.
EXERCISES
16.1 INTRODUCTION
T
he entire existing terrestrial energy resources owe their origin directly
or indirectly to the Sun. However, the most vivid demonstrations of the
terrestrial response to solar transient phenomena come from magnetic
storms, Short Wave Fade out (SWF), ozone variability, weather connection, etc.
The Sun’s influence on the Earth’s atmosphere manifests itself in various ways
such as radiation (UV, IR and x-rays) and gravitational tides. Variations in
217
218 ASTRONOMY AND ASTROPHYSICS
solar activity may be transients (e.g., solar proton events, solar magnetic sector
boundary passages, solar flares), long period variations (e.g., the 27-day synodic
rotation, 11- or 22-year solar cycle), and short period variations. Several workers
have been taken these variations into account in solar activity and correlated them
with atmospheric phenomena like atmospheric electricity, vorticity, ozone content,
and rainfall. Satisfactory physical mechanisms connecting solar activity with
weather has not yet been evolved and the subject remains highly controversial.
The vital problem is that the energy brought into the atmosphere by processes
associated with variable solar activity is much less than the energy of tropospheric
circulation. The physical mechanisms involved have been exhaustively examined
on ‘Response of the atmosphere to solar activity variations’ with an emphasis on
the role of atmospheric electricity in Sun weather relations. The solar controlled
conductivity variations in the lower stratosphere over thunderstorms control the
current flow in the global circuit. This may ultimately influence cloud physical
processes and thereby the atmospheric energy. The solar modulation of the Earth’s
electric field takes place through regulation of ionizing radiation. The lower
mesosphere and stratosphere appear to be regions where the solar signal interacts
with the lower atmosphere.
Solar radiations can be broadly divided into three categories: direct, diffuse, and
global. Diffuse solar radiation also known as sky radiation represents the short
wave energy of solar origin, scattered and diffused downwards by (i) gas molecules,
(ii) dust particles, (iii) water vapor, and (iv) clouds within the atmosphere. Clouds
are visible symbols of atmospheric activity and these are the seats of atmospheric
electricity. No detailed analysis has yet been made between solar diffuse radiation
and cloud discharges. In fact, the phenomena of scattered radiation in the
SOME UNSOLVED PROBLEMS OF SOLAR METERWAVE RADIO BURSTS... 219
atmosphere are very complex and varied. Numerous factors determining the
quantity of scattered radiation change continuously in time and space thereby
making computation difficult. The amount of diffuse radiation received for a
particular geophysical location on the Earth’s surface depends mainly on the solar
altitude, turbidity in the atmosphere, and cloudiness. For a proper understanding
of the radiation climatology of a region, apart from information on direct and
global solar radiation, a knowledge regarding the distribution of diffuse or sky
radiation in the region is necessary. An assessment of the radiation received from
the sky assumes great importance in problems concerning the utilization of solar
energy, in illumination studies, and in architectural research, especially in the
tropics where the diffuse radiation reaches very high levels. Diffuse solar radiation
is related to electrical discharge from clouds.
correlated with the solar microwave radiations at different frequencies. Thus, solar
radiations, which control the electron density in the D-region of the ionosphere
are related to 10.7 cm solar flux, which is the ionospheric index. The geomagnetic
field plays a remarkable role in controlling the electron D-region of the ionosphere.
It has been reported that the ionospheric region below 85 to 90 km holds extra
ionization during the period of geomagnetic storms. The signal transmission is
affected by this change. Magnetic fields also exist on the Sun. The correlation
between the signal amplitude and the mean solar magnetic field is also a subject
of great interest. Since the Sun is rotating about its axis, its magnetic field strength
observed from the Earth’s surface is likely to change gradually with certain
periodicity. Mention may be made here that the slowly varying component of the
solar radio emission, abbreviated as S-component, emanates from the active region
surrounding the sunspots on the solar disk. Correlations of radio signal with the
solar microwave radiation, geomagnetic activity, and solar magnetic field should
be studied in greater detail to examine many questions yet unanswered.
at about +250 kV relative to the Earth. In the return path of the circuit (i.e., over
the fair-weather portions of the globe), a conduction current flows through the
load resistor part of the circuit. Fig. 16.1 reveals how estimates of resistances in
the global circuit elements are made. An order of magnitude approximation for
the thunderstorm generator is 1500 storms at l A per storm for a total current of
3 8 9
(1.5 × 10 A). The thundercloud tops are at potentials of 10 to 10 v relative to
5 6
Earth and ionosphere, which gives a resistance of 10 to 10 ohm for the charging
resistor over the thunderstorm generator (all thunderstorms in parallel). In the
return flow over most of the Earth’s surface, the load resistor is the ionospheric
potential divided by the total current giving a value of about 165 ohm. It has been
estimated (Markson, 1978) that the resistance under the thunderstorm generator is
104 to 105 ohm. Therefore, the resistance over the thunderstorm generator contains
90% of the total circuit resistance and the one under it contains 10%. The top load
resistance is 0.001% and the bottom load resistance is 0.01% of the total circuit
resistance. The actual values of the resistances so obtained are also indicated in
Fig. 16.1.
Figure 16.1: The electrical global circuit of the atmosphere showing how solar modulation of thunderstorm
electrification might influence the atmospheric circulation. The thunderstorm depicted
represents the totality of all global thunderstorms. Estimated resistances in the global circuit
are also shown.
222 ASTRONOMY AND ASTROPHYSICS
16.6.1 Instrumentation
The equipment used for the continuous registration of diffuse solar radiation is
the same as that for global solar radiation, a Moll-Gorczynski pyranometer and
a Cambridge recording milli-voltmeter with the addition of a Schuepp shading
ring to shade the thermopile element with its two hemispherical glass domes from
the direct Sun. The shading ring arrangement, the installation of the instrumental
equipment, the correction to be applied for the sky inevitably obstructed by the
ring, the standardization of the instruments, and the evaluation records have been
carefully done throughout the work. Regular daily observations with Angstrom
pyrgeometers have been taken at the four stations with some limitations. They
can only be used at night when there is no precipitation and the winds are light.
Precautions have to be taken to avoid the effects of the temperature differences.
The instrument has to be calibrated regularly and despite various practical
SOME UNSOLVED PROBLEMS OF SOLAR METERWAVE RADIO BURSTS... 223
difficulties regular nightly measurements have to be taken at all the stations. Three
Tuned Radio Frequency (TRF) receivers at 10, 20, and 30 kHz were used to record
the atmospherics round-the-clock. The receivers used were designed suitably to
ensure a large dynamic range with a view of handling a fairly wide range of field
intensities.
Figure 16.2: Radiation balance of the Earth-atmosphere system for July over the entire path. The figure
indicates the net radiation flux (× 4.186 × 104) at the top of the atmosphere in J m–2 min–1.
The global average values of albedo, absorbed solar, outgoing long-wave and
net radiation are of the order as shown below in Table 16.1:
Net radiation for July –0.012 × 104 to 0.029 × 104 J m–2 min–1
224 ASTRONOMY AND ASTROPHYSICS
The global distribution of the albedo, the outgoing long-wave radiation and
the radiation balances agree fairly well with the ground-based and air-borne
measurements of the same parameters.
Net radiation Q is the difference between the total incoming and outgoing
radiation fluxes, both shortwave and long-wave and is expressed as,
Figure 16.3: Distribution of sunshine hours during (a) pre-monsoon (b) monsoon and (c) post-monsoon
seasons over India.
sky radiation is found by dividing the individual values by the maximum value
obtained. Similar normalized curve for IFIA is further superimposed by considering
the median noise level obtained from the means at 10, 20 and 30 kHz. The figure
reveals that the trends for the two normalized curves are similar with a little more
fluctuation for the D values. It is also interesting to note that the curve for sunshine
hours is almost a mirror image of the other two curves with a minimum during the
cloudy monsoon season and a maximum during the remaining clear months.
Figure 16.4: Curves showing the sunshine hours, and normalized values of (i) diffuse radiation and (ii) IFIA
for Calcutta.
226 ASTRONOMY AND ASTROPHYSICS
Fading of radio waves due to drifting irregularities in the ionosphere has long been
a subject of great interest. It has been shown that the period of such fading would
increase with scale-size of irregularities, while the period decreases with increasing
velocity of the drift involved. Naturally, if drifting irregularities of different scale-
size ranges exist at the relaxation level, one would expect distinctive periods of
fading. Geomagnetic storms are known to be followed by long-lasting storm after-
effects in the ionospheric D-region at high altitudes. Research has shown strong
enhancement in low frequency ionospheric absorption during and after very strong
magnetic disturbances in the mid-latitudes. The duration of the storm after-effects
were found to correlate with the enhancement of high energy (< 40 keV) particles.
VLF and LF data of such absorptions at low latitude may be examined to see how
they are dependent on season, sunspot activity, and geomagnetic disturbances.
The propagation characteristics will help to suggest an origin of the fading and the
role of atmospheric gravity waves at such times.
Sprites last only a few thousandths of a second and extend from above storm clouds
up to about 60 miles high, reaching the bottom of the ionosphere. Radio noise in
the VLF/LF, MF/HF, and VHF/UHF bands associated with the sprites flashes
makes a distinct pop when played back through a speaker producing a sound that
differs from normal lightening discharge signals. The sprites appearing as blood
red flashes with blush tendrils, appear to have some association with cosmic
radiations including the role of solar radiations. As a result sprites have become
a subject of interest to astrophysicists and space-physicists. A sprite flash can be
one of the strong impulsive radio noise sources that pollute radio-astronomical
observations and it may also create problems in space missions by affecting the
space weather.
in two distinctly different forms: (a) blood red flashes may be called sprites and
(b) blue jets appearing in narrow beams may also be called sprites. The flashes last
only a few thousandths of a second and extend up to about 90 km in height (i.e., at
the bottom of the ionosphere).
Figure 16.5: Photograph of some typical unusual flashes as reported both in the stratosphere and
mesosphere, much above the tropospheric thunder clouds.
Recently, the space shuttle in one of its flights recorded 18 such flashes over
Australia, Africa, South and North America. To capture the position and altitude
of the flashes, special low-light-level cameras were used in aircraft and were
recorded on color video.
but frequently exhibit large peak return stroke currents and may lower tens
to hundreds of coulombs of charge to ground. In most of the cases they are
associated with horizontal extensive spider lightening discharges, long dendritic
channels that may linger for many tens of kilometers along cloud base.
4. Positive CG strokes are disproportionately associated with electromagnetic
Q-burst events in the ELF Schumann Resonance (SR) band, the large amplitude
ringing of the entire Earth-ionosphere cavity. The SR continuous spectra are
believed to be driven by the integrated effects of all global lightening. The
Q-bursts, separated in time by minutes to hours, are assumed to be excited by
the largest lightening events on the planet.
5. The combined data of sprites and ELF transients provide confirmation of
their overall coincidence. For example, in Fig. 16.6, the time of the first video
detection of a sprite is coincident with the transient’s onset to within a few
milliseconds. The time of an associated positive ground stroke, having a peak
current of 327 KA, is also shown. The coincidence of sprites with ELF transients
implying extraordinarily large charge transfer favors an electrostatic over an
electromagnetic trigger mechanism.
Figure 16.6: Record showing sprite related ELF transient. The electric and magnetic fields are represented
by thin and bold markings, respectively. The times of positive CG strokes and sprites are
overlaid.
6. Studies have suggested that sprites are caused by the radiated electromagnetic
fields associated with charge acceleration in intra-cloud flashes or return
strokes (Fig. 16.7).
SOME UNSOLVED PROBLEMS OF SOLAR METERWAVE RADIO BURSTS... 229
Figure 16.7: Schematic diagram showing the proposed connection among sprites, positive CG strokes, and
Q-bursts. The positive CG stroke assumed to be the electrostatic source for the sprites and
electromagnetic source for the Q-burst.
The spectrograph usually consists of three parts: (i) an array of four Log Periodic
13
Dipole Array (LPDA) antenna in the frequency band of 40 to 1000 MHz, (ii) a
spectrum analyzer with a 30 dB preamplifier, and (iii) a PC-based data acquisition
system. Each of the LPDAs essentially consists of two booms and the elements are
attached alternately with the two booms. The separation as well as the length of
the elements decrease in the logarithmic scale. The receiver of the spectrograph
consists of a preamplifier having a noise figure of 3 dB, and a spectrum analyzer.
The spectrum analyzer should have a facility of being operated at various IF
bandwidths, integration time, and scan rate. The frequency range may be chosen
as 40 to 215 MHz, and the spectrum may be available in 400 or 600 data points
on a GPIB port, which is to be connected to a PC having an HPIB card for data
acquisition. The computer takes approximately 10 seconds to access the data
through GPIB port and to write in a data file on the hard disk.
The dynamic radio spectrum resulting from the solar bursts consists of three
important parameters: (i) time, (ii) frequency, and (iii) intensity of the spectrum.
These can be intensity mapped on a grayscale and plotted as functions of time
13
The LPDA is a wide-band antenna, which provides more or less uniform gain over a wideband
of frequency. These are very useful for capturing meter-wave radio bursts.
230 ASTRONOMY AND ASTROPHYSICS
14
and frequency. The data browsing software GRAYLEV is used for this purpose
and divides the screen into an array of cells consisting of nine points. The intensity
is quantized into n + 1 = 10 (0 ≤ n ≤ 9) levels and accordingly any particular cell
is weighted in gray. The steady interferences like TV, radio, wireless, etc., appear
as lines parallel to the time axis. The solar bursts may be roughly visualized as
curved, diagonal, or randomly varying lines across the plot. The software provides
flexibility in intensity scaling, which aids in enhancing the visibility of the trace of
a solar-burst on the plot.
The study of solar meter-wave radio bursts dates back nearly half a century.
During this period the impressive series of advances in radio-astronomical
techniques as well as in theoretical concepts in plasma physics has led to a better
understanding of the different aspects of the solar radio emission in meter-wave
regions, designated primarily as Type I, II, III, IV, and V bursts. There are numerous
unresolved issues about the aforesaid radio bursts, which draw the attention of
various groups of uncompromising researchers in this field. Moreover, due to the
advent of different kinds of space observatories, particularly in the x-ray and γ-ray
bands, the whole scenario about these bursts is likely to be changed. The solar
mission Yohkoh with soft x-ray and hard x-ray imaging instruments launched in
1992, the SOHO (Solar and Heliospheric Observations) space mission with the
white-light/UV coronagraph, and the Compton Gamma Ray observatory will be
useful for comparing the data available from these missions with those of metric
bursts that convey information about the coronal features. Before going into the
different kinds of unsolved problems about metric bursts, let us summarize the
important aspects of these bursts, which are quite well known to us.
1. Type I. These are long-lasting (a few hours to several days) radio emissions,
consisting of thousands of short-lived spikes superimposed on a slowly
varying continuum.
The emission mechanism is thought to be a fundamental plasma emission
and is due to coalescence of Langmuir waves with low-frequency waves
(e.g., ion-sound waves); the exciting agency for the Langmuir waves is a
population of energetic particles trapped in a closed magnetic structure
over an active region.
14
The time comes along x-axis, frequency along y-axis, and intensity can be called grayness of a
particular cell.
SOME UNSOLVED PROBLEMS OF SOLAR METERWAVE RADIO BURSTS... 231
2. Type II. These bursts drift slowly from high to low frequencies, with a
drift rate less than 1 MHz/sec. They often give rise to a double structure,
signifying fundamental and second harmonic emission.
The drift is thought to be due to the outward motion of a magneto-
hydrodynamic shock wave at the leading edge of material ejected as
a result of catastrophic explosion of a flare. The emission is due to the
electrons accelerated during the propagation of shock front.
3. Type III. These bursts are intense and highly drifting from high to low
frequencies, with a drift rate of the order of 100 MHz/sec, which decreases
with the decrease of frequency. These bursts occur over a frequency range
of 1 GHz ≥ f ≥ 10 kHz and are generally repeated at the second harmonic.
They are caused by streams of electrons moving out with a velocity
0.33 c, and their passage sets up plasma oscillation—Langmuir waves—
which then radiate at their characteristic frequency.
4. Type IV. These are long-lasting, non-spiky bursts, which may originate in
stationary or moving sources. The stationary type Type IV are broadband
continua with fine structures, and moving type Type IV are broadband
and slowly drifting with smooth continua.
The evolution of such bursts can be explained with the help of the Gyro-
synchrotron emission mechanism.
5. Type V. On dynamic spectrograms they appear as smooth and short-lived
continua following certain Type-III bursts. They are pronounced in the long-
meter wave region and have longer durations at the lower frequencies.
They are generated due to plasma radiation from electrons trapped in high
coronal loops adjacent to open field lines traced by Type-III electrons.
6. Type U and Type J. In Type U, the observed radiation drifts first downward
and then upward in frequency as time proceeds. U-bursts come from the
leading edge of the CME, which comprises an ejected (or rising magnetic
loop) in which those energetic electrons travel with velocity 0.1 to 0.2 c
that gives rise to radio emission.
Sometimes the bursts turn over (i.e., the drift rate goes to zero) but there is
no emission on the descending legs, and they are called J-bursts.
Apart from the above mentioned types, there are other forms of radio bursts,
such as stria bursts, reverse drift pairs, etc., which have been studied little.
In order to clarify issues arising out of different theories the following points in
regard to the observations are to be emphasized.
Type I
(a) Quantitative data on the dynamic spectra of storm bursts obtained from
high-resolution radio spectrography of good sensitivity are needed. The time
profile, frequency, and frequency drift should be measured accurately by the
spectrographs operating in some selected frequency bands and at suitable
locations.
(b) Simultaneous observations in optical, ultraviolet, x-ray, and radio bands are
needed for getting information on the coronal structures, which are the sources
of noise storms.
(c) As the physical processes involved in stationery Type-IV continuum and Type-I
noise storms are thought to be same, the existence of the link between the two
phenomena should be clarified from observational point-of-view.
(d) The inter relationship between metric and decametric storms should be
studied.
(e) Clarification of the degree to which scattering influences observations of Type-I
emissions is needed.
Type II
(a) The details of acceleration mechanism for electrons at the shock front are not
known till today.
(b) It is still not clear whether the emission comes from both the front and back of
the shock.
(c) The inter relationship between the Type-II exciter and the outwardly moving
mass ejected in the coronal transient is unclear.
(d) It is not clear why some shocks in the interplanetary medium lead to Type-II
emissions, while others do not.
Type III
(a) Quasi-linear acceleration of the electron distribution under the influence of the
Langmuir waves would cause the beam to stop in a very short distance (1 km)
and time (nanosecond). But these values are much too small to account for the
observed Type-III behavior.
SOME UNSOLVED PROBLEMS OF SOLAR METERWAVE RADIO BURSTS... 233
Type V
(a) The hypothesis that Type-V emission is due to plasma radiation from electrons
trapped in high coronal loops adjacent to the open field lines traced by
Type-III electron, has been questioned mainly because Type-V radiation
extends to frequencies of 10 MHz or lower.
(b) It is not clear why Type-V source positions in many cases differ from those of
the preceding Type-III.
(c) It is not clear why Type-V radiation is prolonged.
(d) The cause of the polarization reversal from Type-III to Type-V is not clear.
18. Solarly Submillimetric Telescope 405 MHz & 212 GHz 300 s
(SST) Argentina
Culgoora Radio-spectrograph
The Culgoora radio-spectrograph is being operated in the four ranges of
frequency 18 to 57 MHz, 180 to 570 MHz, and 570 to 1800 MHz designated as
236 ASTRONOMY AND ASTROPHYSICS
A-, B-, C-, and D-bands, respectively. The antenna for the A-band is a dual crossed
log-periodic array with two horizontal and vertical feeds. In the B-band, an open
13 m parabolic dish, with a crossed log-periodic feed, which gives horizontal as
well as vertical components over this band. For the C-band, a solid 5 m parabolic
dish, with a crossed log-periodic feed provides both horizontal and vertical
components. The signals are amplified by using Anzac AMC-162 amplifiers in the
A-band, Anzac AMC-117 amplifiers in the B-band, and Anzac AMC-178 amplifiers
in the C and D bands and are then fed to respective spectrum analyzers (Tektronix
2710) meant for the aforesaid bands. The spectrum analyzers are controlled by a
PC through GPIB. The control PC initiates a sweep of each analyzer in sequence
where one sweep consists of 512 samples, with an 8-bit amplitude resolution, and
takes 0.75 seconds. In this way the full frequency range is sampled once in every
3 seconds. After some processing, the signals are displayed on the main observatory
computer, which is a Sun Microsystems Sparc workstation fitted with an 8-bit
color display card. In order to track the Sun across the sky the antenna drives
are controlled by utilizing an XT PC, which runs a program providing the right
ascension and declination of the Sun. Fig. 16.9 shows some of the antennas of this
spectrograph.
90-degree hybrids for the ranges above and below 1 GHz, a preamplifier with a
noise figure less than 4.6 dB, and several PIN diodes and mechanical switches
controlled by the receiver. In the receiver, the signal is amplified and then up-
converted by a 1000 channel direct synthesizer having the range of frequency
11.2 to 12.2 GHz and the switching speed of 0.45 μs. In order to suppress the
image frequencies in the four GHz bands, four 60 MHz band-pass filters with
center frequencies of 8.17, 9.17, 10.17, and 11.17 GHz are used and one of four
corresponding quartz-stabilized oscillators at 8.1, 9.1, 10.1, and 11.1 8Hz is then
selected to mix the signal down to the second intermediate frequency at 70 MHz,
where it is amplified and reduced to 1, 3, or 10 MHz bandwidth. Fig. 16.10 shows
some antennas of this spectrograph.
The duty cycle of 500 µs is divided into four integrating phases of 106 μs each
and several switching times are used to select frequency, band-width, and signal
source in the front end and A/D conversion. Thus, a scan in frequency containing
1 to 500 observing channels is obtained with a time resolution of 0.5 to 250 μs at
one frequency. The timing of the spectrometer is controlled by an atomic clock.
Nancay Radio-spectrograph
This is a multi-channel radio-spectrograph operating in the frequency range of
110 to 470 MHz, which has been subdivided into 36 bands of 10 MHz band-width.
Each of the 10 MHz band-width amplifiers feeds a set of 10 mixers, with local
oscillator frequencies differing from each other in steps of 1 MHz. The amplifiers
with 1 MHz band-width are connected with 120 mixers out of the 360 signal
outputs which give rise to 120 points in spectrum with 1 MHz resolution. In
order to get higher frequency resolution a further division of frequency can be
made by the simultaneous use of 100 more amplifiers of 200 kHz band-width.
The radio-heliograph presently consists of a cross shaped multi-antenna array.
It has a branch running along the East/West line consisting of 19 antennas. The
238 ASTRONOMY AND ASTROPHYSICS
baseline is 3200 m long. It has a branch along the North/South line consisting of
24 antennas over a length of 1250 m. Fig. 16.11 shows some of the antennas of this
spectrograph.
Weissenau Spectrograph
This spectrograph consists of four parts: (a) log-periodic antennas and a
parabolic dish, (b) structural steel bridge, (c) mechanical drive mount system, and
(d) electronic equipment. The antennas system is operated in the frequency ranges
30 to 46, 46 to 86, 86 to 160, and 160 to 1000 MHz with the help of 2 LPDAs, 4 LPDAs,
8 LPDAs, and 7 m parabolic dish with an LP-feeder made of a pyramidal LP plate
structure, respectively. The antennas are mounted on a steel bridge, supported by
4 legs running on 2 profile wheels on a circular steel bridge. The radiometer is a
swept radiometer. Six radiometer trains operate in parallel over bands identical
with the antenna intervals.
Considering the different frequency scales one realizes that the Weissenau
record does show only the red features of the Tremsdorf record, while the green
features and part of the yellow features do not appear on the WEIS record. Indeed,
the band 40 to 90 MHz seems to reproduce weaker than the corresponding band
on the WEIS record contrary to the band 100 to 170 MHz. So the fundamental band
is patchy on the color plot.
SOME UNSOLVED PROBLEMS OF SOLAR METERWAVE RADIO BURSTS... 239
Gauribidanur Radio-spectrograph
This instrument operates in the band of 40 to 150 MHz. It consists of a T shaped
array stretching 1.2 km along the East/West line and 0.45 km along the South.
The signals from individual antennas are high pass filtered with a cut-off at
40 MHz and fed to broadband amplifiers of a 30 dB gain and a noise figure of 3 dB.
There are 128 log-periodic antennas appearing in 16 groups of 8 antennas along
the East/West direction. The internal spacing between the antennas is kept as
10 m. The Southern arm consists of 64 antennas arranged in 16 groups each
consisting of 4 antennas. The internal spacing in this case is 7 m. All the 32 outputs
from the antennas (16 from the East/West and 16 from the South) are taken into the
receiver room where it is down converted to 10.7 MHz having 1 MHz bandwidth.
This is then digitized and correlated in a 1024 channel correlator system. The
resulting output is stored in digital media. The spatial resolution is 5 × 8 arc-minute
at 150 MHz and the sensitivity is nearly 200 jansky for a 1 MHz bandwidth with
5 sec integration. Fig. 16.12 shows the Western arm antennas of this spectrograph.
EXERCISES
5. Explain how the ionizing radiation emitted from the flare region affect the ionosphere
and the radio communication system.
6. Discuss how the D-region of the ionosphere is affected by solar radiation.
7. Explain with a suitable diagram how the atmospheric electrification is modulated by
solar radiation, pointing out the influence of ionizing radiation.
8. Briefly discuss the various associated parameters of diffuse radiation.
9. Describe a Moll-Gorczynski pyranometer for continuous registration of diffuse solar
radiation.
10. What are the different factors controlling the net radiation? Explain them briefly.
11. What is meant by facing of radio waves? Discuss the effects of geomagnetic disturbances
of solar origin.
12. What are sprite flashes? Discuss various visual reports of sprite flashes.
13. Summarize the observational reports of different type of sprites.
14. Draw a schematic diagram showing the connection among sprites, positive CG stokes,
and Q-bursts and explain.
15. Summarize the important aspects of meter-wave radio bursts.
16. Why does the HiRAS consist of three antennas? Explain the role of each antenna.
17. Describe the operation of the new Culgoora radio-spectrograph and the Zurich ETH
spectrograph.
18. Explain the multi-channel radio-spectrograph of Nancay.
19. What are the four different parts of the Weissenau spectrograph? Discuss.
20. Discuss the radio-spectrograph of Gauribidanur.
Chapter 17
THE LIFE OF STARS IN THE
MAIN SEQUENCE
17.1 INTRODUCTION
F
rom the Hertzsprung-Russell diagram in Chapter 5, it is seen that the majority
of stars are found in the main sequence within the temperature/spectral
classification ranging from Type O to Type M. It is, therefore, understood
that stars spend a large fraction of their time of existence as a main sequence star.
Our Sun, a very ordinary G2 star, is one of them. The Sun must have maintained
its luminosity more or less constant until today. Otherwise, the existence of life
on Earth would not have been possible since the temperature of water on Earth
has to be maintained within limits for evolution of life. The formation of the stars
has also been described in Chapter 5. With the knowledge we have gathered
mainly from the Sun and other stars, we shall examine the physical properties of a
main sequence star. Once we understand the processes that dictate the structure of
main sequence stars, it will be easier for us to investigate more about their
evolution.
241
242 ASTRONOMY AND ASTROPHYSICS
Figure 17.1: A simple diagram for deriving the equations for hydrostatic equilibrium in a star. The star
consists of a small core (where the nuclear reactions take place), which is surrounded by a
sphere of hot gases. The core is much smaller than shown in the diagram.
Since in reality the gas is ionized on account of high temperature, and the
ionized particles and electrons behave as ideal particles, we can use simple
equations. Thus, the gas pressure P can expressed based on the ideal gas law as
244 ASTRONOMY AND ASTROPHYSICS
k ρT
P = . (17.1)
m
The quantity ρ is the gas density, m is the mass, T is the temperature, and k is
the Boltzmann constant. Let us now find the variation of the pressure at various
distance r from the center of the star. The balance between the gravitational pull
and the pressure for ejecting the gas outward is called the hydrostatic equilibrium.
Consider a star of mass M and radius R composed of ideal gas. If Pr is the pressure
at a distance r from the center of the star, and dr is a small elemental length of r as
shown in Fig. 17.1, then the change of pressure with distance from the center of the
star can be expressed as
dPr GMr ρr
= 2
. (17.2)
dr r
The quantity Mr is the mass inside the radius r, G is the gravitational constant,
and ρr is the matter density at distance r. Also, the equation of mass continuity can
be expressed as
dMr 2
= 4π r ρr. (17.3)
dr
It is required to complement the above equations with an equation of state,
which expresses a relationship between pressure, temperature, and density. Since
the gas is assumed to be ideal, the pressure can be expressed as
ρ
P = nkT kT . (17.4)
μmH
The quantity μ is the mean molecular weight of the ionized gas, which is the
average atomic mass per particle, n is the number of particles per unit volume
(n = N/V, where N is the total number of particles in a volume V), and mH is the
mass of a neutral hydrogen atom. A hydrogen atom shall yield one electron and
one proton (two particles). It can thus be written
ρH ρX
nH = 2 2 , (17.5)
mH mH
where X is the mass fraction of hydrogen. Similarly, a helium atom shall yield one
atomic nucleus and two electrons (three particles). The helium atom has a mass of
4mH and thus we get
3 ρH 3 ρY
nHe = , (17.6)
4 mH 4 mH
THE LIFE OF STARS IN THE MAIN SEQUENCE 245
where Y is the mass fraction of helium. Finally, all the elements heavier than helium
are considered metals by the astronomers. In general, a metallic atom shall give off
A electrons as it is ionized and its mass is 2 AmH roughly. Thus, we obtain
ρH 1 ρZ
nm = A . (17.7)
2 AmH 2 mH
The total pressure can be now expressed as
P nH kT nHe kT nm kT
(17.8)
3 1 ρkT ρkT ,
2X Y Z
4 2 mH μmH
where the molecular weight is written as
1 = 3 1
2X Y Z. (17.9)
μ 4 2
The typical values of the mass fractions are X = 0.7, Y = 0.28, and Z = 0.02.
However, the value of Z can vary significantly depending on the class of the
star. The velocity distribution of the particles is generally described by Maxwell-
Boltzmann distribution. In certain stars, the temperature is high enough so that the
random motion of the photons can significantly contribute to pressure. Assuming
that the radiation from a star follows the black-body spectrum, the radiation
pressure Prad can be expressed as
1 4
Prad = αT . (17.10)
3
The quantity α is a constant and is equal to 7.55 × 10–16 N m–2 K–4. Since the
core of the star is much smaller compared to its entire volume, its temperature can
be expressed as
Gmav M
T = . (17.11)
3 kR
The quantity mav represents the total mass M of the star divided by the total
number of particles N in the star (mav = M/N). If we substitute the values for the
Sun in Eq. (17.11) we find that the temperature of the core of the Sun is T ≈ 4 ×
106 K. This is only an approximate value. A more detailed analysis shows that the
temperature of the core of the Sun is T ≈ 1.56 × 107 K.
Figure 17.2: The cross-sections of stars having different masses. (a) A star like our Sun. (b) A star of five
solar masses.
Table 17.1: The life periods of main sequence stars having different masses/radii.
It is believed that after the Big Bang matter consisted mainly of hydrogen and
helium. These elements were converted into different elements in the nuclear
factories of the stars by a process called nucleosynthesis. Unlike the nuclear plants
we have on Earth, stars do not produce energy by splitting atoms through nuclear
fission. On the contrary, they join the atoms together through nuclear fusion. The
requirement for the fusion process to start is an extremely high temperature in
millions of K. In the cores of freshly born stars, the hydrogen atoms fuse to form
the next heaviest element, helium. After the hydrogen gets exhausted the helium
atoms fuse together to form carbon. And when the helium runs out the carbon
atoms fuse to form oxygen. This manufacturing process of heavier and heavier
elements continues until the star finally dies. Our Sun is likely to halt during the
carbon-oxygen stage of fusion. Eventually, the stars that are a few times bigger
than the Sun will have cores made entirely of iron. Elements heavier than iron are
created in supernova explosions.
THE LIFE OF STARS IN THE MAIN SEQUENCE 249
Let us now look more deeply into the nuclear reactions taking place in a main
sequence star like our Sun. There are two basic reactions taking place, namely
the Proton-Proton reaction and the Carbon-Nitrogen-Oxygen cycle as shown in
Fig. 17.3.
Figure 17.3: The temperature-dependent basic reactions taking place in a star. (a) Proton-Proton reaction.
(b) Carbon-Nitrogen-Oxygen cycle.
It must be noted that Eqs. (17.12 a) and (17.12 b) must be done twice for every
(17.12 c). 3He and 4He are stable isotopes. To fuse 4He higher energies are required
since the repulsion between two pairs of protons in helium is even more than the
repulsion between only two hydrogen nuclei.
The carbon acts as the catalyst since it initiates the chain of reactions but does
not get consumed. Nitrogen and oxygen can also start similar reactions. According
to the current estimate, the Sun produces 98% to 99% of its energy by means of
the proton-proton reaction and 1% to 2% from the CNO cycle. It is presumed that
CNO would be the major energy producer if our Sun was 50% more massive.
As said earlier, neutrinos interact weakly with matter and therefore escape
out of the star. Detection of neutrinos would allow us to visualize a portion of
the processes taking placing within the core of a star. Until today, the neutrinos
detected by experiments are much less than predicted in theory. This only suggests
that our understanding of Sun is incomplete. One possibility is that there could be
three type of neutrinos and as they leave the core they may possibly change into
a type that is undetectable using today’s instruments. But definite proof is yet to
come.
It can be seen from Fig. 17.4 that the relative amount of energy release from
the two types of reactions described are dependent on the temperature of the core
of the star. For temperatures below 1.8 × 107 K, the proton-proton chain reaction
releases relatively more amount of energy than the CNO cycle. At a temperature of
1.8 × 107 K, both reactions release an equal amount of energy. Beyond this
THE LIFE OF STARS IN THE MAIN SEQUENCE 251
Figure 17.4: In a main sequence star, the proton-proton chain reaction generates a relatively more energy
for core temperatures below 18 × 106 K. Above this temperature the carbon-nitrogen-oxygen
cycle dominates in releasing energy. Thus, the major energy release mechanisms of the massive
stars are the CNO cycles.
EXERCISES
1. State the Russell-Vogt theorem about the structure of a star and explain.
2. Outline the physical conditions existing in a typical star.
3. What is provided by the equations of the stellar structure?
4. Assuming that the star is composed of ideal hot gases, establish the equations of
hydrostatic equilibrium in a star.
5. Explain the energy transport mechanism inside a star.
6. What is radiation pressure? Discuss briefly the main sequence and life span of the stars,
with a special emphasis on the nuclear reactions inside a main sequence star.
7. Explain the temperature dependent basic reactions in a star considering:
(a) Proton-proton chain reaction.
(b) Carbon-Nitrogen-Oxygen cycle.
252 ASTRONOMY AND ASTROPHYSICS
Chapter 18
THE POST-MAIN SEQUENCE
PERIOD OF THE STARS
18.1 INTRODUCTION
I
n this chapter, we will study what happens to a star immediately after its main
sequence life comes to an end (i.e., immediately after the hydrogen fuel comes
to an end). Some part of it may be scientifically guessed assuming that nuclear
reactions were suddenly suspended. It takes about 1,000,000 years for the energy
of a photon to reach the surface of the Sun. In its way, the original gamma ray
photon loses its energy through interaction with the gases and is finally emitted at
the surface as many ultraviolet and optical photons. Thus, if the nuclear furnace is
suspended today, the luminosity of the Sun would probably remain more or less
constant over a long time by human standards. But that is not the complete story.
The star continues to undergo its next life sequence called the post-main sequence.
This process of evolution is also very much mass dependent. We will learn about
them in the following sections.
253
254 ASTRONOMY AND ASTROPHYSICS
Immediately after the hydrogen in the core gets converted to helium, the main-
sequence life period of a star ends and the first part of the post-main sequence
begins through the red giant phase. During a star’s life, ashes of one nuclear-burning
phase provide fuel for next burning phase. The star expands to a extremely large
diameter and may possibly encroach the nearby orbiting planets. The temperature
of the outer layer falls and the color of the star tend toward red. The star evolves
into the spectral class K and M and their position in the H-R (Hertzsprung-Russell)
diagram changes. This red-giant life of a star mainly consists of (a) a red giant
branch phase, (b) horizontal branch phase, and (c) a asymptotic giant branch phase.
During these processes, the helium in the core of the star and the surrounding
hydrogen shell layers get into fusion with different combinations depending on
the mass, pressure, and temperature, which characterize the temporal behavior of
the star. The surface luminosity, temperature, and diameter of the star fluctuate.
New elements are produced, some of which are thrown into space resulting in
mass loss of the star. We shall study them in detail in the following sections.
Figure 18.3: The giant star expands once again. This time the envelope may possibly increase up to 1.5 AU,
equivalent to the orbit of Mars. The star is now an asymptotic giant branch star (AGB), sitting
on the upper right portion of the H-R diagram (see its phases in Fig. 18.4). The luminosity of
a 1 MΘ AGB is about 10,000 times that of our Sun. An example of an AGM star is Mira, also
known as o Ceti.
Figure 18.4: A theoretical picture of red giant phases of 1 MΘ, 2 MΘ, and 5 MΘ stars. The dark lines indicate
major phases of fusion. The dashed lines indicate uncertain phases. (Adapted from a diagram
by I. Iben, Jr.)
in the thin helium-burning shell up to the surface of the star. These heavier nuclei
can be detected in the star’s spectrum, which provides insight into the undergoing
processes deep within the star. The outer layer of an AGB star loses a mass of
about 10–4 MΘ per year due to increase in pulsation (this is about 100 times more
than in the case of RGB stars). Mass loss is an important mechanism that enriches
the interstellar medium with matter processed in the interior of AGB stars. This
is responsible for the chemical evolution of galaxies and the Universe as a whole.
The production of stardust is the starting point of the formation of all solid bodies
ranging from planetesimals to asteroids and planets. The ejected material is
comprised of mixtures of carbon, oxygen, and other elements. The carbon rich
molecules from the stardust tend to shroud the star in the form of clouds. As the
dust clouds expand, they absorb ultraviolet, visual, and near-infra-red radiation
very efficiently and re-emit them at longer wavelengths (mainly in the range of
5 to 100 μm). Thus, AGB stars are more luminous in the infra-red than at visible
wavelengths. The expanding cloud can also be observed at radio wavebands.
Do all stars become a red giant? The answer is likely to be yes. Up to the red
giant phase, stars with any mass evolve in a more or less similar fashion. The only
Figure 18.5: Stars of different masses reach the red giant (or red super-giant) phase in a similar fashion.
(Adapted from a diagram by I. Iben, Jr.)
260 ASTRONOMY AND ASTROPHYSICS
difference is the size of the giant. In the case of a massive star it evolves as a
red super-giant instead of a red giant. It may be noticed in Fig. 18.4, a star with
1 MΘ enters into the giant area, whereas stars with 5 MΘ or more enter into the
super-giant area of the H-R diagram. Fig. 18.5 shows the phases of stars having
different masses reaching the RGB in a similar fashion. The basic physics is same.
However, the HB phase and AGB phase may differ to some extent. The fate of the
star after the completion of the red giant phase could be a white-dwarf, neutron
star, or black hole depending on its mass.
Let us categorize the stars into three groups based on their Mass: small,
medium, and massive, and summarize their process of post-main sequence
evolution toward the red giant stage. We shall also see some of their properties
in the following sub-sections. Table 18.1 lists the elements involved/produced
through thermo-nuclear burning phases inside stars of different masses.
<< MΘ Hydrogen
temperature. Thus, the core cannot cool. As more and more helium fuses, the more
the temperature rises in the core causing more helium to ignite. Within a few hours,
the temperature shoots up to tremendous amount resulting in the generation of
luminosity of 100 billion stars. But even this energy is not sufficient to blow the
star. The outer shell of hydrogen surrounding the helium also burns. The fusion of
helium leads to the production of carbon. About 5% to 10% of the lifetime of the
star is spent in burning the helium and hydrogen out of the main sequence.
After several 100,000 years, the low mass stars cease the helium-burning
process and expel the hydrogen rich envelope resulting to a birth of nebula. These
are also called planetary nebula since when viewed through telescopes their size
appears to resemble the size of planets like Jupiter. Some of the expelled matter
is chemically evolved, particularly carbon. Remnants of the stars contract and get
heated up to 100,000 K, which radiate ultraviolet photons that are absorbed by the
surrounding gaseous shell.
and subsequently burn carbon to form oxygen, nitrogen, and other heavy elements
down to iron. Since iron cannot be fused to produce still heavier elements since
the process absorbs energy instead of releasing, the nuclear burning comes to an
end. Thus, there is no other force left out to counteract gravitational pull and the
enormous mass collapses toward the core. The nucleus collapses, the intermediate
layers burn out instantaneously, high heat is generated in a short period, and a
catastrophic explosion called a supernova takes place. Some supernova remnants
of stars are listed in Table. 18.2. The gases in the supernova attain temperatures
in billions of K. The explosion also produces very rare heavy materials like gold,
lead, uranium, etc. The star becomes exceptionally bright and is visible for days or
weeks. The ejected materials from the outer layers of the star continue to travel and
expand injecting the chemical elements (which it produced in its lifetime of few
million years) into the interstellar medium. New clouds are formed from which
new stars and planets may be born.
EXERCISES
1. What happens to a star immediately after its main sequence life comes to an end?
2. Discuss the Red Giant Phase of the stars.
THE POST-MAIN SEQUENCE PERIOD OF THE STARS 263
19.1 INTRODUCTION
T
he Universe consists of mysterious objects like white-dwarfs, neutron stars,
black holes, pulsars, quasars, dark matter, and the galaxies themselves. Many
of the rare phenomena occur such as the supernova when a star explodes.
Much of the existing knowledge makes us believe that the stars are the major
driving force in the nature and the activities of the Universe. The majority of the
visible Universe is due to stars. Stars are the factories which generate ingredients
from basic hydrogen and helium. The stars also support life in the Universe like
our Sun, which supports life on Earth. Billions of stars form galaxies, clusters, etc.
The remnants of the stars form white-dwarfs, pulsars, quasars and black holes,
etc., which play important roles in the Universe. But a large portion of the Universe
also consists of invisible dark matter, which has not been clearly understood until
today. In this chapter we will begin with the discussion on the objects produced
out of the remnants of the stars and then study galaxies and dark matter.
265
266 ASTRONOMY AND ASTROPHYSICS
19.2 SUPERNOVA
One of the most brilliant phenomena observed in the sky is the supernova, which is
an explosion producing brilliant light visible for weeks in the sky. It may be thought
as astrophysical fireworks. Supernova is broadly classified into two categories:
Type I and Type II. Type I occurs in a binary system of stars where one is a red
16
giant and the other is a white-dwarf. The latter accretes hydrogen and helium
from the former, generates heat by gravitational compression, and starts a rapid
fusion resulting in an explosion. Type II supernovas occur when a red giant or red
super-giant explodes. The supernova types are further sub-categorized based on
the elements found in the explosion through spectral observation. In general, it has
been observed that in a Type-I supernovas the hydrogen Balmer lines are absent,
whereas in a Type-II supernovas they are present. These are further sub-classified
based on the presence of other spectral lines due to different chemical constituents
that appear in the spectra.
16
We shall read about white-dwarfs later in this chapter, but for the time being assume that a
white-dwarf is an extremely small but dense hot ball whose gravity can pull the gases from
nearby and develop a highly pressurized atmosphere around itself. It is the remnant of a star
like our Sun after all its fuel has been consumed.
17
The region of space around a star in a binary system within which orbiting material is
gravitationally bound to that star. If the star expands past its Roche lobe, then the material
outside of the lobe will fall into the other star.
BRILLIANT PHENOMENA: BLACK HOLES, PULSARS, QUASARS, DARK MATTERS, GALAXIES, AND CLUSTERS 267
Figure 19.1: A re-occur of the nova on RS Ophiuchi as seen on February 23, 2006 from Mt. Laguna,
California.
268 ASTRONOMY AND ASTROPHYSICS
Figure 19.2: Multi-wavelength x-ray image of the remnant of Kepler’s Supernova, 1604. (Chandra x-ray
Observatory.)
White-dwarfs, neutron stars, and black holes come under the category of compact
objects in the sky, which are formed from stars after they have passed through
the red giant phase. The white-dwarf is a hot ball that remains when a star uses
all its nuclear fuel. The mass of the white-dwarf is close to that of the Sun and in
any case can’t exceed the Chandrasekhar limit of 1.4 MΘ. It is a compressed object
containing high density material. However, a neutron star is the more compact with
its density about 109 of that of a white-dwarf, and a much smaller volume radius
of around 12 km. The neutron star could be a rotating or non-rotating type. If the
mass of the compact object exceeds more than 2.5 MΘ, invariably a gravitational
collapse takes place and the star ends up as a black hole.
19.3.1 White-Dwarfs
After the triple-alpha process in red giant stars is complete for those which evolved
with a mass less than 4 MΘ, not enough energy is left to start the carbon fusion
270 ASTRONOMY AND ASTROPHYSICS
process. The cores of these stars finally contain mainly carbon and oxygen. The
core shall gradually radiate away its energy and cool down. No energy is left to
hold the gaseous structure against the gravitational collapse. The star thus becomes
an extremely dense object (typically 109 kg/m3) in the sky with a typical mass of
that of the Sun contained within a volume about equal to that of the Earth, which
is called a white-dwarf. Typically, one teaspoon from a white-dwarf would weigh
5 tons. The higher the mass of the white-dwarf the smaller its size or more densely
packed. In the process of the collapse the star gets heated up and glows brilliantly.
The only support against further compression is the electron degeneracy pressure.
The mass of a white-dwarf cannot exceed the Chandrasekhar limit of 1.4 MΘ.
It would take hundreds of billion years before the white-dwarfs would cool to
temperatures at which they would no longer be visible. Our old Sun is expected to
follow the same pattern to becoming a white-dwarf.
Since white-dwarfs are extremely small they are very hard to detect. However,
in a certain case of the binary system of star, the detection could be eased. There are
other methods also to locate white-dwarfs. The first white-dwarf discovered was
Sirius B (Fig. 19.1). In 1844 Friedrich Bessel noticed that the original Sirius (Sirius
A) had a small back and forth motion, which is more likely if being orbited by an
Figure 19.3: An example of a white-dwarf (picture taken from Hubble space telescope, NASA). The
Sirius A and Sirius B form a binary system.
unseen gravitating object. In 1863, Alvan Clark an optician, would resolve this
mysterious object, which was a very small illuminating object compared to Sirius
A, which was none other than a white-dwarf. Thus, it was named as Sirius B. The
orbital period of this binary system is about 50 years. The Hubble space telescope
has detected more than 75 faint (due to the long distance) white-dwarfs in the
globular cluster M4 situated in the constellation Scorpius. Another white-dwarf,
HZ 43, was observed using the x-ray satellite ROSAT. Many of the planetary
nebulae that can be seen using a moderate-sized telescope contain white-dwarfs
near the center. Some of the bright dwarfs are listed in Table 19.2.
S.N. Name of the star Right Asc. Dec. Mag. Distance Diameter
(light yr) (km)
dense. The electrons combine with the protons to form neutrons under influence
of terrific gravitational pressure due to the mass of the star. Typically, the mass of
a neutron star lies between 1.4 MΘ to 2.1 MΘ and its radius lies between 20 km to
10 km. In comparison to the size of the Sun, it is 30,000 to 70,000 times smaller. The
typical density of neutron stars varies between 8 × 1013 gm/cm3 to 2 × 1015 gm/
3
cm . Since the size of a neutron star reduces drastically from the size it had when
it was in the main sequence, its rotation increases by a terrific amount due to high
reduction in the moment of inertia provided it was rotating in the main sequence
period of its life. Typically, the rotation could vary from 1/1000th of a second to
30 seconds. The escape velocity from a neutron star is roughly half the velocity of
light. Neutron stars not only emit like a black body but also have additional ways
of electromagnetic emission. One variety of neutron star is the pulsar, which is a
rapidly rotating neutron star with an extremely high magnetic field with an axis
inclined away from the axis of rotation. Some of these are responsible for sending
pulsed electromagnetic radiations to Earth.
The structural understanding of the neutron star is purely based on existing
mathematical models, which may undergo revisions in the near future. We really
don’t know the true interior structure. Fig. 19.4 shows the interior structure of
a neutron star based on these models. The atmosphere of the star is about 1 m
thick below which lies a layer of crystalline solid iron crust of about 1 km to 2 km
thickness. As we proceed below we find nuclei with ever-increasing numbers of
neutrons under tremendous pressures, which prevent them from decaying. This
is called in the figure, heavy liquid. Further deeper arrives a point called neutron
Figure 19.4: A model of the interior of a neutron star of 1.5 solar mass and 10 to 20 km radius.
BRILLIANT PHENOMENA: BLACK HOLES, PULSARS, QUASARS, DARK MATTERS, GALAXIES, AND CLUSTERS 273
drip from where onwards there are nuclei, free electrons, and free neutrons until
the core is reached. The material in the super-dense core is not well understood.
Possibilities are that the core material could be a super-fluid mixture of neutrons
with few protons and electrons, or it could contain high-energy particles like pions
and kaons in addition to neutrons, or it could even be composed of some strange
matter incorporating quarks heavier than its up and down varieties, or it could
even be quark matter not bound into hadrons.
As shown in Fig. 19.5, the rotation of the strong magnetic field in a neutron
star converts it to an electricity generator. The electric field thus generated is strong
enough to rip the charged particles such as electrons away from the surface of the
star. These particles move along the magnetic field lines and reach the magnetic
poles at accelerated velocities, which produce narrow but intense beams of
radiation, pointing away from the two magnetic poles. When one of these beams
sweeps our Earth, we receive an electromagnetic pulse, very much like we see
the light from a flashlight only when it is pointing toward us. This is sometimes
called the lighthouse effect. Neutron stars whose electromagnetic beams of radiation
happen to sweep across us are called pulsars. A good example is the pulsar situated
at the center of the Crab Nebula, which rotates 33 times per second. This is shown
in Fig. 19.6.
Figure 19.6: This image of Crab Nebula has been prepared with the superposition of x-ray (blue) and
optical (red) images. The surrounding nebular gases are being stirred by the pulsar’s magnetic
field and radiation. (Image created by NASA.)
The most commonly thought black hole is produced out of dying stars having
a mass greater than 20 times that of our Sun. In a main sequence star, there is a
constant tug-of-war between gravity pulling in the gaseous matter and pressure
pushing them out. The equilibrium between the gravity and pressure is maintained
until the star remains stable. The support for this comes from the nuclear furnace
of the star. But when the nuclear fuel comes to an end, gravity gains the upper
hand and the material in the core is compressed even further. The more massive
the core of the star, the more the force of gravity. Thus, the star collapses under its
own weight to a very small volume. The future of the star depends on its mass. For
small stars like the white-dwarf, after the nuclear fuel is burnt out, the repulsive
forces among electrons within the star eventually create enough pressure to halt
further gravitational collapse. The star then cools and dies peacefully. But in the
case of a very massive star (more than 20 MΘ), after its nuclear fuel has exhausted,
it explodes as a supernova. The outer regions of the star are expelled violently into
space, whereas the core, which if more massive than 2.5 MΘ, completely collapses
under its weight and no known repulsive force within the star can push back hard
enough to prevent gravity from completely collapsing the core into a black hole.
The core then compacts into a mathematical point with virtually zero volume and
is said to posses an infinite density. The escape velocity from a black hole is greater
than the speed of light. Since no object, not even light can attain this speed, it is
said that nothing can escape from a black hole. The distance from the black hole
at which the escape velocity is just equal to the speed of light is called the event
horizon. Anything including light, which passes across the event horizon into the
black hole, is trapped forever. The theoretical radius collapse of a star with a mass
of 20 MΘ as it becomes a black hole is shown in Fig. 19.1. As soon as the radius
reduces below 59 km, no radiation can escape from the star.
276 ASTRONOMY AND ASTROPHYSICS
Figure 19.7: Theoretical radius collapse of a star having 20 times the mass of the Sun. As soon as the radius
reduces below 59 km, even radiation cannot escape from the star, and the object becomes
invisible.
The event horizon, also known as the Schwarzchild radius, rs, is a characteristic
radius associated with any non-rotating mass and is given as
rs = 2 G m . (19.1)
c2
Here, G is the gravitational constant, m is the mass of the gravitating object,
–27
and c is the velocity of light. If m is in kilograms, rs = m × 1.48 × 10 meters. A
non-rotating object, which is smaller than its Schwarzchild radius is a black hole.
The Schwarzchild radius of the Sun is only 3 km. Thus, if the radius of the Sun is
somehow brought below 3 km it would turn out to be a black hole. The surface at
the Schwarzchild acts as an event horizon in a non-rotating body. Fig. 19.8 shows
a simulated picture of a black hole.
Figure 19.8: A simulated picture of a black hole possessing ten solar masses as would be seen from a
distance of 600 km with the Milkyway in the background.
galaxy. Their mass could be equivalent to that of billions of Suns. Our own
super-massive black hole is about 2.5 million solar masses. These are thought
to have formed gradually and not in one collapse of a cluster of stars. The
beginning could be with a stellar-sized black hole, and they grow larger
by accretion of matter and other black holes. As evident in the outer parts
of galaxies, the distances between the stars are huge, whereas in the central
region of the galaxies, the stars are very closely packed. Since everything in
the central region of a galaxy is tightly packed, a black hole in the center of a
galaxy becomes more and more massive as stars orbiting the event horizon
are ultimately captured by gravitational pull and added to the mass of the
black hole. The measured velocity of the orbiting stars close to the center of the
galaxy suggests the presence of a super-massive black hole. Sometimes there
are two jets of hot gases perpendicular to the accretion disk of a super-massive
black hole. The length of these jets could be millions of light years. These
are probably caused by the interaction of gas particles with strong, rotating
magnetic fields surrounding the black hole. The Hubble space telescope
provided the best evidence to date that super-massive black holes exist.
Stellar Black Hole: If 3 MΘ of matter of nuclear density (about 1018 kg/m3)
is accumulated, it will fall inside its own Schwarzschild radius, thus forming
a stellar black hole on its own. This happens in the case of massive stars of
more than 20 MΘ. When such a star runs out of nuclear fuel, it can no longer
sustain its own weight and begins to collapse. The star heats up due to the
collapse and some fraction of its outer layer, which may still contain some
278 ASTRONOMY AND ASTROPHYSICS
fresh nuclear fuel, and activates the nuclear reaction once again and explodes.
This is called a supernova. The remaining innermost portion of the star (i.e.,
the core) continues to collapse. If the mass of the core exceeds 2.5 MΘ, a stellar
black hole is formed for it would end up a singularity with theoretical infinite
density.
Primordial Black Hole: When the mass is small, the formations of miniature
black holes require extremely high density material so as to collapse within
the Schwarzchild radius. The exact mechanisms that result in this type of
black hole have not been precisely identified, but a number of hypotheses
have been proposed. It is thought that miniature black holes might have been
formed shortly after the Big Bang when densities were extremely high. The
rapid expansion of some matter might have compressed some slower-moving
matter enough to contract into black holes. Some hypothesize that black holes
can theoretically evaporate and explode. The mass of the black hole determines
the time required for evaporation. In the case of massive black holes, the time
needed is longer than the current accepted age of the Universe. For a black
hole formed at the time of the Big Bang to evaporate today its mass must be
about 1015 g. During its final phase of the evaporation, such a black hole would
explode with a force of several trillion times that of our most powerful nuclear
weapon. So far there has been no observational evidence for miniature black
holes.
19.3.5 Quasars
Quasars can be thought of as an engineering application of black holes to build
powerful ultra wide-band transmitters in the Universe. In other words, black
holes are the dynamos of quasars. In 1960, Thomas Matthews and Allan Sandage
discovered a faint star-like object of about 16th magnitude, which was named
quasar (quasi-stellar). In the Third Cambridge Catalog, this object is recorded
as 3C 48 and is usually used for various calibrations of radio telescopes. Some
quasars emit extremely broad spectrum starting from radio to x-rays, but many
of them do not emit in radio wavelengths. These objects show a high amount of
18
red-shift, which indicates they are extremely distant objects. These are extremely
18
Redshift is defined as an increase in the wavelength of electromagnetic radiation received by
a detector compared with the wavelength emitted by the source, which usually occurs when
there is a relative motion between the source and the observer with increasing distance of
separation with time like Doppler shift. The opposite of red-shift is blue-shift where the received
wavelength is less than the emitted wavelength. The distance between the source object and
the observation can be determined from the measured red-shift or blue-shift. The increase or
decrease of wavelengths can also occur due to the presence of high gravitational fields, which
are called gravitational red-shifts and gravitational blue-shifts.
BRILLIANT PHENOMENA: BLACK HOLES, PULSARS, QUASARS, DARK MATTERS, GALAXIES, AND CLUSTERS 279
40
bright objects producing luminosity of the order of 10 watts. The radiation is
usually thought to be generated by a giant black hole devouring its surrounding
matter.
The study of quasars has provided some clues to scientists about what
happened to the end of Big Bang’s re-ionization. The oldest quasars (largest red-
shift) show absorption regions in front of them, which suggests that during that
time the intergalactic medium was a neutral gas, whereas the younger quasars
show no absorption region but on the contrary their spectra contain a spiky area
known as the Lyman-alpha forest, which suggests that the intergalactic medium
has undergone re-ionization and formed plasma, and that neutral gas exists only
in small clouds. Another interesting characteristic of quasars is that they show
evidence of elements heavier than helium, which suggests that galaxies underwent
a massive phase of star formation between the time of the Big Bang and the first
observed quasars.
19
In the presence of a magnetic field, accelerated electrons near the speed of light move along
helical paths and as a result, they emit radio waves. It is called synchrotron radiation since these
waves are observed in laboratories when physicists send high energy electrons around in circles
using magnetic fields in particle accelerators call synchrotrons.
280 ASTRONOMY AND ASTROPHYSICS
Figure 19.9: The huge disk in the picture represents an active galaxy. The central region contains the super-
massive black hole. A small accretion disk is shown. The relativistic jet is shown to propagate
out of the center in directions perpendicular to the accretion disk.
sight to Earth, relativistic beaming will change its apparent brightness. Since the
calculated speed of these jets is about 99% of the velocity of light, these are called
relativistic jets. It is hypothesized that these jets carry as much mass as in Jupiter.
An actual picture taken by NASA’s Hubble space telescope of a relativistic jet from
the quasar situated in the M87 radio galaxy is shown in Fig. 19.10.
Many stars along with planets, gases, and dust live together in a gravitationally-
bound system called a galaxy. For example, our solar system along with hundreds
of billions of stars and many other forms of matter make the Milkyway galaxy. It has
been observed from extremely large distances that even galaxies reside in groups
forming a gravitating system called the cluster. The clusters contain the galaxies
as well as many other types of material, which are within the intergalactic space.
All the matter in the Universe is held together in groups by gravitational forces,
which are generated by the matter themselves. The intracluster 20 medium space
20
Intracluster medium (ICM) is the superheated gas present at the center of any galaxy cluster.
BRILLIANT PHENOMENA: BLACK HOLES, PULSARS, QUASARS, DARK MATTERS, GALAXIES, AND CLUSTERS 281
Figure 19.10: Quasar in the M87 radio galaxy emitting a relativistic jet of electrons and other sub-atomic
particles traveling at nearly the speed of light from its center. The blue of the jet contrasts with
the yellow glow from the combined light of billions of unseen stars and the yellow, point-like
globular clusters that make up this galaxy. (Image: NASA Hubble Space Telescope.)
within the clusters is filled with a hot gas that shines in x-rays instead of visible
light due to extremely high temperatures of the order of 107 K to 108 K. From the
temperature distribution, it is possible to estimate how much the gases are being
squeezed by the force of gravity from all the material within the cluster. From
these studies the quantity of matter in that part of space can be estimated.
The results are quite astonishing. It turns out that there are five times more
material in clusters of galaxies than we would expect from the galaxies and hot
gas we can see. The reason is quite logical. All matter does not emit light and
so we can’t see them. But their presence is felt from the indirect measurement.
This invisible matter contributing five times more gravity is called dark matter.
Evidence for the existence of dark matter is also found in the galaxies as their
dynamics of rotation does not tally with the visible mass. We will learn about this
in the following sections.
282 ASTRONOMY AND ASTROPHYSICS
19.5 GALAXIES
Figure 19.11: A majestic spiral galaxy NGC 4414 imaged by the Hubble Space Telescope in 1995. It is
situated in the constellation Coma Berenices, having a diameter of nearly 56,000 light years,
and is roughly at a 60 million light years distance. (Image created by NASA.)
BRILLIANT PHENOMENA: BLACK HOLES, PULSARS, QUASARS, DARK MATTERS, GALAXIES, AND CLUSTERS 283
After the elliptical section the diagram splits into two arms of the tuning fork.
The upper arm contains the spiral galaxies, which start from S0 at the junction of
the tuning fork. The arms of spiral galaxies seem to emerge directly from the core.
The letter S stands for spiral, while the number 0 attached to it means no arms. In
the other cases we do not have numbers but a lowercase letter that denotes how
wound up the arms are. The range of letters is explained below:
(i) Sa. A bright central disk with smooth but tightly wound arms.
(ii) Sb. Slightly better defined spiral arms than Sa.
(iii) Sc. Much looser wound spiral arms than Sb.
(iv) Sd. Very loose arms generating most part of the luminosity.
The lower branch is for barred spiral galaxies. A barred spiral galaxy is defined
as a spiral galaxy with a band of bright stars emerging from the center and running
across the middle of the galaxy. The spiral arms appear as if having emerged from
the ends of the bar in these galaxies. The letter B stands for barred. The subscript
lowercase letters indicate how heavily defined the bar is as given below:
(i) SBa. A bright central region and tightly spaced spirals.
284 ASTRONOMY AND ASTROPHYSICS
(ii) SBb. Relatively better defined arms than SBa, which are more loosely
wound.
(iii) SBc. A much dimmer central portion having even looser arms than SBb.
As was mentioned, our Milkyway Galaxy was earlier believed to be a Sb
galaxy like the Andromeda Galaxy, but today it is believed to be a SBb galaxy. Note
that Hubble’s classification was based on photographs of the galaxies, which made
him initially think that elliptical galaxies were an early form of galaxies, which
might have later evolved into spiral galaxies. Although today’s understanding
suggests that the situation is more or less opposite his early belief left its imprint
in the astronomical world where people still speak of an early type or a later type of
galaxy based upon whether a galaxy’s type appears to the left or to the right in the
diagram. We shall describe the galaxies in more detail in the next sections.
n = 10 1 b (19.2)
a
(b) Lenticular galaxy: Only the S0 and SB0 fall into this category. These
galaxies appears to have disk-type structures with a spherical bulge at the
center projecting out of it, and not show any spiral arms.
(c) Spiral galaxy: These galaxies range can be from Sa to Sd. These galaxies
have a central bulge and an outlying disk consisting of the spiral arms.
The arms are centered about the bulge, and vary from tightly wound (Sa)
to very loosely wound (Sc and Sd). Sc and Sd have less pronounced central
bulges.
(d) Barred spiral galaxy: These galaxies range from SBa to SBd. These galaxies
have a spiral structure similar to the spiral galaxies, but the arms instead of
oozing out from the bulge are projected out from the ends of a bar running
through the bulge. SBa have the most tightly wound arms. Here again,
SBa through SBd refer to how tightly these arms are wound.
BRILLIANT PHENOMENA: BLACK HOLES, PULSARS, QUASARS, DARK MATTERS, GALAXIES, AND CLUSTERS 285
(e) Irregular galaxy: These galaxies are symbolized as Irr and classified into
two categories: Irr-I and Irr-II. The Irr-I type of galaxies show a spiral
structure but are deformed in some way, whereas the Irr-II are those
galaxies which do not fit into any other said category.
The known properties of the galaxies are listed in Table 19.3.
1. Spiral Halo:
&
9
10 – 10
11 8
10 – 10 10
5–250 Population-II 77.0
Barred Spiral Disk:
Population-I
5 13
2. Elliptical 10 – 10 105 – 1011 1–205 Population-II 20.0
8 10
3. Irregular 10 – 10 107 – 109 1–10 Population-I 3.0
Elliptical Galaxies
The size and mass of these galaxies widely varies. They can be as small as
0.1 kpc to over 100 kpc, while their masses may vary from 107 to about 1013 MΘ.
The smallest types among these are called dwarf elliptical galaxies, which may not
be larger than even a typical globular cluster but may contain considerable mass
in the form of dark matter. The motion of the stars within the frame of the galaxy
is dominated by randomness. These galaxies contain few young stars, very little
interstellar matter, and a few open star clusters. They also contain old stars, the so-
called Population-II stars. It has been found that the larger elliptical galaxies have
a system of globular clusters, which indicates an old population. It is believed
286 ASTRONOMY AND ASTROPHYSICS
that in these galaxies the star formation has come to an end. If there is any star
formation, the rate seems to be very low. These galaxies appear yellow-red in
color, which is in contrast to the distinct blue tinge of typical spiral galaxies (which
may be largely from young stars).
As seen from the tuning fork diagram, elliptical galaxies were originally
thought to vary from spherical to highly elongate. E0 is considered to be the most
spherical. As the number increases from E0 to E7 the corresponding galaxies
become thinner and thinner. It is now understood that the majority of elliptical
galaxies are moderately thin and the classification by Hubble is a result of the
angle of observation of the galaxy. The elliptical galaxies have two distinct physical
structures, boxy or giant elliptical galaxy and disky or normal and low luminosity
elliptical galaxies. The boxy structure is a result of the random motion of the stars
within the galaxy, which could be large in a particular direction, whereas in the
disky structured elliptical galaxies, the random motions are nearly isotropic but
appear flattened due to rotation. The dwarf elliptical galaxies are probably not
true ellipticals and many astronomers refer to them as dwarf spheroidal galaxies,
a debatable issue. Elliptical galaxies and the bulges of disk-type galaxies have
similar properties, and are generally regarded as the same physical phenomenon.
Fig. 19.13 shows an image of an elliptical galaxy from NASA.
Figure 19.13: An image of an elliptical galaxy, Messier 105 (also known as NGC 3379), in the constellation
of Leo, from NASA.
The elliptical galaxies are generally found in the cores of galaxy clusters and
in compact groups of galaxies. Recent observations show blue star clusters within
some elliptical galaxies together with other structures. These evidences have
made scientists today think that an elliptical galaxy is formed when two galaxies
BRILLIANT PHENOMENA: BLACK HOLES, PULSARS, QUASARS, DARK MATTERS, GALAXIES, AND CLUSTERS 287
of any type possessing comparable masses collide or merge with each other.
This phenomenon is called as galaxy mergers. It is also thought that many such
galaxy mergers have taken place in earlier times and are still occurring today but
infrequently. Minor galaxy mergers involve two galaxies of different masses and
are not limited to the giant ellipticals. Even our Milkyway galaxy is known to be
digesting a couple of small galaxies currently. Examples of elliptical galaxies are
M32, M49, M59, M60, M87, M89, M105, and M110.
Spiral Galaxies
These galaxies consist of a dense central part (a bulge) attached to arms, which
rotate around its center. The galaxy appears like a spiral from a distance. The arms
form a disk-like structure with a bulge at its center. Fig. 19.11 shows a majestic
spiral galaxy, NGC 4414, imaged by the Hubble Space Telescope in 1995. The bulge
resembles an elliptical galaxy, containing many old stars (Population-II stars),
and possibly a super-massive black hole at its center. The disk is a flat rotating
assembly consisting of interstellar matter, young stars (Population-I stars), and
open star clusters.
In 1964 the first acceptable theory came in from C. C. Lin and Frank Shu. Their
suggestion was that the spiral arms were manifestations of spiral density waves.
They presumed that stars travel in slightly elliptic orbits and the orientations of
their orbits are correlated. In other words, the ellipses vary in their orientation with
respect to one another in a smooth way with increasing distance from the galactic
center. This is shown in Fig. 19.14. It can be understood that the elliptical orbits
come close together in certain areas to give the look of the arms. Thus, the stars
do not remain in a fixed position, which we now see, but they all pass through the
arms as they travel in their own orbits around the galactic center. There are also
alternative theories based on waves of new star formations.
Figure 19.14: Different elliptical orbits of stars around the galactic center are shown, which makes the
galaxy look like a spiral from a distance.
288 ASTRONOMY AND ASTROPHYSICS
The new stars are formed in their bright arms within the disk, which extends
in a logarithmic fashion from the bulge. The disks of these galaxies are generally
surrounded by large spheroid halos of Population-II stars of which many are
concentrated in the globular cluster (orbiting the galactic center). For a long time our
own Milkyway galaxy was thought to be a pure spiral galaxy, but recent research
from the Spitzer Space Telescope suggests that it is a barred spiral. Examples of
spiral galaxies are Triangulum, Whirlpool, Milkyway, Andromeda, and Sunflower
galaxies.
Irregular Galaxies
The irregular galaxies are those which do not fall into the Hubble classification
for galaxies. In other words, these are neither spiral nor elliptical and often have
a chaotic appearance. Neither a nuclear bulge nor any spiral arm can be traced
out from the structure. Collectively, they form about one-fourth of all the galaxies.
It is believed that most irregular galaxies were once spiral or elliptical in nature
but were deformed by the action of gravitation. Fig. 19.15 shows an image of an
irregular galaxy, Messier 82 (also known as NGC 3034 or the Cigar Galaxy), from
NASA.
Examples of irregular galaxies are IC 1613, Irregular Galaxy IC 10, Leo A,
Messier 82 (The Cigar Galaxy), NGC 1569, NGC 3109, Pegasus Dwarf Irregular
Galaxy, Phoenix Dwarf, Sculptor Dwarf Irregular Galaxy, Sextans A, and Wolf-
Lundmark-Melotte.
Figure 19.15: An image of an irregular galaxy Messier 82 (also known as NGC 3034 or the Cigar Galaxy)
from NASA.
BRILLIANT PHENOMENA: BLACK HOLES, PULSARS, QUASARS, DARK MATTERS, GALAXIES, AND CLUSTERS 289
Figure 19.16: An artist sketch of the Milkyway galaxy as viewed from Earth in a dark night with a clear sky.
ecliptic relative to the galactic plane. The Milkyway divides the night sky into two
roughly equal hemispheres, which indicates that our Solar System lies close to the
galactic plane. The diameter of the disk of the Milkyway galaxy spans between
80,000 to 100,000 light years. The thickness of the galactic core is about 1000 light
21
years. It is composed of 200 to 400 billion stars. The Galaxy is thought to be a
large barred spiral galaxy of Hubble type SBbc (loosely wound barred spiral arms)
having a mass of 5.8 × 1011 MΘ. Its integrated absolute visual magnitude has been
21
Note that the term Galaxy beginning with a capital G is generally used as a substitute name for
the Milkyway galaxy, whereas the term galaxy beginning with a lowercase g is used with any
galaxy. This convention has been followed in many textbooks and is also adapted by us except
in headings where Galaxy doesn’t necessarily substitute the Milkyway.
290 ASTRONOMY AND ASTROPHYSICS
estimated to be –21.3. The currently estimated age of the Galaxy based on research
performed in 2004 is 13.6 × 109 years, which is nearly the age of the Universe itself.
Current estimates tell the distance from the Sun to the galactic center as 26,000
± 1400 light years.
The central object of the Galaxy is strongly suspected to be a super-massive
black hole called Sagittarius A*, which is a bright and very compact source of radio
emission. Like many other galaxies, the distribution of mass in the Galaxy is such
that the orbital speed of most of the stars within the Galaxy does not depend
strongly on its distance from the center. Away from the central bulge, the typical
stellar velocity ranges between 210 to 240 km/s. The orbital period of the stars
follows a different law compared to those in our Solar System. The theoretical
explanation suggests the presence of dark matter inside the Galaxy.
The bar of the Galaxy is estimated to be about 27,000 light years in length. It
primarily consists of red stars, which are believed to be very old. A ring called the
5-kpc ring surrounds the bar. It contains a large fraction of the molecular hydrogen
within the Galaxy. Most of the star formation takes place there. If viewed from an
extremely large distance such as from the neighboring Andromeda galaxy, this
ring would be the brightest feature of the Milkyway. The galactic center is shown
in Fig. 19.17.
Figure 19.17: The galactic center of the Milkyway galaxy along the direction of Sagittarius. The primary
stars of Sagittarius are marked in red circles.
It is believed that there are four major spiral arms, all of which begin at the
Galaxy’s center. Each of the arms describes a logarithmic spiral nature having
a pitch of roughly 12°. The arms have been named Perseus-Arm, Norma-Arm,
BRILLIANT PHENOMENA: BLACK HOLES, PULSARS, QUASARS, DARK MATTERS, GALAXIES, AND CLUSTERS 291
Scutum-Arm, and Sagittarius Arm. Additionally, there are smaller arms like the
Orion-Arm containing the Sun. The Norma-Arm also has several other names:
Cygnus-Arm, Outer-Arm, Cygnus-Norma-Arm, and kpc-Arm (only for the inner
region). A sketch of the measured and extrapolated arms of the Galaxy is shown in
Fig. 19.18. Note the position of our Sun in the Orion-Arm. The Solar System takes
nearly 225 to 250 million years to complete one orbit rotation. This amount of time
is called a galactic year. The orbital velocity of our Solar System is 217 km/s.
Figure 19.18: A rough map of the observed (solid lines) and extrapolated (dotted lines) structure of the
spiral arms of the Milkyway galaxy. The position of the Sun is shown in the Orion-Arm.
22
It is thought to have close resemblance with the outlook of our Milkyway galaxy at a very large
distance outside the Galaxy. The NGC 7331 galaxy is often referred to as the Milkyway’s twin.
292 ASTRONOMY AND ASTROPHYSICS
Figure 19.19: A picture of the Galaxy provided by NASA: an artist’s conception of the Milkyway Galaxy as
viewed along its axis.
halo is not understood extensively. Extensive studies of the rotation of the Galaxy
indicate that the halo dominates the mass of the Galaxy, but the material is not
visible, which is called dark matter. As can be observed from the figures, the disk
of the Galaxy is more or less a flat rotating system that contains the Sun and other
young and intermediate aged stars. The formations of stars take place in the disk,
Figure 19.20: The NGC 7331 galaxy, which is often referred to as the Milkyway’s twin. This is how our
Milkyway galaxy might look like to an observer positioned outside our home Galaxy.
BRILLIANT PHENOMENA: BLACK HOLES, PULSARS, QUASARS, DARK MATTERS, GALAXIES, AND CLUSTERS 293
especially in the spiral arms. The Sun is nearly 2/3 of the distance from the center
to the edge of the disk, which is about 25,000 light years by most modern estimates.
The disk also contain atomic hydrogen (H-I regions) and molecular hydrogen (H-
II regions) gas and dust.
Note that today’s astronomical studies suggest that our neighboring
Andromeda galaxy is approaching us at a speed of 300 km/s. This means that the
Milkyway may collide with it in three to four billion years of time. However, the
collision could possibly destroy some of the elements in both the galaxies, but both
of these galaxies will finally merge together to form an elliptical galaxy.
GM
v = 2 1 (19.3)
r a
v = GM (19.4)
r
The center of mass of a galaxy more or less resides in the nuclear bulge of the
galaxy. All the elements rotate around this center. The stars revolve around the
center at a constant speed over a large range of distances from the center of the
galaxy. Fig. 19.21 shows the variation of rotation speed (orbital velocity) of the
elements, stars, gases, etc., as a function of distance from the center of the galaxy
for a general case and the Milkyway. The theoretical curve does not match the
observed one. Instead of reduction in speed with increase in radii, the rotational
velocity stabilizes and remains roughly constant around 220 km/s. It is also not
understood whether the dip in the graph at about 2.9 Kpc is real or just an artifact
23
The center of mass of a system constituted by distributed masses present int. It is a specific point
at which the system’s entire mass behaves as if it were concentrated.
294 ASTRONOMY AND ASTROPHYSICS
introduced different motions in the gases. The flat behavior is also true in the case
of many other galaxies. A typical rotational behavior is shown for some other
galaxies.
There is no direct way to measure the mass of any distant object. All we have
to do is measure some parameter related to the mass, which can be observed from
Earth. Thus, most of the information that reaches us is only from luminous objects
like stars and hot gases. As seen in the graph the rotational velocity does not
decrease in the expected manner, which makes us guess that there are substantial
extra amounts of matter (could be in the form of gases), which we cannot observe
directly, whose gravity is responsible for the behavior in the curve. This is some of
the strongest evidence for dark matter. There are controversies on the issue of dark
matter like the MOND theory, Modification of Newtonian Dynamics. This theory
suggests that the physics of gravity changes on a large scale but, until recently, was
not a relativistic theory.
The early stage of the Universe is believed to have been almost smooth and
featureless. With its evolution, with time, it gradually developed itself and became
organized. The planets revolve around the Sun and form the Solar System. Similarly,
the Sun, along with hundreds of billions of stars, gases, and dust move around the
galactic center of our Milkyway galaxy. As we move up in the scale we find that
galaxies organize themselves in groups both small and large. The groups include
the galaxies and any other material, which is within the intergalactic space. These
are considered to be the densest parts when the Universe is studied as a large
BRILLIANT PHENOMENA: BLACK HOLES, PULSARS, QUASARS, DARK MATTERS, GALAXIES, AND CLUSTERS 295
scale structure and are the largest gravitationally-bound objects. Depending on the
mass, these groups of gravitationally-bound systems can be classified into three
groups: galaxy groups, clusters, and super clusters. These are explained below.
(i) Groups of galaxies. These are the smallest aggregates of galaxies and
typically contain less than 50 galaxies. The estimated mass is about
1013 MΘ within a diameter of 1 to 2 Mpc. The group to which our Milkyway
belongs is called the local group and contains more than 40 galaxies. The
spread of velocities for the individual galaxies is approximately 150 km/s.
Note that this definition should be used as a guide only, since larger and
more massive galaxy systems are sometimes classified as galaxy groups.
(ii) Clusters of galaxies. Still larger than the groups are clusters containing 50
to 1000 galaxies, hot x-ray emitting gas, and large amounts of dark matter.
The distribution of these three components is roughly the same within
the cluster. The estimated mass is between 1014 MΘ to1015 MΘ within a
diameter of 2 to 10 Mpc. The spread of velocities for the individual galaxies
is approximately 800 to1000 km/s. Fig. 19.22 shows a picture of a massive
galaxy cluster, Abell 1689.
Figure 19.22: One of the most massive galaxy clusters observed is the Abell 1689 in the constellation Virgo.
The yellow galaxies belong to the cluster itself, but the red and blue distorted streaks are
background galaxies, which are being gravitationally lensed by the cluster. Some of the lensed
galaxies are more than 13 billion light years away. The lensing zone itself is 2 million light
years across.
296 ASTRONOMY AND ASTROPHYSICS
(iii) Super clusters. An even larger scale of visualization is called a super cluster
consisting of groups, clusters, and even some isolated galaxies. It was
once thought of as the largest structure in nature, but today it is known be
subordinate to the enormous walls or sheets sometimes called super cluster
complexes. Super clusters range between 10 Mpc to 150 Mpc. Fig. 19.23
shows a map of the nearest super clusters.
Fig. 19.24 shows light from a distant galaxy reaching the Earth after getting
bent by the gravity of a massive object. The orange arrows illustrate the apparent
position of the background source. The white arrows illustrate the path of the
light from its true source location. The path of the light from the source is curved,
which distorts the image, and hence, the apparent location of the source could be
different from its true location. Additionally, multiple images of a single source
might be observed. Note that if the observer, source, and massive object join in a
straight line, the image of the source will appear like a ring behind the massive
object. This phenomenon, usually referred to as Einstein ring, was first mentioned
in 1924 by Orest Chwolson and later quantified by Einstein in 1936.
Figure 19.24: Light from a distant galaxy reaching the Earth after getting bent due to a massive object. The
orange arrows illustrate the apparent position of the background source. The white arrows
illustrate the path of light from its true source location.
The effects of gravitational lensing may be noted in Fig. 19.23 and more clearly
in Fig. 19.25. In the latter figure, a cluster of galaxies called Abell 2218 is shown,
which acts as a powerful gravitational lens, magnifying all galaxies lying behind its
core. The gravitationally lensed galaxies appear to be stretched along the cluster’s
center and some of them are imaged multiple times. The multiple images usually
appear in pairs along with a generally faint third counter image, as is the case for
very distant objects.
298 ASTRONOMY AND ASTROPHYSICS
Figure 19.25: Abell 2218, a cluster of galaxies nearly 3 billion light years away in the constellation Draco. It
acts as a powerful gravitational lens, magnifying all galaxies lying behind its core. The lensed
galaxies appear to be stretched along the cluster’s center and some of them are imaged
multiple times.
The gravitational lensing may be classified into three categories: strong lensing,
weak lensing, and micro lensing depending on the characteristic properties of the
image. These are listed in Table 19.4.
Strong The distortion of the background object is easily visible, i.e., formation of
Einstein rings, arcs, and multiple images.
Weak The distortions are much smaller and can only be detected by analyzing large
numbers of objects. The lensing action shows up statistically as a stretching
of the background objects located perpendicular to the direction to the center
of the lens.
Micro The distortions are not visible in shape, but the amount of light received from
a background object varies with time. The background source can be star, a
quasar, or any distant emitting object within or outside our Galaxy, while the
lens can also be star, a quasar, or any massive object coming in between.
BRILLIANT PHENOMENA: BLACK HOLES, PULSARS, QUASARS, DARK MATTERS, GALAXIES, AND CLUSTERS 299
19.8 RED-SHIFTS
When the wavelength of visible light increases due to some physical phenomena it
is called red-shift. In other words, if the frequency of the light observed is less than
that of the transmitting source, it is called a red-shift. This is shown in Fig. 19.26.
The opposite of red-shift is blue-shift where the observed frequency is more than
the source frequency (wavelength decreases). The terms red-shift and blue-shift are
not restricted to use for visible light alone but are also used for any electromagnetic
waves. These may be quantified by the relative difference between the observed
and emitted wavelengths (or frequency) of an object. In astronomy this change is
measured relatively such that it is dimensionless. If λ and λ0 are the measured and
emitted wavelengths then the red-shift z is given as in Eq. (19.5). The red-shift z
can also be expressed in terms of frequency as shown in Eq. (19.6), where ν and ν0
are the measured and emitted wavelengths. If z is positive it is a red-shift and if z
is negative it is a blue-shift.
λ λ0 λ
z = 1 (19.5)
λ0 λ
ν0 ν ν0
z = 1 (19.6)
ν ν
Figure 19.26: The top portion shows a spectrum of light containing absorption lines (dark lines). The
bottom is the spectrum obtained after undergoing a change (increase) in wavelength. Note
that the lines have moved to the right.
z = 1 v γ 1 (19.7)
c
γ = v2
1 (19.8)
c2
v cos θ
z = 1 γ 1 (19.9)
c
(ii) Cosmological Red-Shift. If we consider the Universe to be expanding
then the distance between any two points for a particular instant will be
different at a different instant. This suggests that the space between these
two points is expanding with time. If we take multiple points, the same
principle applies and the distances between individual points will increase
with time. If the velocity of light were fixed, then due to the continuous
expansion with time a light source at one point as seen from any other
point will constantly change its frequency with time. The frequency will
be reducing with time. In other words, a red-shift will be observed. Since
the distance is a function of time in this case, if we know the red-shift
then we may also obtain the distance of separation. This type of red-shift
is also known as cosmological red-shift or Hubble red-shift. On the other
hand, if the Universe was contracting then the blue-shift would have been
observed. The expression for cosmological red-shift is given in Eq. (19.10),
where aT1 and aT2 are the values of the cosmic scale factor at times T1 and
T2, respectively. A red-shift occurs if T2 > T1 and a blue-shift occurs if
T2 < T1.
aT
z = 2
−1 (19.10)
aT
1
24
Speed comparable to the velocity of light.
BRILLIANT PHENOMENA: BLACK HOLES, PULSARS, QUASARS, DARK MATTERS, GALAXIES, AND CLUSTERS 301
1
z = 1 (19.11)
2 GM
1 2
rc
v = H0D. (19.12)
EXERCISES
3. What is a supernova? What are its two categories and what are their characteristic
differences?
4. Discuss the properties of a Type I supernova.
5. Give the names of different novas detected from 1970 onwards.
6. Discuss the properties of a Type II supernova and point out the characteristic differences
with Type I.
7. What is the order of energy released by a Type II supernova as compared to that of the
Sun?
8. What are compact objects in the sky?
9. Write a few lines about white-dwarfs, neutron stars, and black holes.
10. Write down the important properties of bright white-dwarfs.
11. How is a neutron star formed? What is the value of its typical density?
12. Give a model of the interior of a neutron star.
13. What are pulsars? Write down their important characteristics.
14. Give the basic structural model of a pulsar.
15. From the available theory, write down the basic properties and structures of black
holes.
16. How are the black holes formed? Give a clear interpretation of the same.
17. What is the ‘Schwarzchild radius’? What would be the value of this for our Sun?
18. Depending on the Schwarzchild radius, categorize the black holes into three types and
discuss them.
19. What are quasars? Write down the characteristics of quasars.
20. Explain the mechanism for emissions from quasars.
21. What is dark matter? Give evidence for the existence of dark matter.
22. What is the typical number of stars in a galaxy? What is intergalactic space?
23. Discuss the different physical structures of galaxies. Explain the characteristics of
elliptical and spiral galaxies.
24. Classify galaxies based on the Hubble sequence.
25. Make a general classification of galaxies and explain the different types.
26. Write down the known properties of the galaxies.
27. Mention the important characteristics of elliptical galaxies. Compare them with spiral
galaxies.
28. What are irregular galaxies? Describe briefly our Milkyway, and then galaxies of
different types.
29. What is meant by a rotation of galaxies?
30. Depending on the mass, classify the galaxy groups and explain them.
31. What is meant by a ‘red-shift’? Write down its properties.
32. Explain the various mechanisms for red-shifts and blue-shifts.
33. Write down the Hubble’s law in relation to the expanding universe.
Chapter 20
PROBLEMS OF RADIO INTERFERENCE
AND NEW DEVELOPMENTS IN
RADIO ASTRONOMY
20.1 INTRODUCTION
I
n this chapter we discuss the problems of radio interference, spectrum
monitoring instruments, methods of terrestrial spectrum analysis, and some
new developments taking place in the field of radio astronomy. Spectrum
management is generally not covered in other texts of radio astronomy. We felt its
importance and have treated the subject seriously. Finally, we describe an ongoing
international project for futuristic radio astronomy involving several countries
called the Square Kilometer Array.
303
304 ASTRONOMY AND ASTROPHYSICS
our nearest star, the Sun. Its flux is measured in units of SFU (Solar Flux Units),
where 1 SFU is equal to 10–22 W m–2 Hz–1. A distant radio station may produce a
signal flux density of the order of several millions of SFUs. The radio power flux
densities from distant radio galaxies received on Earth are generally of the order
25 –26
of several micro-jansky to few milli-jansky, where 1 jansky is equal to 10 W
m–2 Hz–1. The astronomical sources produce signals of the nature of a random
noise (Gaussian noise) generally in both polarizations and spread over a wide
band of frequencies. When seen at the output of a radio telescope receiver, the
spectrum appears to be more or less flat over a substantial amount of bandwidth
(unless absorption lines are present). It appears like a locally terminated spectrum
analyzer with a matched load (without connecting any signal source). A locally
generated wideband radio noise like that produced from an electrical arc welding
may be quite high in magnitude to obscure the radio astronomy signal. Any kind
of signal or noise (other than that produced from astronomical sources) interfering
or obscuring the astronomical signals is called RFI (Radio Frequency Interference)
by radio astronomers.
RFI can enter into the radio telescope data through the side-lobes of the
telescope-antenna or as a leakage through the cables/shielding of the electronic
receiver system. This is one of the reasons that radio telescopes are located in
very remote places where human population is rarely present. Fig. 20.1 shows
the various possible ways by which the RFI can directly enter into the telescope
Figure 20.1: Methods by which the RFI gets coupled to the radio astronomy data through the antennas.
(a) RFI signals enter through side-lobes when the antenna points to zenith. (b) RFI signals enter
through spill-over and side-lobes when the antenna points at a small angle to zenith. (c) RFI
signals enter through spill-over, mesh leakage, and side-lobes when the antenna points at a large
angle to zenith.
25
This unit is named after the discoverer of radio emission from extra terrestrial objects by Sir Karl
Jansky in 1933.
PROBLEMS OF RADIO INTERFERENCE AND NEW DEVELOPMENTS IN RADIO ASTRONOMY 305
data through the antenna in low frequency radio telescopes (VHF/UHF). In part
(a), RFI signals enter through side-lobes when the antenna points to zenith. In part
(b), RFI signals enter through spill-over and side-lobes when the antenna points
at a small angle to zenith. In part (c), RFI signals enter through spill-over, mesh
leakage, and side-lobes when the antenna points at a large angle to zenith.
RFI can be broadly classified into two categories: narrowband RFI and
broadband RFI. Fig 20.2 shows the spectrum of narrow and broadband RFI.
Narrowband RFI is generally produced from communication transmitters. In
general, the communication signal may be represented in terms of weighted sine
and cosine functions of frequency (by applying a Fourier transform on the time
domain signal). The spectral power distribution of these signals lie concentrated
around some discrete frequencies represented by sine and cosine terms. Thus,
the overall spectrum appears as tall narrowband lines (communication signals)
standing over a wide uniformly distributed power over the entire frequency band
(radio astronomy signals). In order to process the acquired astronomical data, a
simple technique generally adopted is to remove the data containing these
narrowband lines and process the rest of the data.But this process is not very
Figure 20.2: Examples of narrow and broadband RFI. The power spectrum is plotted against time and
frequency. (a) Three-dimensional plot of raw RFI spectrum. (b) Contour plot of the same
data.
306 ASTRONOMY AND ASTROPHYSICS
efficient if a large portion of the data is wasted. There are other techniques to
eliminate these lines, which are called RFI mitigation. The process mathematically
separates the lines from the original data. A separate antenna is used to collect
the signature of the RFI, which is properly scaled and vectorially (amplitude and
phase) subtracted from the original data. Thus, the remaining data is free from
RFI. However, not all types of RFI can be eliminated by this process. Moreover,
proper signature of the RFI in quality and quantity are needed for this process.
The other type of RFI, as stated earlier, is broadband RFI. These are generated
from man-made electrical processes where a large amount of current in electrical
circuits is switched. In other words, di/dt is large, where i represents the current
and t the time. Therefore, if any amount of inductance is present in the circuit
the induced voltage will be very large and the system can radiate. Examples of
such processes are electrical arc wielding, electric ignition in motor vehicles, SMPS
(Switched Mode Power Supply), etc. The spectral distribution of these signals is
extremely large, which may be understood by looking at the Fourier transform of a
pulse. Thus, the presence of such signals can radio-jam the telescope and the entire
extra-galactic data may get spoiled. This kind of RFI could be more dangerous
than the narrowband type since these do not leave any scope for data usage.
26
Giant Meterwave Radio Telescope. It is a low-frequency radio telescope (50 MHz to 1620 MHz)
located in Narayangaon, Pune, India.
PROBLEMS OF RADIO INTERFERENCE AND NEW DEVELOPMENTS IN RADIO ASTRONOMY 307
Figure 20.3: A picture of the antennas of the RFI monitoring system of GMRT. Four log-periodic antennas
are mounted on a tower of 20 m height in the E-plane directed along East, West, North, and
South directions.
amplified suitably by an LNA (Low Noise Amplifier) and fed to the spectrum
analyzer kept inside the electronics laboratory building through an RF-cable. The
traces of the spectrum analyzer are re-displayed and saved by the computer with
time and directional stampings using its parallel port to GPIB interface. Settings of
the spectrum analyzer and observation period are user defined and set using the
computer. The computer cycles the antennas over the entire period of observation
and keeps recording the data in five different data files.
27
Graphical User Interface of the software through which the operational parameters are set and
the results are displayed.
308 ASTRONOMY AND ASTROPHYSICS
Figure 20.4: Block diagram of the omni-directional RFI Spectrum Monitoring System of GMRT.
holding all the four spectrum data. The software is accessible from the local area
network. Fig. 20.5 (a) shows the GUI objects of the system running the software.
At the back is the system settings window. The details of the spectrum analyzer’s
traces are visible on the top GUI. It also displays the antenna selected, time
progress bar, remaining time, etc. Fig. 20.5 (b) shows the GUI of a data visualizing
software, which can be operated from any terminal connected to the local area
Figure 20.5: GUI of the software used. (a) GUI of the system operating the software. (b) Data visualizing
and radio direction finding software.
PROBLEMS OF RADIO INTERFERENCE AND NEW DEVELOPMENTS IN RADIO ASTRONOMY 309
network. This software tool is used for viewing the completed data files and can
also be directly connected online to the system operating software for visualizing
current data. The East, West, North, and South spectrum (one cycle) is displayed
in four colors. A crosswire can be positioned by the mouse on the spectrum to
visualize the power levels from all four directions for any particular frequency.
The corresponding radio direction appears in a vector plot GUI sitting on top.
This innovative method compares the RF-power levels with the antenna radiation
patterns and finds the radio direction. It can generate time-averaged plots from
completed data and also visualize the data trace-by-trace over the observed time.
Frequency Span
No. of Spectral Channels = (20.1)
Δν res
TspecAna = Noise Floor TSpecAnaPhy (20.2)
k Δν res
28
This is done in the data processing after the observation is complete.
310 ASTRONOMY AND ASTROPHYSICS
Table 20.1: Used specifications and characteristics of the devices and instruments used.
Log-Periodic Dipole Array GAnt = 4.5 dBi, Frequency Range: 50 MHz to 1800 MHz (see
Table 20.2 for antenna off-source temperatures)
Spectrum Analyzer HP-8590L, Δνres=1 kHz, Δνvid =1 kHz, Noise Floor = −110 dBm
RF-Cable-2 LMR-400, Length = 120 m (see Table 20.2 for insertion loss)
1
τ = (20.3)
2Δ ν vid
TspecAna
TspecAnaSys = (20.4)
GLNA εcable1 εcable 2 ε RFswitch
1
TRFswitchSys = TRFswitchPhy 1 (20.5)
ε RFswitch
T 1
Tcable1Sys = cable1Phy ε 1 (20.6)
ε
RFswitch cable1
Tcable 2 Phy 1
Tcable2Sys = 1 (20.7)
GLNA ε RFswitch εcable1 εcable 2
T (NFLNA 1)
TLNAsys = LNAphy (20.8)
εcable1 ε RFswitch
Tsys
ΔTmin = (20.11)
τeff Δν res
2
Table 20.2: Sys. Char. at −180 dB [watt/m /Hz], Δνres= Δνvid =1 kHz, Nsweep = 100.
SNR (dB) 70 65 61 56 50 47 43
between the two antennas. The antennas are kept just outside the arc of the circle.
Thus, the arc-length of the semicircle between the two antennas is R αv, which
may be roughly approximated to (L + D) if αv ≤ 90°. Thus, R may be expressed as
in Eq. (20.14). Here, αv has been named here as the viewing angle, since the radio
sources to be resolved should fall within this angle. The maximum value of αv is
αmax and is determined by the maximum side-lobe level in the radiation pattern, at
the operating frequency. As shown in Fig. 20.7, a circle is drawn on the major lobe
with a radius equal to the maximum side-lobe level. The angle subtended by the
cuts on the major-lobe is 2 αmax.
1
L3 2
D ≥ 0.62 (20.13)
λ max
L D
R ≥ (20.14)
α
v
Figure 20.6: The simplest two-antenna system for radio-direction finding. Two identical log-periodic dipole
arrays are positioned on the arc of a circle of radius R. αv is the viewing angle within which the
radio source must be present, αs is the source angle, L is the length of the largest element of
any antenna, and D is the minimum distance of separation.
PROBLEMS OF RADIO INTERFERENCE AND NEW DEVELOPMENTS IN RADIO ASTRONOMY 313
From Fig. 20.7, considering Axisi the reference axis, let the gain patterns of
the adjacent antennas (i.e., Antennai and Antennaj) within the angle αv subtended
by them be Gi(α) and Gj(α), respectively. Let their ratio be denoted by GR(α) as
expressed as in Eq. (20.15). It must be noted that each point of GR(α) within αv has
a unique value. The antennas deliver output powers Pi and Pj corresponding to a
radio source falling within αv at the source frequency. The angle of the source can
be obtained from Eq. (20.16).
Figure 20.7: Determination of the maximum angle of separation αmax between two adjacent identical
antennas. A circle having a radius equal to the maximum side-lobe level cuts the major-lobe
at two points. These points are joined to the center with straight lines. The angle subtended
is 2 αmax.
G (α )
GR (α) = i , 0 α αv (20.15)
G j (α )
ij
P
αS = GR1 i (20.16)
Pj
In order to increase the effective viewing angle, more antennas could be
positioned on the arc of radius R at angles 2 αv, 3 αv, 4 αv, etc., with respect to
Axisi until the entire 360° (2π of arc) is covered. The number of antennas required
is expressed in Eq. (20.17). However, if N is not an integer, αv can be reduced so
314 ASTRONOMY AND ASTROPHYSICS
that N takes the value of the next possible integer. An RF-switch of the type SPNT
can be used to multiplex the antenna outputs into the spectrum analyzer. The RFI
monitoring system of GMRT (see Fig. 20.4) is a practically working system with
N = 4. For resolving a radio source, the RF-switch is swept from Antennai to AntennaN
and correspondingly RF-powers Pi to PN are obtained, respectively, for any
particular frequency of observation. The source lies within the angle subtended by
the two adjacent antennas, which deliver the highest and second highest powers.
Thus, for a multi-antenna system αs can be expressed as in Eq. (20.18).
N = 2π (20.17)
αv
1 P
αs = (i 1) α v GR i , i 1, 2, ... N (20.18)
Pj
The power flux density (power per unit area per unit of bandwidth) available
on the circle of radius R at the angle αs with respect to Axisi in the direction of the
radio source shown in Fig. 20.6 can be expressed as in Eq. (20.19). With an antenna
gain at αs in dBi, power at antenna terminals in dB[watt], frequency of operation
in MHz and resolution bandwidth in kHz, and the flux density in dB[W m–2
Hz–1] can be expressed in Eq. (20.20). For direction finding, the signal from both
the adjacent antennas should be above the noise floor for detection. In the worst
case the signal may be along one particular antenna axis so that the gain ratio of
the antennas is at its maxima. Under the said condition, the signals received from
both the antennas will be visible if the minimum signal-to-noise ratio of the system
at that frequency is expressed as in Eq. (20.21). This is sometimes achieved by
spectrum frame averaging of a large quantity of data.
4 π Pm
S = 2 (20.19)
Gm (χ) λ Δν res
χ α x if m i , χ (α v α s ) if m j
the formula for modeling wide-beam patterns as expressed in Eq. (20.22), where
B0 represents the maximum value of antenna gain and κ is the angle scaling factor.29
Thus, with reference to Fig. 20.4, the power ratio for a source under observation
can be expressed as in Eq. (20.23). From this the source angle is derived as in
Eq. (20.24). The angular accuracy of the computed source angle due to a deviation
in the radiation pattern of one of the antennas from its ideal value is analyzed
here. The relative change in source angle with respect to viewing angle due to
a small change in antenna output power ratio is expressed in Eq. (20.25). This is
obtained by differentiating Eq. (20.24) with respect to ρij. The usable portion of the
averaged E-plane radiation pattern of the four antennas for the system described
earlier is shown in Fig. 20.8. The side-lobes and back-lobe reside within the shaded
area and are of less concern to us since their gains are well below the operating
region. Eq. (20.22) fairly models this pattern with n = 6, κ = 0.69 and B0 = 2.82
(4.5 dBi). Table 20.3 lists the percentage error in source angle relative to viewing
angle αv due to relative deviation of ρij for a gain radiation pattern modeled for
Fig. 20.7 using Eq. (20.22). It is observed that the error at half viewing angle is
nearly one and a half times of that at its minimum and maximum values.
n
G(θ, φ) B0 cos ( κθ);
(20.22)
π
0 θ , 0 φ 2π, n κ
2
n
cos ( κ α s )
ρij
cos κ (α v α s ) (20.23)
P
also ρij i
Pj
1
1 1 1
αs = tan ρij n cot ( κ α v ) (20.24)
κ sin ( κ α v )
1 dρ
1
ρ ji n
ji
n α v κ sin ( κ α v ) ρ ji
α 2
1
1
1 ρij n cot ( κ α v ) (20.25)
sin ( κ α v )
dα s
also α
αv
and d ρ ji ρ ji Actual
ρ ji
29
For low-gain antennas having wide major-lobes, (first nulls appear after 90° on either side).
316 ASTRONOMY AND ASTROPHYSICS
Figure 20.8: Usable region of the averaged E-plane power radiation pattern of the antennas used in the
system shown in Fig. 20.4.
Table 20.3: Percentage error in αs relative to αv due to relative deviations of ρji for a gain
radiation pattern model with n = 6, B0 = 2.82, κ = 0.69, αv = 90° using Eq. (20.25).
Fig. 20.9. The narrowband spectral lines are separated as shown in the bottom half
of the figure using some computer algorithms. Note that, except where the spectral
lines are present, all other portions are removed. For each spectral line, the largest
and the second largest powers are determined and the corresponding antennas are
identified. If they come from adjacent antennas, the directions and absolute fluxes
are determined using Eq. (20.18) in combination with Eqs. (20.24) and (20.20),
respectively, otherwise the data is rejected. This is done using software specially
Figure 20.9: Averaged spectrums from the four-antenna system of Fig. 20.4. Top: Raw spectrum data from
different directions. Bottom: Narrowband signal spectrums separated from raw data, Δνres =
1 kHz, Δνvid = 1 kHz, Nsweep = 100.
prepared based on the equations. The resolved radio source angles and flux
densities over the entire frequency band are shown in Fig. 20.10. The impedance,
radiation pattern, etc., of the log-periodic dipole antennas repeats with a logarithm
of frequency. For the system shown in Fig. 20.4, the maximum antenna gain
deviation from the averaged pattern (shown in Fig. 20.8) is within 1.5 dB. If the
deviations in both the adjacent antenna patterns are opposite, the maximum error
that can occur is 12.77%, which corresponds to the row with serial number 4 of
Table 20.3. The maximum error corresponding to that is 11.5° from the actual radio
direction if the model pattern is used. Table 20.4 lists the actual errors measured
318 ASTRONOMY AND ASTROPHYSICS
for some of the FM radio and TV spectrums. Note that the percentage errors are
calculated with αv = 90°. Since the system is omni-directional, the effective viewing
angle is 360°, in which case the tabulated errors should be divided by 4.
Figure 20.10: Resolved terrestrial radio sources. Top: Angular positions of the sources in a Horizon
coordinate system. Bottom: Flux density at the vicinity of the antennas.
x11 x1n
X
(20.26)
xm1 xmn
where,
m no. of time stamped traces,
n no. of data points/trace
30
It is an optional choice that may help in improving the signal-to-noise ratio. The maximum
number of traces that can be integrated/averaged is decided based on the number of traces
occupied by the spectrum feature, which is under study.
320 ASTRONOMY AND ASTROPHYSICS
Δν span
χ = , χ 10 (20.28)
Δν MaxNarrowBand
For further processing, the data in Eq. (20.26) is now reorganized in groups of
sub-matrices of Xij as in Eq. (20.30). There are k numbers of sub-matrices in each
row. Each of the Xij sub-matrix is a row matrix containing (n/k) columns. The ratio
(n/k) and k must be integers. Strictly speaking, the value of k should be chosen
based on the definition set for the minimum bandwidth of the broadband spectrum
under study in terms of the number of data points it occupies in any single time
stamped trace or in any row of the matrix X. But for all practical purposes, the
value of k can be chosen greater than or equal to 10 times of nMaxNarrowBand, which is
the maximum number of data points occupied by a narrowband spectrum.
X11 X1k
X ,
X m1 X mk
where
n n n (20.30)
Xij xi ( j 1) 1 xi ( j 1) 2 xi j
k k k
n
is an integer, k 10 nMaxNarrowBand
k
Eq. (20.32) and is saved as matrix X1. This process of chopping may be repeated
several times by substituting the original data matrix X of Eq. (20.27) with X1 of
Eq. (20.32) and process the Eqs. (20.27) through (20.32) with smaller values of p
until the narrowband spectrum features diminish toward acceptable limits.
x111 x11n
X1
x1m1 x1mn
(20.32)
where, x1ij mij if xij ρ mij
x1ij xij otherwise
where, ρ 1
If data in X1 is plotted, it may be observed that the slope of each of the power
versus frequency curve in each row changes abruptly on either side of the locations
where the major power of the narrowband spectrum features existed earlier. In
other words, traces of the narrowband spectrum features still exist however less in
quantity. In order to further smoothen, a discrete/fast Fourier transform operation
is performed first. Next, the high frequency components in the Fourier transform
output are attenuated substantially and then a discrete inverse Fourier transform
is taken. The magnitude of the inverse Fourier transform is smooth and contains
the broadband spectrum features.
However, if the transformations are directly applied to the data, undesirable
effects near the edges (channel 1 and channel n) due to the edge discontinuities
takes place. In order to prevent them, two matrices, namely P and Q, containing
the extreme end data values are padded on either side of the matrix X1 and saved
as matrix Y1 as shown in Eq. (20.33). The Fourier transform is applied to every
individual time stamped row of matrix Y1, and the results are stored in matrix F
as shown in Eq. (20.34). Thereafter, attenuation is applied to the high-frequency
components by suitably selecting the values of αij and storing products in matrix
F1 as shown in Eq. (20.35). Following this, an inverse Fourier transform is applied
322 ASTRONOMY AND ASTROPHYSICS
on each of the time stamped data rows of F1, and the results are saved in matrix
Y2 as expressed in Eq. (20.36). The values in the extra padded sections in Y2 are
removed and the magnitudes of the remaining complex data are saved in X2 as
expressed in Eq. (20.37). The entire smoothening operation involving Eqs. (20.33)
through (20.37) are called second-level broadband processing. It is observed that
the operations involving Eqs. (20.35) through (20.37) may produce small dips at
certain regions. These may be safely corrected by comparing with the raw data
thereby saving the final values in matrix X3 as expressed in Eq. (20.38). This
operation is called third-level broadband processing.
x111 x111
Y1 [P X1 Q], where P ,
x1m1 x1m1
(20.33)
11
1n 2 k
x11n x11n f f
Q , F
x1mn x1mn f
m1 f mn 2 k
y 211 y 21n 2 k
Y2
y2
m1 y 2mn 2 k
(20.36)
n 2k 2π 10.5 l 1j 1
y 2il j 1
dij exp
n 2k
,
1 i m, l 1, 2, ...(n 2k )
PROBLEMS OF RADIO INTERFERENCE AND NEW DEVELOPMENTS IN RADIO ASTRONOMY 323
x 211 x 21n
X2 (20.37)
x 2m1 x 2mn
x 2ij y 2i( j k ) where 1 i m, 1 j n
x 311 x 31n
X3 ,
x 3m1 x 3mn
(20.38)
x 3ij xij if x 2ij κ(smallest xij ) xij ,
0 κ 1,
x 3ij x 2ij otherwise
Figure 20.11: A three-dimensional plot and contour plot of the extracted broadband spectrum features in
time-frequency-power space for the raw data displayed in Fig. 20.2.
324 ASTRONOMY AND ASTROPHYSICS
Figure 20.12: A three-dimensional plot and contour plot of the extracted narrowband spectrum features in
time-frequency-power space for the data shown in Fig. 20.2.
31
It is assumed that the noise-floor of the instrument is a constant across the entire observed band.
An ideal instrument is noise-free.
PROBLEMS OF RADIO INTERFERENCE AND NEW DEVELOPMENTS IN RADIO ASTRONOMY 325
x 411 x 41n
X4
x 4 m1 x 4 mn
(20.40)
x 4ij xij if xij γ(smallest x 3ij ) x 3ij ,
0 γ1
otherwise x 4ij Noise-floor,
1 i m, 1 j n
S X 3 or X 4
MaskLevel δ(smallest sij ) , 1 δ 2
(20.41)
e11 e1n
E
em1 emn (20.42)
if sij MaskLevel, eij 1
else eij 0, 1 i m, 1 j n
326 ASTRONOMY AND ASTROPHYSICS
Figure 20.13: Percentage spectral occupancy over entire observed time for the data shown in Fig. 20.2.
who wishes to write his/her own computer code. It is felt that these algorithms
can also be extended for other EMI/EMC analysis apart from radio astronomy.
Figure 20.14: Categorically analyzed of terrestrial spectrum data collected at the GMRT site.
Table 20.5: Parameter settings for the various results of data analysis presented.
Parameter Fig. 20.11 Fig.20.14 Fig. 20.14 (c) Fig. 20.14 (d)
Fig. 20.12 (a) & (b) 450–850 MHz 850–1250 MHz
120–160 MHz 50–450 MHz
χ 10 10 10 10
k 50 50 50 50
ς 20 20 20 20
In general, the allocated bandwidth for radio astronomy is only a few MHz centered
at certain frequencies important from radio astronomy’s point of view. There is a
huge ongoing battle between the communication industries and radio astronomers
for capturing unallocated spectrum bands. There is very little bandwidth left for
allocation at higher frequencies. As the atmosphere is not transparent to the whole
electromagnetic spectrum, the spectrum use for radio astronomy is very limited.
Out of this, the majority of the bandwidth is used for commercial purposes due
to which radio astronomers are taking refugees in deserted places where the
communication signals do not reach the radio telescopes. There is also some
thinking among astronomers to launch satellites and also possibly use the Moon
for setting up radio telescope antennas. But these are an extremely adventurous
and costly affair. Radio astronomers have formed international committees to
present their claims for spectrum usage globally.
Certain frequency bands are not always occupied by communication users.
Radio astronomers feel that these frequency bands could be used efficiently for
exploring the space when they are free. The major requirement of such technologies
is an ultra wideband feed fitted on a parabolic dish, which can cover decades of
bandwidths. The free spectrum can therefore be extracted by setting a selectable
bandwidth band-pass filter. This futuristic view of acquiring increased bandwidth
of the front-end system of radio telescopes has led many engineers to develop
highly sensitive antenna feeds covering decade bandwidths. Some are reported to
be partially successful and have shown the possibility of attaining perfection. One
such feed has also been reported from the GMRT. This feed has already attained
the perfection for use in communication systems.
PROBLEMS OF RADIO INTERFERENCE AND NEW DEVELOPMENTS IN RADIO ASTRONOMY 329
The detail of the antenna feed is shown in Fig. 20.15. It consists of two log-periodic
planar toothed trapezoidal structures (for orthogonal polarizations) with a Step-
Lane reflector (thick lines) or a Pyramidal reflector (dotted lines).32 An optional
secondary reflector can also be seen, which boosts the gain at the lowermost
frequencies. The antenna has been fabricated first on printed circuit board and
later again on a metallic sheet cut, shaped, and fixed on a low permittivity thin
dielectric surface (for mechanical support). In order to further reduce the effect
of the dielectric, slots were made by removing the dielectric between the dipoles.
This has also reduced the resistance to air flow. For any particular design, only one
of the frequency independent reflectors, Step-Lane or Pyramidal, is positioned
below the antennas. The step widths of the Step-Lane reflector are kept identical
Figure 20.15: Elevation and plan of the dual polarized planar log-periodic antenna structures associated
with a Step-Lane reflector (dark lines) or a Pyramidal reflector (dotted lines). An optional
secondary reflector is provided for improving the gain at the lowermost frequencies.
32
The names Step-Lane and Pyramidal has been given by the inventors of the frequency independent
reflectors for use with frequency independent antennas.
330 ASTRONOMY AND ASTROPHYSICS
to the dipole widths of the antennas, such that the former hide themselves behind
the antennas when viewed from the top. For any particular frequency, the Z-
axis distance between a dipole and its corresponding reflecting step is about a
quarter of an wavelength. Similarly, for the case of the Pyramidal reflector, the
Z-axis distance between any dipole and the reflector surface is maintained as a
quarter of a wavelength. Thus, for any radiation falling directly perpendicular to
the antenna plane (along Z-axis) within the designed frequency range, there exists
a dipole and a reflecting surface to enhance the gain at that frequency. Essentially,
within the designed spectrum band, the flux escaping the respective frequency
dipoles get reflected back on them with an identical phase for every frequency,
and hence named frequency independent reflector. Effectively, the reflector converts
the bidirectional radiation pattern into a unidirectional pattern with improved
average antenna gain approximately by a factor of 2.0.
The gain simulations were done using an antenna and electromagnetic
modeling software script on MATLAB. Two computational models were prepared;
one with the Step-Lane reflector and other with the Pyramidal reflector. The
secondary reflector was included in both the models. Each of these models was
broken into more than 4000 finite element triangles. The modeling and meshing
were done using the GNU GMSH software. Fig. 20.16 shows the simulated and
measured33 gains over the 1:10 frequency bandwidth. The quantity F stands for the
minimum cut-off frequency of the ultra-wide band. The average gain over the ultra-
wide band is 8.9 and 9.1, respectively, for a Step-Lane reflector and a Pyramidal
Figure 20.16: Gain variation of the antenna-feed over the ultra-wide band. Please note that the measured
gain is for an antenna-feed designed for a frequency range of 200 MHz to 2000 MHz for use
on the GMRT dishes.
33
The measured gains stand for a design specially done for GMRT using a Step-Lane reflector. The
frequency range of this antenna-feed is from 200 MHz to 2000 MHz. The Pyramidal reflector
version of the antenna-feed is under the process of fabrication and yet to be tested.
PROBLEMS OF RADIO INTERFERENCE AND NEW DEVELOPMENTS IN RADIO ASTRONOMY 331
The sensitivity of a radio telescope is dependent on the power it can collect from the
astronomical sources. Thus, the more light (radio wave) collecting area a telescope
can have, the more will be its sensitivity. The angular resolution of any telescope
is measured roughly as λ/D, where λ is the wavelength and D is the diameter
of the reflector. The interferometers can further increase the angular resolution,
which is a function of the baselines. Unlike optical telescopes that operate in
wavelengths of angstrom units, radio telescopes operate in wavelengths from
a few mm to a few meters (the lower part of the spectrum of radio telescopes
are of great importance). The major problem of building large radio telescopes is
the huge size of the paraboloid reflectors. The sensitivity of a telescope depends
on the total light collecting area. In radio wavelengths, the signals are extremely
weak and require a huge collecting area. Thus, radio telescopes are fabricated in
pieces—many antennas are placed in an array. The array powers may be combined
in interferometer mode or the phased array mode depending on applications. For
example, for observing distant sources like a radio galaxy which is an extended
radio source, the interferometer mode would be useful since one requires high
angular resolution for producing an intensity map of the extended radio source.
On the other hand, for observing a pulsar, one may use the phased array mode
since its electromagnetic properties are of interest to the observer and not the
shape of the source. But in either case, the total collecting area plays a major role in
improving sensitivity.
The Square Kilometer Array, or SKA, is a proposed integrated radio telescope,
which intends to produce a light collecting area in squares of kilometers over the
frequency range roughly covering 100 MHz to 25 GHz. The frequency range may
further be extended below and above the said limits. The super collecting area
will make it 50 times more sensitive than any other radio telescope array that
exists today. Various engineering techniques such as electronic beam switching
and antenna band sharing are proposed to be implemented so that several radio
astronomers may use the instrument simultaneously for looking into different
parts of the sky and at different frequencies of observations. Interferometer
332 ASTRONOMY AND ASTROPHYSICS
techniques will be used to make images of much distant sources, which are not
possible today with present instruments. An international consortium consisting
of representatives from various countries is working on the design aspects of the
SKA, which is to be finalized in 2008.
The array is expected to comprise many antennas spread over an area of several
thousands of kilometers. Half of its collecting area will be from a compact core of
antennas placed within 150 km. The remaining antennas will be spread up to a few
thousand kilometers to increase its angular resolution by increasing the baseline.
A choice has to be made between the two sites short listed for building the array,
Mileura station in Western Australia and Karoo Basin in the Northern Cape area
of South Africa. After the site is selected, construction of the SKA will begin in
2010, with observations to start from 2015 onwards. The complete functionality
is expected to be available from 2020 onwards. The instrument will be sensitive
to detect every AGN out to a red-shift of 6, when the Universe was less than one
billion years old.
EXERCISES
BACKGROUND INFORMATION,
FORMULAE, RELATIONSHIPS
I
n the Sun the gases that are at high temperatures produce electrically charged
particles. Due to their collisions, electromagnetic waves of different wavelengths
are emitted by the Sun, which extend from gamma rays (λ = 10–13 m) to radio
3
waves (λ = 10 m). However, about 99% of solar radiation lies between 0.15 to
4.0 μm. Out of this: (i) about 46% lies in IR (0.7 μm to 2.3 μm), which causes radiant
heat with some chemical effect; (ii) about 44% lies in the visible light (0.38 μm to
0.7 μm), violet to red; (iii) about 9% lies in the UV radiation (0.15 μm to 0.38 μm),
which causes photochemical effects, sunburn, bleaching, etc. All these radiations
8
travel with 3 × 10 m/s in an average surface temperature of 6,000° K, while the
Earth emits at its average surface temperature of 15° C or 288° K.
From Wein’s law we can write
constant
λmax = ,
T
where the constant = 2897 μmK ≈ 3000 μmK.
333
334 ASTRONOMY AND ASTROPHYSICS
3000 μ mK
λmax = 0.5 μm.
6000 K
This shows that the Sun emits a maximum amount of radiation energy
corresponding to the wavelength of 0.5 μm.
On the other hand, we can write for the Earth:
3000 μ mK
λmax = 10 μm [considering T = 300 K].
300 K
This shows that the Earth emits a maximum amount of radiant energy
corresponding to the wavelength of 10 μm.
From the above calculation, we understand why solar radiation is called
shortwave radiation and terrestrial radiation is called long-wave radiation.
The mean Earth-Sun distance is a⊕= 1 AU ≈ 1.5 × 1011 m. Also, the angular
diameter of the Sun as observed from the Earth = 9.3 × 10–3 radian, which is slightly
more than half a degree. Hence,
arc
= radian; i.e., arc = radius × radian.
radius
Now, the Sun’s diameter = Earth-Sun distance × angular diameter of the Sun
11 –3 8
= 1.5 × 10 × 9.3 × 10 m = 13.95 × 10 m.
8
The radius of the Sun ≈ 6.97 × 10 m.
We have
the mass of the Sun, MΘ = 1.99 × 1030 kg, and
the radius of the Sun, RΘ = 6.97 × 108 m.
Now using the formula,
G MΘ
gΘ = [where, G = 6.67 × 10–1 Nm2/kg],
RΘ2
MΘ 3 MΘ 3 1.99 10 30
ρΘ = 1.41 10 3 kg/m 3
8 3
VΘ 4 πRΘ 3
4 3.14[6.96 10 ]
The amount of solar radiation (i.e., energy incident on a unit area in unit time
on a surface held at right angles to the solar beam at the outer boundary of the
atmosphere) is called the solar constant. If it is denoted by SΘ then,
LΘ 3.83 10 26 7 2
S = 6.29 10 w/m
8 2
4 π RΘ 2
4 3.14[6.96 10 ]
It has been noted that mass of the Sun consists of about 80% hydrogen,
19% helium, and 1% other elements. The Sun is a burning gaseous globe with a
15 30
radius of about 7 × 10 km and a mass of about 2.0 × 10 kg.
It has been estimated the solar energy received per year is 1.5 × 1015 Mwh. On a
clear sunny day, 1 acre of land in the daytime receives solar energy equivalent to
16 bbl (bbl = billion barrels) of crude oil. The total solar energy that can be
intercepted by the Earth = 2.55 × 1018 cal min–1 = 3.67 × 1021 cal day–1 .
The deposition of solar radiation in the Earth atmosphere is given by
Qs = Cr + Ar + Ca + Aa + (Q + q) (1 – α) + (Q + q)α.
Here, Qs = amount of solar radiation incident at the top of the atmosphere
= 263 kly/yr.
Q = Direct beam of solar radiation incident on horizontal surface at
the ground = 82 kly/yr.
q = Diffused solar radiation incident on horizontal surface at the
ground = 58 kly/yr.
Cr = Solar radiation reflected and scattered back to space by clouds.
Ar = Solar radiation reflected and scattered back to space by aerosol
(i.e., air molecules, dust, and water vapor of the atmosphere).
Ca = Solar radiation absorbed by the clouds.
Aa = Solar radiation absorbed by atmospheric aerosols.
α = Albedo of the surface ≈ 0.1143.
Solar energy per unit area is usually expressed in langely (ly) or kilolangely
(kly), where
2
1 ly = 1 cal/cm
1kly = 1000 ly
APPENDIX A 337
26
It is assumed that in each minute the Sun radiates about 56 × 10 cal of energy
in all. Now the solar constant
56 10 26 56 10 26
SΘ = cal/min
4 π( a 2 ) 4 3.14 [1.5 1013 ]2
2
= 2.0 ly/min = 1359 w/m .
2
The total solar energy intercepted by the Earth in unit time = π R⊕ SΘ
The amount of solar energy received at the top of the atmosphere in unit area
in unit time is given by,
πR 2SΘ SΘ
Qs = 0.5 ly/min
4 π R 2
4
6 –2
Internal energy = 1674.8 × 10 J m .
6 –2
Total PE = 2242.3 × 10 J m .
3 –2
KE = 1153.4 × 10 J m .
1153.4 10 3 1
Therefore, an average amount of KE = of total PE.
6
2242.3 10 2000
It is seen that the stars take about 3 minutes and 56 seconds less than the Sun to go
once around the Earth. So we can write,
1 civil day = 24 hours.
h m s
1 sidereal day = 23 56 4 .
The time taken by the Sun to complete a revolution around the Earth, as observed
from the Earth, is defined as a solar year. Depending on the points chosen we have
different types of solar years, out of which we consider here the sidereal solar year
and the tropical solar year.
The duration of the sidereal year = 365.256364 days.
The duration of a tropical year = 365.242190 days.
It is to be noted that
T 365.25 365.25 24
1 sidereal day = civil day = hours
T 1 366.25 366.25
The angular extent of the ecliptic is 360° and the Sun takes about 365.25 civil days
to complete one revolution. A solar month is equal to about 365.25/12 days. To be
precise, we have
APPENDIX A 339
If H is the hour angle, then cos H = (sin a – sin φ sin δ)/cos φ cos δ, where a is the
altitude, δ is the declination, and φ is the latitude at the time of sunrise and sunset,
and the Sun’s altitude a = 0.
Therefore, cos H = –tan φ tan δ.
This is the required relations to be used. We give an example below.
Find the times of sunrise and sunset and also the duration of day and night at
Bangalore (φ = 13° N) on May 15, given that the Sun’s declination is 18°47’N.
Solution:
STANDARD SYMBOLS
n Symbol of a neutron
+
p Symbol of a proton
–
e Symbol for an electron
+
e Symbol for a positron
λ Symbol of wavelength
μ Symbol for refractive index
γ Symbol of gamma ray
ν Symbol of frequency
12
6C Symbol of nuclei for carbon isotope. The superscript represents the
atomic mass (in amu) and the subscript is the atomic number.
G Gravitational constant = 6.670 × 1011 N m2 kg–2
h Planck’s constant = 6.6262 × 10–34 J s
k Boltzmann’s constant = 1.3806 × 10–23 J K–1
me Mass of an electron = 9.1096 × 10–31 kg
amu Atomic mass unit = 1.6605 × 10–27 kg
c Velocity of light in vacuum = 2.9979 × 108 m s–1
e Charge of an electron 1.6022 × 10–19 C
eV Electron Volt = 1.6022 × 10–19 J
U,B,V Magnitudes through the standard U B V photometric system
ε0 Permittivity of free space = 8.854 × 10–12 F m–1
μ0 Permeability of vacuum = 4π × 10–7 H m–1
341
342 ASTRONOMY AND ASTROPHYSICS
APPENDIX C
Included on the CD-ROM are figures from the text(most in color) and other
files related to topics for astronomy and astrophysics
See the “README" files for any specific information/system requirements
related to each file folder, but most files will run on Windows XP or higher
343
344 ASTRONOMY AND ASTROPHYSICS
REFERENCES
345
346 ASTRONOMY AND ASTROPHYSICS
di Gesu V., et al. (1992). Data Analysis in Astronomy IV (New York: Plenum).
Edward W. K. and Michael S. T. (1990). The Early Universe (Westview Press).
Emerson D. T. and Payne J. M. (ed) (1995). Multi-Feed Systems for Radio Telescopes
(Astronomical Society of the Pacific Conference Series vol. 75).
Evans J. V. and Hagfors T. (1968). Radar Astronomy (Dordrecht: Reidee).
Goddard D. E. and Milne D. K. (ed) (1994). Parkes: Thirty Years of Radio Astronomy
(CSIRO).
Golay M. (1974). Introduction to Astronomical Photometry (Dordrecht: Reidel).
Guyenne T. (ed) (1995). Future Possibilities for Astronomy in Space (European Space
Agency).
Harwit M. (1973). Astrophysical Concepts (New York: Wiley)
Hearnshaw J. B. (1996). Measuring of Starlight (Cambridge: Cambridge University Press).
Hiltner W. A. (ed) (1962). Astronomical Techniques (Stars and Stellar Systems vol. II)
(Chicago: University of Chicago Press).
Hopkins J. (1976). Glossary of Astronomy and Astrophysics (Chicago: University of
Chicago Press).
James B. and Scott T. (1987). Galactic Dynamics (Princeton University Press).
Jeschek C. and Murtagh F. (ed) (1990). Errors, Bias and Uncertainties in Astronomy
(Cambridge: Cambridge University Press).
Joardar S. and Bhattacharya A. B. (2006). Journal of Electromagnetic Waves and Applications,
vol. 20, p. 1195 & p. 1465.
Joardar S. and Bhattacharya A. B. (2006). Journal of Electromagnetic Waves and Applications,
vol. 20, No. 11, 1465–1479, 2006.
Joardar S. and Bhattacharya A. B. (2007). Journal of Electromagnetic Waves and Applications,
vol. 21, p. 441.
Joardar S., Datta A., Bhattacharya A. B. and Bose S. K. (2007). IETE Journal of Research,
vol. 53, no. 1, p. 83.
Joardar S. and Bhattacharya A. B. (2007). Progress In Electromagnetics Research, PIER 72,
p.127.
Kitchin C. R. (1995). Telescopes and Techniques (Berlin: Springer).
Kitchin C. R. and Ferrest R. F. (1997). Seeing Stars (Berlin: Springer).
Kitchin C. R. (2005). Astrophysical Techniques (New Delhi: Overseas Press India
Pvt. Ltd.)
Kraus J. D. (1966). Radio Astronomy (New York: McGraw Hill).
Lang K. R. (1978). Astrophysical Formula (Berlin: Springer).
Mackintosh A. (1986). Advanced Telescope Making Techniques, vol. I & II (William Bell).
REFERENCES 347
349
350 ASTRONOMY AND ASTROPHYSICS
V
X
Varahamihira, 2
Vedas, 1 X-ray telescopes, 80
Vega, 55
Venus, 15-21, 24, 26-30
Y
Very Large Array (VLA), 86
Very Long Baseline Interferometry (VLBI), Yajnavalkya, 1
9
Yojan, 2
Viaduct, 119
Viking orbiter, 18
Virgo, 295 Z
Virgo-A, 148
Zeeman effect, 73
Vogt, Heinrich, 242
Zeta Puppis, 54
Volcanism, 35
Zones, 15-21
Voltage polar diagrams, 167
Zurich ETH spectrograph, 236-237
Voyager-1, 17
Zwicky, Fritz, 271