Spectral lines are narrow () emission or absorption features in the spectra of gaseous and ionized sources. Examples of radio spectral lines include recombination lines of ionized hydrogen and heavier atoms, rotational lines of polar molecules such as carbon monoxide (CO), and the cm hyperfine line of interstellar Hi.
Spectral-line emission and absorption are intrinsically quantum phenomena. Classical particles and waves are idealized concepts like infinitesimal points or perfectly straight lines in geometry; they don’t exist in the real world. Some things are nearly waves (e.g., radio waves) and others are nearly particles (e.g., electrons), but all share characteristics of both particles and waves. Unlike idealized waves, real radio waves do not have a continuum of possible energies. Instead, electromagnetic radiation is quantized into photons whose energy is proportional to frequency: . Unlike idealized particles, real particles of momentum are associated with waves whose De Broglie wavelength is . An electron’s stable orbit about the nucleus of an atom shares a property with standing waves: its circumference must equal an integer number of wavelengths. Planck’s constant erg s in these two equations is a quantum of action whose dimensions are (masslengthtime), the same as (energytime) or (angular momentum) or (lengthmomentum). Although has dimensions of energytime, physically acceptable solutions (the wave functions and their derivatives must be finite and continuous) to the time-independent Schrödinger equation exist only for discrete values of the total energy, so spectral lines have definite frequencies resulting from transitions between discrete energy states. A second quantum effect important to spectral lines, particularly at radio wavelengths where , is stimulated emission (Section 7.3.1). Fortunately, the fundamental characteristics of radio spectral lines from interstellar atoms and molecules can be derived from fairly simple applications of quantum mechanics and thermodynamics.
Spectral lines are powerful diagnostics of physical and chemical conditions in astronomical objects. Their rest frequencies identify the specific atoms and molecules involved, and their Doppler shifts measure radial velocities. These velocities yield the redshifts and Hubble distances of extragalactic sources, plus rotation curves and radial mass distributions for resolved galaxies. Collapse speeds, turbulent velocities, and thermal motions contribute to line broadening in Galactic sources. Temperatures, densities, and chemical compositions of Hii regions, dust-obscured dense molecular clouds, and diffuse interstellar gas are also constrained by spectral-line data. Radio spectral lines have some unique characteristics:
Their “natural” line widths are much smaller than Doppler-broadened line widths, so gas temperatures and very small changes in radial velocity can be measured.
Stimulated emission is important because . This causes line opacities to vary as and favors the formation of natural masers.
The ability of radio waves to penetrate dust in our Galaxy and in other galaxies allows the detection of line emission emerging from dusty molecular clouds, protostars, and molecular disks orbiting AGNs.
Frequency (inverse time) can be measured with much higher precision than wavelength (length), so very sensitive searches for small changes in the fundamental physical constants over cosmic timescales are possible.
Although the interstellar medium (ISM) of our Galaxy is dynamic, it tends toward a rough pressure equilibrium because mass motions with speeds up to the speed of sound try to reduce pressure gradients. Temperatures equilibrate more slowly, so there are wide ranges of temperature and particle number density consistent with a given pressure and the ideal gas law
Typical ISM pressures lie in the range – [jt11]. Radiative cooling by spectral-line emission depends strongly on temperature, so most of the ISM exists in several distinct phases having comparable pressures but quite different temperatures:
cold (10s of K) dense molecular clouds
cool ( K) neutral Hi gas
warm ( K) neutral Hi gas
warm ( K) ionized Hii gas
hot ( K) low-density ionized gas (in bubbles formed by expanding supernova remnants, for example)
All but the hottest phase are sources of radio spectral lines.
The semiclassical Bohr atom (Figure 7.1) contains a nucleus of protons and neutrons around which one or more electrons move in circular orbits. The nuclear mass is always much greater than the sum of the electron masses , so the nucleus is nearly at rest in the center-of-mass frame. The wave functions of the electrons have De Broglie wavelengths
where is the electron’s momentum and is its speed. Only those orbits whose circumferences equal an integer number of wavelengths correspond to standing waves and are permitted. Thus the Bohr radius of the th permitted electron orbit satisfies the quantization rule
where the number is called the principal quantum number. The requirement that
implies that the orbital angular momentum is an integer multiple of the reduced Planck’s constant . The relation between and is determined by the balance of Coulomb and centrifugal forces on electrons in circular orbits. For a hydrogen atom,
Numerically, the Bohr radius of a hydrogen atom whose electron is in the th electronic energy level is
The Bohr radius of a hydrogen atom in its ground electronic state () is only cm, but the diameter of a highly excited () radio-emitting hydrogen atom in the ISM can be remarkably large: m, which is bigger than most viruses!
The electron in a Bohr atom can fall from the level to , where and are any natural numbers (), by emitting a photon whose energy equals the energy difference between the initial and final levels. Such spectral lines are called recombination lines because formerly free electrons recombining with ions quickly cascade to the ground state by emitting such photons. Astronomers label each recombination line using the name of the element, the final level number , and successive letters in the Greek alphabet to denote the level change : for , for , for , etc. For example, the recombination line produced by the transition between the and levels of a hydrogen atom is called the H91 line.
The total electronic energy is the sum of the kinetic () and potential () energies of the electron in the th circular orbit:
The electronic energy change going from level to level is equal to the energy of the emitted photon:
so the photon frequency is
The factor in large parentheses is called the Rydberg constant , where the subscript refers to the limit of infinite nuclear mass :
The dimensions of are length, and the product is the Rydberg frequency:
Allowing for the relatively large but finite nuclear mass and repeating the analysis above in the atomic center-of-mass frame yields the same frequency formula with replaced by :
The hydrogen nucleus is a single proton of mass so and for a hydrogen atom is
Thus the frequency of the photon produced by the H109 transition ( to ) is
The mass of a neutron is about equal to the mass of a proton, so the He nucleus consisting of two protons and two neutrons has mass , the isotope of carbon with six protons and six neutrons has , and so on. Electrons recombining onto singly ionized atoms with any number of protons and electrons orbit in the potential produced by a net charge of one proton, so the recombination lines of heavier atoms are very similar to those of hydrogen, but at the slightly higher frequencies (Figure 7.2) given by Equation 7.12. For example, the primordial abundance of the rare helium isotope He is important because it depends on the photon/baryon ratio in the early universe. The abundance of He in Galactic Hii regions has been measured via radio recombination-line emission and indicates that baryons account for only a few percent of the critical density needed to close the universe.
The strongest radio recombination lines are produced by transitions with , so the approximation
yields simpler (but not extremely accurate) approximations for radio recombination line frequencies
and for the frequency separation between adjacent lines
Adjacent high- (low-) radio recombination lines have such small fractional frequency separations (Figure 7.3) that two or more transitions can often be observed simultaneously and averaged, to reduce the observing time needed to reach a given signal-to-noise ratio.
The H109 line was first detected by P. Mezger in 1965, despite (incorrect) theoretical predictions that pressure broadening would smear out the lines in frequency and make them undetectable. It is true that atomic collisions in the interstellar medium significantly disturb the energy levels of large atoms, but this disturbance is about the same for adjacent energy levels, so the differential disturbance that alters the line frequency is actually much smaller. His advice: “Don’t abandon an observation just because you have been told that it will fail.”
The spontaneous emission rate is the average rate at which an isolated atom emits photons. Rigorous quantum-mechanical calculations of spontaneous emission rates are complicated, but a fairly good classical approximation can be derived by noting that radio photons are emitted by atoms with and invoking the correspondence principle, Bohr’s hypothesis that systems with large quantum numbers behave almost classically. The time-averaged radiated power for classical transitions is given by Larmor’s formula (Equation 2.143) for an electric dipole with dipole moment :
The photon emission rate (s) equals the average power emitted by one atom divided by the energy of each photon. The spontaneous emission rate for transitions from level to level is denoted by :
in the limit (Equation 7.14). In that limit, also.
The atomic radius (Equation 7.6) is
Evaluating the constants yields the spontaneous emission rate of hydrogen atoms:
For example, the 5.0089 GHz H109 transition rate is .
The associated natural line width or intrinsic line width follows from the uncertainty principle: . Substituting for and for for each energy level involved in the transition and summing these two uncertainties yields
Natural broadening is negligibly small at the large that produce radio-frequency photons. Collisions of the emitting atoms cause collisional broadening, where the amount of collisional broadening is a small fraction of the collision rate when . Except for very large , collisional broadening is also small and the actual line profile (normalized intensity as a function of frequency) is primarily determined by Doppler shifts reflecting the radial velocities of the emitting atoms. The sign convention is for sources moving away from the observer. Radial velocities may be microscopic (from the thermal motions of individual atoms) or macroscopic (from large-scale turbulence, flows, or rotation). In the nonrelativistic limit , the Doppler equation (Equation 5.142) relating the observed frequency to the line rest frequency reduces to
so nonrelativistic radial velocities can be estimated from
The thermal component of the line profile from a recombination-line source in LTE is determined by the Maxwellian speed distribution (Equation B.49) of atoms with mass and temperature . The speed in any one coordinate of an isotropic distribution is of the total speed in three dimensions, so the Gaussian
is the normalized () radial velocity distribution. The normalized line profile (Figure 7.4) for thermal emission is
This is a Gaussian line profile. Its full width between half-maximum points (FWHM) is the solution of
For example, the FWHM of the H109 line ( GHz) in a quiescent (no macroscopic motions) Hii region with temperature K is
Notice that the thermal line width is much larger than the natural line width Hz.
Normalization (requiring ) implies that the value of at the line center () is
For a given integrated (over frequency) line strength, the line strength per unit frequency at any one frequency (e.g., at ) is inversely proportional to the line width . Integrated line strengths are frequently specified in the astronomically convenient units of , where .
The spontaneous emission coefficient is the average photon emission rate (s) for an “undisturbed” atom or molecule transitioning from an upper (U) to a lower (L) energy state. The spectral-line radiative transfer problem also involves the absorption coefficient and the stimulated emission coefficient (Figure 7.5). Einstein showed that both the absorption and stimulated emission coefficients can be calculated from the spontaneous emission coefficient.
Consider any two energy levels and of a quantum system such as a single atom or molecule. The photon emitted or absorbed during a transition between the upper and lower states will have energy
and contribute to a spectral line with rest frequency . The energy levels actually have small but finite widths, so the spectral line has some narrow line profile centered on and conventionally normalized such that . A system in the lower energy state may absorb a photon of frequency and transition to the upper state. The rate (s) for this process is proportional to the profile-weighted mean radiation energy density
of the surrounding radiation field, so the Einstein absorption coefficient is defined to make the product
equal the average rate (s) at which photons are absorbed by a single atomic or molecular system in its lower energy state.
Einstein realized that there must be a third process in addition to spontaneous emission and absorption. It is stimulated emission, in which a photon of energy stimulates the system in the upper energy state to emit a second photon with the same energy and direction. The rate for this process is also proportional to , so by analogy with Equation 7.40, the Einstein stimulated-emission coefficient is defined to make the product
equal the average rate (s) of stimulated photon emission by a single quantum system in its upper energy state. Beware that some authors use instead of in Equation 7.41 to define a stimulated emission coefficient that is times the given by Equation 7.41 [rl79].
Stimulated emission is sometimes called negative absorption. Negative absorption is not familiar in everyday life because it is much weaker in room-temperature objects at visible wavelengths, where , but negative absorption competes effectively with ordinary absorption at radio wavelengths where .
Consider a macroscopic collection of many atoms or molecules in full thermodynamic equilibrium (TE) with the surrounding radiation field. TE is a stationary state. If there are atoms or molecules per unit volume in the (upper, lower) energy states, then the average rate of photon creation by both spontaneous emission and stimulated emission must balance the average rate of photon destruction by absorption:
In TE, the ratio of to is fixed by the Boltzmann equation
where and are the numbers of distinct physical states having energies and , respectively. The quantities and are called the statistical weights of those energy states. Examples of statistical weights include the following:
Hydrogen atoms have , where is the principal quantum number. The number is the product of the 2 electron spin states and orbital angular momentum states in the th electronic energy level.
Rotating linear molecules (e.g., carbon monoxide, CO) have , where is the angular-momentum quantum number. For each , there are possible values of the -component of the angular momentum: .
Hydrogen atoms have two hyperfine energy levels whose difference yields the cm ( MHz) Hi line: and .
Solving Equation 7.42 for the profile-weighted mean energy density of blackbody radiation connects the properties of the quantum system (atom or molecule) to the radiation, just as Kirchhoff’s law (Equation 2.30) did for continuum radiation:
for matter and Equation 2.92 implies
for radiation. Inserting the Planck radiation law (Equation 2.86) for near gives
for all temperatures , so the equation
Equations 7.50 and 7.51 are called the equations of detailed balance. The two equations relate the three quantities , , and , so all three can be computed if only one (e.g., the spontaneous emission coefficient ) is known. Equations 7.50 and 7.51 also prove that cannot be zero; that is, stimulated emission must occur.
Note that these Equations 7.50 and 7.51 are valid for any microscopic physical system because they relate the coefficients , , and , characteristic of individual atoms or molecules for which the macroscopic statistical concepts of TE or LTE are meaningless. Although TE was assumed for their derivation, the dependences on temperature and frequency dropped out for a line at a single frequency . Thus Equations 7.50 and 7.51 also apply to all macroscopic systems, whether or not they are in TE or even LTE. (Recall the derivation of Kirchhoff’s law (Equation 2.30), which also made use of full TE but also relates the emission and absorption coefficients of any matter in LTE at temperature , independent of the actual ambient radiation field.)
The two Equations 7.50 and 7.51 for detailed balance relate the three Einstein coefficients and allow the spectral-line radiative transfer problem to be solved in terms of the spontaneous emission coefficient alone. The radiative transfer equation (2.27) is
where is the specific intensity, is the net fraction of photons absorbed (the difference between ordinary absorption and negative absorption) per unit length, and is the volume emission coefficient.
The ordinary opacity coefficient is the fraction of spectral brightness removed per unit length by absorption from the lower level to the upper level. At frequencies near the line center frequency , the photon energy is , the number of absorbers per unit volume is , the number of absorptions per unit area per unit time is , the fraction of photons absorbed per unit area per unit length is , and the photon energy loss per unit length at frequency is
Stimulated emission is best treated as negative absorption because, like ordinary absorption and unlike spontaneous emission, its strength is proportional to . The derivation of Equation 7.53 can be repeated to give
The spontaneous emission coefficient is the spectral brightness (power per unit frequency per steradian) added per unit volume by spontaneous transitions from the upper to lower energy levels. As above, the line photon energy is approximately , the number density in the upper energy level is , and the photon emission rate per unit volume is . These photons are emitted isotropically over , so
The ratio of the emission coefficient to the net absorption coefficient is
Equation 7.50 can be used to eliminate from this ratio also:
Finally, Equation 7.51 can be used to eliminate both and :
In LTE, Kirchhoff’s law independently implies
recovering the Boltzmann equation (Equation 7.43) for LTE (not just for full TE):
Equation 7.58 and the assumption of LTE allows the substitution of
to yield the net line opacity coefficient in LTE:
in terms of the spontaneous emission rate only; the stimulated emission coefficient and absorption coefficient have been eliminated.
in Equation 7.67 is the sum of two terms, the first representing ordinary absorption and the second accounting for the negative absorption of stimulated emission. In the Rayleigh–Jeans limit ,
Thus stimulated emission nearly cancels pure absorption and significantly reduces the net line opacity at radio frequencies. Because and , the product is independent of temperature. The brightness of an optically thin () radio emission line is proportional to the column density of emitting gas but can be nearly independent of the gas temperature. Thus the Hi line flux () of an optically thin galaxy is proportional to the total mass of neutral hydrogen in the galaxy but says nothing about its temperature.
Even if a macroscopic two-level system is not in LTE, its excitation temperature can be defined by
The excitation temperature is not a real temperature; it only measures the ratio of to . In a two-level system, the excitation temperature is determined by a balance between radiative and collisional excitations and de-excitations. If collisions cause excitations per unit volume per unit time from the lower level to the upper level and de-excitations per unit volume per unit time from the upper level to the lower level, then Equation 7.42 becomes
and detailed balance requires
where is the ambient radiation brightness temperature, in favor of . Equation 7.50 eliminates
where is the kinetic temperature of the gas. The numerator of Equation 7.73 becomes
and the denominator is
If the spontaneous emission rate is much larger than the collision rate, Equation 7.80 yields ; if the collision rate is much higher than the spontaneous emission rate, . For any and , lies between and .
If the upper energy level is overpopulated, that is,
then is actually negative,
is negative, and Equation 7.67 gives a negative net opacity coefficient . Negative net opacity implies brightness gain instead of loss; the intensity of a background source at frequency will be amplified. At radio wavelengths this phenomenon is called maser (an acronym for microwave amplification by stimulated emission of radiation) amplification. Astronomical masers are common at radio frequencies because and hence even in TE. These sources can have line brightness temperatures as high as K, which is much higher than the kinetic temperature of the masing gas. For a clear presentation covering the basics of astronomical masers, written by Reid & Moran, see ver88.
Our model for an astronomical maser starts with the radiative transfer equation for a two-level system. Assume for simplicity that so that Equation 7.50 implies and Equation 7.51 implies . Then Equation 7.57 simplifies to
Then at the line center frequency ,
The maser optical depth
is called the maser gain over the path of integration, and the maser amplifies the intensity of background radiation by the factor . In a laboratory maser, radiation is trapped in a high-Q resonant cavity to create effective path lengths up to times the cavity length. In an astronomical maser there is no cavity, so the radiation makes only a single pass and the physical path must be much longer () for significant gain to occur.
Maser emission quickly depopulates the upper energy level, so masers have to be “pumped” to emit continuously. Typically one or more higher energy levels absorb radiation from a pump source (e.g., infrared continuum from a star or an AGN), and radiative decays preferentially repopulate the upper energy level. This radiative pumping process produces no more than one maser photon per pump photon, so the pump energy required is proportional to the frequency of the pump photon. If the maser photon emission rate is limited by the pump luminosity, the maser is described as being saturated; if the pump power is more than adequate, the maser is unsaturated.
Strong, compact maser sources are powerful tools for high-resolution imaging and precision astrometry. They are being used to measure accurate trigonometric distances to individual stars in our Galaxy, the size and structure of our Galaxy, black-hole masses in AGNs, and distances to galaxies up to [rei14].
The most spectacular example is the circumnuclear 22 GHz megamaser disk surrounding the Seyfert nucleus of NGC 4258 [mor08] and illustrated in Figure 7.6. The observed maser spectrum has hundreds of narrow lines in three groups clustered around the systemic recession velocity of NGC 4258 and at velocities further redshifted and blueshifted by the rotation velocities in the disk.
The nearly edge-on maser disk was imaged with high spectral and angular resolution by the VLBA, and its rotation curve is shown in Figure 7.7. The systemic lines are concentrated within a region wide and centered on the line of sight to the supermassive black hole (SMBH) at the center of the disk. Photons from the central radio continuum source are amplified and beamed in our direction, so systemic maser clouds are visible only when they are within a small angle from the line of sight. For a constant , their yields the sloped straight line in Figure 7.7. This slope corresponds to the gravitational acceleration of the SMBH, and long-term monitoring the velocities of individual systemic-maser lines independently yields . The redshifted and blueshifted maser clouds are visible only near the tangent points of the disk near , which have the longest gain paths at nearly constant velocities that follow the perfectly Keplerian curves in Figure 7.7. Combining the acceleration, velocity, and angular-size measurements yields a unique solution for both the SMBH mass and the distance to NGC 4258.
Astronomical sources of radio recombination lines are often in local thermodynamic equilibrium (LTE). The spontaneous emission rate from atomic physics and the laws of radiative transfer for spectral lines can be combined to model sources in LTE. However, LTE does not apply to all recombination lines, so departures from LTE must be recognized and treated differently.
Equation 7.67 gives the absorption coefficient at the center frequency of the electronic transition of hydrogen in an Hii region in local thermodynamic equilibrium (LTE) at electron temperature :
and is the number density of atoms in the th electronic energy level. At radio frequencies, it can be assumed that , , and . Equation 7.23 gives the spontaneous emission rate
and Equation 7.37 parameterizes the normalized line profile
The number density of atoms in the th electronic energy level is given by the Saha equation, a generalization of the Boltzmann equation (for a derivation of the Saha equation, see rl79):
where is the ionization potential of the th energy level. For large , and the exponential factor can be ignored. Combining the results from Equations 7.87 through 7.92 yields the opacity coefficient at the line center frequency :
Some algebra reduces this to
Notice that the electronic energy level has dropped out; Equation 7.94 is valid for all radio recombination lines with . The optical depth at the line center frequency can be expressed in terms of the emission measure defined by Equation 4.57:
In astronomically convenient units the line center opacity is
Because in all known Hii regions, the brightness temperature contributed by a recombination emission line at its center frequency is
At frequencies high enough that the free–free continuum is also optically thin, the peak line-to-continuum ratio (which occurs at frequency ) in LTE is
where is the line FWHM expressed as a velocity and the typical He/H ion ratio is . The term in square brackets is necessary because He contributes to the free–free continuum emission but not to the hydrogen recombination line. The line-to-continuum ratio yields an estimate of the electron temperature that is independent of the emission measure so long as the frequency is high enough that the continuum optical depth is small.
Recombination lines can be used to find the electron temperatures of Hii regions in LTE. Solving Equation 7.98 explicitly for gives the useful formula
By mapping the recombination line-to-continuum ratios in a number of H transitions, Lockman and Brown determined the temperature distribution in the Orion Nebula (Figure 7.8), a nearby Hii region.
Differences between the rest and observed frequencies of radio recombination lines are attributed to Doppler shifts from nonzero radial velocities. With a simple rotational model for the disk of our Galaxy, astronomers can convert radial velocities to distances, albeit with some ambiguities, and map the approximate spatial distribution of Hii regions in our Galaxy (Figure 7.9). They roughly outline the major spiral arms.
A plot showing the observed electron temperatures of Galactic Hii regions (Figure 7.10) reveals that temperature increases with distance from the Galactic center.
The explanation for this trend is the observed decrease in metallicity (relative abundance of elements heavier than helium) with galactocentric distance. Power radiated by emission lines of “metals” is the principal cause of Hii region cooling.
Radio recombination line strengths are much less affected by dust extinction than optical lines (e.g., the H and H lines) are, so radio recombination lines are useful quantitative indicators of the ionization rates and hence star-formation rates in dusty starburst galaxies such as M82 (Figure 7.11).
A molecule is called polar if its permanent electric dipole moment (Equation 7.120) is not zero. Symmetric molecules (e.g., the diatomic hydrogen molecule H) have no permanent electric dipole moment, but most asymmetric molecules (e.g., the carbon monoxide molecule CO) do have asymmetric charge distributions and are polar. The electric dipole moments of polar molecules rotating with constant angular velocity appear to vary sinusoidally with that angular frequency, so polar molecules radiate at their rotation frequencies. The intensity of this radiation can be derived from the Larmor formula expressed in terms of dipole moments instead of charges and charge separations.
The permitted rotation rates and resulting line frequencies are determined by the quantization of angular momentum. The quantization rule for the permitted electron orbital radii
in the Bohr atom quantizes the orbital angular momentum in multiples of :
The rule that angular momentum is an integer multiple of is universal and applies to the angular momentum of a rotating molecule as well.
Consider a rigid diatomic molecule (Figure 7.12) whose two atoms have masses and and whose centers are separated by the equilibrium distance . The individual atomic distances and from the center of mass must obey
In the inertial center-of-mass frame,
where is the moment of inertia and is the angular velocity of the rotation. Nearly all of the mass is in the two compact (much smaller than ) nuclei, so and . It is convenient to rewrite this as
is the reduced mass of the molecule.
The rotational kinetic energy associated with this angular momentum is
The quantization of angular momentum to integer multiples of implies that rotational energy is also quantized. The corresponding energy eigenvalues of the Schrödinger equation are
Note the inverse relation between permitted rotational energies and the moment of inertia . If the upper-level rotational energy is much higher than , few molecules will be collisionally excited to that level and the line emission from molecules in that level will be very weak. For example, the minimum rotational energy of the small and light H molecule is equivalent to a temperature K, which is much higher than the actual temperature of most interstellar H. Only relatively massive molecules are likely to be detectable radio emitters in very cold (tens of K) molecular clouds.
Quantization of rotational energy implies that changes in rotational energy are quantized. The energy change of permitted transitions is further restricted by the quantum-mechanical selection rule
Going from to releases energy
The frequency of the photon emitted during this rotational transition is
where is the angular-momentum quantum number corresponding to the upper energy level. In terms of the molecular reduced mass and equilibrium nuclear separation ,
Thus a plot of the radio spectrum of a particular molecular species in an interstellar cloud will look like a ladder (Figure 7.13) whose steps are all harmonics of the fundamental frequency that is determined solely by the moment of inertia of that species. The relative intensities of lines in the ladder depend on the temperature of the cloud. Since , the lowest frequency of line emission depends on the mass and size of the molecule. Large, heavy molecules in cold clouds may be seen at centimeter wavelengths, but smaller and lighter molecules emit only at millimeter wavelengths.
For example, the laboratory spectrum of the CO carbon-monoxide molecule shows that the fundamental transition emits a photon at GHz. (See the online spectral-line catalog called Splatalogue88http://www.splatalogue.net/. for accurate frequencies of radio spectral lines.) The distance between the C and O nuclei can be estimated from
where the reduced mass is
Thus the equilibrium distance between the C and O nuclei is
The centrifugal forces acting on the nuclei increase as the molecule spins more rapidly, so a nonrigid bond will stretch and will increase slightly with . Spectral lines emitted by more rapidly rotating CO molecules will have frequencies slightly lower than the harmonics of the line: GHz and the actual frequency is GHz. Chemists use these line frequencies to determine , and the difference between and is a measure of the stiffness of the carbon–oxygen chemical bond. Since the actual frequency is only slightly less than the harmonic frequency, the stiffness of the “spring” connecting the atoms is quite high. Consequently the fundamental vibrational frequency of the CO molecule is much higher than the fundamental rotational frequency, and CO emits vibrational lines at mid-infrared wavelengths m.
Equation 7.111 can be used to calculate frequencies for molecules containing rare isotopes (e.g., CO) that might be more difficult to measure in the lab:
so we expect
The actual CO frequency is 110.201354 GHz.
Polar diatomic molecules emit a harmonic series of radio spectral lines at millimeter wavelengths. Bigger and heavier linear polyatomic molecules have ladders of lines starting at somewhat lower frequencies. Nonlinear molecules such as the symmetric-top ammonia (NH) with two distinct rotational axes have more complex spectra consisting of many parallel ladders (Figure 7.14).
Molecules are excited into states by ambient radiation and by collisions in a dense gas. The minimum gas temperature needed for significant collisional excitation is
where is the rotational energy of the upper energy level for the transition. Thus a minimum gas kinetic temperature
is required to excite the transition at frequency . For example, the minimum gas temperature needed for significant excitation of the CO line at GHz (Figure 7.15) is
The values for many molecular lines may be found in the online spectral-line catalog Splatalogue. If K, then radiative excitation by the cosmic microwave background is ineffective.
Larmor’s formula for a time-varying dipole can be applied to estimate the average power radiated by a rotating polar molecule. The electric dipole moment of any charge distribution is defined as the integral
over the volume containing the charges. In the case of two point charges and with separation ,
When a molecule rotates with angular velocity , the projection of the dipole moment perpendicular to the line of sight varies with time as . Equation 7.101 states that
Equation 2.133 from the derivation of Larmor’s formula states that each charge contributes
to the radiated electric field at distance from the source. The fields from both charges add in phase because , so the total radiated field is
Thus the instantaneous power emitted is
and the time-averaged power is
This can be expressed as
defines the mean electric dipole moment. For the radiative transition between upper and lower energy levels U and L, the spontaneous emission coefficient is
The value of for the transition of a linear rotating molecule with dipole moment is
(This equation reflects the complex angular wave functions involved, and we won’t derive it here.)
Dipole moments are often expressed in debye units defined by , where is a typical charge imbalance for a polar molecule and is a typical value for the interatomic separation . For example, the CO molecule has dipole moment .
For example, the spontaneous emission coefficient for the CO line at GHz is
This is close to the more accurate Splatalogue value, s.
The typical time s for a CO molecule to emit a photon spontaneously may be longer than the average time between molecular collisions in an interstellar molecular cloud, so CO can approach LTE with the excitation temperature of the line being nearly equal to the kinetic temperature of the molecular cloud. For any molecular transition, there is a critical density defined by
at which the radiating molecule suffers collisions at the rate equal to the spontaneous emission rate . Typical collision cross sections are cm, and the average velocity of the abundant H molecules is cm s if K. Thus the critical density of the CO transition is
Many Galactic molecular clouds have higher densities than this, so Galactic CO emission is strong and widespread. Also, photons from a particular transition may be repeatedly absorbed and reemitted within the molecular cloud. Such line trapping lowers the effective emission rate and reduces the effective value of needed for LTE.
Whether or not a molecular cloud is in LTE, Equation 7.67 with the temperature replaced by excitation temperature can be used to calculate the line opacity coefficient:
At the line center frequency , Equation 7.37 gives
where is the FWHM line width expressed as a frequency and is the line width in velocity units (e.g., km s). The line-center optical depth is along the line of sight, and the column density is along the line of sight, so
In the Rayleigh–Jeans approximation , the brightness-temperature difference between the line center and the nearby off-line continuum is
where is the brightness temperature of any background continuum emission (e.g., from the CMB). In the limit of low line optical depth, the line brightness
is proportional to the column density . The spontaneous emission coefficient is proportional to (Equation 7.131), so spectral lines tend to have higher brightness temperatures and become more prominent at higher radio frequencies (Figure 7.16).
The hydrogen molecule H is by far the most abundant molecule in interstellar space. Unfortunately, it is symmetric so its dipole moment is zero. Observable but comparatively rare polar molecules such as CO are only tracers, and total column densities of molecular gas must be estimated indirectly by the use of a fairly uncertain CO to H conversion factor relating column density in cm to CO velocity-integrated line brightness in . The best current value for our Galaxy is
is probably higher in galaxies with low metallicity and lower in starburst galaxies [bol13]. Figure 8.12 shows the distribution of CO emission tracing molecular gas and obscured star formation in the interacting starburst galaxies NGC 4038/9.
Isotopologues are molecules that differ only in isotopic composition; that is, only in the numbers of neutrons in their component atoms. Lines of the most abundant isotopologue of carbon monoxide, , are often optically thick, so the line brightness temperature approaches the molecular gas kinetic temperature and is nearly independent of column density. Lines of rarer isotopologues such as or are usually optically thin and can be used to measure the column densities needed to estimate the total mass of molecular gas in a source. Intensity ratios of optically thin lines from different levels can be used to measure excitation temperature, which is close to the kinetic temperature in LTE.
Transitions with high emission coefficients (e.g., the HCN (hydrogen cyanide) line at GHz has s) are collisionally excited only at very high densities ( cm for HCN ). They are valuable for highlighting only the very dense gas directly associated with the formation of individual stars.
The discovery of ammonia (NH) in the direction of the Galactic center by crt+68 immediately led to the realization that the interstellar medium must contain regions much denser than previously expected because the critical density needed to excite the NH line is cm.
Hydrogen is the most abundant element in the interstellar medium (ISM), but the symmetric H molecule has no permanent dipole moment and does not emit detectable spectral lines at radio frequencies. Neutral hydrogen (Hi) atoms are abundant and ubiquitous in low-density regions of the ISM. They are detectable in the cm ( MHz) hyperfine line. Two energy levels result from the magnetic interaction between the quantized electron and proton spins. When the relative spins change from parallel to antiparallel, a photon is emitted.
The Hi line center frequency is
where is the nuclear -factor for a proton, is the dimensionless fine-structure constant, and is the hydrogen Rydberg frequency (Equation 7.12).
By analogy with the emission coefficient of radiation by an electric dipole
the emission coefficient of this magnetic dipole is
where is the mean magnetic dipole moment for Hi in the ground electronic state (). The magnitude is called the Bohr magneton, and its value is
Thus the emission coefficient of the 21-cm line is only
so the radiative half-life of this transition is very long:
Such a low emission coefficient implies an extremely low critical density (Equation 7.135) , so collisions can easily maintain this transition in LTE, even in the outermost regions of a normal spiral galaxy and in tidal tails of interacting galaxies.
Regardless of whether the Hi is in LTE or not, we can define the Hi spin temperature (the Hi analog of the molecular excitation temperature defined by Equation 7.70) by
where the statistical weights of the upper and lower spin states are and , respectively. Note that
is very small for gas in LTE at K, so in the ISM
Inserting these weights into Equation 7.67 gives the opacity coefficient of the line:
where is the number of neutral hydrogen atoms per cm. The neutral hydrogen column density along any line of sight is defined as
The total opacity of isothermal Hi is proportional to the column density. If , then the integrated Hi emission-line brightness is proportional to the column density of Hi and is independent of the spin temperature because and in the radio limit . Thus can be determined directly from the integrated line brightness when . In astronomically convenient units it can be written as
where is the observed 21-cm-line brightness temperature at radial velocity and the velocity integration extends over the entire 21-cm-line profile. Note that absorption by Hi in front of a continuum source with continuum brightness temperature , on the other hand, is weighted in favor of colder gas (Figure 7.17).
The equilibrium temperature of cool interstellar Hi is determined by the balance of heating and cooling. The primary heat sources are cosmic rays and ionizing photons from hot stars. The main coolant in the cool atomic ISM is radiation from the fine-structure line of singly ionized carbon, Cii, at m. This line is strong only when the temperature is at least
so the cooling rate increases exponentially above
The actual kinetic temperature of Hi in our Galaxy can be estimated from the Hi line brightness temperatures in directions where the line is optically thick () and the brightness temperature approaches the excitation temperature, which is close to the kinetic temperature in LTE. Many lines of sight near the Galactic plane have brightness temperatures as high as 100–150 K, values consistent with the temperature-dependent cooling rate.
Neutral hydrogen gas in the disk of our Galaxy moves in nearly circular orbits around the Galactic center. Radial velocities measured from the Doppler shifts of Hi cm emission lines encode information about the kinematic distances of Hi clouds, and the spectra of Hi absorption in front of continuum sources can be used to constrain their distances also. Hi is optically thin except in a few regions near the Galactic plane, so the distribution of hydrogen maps out the large-scale structure of the whole Galaxy, most of which is hidden by dust at visible wavelengths.
Figure 7.18 shows a plan view of the Galactic disk. The Sun () lies in the disk and moves in a circular orbit around the Galactic center. The distance to the Galactic center kpc and the Sun’s orbital speed km s have been measured by a variety of means [rei93]. All Hi clouds at galactocentric distance are assumed to be in circular orbits with angular velocity , where is a monotonically decreasing function of . For cloud 1 at galactocentric azimuth on the line of sight at Galactic longitude , the observed radial velocity relative to the Sun is given by
Using the trigonometric identities and we obtain
To apply this equation, we need to determine the rotation curve . The maximum radial velocity on the line of sight at longitude is called the “terminal velocity” . Since decreases with , this velocity occurs at the minimum where the orbit is tangent to the line of sight:
We can determine the rotation curve from measurements of spanning a wide range of and thus of .
The 1420 MHz Hi line is an extremely useful tool for studying gas in the ISM of external galaxies and tracing the large-scale distribution of galaxies in the universe because Hi is detectable in most spiral galaxies and in some elliptical galaxies.
Because cm is such a long wavelength, many galaxies are unresolved by single-dish radio telescopes. For example, the half-power beamwidth of the 100-m GBT is about 9 arcmin at cm. Thus a single pointing is sufficient to obtain a spectral line representing all of the Hi in any but the nearest galaxies.
The observed center frequency of the Hi line can be used to measure the radial velocity of a galaxy. The radial velocity of a galaxy is the sum of the recession velocity caused by the uniform Hubble expansion of the universe and the “peculiar” velocity of the galaxy. The radial component of the peculiar velocity reflects motions caused by gravitational interactions with nearby galaxies and is typically km s in magnitude. The Hubble velocity is proportional to distance from the Earth, and the Hubble constant of proportionality has been measured as [pl15]. If the radial velocity is significantly larger than the radial component of the peculiar velocity, the observed Hi frequency can be used to estimate the distance to a galaxy.
Beware that astronomers still use inconsistent radial velocity conventions that were established when most observed radial velocities were much less than the speed of light. The approximation
where is the line frequency in the source frame and is the observed frequency, was used to define the radio velocity for any as
because radio astronomers measure frequencies, not wavelengths. Optical astronomers measure wavelengths, not frequencies, so the nonrelativistic approximation
is the basis for the optical velocity defined for any by
where is the redshift defined by Equation 2.127. The optical and radio velocity conventions are not exactly the same, and neither agrees with the relativistically correct radial velocity calculated from Equation 5.142. Occasionally an observer confuses velocity conventions, fails to center the observing passband on the correct frequency, and ends up with only part of the Hi spectrum of a galaxy. Outside the local universe () the concept of distance itself becomes more complicated. To calculate distances to astronomical objects with larger redshifts, see hog99.
For example, the Hi emission-line profile (Figure 7.19) of the galaxy UGC 11707 can be used to estimate its distance . The observed line center frequency is MHz, so the “radio” and “optical” velocities are
Using the “optical” velocity gives
If the Hi emission from a galaxy is optically thin, then the integrated line flux is proportional to the mass of Hi in the galaxy, independent of the unknown Hi temperature. It is a straightforward exercise to derive from Equation 7.155 the relation
for the total Hi mass of a galaxy. The integral over the line is called the line flux and is usually expressed in units of Jy km s. For example, to estimate the Hi mass of UGC 11707, assume . The single-dish Hi line profile of UGC 11707 (Figure 7.19) indicates a line flux
Small statistical corrections for nonzero can be made from knowledge about the expected opacity as a function of disk inclination, galaxy mass, morphological type, etc.
A well-resolved Hi image of a galaxy yields the total mass enclosed within radius of the center if the gas orbits in circular orbits.
If the mass distribution of every galaxy were spherically symmetric, the gravitational force at radius would equal the gravitational force of the enclosed mass . This is not a bad approximation, even for disk galaxies. Thus for gas in a circular orbit with orbital velocity ,
where is the mass enclosed within the sphere of radius and is the orbital velocity at radius , so
Note that the velocity is the full rotational velocity, not just its radial component , where is the inclination angle between the galaxy disk and the line of sight. The inclination angle of a thin circular disk can be estimated from the axial ratio
where and are the minor- and major-axis angular diameters, respectively. Converting from CGS to astronomically convenient units yields
and we obtain the total galaxy mass inside radius in units of the solar mass:
Thus the total mass of UGC 11707 (Figure 7.20) can be estimated from
UGC 11707 is a relatively low-mass spiral galaxy.
This “total” mass is really only the mass inside the radius sampled by detectable Hi. Even though Hi extends beyond most other tracers such as molecular gas or stars, it is clear from plots of Hi rotation velocities versus radius that not all of the mass is being sampled, because we don’t see the Keplerian relation which indicates that all of the mass is enclosed within radius . Most rotation curves, one-dimensional position-velocity diagrams along the major axis, are flat at large , suggesting that the enclosed mass as far as we can see Hi. The large total masses implied by Hi rotation curves provided some of the earliest evidence for the existence of cold dark matter in galaxies.
Because detectable Hi is so extensive, Hi is an exceptionally sensitive tracer of tidal interactions between galaxies. Long streamers and tails of Hi trace the interaction histories of pairs and groups of galaxies. See Figure 8.11 showing Hi in the M81 group of galaxies and Figure 8.12 revealing long tidal tails of Hi from the “Antennae” galaxy pair NGC 4038/9.
Another application of the Hi spectra of galaxies is determining departures from smooth Hubble expansion in the local universe via the Tully–Fisher relation. Most galaxies obey the empirical luminosity–velocity relation [tf77]:
where is the maximum rotation speed. Arguments based on the virial theorem can explain the Tully–Fisher relation if all galaxies have the same central mass density and density profile, differing only in scale length, and also have the same mass-to-light ratio. Thus a measurement of yields an estimate of that is independent of the Hubble distance . The Tully–Fisher distance can be calculated from this “standard candle” and the apparent luminosity. Apparent luminosities in the near infrared (m) are favored because the near-infrared mass-to-light ratio of stars is nearly constant and independent of the star-formation history, and because extinction by dust is much less than at optical wavelengths. Differences between and are ascribed to the peculiar velocities of galaxies caused by intergalactic gravitational interactions. The magnitudes and scale lengths of the peculiar velocity distributions are indications of the average density and clumpiness of mass on megaparsec scales.
Most of the baryonic matter in the early universe was fully ionized hydrogen and helium gas, plus trace amounts of heavier elements. This smoothly distributed gas cooled as the universe expanded, and the free protons and electrons recombined to form neutral hydrogen at a redshift when the age of the universe was about years. The hydrogen remained neutral during the dark ages prior to the formation of the first ionizing astronomical sources—massive () stars, galaxies, quasars, and clusters of galaxies—by gravitational collapse of overdense regions. These astronomical sources gradually started reionizing the universe when it was several hundred million years old () and completely reionized the universe by the time it was about years old (). This era is called the epoch of reionization.
The highly redshifted Hi signal was nearly uniform during the dark ages, and it developed structure on angular scales up to several arcmin when the first astronomical sources created bubbles of ionized hydrogen around them. As the bubbles grew and merged, the Hi signal developed frequency structure corresponding to the redshifted Hi line frequency. The characteristic size of the larger bubbles reached about 10 Mpc at and produced Hi signals having angular scales of several arcmin and covering frequency ranges of several MHz. These Hi signals encode unique information about the formation of the earliest astronomical sources.
The Hi signals produced by the EOR will be very difficult to detect because they are weak (tens of mK), relatively broad in frequency, redshifted to low frequencies ( MHz) plagued by radio-frequency interference and ionospheric refraction, and lie behind a much brighter (tens of K) foreground of extragalactic continuum radio sources. Nonetheless, the potential scientific payoff is so great that several groups around the world are developing instruments to detect the Hi signature of the EOR. Two such instruments are PAPER (Precision Array to Probe the Epoch of Reionization) [Par14] and the Murchison Widefield Array [tin13] shown in Figure 8.7.