The Hydrogen 21-cm Line and Its Applications To Radio Astrophysics
The Hydrogen 21-cm Line and Its Applications To Radio Astrophysics
The Hydrogen 21-cm Line and Its Applications To Radio Astrophysics
MIT Undergraduate
(Dated: March 9, 2008)
We explore some applications of the 21-cm neutral hydrogen line in astrophysics. From observa-
tions in a narrow band around 21 cm we nd estimates for the peak temperature of the chromosphere
of the sun, as well as the angular extent of a portion of the galactic plane around l = 30
. We de-
termine the rotation curve of the Milky Way galaxy (for r < R) and explore its spiral structure
for galactic longitudes 30
< l < 80
=
p
B
e
. When a hydrogen
atom transitions from the excited state into the ground
state, a photon is emitted carrying away the energy dif-
ference. This photon has a wavelength of 21 cm.
The spin-ip transition, as it is termed, is a highly for-
bidden process, with a mean-life of approximately 10
7
years. The Heisenberg Uncertainty relationship then
leads us to expect a very sharp emission line with small
E energy dispersion, or line width, in frequency. This
obs
=
source
1
1 +
(1)
we can obtain the recessional velocity, v
rec
, by solving
the equation for .
3. EQUIPMENT AND TEMPERATURE
CALIBRATION
For our experiment we have access to the Small Ra-
dio Telescope on the roof of Building 26 developed by
Haystack Laboratories. A motor and feedback loop con-
trols the pointing of the telescope. Radio signals received
are fourier transformed and sorted into frequency bins
in a multi-channel accumulated spectrum of counts con-
verted to power (in K) against frequency.
A temperature calibration is performed internally by
the electronics of the system. A noise diode emitting
black body radiation at 115 K is activated at the vertex
of the parabolic dish. The eect is a signal that lls up
the entire antenna and is read by the feed horn. Subse-
quent signals are compared against this calibration and
converted into an eective temperature reading that is
proportional to the power received. The dish has a half-
power beamwidth of approximately 7.5
[5]. We point
our receiver to a patch of the cold sky (elevation > 50
)
in order to perform our calibration and obtain a baseline
system noise temperature, T
sys
.
4. TELESCOPE POINTING (25-POINT SCAN)
We use the sun as a means of correcting pointing errors
on the telescope. An Npoint scan of the sun creates a 5
5 grid around a central point and takes a temperature
measurement at each point in the grid. 25 temperatures
2
are recorded in all and can be used to create a tempera-
ture contour plot of the sun, see Figure 1. The telescope
recenters on the point of maximal temperature; the o-
sets are recorded. By using this method, we nd a general
repeatable pointing error of approximately 2
.
FIG. 1: A 25-point scan of the sky surrounding the sun. Contour
curves are shown.
5. BEAM WIDTH DETERMINATION AND
PEAK TEMPERATURE OF THE SUN
We point the telescope approximately 45 minutes
ahead of the sun. By turning o the tracking mechanism,
we allow the sun to drift through our antenna beam, tak-
ing integrated readings approximately 10 seconds apart.
We then sum over the power recorded in each bin in or-
der to obtain an average power per bin as a function of
time. This manipulation, as well as the integration time,
is factored into the error on each data point. Our results
from the scan is given in Figure 2.
The Fraunhofer Diration pattern, intensity as a func-
tion of angle , is described in the following relation in-
volving the Bessel function J
1
[4]:
I() = I
0
_
2J
1
(kasin)
kasin
_
2
(2)
The theoretical rst minimum occurs at sin =
1.22/D. Here, = 21 cm and D = 8.2 ft. We nd
it to be 5.87
)
FIG. 2: Drift scan over the sun. The Bessel diraction pattern
is smeared out by the nite angular extent of the sun and the
resulting curve is a reasonable gaussian t. We use 2 on the
gaussian distribution as an approximation for the rst minimum of
the Bessel function.
FIG. 3: Drift scan over the galactic plane around l = 30
. A
similar procedure is used to extract the rst minimum.
and the reduction of a circularly polarized radiation to a
linearly polarized receiver. We nd a T
sun
of 3.510
4
K.
T
sun
= 2
_
7.5
2
.5
2
_
77.5K = (3.5 .3) 10
4
K (3)
Here, we use the half-power beamwidth 7.5
instead of
the full beamwidth in order to account for the non-
uniform receiving power of the antenna. Since this re-
sult is now far above the known surface temperature
of the sun ( 6000K), we are led to believe that the
bulk of the 21 cm radiation from the sun has its ori-
gins in its chromosphere as opposed to its surface. In
another trial, we allow a portion of the galactic plane
around l = 30
. We use a simpli-
ed model and treat both solar and galactic emissions as
3
step functions to be convolved with the diracting beam
pattern. We estimate the bulk emission from this part of
the galaxy to be coming from an angular region approx-
imately 8.7 5.48 0.5 = (2.7 0.2)
across.
6. EXPLORING GALACTIC STRUCTURE AND
DYNAMICS
Recalling that the sharpness of the 21-cm hydrogen
line allows for very precise doppler spectroscopy, we make
use of the abundance of neutral hydrogen in the Milky
Way to derive a galactic rotation curve (a v(r) for every
r). We model the galaxy as a dierentially rotating thin
disk with stellar and interstellar material in circular or-
bits around a shared origin. This symmetry assumption
allows the average velocity to be independent of polar
angle ,
v
r
, r > r
0
(4)
This is known as a Keplerian velocity curve, which is
obeyed by the planets in our solar system. For the Milky
Way galaxy, however, this model is unrealistic, as we well
know the mass is not concentrated at the galactic center
but spread out over many kpc. From observational evi-
dence, we know the Milky Way to consist of two visible
parts: a spherical central bulge with radius of approx-
imately r = 1.5 kpc with a constant mass density
B
[3], and an outer extent with luminosity falling o ex-
ponentially with r. Considering rst the central bulge,
integrating the constant density we nd a total enclosed
mass M(r) proportional to r
3
and a v(r) that increases
linearly in r.
Treating the outer regions is only slightly more in-
volved. We begin with the empirically determined sur-
face luminosity density curve for the galactic disk, (r) in
units of [L
/kpc
2
] [2]. Using two-dimensional geometry
now, we assign an average mass-to-light ratio, 0.6,
to convert to surface mass density in units of [M
/kpc
2
].
Finally, we integrate over r to nd total enclosed mass
at a given radius:
M(r) =
_
r
0
[(r)](2r)dr =
_
r
0
[
0
e
r/b
](2r)dr
M(r) =
0
_
b
2
_
1 e
r/b
_
rbe
r/b
_
(5)
Here, b is a characteristic radius, taken to be about 3 kpc,
and
0
is the maximum surface mass density. From there
we can quickly derive a v(r) using the relation v(r) =
_
GM(r)/r. We nd that the limiting behavior for v(r)
at large r is also Keplerian (goes as 1/
r). Since in
reality, we will be observing clumps of gravitating matter
as opposed to a single point mass, the v(r) from here on
will represent an average velocity.
6.2. Experimental Methods
We point our telescope to various galactic longitudes l
and collect doppler shifted 21-cm spectra. However, there
is no way to separate the accumulated received radiation
by source distance. In order to derive a rotation curve, we
will need to be able to associate each measured velocity
with a non-ambiguous radius. We do this by considering
only the maximum red/blue shift in each direction. We
propose that within 0 < l < 90, the maximum recessional
velocity corresponds to the matter traveling directly away
from us in a circle of radius R that is tangent to our line
of sight, assuming a relatively at rotation curve, due to
the maximization of the dot product v s, where s is the
line of sight. This is best understood with a diagram
(Figure 4).
FIG. 4: The maximum recessional velocity along a line of sight
corresponds to matter moving in a circle with radius R around the
galactic center.
The only visible portion of the galactic plane suitable
for this analysis has 0 < l < 90. We sweep this range to
nd maximal v
rec
in each direction and use trigonometry
to nd its associated radius R = R
sin(l).
4
6.3. Calibration
We use a frequency generator to conduct a 3-point cal-
ibration of our output. This procedure converts each bin
in our accumulated spectrum into a calibrated frequency
value. We obtain calibration errors through a linear t.
FIG. 5: Calibration of spectrum. Frequency vs. bin number.
6.4. Corrections - VLSR and Dispersion
Two key corrections were made to our raw data in
order to obtain the velocity prole. First, in order to ac-
count for the velocity dispersion due to relative motion
of gas clusters in a given region of the galaxy, we point
our telescopes exactly anti-center, at l = 180
. Since we
expect no net red- or blueshift in this direction (average
motion perpendicular to our line of sight), the width of
the spectrum can be used to determine a representative
velocity disperson. See Figure 6. We subtract this ex-
pected dispersion from the maximum redshift along our
line of sight to nd the maximum average recessional ve-
locity.
FIG. 6: 21cm power spectrum obtained from observation along l =
180. We nd the velocity dispersion to be approximately 20km/s.
A second correction of importance removes the eect
of the receivers motion relative to the galactic center.
This includes earths 30 km/s orbit around the sun and
the suns 20 km/s motion in the direction of RA 18h
Dec 30
and lattitude
0
.
Comparison to expected velocity functions.
6.6. Discussion
We see clearly non-Keplerian behavior in our rotation
curve. For small R, the rotational velocity rises linearly
in r as predicted by the bulge model, however for large
R the observed behavior deviates signicantly from the
expected Keplerian 1/
r
v(r) v
LSR
_
(7)
From r, we can now nd the angle :
(r) = G(r) =
_
l + sin
1
_
r
r
sin(l)
__
(8)
FIG. 10: Plot of the sources of maximum intensity 21 cm emission
between 30 < l < 80. We see possible evidence of two distinct
density waves (spiral arms) of high H1 concentration.
Performing this analysis for each redshifted peak along
30 < l < 80, we obtain the following polar plot showing
6
the locations of highest H1 density in the viewable por-
tion of the galactic plane (Figure 10). There is evidence
of some spiral structure as indicated by two possible arms
of high density H1 gas, whose locations are approximated
by superimposed curves to aid viewing. The outer curve
is likely the Orion arm, of which the solar system is a
part, and the inner may be a short portion of the Saggi-
tarius Arm.
7. ERRORS
The empirically determined parameters used in this
experiment, as a whole, were a surprisingly large contrib-
utor to general error. The orbital radius of the sun, R
,
for example, is only known to approximately 6-10% ac-
curacy, and was the largest contributor to uncertainty in
r in our galactic rotation curve. Similarly, v
LSR
and v
and were
corrected occasionally with osets to eliminate its sys-
tematic eects.
The need for a radio-quiet environment for high-
accuracy radio astrophysics became apparent after only a
few trials. Aside from noise within the electronics (tem-
perature calibrations determine T
sys
), radio towers, cell
phones, even satellites distorted our spectra immensely.
Although mostly contributing a noise oor, as we moved
to lower galactic longitudes and therefore nearer to the
horizon, not only were there occasionally loud signals
that blotted out the entire spectrum, but the noise oor
also shifted and began to take on some slope. This
made determination of maximum redshift dicult and
although we attempted to correct for this eect by tak-
ing a noise background at a similar elevation away from
the galactic plane, contributed signicantly to our error
in the lower longitudes.
Non-validity ranges in some of our assumptions is a
matter of greater concern. These errors contribute sys-
tematically to our results and may result in signicant
corrections and shifts. The assumption that v(r) is
mostly at and therefore the highest recessional velocity
corresponds to the radius of closest approach is obviously
called into question by our velocity curve itself. Wrong
attributions of maximal velocity would result in an over-
estimate of v(r) at small r. In addition, the assumption
of perfectly circular rotation around the galactic center
is refuted by evidence of spiral arms of high mass den-
sity with radial components. To some extent, this will
aect our velocity curve, but to a greater extent, it will
change the locations of the maxima in our plot of density
maxima and warp our spiral arms.
We would like to also point out certain limitations
to our method in estimating velocity dispersion. Since
stronger gravitational interactions would likely result in
a wider distribution of velocities, regions of lower stellar
density would not necessarily suer the same dispersive
forces as regions of higher stellar density. This would
result in an underestimate of certain velocities and an
overestimate of others.
8. CONCLUSIONS
In the earlier portions of our lab, we determined the
chromospheric temperature of the sun to be approxi-
mately (3.5 .3) 10
4
K, and the full beamwidth of
the telescope to correspond to a the circular diraction
theorys prediction of 1.22/D. Additionally, we mea-
sured the thickness of the peak 21cm emitting hydrogen
in the region of galactic longitude 30
and found it to be
approximately (2.7 0.2)
in angular extent.
In the latter parts of our experiment we went into some
detail investigating the dynamics and structure of the
Milky Way galaxy. We derived a rotation curve from
doppler spectroscopy performed on 21-cm H1 emissions
and found it to be non-Keplerian in nature, indicating
some misalignment between the expected gravitational
eects due to visible matter and the actual dynamics of
our galaxy. We oered the existence of a dark matter
halo as a possible explanation of this phenomenon, sup-
ported by our experimental results. Additional analysis
of the 21-cm spectra in various directions in the galactic
plane indicated the presence of 2 possible spiral arms.
Due to the large uncertainties in this particular analysis
the results are interesting but inconclusive and denitely
warrants further investigation. In the end, we believe we
have demonstrated the 21-cm neutral hydrogen line to be
a powerful tool for exploring our galaxy and beyond.
[1] Drees, M. Dark Matter. Technical University, Munich.
[2003]
[2] Mihos, Chris. Rotation Curves. Case Western Reserve
University. http://burro.astr.cwru.edu/ [2001]
[3] University of Oregon Physics.
http://physics.uoregon.edu/ jimbrau/ [2001]
[4] Diraction from a Circular Aperture. Paul Padley, Con-
nexions. [2005]
[5] Sewell, 21 cm Radio Astrophysics, 8.14 Course Reader,
[2007]
[6] van de Hulst, H.C., Muller, C.A., Oort, J.H. Spiral Struc-
ture... Derived from Hydrogen Emission at 21cm Wave-
length. Bulletin of Astronomical Institutes of Nether-
lands. Volume XII. Number 452. [1954]
7
[7] Shu, F.H. Our Galaxy: The Milky Way System. The
Physical Universe - Chapter 12. University Science Books.
[1982]
Acknowledgments
All non-linear ts were made with the MATLAB
scripts made available to us on the Junior Lab website:
http://web.mit.edu/8.13/www/jlmatlab.shtml.
We would like to thank the Junior Lab sta for their
assistance and Scott Sewell in particular for xing the
telescope repeatedly and in many ways so that this ex-
periment can be possible.