= Solution
The <Boltzmann distribution> for two atomic levels at thermodynamic equilibrium is
$$
\boxed{\frac{N_a}{N_b}=\frac{\omega_a}{\omega_b}
\exp\left[-\frac{E_a-E_b}{kT}\right],}
$$
where $\omega$ is the <statistical weight of an atomic level>. For a resolved level of total angular momentum $J$, $\omega=2J+1$ when magnetic substates are unresolved.
To extend this counting to the continuum, assume an ideal, nondegenerate <plasma> with classical free <Electrons> and neglect the small difference between the translational masses of the two heavy ionic species. Let $z_e=e^{\mu_e/kT}$ be the <electron fugacity in the Saha equation>, with the <Electron> <chemical potential> measured relative to its rest energy. The two <Electron> spin states and the phase-space cell volume $h^3$ give
$$
N_e=\frac{2z_e}{h^3}\int_{\mathbb R^3}e^{-p^2/(2mkT)}\,d^3p
=2z_e\left(\frac{2\pi mkT}{h^2}\right)^{3/2}.
$$
Chemical equilibrium for capture equates the bound species' <chemical potential> to the sum for its parent <ion> and free <Electron>. Thus applying <Boltzmann factors> to a bound level at energy $-I_n$ gives $N_n/N_+=(\omega_n/\omega_+)z_e e^{I_n/kT}$. Eliminating the <Electron> fugacity proves the level-specific <Saha ionization equation>:
$$
\boxed{\frac{N_n}{N_eN_+}
=\frac{\omega_n}{2\omega_+}
\left(\frac{h^2}{2\pi mkT}\right)^{3/2}e^{I_n/kT}.}
$$
The factor two is the <Electron> spin multiplicity; the factor with power $3/2$ is the inverse free-electron translational state density. A total-ionization Saha equation sums the bound states into an internal partition function, whereas this expression refers to one specified level.
For <dielectronic recombination>, write the entrance <ion> as $X^{(z+1)+}$ in ground state $i$, and the intermediate state of $X^{z+}$ as $d=(j,n\ell)$. Capture excites the core from $i$ to $j$ while binding the incident <Electron>. The resonance energy relative to the entrance continuum is $\bar E=E_{ij}-I_{n\ell}>0$. The doubly excited state can eject the <Electron> by <autoionization>, or undergo <radiative stabilization> to a bound state.
Let $A^a_{di}$ be the partial <autoionization> rate back to the entrance <ion>, $A_a$ the total <autoionization> rate, $A_r$ the total radiative decay rate, and $A_s$ the rate of radiative decays that stabilize into bound levels. For an isolated narrow resonance, the Saha-Boltzmann equilibrium ratio has the bound-state energy replaced by the positive resonance energy:
$$
\frac{N_d^{\rm eq}}{N_eN_i}
=\frac{g_d}{2g_i}\left(\frac{h^2}{2\pi mkT}\right)^{3/2}e^{-\bar E/kT}.
$$
At equilibrium, inverse <dielectronic capture> and <autoionization> fluxes balance: $N_eN_iC_d=N_d^{\rm eq}A^a_{di}$. This establishes the capture coefficient even when the actual <plasma> is not in local thermodynamic equilibrium. In the low-density, weak-radiation recombining <plasma>, the intermediate level instead has the stationary balance
$$
N_eN_iC_d=N_d(A_a+A_r).
$$
Multiplying its population by the successful stabilization rate gives the actual recombination coefficient
$$
\boxed{\alpha_d(d)=\frac{g_d}{2g_i}
\left(\frac{h^2}{2\pi mkT}\right)^{3/2}
e^{-\bar E/kT}\frac{A^a_{di}A_s}{A_a+A_r}.}
$$
Here the equilibrium population was used only to determine the capture rate by detailed balance; it was not assumed to equal the intermediate population when <photons> escape. This is precisely why the competing decay probabilities appear.
Use the supplied core-transition relation between the <Einstein coefficients> and <atomic oscillator strength>,
$$
A_{ji}=\frac{\alpha^4c}{2a_0}\left(\frac{E_{ij}}{I_H}\right)^2\frac{g_i}{g_j}f_{ij},
\qquad \frac{h^2}{2\pi m}=4\pi a_0^2I_H.
$$
Here $\alpha$ is the <fine-structure constant> and $a_0$ the <Bohr radius>, not a recombination coefficient. Factor out $A_{ji}$ and define
$$
\boxed{\beta_{j,n\ell}=\frac{g_d}{2g_j}\,
\frac{A^a_{di}A_s}{A_{ji}(A_a+A_r)}.}
$$
Substitution yields
$$
\begin{aligned}
\alpha_d(j,n\ell)
&=4\pi^{3/2}\alpha^4a_0^2c\,
\frac{E_{ij}^2}{I_H^{1/2}(kT)^{3/2}}
e^{-\bar E/kT}f_{ij}\beta_{j,n\ell}\\
&=\boxed{4\pi^{3/2}\alpha^4a_0^2c
\left(\frac{E_{ij}}{I_H}\right)^{1/2}
\left(\frac{E_{ij}}{kT}\right)^{3/2}
e^{-\bar E/kT}f_{ij}\beta_{j,n\ell}.}
\end{aligned}
$$
This definition accommodates several entrance and decay channels. In the common spectator-electron approximation with one <autoionization> channel and stabilization by the core transition, $A_s=A_r=A_{ji}$ and $A^a_{di}=A_a$. If all coupled substates of the captured $n\ell$ <Electron> are summed, $g_d=2(2\ell+1)g_j$, giving
$$
\boxed{\beta_{j,n\ell}=(2\ell+1)\frac{A_a}{A_a+A_{ji}}.}
$$
For an intermediate state resolved by <fine structure>, retain its actual $g_d$ instead of that summed weight. The stabilization probability itself is $A_s/(A_a+A_r)$; \b[the requested $\beta_{j,n\ell}$ is not just this branching probability, because the prefactor has already factored out the core radiative rate and <statistical weights of atomic levels>.] Radiative transitions to another still-autoionizing state require that state's eventual survival probability, rather than automatically counting every emitted <photon> as stabilization.
The total <dielectronic recombination> coefficient sums all accessible positive-energy resonances. In this isolated-resonance approximation it has the form
$$
\alpha_d(T)=T^{-3/2}\sum_d c_d e^{-\bar E_d/kT},\qquad c_d\geq0,
$$
with atomic factors included in $c_d$. The <resonance-temperature dependence of dielectronic recombination> follows by differentiating a single term's logarithm:
$$
\frac{d}{dT}\log\alpha_d(d)=-\frac{3}{2T}+\frac{\bar E_d}{kT^2},
\qquad \boxed{kT_{\rm peak}=\frac23\bar E_d.}
$$
Thus low temperatures suppress a positive-energy resonance, although resonances very near threshold can still matter; at sufficiently high <temperature> a fixed set of resonances falls as $T^{-3/2}$. Several core-excitation series can produce several peaks. When <autoionization> is much faster than radiation, the capture-and-survival product is limited by the stabilizing radiative rate; when capture is weak, <autoionization> is the limiting inverse-capture rate. For a fixed Rydberg angular channel, <autoionization> typically decreases at high principal quantum number while the core radiative rate changes little. Capture into very high levels is vulnerable to subsequent collisional or field ionization, so finite-density survival can suppress the total rate. \b[<dielectronic recombination> is resonant and strongly temperature-selective; its magnitude can be important in coronal ionization balance.]
Back to article page
