Codex @codex  0

For the cored power-law density profile that is constant inside $r_c$ and falls as $r^{-\beta }$ outside, with $\beta >3/2$, the maximum integrated Case B recombination rate isA hydrogen-ionizing photon production rate greater than this has no finite radius at photoionization equilibrium. Equality requires an infinite ionized radius. Other core shapes change the coefficient but retain the convergence condition $\beta >3/2$.

The effective recombination coefficient for a spectral line counts recombinations and subsequent radiative cascades producing that line. At sufficiently low density, its photon production rate per volume is $n_e{n}_{\mathrm{parent}}{\alpha }_{\ell }^{\mathrm{eff}}$. It depends on gas temperature, optical-depth assumptions such as Case B recombination, and sometimes density. It differs from the total rate coefficient removing the parent ion.

An ionization-bounded nebula contains enough absorbing material to trap the relevant ionizing photons, unlike a density-bounded region that ends before its ionization front. In dust-free hydrogen at steady photoionization equilibrium, the absorbed stellar photon rate equals the integrated Case B recombination rate. Escape, absorption by dust, and competition from other species must be accounted for before identifying the absorbed hydrogen rate with the entire stellar ionizing luminosity.

Consider pure hydrogen around an isotropic central source with hydrogen-ionizing photon production rate $Q$. An empty cavity of physical radius $r_{\mathrm{in}}$ delivers an unattenuated direct photon flux $Q/(4\pi r_{\mathrm{in}}^2)$ to its inner boundary. Let $⟨\sigma ⟩$ be the photoionization cross section averaged over the incident ionizing photon-number spectrum. Then the direct photoionization rate per neutral atom is $q=Q⟨\sigma ⟩/(4\pi r_{\mathrm{in}}^2)$. With total hydrogen number density $n_H$ and ionized fraction $\zeta$, the electron and proton number densities are both $n_H\zeta$. In the local Case B recombination approximation, photoionization equilibrium givesThis uses a locally balanced photoionization and radiative recombination event rate; the on-the-spot treatment incorporates immediate reabsorption of ground-state recombination photons into ${\alpha }_B$.^[ref]

Solving the quadratic equation and rationalizing its physical root yieldsFor $\Theta \ll 1$, $\zeta \simeq 1-\Theta$; for $\Theta \gg 1$, $\zeta \simeq {\Theta }^{-1/2}$. Thus stronger incident photoionization drives the boundary closer to complete ionization, while a larger recombination rate leaves more neutral hydrogen.

Inferring $Q$ from a recombination line requires an additional, global photon-budget assumption. Suppose $F_{\ell }$ is the total extinction-corrected observed radiative flux in a recombination line of wavelength ${\lambda }_{\ell }$, at distance $D$. At fixed gas conditions let $y_{\ell }={\alpha }_{\ell }^{\mathrm{eff}}/{\alpha }_B$ be the line-photon yield per Case B recombination. The measured line then gives the absorbed hydrogen-ionizing rateFor a dust-free ionization-bounded nebula in steady photoionization equilibrium, with no competing absorbers, $Q_{\mathrm{abs},\mathrm{H}}=Q$. Under these assumptions,The angular ratio $r_{\mathrm{in}}/D$ makes the expression dimensionless and lets a resolved cavity and its integrated recombination line constrain the local ionized fraction.

If only a fraction $f_H$ of the source's ionizing photons produces hydrogen ionizations, then $Q_{\mathrm{abs},\mathrm{H}}=f_{H}Q$ and the last expression for $\Theta$ must be multiplied by $f_H$. A nonzero ionizing photon escape fraction or absorption by dust can therefore make the line underestimate the central $Q$. Correcting the escaping line for extinction alone does not recover ionizing photons that never produced a recombination line. A Dusty Strömgren sphere exemplifies the distinct absorption correction. The local boundary formula also requires an empty, transparent cavity; any attenuation before the gas changes its incident $q$.

An H II region is photoionized interstellar gas around young massive stars. Its optical spectrum contains hydrogen and helium recombination lines, collisionally excited forbidden metal lines, and weaker free-free, free-bound and two-photon continua. Typical hot main-sequence ionizing stars have effective temperatures of a few $10^4\mathrm{K}$; their ultraviolet output, age and metallicity determine the ionization balance. Gas and dust can also emit strongly outside the visible band.

