Neutral Hydrogen

A little more than a year ago, while examining a newly made ALFALFA grid, Riccardo noticed a rather bright source in the constellation of Leo, moving away from us with a velocity of only 264 km/sec. It was not (or only barely) resolved by ALFA’s 4 arcminute beam, and its HI line velocity width was very narrow, indicating either a fully face-on and/or a very low mass object. In fact, it met the criteria of a “ultra compact high velocity cloud”, the targets for which Betsey Adams has been hunting. With such a low redshift, it was not clear whether the signal arose from a cloud in the Milky Way’s halo or a previously unidentified, tiny nearby galaxy. Quick checks of the public Sloan Digital Sky Survey (SDSS) and Digital Sky Survey (DSS) images showed no obvious associated stellar population but a suggestion of some faint, blue-ish emission. Could this really be a very faint, very small Milky Way neighbor, a bona fide (almost) optically-invisible (dark) galaxy? We needed to determine its distance and to look for evidence of rotation (which would suggest the presence of dark matter), so the quest to obtain the required additional observations began.

Optical image of the starlight in Leo P.

Optical image of Leo P showing its starlight

Being able to respond quickly to potentially exciting discoveries is one of the reasons ALFALFA is a team effort. So right away, we contacted ALFALFA team members Kathy Rhode and John Salzer at Indiana University, because IU has access to good imaging instruments on the WIYN (Wisconsin-Indiana-Yale-NOAO) 3.5 meter telescope in Arizona, and John Cannon at Macalester College who has been undertaking the SHIELD (Survey of HI in Extremely Low mass Dwarfs) program with the VLA. We made a special plea to the director of the VLA for “director’s discretionary time” to take a quick peek at the HI source to look for rotation. Knowing us not to ask without good reason, the director approved our last minute request. Although we were sure of the reality of the signal, the (awesome) Undergraduate ALFALFA team, during one of the ALFALFA followup runs took a spectrum centered on the optical object to confirm the position and radio characteristics of the ALFALFA signal. Kathy and John S. were able to take some quick images during an already-scheduled observing run. Indeed there were stars, and even more importantly, not very many! And a single HII region, proving that star formation is taking place. The VLA observations were made a few months later and the map that John C. and his student Elijah Bernstein-Cooper (read Elijah’s comments in an earlier post) resolved and localized the HI gas and confirmed that the object is rotating. Rotation signifies the presence of a significant amount of dark matter proving its extragalactic nature. A truly tiny object, Leo P contains only a few hundred thousand stars, in contrast to the Milky Way’s tens of billions, but Kathy and John were able to tease out an H-R diagram of the stars, yielding a distance of about 1.75 Mpc (or 5 million light years). So, while Leo P meets Betsey’s criteria to be an ultra compact high velocity cloud, it is also a bona fide galaxy, discovered because of its hydrogen gas, not its starlight. In fact, it contains more mass in gas than in stars. Most recent spectroscopic observations made by another ALFALFA team member Evan Skillman of the University of Minnesota confirm its pristine nature as an object that has undergone very little enrichment in heavy elements due to nucleosynthesis in stars, earning it the designation “P” for “pristine”. We believe that Leo P has managed to retain its gas without forming stars because, in contrast to most dwarf galaxies which reside near large ones, it lives virtually isolated in the local universe, just outside the Local Group.

Leo P is the first example of the class of gas-bearing tiny galaxy for which ALFALFA was specifically designed to look. Betsey’s thesis has already produced a catalog of similar “dark galaxy” candidates even though the survey data processing is not yet complete. As in the case of Leo P, we are pursuing the required detailed observations of the very best candidates (see her post on her March 2013 observing run at WIYN with its new pODI camera). The ALFALFA hunt for (almost) dark galaxies continues, but now we have shown that they do exist and that we can find them!

As Betsey has mentioned, ALFALFA is a great way to find hydrogen-rich galaxies regardless of their stellar content. With ALFALFA, we can find lots of galaxies, including a lot that have very few stars and are therefore incredibly dim. This is just another way to learn about the incredible diversity and structure in the Universe around us, and when we’re done we’ll have a big pile of about 30,000 galaxies to work with. There’s a lot that we can learn about the Universe in this way, but I’d like to tell you about just one: the HI Mass Function. We call it the HIMF for short.


I talked in an earlier post about what the 21 cm line is but you might be wondering why we are interested in observing this line. Why spend so much time on a survey focused on detecting this transition line?

One of the main reasons is that it allows a different method of selecting galaxies. To date, most large surveys have been optical (or near infrared) surveys. This means that galaxies are identified based on the stars contained within them that shine at optical wavelengths. With ALFALFA, we are instead looking for galaxies based on their 21 cm line emission, which is a transition line of neutral hydrogen. So we are identifying galaxies based on their gas, rather than stellar, content. As a result, ALFALFA is able to easily find galaxies that don’t have a lot of stars (low surface brightness galaxies) but that have lots of gas; these are the types of galaxies that are often overlooked by optical surveys. Of course, ALFALFA misses galaxies that have no gas and lots of stars (such as ellipticals) that are detected in optical surveys, so you want surveys at all different wavelengths.

