Collisional ionization equilibrium Created 2026-10-06 Updated 2026-10-14
In collisional ionization equilibrium, collision-driven ionization and recombination balance for the populated ion stages, determining their relative abundances at a given temperature. The resulting cooling curve includes the loss of line-emitting bound states as temperature increases. Rapidly cooling gas or external photoionization can invalidate this equilibrium assumption.
The solar corona is a hot, dilute plasma in which collisional ionization mainly proceeds by direct electron impact and by excitation-autoionization. The principal inverse ion-changing processes are radiative recombination and dielectronic recombination. Three-body recombination is negligible at ordinary coronal densities. Photoionization, stimulated recombination and charge exchange are normally secondary in this fully ionized, collision-dominated regime, although they should be reconsidered near irradiated or neutral interfaces.
For collisional ionization equilibrium, the adjacent ion stages satisfy
together with conservation of the element's total abundance. Thus the low-density ion fraction depends primarily on . This fixes the amount of Si X available for line emission, but not the populations of its excited levels: a coronal approximation does not impose a Boltzmann level distribution.
Within an ion, electrons collisionally excite, de-excite and redistribute levels. Their collision rate coefficient for an atomic transition is
The rate per ion is . Maxwellian detailed balance gives, for ,
Strong allowed transitions undergo rapid spontaneous emission; forbidden magnetic-dipole and electric-quadrupole transitions give much smaller decay rates and can leave a metastable atomic level appreciably populated. Radiative cascades feed levels from higher states. Absorption and stimulated emission add rates proportional to the local radiation field through the Einstein coefficients when significant. Proton collisions can also mix closely separated levels separated by fine structure and should be included when their rates matter.
Write . A stationary collisional-radiative model solves
where, in the simplest electron-plus-radiative model, with for an upward transition. Add proton and radiation rates when appropriate. Include enough levels and radiative branches to capture population feeding. The thermal collision data, radiative probabilities and wavelengths then determine each line's optically thin atomic line intensity:
The factor distributes isotropic emission over solid angle. photon intensities omit . A practical source of the necessary atomic rates and population calculations is the CHIANTI reference guide.
Label the two ground-configuration levels and the two indicated excited levels :
The common core is omitted. Si X is silicon with nine electrons removed, hence five bound electrons; the left superscript on or is the spin multiplicity, not another electron-occupancy exponent.
The observed lines correspond to decays and . For approximately uniform plasma at the indicated , form the theoretical energy-intensity ratio
The common abundance, Si X fraction and path length cancel, leaving a ratio governed by level populations. The frequency factor is for the rounded wavelengths. With observed photon counts, first apply the detector calibration and use the corresponding photon-ratio convention.
The electron-density diagnostic from metastable populations can be seen in a reduced rate model. The upper ground level has only a weak radiative decay to . Ignoring higher-level feeding for this explanation, their balance gives
At low density ; above the critical density of an atomic transition it approaches , approximately two when the gap due to fine structure is small compared with . This saturation can occur while the strong upper-level transitions still remain radiatively dominated.
Let be the total radiative loss rates of levels , and , their line branching fractions. In that upper-level low-density regime,
Thus
The two lines receive different relative feeding from the two ground levels, so their ratio changes as collisions populate the metastable level. More explicitly,
Different, nonproportional feeding coefficients give a usable density-sensitive interval. With , this reduced model even permits explicit inversion:
Physical solutions require positive density and a ratio within the allowed calibration range. The actual Si X calibration should use the full multilevel population system, not discard cross-excitation, cascades or all competing branches simply to use this illustrative formula.
Operationally, evaluate the model ratio at the known temperature over a density grid and match the calibrated observed ratio to it, propagating line-flux and atomic-data uncertainties. No observed ratio or numerical atomic rate set is supplied here, so the result is this inference procedure rather than a numerical density. At low- or high-density saturation, the ratio constrains density poorly. Deconvolve blends or sum every unresolved transition in the theoretical feature; for example, the neighboring branch can contribute near the 356-Angstrom feature when unresolved. Atomic line identifications and that nearby branch are tabulated in Del Zanna's coronal-spectroscopy thesis.
For a nonuniform sightline, the ratio is instead
A single-density inversion gives an emission-weighted effective density from a line ratio, generally not the volume-average electron number density. The density estimate is meaningful within its temperature, optical-depth, line-blending and spatial-uniformity assumptions. This is the precise interpretation of an average density obtained from the observed pair.
For an optically thin low-density gas in collisional ionization equilibrium, the astrophysical cooling function collects losses per pair of hydrogen nuclei into . Its units are ; the time unit is missing from the printed guide to the vertical-axis labels. For the primordial atomic cooling curve, only hydrogen and helium provide bound-state coolants.
Figure 1.
Schematic primordial cooling function with hydrogen and helium features, metal enrichment and photoionization effects
. Details.
Near , collisional excitation of hydrogen becomes effective and subsequent photon emission removes thermal energy. The excitation rate is exponentially suppressed below the atomic energy threshold. Atomic line cooling and collisional ionization generate a strong hydrogen feature around a few ; helium excitation and ionization produce further structure around . Radiative recombination also contributes. Once the gas is highly ionized, these bound-state losses weaken. At –, thermal bremsstrahlung dominates, with approximately apart from slowly varying factors. The figure is an original qualitative sketch anchored to the two supplied values, not an atomic-rate calculation.
Increasing galactic metallicity adds many ionic transitions, producing metal-line cooling and substantially enhancing cooling, particularly in the – range. Metals also provide low-energy fine-structure transitions below . The change is not a uniform vertical shift: the positions and strengths of features depend on the elemental abundances and ionization state.
Photoionization removes bound electrons even where collisions alone would leave atoms neutral, often suppressing the hydrogen and helium line peaks and changing the metal-ion population. It also adds photoionization heating. The net thermal loss is then , with possible heating-cooling equilibrium near . A single density-independent no longer describes every irradiated cloud: the answer also depends on radiation intensity, spectrum, shielding and density. The illustrative photoionized curve denotes altered cooling alone, not the net loss including heating.
Below , primordial gas can use molecular line cooling, especially rotational and vibrational transitions of molecular hydrogen and, in suitable chemical conditions, hydrogen deuteride. Molecule formation and protection against photodissociation are essential. Enriched gas can additionally use metal-line cooling, other molecular species and dust thermal emission, with gas-dust energy exchange important at high density. At high cosmological redshift, residual free electrons can transfer energy to the cosmic microwave background through Compton cooling by the cosmic microwave background if the gas is hotter than the radiation. Expansion can also cool gas adiabatically, but neither that process nor Compton cooling by the cosmic microwave background is an atomic loss law. The atomic threshold is a cooling-channel limitation, not a universal minimum gas temperature.
For an optically thin, low-density plasma, the astrophysical cooling function packages collisional radiation losses and temperature-dependent ion fractions into the coefficient multiplying . Because the volume loss rate has units , the coefficient has units . The energy-density unit printed in the numerical hint cannot be the unit of this coefficient. Interpret the quoted logarithmic values in the dimensionally consistent cooling-coefficient unit.
Assume collisional ionization equilibrium, a primordial hydrogen-helium mixture and no external photoheating. In the specified range, the primordial atomic cooling curve has the following features. Just above , thermal electrons begin to excite neutral hydrogen efficiently; subsequent line emission, especially Lyman-alpha emission, causes a steep rise. The excitation rate contains a threshold factor of order . Hydrogen line cooling is strong near a few times ; the supplied value at provides a useful low-temperature label.
As hydrogen becomes ionized, neutral-hydrogen line cooling declines. Helium excitation and ionization produce a further shoulder or peak around . Collisional excitation, collisional ionization and radiative recombination all contribute: excitation photons remove electron kinetic energy, ionization consumes it, and recombination produces free-bound radiation. Once hydrogen and helium are almost fully stripped, their bound-state cooling disappears and the curve falls into a relatively inefficient interval. At high temperature, thermal bremsstrahlung dominates, with an approximate tail and a weak Gaunt-factor correction.
The sketch uses the two supplied numerical labels and a qualitative hydrogen-helium interpolation; it is not a tabulated atomic-rate calculation.
Figure 1.
Qualitative primordial atomic cooling curve with hydrogen and helium line features, a bremsstrahlung tail, and an illustrative metal-enriched comparison
. Details.
Metal-line cooling raises the cooling coefficient markedly over much of –, because heavier elements supply many ions and excitation transitions after hydrogen and helium have lost their bound electrons. It also broadens and reshapes the line-cooling peaks. Fine-structure lines can permit cooling below the hydrogen atomic threshold; molecular hydrogen can likewise cool metal-free gas below that threshold, but lies outside the requested temperature range. At sufficiently high temperatures thermal free-free emission again dominates the continuum. The enhanced curve in the figure is a schematic comparison, not a numerical claim about a specified metallicity. Thus metals generally shorten the cooling time and extend the temperature range of efficient cooling; ionization state, abundance and radiation field determine the actual curve.