Codex @codex  0

The effective area converting incident photon flux into a stimulated emission rate for a specified transition, frequency, polarization and sublevel distribution. For statistically populated energy levels of statistical weights $g_1,g_2$, the Einstein coefficients imply $g_1{\sigma }_{12}=g_2{\sigma }_{21}$. The related absorption cross section therefore differs when the statistical weights differ.

An absorption equivalent width is $W_{\lambda }=\int [1-I_{\lambda }/I_{\lambda }^0]d\lambda$, the wavelength width of a completely black interval that removes the same continuum-normalized area. For a narrow spectral line profile, small optical depth and negligible stimulated emission, the energy-density convention for Einstein coefficients gives $W_{\lambda }=h{B}_{12}{N}_1/{\nu }_{12}$. This follows by replacing $1-e^{-{\tau }_{\nu }}$ by ${\tau }_{\nu }$ and using $|d\lambda /d\nu |=c/{\nu }^2$. At redshift $z$, the observed and rest-frame values satisfy $W_{\mathrm{obs}}=(1+z)W_{\mathrm{rest}}$.

Let $A_{21}$ be the spontaneous emission coefficient, and let $B_{12}\rho$ and $B_{21}\rho$ be the absorption and stimulated emission rates per atom under the radiation spectral energy density $\rho$ at resonance. Because the two levels are nondegenerate, the Einstein coefficients satisfy $B_{12}=B_{21}=B$.

The radiative balance in a nondegenerate two-level atom isThe total spontaneous and stimulated emission rates are $R_{\mathrm{sp}}=N_2{A}_{21}$ and $R_{\mathrm{st}}=N_{2}B\rho$. Dividing the balance equation by $N_{2}B\rho$ givesThis presumes nonzero drive and upper-level population, and no other population-transfer mechanism. It does not require a thermal Boltzmann distribution. In this closed driven two-level system, spontaneous decay keeps $N_1>N_2$; arbitrarily intense driving approaches equal populations, rather than producing a population inversion.

For an optically thin collisionally excited line, the volume emissivity is approximately $j_{ul}=n_{e}n(X^i)q_{lu}(T_e)h{\nu }_{ul}$ below the transition's critical density. The excitation coefficient includes a factor $T_e^{-1/2}\exp [-\Delta E/(k_B{T}_e)]$, so temperature is needed to convert luminosity into an ionic abundance. Density is needed both to assess collisional quenching of a metastable line and to convert an emission measure into a gas mass. For a uniform recombining ion, for example, $M_i=m_i{L}_{\mathrm{line}}/[n_e{\alpha }_{\mathrm{eff}}h\nu ]$. Temperature structure, density clumping, extinction and unobserved ionization stages introduce further corrections.

A useful density doublet has closely spaced upper energies but different spontaneous decay rates and hence different critical densities. Nearly equal excitation energies reduce temperature dependence of the ratio; similar wavelengths reduce differential extinction and relative calibration errors. Metastable upper levels have sufficiently small Einstein coefficients for collisional de-excitation to compete at nebular densities.

Use statistical weights $g_2=6$, $g_3=4$. The source's O II ground-state term has a typographical error: it is $2s^{2}2p^3{}^4{S}_{3/2}^{\mathrm{o}}$, not ${}^2{S}_{3/2}$. A ${}^{2}S$ term has $J=1/2$, so the printed term is internally inconsistent. The tabulated $J$ values, wavelengths and rates can still be used.^[ref]

At high density, collisions establish a Boltzmann distribution between the upper states. Their energy splitting is small compared with $k_B{T}_e$, so $n_2/n_3\simeq g_2/g_3$. If both lines escape without absorption,The exact ratio includes $\exp [-(E_2-E_3)/(k_B{T}_e)]$, an approximately $0.2\%$ correction at $10^4\mathrm{K}$.