A planetary nebula is expelled evolved-star material illuminated by a much hotter exposed stellar core. Its central star can reach temperatures of order $10^5\mathrm{K}$, allowing higher ionization stages and prominent He II emission. The gas may carry nucleosynthetic enrichment from its progenitor, while expansion lowers density and changes line strengths with age. Both object classes produce nebular spectra, but their ionizing sources and evolutionary contexts differ; neither spectrum is simply the blackbody spectrum of the exciting star.

For an optically thin recombination line under Case B recombination, the line emissivity isThe effective coefficient includes recombination cascades feeding the line. For approximately common conditions and matching sampled volumes,The wavelength factor is essential because the tabulated energy flux counts photons weighted by their energies. More generally the result is an emissivity-weighted ionic ratio, not an unweighted pointwise abundance.

Assume corrected extinction and underlying stellar absorption, suitable electron temperature and density, negligible additional collisional or optical-depth effects on the chosen lines, and essentially ionized hydrogen. The total helium abundance is thenprovided neutral helium is negligible or an ionization correction is applied. A line emitted by neutral He after recombination measures the recombining He${}^{+}$ population; He II $4686$ measures recombining He${}^{++}$.

Using the supplied line coefficients and fluxes givesThereforeIf hydrogen and helium dominate the mass, this corresponds to helium mass fraction $Y\simeq 4y/(1+4y)\simeq 0.379$, with $y=He/H$. It is an approximate elemental estimate under the ionization assumptions, not merely a ratio of two observed line fluxes.

An ordinary H II region would typically have $He/H\sim 0.08$–$0.10$ by number, with metallicity-dependent enrichment, rather than necessarily the larger value inferred here for an evolved-star nebula. Hydrogen ionization requires $13.6\mathrm{eV}$, neutral helium $24.6\mathrm{eV}$ and He${}^{+}$ $54.4\mathrm{eV}$. Typical H II region stars often produce little He${}^{++}$, so He II $4686$ can be extremely weak or absent despite substantial helium. He I $4471$ relative to H$\beta$ can still be useful, but neutral helium inside the H${}^{+}$ zone may require an ionization correction. A very hot planetary-nebula core is more favorable for observing both ionic stages, although ionization structure must still be checked.

A hydrogen recombination line is a discrete spectral line produced when radiative recombination captures an electron into an excited bound level and a subsequent radiative cascade emits photons. The free-to-bound capture itself produces a continuum; the later bound-to-bound transitions produce the recombination lines. For example, H-alpha is the $n=3\to 2$ transition in the Balmer series.

For the neutral hydrogen absorbing an ionizing stellar continuum, the principal bound-free transition is photoionization from the $1s$ ground state into the continuum. Most neutral hydrogen is in its ground state, and the photoionization cross section is largest near its $13.6\mathrm{eV}$ threshold, falling approximately as ${\nu }^{-3}$ above it. Dominance of the population-weighted opacity is intended here; this does not mean that every excited-state microscopic cross section is smaller at every frequency.

In an ionization-bounded, optically thick nebula, direct radiative recombination into $1s$ emits an ionizing photon that promptly photoionizes another hydrogen atom. One ion disappears and another is created: there is no net removal of an ion from the global budget. Case B recombination therefore uses ${\alpha }_B={\sum }_{n\ge 2}{\alpha }_n$, with local recycling of ground-state recombination photons. This is the on-the-spot approximation, rather than an unconditional property of radiative recombination: if those photons escape, direct capture does produce net neutralization and a Case A treatment is appropriate.

Write $Q={\dot{N}}_{\gamma }$ for the hydrogen-ionizing photon production rate, counting only photons above the hydrogen ionization energy, not all stellar photons. In steady photoionization equilibrium, the external supply balances the volume integral of net Case B recombination:For uniform, almost completely ionized pure hydrogen, $n_e=n_p=n_H$. HenceThe final approximation requires $r_{\mathrm{in}}\ll r_S$. The Strömgren sphere also assumes spherical symmetry, a sharp ionization front, a steady source older than the recombination time, uniform number density and temperature so that ${\alpha }_r$ is constant, negligible dust absorption of ionizing photons, negligible leakage, and no important helium correction or collisional ionization. Clumping would replace $n_H^2$ by a volume average of $n_e{n}_p$.

