IOAA 2016 Problems
IOAA 2016 Problems
IOAA 2016 Problems
PROBLEMS
Table of Constants
Page 1 of 1
Fundamental Constants
Speed of light in vacuum = 2.998 108 m s 1
Planck Constant = 6.626 1034 J s
Boltzmann Constant B = 1.381 1023 J K 1
Stefan-Boltzmann Constant = 5.670 108 W m2 K 4
Charge of electron = 1.602 1019 C
Universal Gravitational Constant = 6.674 1011 N m2 kg 2
Universal Gas Constant = 8.315 J mol1 K 1
Avogadro Constant A = 6.022 1023 mol1
Wien's displacement law m = 2.898 103 m K
Mass of electron e = 9.109 1031 kg
Mass of proton p = 1.673 1027 kg
Mass of neutron n = 1.675 1027 kg
Atomic Mass Unit (a.m.u.) = 1.661 1027 kg
Astronomical Data
1 parsec (pc) = 3.086 1016 m
1 astronomical unit (AU) = 1.496 1011 m
Solar Mass = 1.989 1030 kg
Solar Radius = 6.955 108 m
Solar Luminosity = 3.826 1026 W
Apparent magnitude of the Sun at mid-day = 26.72 mag
Solar Constant (at Earth) = 1366 W m2
Apparent angular diameter of Sun = 30
Earth Mass = 5.972 1024 kg
Earth Radius = 6.371 106 m
1 tropical year = 365.242 solar days
= 3.156 107 s
Theoretical Examination
Page 1 of 1
(T4) Shadows 10
An observer in the northern hemisphere noticed that the length of the shortest shadow of a 1.000 m
vertical stick on a day was 1.732 m. On the same day, the length of the longest shadow of the same
vertical stick was measured to be 5.671 m.
Find the latitude, , of the observer and declination of the Sun, , on that day. Assume the Sun to be a
point source and ignore atmospheric refraction.
A massive body thus behaves somewhat like a focusing lens. The light rays coming from infinite distance
beyond a massive body, and having the same impact parameter , converge at a point along the axis, at a
distance from the centre of the massive body. An observer at that point will benefit from huge
amplification due to this gravitational focusing. The massive body in this case is being used as a
Gravitational Lensing Telescope for amplification of distant signals.
(T10.1) Consider the possibility of our Sun as a gravitational lensing telescope. Calculate the shortest 6
distance, min , from the centre of the Sun (in A. U.) at which the light rays can get focused.
(T10.2) Consider a small circular detector of radius , kept at a distance min centered on the -axis and 8
perpendicular to it. Note that only the light rays which pass within a certain annulus (ring) of
width (where ) around the Sun would encounter the detector. The amplification factor
at the detector is defined as the ratio of the intensity of the light incident on the detector in the
presence of the Sun and the intensity in the absence of the Sun.
Express the amplification factor, m , at the detector in terms of and .
(T10.3) Consider a spherical mass distribution, such as dark matter in a galaxy cluster, through which 6
light rays can pass while undergoing gravitational bending. Assume for simplicity that for the
gravitational bending with impact parameter, , only the mass () enclosed inside the radius
is relevant.
What should be the mass distribution, (), such that the gravitational lens behaves like an ideal optical
convex lens?
Theoretical Examination
Page 4 of 6
The test mass loses energy due to the emission of gravitational waves. As a result the orbit keeps on
shrinking, until the test mass reaches the surface of the object, or in the case of a black hole, the innermost
stable circular orbit ISCO which is given by ISCO = 3sch , where sch is the Schwarzschild
radius of the black hole. This is the epoch of merger". At this point, the amplitude of the gravitational
wave is maximum, and so is its frequency, which is always twice the orbital frequency. In this problem,
we will only focus on the gravitational waves before the merger, when Keplers laws are assumed to be
valid. After the merger, the form of gravitational waves will drastically change.
(T11.1) Consider the observed gravitational waves shown in the figure above. Estimate the time period, 3
0 , and hence calculate the frequency, 0 , of gravitational waves just before the epoch of
merger.
(T11.2) For any main sequence (MS) star, the radius of the star, MS , and its mass, MS , are related by 10
a power law given as,
MS (MS )
where = 0.8 for < MS
= 1.0 for 0.08 MS
If the central object were a main sequence star, write an expression for the maximum frequency
of gravitational waves, MS , in terms of mass of the star in units of solar masses (MS / )
and .
(T11.3) Using the above result, determine the appropriate value of that will give the maximum 9
possible frequency of gravitational waves, MS,max for any main sequence star. Evaluate this
frequency.
(T11.4) White dwarf (WD) stars have a maximum mass of 1.44 (known as the Chandrasekhar limit) 8
and obey the mass-radius relation 1/3 . The radius of a solar mass white dwarf is equal
Theoretical Examination
Page 5 of 6
to 6000 km. Find the highest frequency of emitted gravitational waves, WD,max , if the test
mass is orbiting a white dwarf.
(T11.5) Neutron stars (NS) are a peculiar type of compact objects which have masses between 8
1 and 3 and radii in the range 10 15 km. Find the range of frequencies of emitted
gravitational waves, NS,min and NS,max , if the test mass is orbiting a neutron star at a distance
close to the neutron star radius.
(T11.6) If the test mass is orbiting a black hole (BH), write the expression for the frequency of emitted 7
gravitational waves, BH , in terms of mass of the black hole, BH , and the solar mass .
(T11.7) Based only on the time period (or frequency) of gravitational waves before the epoch of merger, 5
determine whether the central object can be a main sequence star (MS), a white dwarf (WD), a
neutron star (NS), or a black hole (BH). Tick the correct option in the Summary Answersheet.
Estimate the mass of this object, obj , in units of .
(T12) AstroSat
India astronomy satellite, AstroSat, launched in September 2015, has five different instruments.
In this question, we will discuss three of these instruments (SXT, LAXPC, CZTI), which point in the
same direction and observe in X-ray wavelengths. The details of these instruments are given in the table
below.
You should note that LAXPC energy range is divided into 8 different energy band counters of equal
bandwidth with no overlap.
(T12.1) Some X-ray sources like Cas A have a prominent emission line at 0.01825 nm corresponding 13
to a radioactive transition of 44Ti. Suppose there exists a source which emits only one bright
emission line corresponding to this transition. What should be the minimum relative velocity
() of the source, which will make the observed peak of this line to get registered in a different
energy band counter of LAXPC as compared to a source at rest?
These instruments were used to observe an X-ray source (assumed to be a point source), whose energy
spectrum followed the power law,
( ) = 23 [in units of counts/keV/m2 /s ]
where is the energy in keV, is a constant and () is photon flux density at that energy. Photon flux
density, by definition, is given for per unit collecting area (m2 ) per unit bandwidth (keV) and per unit
Theoretical Examination
Page 6 of 6
time (seconds). From prior observations, we know that the source has a flux density of 10
counts/keV/m2 /s at 1 keV, when measured by a detector with 100% photon detection efficiency. The
counts here mean the number of photons reported by the detector.
As the source flux follows the power law given above, we know that for a given energy range from
1 (lower energy) to 2 (higher energy) the total photon flux (T ) will be given by
1 3 1 3
T = 3 (2 1 ) [in units of counts/m2 /s ]
(T12.2) Estimate the incident flux density from the source at 1 keV, 5 keV, 40 keV and 100 keV. Also 8
estimate what will be the total count per unit bandwidth recorded by each of the instruments at
these energies for an exposure time of 200 seconds.
(T12.3) For this source, calculate the maximum exposure time (S ), without suffering from saturation, 4
for the CCD of SXT.
(T12.4) If the source became 3500 times brighter, calculate the expected counts per second in LAXPC 8
counter 1, counter 8 as well as total counts across the entire energy range. If we observe for
longer period, will the counter saturate due to any individual counter or due to the total count?
Tick the appropriate box in the Summary Answersheet.
(T12.5) Assume that the counts reported by CZTI due to random fluctuations in electronics are about 10
0.00014 counts per pixel per keV per second at all energy levels. Any source is considered as
detected when the SNR (signal to noise ratio) is at least 3. What is minimum exposure time,
, needed for the source above to be detected in CZTI?
Note that the noise in a detector is equal to the square root of the counts due to random
fluctuations.
(T12.6) Let us consider the situation where the source shows variability in number flux, so that the 7
factor increases by 20%. AstroSat observed this source for 1 second before the change and
1 second after this change in brightness. Calculate the counts measured by SXT, LAXPC and
CZTI in both the observations. Which instrument is best suited to detect this change? Tick the
appropriate box in the Summary Answersheet.
Theoretical Examination
Page 1 of 27
Solution:
T 2.0
The colour of the clear sky during night is the same as during daytime, since the
spectrum of sunlight reflected by the Moon is almost the same as the spectrum of
sunlight. Only the intensity is lower.
(T1.2) An astronomer at Bhubaneswar marks the position of the Sun on the sky at 05:00 UT 2
every day of the year. If the Earths axis were perpendicular to its orbital plane, these
positions would trace an arc of a great circle.
Solution:
T 2.0
If the Earths axis were perpendicular to its orbital plane, the celestial equator will
coincide with ecliptic and the Sun will remain along the celestial equator every day.
However, as the Earths orbit is elliptical, the true sun would still lead or lag mean
sun by a few minutes on different days of year.
(T1.3) If the orbital period of a certain minor body around the Sun in the ecliptic plane is less 2
than the orbital period of Uranus, then its orbit must necessarily be fully inside the orbit
of Uranus.
Solution:
F 2.0
The semi-major axis of the orbit of the body will be less than that of Uranus. However
the minor bodys orbit may have a high eccentricity, in which case it may go outside
that of Uranus.
(T1.4) The centre of mass of the solar system is inside the Sun at all times. 2
Solution:
F 2.0
The centre of mass of Sun-Jupiter pair is just outside the Sun. Thus, if all gas giants
are on same side of the Sun, the centre of mass of Solar system is definitely outside
the Sun.
(T1.5) A photon is moving in free space. As the Universe expands, its momentum decreases. 2
Solution:
T 2.0
For photons the wavelength increases when the Universe expands.
Theoretical Examination
Page 2 of 27
Solution:
As the gas is monatomic,
3 1 2
kB TT mg vrms
2 2
Mg 2
3kB TT v
N rms
sA
3kB NA TT
vrms 4.0
Mg
50% deduction if 3/2 pre-factor is not used and 1/2 or 1 are used instead.
Full credit if students writes the relation for vrms directly.
To remain in atmosphere,
r
vesc 1 2GMT
vrms < =
6 6 RT
s r
3kB NA TT GMT
<
Mg 18RT
54kB NA TT RT
Mg > 4.0
GMT
54 1.381 1023 6.022 1023 93.7 2.575 106
> g
6.6741 1011 1.23 1023
> 13.2 g 1.5
Thus, all gases with atomic weight more than Amin = 13.2 will be retained in the atmosphere
of Titan. 0.5
Half mark for understanding that atomic mass has no units.
Alternative solution
3 1 2
kB TT mg vrms
2 2
s
3kB TT
vrms 4.0
mg
54kB TT RT
mg > 4.0
GMT
mg > 2.19 1026 kg 1.0
mg 2.19 1026 kg
Amin = =
atomic mass unit 1.66 1027 kg
Theoretical Examination
Page 3 of 27
Solution:
m0 /c m0 0.3
= = 4 = 3000
r0 /c r0 10
At ze , both matter density and radiation density were equal.
r = m
r0 (1 + ze )4 = m0 (1 + ze )3
m0
1 + ze = = 3000 2.0
r0
ze ' 3000 1.0
Only ze = 2999 and ze = 3000 are acceptable answers.
(T3.2) Assuming that the radiation from the early Universe has a blackbody spectrum with a 4
temperature of 2.732 K, estimate the temperature, Te , of the radiation at redshift ze .
Solution:
As the Universe behaves like an ideal black body, the radiation density will be pro-
portional to the fourth power of the temperature (Stefans law).
4
Te r
= e
T0 r0
r (1 + ze )4
= 0 2.0
r0
4
Te
= (1 + ze )4 1.0
2.732
Te
= 1 + ze = 3000
2.732
Te = 3000 2.732
Te = 8200 K 1.0
8100 Te 8200 gives 1.0; 8200 < Te 9000 gives 0.5; else 0.
(T3.3) Estimate the typical photon energy, E (in eV), of the radiation as emitted at redshift 3
ze .
Theoretical Examination
Page 4 of 27
Solution:
Wiens law:
0.002898 m K
max =
Te
0.002898
= m = 354 nm 1.0
8200
hc
E =
max
6.62 1034 3 108 5.62 1019
= J = eV 1.0
354 109 1.602 1019
E = 3.5 eV 1.0
Alternative solution:
E = kB Te
1.13 1019
= 1.38 1023 8200 J = eV 2.0
1.602 1019
E = 0.71 eV 1.0
Use of either Wiens law or E = kB T gets full credit. E = 3kB T /2 or
similar gets no credit. Answers with E = 3kB T or E = 2.7kB T also get
full credit.
(T4) Shadows 10
An observer in the northern hemisphere noticed that the length of the shortest shadow of a 1.000 m
vertical stick on a day was 1.732 m. On the same day, the length of the longest shadow of the
same vertical stick was measured to be 5.671 m.
Find the latitude, , of the observer and declination of the Sun, , on that day. Assume the Sun
to be a point source and ignore atmospheric refraction.
Solution:
As the longest shadow of the Sun on the given day is of finite length, the Sun is circumpolar
for this observer on this day. 2.0
90
90
A
1
2
O S
In the figure above, the left panel shows the shadow OS formed by stick OA (of length 1.000 m),
and the right panel shows the Suns location in two cases.
Theoretical Examination
Page 5 of 27
2 = (90 ) = 90 + 1.0
cot( 90 + ) = 5.671
1
tan( 90 + ) =
5.671
1 1
+ = tan + 90 = 100 = 1.745 rad 1.5
5.671
Solving,
= 80 = 1.396 rad 1.0
= 20 = 0.349 rad 1.0
Given high accuracy of shadow length, only 0.5 is allowed.
One can also solve the question by manipulating tan( ) and tan( + ), to get tan()
and tan( ).
Solution:
As the dish is pointed towards northern meridian at zenith angle of 39.7 , altitude of the
centre of the beam is
a = 90.00 z = 90.00 39.70 = 50.30 1.0
Theoretical Examination
Page 6 of 27
(T6.1) Find the ratio of radii of the star in its most contracted and most expanded states 7
(R1 /R2 ).
Solution:
We first find flux ratio and then use Stefans law to compare the fluxes.
F1
m1 m2 = 2.5 log 1.0
F2
F1
= 100.4(m1 m2 ) = 100.4(3.464.08)
F2
= 1.77 1.0
Li = 4Ri2 Ti4 1.0
4R12 T14 4R22 T24
F1 = F2 = 1.0
4D2 4D2
Theoretical Examination
Page 7 of 27
F1 R2 T 4
= 12 14
F2 R T2
r2 2
R1 F1 T2
=
R2 F2 T1
T2 1
From Wiens displacement law, = . 1.0
T1 2
r 2
R1 F1 1
= 1.0
R2 F2 2
531.0 2
= 1.77
649.1
R1
= 0.890 1.0
R2
Acceptable range: 0.010.
(T6.2) Find the radii of the star (in metres) in its most contracted and most expanded states 3
(R1 and R2 ).
Solution:
R2 R1 = v P/2 2.0
9.84
R2 R1 = 12.8 103 86 400 m
2
(1 0.890)R2 = 5.441 109 m
R2 = 4.95 1010 m 0.5
(T6.3) Calculate the flux of the star, F2 , when it is in its most expanded state. 5
Solution:
To get the absolute value of flux (F2 ) we must compare it with observed flux of the
Sun.
F2
m2 m = 2.5 log
F
F2 = F 100.4(m2 m ) 3.0
L
= 2 100.4(4.08+26.72)
4a
3.826 1026 4.7863 1013
= 11 2
W m2
4(1.496 10 )
F2 = 6.51 1010 W m2 2.0
Acceptable range: 0.04 1010 W m2 .
Theoretical Examination
Page 8 of 27
Solution:
r s
R22 T24
r
L2
Dstar = = = R2 T22 2.0
4F2 F2 F2
2.898 103 m K
Wiens law: T2 = 1.0
2
2 s
2.898 103 5.670 108
Dstar = 4.95 1010
649.1 109 6.51 1010
Dstar = 9.208 1018 m = 298 pc 2.0
Acceptable range: 298 2pc (depends on truncation).
Solution:
The magnification will be given by,
fo
m0 = 1.0
fe
100
= = 100 1.0
1
The magnification is m0 = 100
Length of the telescope will be
L0 = fo + fe 1.0
= 100 + 1 = 101 cm 1.0
The telescope length will be L0 = 101 cm
Exact answer required for credit.
An introduction of a concave lens (Barlow lens) between the objective lens and the prime focus is a
common way to increase the magnification without a large increase in the length of the telescope.
A Barlow lens of focal length 1 cm is now introduced between the objective and the eyepiece to
double the magnification.
(T7.2) At what distance, dB , from the prime focus must the Barlow lens be kept in order to 6
obtain this desired double magnification?
Solution:
We use the following sign convention. Lens is the origin. Direction along the direction
Theoretical Examination
Page 9 of 27
1 1 1
of light is taken as positive. The lens formula is = (f is positive for convex
f v u
v
lens and negative for concave lens). Magnification is m = . Solutions using other
u
sign conventions are acceptable.
Let v be image distance from the Barlow lens.
1 1 1
= 1.0
fB v u
Distance of Barlow lens, dB , before the prime focus is same as the object distance,
u, in this case. 1.0
1 1 1
=
fB v dB
v v
Also, mB = 2 = = 0.5
u dB
1 2
= 1.0
dB v
1 1 2
=
1 v v
1
1 =
v
v = 1 cm
v 1 cm
dB = = = 0.5 cm 2.5
2 2
The positive sign for dB indicates that the Barlow lens was introduced 0.5 cm before
the prime focus.
Solution:
The increase in the length will be,
L = v dB 2.0
= 1.0 0.5 = 0.5 cm 2.0
Thus, the length will be increased by L = 0.5 cm
A telescope is now constructed with the same objective lens and a CCD detector placed at the
prime focus (without any Barlow lens or eyepiece). The size of each pixel of the CCD detector is
10 m.
(T7.4) What will be the distance in pixels between the centroids of the images of the two stars, 6
np , on the CCD, if they are 2000 apart on the sky?
Solution:
Plate scale at prime focus is given by,
1 1 rad
s= = = 0.206 265 arcsec/m 2.0
fo 1m
Theoretical Examination
Page 10 of 27
Solution:
The U -band is defined as (360 40) nm. Thus, the maximum, minimum and average
frequencies of the band will be,
c
max = = 9.369 1014 Hz
max
min = 7.495 1014 Hz
avg = 8.432 1014 Hz
= max min
= 1.874 1014 Hz 2.0
fst1 = 3631 100.415
= 3.631 mJy = 3.631 1029 W Hz1 m2 2.0
Now, N0 havg = fst1 A t 2.0
where, A = 1 m2 & t = 1 s
1.874 1014 3.631 1029
N0 =
6.626 1034 8.432 1014
' 12180 2.0
Exact calculation including integration is accepted with full credit (exact
answer: 12190).
Accepted range: 12180 200.
Using flat spectrum for instead of is considered a major conceptual
error, and will incur penalty of 2.0 marks.
This star is being observed in the U -band using a ground based telescope, whose primary mirror
has a diameter of 2.0 m. Atmospheric extinction in U -band during the observation is 50%. You
may assume that the seeing is diffraction limited. Average surface brightness of night sky in
U -band was measured to be 22.0 mag/arcsec2 .
Theoretical Examination
Page 11 of 27
(T8.2) What is the ratio, R, of number of photons received per second from the star to that 8
received from the sky, when measured over a circular aperture of diameter 200 ?
Solution:
Let us call sky flux per square arcsec as and total sky flux for the given aperture
as sky . Let total star flux be st .
sky = A = (1 arcsec)2 = 3.0
msky = 22.0 + 2.5 log10 1.0
sky
= 22.0 + 2.5 log10
= 22.0 2.5 log10 ()
msky = 20.76 mag 1.0
As extinction is 50% 1.0
st2 0.5st1
R= = = 0.5 10(20.7615)/2.5
sky sky
' 100 2.0
Accepted range: 5.
A student may also calculate number of photons incident per second per
metre square in this case and then compare it with the answer in the first
case to get the correct ratio.
(T8.3) In practice, only 20% of U -band photons falling on the primary mirror are detected. How 4
many photons, Nt , from the star are detected per second?
Solution:
Nt 1 m2 = N0 0.5 0.2 At 2.0
2
2.0
Nt = 12180 0.5 0.2 = 1233 1.0
2
Nt ' 3813 1.0
Accepted range (3813 50).
engine burn?
Solution:
Let apogee and perigee distances be ra and rp respectively.
rp = R + hip = (6371 + 264) km = 6635 km 1.0
ra = R + hia = (6371 + 23904) km = 30 275 km 1.0
Conservation of energy and angular momentum gives
GmM
E=
rp + ra
Total energy at perigee
1 GmM GmM
mvp 2 =E=
2 rp rp + ra
1 2 GM rp
vp = 1
2 rp rp + ra
s
GM ra
vp = 2 2.0
ra + rp rp
s
2 6.674 1011 5.972 1024 3.0275 107
=
6.635 106 (3.0275 + 6.635 106 )
= 9.929 km s1 1.0
As the engine burn is just 41.6 s, we assume that the entire impulse is applied instan-
taneously at perigee. The impulse is J = 1.73 105 kg m s1 . Note that the total
mass of MOM must include the fuel, so we have to use
m = 500 + 852 = 1352 kg. 1.0
Change in velocity due to impulse at perigee is
J 1.73 105
v = = = 128.0 m s1 1.0
m 1352
The new velocity will be given by (we use 0 symbol to denote quantities after the first
orbit-raising maneuvre)
vp0 = vp + v = 10.06 km s1 1.0
The perigee remains unchanged. So we get rp0 = rp . 1.0
Since the satellite is moving faster, the new apogee will be higher.
s
r0
vp0 = 2GM 0 0 a 0
rp (ra + rp )
rp0 2GM
1+ 0
= 0 2 2.0
ra (vp ) rp0
2 6.674 1011 5.972 1024
= = 1.188
(10.06 103 )2 6.635 106
6635
ra0 = = 35 380 km 2.0
0.188
ha = 35380 6371
= 29 009 km 1.0
Theoretical Examination
Page 13 of 27
(T9.2) Find eccentricity (e) of the new orbit after the burn and new orbital period (P ) of MOM 6
in hours.
Solution:
As seen above,
rp0 1e
0
= 0.188 = 1.0
ra 1+e
1 0.188
e= = 0.683 1.0
1 + 0.188
Acceptable range: 0.002
The new orbital semi-major axis and orbital period will be,
ra0 + rp0
a0 = 1.0
2
35380 + 6635
= = 20 933 km 1.0
s 2
a0 3
P = 2 1.0
GM
s
(2.0933 107 )3
= 2 = 30 136 s = 8.37 h 1.0
6.674 1011 5.972 1024
Acceptable range: 0.1 h
b
r
x-axis
A massive body thus behaves somewhat like a focusing lens. The light rays coming from infinite
distance beyond a massive body, and having the same impact parameter r, converge at a point
along the axis, at a distance fr from the centre of the massive body. An observer at that point
Theoretical Examination
Page 14 of 27
will benefit from huge amplification due to this gravitational focusing. The massive body in this
case is being used as a Gravitational Lensing Telescope for amplification of distant signals.
(T10.1) Consider the possibility of our Sun as a gravitational lensing telescope. Calculate the 6
shortest distance, fmin , from the centre of the Sun (in A.U.) at which the light rays can
get focused.
Solution:
The rays travelling closer to the gravitational body will bend more. Thus, we get
shortest convergence point where the rays just grazing the solar surface will meet
each other.
b
R
b
fmin
2Rsch R
b = ' 2.0
R fmin
R 2 R2 c2
fmin = = 2.0
2Rsch 4GM
(6.955 108 2.998 108 )2
= m
4 6.674 1011 1.989 1030
8.188 1013
= 8.188 1013 m = AU
1.496 1011
fmin = 547.3 AU 2.0
(T10.2) Consider a small circular detector of radius a, kept at a distance fmin centred on the 8
x-axis and perpendicular to it. Note that only the light rays which pass within a certain
annulus (ring) of width h (where h R ) around the Sun would encounter the detector.
The amplification factor at the detector is defined as the ratio of the intensity of the light
incident on the detector in the presence of the Sun and the intensity in the absence of
the Sun.
Express the amplification factor, Am , at the detector in terms of R and a.
Solution:
The following figure needs to be drawn. 1.0
h 2
R Detector
a
fmin
f2
Theoretical Examination
Page 15 of 27
The light bending from the surface of the Sun (r = R ) will intersect the detector
at its centre, as it is kept at fmin .
The boundary of the detector will be intersected by a light ray with r = R + h.
This ray will intersect the x-axis at a distance f2 .
(R + h)2
f2 = 2.0
2Rsch
Same argument as in Part 1
For small angles,
a = (f2 fmin )2
2
(R + h)2 2R h + h2
R 2Rsch
= =
2Rsch 2Rsch (R + h) R + h
' 2h 2.0
Let original intensity of the incoming radiation be I0 .
The flux at the detector in the presence of Sun is = I0 2 R h 1.0
The flux at the detector in the absence of Sun is 0 = I0 a2 1.0
The amplification is therefore
I0 2R h R
Am = = = 1.0
0 I0 a2 a
(T10.3) Consider a spherical mass distribution, such as a dark matter cluster, through which 6
light rays can pass while undergoing gravitational bending. Assume for simplicity that
for the gravitational bending with impact parameter, r, only the mass M (r) enclosed
inside the radius r is relevant.
What should be the mass distribution, M (r), such that the gravitational lens behaves
like an ideal optical convex lens ?
Solution:
2
1
r1
r2
f1 = f2
All rays should focus at the same spot. This should be evident from figure drawn on
answersheet or otherwise. 2.0
Let there be two rays with impact parameters r1 and r2 . The corresponding distances
of focus will be
ri2 ri2 c2
fi = = 2.0
2rschi 4GM (ri )
Same argument as in Part 1
Theoretical Examination
Page 16 of 27
The test mass loses energy due to the emission of gravitational waves. As a result the orbit keeps
on shrinking, until the test mass reaches the surface of the object, or in the case of a black hole,
the innermost stable circular orbit ISCO which is given by RISCO = 3Rsch , where Rsch is
the Schwarzschild radius of the black hole. This is the epoch of merger. At this point, the
amplitude of the gravitational wave is maximum, and so is its frequency, which is always twice
the orbital frequency. In this problem, we will only focus on the gravitational waves before the
merger, when Keplers laws are assumed to be valid. After the merger, the form of gravitational
waves will drastically change.
(T11.1) Consider the observed gravitational waves shown in the figure above. Estimate the time 3
period, T0 , and hence calculate the frequency, f0 , of gravitational waves just before the
epoch of merger.
Theoretical Examination
Page 17 of 27
Solution:
From the graph, just before the peak of emission, the time period of gravitational
waves is approximately (0.007 0.004) s.
T0 0.007 s 2.0
Acceptable range: 0.003 to 0.011 s
That is, the frequency of the gravitational waves is f0 142.86 Hz . 1.0
Acceptable range: 333.33 to 90.91 Hz
Answer given in terms of angular frequency with correct value and units
gains full credit.
(T11.2) For any main sequence (MS) star, the radius of the star, RMS , and its mass, MMS , are 10
related by a power law given as,
RMS (MMS )
where = 0.8 for M < MMS
= 1.0 for 0.08M MMS M
If the central object were a main sequence star, write an expression for the maximum
frequency of gravitational waves, fMS , in terms of mass of the star in units of solar masses
(MMS /M ) and .
Solution:
Since m M then, by Keplers third law
r
1 GM
forbital = 4.0
2 r3
Hence the frequency of the gravitational waves is
r
1 GM
fgrav = 2forbital = 1.0
r3
The frequency will be maximum when r = RMS . 1.0
(T11.3) Using the above result, determine the appropriate value of that will give the maximum 9
possible frequency of gravitational waves, fMS,max for any main sequence star. Evaluate
Theoretical Examination
Page 18 of 27
this frequency.
Solution:
For possible values of given in the question, the exponent 13
2 is negative. Thus,
if MMS > M , the frequency will be smaller. Thus, for highest possible frequency
coming from a main sequence star, you should take lowest possible mass i.e. = 1.0
4.0
131
MMS ( 2 )
s
1 GM
fMS,max = 3
R M
s
1 GM M
= 3 2.0
R MMS
The frequency of gravitational waves will be given by,
s
1 6.674 1011 1.989 1030 1
fMS,max = 8 3
(6.955 10 ) 0.08
fMS,max = 2.5 mHz 3.0
Answer given in terms of angular frequency with correct value and units
gains full credit.
(T11.4) White dwarf (WD) stars have a maximum mass of 1.44M (known as the Chandrasekhar 8
limit) and obey the mass-radius relation R M 1/3 . The radius of a solar mass white
dwarf is equal to 6000 km. Find the highest frequency of emitted gravitational waves,
fWD,max , if the test mass is orbiting a white dwarf.
Solution:
The maximum frequency would be when r = RWD . 1.0
We use the notation RW D for the radius of a solar mass white dwarf. Then for
white dwarfs
3 3 MW D
RW D = RW D 1.0
M
s
1 GMWD
fWD = 3
RW D
s
1 GM MW D
= 3 2.0
RW D M
For maximum frequency, we have to take highest white dwarf mass. 2.0
s
1.989 1030
fWD,max = 2.600 106 1.44
(6000 103 )3
= 2.600 106 95.96 103 1.44
fWD,max = 0.359 Hz 2.0
Answer given in terms of angular frequency with correct value and units
gains full credit.
(T11.5) Neutron stars (NS) are a peculiar type of compact objects which have masses between 8
Theoretical Examination
Page 19 of 27
1 and 3M and radii in the range 10 15 km. Find the range of frequencies of emitted
gravitational waves, fNS,min and fNS,max , if the test mass is orbiting a neutron star at a
distance close to the neutron star radius.
Solution:s
1 GMNS
fgrav = 3
RNS
The lowest possible frequency is when MNS is lowest and RNS is the highest. 2.0
For MNS = M and R = 15 km, we get
fNS,min = 1.996 kHz 2.0
Similarly, the largest possible frequency is when MNS is the largest and RNS is the
smallest. 2.0
For MNS = 3M and R = 10 km, we get
fNS,max = 6.352 kHz 2.0
Answer given in terms of angular frequency with correct value and units
gains full credit.
(T11.6) If the test mass is orbiting a black hole (BH), write the expression for the frequency of 7
emitted gravitational waves, fBH , in terms of mass of the black hole, MBH , and the solar
mass M .
Solution:
For black holes, we have to consider RISCO . 2.0
Hence the equation will be,
s
1 GM M
fBH = 3 2.0
27Rsch MBH
M
fBH = 4.396 kHz 3.0
MBH
Answer given in terms of angular frequency with correct value and units
gains full credit.
(T11.7) Based only on the time period (or frequency) of gravitational waves before the epoch of 5
merger, determine whether the central object can be a main sequence star (MS), a white
dwarf (WD), a neutron star (NS), or a black hole (BH). Tick the correct option in the
Summary Answersheet. Estimate the mass of this object, Mobj , in units of M .
Solution:
We found the frequency of the LIGO-detected wave to be 166.67 Hz just before
merger. As per our analysis above, only black holes can lead to emission in this
frequency range.
Black Hole 2.0
By using corresponding expression,
4396
Mobj = M 31M 3.0
142.86
Any answer between 13 to 50 will get full credit.
Theoretical Examination
Page 20 of 27
(T12) Exoplanets
Two major methods of detection of exoplanets (planets around stars other than the Sun) are the
radial velocity (or so-called wobble) method and the transit method. In this problem, we find
out how a combination of the results of these two methods can reveal a lot of information about
an orbiting exoplanet and its host star.
Throughout this problem, we consider the case of a planet of mass Mp and radius Rp moving in a
circular orbit of radius a around a star of mass Ms (Ms Mp ) and radius Rs . The normal to the
orbital plane of the planet is inclined at angle i with respect to the line of sight (i = 90 would
mean edge on orbit). We assume that there is no other planet orbiting the star and Rs a.
Wobble Method:
When a planet and a star orbit each other around their barycentre, the star is seen to move
slightly, or wobble, since the centre of mass of the star is not coincident with the barycentre of
the star-planet system. As a result, the light received from the star undergoes a small Doppler
shift related to the velocity of this wobble.
The line of sight velocity, vl , of the star can be determined from the Doppler shift of a known
spectral line, and its periodic variation with time, t, is shown in the schematic diagram below.
In the diagram, the two measurable quantities in this method, namely, the orbital period P and
maximum line of sight velocity v0 are shown.
vl
P
v0
(T12.1) Derive expressions for the orbital radius (a) and orbital speed (vp ) of the planet in terms 3
of Ms and P .
Solution:
Keplers law:
GMs 2 1/3
a= P 1.0
4 2
Gravitational force provides centripetal acceleration:
r
GMs
vp = 0.5
a
2GMs 1/3
vp = 1.5
P
Theoretical Examination
Page 21 of 27
(T12.2) Obtain a lower limit on the mass of the planet, Mp, min in terms of Ms , v0 and vp . 4
Solution:
Momentum conservation:
Mp v p = M s v s 1.5
Observed quantity is v0 = vs sin i. Thus,
Ms vs sin i Ms v 0
Mp, min = Mp sin i = =
vp vp
Transit Method:
As a planet orbits its host star, for orientations of the orbital plane that are close to edge-on
(i 90 ), it will pass periodically, or transit, in front of the stellar disc as seen by the observer.
This would cause a tiny decrease in the observed stellar flux which can be measured. The schematic
diagram below (NOT drawn to scale) shows the situation from the observers perspective and the
resulting transit light curve (normalised flux, f , vs time, t) for a uniformly bright stellar disc.
Rs
bRs
1 2 3 4
Rp
f
1
tF
tT
If the inclination angle i is exactly 90 , the planet would be seen to cross the stellar disc along a
diameter. For other values of i, the transit occurs along a chord, whose centre lies at a distance
bRs from the centre of the stellar disc, as shown. The no-transit flux is normalised to 1 and the
maximum dip during the transit is given by .
Theoretical Examination
Page 22 of 27
The four significant points in the transit are the first, second, third and fourth contacts, marked
by the positions 1 to 4, respectively, in the figure above. The time interval during the second and
third contacts is denoted by tF , when the disc of the planet overlaps the stellar disc fully. The
time interval between the first and fourth contacts is denoted by tT .
These points are also marked in the schematic diagram below showing a side-on view of the
orbit (NOT drawn to scale).
i
4
3
a Observer
1 2
(T12.3) Find the constraint on i in terms of Rs and a for the transit to be visible at all to the 2
distant observer.
Solution:
bRs = a cos i 1.0
Therefore, for visibility, 0 b 1 i cos1 (Rs /a) 1.0
Solution:
Blackbody brightness is proportional to area. Since the observer is far away from
the star-planet system, size of silhouette of planet on stellar disc is independent of a.
2
Rp
= 1.0
Rs
Solution:
Circular orbit uniform orbital speed
t a
= = 1.0
P 2a 2
where is the angle subtended by the planet at the centre of the star during transit
(over time t).
Theoretical Examination
Page 23 of 27
14 l14 3
Rs + Rp Rs Rp 23 l23
a
bRs
1 2 3 4
1 2
l14 /2 l23 /2
l23 /2
sin(23 /2) = 2.0
a
s
2
l23 Rs Rp
23 = 2 sin1 = 2 sin1 1 b2 1.0
2a a Rs
s
2
P P Rs Rp
tF = 23 = sin1 1 b2 1.0
2 a Rs
Similarly,
s
2
P Rs Rp
tT = sin1 1+ b2 1.0
a Rs
(T12.6) In the approximation of an orbit much larger than the stellar radius, show that the 5
parameter b is given by
2 1/2
tF
1+
tT
b = 1 + 2
2
tF
1
tT
Solution:
Since Rs a, use sin1 x x. 2.0
s
Rp 2
P Rs
tT 1+ b2
a Rs
s
2
P Rs Rp
tF 1 b2
a Rs
Theoretical Examination
Page 24 of 27
Dividing, and putting Rp /Rs = ,
" #1/2
tF (1 )2 b2
= 1.0
tT (1 + )2 b2
2 1/2
2 1/2
tF 2
tF
(1 )2 (1 + ) 1+
tT tT
b= = 1 + 2 2.0
2 2
tF tF
1 1
tT tT
Expressions lacking the use of approximation Rs /a 1, but otherwise
correct will get a penalty of 2.0.
Use of approximation with proper justification at a later stage than at
the first step will get full credit.
(T12.7) Use the result of part (T12.6) to obtain an expression for the ratio a/Rs in terms of 3
measurable transit parameters, using a suitable approximation.
Solution:
s
2
P Rs Rp
tT = 1+ b2 1.0
a Rs
a 2P 1/4
= 2.0
Rs (t2T t2F )1/2
Expressions lacking the use of approximation Rs /a 1, but otherwise
correct will get a penalty of 1.0.
If penalty has already been imposed in part (T12.6), no further penalty
for lack of approximation.
(T12.8) Combine the results of the wobble method and the transit method to determine the 6
Ms
stellar mean density s in terms of tT , tF , and P .
4Rs3 /3
Theoretical Examination
Page 25 of 27
Solution:
From part (T12.7)
2P 1/4
a = Rs
(t2T t2F )1/2
From part (T12.1)
GMs 2 1/3
a= P
4 2
Combining,
" #3
GMs 2 2P 1/4
P = Rs 3.0
4 2 (t2T t2F )1/2
Identifying the two equations for combining gets credit. No credit for
writing only one equation or irrelevant equations.
Ms 3 4 2 8P 3 3/4
s = 2.0
4Rs3 /3 4 P 2 G 3 (t2T t2F )3/2
24 P ()3/4
s = 1.0
2 G (t2T t2F )3/2
Rocky or gaseous:
Let us consider an edge-on (i = 90 ) star-planet system (circular orbit for the planet), as seen
from the Earth. It is known that the host star is of mass 1.00M . Transits are observed with
a period (P ) of 50.0 days and total transit duration (tT ) of 1.00 hour. The transit depth () is
0.0064. The same system is also observed in the wobble method to have a maximum line of sight
velocity of 0.400 m s1 .
(T12.9) Find the orbital radius a of the planet in units of AU and in metres. 2
Solution:
From Keplers third law (with same mass of host star):
P 2/3
a
= 1.0
a P
2/3
50.0
a= 1 AU = 0.266 AU 0.5
365.242
= 0.266 1.496 1011 m = 3.97 1010 m 0.5
Solution:
Edge-on b = 0 1.0
" #1/2
tF (1 )2 b2 1
= = = 0.8519 1.0
tT (1 + )2 b2 1+
Theoretical Examination
Page 26 of 27
(T12.11) Obtain the mass Mp and radius Rp of the planet in terms of the mass (M ) and radius 8
(R ) of the Earth respectively. Is the composition of the planet likely to be rocky or
gaseous? Tick the box for ROCKY or GASEOUS in the Summary Answersheet.
Solution:
From parts (T12.1) and (T12.9)
r s
GM 6.674 1011 1.989 1030
vp = = = 57.798 km s1 1.0
a 3.97 1010
Assumption of small planet (Mp Ms ) is valid because is very small; less dense
planet would make the assumption stronger!
Solution:
f Transit no. 1 f Transit no. 2
t t
t t
f Transit no. 5
(T12.13) Throughout the problem we have considered a uniformly bright stellar disc. However, 2
real stellar discs have limb darkening. Draw a schematic transit light curve when limb
darkening is present in the host star.
Solution:
f
Non-flat bottom with central minimum gets 2.0. Curvature of ingress and
egress are tolerated.
Data Analysis Examination
Page 1 of 4
4 2
() = t where t =
B2
where B is the orbital period of the binary, 0 is the intrinsic spin period of the pulsar and is the radius
of the orbit.
The following table gives one such set of measurements of and at different heliocentric epochs, ,
expressed in truncated Modified Julian Days (tMJD), i.e. number of days since MJD = 2,440,000.
No. T P a
(tMJD) (s) (m s-2)
1 5740.654 7587.8889 - 0.92 0.08
2 5740.703 7587.8334 - 0.24 0.08
3 5746.100 7588.4100 - 1.68 0.04
4 5746.675 7588.5810 + 1.67 0.06
5 5981.811 7587.8836 + 0.72 0.06
6 5983.932 7587.8552 - 0.44 0.08
7 6005.893 7589.1029 + 0.52 0.08
8 6040.857 7589.1350 + 0.00 0.04
9 6335.904 7589.1358 + 0.00 0.02
By plotting () as a function of (), we can obtain a parametric curve. As evident from the relations
above, this curve in the period-acceleration plane is an ellipse.
In this problem, we estimate the intrinsic spin period, 0 , the orbital period, B , and the orbital radius, ,
by an analysis of this data set, assuming a circular orbit.
(D1.1) Plot the data, including error bars, in the period-acceleration plane (mark your graph as D1.1). 7
(D1.2) Draw an ellipse that appears to be a best fit to the data (on the same graph D1.1). 2
(D1.3) From the plot, estimate 0 , t and t , including error margins. 7
(D1.4) Write expressions for B and in terms of 0 , t , t . 4
(D1.5) Calculate approximate value of B and based on your estimations made in (D1.3), including 6
error margins.
(D1.6) Calculate orbital phase, , corresponding to the epochs of the following five observations in 4
the above table: data rows 1, 4, 6, 8, 9.
(D1.7) Refine the estimate of the orbital period, B , using the results in part (D1.6) in the followin g
way:
(D1.7a) First determine the initial epoch, 0 , which corresponds to the nearest epoch of zero 2
orbital phase before the first observation.
(D1.7b) The expected time, calc , of the estimated orbital phase angle of each observation is 7
given by,
calc = 0 + ( + ) B ,
360
Data Analysis Examination
Page 2 of 4
where n is the number of full cycle of orbital phases that may have elapsed between 0
and (or calc ). Estimate n and calc for each of the five observations in part (D1.6).
Note down difference OC between observed and calc . Enter these calculations in
the table given in the Summary Answersheet.
(D1.7c) Plot OC against (mark your graph as D1.7). 4
(D1.7d) Determine the refined values of the initial epoch, 0 ,r, and the orbital period, B,r. 7
The composite graphic1 below shows multiple snapshots of the Moon taken at different times during the
total lunar eclipse, which occurred in this month. For each shot, the centre of frame was coinciding with
the central north-south line of umbra.
For this problem, assume that the observer is at the centre of the Earth and angular size refers to angular
diameter of the object / shadow.
1
Credit: NASAs Scientific Visualization Studio
Data Analysis Examination
Page 3 of 4
Calculate the angular size of the penumbra, penumbra , in terms of Moon . Assume the observer
to be at the centre of the Earth.
(D2.6) Let Earth be angular size of the Earth as seen from the centre of the Moon. Calculate the angular 5
size of the Moon, Moon , as would be seen from the centre of the Earth on the eclipse day in
terms of Earth .
(D2.7) Estimate the radius of the Moon, Moon , in km from the results above. 3
(D2.8) Estimate the shortest distance, perigee , and the farthest distance, apogee , to the Moon. 4
(D2.9) Use appropriate data from September 10 to estimate the distance, Sun , to the Sun from the 10
Earth.
Light curves of all type Ia supernovae can be fit to the same model light curve, when they are scaled
appropriately. In order to achieve this, we first have to express the light curves in the reference frame of
the host galaxy by taking care of the cosmological stretching/dilation of all observed time intervals, obs ,
by a factor of (1 + ). The time interval in the rest frame of the host galaxy is denoted by gal .
Data Analysis Examination
Page 4 of 4
The rest frame light curve of a supernova changes by two magnitudes compared to the peak in a time
interval 0 after the peak. If we further scale the time intervals by a factor of s (i.e. = gal ) such
that the scaled value of 0 is the same for all supernovae, the light curves turn out to have the same
shape. It also turns out that is related linearly to the absolute magnitude, peak , at the peak luminosity
for the supernova. That is, we can write
= + peak ,
where and are constants. Knowing the scaling factor, one can determine absolute magnitudes of
supernovae at unknown distances from the above linear equation.
The table below contains data for three supernovae, including their distance moduli, (for the first two),
their recession speed, , and their apparent magnitudes, obs , at different times. The time obs
peak shows number of days from the date at which the respective supernova reached peak brightness.
The observed magnitudes have already been corrected for interstellar as well as atmospheric extinction.
(D3.1) Compute gal values for all three supernovae, and fill them in the given blank boxes in the 15
data tables on the BACK side of the Summary Answersheet. On a graph paper, plot the points
and draw the three light curves in the rest frame (mark your graph as D3.1).
(D3.2) Take the scaling factor, 2 , for the supernova SN2006IS to be 1.00. Calculate the scaling 5
factors, 1 and 3 , for the other two supernovae SN2006TD and SN2005LZ, respectively, by
calculating 0 for them.
(D3.3) Compute the scaled time differences, s , for all three supernovae. Write the values for s in 14
the same data tables on the Summary Answersheet. On another graph paper, plot the points and
draw the 3 light curves to verify that they now have an identical profile (mark your graph as
D3.3).
(D3.4) Calculate the absolute magnitudes at peak brightness, peak ,1, for SN2006TD and peak ,2 , for 6
SN2006IS. Use these values to calculate and .
(D3.5) Calculate the absolute magnitude at peak brightness, peak,3 , and distance modulus, 3, for 4
SN2005LZ.
(D3.6) Use the distance modulus 3 to estimate the value of Hubble's constant, 0 . Further, estimate 6
the characteristic age of the universe, H .
Data Analysis Examination
Page 1 of 20
For a pulsar in a binary orbit, the measured pulsar spin period (P ) and the measured line-of-sight
acceleration (a) both vary systematically due to orbital motion. For circular orbits, this variation
can be described mathematically in terms of orbital phase (0 2) as,
2P0 r
P () = P0 + Pt cos where Pt =
cPB
4 2 r
a() = at sin where at =
PB2
where PB is the orbital period of the binary, P0 is the intrinsic spin period of the pulsar and r is
the radius of the orbit.
The following table gives one such set of measurements of P and a at different heliocentric epochs,
T , expressed in truncated Modified Julian Days (tMJD), i.e. number of days since MJD =
2,440,000.
T P a
No. (tMJD) (s) (m s2 )
1 5740.654 7587.8889 0.92 0.08
2 5740.703 7587.8334 0.24 0.08
3 5746.100 7588.4100 1.68 0.04
4 5746.675 7588.5810 +1.67 0.06
5 5981.811 7587.8836 +0.72 0.06
6 5983.932 7587.8552 0.44 0.08
7 6005.893 7589.1029 +0.52 0.08
8 6040.857 7589.1350 +0.00 0.04
9 6335.904 7589.1358 +0.00 0.02
By plotting a() as a function of P (), we can obtain a parametric curve. As evident from the
relations above, this curve in the period-acceleration plane is an ellipse.
In this problem, we estimate the intrinsic spin period, P0 , the orbital period, PB , and the orbital
radius, r, by an analysis of this data set, assuming a circular orbit.
(D1.1) Plot the data, including error bars, in the period-acceleration plane (mark your graph as 7
D1.1).
Solution:
Graph Number : D1.1
Data Analysis Examination
Page 2 of 20
(D1.2) Draw an ellipse that appears to be a best fit to the data (on the same graph D1.1). 2
Solution:
See above
Solution:
Values are determined from lengths of axes of ellipse and mid-point of P -axis.
Error margins may be determined by estimating extreme ellipses covering the points
with errorbars. Any reasonable method of estimating error margins will be accepted.
Marking table:
Parameter Half credit Full credit Half credit
Minimum Minimum Maximum Maximum
Pt (s) 0.59 0.63 0.71 0.75
Pt (s) 0.01 0.02 0.04 0.05
P0 (s) 7588.38 7588.43 7588.53 7588.58
P0 (s) 0.01 0.02 0.04 0.05
at (m s2 ) 1.61 1.65 1.77 1.81
2
at (m s ) 0.04 0.05 0.07 0.08
Wrong values due to wrong/poor plot/fit in (D1.1) and (D1.2) WILL
BE penalised.
Error estimation is based on graph drawing. Quoted values corre-
spond to the envelope of possible ellipses drawn to include all points
with errorbars. Any reasonable method to estimate error to be given
credit.
Solution:
We can easily recover the orbital period (PB ) and the radius of the orbit (r) in a
circular orbit:
4 2
at = 2 r
PB
PB2 at
r= 1.0
4 2
2P0 r
Pt =
PB c
2P0 PB2 at P0 PB at
= 2
=
PB c 4 2c
Data Analysis Examination
Page 4 of 20
Pt 2c
PB = 1.0
P0 at
Pt 2c 2
at
r= 2 1.0
4 P0 at
2 2
Pt c
r= 1.0
P0 at
Alternative algebraic routes accepted. Each of PB and r carry 2.0 marks.
(D1.5) Calculate approximate value of PB and r based on your estimations made in (D1.3), 6
including error margins.
Solution:
Pt 2c
PB =
P0 at
0.67 2 2.998 108
= s
7588.48 1.71
= 96 260 s = 1.125 70 d 1.0
s
Pt 2 P0 2 at 2
PB = PB + + 1.0
Pt P0 at
s
0.02 2
2
0.06 2
0.02
= 1.12570 + + d
0.67 7588.48 1.71
= 1.12570 0.0461 ' 0.052 d
PB = (1.13 0.05) d 1.0
2
c2
Pt
r=
P0 at
2 2
2.998 108
0.67
= m
7588.48 1.71
r = 4.097 39 108 m = 2.738 90 103 AU 1.0
s
2Pt 2 2P0 2 at 2
r = r + + 1.0
Pt P0 at
s 2
2 0.02 2 0.06 2
3 2 0.02
= 2.738 90 10 + + AU
0.67 7588.48 1.71
= 2.738 90 103 0.069 25 AU ' 0.19 103 AU
r = (2.74 0.19) 103 AU 1.0
Errors in PB and r can be also estimated as maximum possible (worst case) error.
In such case, errors would be about 1.5 times the standard error calculated above
(PB = 0.07, r = 0.25).
(D1.6) Calculate orbital phase, , corresponding to the epochs of the following five observations 4
Data Analysis Examination
Page 5 of 20
Solution:
Using these newly determined orbital parameters, we can calculate the angular orbital
phase for each data point, i.e., for each pair of acceleration and period measured (P ,
a).
1 a Pt
= tan
at P P0
Care has to be taken to choose the value of the phase from among , , 2 ,
depending on the sign of cos and sin .
Sr. T P a
no. (tMJD) (s) (m s2 )
1 5740.654 7587.8889 0.92 148.62
4 5746.675 7588.5810 +1.67 278.77
6 5983.932 7587.8552 0.44 164.57
8 6040.857 7589.1350 +0.00 0.00
9 6335.904 7589.1358 +0.00 0.00
Credit for each correct value: 1.0 for first three, 0.5 for last two.
Credit for or 2 is 0.5 per value.
All values wrong due to wrong expression for gets a maximum of
1.0 mark.
Values in radians accepted.
(D1.7) Refine the estimate of the orbital period, PB , using the results in part (D1.6) in the
following way:
(D1.7a) First determine the initial epoch, T0 , which corresponds to the nearest epoch 2
of zero phase before the first observation.
Solution:
T1 T0 1 1
= T0 = T1 PB 1.0
PB 2 2
148.62
T0 = 5740.654 1.12570 tMJD
360
T0 = 5740.189 tMJD 1.0
Tolerance: 0.002 tMJD. Using P0 instead of PB gets zero.
(D1.7b) The expected time, Tcalc , of the estimated phase of each observation is given by 7
Tcalc = T0 + n + PB ,
360
where n is the number of full cycle of orbital phases elapsed between T0 and
Tcalc . Estimate n and Tcalc for each of the five observations in part (D1.6). Note
down difference TOC between observed T and Tcalc . Enter these calculations
in the table given in the Summary Answersheet.
Solution:
Data Analysis Examination
Page 6 of 20
Tcalc = T0 + n + PB
360
where n = Integer part of [(T T0 )/PB ].
Solution:
Graph Number : D1.7
(D1.7d) Determine the refined values of the initial epoch, T0,r , and the orbital period, 7
PB,r .
Data Analysis Examination
Page 7 of 20
Solution:
A linear fit to the plot of TOC vs n gives the offset of period per cycle
(slope) and the shift in the zero-phase point (intercept). 2.0
This concept, which may be evident in the subsequent calculation,
gains the credit, explicit statement is not necessary.
From a linear fit,
Slope = 0.000 43 d/n Intercept = 0.010 d 3.0
The composite graphic1 below shows multiple snapshots of the Moon taken at different times
during the total lunar eclipse, which occurred in this month. For each shot, the centre of frame
was coinciding with the central north-south line of umbra.
For this problem, assume that the observer is at the centre of the Earth and angular size refers to
angular diameter of the relevant object / shadow.
Solution:
From the table we see that the angular size of Moon is smallest close to the New
Moon day. Thus, the answer is New Moon . 3.0
Justification is NOT necessary for full credit.
(D2.2) In September 2015, the ascending node of lunar orbit with respect to the ecliptic is closest 4
to
New Moon / First Quarter / Full Moon / Third Quarter.
Tick the correct answer in the Summary Answersheet. No justification for your answer
is necessary.
Solution:
As there is an eclipse happening in this month, the lunar nodes are close to Full
Moon day and New Moon day. Next we notice that lowest declination of Moon is
just 18 . This means that after the New Moon day, the orbit of Moon is above the
ecliptic. In other words, the ascending node is near the New Moon . 4.0
Justification is NOT necessary for full credit.
(D2.3) Estiamte the eccentricity, e, of the lunar orbit from the given data. 4
Solution:
The largest angular size of the Moon in the ephemerides is 2008.300 and the smallest
angular size is 1763.700 . The distance is inversely proportional to the angular size.
Hence ratio of distance at perigee to the distance at apogee is:
rperigee a0 1e
Ratio = = 2.0
rapogee 1+e a
0
1
Credit: NASAs Scientific Visualization Studio
Data Analysis Examination
Page 10 of 20
1e 1763.7
= = 0.87821
1+e 2008.3
1 0.87821
e= = 0.064846
1 + 0.87821
e ' 0.065 2.0
Rounding off is done to account for the fact that our data are not contin-
uous, hence exact angular sizes at perigee and apogee are not known.
A non-rounded answer will also receive full credit.
(D2.4) Estimate the angular size of the umbra, umbra , in terms of the angular size of the 8
Moon, Moon . Show your working on the image given on the backside of the Summary
Answersheet.
Solution:
The following construction is shown. 5.0
(D2.5) The angle subtended by the Sun at Earth on the day of the lunar eclipse is known to be 9
Sun = 1915.000 . In the figure below, S1 R1 and S2 R2 are rays coming from diametrically
opposite ends of the solar disk. The figure is not to scale.
S1
R2
Sun
R1
S2
Earth
Moon
Calculate the angular size of the penumbra, penumbra , in terms of Moon . Assume the
observer to be at the centre of the Earth.
Solution:
The following diagram needs to be drawn.
S1
1 E
P Sun
2
S2
B R
J O
Q
Moons path
Angular size of umbra is umbra = 2]BOQ
Angular size of penumbra is penumbra = 2]AOQ
We have
QA = QP + P A
= OE + P E tan 1
Data Analysis Examination
Page 12 of 20
R + dMoon 1 since P A P E
Sun
R + dMoon since 1 2 Sun /2 2.0
2
and
QB = QP P B
= OE P E tan 2
R dMoon 2
Sun
R dMoon 2.0
2
1 QA QA
penumbra = 2]AOQ = 2 tan 2 since QA OQ
OQ OQ
Sun
R + dMoon
=2 2 = 2R + 1.0
Sun
dMoon dMoon
and
1 QB QB
umbra = 2]BOQ = 2 tan 2
OQ OQ
Sun
R dMoon
=2 2 = 2R 1.0
Sun
dMoon dMoon
Subtracting,
penumbra umbra = 2Sun penumbra = umbra + 2Sun 1.0
We have,
umbra = 2.76Moon and Sun = 1915.000
From the given data, Moon = 2008.300 . 1.0
Therefore,
1915.0
penumbra = 2.76Moon + 2 Moon
2008.3
Alternative solution:
In the figure below, rays HEA and IF C are coming from one edge of solar disk
and rays HF D and GEB are coming from the opposite edge. The observer (O) is
assumed to be at the centre of the Earth. The Moon travels along the path ABCD
during the course of eclipse.
Data Analysis Examination
Page 13 of 20
G
A
E
Sun Sun
B
J Sun O
H
C
Sun Sun
F
D
I
Moons path
A
E
Sun
B
J Sun O
C
Sun
F
D
Moons path
From figure,
]AEB = ]GEH = ]HF I = ]DF C = ]EJF = Sun 3.0
umbra = ]BOC = 2.76Moon 1.0
penumbra = ]AOD 1.0
]AOD = ]AOB + ]BOC + ]COD
]AOB = ]AEB 1.0
]COD = ]CF D 1.0
penumbra ' ]AEB + umbra + ]CF D
= 2Sun + 2.76Moon = 2 1915.000 + 2.76 2008.300
penumbra = 9372.900 = 4.67Moon 2.0
(D2.6) Let Earth be angular size of the Earth as seen from the centre of the Moon. Calculate 5
the angular size of the Moon, Moon , as would be seen from the centre of the Earth on
the eclipse day in terms of Earth .
Solution:
From the Moon,
2R
Earth = 1.0
dMoon
Data Analysis Examination
Page 14 of 20
Alternative solution:
Let us say that the Moon is at position of B. Thus, angular size of Earth as seen
from this position will be, (see figure in the previous part)
Earth = ]EBF = ]BF D 2.0
= ]BF C + ]CF D
' umbra + Sun 1.0
The angular size of the Full Moon on 28 September as seen in the table is 2008.300 .
Earth = 2.76 2008.300 + 1915.000 = 7453.000
Earth
Moon = 0.269Earth = 2.0
3.72
(D2.7) Estimate the radius of the Moon, RMoon , in km from the results above. 3
Solution:
Thus, the radius of Moon will be,
R
RMoon =
3.72
6371
RMoon = 1.0
3.72
RMoon ' 1713 km 2.0
Acceptable range: 20 km.
(D2.8) Estimate the shortest distance, rperigee , and the farthest distance, rapogee , to the Moon. 4
Solution:
The shortest and longest distances will be,
2 1713 206265
rperigee =
2008.3
rperigee = 3.52 105 km 2.0
2 1713 206265
rapogee =
1763.7
rapogee = 4.01 105 km 2.0
(D2.9) Use appropriate data from September 10 to estimate the distance, dSun , to the Sun from 10
the Earth.
Data Analysis Examination
Page 15 of 20
Solution:
Moon
M
M
S S E E
Sun Earth
Moons orbit
On September 10, phase of Moon is 0.097 and elongation of Moon is 36.2 . Angular
size of the Moon on this day is 1792.000 . Therefore, distance to Moon (from Earth)
on September 10 is
2 1713 206265
dMoon,10 =
1792.0
= 3.94 105 km 2.0
Let ]EM S = M
]ESM = S
]SEM = E
E = 36.2 1.0
1 + cos M
phase = 2.0
2
1
M = cos (2 phase 1)
= cos1 (2 0.097 1) = cos1 (0.806)
= 143.71 2.0
S = 180 E M
= 180 36.2 143.71
= 0.09 1.0
Now using sine rule,
dSun sin M
= 2.0
dMoon,10 sin S
3.94 108 sin 143.71
dSun =
sin 0.09
dSun = 1.48 1011 m 1.0
Data Analysis Examination
Page 16 of 20
(D3.1) Compute tgal values for all three supernovae, and fill them in the given blank boxes 15
in the data tables on the BACK side of the Summary Answersheet. On a graph paper,
plot the points and draw the three light curves in the rest frame (mark your graph as
D3.1).
Solution:
Redshifts for the three supernovae are z1 = 0.0151, z2 = 0.0314 and z3 = 0.0402. 1.5
(d) (mag) (d) (d) (d) (mag) (d) (d) (d) (mag) (d) (d)
15.00 19.41 14.78 20.00 15.00 18.35 14.54 14.54 15.00 20.18 14.42 17.03
10.00 17.48 9.85 13.34 10.00 17.26 9.70 9.70 10.00 18.79 9.61 11.35
5.00 16.12 4.93 6.67 5.00 16.42 4.85 4.85 5.00 17.85 4.81 5.68
0.00 15.74 0.00 0.00 0.00 16.17 0.00 0.00 0.00 17.58 0.00 0.00
5.00 16.06 4.93 6.67 5.00 16.41 4.85 4.85 5.00 17.72 4.81 5.68
10.00 16.72 9.85 13.34 10.00 16.82 9.70 9.70 10.00 18.24 9.61 11.35
15.00 17.53 14.78 20.00 15.00 17.37 14.54 14.54 15.00 18.98 14.42 17.03
20.00 18.08 19.70 26.67 20.00 17.91 19.39 19.39 20.00 19.62 19.23 22.70
25.00 18.43 24.63 33.34 25.00 18.39 24.24 24.24 25.00 20.16 24.03 28.38
30.00 18.64 29.56 40.01 30.00 18.73 29.09 29.09 30.00 20.48 28.84 34.06
(D3.2) Take the scaling factor, s2 , for the supernova SN2006IS to be 1.00. Calculate the scaling 5
factors, s1 and s3 , for the other two supernovae SN2006TD and SN 2005LZ, respectively,
by calculating t0 for them.
Solution:
From the graph D3.1, SN2006IS took 22.0 d to fade by 2 magnitudes.
That is, t0 (SN2006IS) = 22.0 d.
Similarly, t0 (SN2006TD) = 16.4 d.
And t0 (SN2005LZ) = 18.8 d.
Acceptable range: 1.0 days 3.0
Thus, stretching factors for these two supernovae are
22.2
s1 = = 1.354
16.4
22.2
s3 = = 1.181 2.0
18.8
(D3.3) Compute the scaled time differences, ts , for all three supernovae. Write the values for 14
ts in the same data tables on the Summary Answersheet. On another graph paper,
plot the points and draw 3 light curves to verify that they now have an identical profile
(mark your graph as D3.3).
Solution:
Filling the scaled values in the fourth column of the table (ts in table above) 3.5
10.5
(D3.4) Calculate the absolute magnitudes at peak brightness, Mpeak,1 , for SN2006TD and 6
Mpeak,2 , for SN2006IS. Use these values to calculate a and b.
Solution:
To get a and b,
Mpeak,1 = mpeak,1 1 = 15.74 34.27 mag
= 18.53 mag
Mpeak,2 = mpeak,2 2 = 16.17 35.64 mag
= 19.47 mag 2.0
s1 s2 1.354 1 0.354
b= = mag1 = mag1
Mpeak,1 Mpeak,2 18.53 (19.47) 0.94
Data Analysis Examination
Page 20 of 20
(D3.5) Calculate the absolute magnitude at peak brightness, Mpeak,3 , and distance modulus, 3 , 4
for SN2005LZ.
Solution:
s3 = a + bMpeak,3
s3 a 1.181 8.325 7.144
Mpeak,3 = = mag = mag
b 0.3762 0.3762
Mpeak,3 = 18.99 mag 2.0
(D3.6) Use the distance modulus 3 to estimate the value of Hubbles constant, H0 . Further, 6
estimate the characteristic age of the universe, TH .
Solution:
Distance to SN2005LZ is
3 3
d = 10( 5 +1) pc = 10( 5 +16) Mpc
3
36.57
= 10( 5 5) Mpc = 102.314 Mpc
' 206 Mpc
cz3 12060
H0 = = km s1 Mpc1
d3 206
H0 = 58.5 km s1 Mpc1 4.0
1 3.086 1022
TH = = yr
H0 58.5 103 3.156 107
TH = 16.7 Gyr 2.0
Extra factor of 2/3 allowed in the value of TH .
Observational Examination (OM)
Page 1 of 1
Instructions
(OM1) Mark any 5 (five) of the following stars on the map by putting a circle (O) around the appropriate star and 20
writing its code next to it. If you mark more than 5 stars, only the first 5 in serial order will be considered.
(OM2) Mark location of any 3 (three) of the following galaxies on the map by putting a + sign at appropriate 15
place in the map and writing its code next to it. If you mark more than 3 galaxies, only the first 3 in serial
order will be considered.
(OM4) Show position of Autumnal Equinox (descending node of the ecliptic) on the map by a + sign and label 5
it as A.
(OM5) Draw local meridian for Bhubaneswar on Winter Solstice day (22nd December) at local midnight and label 5
it as M.
Observational Examination (OP)
Examination Page 1 of 1
Instructions
A. Logistical
(1) Total duration of the examination is 30 minutes.
(2) After you take your designated seat inside planetarium dome, you will be given 5 minutes to read all
the questions. Planetarium will be turned on only after this time. You will be given 5 minutes for
familiarisation with the sky after it has been turned on.
(3) At your seat you will find a writing board, a torch and a pen. Please use this pen to mark final answers
on the maps provided (Map 1 and Map 2). Answers marked by other writing instruments will not be
considered.
(4) At any time during the examination, you may stand up at your place or turn around. However, you
are not allowed to move out of your place.
(5) You must maintain strict silence throughout the examination.
(6) After your examination is over, keep writing board, torch and pen at the same place. Fold your map,
put it back inside the envelope and hand it over to the volunteers.
B. Academic
(1) Your Summary Answersheet consists of an A4 size sheet on which sky maps are printed on both
sides (Map 1 and Map 2). Please make sure that you write your contestant code on both sides of the
sheet in the boxes provided for the purpose.
(2) To mark the answers put a + sign at appropriate place on the map and write code as given in the
question next to it.
(3) Some inaccuracy in marking of a position will be tolerated, with an appropriate penalty.
(4) Start of a new part of each question will be announced. Simultaneously, a corresponding code will
be displayed on the dome. Please pay attention.
Observational Examination (OP)
Page 1 of 1
(OP1) Eight well known historical supernovae will appear in the projected sky one at a time (not necessarily in 40
chronological order). You have to identify the appropriate map (Map 1 / Map 2) where a particular
supernova belongs and mark it in the corresponding map with + sign and write codes S1 to S8
besides it.
Each supernova code will be projected on dome for 10 seconds, followed by appearance of supernova for
60 seconds and then 20 seconds for you to mark the answers.
(OP1.1) For S1, S2, S3, S4 and S5, the projected sky corresponds to the sky as seen from Rio de Janeiro
on the midnight of 21st May.
(OP1.2) For S6, S7 and S8, the projected sky corresponds to the sky as seen from Beijing on the midnight
of 20th November. There will be a gap of two minute after S5 for change over and adaptation
to new sky.
(OP2) We are now projecting sky of another planet. The sky will be slowly rotated for 5 minutes . Identify the 10
visible celestial pole of this planet and mark it with a + sign and label it as P on the appropriate map
(Map 1 / Map 2).
Observational Examination (OT)
Page 1 of 1
Instructions
(1) Total duration of the examination is 25 minutes. There will be a bell every five minutes.
(2) The examiner at your station will give you a writing board, a torch and pen. Please use only the given pen
to mark final answers on the map provided in the Summary Answersheet. Answers marked using any
other writing instrument will not be considered.
(3) After your examination is over, return the writing board, torch, pen and the envelope including your
answersheet to the examiner.
Observational Examination (OT)
Examination Page 1 of 1
When you arrive at your observing station, DO NOT disturb the telescope before attempting the first
question (OT1).
(OT1) The telescope is already set to a deep sky object. Identify the object and tick the correct box in the 10
Summary Answersheet.
Note: You can use any technique to identify the object. However, if you disturb the telescope, you will
NOT be helped to bring it back to the original position.
(OT2)
(OT2.1) Point the telescope to M45. Show the object to the examiner. 5
Note: 1. After 5 minutes, 1 mark will be deducted for a delay of every minute (or part thereof)
in pointing the telescope.
2. You have a single chance to be evaluated. If your pointing is incorrect the examiner
will change the pointing to M45 for the next part of the question.
(OT2.2) Your Summary Answersheet shows telescopic field of M45. In the image, seven (7) brightest 15
stars of the cluster are replaced by + sign. Compare the image with the field you see in the
telescope and number the + marks from 1 to 7 in the order of decreasing brightness (brightest
is 1 and faintest is 7) of the corresponding stars.
(OT3) The examiner will give you a moon filter, an eyepiece with a cross-wire and a stopwatch. Point the 20
telescope towards the Moon. Attach the filter to the telescope. On the surface of the Moon, you will see
several seas (maria) which are nearly circular in shape. Estimate the diameter of Mare Serenitatis, MSr ,
labelled as 1 in the figure below, as a fraction of the lunar diameter, Moon , by measuring the telescope
drift times, Moon and MSr , for the Moon and the mare, respectively.
Group Examination
Page 1 of 2
(G1) A spacecraft of mass m and velocity approaches a massive planet of mass M and orbital velocity , as
measured by an inertial observer. We consider a special case, where the incoming trajectory of the
spacecraft is designed in a way such that velocity vector of the planet does not change direction due to
the gravitational boost given to the spacecraft. In this case, the gravitational boost to the velocity the
spacecraft can be estimated using conservation laws by measuring asymptotic velocity of the spacecraft
before and after the interaction and angle of approach of the spacecraft.
)
(G1.1) What will be the final velocity ( f of the spacecraft, if
and
are exactly anti-parallel (see 3
Figure 1).
(G1.2) Simplify the expression for the case where << . 1
(G1.6) On the day of the encounter, around what standard time (std ) had the Jupiter transited the 6
meridian in the sky of Bhubaneswar (20.27 N; 85.84 E; UT + 05: 30)?
(G1.7) Speed of the spacecraft (in km s 1) as measured by the same observer on some dates before the 12
encounter and some dates after the encounter are given below. Here day n is the date of
encounter. Use these data to find the orbital speed of Jupiter () on the date of encounter and
angle .
(G1) A spacecraft of mass m and velocity ~v approaches a massive planet of mass M and orbital velocity
~u, as measured by an inertial observer. We consider a special case, where the incoming trajectory
of the spacecraft is designed in a way such that velocity vector of the planet does not change
direction due to the gravitational boost given to the spacecraft. In this case, the amount of
gravitational boost to the velocity the spacecraft can be roughly estimated using conservation
laws by measuring asymptotic velocity of the spacecraft before and after the interaction and angle
of approach of the spacecraft.
Planet Planet
u u
v
spacecraft
v
spacecraft
Figure 1 Figure 2 Figure 3
(G1.1) What will be the final velocity (v~f ) of the spacecraft, if ~v and ~u are exactly anti-parallel 3
(see Figure 1).
Solution:
Let v~f and u~f be the final velocity of the spacecraft and the planet respectively. As
the planet For anti-parallel case, using conservation of linear momentum,
M~u + m~v = M u~f + mv~f
M u mv = M uf + mvf
m
uf = u (vf + v)
M
Now, using conservation of energy,
M u2 + mv 2 = M u2f + mvf2
m 2 m 2 m
u2 + v = u (vf + v) + v2
M M M f
m 2 m m m m 2
u2 + v = u2 + (vf + v)2 2u (vf + v) + vf
M
M M
M
M
m
0= (vf + v)2 2u(vf + v) + (vf2 v 2 )
M
m
0= (vf + v) f + v) 2u
(v (v
f + v) + f + v)(vf v)
(v
M
m
0= (vf + v) 2u + vf v
M
m m
vf 1 + = 2u + 1 v
M M
m
2u + 1 M v
vf = m
1+ M
Alternative solution in COM frame
Solution:
If m M ,
vf 2u + v
(G1.3) If angle between ~v and ~u is and m M (see Figure 2), use results above to write 3
expression for the magnitude of final velocity (vf ).
Solution:
As velocity vector of the planet is not changing direction, there is no momentum
transfer in direction perpendicular to ~u. We will resolve ~v and v~f into components
parallel and perpendicular to ~u.
vx = v cos vy = v sin
vfx = 2u + v cos vfy = v sin
vf2 = vf2x + vf2y = (2u + v cos )2 + (v sin )2
= 4u2 + 4uv cos + v 2 cos2 + v 2 sin2
= 4u2 + 4uv cos + v 2
p
vf = 4u2 + v 2 + 4uv cos
(G1.4) Table on the last page gives data of Voyager-2 spacecraft for a few months in the year 8
1979 as it passed close to Jupiter. Assume that the observer is located at the centre of
the Sun. The distance from the observer is given in AU and is heliocentric ecliptic
longitude in degrees. Assume all objects to be in the ecliptic plane. Assume the orbit
of the Earth to be circular. Plot appropriate column against the date of observation to
find the date at which the spacecraft was closest to the Jupiter, and label the graph as
G1.4.
Solution:
From the graph, it can be inferred that the encounter with Jupiter occured on
10th July (day 191) and its distance from the Sun on that day is 5.33121 AU
(G1.5) Find the Earth-Jupiter distance, (dEJ ) on the day of the encounter. 4
Group Examination
Page 3 of 8
Solution:
The day number of Vernal Equinox is 80. Thus, ecliptic longitude of the Sun as seen
from the Earth on the day of encounter will be,
= (191 80) 360 /365.25 = 109.4045
Thus, the ecliptic longitude of the Earth as seen from the Sun on the day of encounter
will be,
= 180 + 109.4045 = 289.4045
Applying cosine rule,
q
dJ = d2 + d2J 2d dJ cos
= 12 + 5.33122 2 1 5.3312 cos(289.4045 137.5628 )
p
= 6.2308 AU
i.e. the Earth is 6.2308 AU from Jupiter on that day.
(G1.6) On the day of the encounter, around what standard time (tstd ) had the Jupiter transited 6
the meridian in the sky of Bhubaneswar (20.27 N; 85.84 E; UT + 05:30)?
Solution:
Thus, the angle of eastern elongation for Jupiter (]SEJ) on that day would be,
1 5.3312 sin(289.4045 137.5628 )
= sin
6.2308
= 23.8146
It would rise 95 minutes after the Sun rise, i.e. around 7:35am. It would transit the
meridian after around 6 hours i.e. around 13:35 local time or 13:22 IST .
For more precise answer, R.A. of Jupiter on the day of encounter is approximately,
Jgeocentric = 109.4045 + 23.8146 = 133.2191
tan J = tan Jgeocentric cos
= tan 133.2191 cos 23 260
J = 135.68 = 9h 3m
Thus, it will culminate at that sidereal time. On that day, sidereal time at noon
is 07:24 (111 days from V.E. times 4 minutes). Thus, it will culminate 1 hour 39
minutes after the local noon i.e. at 13:39 local time or at about 13:26 IST .
Group Examination
Page 4 of 8
To V.E. S
S E
1
A
J J D
u
B 2
Figure 4 C
Figure 5
(G1.7) Speed of the spacecraft (in km s1 ) as measured by the same observer on some dates 12
before the encounter and some dates after the encounter are given below. Here day n is
the date of encounter. Use these data to find the orbital speed of Jupiter (u) on the date
of encounter and angle .
date n 45 n 35 n 25 n 15 n5 n
vtot 10.1408 10.0187 9.9078 9.8389 10.2516 25.5150
date n+5 n + 15 n + 25 n + 35 n + 45
vtot 21.8636 21.7022 21.5580 21.3812 21.2365
Solution:
In Figure 5, path of Voyager-2 is shown as A-B-C. The Sun is shown as S and the
Jupiter is shown as J. From the data we note that r is increasing continuously. The
same should be reflected in the diagram. For practical purpose, J and B are the same
points. The direction of velocity vector of Jupiter is given by JD. In the figure,
]ASB = 1 ]ASB = 2
]ASC = ]ABD =
]ABC = 1 ]DBC = ]ABC ]ABD = 1
]SAB = 1 ]SCB = 2
Now the lines originating from the Sun indicate radial direction on the respective
dates. Let us take speed of the spacecraft sufficiently far from the day 190, to avoid
any influence of Jupiter in initial and final velocity estimation. We can choose dates
35 days on either side of July 10 i.e. June 5 and August 14.
1 = 137.5628 136.0736 = 1.4892
p
l(AB) = l(SA)2 + l(SB)2 2 l(SA) l(SB) cos 1
p
= 5.174872 + 5.331212 2 5.17487 5.33121 cos 1.4892
= 0.207 55 au
5.33121 sin 1.4892
l(SB) sin 1
1 = sin1 = sin1
l(AB) 0.20755
Group Examination
Page 5 of 8
= sin1 (0.66755)
1 = 41.8783 or 138.1217
2 = 141.2007 137.5628 = 3.6379
p
l(BC) = l(SC)2 + l(SB)2 2 l(SC) l(SB) cos 2
p
= 5.450852 + 5.331212 2 5.45085 5.33121 cos 3.6379
= 0.362 53 au
1 l(SB) sin 2 1 5.33121 sin 3.6379
2 = sin = sin
l(BC) 0.36253
= sin1 (0.66755)
2 = 68.9199 or 111.0801
from the figure, 1 should be obtuse and 2 may be acute. In SABC
= 2 1 = 141.2007 136.0736
= 5.1271
1 = 360 1 2
= 360 5.1271 138.1217 68.9199
1 = 147.8313
In 4SBC, we notice
]SBC = 180 2 2
= 180 68.9199 3.6379
= 107.4422
vy v sin
tan ]DBC = f =
v xf v cos + 2u
sin
tan(1 ) =
cos + 2 uv
2u sin
= cos
v tan(1 )
We use this expression to find |~u|.
vf2 = 4u2 + v 2 + 4uv cos
vf2 4u2 4u
2
= 2
+1+ cos
v v v
2u 2
v 2
f 2u
= +1+2 cos
v v v
2
sin sin
= cos + 1 + 2 cos cos
tan(1 ) tan(1 )
2
sin 2 sin cos
2 2 sin cos
= + cos + 1 + 2 cos2
tan(1 ) tan( ) tan( )
1 1
sin2
= + 1 cos2
tan2 (1 )
sin2
+ sin2 = sin2 cot2 (1 ) + 1
= 2
tan (1 )
Group Examination
Page 6 of 8
sin2
=
sin2 (1 )
vf sin sin
= =
v sin(1 ) sin 1 cos cos 1 sin
v
= sin 1 cot cos 1
vf
sin 1
tan = v
vf + cos 1
sin 147.8313
= 10.0187
= 1.4088
21.3812 + cos 147.8313
= 180 54.6328 = 125.3672
vf2 = 4u2 + v 2 + 4uv cos
21.38122 = 4u2 + 10.01872 + 4u 10.0187 cos 125.3672
(457.1557 100.3743)
0 = u2 5.7990u
4
2
0 = u 5.7990u 89.1953
5.7990 + 5.79902 + 4 89.1953
u=
2
= 12.7789 km s1
Jupiters orbital velocity on the day of encounter is 12.779 km s1 and the angle
between the initial velocity of the spacecraft and Jupiters velocity vectors is 125 220 .
Solution:
The angle between ~r and ~u on the day of encounter will be,
= ]SBC (1 )
= 107.4422 147.8313 + 125.3672
= 84.9781
Now we use angular momentum conservation to estimate eccentricity. If up and rp
represent perihelion velocity and perihelion distance of Jupiter,
s
GM 1 + e
rp up = aJ (1 e)
aJ 1e
p
= GM aJ (1 e2 )
rp up = ru sin
r2 u2 sin2
1 e2 =
GM aJ
2
5.331212 1.496 1011 12.7789 103 sin2 84.9781
=
6.6741 1011 1.9891 1030 5.202 60
= 0.99761
e = 1 0.99761 = 0.0489
Group Examination
Page 7 of 8
Solution:
To estimate longitude of perihelion, one should estimate true anomaly of Jupiter on
that day.
a(1 e2 )
r=
1 + e cos
a(1 e2 ) 5.20260 0.99761
0.0489 cos = 1= 1
r 5.33121
= 0.02646
= 122.754
Thus, the longitude of perihelion of Jupiter is,
p = J
= 137.5628 122.754
p = 14.809
Group Examination
Page 8 of 8