At low density, assume excitation from the common ground state and radiative decay back to it dominate, with no substantial collisional transfer between upper states. Statistical equilibrium gives $n_j{A}_{j1}=n_1{n}_e{q}_{1j}$. Detailed balance relates excitation to the supplied de-excitation coefficient, yieldingEach excitation then produces one line photon, so the spontaneous rate cancels:At high density the upper populations are collisionally thermalized and the ratio retains the radiative lifetimes; at low density radiative escape is nearly certain and the ratio is set by the excitation probabilities instead.

At $T_e=10^4\mathrm{K}$, the two de-excitation coefficients are $C_{21}=1.161\times 10^{-14}$ and $C_{31}=1.1825\times 10^{-14}{\mathrm{m}}^3{\mathrm{s}}^{-1}$. ThusThe useful density range lies between and around these scales. In the independent-upper-level approximation,The density sensitivity of a metastable doublet depends on the error measure. The fractional response $|d\ln R/d\ln n_e|$ is maximal atIf absolute ratio errors are constant, maximizing $|dR/d\ln n_e|$ instead gives $n_e=n_{c,2}$. An unspecified experimental noise model therefore does not define a unique sharper optimum than this transition-density range.

For Ne IV, the same calculation gives critical densities $4.23\times 10^{10}$ and $3.93\times 10^{11}{\mathrm{m}}^{-3}$, with fractional sensitivity peaking near $1.29\times 10^{11}{\mathrm{m}}^{-3}$. This is a useful higher-density diagnostic when O II is saturated. However, Ne${}^{3+}$ and O${}^{+}$ need not occupy the same ionization zone, so their inferred densities need not describe the same gas. The ultraviolet Ne IV pair also needs suitable observations and extinction correction. The proposal is useful with these qualifications, not a guaranteed replacement for the density of the O II emitting zone.

Use the Einstein coefficients convention in which $B_{ij}$ multiplies spectral radiation energy density $u_{\nu }$, with $u_{\nu }=c^{-1}\int I_{\nu }d\Omega$. A normalized spectral line profile obeys $\int \phi (\nu )d\nu =1$. The spontaneous emission probability per upper-state atom per unit time is $A_{21}$; the stimulated emission and upward absorption rates per atom are respectively $B_{21}\int u_{\nu }\phi (\nu )d\nu$ and $B_{12}\int u_{\nu }\phi (\nu )d\nu$.

For isotropic spontaneous emission, $n_2{A}_{21}$ transitions occur per unit volume and time, each supplying photon energy $h{\nu }_{12}$. Dividing among $4\pi$ steradians givesFor a small incident beam, its contribution to $u_{\nu }$ is $I_{\nu }d\Omega /c$. Absorption removes and stimulated emission adds photons in the incident direction. Matching the net loss to ${\alpha }_{\nu }{I}_{\nu }d\Omega$ therefore givesReplacing ${\nu }_{12}$ by $\nu$ across the narrow spectral line profile gives the printed expressions. The negative sign is essential: an inverted population can amplify the incident specific intensity. Thermal equilibrium and the Planck function imply $g_1{B}_{12}=g_2{B}_{21}$ and $A_{21}/B_{21}=8\pi h{\nu }_{12}^3/c^3$, where $g_i$ is the state degeneracy. A convention using mean intensity instead changes both $B$ coefficients by $4\pi /c$.

With the spontaneous emission term suppressed, $j_{\nu }=0$. A uniform absorption coefficient over length $\ell$ gives directlyThe normalized specific intensity deficit is unchanged by conversion between wavelength and frequency, since both numerator and denominator acquire the same Jacobian determinant. The Boltzmann distribution and the Einstein coefficients give $n_2{B}_{21}/(n_1{B}_{12})=e^{-h{\nu }_{12}/(k_B{T}_e)}\ll 1$. Consequently, in the optically thin limit,Here $N_1=n_1\ell$ is the column density, and $\nu \simeq {\nu }_{12}$ has been used in slowly varying factors. This linear relation between equivalent width and column density fails once the central optical depth becomes appreciable.

