Saha Equation Calculator

The Saha equation, derived by Meghnad Saha in 1920, gives the equilibrium ratio of ionized to neutral atoms in a gas at a given temperature and density, assuming the gas is in local thermodynamic equilibrium. It combines statistical mechanics (counting the available quantum states) with the Boltzmann factor for the ionization energy.

For hydrogen, n_e n_II / n_I = (2πm_e k_B T/h²)^(3/2) e^(−χ/k_B T), where n_e, n_II, and n_I are the electron, ion, and neutral densities, χ = 13.6 eV is the hydrogen ionization energy, and the degeneracy-ratio prefactor for ground-state hydrogen simplifies to exactly 1.

Before Saha's 1920 work, astronomers could not explain why stars of very similar composition showed wildly different absorption line strengths. Saha showed that a star's spectral type is primarily determined by its surface temperature, through the ionization state of the elements in its atmosphere, not by fundamentally different chemical composition, resolving a major puzzle in stellar astrophysics.

The Boltzmann factor e^(−χ/k_B T) changes exponentially with temperature, so even a modest change in T can shift the exponent by several units, swinging the ionization fraction by many orders of magnitude. This is why stars with only slightly different surface temperatures can show dramatically different spectral line strengths for the same element.

The ionization fraction x is the proportion of hydrogen atoms that are ionized, with n_e = n_II = x·n_total and n_I = (1−x)·n_total, where n_total is the combined density of neutral and ionized hydrogen entered into the calculator.

The photosphere (about 5,800 K, or roughly 0.5 eV) is hot by everyday standards but still far below hydrogen's 13.6 eV ionization energy, so the exponential Boltzmann suppression keeps the ionization fraction extremely low, around one part in ten thousand. Only trace elements with much lower ionization energies (like sodium or calcium) become significantly ionized at photospheric temperatures.

LTE means the gas is dense enough and collisions frequent enough that the local ionization state matches what a Boltzmann/Saha equilibrium at the local temperature and density would predict, even if the gas as a whole is not in complete equilibrium (like radiating starlight). It is an excellent approximation for stellar interiors and photospheres, but breaks down in tenuous plasmas like the solar corona or the interstellar medium, where radiative and collisional rates fall out of balance.

The equation as usually written gives the ratio n_e n_II/n_I, but since n_e = n_II = x·n and n_I = (1−x)·n for hydrogen, substituting turns it into a quadratic equation in x, n·x² + RHS·x − RHS = 0, which this calculator solves exactly using the quadratic formula.

Yes, in general it applies to any element, with the ionization energy χ and the degeneracy ratio of the relevant states adjusted accordingly. This calculator is scoped specifically to hydrogen, since it has the simplest and most astrophysically important form, with the degeneracy prefactor working out to exactly 1.

The Saha equation determines the initial ionization state of a plasma formed by heating a gas, connecting classical statistical mechanics to the birth of the fourth state of matter. It remains a standard tool for estimating ionization balance in stellar atmospheres, laboratory discharges, and any thermal plasma source.