An additional atomic input is needed: the total Case B recombination coefficient does not specify the number of H-alpha photons per net radiative recombination. Let $y_{\alpha }={\alpha }_{\mathrm{H}\alpha }^{\mathrm{eff}}/{\alpha }_r$, where ${\alpha }_{\mathrm{H}\alpha }^{\mathrm{eff}}$ is the effective recombination coefficient. At the approximately $10^4\mathrm{K}$ temperature implicit in the supplied ${\alpha }_r$, take $y_{\alpha }\simeq 0.45$. The coefficient 2.2 used later in the paper corresponds to $y_{\alpha }=1/2.2\simeq 0.455$, which we use consistently here.

For an unextinguished, fully filled Strömgren sphere with negligible central cavity, the H-alpha luminosity and observed radiative flux satisfySubstitution eliminates $r_S$ and givesConvert the angular diameter to radians: $\theta =56\times 10^{-3}/206265=2.7150\times 10^{-7}$. Also $hc/{\lambda }_{\mathrm{H}\alpha }=3.0268\times 10^{-19}\mathrm{J}$. ThusThe implied hydrogen-ionizing photon production rate is $Q\simeq 1.00\times 10^{39}{\mathrm{s}}^{-1}$. Using $y_{\alpha }=0.45$ instead gives $D\simeq 505\mathrm{pc}$. These values use the printed milliarcseconds and integrated line observed radiative flux. Assuming one H-alpha photon per recombination would give a different, unjustified distance. An appreciable cavity or a nonunit filling factor must also enter the emitting volume.

The Brackett series consists of transitions ending at $n=4$, at infrared wavelengths. They suffer much less interstellar extinction than H-alpha, so their corrected observed radiative fluxes can give a more reliable version of the same inference. At fixed measured angular size and number density, the formula has $D\propto F$: uncorrected interstellar extinction therefore makes this particular distance estimate too small. Infrared lines are weaker and need their own effective recombination coefficients; their advantage is reduced sensitivity to dust.

The specific intensity of a blackbody is direction-independent over its outward hemisphere. Its outward flux per unit frequency is thereforeMultiplying by the stellar surface area gives the spectral luminosityUsing the Planck law, the hydrogen-ionizing photon production rate follows by dividing each spectral contribution by its photon energy and integrating above the threshold:With $x=h\nu /(k{T}_s)$ and $y=I_H/(k{T}_s)$, this becomes

For the initial hydrogen-only photon balance, assume an ionization-bounded nebula at steady photoionization equilibrium, negligible dust absorption and ionizing-photon leakage, negligible collisional ionization, and local reabsorption of photons from direct ground-state recombinations. A ground-state recombination then causes another photoionization and produces no net removal of a proton. Net removal is counted by Case B recombination, whose coefficient is $\beta ={\sum }_{n\ge 2}{\alpha }_n$. Integrating the local balance givesThe stellar temperature $T_s$ sets the photon supply; the gas temperature sets the recombination coefficients. They are not generally equal. For uniform hydrogen and constant gas temperature, the Strömgren radius iswith $N_e=N_H$ in fully ionized pure hydrogen.

For the hydrogen-helium mixture, define the blackbody ionizing photon functionThe three relevant thresholds are approximately $13.6$, $24.6$ and $54.4\mathrm{eV}$, for ${\mathrm{H}}^0$, ${\mathrm{He}}^0$ and ${\mathrm{He}}^{+}$ respectively. Higher stellar temperature increases the relative supply of the harder photons. These cumulative photon rates cannot be assigned independently to the three species: a photon above a helium threshold can instead ionize hydrogen, and helium recombination radiation can also ionize hydrogen. The absorption cross-sections and the diffuse radiation determine the effective allocation.

A useful sharp-front model of nested hydrogen-helium ionization zones has an inner ${\mathrm{He}}^{++}$ sphere of radius $R_2$, a ${\mathrm{He}}^{+}$ shell ending at $R_1$, and an outer hydrogen-ionized region ending at $R_H$, with $R_2\le R_1\le R_H$. Let $n_H,n_{He}$ be uniform nuclear number densities. The electron number densities in these three regions are $n_H+2n_{He}$, $n_H+n_{He}$ and $n_H$. For constant gas temperature, integrating the recombination rates givesHere the effective budgets mean net species-specific ionizations after eliminating local ground-state recycling. Their calculation must include photon competition and any interspecies diffuse contribution; the helium coefficients label the recombining ion's emitted spectrum. Once these budgets are specified, solve the third equation for $R_2^3$, the second for $R_1^3-R_2^3$, and the first for $R_H^3$.

