Processing math: 100%
Skip to main content
Library homepage
 

Text Color

Text Size

 

Margin Size

 

Font Type

Enable Dyslexic Font
Physics LibreTexts

8.6: Saha's Equation

( \newcommand{\kernel}{\mathrm{null}\,}\)

Consider the reversible reaction

AA++e,

where A is a neutral atom and A+ its first ion. Let N0, Ni, Ne be the numbers of neutral atoms, ions and electrons respectively, held in a box of volume V. Then

NeNiN0=S(T,P)

is the Saha function. It is a function of temperature and pressure, high temperature favouring ionization and high pressure favouring recombination. (Students of chemistry will recognize this as an application of Le Chatelier's principle.) The equation tells us the relative numbers of the three types of particle (i.e. the degree of ionization) in an equilibrium situation when the number of ionizations per second is equal to the number of recombinations per second.

Now the number of particles in a given energy level is proportional to the Boltzmann factor for that level, and the total number of particles is proportional to the sum of the Boltzmann factors for all the levels - i.e. to the partition function. Thus, using the symbol Q to denote partition functions, we have, for Saha's equation:

NeNiN0=QeQiQ0.

The partition function is the sum of the Boltzmann factors over all the states, translational and internal (electronic). The total energy of a particle is the sum of its translational and internal energies, so that total partition function is the product of its translational and internal partition functions, for which we shall use the symbol u. Therefore

NeNiN0=(2πmkTh2)32Vueuiu0.

Here m=memim0, which is very little different from me.

The internal partition function of the electron is merely its statistical weight 2s+1=2. In calculating the internal partition functions of the atom and ion, it is, of course, essential to measure all energies from the same level - the ground level of the neutral atom, for example. You would not be so foolish, of course, as to calculate the partition function of the atom by using, in the Boltzmann sum, the energies of each level above the ground level of the atom, and then calculate the partition function of the ion by using the energies above some quite different level, such as the energy of the ground stage of the ion, would you? Most of us probably would be so foolish; we should, of course, add χi (the ionization energy) to each level of the ion. If we do calculate the partition function of the ion using energies above the ground level of the ion, we must then multiply the result, ui, by eχi/kT. Being even more careful, we should remember, from section 8.5, that the ionization energy is in practice lowered by a small amount Δχi. Finally, let us now use the symbol n to denote number of particles per unit volume (so that n=N/V, and on no account to be confused with the quantum numbers n used in section 8.5) and we arrive at the usual form for Saha's equation:

nenin0=(2πmkTh2)322uiu0exp((χiΔχi)kT).

It might be noted that Δχi is a function of ne, which leads to a slight complication in the computation of Saha's equation, which we shall encounter in one of the problems that follow.

It should be remarked that Saha's equation played an extremely important role in the understanding of the stellar spectral sequence. As is well known, the sequence of spectral types OBAFG… is a result of the degree of ionization and excitation of the elements as a function of temperature, and the difference in the degree of ionization between main sequence stars and giants of a given temperature is the result of the higher degree of ionization in the relatively low pressure atmospheres of the giant stars. This could not be understood at a quantitative level until the development of Saha's equation and its application to stellar atmospheres by Saha, Fowler, Milne and Payne (later Payne-Gaposhkin) in the 1920s.

Problems.

Problem 1. Verify that equation ??? balances dimensionally.

Problem 2. In the bad old days, we did calculations using logarithm tables! You probably will never have to do that, but this exercise will nevertheless turn out to be useful. Show that Saha's equation can be written as

lognenin0=A+32logT+loguiu0B(ViΔViT)

and determine the numerical value of the constants A and B if n is expressed in m3, T is in K, Vi is the ionization potential in volts, and logarithms are common logarithms (i.e. to the base 10). I make them A=21.68 and B=5039.7. It is common also to write equation ??? in the form

lognenin0=27.2432logθ+loguiu0θ(ViΔVi),

in which θ=5039.7/T.

Problem 3. You are going to calculate the Saha function for hydrogen, and so you need the partition functions for the electron, the neutral hydrogen atom and the hydrogen ion. The electron is easy. Its spin is 1/2, so its partition function is 2, as already discussed and indeed already incorporated into Saha's equation. The partition function for neutral hydrogen can be taken to be the statistical weight of its n=1 shell, which is 2n2=2. Now what about the statistical weight of the hydrogen ion? The hydrogen ion is a proton, which has spin =1/2. Therefore, as for the electron, should the partition function be 2? The answer is no! When calculating the Saha equation for hydrogen, you should take the partition function of the proton to be 1. This probably seems entirely illogical and you are probably quite sure that I am wrong. But before coming to this conclusion, read again the last paragraph of section 8.4 and the first paragraph of section 8.5i. You can include the effect of the proton spin provided that you take it into account for both H and H+. If you insist that the statistical weight of the proton is 2, you must also insist that the statistical weight of the n=1 shell of H is 4n2=4. As with the Boltzmann equation, you usually need not be concerned with nuclear spin - it cancels out in ui/u0. There is nothing, however, to cancel out the factor 2 for the electron partition function.

Problem 4. A kilogram of water is contained is a box of volume one cubic metre. The box is made of material that won't melt or vaporize. (!) Draw a graph showing, as a function of temperature, from 10,000 K to 20,000 K, the logarithm of the numbers per unit volume of each of the species e, H, H+, O, O+. (Assume there are no higher ionization stages, and no molecules.)

This is going to be difficult, and will involve some computation. The first thing you are going to need to do is to calculate the Saha functions of the species involved as a function of temperature., and for this, you will need the partition functions . For the electron, the partition function is 2 and is already incorporated in the Saha equation. As discussed in Problem 3, the partition function of H+ should be taken to be 1, and, as for the electron, it is mercifully independent of temperature. For neutral H, the statistical weight of the ground shell is 2, and, as you will already have discovered if you carried out the little calculation in section 8.5ii, you can probably take this to be the partition function at all temperatures between 10,000 and 20,000 K. The ground term of neutral oxygen is 3P, and the first excited term is 1D at a term value of 15868 cm1. Probably only the ground term contributes appreciably to the partition function, but you might want to check that. The ground term of O II is 4So and the first excited term is a long way up. Thus, for all species, you can probably take the partition function to be the statistical weight of the ground term and independent of temperature. You are very lucky that I did not throw a metal into the mix!

The ionization limits of H and O are at 101678 and 109837 cm1 respectively. We don't yet know by how much, Δχi the ionization energy is lowered, so we'll ignore Δχi to begin with. You now have enough information to calculate the Saha functions for hydrogen and for oxygen as a function of temperature.

We have been asked to calculate the numbers of five species and we are therefore going to need four equations. We started with one kilogram of water. We know that 18 grams of water contain A molecules, where A is Avogadro's number. (You may want to look up the exact molar mass of H2O.) Thus, we started with 1000A/18 molecules. Let's call that number N.

Now see if you agree with the following five equations.

NH+NH+=2N

NO+NO+=N

NH++NO+=Ne

NeNH+=VSHNH

NeNO+=VSONO

We do indeed have five equations in the five unknowns, and all other quantities in the equations are known. The first two express the stoichiometry of the water molecule. The third expresses electrical neutrality. The last two are Saha's equations for hydrogen and oxygen. They five equations are not very difficult equations, although the last two are nonlinear, which makes them slightly awkward. We can make them a little easier. I am going to re-write them:

NH+NH+2N=0

NO+NO+N=0

NH++NO+Ne=0

VSHNHNeNH+=0

VSONONeNO+=0

If we now, just for the moment, suppose that Ne is a quantity whose value is known. In that case, equations ??? to ??? would comprise five linear equations in just four unknowns. The condition for these equations to be consistent is that the determinant of the coefficients and the constant terms be zero:

11002N0011N0101Ne=0VSHNe00000VSONe0

If I have done my algebra right (you might like to check this), this is a cubic equation in Ne:

N3e+V(SO+SH)N2e+V[VSHSO(2SH+SO)N]Ne3V2SHSON=0.

You can now solve this for Ne. Remember, however, that, in calculating the Saha functions, you ignored the lowering of the ionization potential ΔVi, so what you have found is but a first approximation for Ne. This is a function of the electron density. One rather old formula, due to Unsöld, is

ΔVi=7×105n1/3e,

where Vi is in volts and ne(=Ne/V) is in m-3 (indicating that the lowering of the potential is proportional to the distance between the charged particles.). If you can find a more recent and perhaps better formula, by all means use it. You can use your first approximation for ne to calculate the Saha functions, and then repeat the entire calculation and continue to iterate until you obtain an acceptably small change in ne. It is probable that only one iteration will be necessary. Having found ne (and Ne) it is then straightforward to find the remaining unknowns from equation ??? to ???. Just make sure that you make use of all five equations as a check for arithmetic mistakes.

Problem 5. A bottle of methyl cyanate CH3CNO is held in a cylinder with a movable piston such that the pressure remains constant at one pascal. Heat it up to some temperature at which you can assume that the only species present are electrons and the neutral atoms and singly ionized stages of C, H, O and N - i.e. no higher ionization stages and no molecules. The cylinder and piston are not allowed to melt - they are required only to provide constant pressure conditions, in contrast to the previous problem in which you had a constant volume. You have nine unknowns, and you will need nine equations. In this problem you are asked merely to write down the nine required equations, not necessarily to solve them unless you particularly want to. There will be Saha equations for each element, an equation expressing electrical neutrality, and four equations reflecting the stoichiometry. You may assume the ideal gas equation P=nkT.


This page titled 8.6: Saha's Equation is shared under a CC BY-NC 4.0 license and was authored, remixed, and/or curated by Jeremy Tatum via source content that was edited to the style and standards of the LibreTexts platform.

Support Center

How can we help?