Hydrogen photoionization in a magnetized medium: the rigid-wavefunction approach revisited

Realistic modeling of stellar spectra requires accurate radiative opacity coefficients. Owing to the fragmentary nature of existing data from rigorous quantum-mechanical calculations, photoionization coefficients based on the rigid-wavefunction approximation remain the only practical option for studies of magnetic white dwarfs. Although variants of this approach have been widely used in spectral analyses for decades, a complete and explicit treatment of degeneracy-level breaking has not previously been presented. In this work, we provide a comprehensive description of this procedure, including explicit expressions for the photoionization probability of individual bound-free transitions as functions of magnetic field strength and radiation polarization. We also evaluate the occupation numbers of bound states in a magnetized gas under ionization equilibrium, enabling the calculation of absolute photoionization opacities. Because high-lying atomic states are strongly perturbed by the magnetic field and ultimately dissolved, substantial modifications of the monochromatic absorption are found even for fields below 10 MG–a regime where fully rigorous quantum calculations are numerically demanding and have not yet been applied. Over a wide range of magnetic field strengths, pronounced dichroic features appear in the hydrogen continuum absorption.

💡 Research Summary

The paper addresses a critical gap in the modeling of magnetic white dwarf (MWD) atmospheres: the lack of comprehensive hydrogen photo‑ionization cross‑section data across the magnetic field strengths that characterize most known MWDs (0.01–1000 MG). While fully quantum‑mechanical calculations exist only for fields above ≈10 MG and for a limited set of transitions, realistic atmosphere models require thousands of bound–free cross sections for many photon energies, polarizations, and field strengths. The authors therefore revisit the rigid‑wavefunction approximation (RWA), originally proposed by Lamb & Sutherland, and provide a complete, explicit formulation that can be applied from the weak‑field regime up to several megagauss.

The theoretical development begins with the electric‑dipole matrix element ⟨nlm| r·e_q |kl′m′⟩. By invoking the Wigner‑Eckart theorem, the authors separate the angular dependence (encoded in Wigner 3j symbols) from the reduced radial matrix element ⟨nl∥r∥kl′⟩. This leads to compact expressions for the partial cross sections (Eqs. 4–10) and introduces two sets of geometrical weights, A_q^lm and B_q^lm, which quantify the contribution of transitions with Δl = +1 and Δl = −1 respectively. Tables 1 and 2 list these weights for a range of orbital quantum numbers, clearly showing how they depend on the magnetic quantum number m and on the photon polarization q (π, σ⁺, σ⁻). The symmetry relations (Eq. 18) guarantee that the total cross section summed over polarizations recovers the field‑free result when B → 0.

To obtain absolute cross sections, the reduced radial matrix elements are expressed through the continuum oscillator‑strength density d f/dE. The authors adopt the analytic continuation method of Menzel & Pekeris, which extends discrete bound‑bound oscillator strengths into the continuum by replacing the principal quantum number n′ with a complex wave number k. Equation 24 and Table 3 provide explicit polynomial forms Q_nl,kl′(k) for transitions up to n = 5, enabling rapid evaluation of d f/dE for any photon energy.

A crucial part of the work is the treatment of ionization equilibrium in a magnetized plasma. The magnetic field shifts and splits both bound and free electron energies, altering the statistical weights of each level. The authors compute occupation numbers n_i by solving a modified Saha equation that includes field‑dependent level energies and accounts for the dissolution of high‑lying states (the “magnetic pressure ionization” effect). This yields a self‑consistent set of bound‑state populations that can be combined with the partial cross sections to produce absolute photo‑ionization opacities.

The results, presented in Section 5, demonstrate several important physical effects. First, even for fields below 10 MG, the high‑n hydrogen states are strongly perturbed and eventually merge into the continuum, leading to a substantial reduction of the monochromatic opacity compared with field‑free calculations. Second, the dichroic nature of the magnetized medium is evident: σ⁺ and σ⁻ polarizations can differ by up to 20 % in cross‑section magnitude, while the π component behaves differently still. Third, the RWA predictions agree well with the limited fully quantum‑mechanical data available (e.g., Zhao & Stancil 2007) after smoothing over the fine resonance structure, which is expected to be washed out by field variations across a stellar surface and by instrumental resolution.

In the discussion, the authors acknowledge that the RWA is formally exact only in the zero‑field limit, yet they argue that its empirical reliability extends into the megagauss regime because the magnetic field primarily modifies the energies, not the spatial overlap of the bound and free wavefunctions in the small region that dominates the dipole matrix element. They also point out that the method is computationally inexpensive, allowing the generation of complete opacity tables for any desired grid of (T, ρ, B) values.

Finally, the paper concludes that the rigid‑wavefunction approach, now fully documented with explicit formulas for geometrical weights, branching fractions, and occupation numbers, provides the only practical route to obtain the extensive set of hydrogen photo‑ionization cross sections required for modern MWD atmosphere modeling. The authors suggest future work to combine RWA‑based opacities with the few available high‑field quantum‑mechanical results to construct hybrid opacity tables that retain the detailed resonance structure where it matters most.