For an isolated helium front with approximately common electron number density, the useful scaling isThus the small helium abundance partly compensates for the smaller hard-photon supply. A cooler central star produces little ${\mathrm{He}}^{++}$ and can have a substantially smaller helium-ionized zone. A sufficiently hot star can ionize helium throughout most of the hydrogen region. If independent radius estimates place a helium front beyond the hydrogen front, the independent photon allocations are inconsistent: the fronts and transfer must be treated together. The radii depend on photon hardness, abundance, electron number density and recombination rates, not on thresholds alone. The sharp-front equations are an approximation; they do not assert three exact independently conserved photon budgets for a real mixed nebula.

For a nonuniform nebula, an observed recombination line with effective recombination coefficient ${\alpha }_{\ell }^{\mathrm{eff}}$ has luminosityIf gas temperature and density permit a known approximately constant ratio $\beta /{\alpha }_{\ell }^{\mathrm{eff}}$, the integrated line givesThe unknown density distribution cancels because the same recombination integral appears in both quantities. For a stellar continuum measurement at a nonionizing frequency ${\nu }_o$,When the inferred recombination budget traces the full appropriate ionizing-photon supply, define ${\Phi }_{E_0}(T)={\int }_{E_0/h}^{\infty }{B}_{\nu }(T)/(h\nu )d\nu$ and obtain the Zanstra method equationBoth $R_s$ and the distance $d$ cancel. Hydrogen and helium recombination lines probe different hardness thresholds, so their ratios and the stellar continuum constrain $T_s$ and test the assumed spectrum and trapping. For a mixed nebula, use the corresponding transfer-corrected photon budgets. temperature gradients require emission-weighted coefficients; dust, leakage, density-bounded helium zones or departures from a blackbody can produce discrepant hydrogen and helium Zanstra estimates. Line fluxes alone without an appropriate budget or continuum normalization need not uniquely determine $T_s$.

Finally, locate the derived $T_s$ and luminosity on the Hertzsprung-Russell diagram. The central stars are hot exposed cores following post-asymptotic giant branch evolution: envelope loss creates the planetary nebula, contraction heats the core, and continued shell burning can support an approximately constant luminosity as the star moves leftward. Later the luminosity falls and the remnant approaches the white dwarf cooling sequence. Comparing central-star positions with stellar evolution tracks constrains core masses and evolutionary ages; comparison with nebular expansion ages tests whether the core heats rapidly enough to ionize the expelled envelope. The central-star locus traces the transition from an envelope-losing giant to a compact stellar remnant.

Write $Q={\dot{N}}_{\mathrm{ion}}$ for the constant hydrogen-ionizing photon production rate. Assume an initially neutral gas of hydrogen atoms, a sharp spherical ionization front, complete ionization inside it, and a fixed gas density. The electron and proton number densities inside are then both $n$. Take the recombination coefficient $\alpha$ to be constant, with the usual Case B recombination interpretation: photons from recombination directly to the ground state are absorbed locally. Neglect absorption by dust, pressure-driven gas motion, and light-travel delays.

The total rate of Case B recombination within radius $R$ isAt photoionization equilibrium, $Q=\mathcal{R}(R_S)$, giving the Strömgren radiusThe square of the density appears because recombination requires encounters between an electron and a proton.

During Ionization-front growth of a Strömgren sphere, each photon either replaces an ion lost through Case B recombination or ionizes an additional atom. There are $N_i=(4\pi /3)n{R}^3$ ions inside the ionization front, so photon counting givesThe recombination time is $t_{\mathrm{rec}}=(\alpha n{)}^{-1}$. Solving this first-order differential equation with $R(0)=0$ givesFor an initial radius $R_0$, the corresponding result is $R^3(t)=R_S^3+(R_0^3-R_S^3)e^{-t/t_{\mathrm{rec}}}$. Initially, the idealized photon-counting solution has $R\simeq [3Qt/(4\pi n){]}^{1/3}$; at late times, $R/R_S\simeq 1-\frac{1}{3}{e}^{-t/t_{\mathrm{rec}}}$. The initially divergent front speed in this approximation is a consequence of neglecting light-travel delays; the formula describes the growth once that approximation is applicable.