Observing the 21 cm line offers a few other advantages that are worth mentioning. If you know the total flux of radiation from the 21 cm line at your telescope and the distance to a galaxy, you can compute exactly the mass of neutral hydrogen in that galaxy. This gas mass can then be used to compare properties across different galaxies and look for correlations. One of the other major advantages is that due to how radio telescopes receive radiation, when observing a source you obtain spectral information naturally. By spectral information, I am referring to the fact that you record the flux as a function of frequency (or velocity from the Doppler Effect). You can see this below with data from an example galaxy – a classic double-horned HI (neutral hydrogen) line profile. The total flux of the galaxy comes from calculating the area underneath the peak. There’s also other information available here, though. The center of the profile gives the redshift, or recessional velocity, of the galaxy which can be used to estimate a distance. The width of the profile corresponds to the rotational velocity of the galaxy (with an inclination factor). The beauty of 21 cm observations is that we record all this data at once. If you wanted to do the same optically, you would first need to find a galaxy in a photometric survey (essentially, take a picture of it with a telescope) and then target it a second time for spectroscopic follow-up to find a redshift for the galaxy. Of course, optical surveys have their advantages too.

Example HI Line Profile

Example HI Line Profile

The ALFALFA Survey is trying to detect galaxies through emission at the ’21cm line’. In fact, you’ll see mention of the 21cm line throughout extragalactic radio astronomy as it is one the most common observations. This might lead you to wonder: What is the 21cm line? And why do we want to observe it? For now, I’ll address the first of those questions; you’ll learn all about why the 21cm line is interesting in a later post.

Alternately, this post could be titled: “A Quick and Dirty Introduction to the Hydrogen Atom”.  The ’21cm line’ is a specific transition in the hydrogen atom, so called because the energy of the transition corresponds to a photon with a wavelength of 21 cm.  Physicists love hydrogen because it’s the simplest atom – one proton and one electron – meaning it can be solved exactly.  Astronomers love hydrogen because it fills the universe and observing it in different states offers lots of information about conditions found in astronomical situations.

Most people are familiar with the simple Bohr model of the atom where electrons orbit the nucleus in quantized orbits.  Of course, like many physics explanations, the Bohr model isn’t accurate but it does provide a convenient picture for much of the physics so we continue to use it.  The transitions in atoms that most people are familiar with are electronic transitions.  The orbits of electrons in an atom have quantized energy states and there is a specific energy associated with a move from one orbit to another.  For hydrogen, these electronic transitions have energies that release photons in the ultraviolet, optical and infrared parts of spectra.  The energies associated with the 21 cm line are much, much lower.  In order to understand where the 21cm line comes from we have to look further at the structure of hydrogen. The first set of corrections made are termed “fine structure” and account for the fact that the base calculation ignores relativistic effects.  These aren’t the corrections that cause the 21 cm line, though.  For that we need to go to a second set of corrections (which have an even smaller change in energy levels) – the “hyperfine structure” of hydrogen.

An illustration of the Bohr model of the atom

An illustration of the Bohr model of the atom

The hyperfine structure refers to the weak interaction between the spins of the proton and electron.  Spin is an extremely quantum idea and hard to describe with a classical analogy.  The best way to think of it is as particles spinning on their own axis – it behaves as an angular momentum term in many ways.  However, it does have some odd properties, which we will ignore for now.  If you keep the simple picture of a charged particle spinning in mind, you can imagine a slight magnetic field resulting as moving charges cause magnetic fields.  Similarly, the proton creates a slight magnetic field.  The interaction of these two magnetic fields results in the 21 cm transition.

The electron and proton can both be thought of as magnetic dipoles generating magnetic fields.  There is then an energy difference between when the dipole moments of the electron and proton are aligned or anti-aligned. The state where the two spins (or magnetic moments) are anti-aligned is a lower energy state; when the atom transitions from a “spin-up” electron to a “spin-down” electron, a photon with a wavelength of 21 cm (corresponding to the energy difference between the levels) is released.

Illustration of the 21 cm transition

Illustration of the 21 cm transition

I’m going to try ending posts with the math so that it’s available if you’re interested and easy to skip if you’re not – hopefully everything I write will be clear without any equations.

The energy of a photon is related to its frequency by E=h\nu where \nu is the frequency and h is Planck’s constant.

If you would like to know the wavelength of this photon, you can find it from the speed of light: c=\nu \lambda where \lambda is the wavelength.

If you would like to calculate the energy of the hyperfine transition a back-of-the-envelope approximation is:

\Delta E \sim \frac{2 \mu_e \mu_p}{a_0^3}

where a_0 is the Bohr radius of an atom and \mu_{e,p} are the magnetic moments of the electron and proton, respectively.

Substituting in for the magnetic moments this becomes,

\Delta E \sim \frac{1}{2} g_p \frac{m_e}{m_p} \alpha^2 \frac{e^2}{a_0}

where g_p is the degeneracy of the proton ( ~2.79), m_{e,p} are the masses of the electron and proton, respectively, \alpha is the fine structure constant (~1/137), and e is the charge of the electron.

This approximation is very close to the true calculation using quantum mechanics:

\Delta E = \frac{8}{3} g_p \frac{m_e}{m_p} \alpha^2 \frac{e^2}{a_0}