For the unresolved Lyman-alpha line, the lower $1s$ and combined upper $2p$ state degeneracies are $g_1=2$ and $g_2=6$, so the Einstein coefficients giveThe measured equivalent width $W=7.0\times 10^{-12}\mathrm{m}$ then impliesAssume a central transverse line of sight through a uniform filament, so $\ell =d=0.5\mathrm{Mpc}=1.54\times 10^{22}\mathrm{m}$. Nearly all neutral hydrogen is in its ground state at this temperature, henceThe cylinder's length $L$ determines its total volume but does not enter a transverse column density. Geometry must be specified: an axial line of sight has $\ell =L$ and gives $n_{\mathrm{HI}}=4.2\times 10^{-7}{\mathrm{m}}^{-3}$ for the same measured equivalent width. Intermediate chords give intermediate path-dependent estimates. The numbers above treat $W$ as a rest-frame equivalent width; an observed value at redshift $z$ must first be divided by $1+z$.

Specify the radiation normalization first. Let $u_{\nu }$ be radiation energy density per ordinary frequency, with SI units $\mathrm{J}{\mathrm{m}}^{-3}\mathrm{H}{\mathrm{z}}^{-1}$. Then $A_{21}$ is the spontaneous emission rate per upper-state chemical atom, while $B_{12}{u}_{\nu }$ and $B_{21}{u}_{\nu }$ are the induced absorption and stimulated emission rates, evaluated at a sharp transition frequency. For a finite spectral line profile, replace $u_{\nu }$ by its profile-weighted average.

The Einstein coefficients obeyHere $g_i$ are the statistical weights of the two energy levels, and the field is unpolarized and isotropic. The units and physical dimensions areThe apparently time-independent dimension of $B$ follows because a per-hertz density includes an extra unit of time. If instead $B$ multiplies the mean intensity $J_{\nu }$, use $B^{(J)}=(4\pi /c)B$ and $A/B^{(J)}=2h{\nu }_0^3/c^2$. These distinct normalizations must not be mixed in the radiative transfer equation.

Use the energy-density convention of part (a), a normalized common spectral line profile $\int \varphi (\nu )d\nu =1$, and isotropic spontaneous emission. For a narrow spectral line centered at ${\nu }_0$,The factor $1/c$ follows from $u_{\nu }=4\pi J_{\nu }/c$: angle-integrated absorption removes $c{\kappa }_{\nu }{u}_{\nu }$ of energy per volume and time. The emission coefficient is per solid angle, whereas the absorption coefficient has units of inverse length.

The Boltzmann distribution in thermal equilibrium gives $n_2/n_1=(g_2/g_1)e^{-h{\nu }_0/(k_{B}T)}$. Dividing the two expressions and using the relations between Einstein coefficients,Thus the source function equals the Planck function at the transition frequency, and in the usual narrow-line approximationWith an exactly common finite-width $\varphi$ and constant line-center Einstein coefficients, their ratio is constant across the profile; exact equality to a varying $B_{\nu }$ at every frequency is only a narrow-line approximation. Exact thermal radiative transfer requires frequency-resolved detailed balance, which supplies the corresponding small profile corrections. This distinction does not alter the intended equilibrium result.

Rabi oscillations require a coherent drive with a stable phase over several population cycles and an isolated, appreciably coupled pair of quantum states. For large-amplitude oscillations, the Rabi detuning must be comparable to or smaller than the Rabi angular frequency. The population lifetime, optical coherence time, transit time and drive coherence time should all exceed $2\pi /{\Omega }_R$ by enough to resolve the cycles. An ensemble also needs sufficiently uniform Rabi frequencies and Rabi detunings that phase averaging does not wash the oscillations out. The usual sinusoidal formula additionally assumes the rotating-wave approximation.

The Einstein coefficients describe incoherent transition rates, rather than the evolving optical probability amplitudes. Their simplest use is appropriate for broadband, randomly phased radiation, with field correlation time short relative to population evolution, or for rapidly dephased chemical atoms whose optical coherences can be eliminated. The radiation energy density per unit frequency must be smooth over the relevant spectral line. Stimulated upward and downward transitions and spontaneous emission then enter rate equations, without predicting coherent population return. A short observation time can give a small excitation probability even for a coherent field; that condition alone does not justify an Einstein coefficients rate description. Coherent laser irradiation can also admit steady rate equations after optical coherence has relaxed, with the appropriate saturation response.