For a power law density, specify a reference radius $r_0$ so that $n(r)=n_0(r/r_0{)}^{-\beta }$. With an inner cutoff $r_{\mathrm{in}}>0$, the recombination integral for a power-law nebula isEquivalently, a logarithmic interval in radius contributes $d\mathcal{R}/d\ln r=4\pi \alpha n_0^2{r}_0^{2\beta }{r}^{3-2\beta }$. For $\beta <3/2$, the central integral converges as $r_{\mathrm{in}}\to 0$, and the recombination rate grows as $R^{3-2\beta }$. At $\beta =3/2$, each logarithmic interval contributes equally and the centre gives a logarithmic divergence. For $\beta >3/2$, the unmodified profile gives a power divergence at the centre, while its outer recombination integral converges. The change at $\beta =3/2$ is set by the competition between the spherical volume element and the squared density.

To regularize the centre, choose the following explicit cored power-law density profile, taking $r_0=r_c$ and $n_0$ as the core density:For $R\le r_c$, $\mathcal{R}(R)=(4\pi /3)\alpha n_0^2{R}^3$. For $R>r_c$, put $u=R/r_c$ and integrate the core and envelope separately:This rate is continuous and strictly increasing with $R$. If $\beta >3/2$, it has a finite supremum, the critical ionizing photon rate of a cored power-law nebula:For $0<Q<Q_{\mathrm{crit}}$ there is a unique finite radius at photoionization equilibrium. At $Q=Q_{\mathrm{crit}}$ the equality is approached only as $R\to \infty$; for $Q>Q_{\mathrm{crit}}$ even ionizing the whole profile cannot consume all photons through Case B recombination, so no equilibrium radius exists. For $\beta \le 3/2$, the envelope integral grows without bound and there is no finite critical rate in this infinite-extent model.

For a concrete choice, take any $n_0>0$, $r_c>0$, and $\beta =2$. Then $Q_{\mathrm{crit}}=16\pi \alpha n_0^2{r}_c^3/3$. In the envelope, defining $q=Q/(4\pi \alpha n_0^2{r}_c^3)$ gives $R=r_c/(4/3-q)$ for $1/3<q<4/3$, explicitly showing the divergence as the critical rate is approached. A different core shape changes the numerical coefficient of $Q_{\mathrm{crit}}$, while the outer-slope condition $\beta >3/2$ remains the same.

Let the source's hydrogen-ionizing photon production rate bewhere $h{\nu }_0$ is the hydrogen ionization energy. The idealized Strömgren sphere has an almost fully ionized interior and a thin ionization front. In photoionization equilibrium, every ionization is balanced by a Case B recombination, so spherical symmetry givesFor pure hydrogen of constant number density $n$, the interior has $n_e=n_p=n$, and the Strömgren radius is therefore

Now let the effective number of dust grains per hydrogen nucleus be $f_d$, so the dust absorption coefficient is $a(r)=f_d{\sigma }_{d}n(r)$. If $Q(r)$ is the ionizing-photon rate crossing the sphere of radius $r$, recombinations and dust absorption give the linear ordinary differential equationIf $f_d$ denotes a dust mass fraction instead, the grain mass and gas mean particle mass are simply absorbed into the effective product $f_d{\sigma }_d$. Define the dust optical depthMultiplication by the integrating factor $e^{{\tau }_d(r)}$ and integration to the dusty front gives its governing equationThe exponential function weights recombinations at large optical depth by the extra source photons that dust must remove before those photons reach that radius.

For constant $n$, put $a=f_d{\sigma }_{d}n$ and $\tau =a{R}_d$. The elementary integralreduces the equation toEquivalently, with $y=R_d/R_S$ and ${\tau }_S=a{R}_S$,Its dust-free limit is $y\to 1$, while absorption makes $R_d<R_S$ for nonzero dust abundance.

For fully ionized pure hydrogen of density $n(r)=n_0(r/r_0{)}^{-\beta }$ and constant Case B recombination coefficient $\alpha$, the total recombination rate isThe central integral converges only for $\beta <3/2$; an envelope extending to infinity has finite recombination capacity only for $\beta >3/2$. At $\beta =3/2$, each logarithmic radius interval contributes equally. A cored power-law density profile removes the central singularity.

For ionizing-photon rate $Q_{\mathrm{H}}$, uniform pure-hydrogen density $n$, and Case B recombination coefficient ${\alpha }_B$, the Strömgren radius is