Hydrogen Neutral vs Ionized Fraction Calculator
Calculate the fractions of neutral hydrogen atoms (H0) and ionized hydrogen (protons, H+) using the Saha ionization equilibrium for pure hydrogen plasma.
Model assumptions: local thermodynamic equilibrium, pure hydrogen gas, ideal behavior, ground-state partition approximation. For non-LTE astrophysical plasmas, use detailed radiative transfer and collisional-radiative models.
How to Calculate the Fractions of Neutral and Ionized Hydrogen (Protons): Complete Expert Guide
If you want to calculate the fractions of neutral hydrogen and ionized hydrogen (protons), the key quantity is the ionization fraction. In many physics and astrophysics applications, this fraction tells you how much hydrogen exists as neutral atoms (H0) versus how much exists as protons (H+) plus electrons. This matters in stellar atmospheres, H II regions, plasma diagnostics, fusion-adjacent modeling, early-universe calculations, and spectroscopy.
The calculator above uses the Saha ionization equation for pure hydrogen under local thermodynamic equilibrium (LTE). In practical terms, if you know the temperature and hydrogen number density, you can estimate the neutral and ionized fractions quickly and with physically meaningful scaling.
Why this calculation is important
- Spectral interpretation: The strength of Balmer and Lyman lines depends on neutral hydrogen populations and excitation states.
- Plasma conductivity: Ionized fractions determine electron density, which controls conductivity and collision rates.
- Astrophysical environment classification: Neutral clouds, warm ionized medium, and H II regions can be distinguished by ionization fraction.
- Thermal evolution: Cooling and heating channels depend strongly on how much hydrogen is ionized.
Core physical idea
Hydrogen ionization is the process:
At equilibrium, the ratio of ionized to neutral populations follows the Saha relation:
where chi is hydrogen ionization energy (13.598 eV), T is temperature, and n_e, n_p, n_H0 are electron, proton, and neutral hydrogen number densities.
For pure hydrogen gas, n_e = n_p, and total hydrogen nuclei number density is:
Define ionized fraction x = n_p / n_total. Then neutral fraction is (1 – x). Substituting gives:
where S is the Saha right-hand side. Solving this equation yields x directly.
Step-by-step method used by the calculator
- Read input temperature T (K).
- Read total hydrogen number density n_total and convert units to m^-3 if needed.
- Compute Saha factor S using physical constants in SI.
- Compute y = S / n_total.
- Solve x from x^2 + yx – y = 0 using the positive physical root.
- Calculate:
- Ionized fraction: x
- Neutral fraction: 1 – x
- Proton density: n_p = x n_total
- Neutral density: n_H0 = (1 – x)n_total
- Electron density: n_e = n_p
Reference constants and physical values
| Quantity | Value | Why it matters |
|---|---|---|
| Hydrogen ionization energy chi | 13.598 eV (2.179 x 10^-18 J) | Sets the exponential barrier for ionization |
| Boltzmann constant k_B | 1.380649 x 10^-23 J/K | Converts thermal energy scale from temperature |
| Electron mass m_e | 9.1093837015 x 10^-31 kg | Appears in translational partition term |
| Planck constant h | 6.62607015 x 10^-34 J s | Controls quantum state density factor |
| eV conversion | 1 eV = 1.602176634 x 10^-19 J | Converts tabulated atomic energies into SI |
Real-world ionization statistics in different environments
The values below are typical literature-scale ranges used in astrophysics and plasma diagnostics. Exact values vary with radiation field, metallicity, density structure, and nonequilibrium effects.
| Environment | Typical Temperature | Typical Hydrogen Density | Approximate Ionized Fraction x |
|---|---|---|---|
| Molecular cloud interiors | 10 to 30 K | 10^2 to 10^6 cm^-3 | 10^-8 to 10^-6 |
| Cold neutral medium (CNM) | 50 to 100 K | 20 to 50 cm^-3 | about 10^-4 |
| Warm neutral medium (WNM) | 6000 to 10000 K | 0.2 to 0.5 cm^-3 | about 0.01 to 0.1 |
| Warm ionized medium (WIM) | around 8000 K | 0.05 to 0.2 cm^-3 | about 0.8 to 0.99 |
| Classical H II regions | 7000 to 10000 K | 10^2 to 10^4 cm^-3 | greater than 0.99 |
| Post-reionization intergalactic medium | about 10^4 K (order of magnitude) | very low | near unity, often greater than 0.999 |
Worked interpretation example
Suppose T = 10,000 K and n_total = 1000 cm^-3. At this temperature, thermal energy is large enough that the Saha factor can become very large compared with density, which tends to drive ionization upward. When y = S/n_total is large, x approaches 1. This means most hydrogen nuclei are protons, and neutral fraction becomes small. In this regime, small changes in temperature can strongly change neutral fraction due to the exponential term exp(-chi/k_B T).
By contrast, at low temperature such as 100 K, exp(-chi/k_B T) is effectively tiny. Then y is tiny and x is near zero. Hydrogen is overwhelmingly neutral. This sharp transition behavior is one reason ionization diagnostics are so sensitive to thermal and radiation conditions.
What controls the neutral versus ionized split the most?
- Temperature: The dominant control in Saha equilibrium because of the exponential Boltzmann factor.
- Total density: At fixed T, higher density generally suppresses ionization fraction in Saha balance.
- Radiation field: Strong UV photons can photoionize hydrogen beyond LTE expectations.
- Recombination timescale: If timescales are long, gas may be out of equilibrium.
- Composition and electrons from other species: In real gas mixtures, metals and helium alter electron budget.
When Saha is accurate and when it is not
Saha equilibrium is best when collisions and local thermodynamic equilibrium assumptions are valid, such as dense stellar photospheric layers. It becomes less reliable in diffuse nebulae dominated by radiation transport, in transient plasmas, in strongly magnetized nonequilibrium media, or where departures from Maxwellian particle distributions are important.
In photoionized nebulae, for example, ionization is often governed by balance between photoionization and recombination coefficients rather than simple LTE Saha statistics. In those situations, codes based on collisional-radiative models and radiative transfer are preferred.
Best practices for accurate calculations
- Use consistent units throughout. This calculator internally uses SI (m^-3).
- Check if your system is near LTE. If not, treat Saha output as a first estimate.
- Use realistic density ranges for your environment.
- Perform sensitivity checks by varying temperature ±10% to see fraction response.
- Compare with observational constraints such as emission line ratios, free-free continuum, or dispersion measures.
Common mistakes to avoid
- Mixing up cm^-3 and m^-3. Remember 1 cm^-3 = 10^6 m^-3.
- Assuming neutral fraction plus proton fraction can differ from 1 in pure hydrogen LTE. They should sum to 1 by definition.
- Applying Saha to strongly photoionized, low-density gas without caution.
- Ignoring partition function corrections at higher precision levels.
Authoritative references and further reading
For high-quality reference data and fundamentals, see:
- NIST Chemistry WebBook: Hydrogen thermochemical and energetic data (.gov)
- HyperPhysics (Georgia State University): Saha equation overview (.edu)
- NASA GSFC science toolbox: ionization concepts (.gov)
Bottom line
To calculate the fractions of neutral and ionized hydrogen (protons), the most direct LTE route is Saha equilibrium with temperature and total hydrogen density. The result is a physically interpretable ionized fraction x and neutral fraction (1 – x), plus number densities for protons, electrons, and neutral atoms. Use this as a robust baseline, then move to non-LTE models when radiation and kinetics dominate.