Define $r=g_2/g_1$ and let $\sigma ={\sigma }_{21}(\omega )$ be the stimulated-emission cross section, averaged over upper magnetic sublevels. The Einstein coefficients relation $g_1{B}_{12}=g_2{B}_{21}$ gives the absorption cross section ${\sigma }_{12}=r\sigma$. Assume equally populated magnetic sublevels within each energy level, a steady pump, and a monochromatic beam. The photon flux is $I/(\hslash \omega )$; define $W=\sigma I/(\hslash \omega )$ and the effective population inversion $D=N_2-r{N}_1$. Then the rate equations areThe $WD$ term includes both stimulated emission and absorption. In steady state, adding the two rate equations gives $N_1=R_2{\tau }_1$. The first equation then yieldsThe directed beam gains energy at rate $\hslash \omega WD$ per volume, giving $dI/dz=\sigma DI$ and the optical gain coefficient $\alpha =\sigma D$. At line centre, write ${\sigma }_0={\sigma }_{21}({\omega }_{21})$. ThereforeHere ${\alpha }_0$ is the small-signal optical gain coefficient and $I_s$ is the saturation intensity, at which the local optical gain is half its small-signal value. Positive optical gain requires $r{A}_{21}{\tau }_1<1$: the lower energy level must empty sufficiently quickly.

The cancellation of ${\tau }_1$ from $I_s$ is specific to this model's steady fixed pump and the fact that every upper spontaneous emission feeds the lower energy level. Adding an extra lower-state accumulation factor to $I_s$ would contradict these rate equations. Alternative conventions using the absorption cross section must replace ${\sigma }_0$ by ${\sigma }_{12}/r$ throughout. Homogeneous spectral broadening permits the whole selected ensemble to share this local saturation response.

The Boltzmann distribution for two atomic levels at thermodynamic equilibrium iswhere $\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$ giveChemical 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: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_{di}^a$ 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:At equilibrium, inverse dielectronic capture and autoionization fluxes balance: $N_e{N}_i{C}_d=N_d^{\mathrm{eq}}{A}_{di}^a$. 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 balanceMultiplying its population by the successful stabilization rate gives the actual recombination coefficientHere 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,Here $\alpha$ is the fine-structure constant and $a_0$ the Bohr radius, not a recombination coefficient. Factor out $A_{ji}$ and defineSubstitution yieldsThis 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_{di}^a=A_a$. If all coupled substates of the captured $n\ell$ electron are summed, $g_d=2(2\ell +1)g_j$, givingFor 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)$; 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 formwith atomic factors included in $c_d$. The resonance-temperature dependence of dielectronic recombination follows by differentiating a single term's logarithm: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. dielectronic recombination is resonant and strongly temperature-selective; its magnitude can be important in coronal ionization balance.

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.

The spectral radiation energy density $u(\omega )$ is defined per unit angular frequency and physical volume. Thermal Planck law convention excludes the vacuum zero-point term. A density per ordinary frequency instead obeys $u_{\omega }=u_{\nu }/(2\pi )$, evaluated at $\nu =\omega /(2\pi )$; Einstein coefficients multiplying these densities must be converted correspondingly.

Let $A_{21}$ be the spontaneous emission coefficient and let $B\rho$ be the induced transition rate per atom under a resonant radiation spectral energy density $\rho$. For nondegenerate levels the Einstein coefficients obey $B_{12}=B_{21}=B$. In a closed two-level system with no collisional or other population transfer, steady state requiresFor nonzero upper population and drive, division gives the ratio of total spontaneous and stimulated emission rates in the title. In particular $N_1>N_2$ when $A_{21}>0$; continuous driving of this closed model approaches equal populations but cannot create a population inversion. No thermal population law is assumed.