advanced30 min studyLesson 888 in Physics

Photoionization equilibrium and H II regions

Distinguish total-field ionization balance from the Case B photon budget and derive a static hydrogen ionized-region scale.

On this page

What you’ll understand

  • State whether diffuse photons are included in a photoionization rate.
  • Derive the Case B Strömgren radius from a photon-number budget.
  • Compare recombination times and density-dependent recombination brightness.

An H II region is hydrogen kept ionized by photons above the ionization threshold. Its ionization state need not follow the Saha equation: radiation from a nearby source can maintain ionization in a relatively cool dilute gas. The simplest model balances ionizations against recombinations at a specified temperature. Careful photon accounting matters because a recombination can itself emit another ionizing photon. We will distinguish the full local radiation field from the commonly used on-the-spot approximation before estimating an ionized region's radius.

Ionization balance counts rates rather than imposing LTE

Let neutral density be n HI, proton density n p and electron density n e, all per cubic metre. The photoionization rate per neutral atom is an integral of photon intensity times the atomic cross section over frequencies above threshold. If J includes both direct and diffuse photons, steady balance uses the total recombination coefficient alpha A. That coefficient has cubic metres per second units, so the recombination rate per volume has inverse cubic metres per second units. Collisional ionization is neglected in this photoionization-dominated example.

The temperature must be supplied or found from a separate heating and cooling calculation. It controls recombination coefficients but does not by itself fix the external ionizing spectrum. A Maxwellian electron distribution can coexist with strongly non-LTE hydrogen populations. Density structure and shielding also alter the neutral fraction even when the ionized volume is nearly fully ionized. Thus an ionization equilibrium is a rate balance with declared radiation and atomic processes; it is not synonymous with thermodynamic equilibrium or a guarantee of uniform temperature.

Γtot=4π∫ν0∞Jν,totσνhν dν,nHIΓtot=nenpαA(T),αB=αA−αground\Gamma_{\mathrm{tot}}=4\pi\int_{\nu_0}^{\infty}\frac{J_{\nu,\mathrm{tot}}\sigma_\nu}{h\nu}\,d\nu,\qquad n_{\mathrm{HI}}\Gamma_{\mathrm{tot}}=n_en_p\alpha_A(T),\qquad \alpha_B=\alpha_A-\alpha_{\mathrm{ground}}
Gamma-total has s to the minus one units, J-total includes diffuse ionizing radiation, and sigma is in square metres. Alpha-ground counts direct recombinations to the ground state. Replacing alpha-A with alpha-B while retaining the full diffuse Gamma would double count recycled photons.↔ Scroll wide equations horizontally.

Case B removes locally recycled photons from the budget

A recombination directly to the hydrogen ground state emits an ionizing photon. In an optically thick ionizing medium, the on-the-spot approximation assumes that photon causes another nearby ionization. These two local events cancel from the budget of stellar photons. The remaining effective recombination coefficient is alpha B, excluding direct ground-state recombinations. Correspondingly the reduced local equation uses the external stellar photoionization rate, rather than the total direct-plus-diffuse rate. This approximation is inappropriate when the recombination ionizing photons escape freely.

For a static, uniform, pure-hydrogen sphere with a sharp ionization front, set n e and n p approximately equal to the hydrogen density inside. If all emitted source ionizing photons are absorbed by hydrogen, integrating the Case B recombination rate gives the Strömgren radius. The source quantity Q H is photons per second, not its bolometric energy luminosity. Dust absorption, leakage, helium, clumps and a nonuniform density change this budget. Optical thickness to ionizing photons also does not require an optically thick radio continuum or visible recombination line.

nHIΓext=nenpαB(on the spot),RS=(3QH4παBnH2)1/3,trec,0=1αBne,0n_{\mathrm{HI}}\Gamma_{\mathrm{ext}}=n_en_p\alpha_B\quad(\text{on the spot}),\qquad R_S=\left(\frac{3Q_H}{4\pi\alpha_Bn_H^2}\right)^{1/3},\qquad t_{\mathrm{rec},0}=\frac1{\alpha_Bn_{e,0}}
The radius assumes an ionization-bounded, fully ionized interior of uniform hydrogen density, fixed temperature, no dust and no photon leakage. The recombination time is an initial characteristic time, not a claim that a freely recombining hydrogen ionization fraction decays exponentially.↔ Scroll wide equations horizontally.

Recombination and hydrodynamic evolution are different clocks

After a source switches off, pure hydrogen supplies its own electrons. At fixed total density and temperature, with the same effective Case B recycling approximation maintained, the ionized fraction satisfies a loss rate proportional to its square. Its solution is a reciprocal linear function of time, not an exponential with constant electron density. An exponential can be appropriate for a trace ion in an externally fixed electron bath, but that is a different physical population assumption. The distinction matters when translating fading emission into an elapsed time.

