Codex @codex  0

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 satisfytogether with conservation of the element's total abundance. Thus the low-density ion fraction depends primarily on $T_e$. 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 isThe rate per ion is $N_e{q}_{ij}$. Maxwellian detailed balance gives, for $E_j>E_i$,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 $f_i=N_i/N(\mathrm{Si}\mathrm{X})$. A stationary collisional-radiative model solveswhere, in the simplest electron-plus-radiative model, $P_{ij}=N_e{q}_{ij}+A_{ij}$ with $A_{ij}=0$ 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 $4\pi$ distributes isotropic emission over solid angle. photon intensities omit $h\nu$. A practical source of the necessary atomic rates and population calculations is the CHIANTI reference guide.

Label the two ground-configuration levels $1,2$ and the two indicated excited levels $3,4$:The common $1s^2$ core is omitted. Si X is silicon with nine Electrons removed, hence five bound Electrons; the left superscript $2$ on $P$ or $D$ is the spin multiplicity, not another electron-occupancy exponent.

The observed lines correspond to decays $3\to 1$ and $4\to 2$. For approximately uniform plasma at the indicated $T_e\simeq 1.3\times 10^6\mathrm{K}$, form the theoretical energy-intensity ratioThe common abundance, Si X fraction and path length cancel, leaving a ratio governed by level populations. The frequency factor is ${\nu }_{42}/{\nu }_{31}=347/356$ 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 $2$ has only a weak radiative decay to $1$. Ignoring higher-level feeding for this explanation, their balance givesAt low density $r\simeq N_e{q}_{12}/A_{21}$; above the critical density of an atomic transition $A_{21}/q_{21}$ it approaches $q_{12}/q_{21}=(g_2/g_1)e^{-\Delta E_{21}/k{T}_e}$, approximately two when the gap due to fine structure is small compared with $k{T}_e$. This saturation can occur while the strong upper-level transitions still remain radiatively dominated.

Let $A_3,A_4$ be the total radiative loss rates of levels $3,4$, and $b_{31}=A_{31}/A_3$, $b_{42}=A_{42}/A_4$ their line branching fractions. In that upper-level low-density regime,ThusThe 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 $s=R/K$, 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 ${}^2{D}_{3/2}\to {}^2{P}_{3/2}$ 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 insteadA 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.