Recombination-line emissivity usually weights electron times proton density with a line-specific effective coefficient. That coefficient is not automatically the total alpha B coefficient. For fixed geometry, temperature and line coefficient, brightness therefore falls as the ionized fraction squared in the simple recombining model. A real heated ionized region can expand as pressure accelerates the surrounding gas; its eventual radius need not equal a static initial Strömgren estimate. Comparing recombination and sound-crossing times helps decide when a fixed-density ionization calculation is a useful approximation.

A synthetic ionization-bounded hydrogen sphere

Take QH=(8π/3)×1045 s−1Q_H=(8\pi/3)\times10^{45}\,\mathrm{s^{-1}}, nH=108 m−3n_H=10^8\,\mathrm{m^{-3}} and a prescribed fixed-temperature αB=2×10−19 m3 s−1\alpha_B=2\times10^{-19}\,\mathrm{m^3\,s^{-1}}. Find the static Strömgren radius and initial recombination time. Use one year = 3.15 × 10^7 s.

  1. The volume recombination rate is αBnH2=2×10−3 m−3 s−1\alpha_Bn_H^2=2\times10^{-3}\,\mathrm{m^{-3}\,s^{-1}}. This is a photon-number budget, so no mean photon energy multiplies Q.

  2. The cubed radius is 3QH/(4παBnH2)=1048 m33Q_H/(4\pi\alpha_Bn_H^2)=10^{48}\,\mathrm{m^3}. Taking the cube root gives RS=1016 mR_S=10^{16}\,\mathrm m.

  3. Initially n e equals n H, giving trec,0=1/(2×10−19×108)=5×1010 st_{\mathrm{rec},0}=1/(2\times10^{-19}\times10^8)=5\times10^{10}\,\mathrm s, or about 1590 years.

The result

The static idealized radius is 10^16 m and its initial recombination time is about 1.59 × 10^3 years. Expansion, dust and leakage have been excluded.

Explore the model

Interactive lab

Interactive lab · function

A recombination emissivity weights a density product

Reset and set BOTH Coefficient a = 1 and Position x = 1; predict y = 1. For a fresh comparison, Reset and set a = 1 and x = 2; predict y = 4. Interpret x as normalized equal electron and proton density at fixed recombination coefficient.

A recombination emissivity weights a density productThis lab uses only y = ax². Playback traces the curve within the fixed integration interval from 0 to x; it does not estimate numerical error or change the chosen point. Current readings: y at chosen x: 0.608; Tangent slope · 2ax: 1.35; Signed area · ax³/3: 0.182. All parameters are adjustable with the labeled controls below.−2.4−1.201.22.4−9−4.504.59xy(x, ax²)
y = ax²TangentSigned integral
y at chosen x
0.608
Tangent slope · 2ax
1.35
Signed area · ax³/3
0.182
t = 0 s

This lab uses only y = ax². Playback traces the curve within the fixed integration interval from 0 to x; it does not estimate numerical error or change the chosen point. Parameter changes pause playback and recompute the model at the current elapsed time. Reset restores both parameters and returns time to zero. Playback starts paused. Adjust either parameter, then choose Play or Step.

Look forThis quadratic illustrates a local recombination-rate density scaling. It has no ionizing spectrum, diffuse photon recycling, ionization front, expanding gas, or radiation-hydrodynamics solver.

Check your understanding

In the on-the-spot Case B balance, how are direct ground-state recombinations handled?

Try it yourself

Switch off the source in a fixed-density hydrogen toy model with initial ionized fraction one. Keep T and alpha B fixed, retain the stated local recycling approximation, and neglect expansion. Find the ionized fraction and a thin recombination-line brightness ratio at one initial recombination time.

Take a moment to make a prediction or work it out on paper. Open a hint when you need one.

Give me a hint

Integrate dx/dt=−αBnHx2dx/dt=-\alpha_Bn_Hx^2 rather than imposing a constant electron density. The line brightness is proportional to x squared under the stated conditions.

Compare with the solution

Integration gives 1/x=1+αBnHt1/x=1+\alpha_Bn_Ht, hence x=1/[1+t/trec,0]x=1/[1+t/t_{\mathrm{rec},0}]. At one initial recombination time x is one half, not e to the minus one. Electron and proton densities each halve, so their product and the thin line brightness become one quarter of the initial values. This assumes unchanged path length and line coefficient as well as continued local absorption of the recycled ionizing photons; a cooling or expanding region requires additional evolution equations.

Take it with you

  • Total-field photoionization balance uses alpha A; the reduced on-the-spot external budget uses alpha B.
  • The Strömgren radius follows a photon-number budget with explicit geometry and density assumptions.
  • Pure hydrogen recombines with a changing electron density, giving reciprocal rather than exponential fading in the fixed-temperature toy model.

Read further

Original sources and deeper explanations to follow at your own pace.

Optional progress, saved in this browser when available. Revisit any lesson at any time.

Lesson not marked complete.