Issue |
A&A
Volume 693, January 2025
|
|
---|---|---|
Article Number | L5 | |
Number of page(s) | 4 | |
Section | Letters to the Editor | |
DOI | https://doi.org/10.1051/0004-6361/202452569 | |
Published online | 03 January 2025 |
Letter to the Editor
Strong signature of right-handed circularly polarized photoionization close to the cyclotron line in the atmosphere of magnetic white dwarfs
Instituto de Ciencias Astronómicas, de la Tierra y del Espacio (CONICET-UNSJ), Av. España 1512 (sur), 5400 San Juan, Argentina
⋆ Corresponding author; rene.rohrmann@gmail.com
Received:
10
October
2024
Accepted:
9
December
2024
Magnetic fields break the symmetry of the interaction of atoms with photons with different polarizations, yielding chirality and anisotropy properties. The dependence of the absorption spectrum on the polarization, a phenomenon known as dichroism, is present in the atmosphere of magnetic white dwarfs. Its evaluation for processes in the continuum spectrum has been elusive so far due to the absence of appropriate ionization equilibrium models and incomplete data on photoionization cross sections. We combined rigorous solutions to the equilibrium of atomic populations with approximate cross sections to calculate the absolute opacity due to photoionization in a magnetized hydrogen gas. We predict a strong right-handed circularly polarized absorption (χ+) formed blueward of the cyclotron resonance for fields from about 14 to several hundred megagauss. In energies lower than the cyclotron fundamental, this absorption shows a deep trough with respect to linear and left-handed circular polarizations that steepens with the field strength. The jump in χ+ is due to the confluence of a large number of photoionization continua produced by right-handed circularly polarized transitions from atomic states with a nonnegative magnetic quantum number toward different Landau levels.
Key words: atomic processes / magnetic fields / opacity / stars: atmospheres / white dwarfs
© The Authors 2025
Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
This article is published in open access under the Subscribe to Open model. Subscribe to A&A to support open access publication.
1. Introduction
About 800 magnetic white dwarfs (MWDs) have been discovered so far (Ferrario et al. 2020); most have hydrogen-rich atmospheres (similar to non-MWDs) and magnetic fields (B) from around 0.01 up to nearly 103 megagauss (MG). Analysis of a volume-limited sample suggests that 20% of all white dwarfs have magnetic fields in their surfaces (Bagnulo & Landstreet 2021). Field intensities in MWDs are mainly detected via displacements of spectral lines and Zeeman splittings (Amorim et al. 2023). They can also be detected via narrow- and broadband circular polarization measurements and spectropolarimetry; the latter is a valuable method for measuring fields when it is difficult to identify spectral lines (Berdyugin et al. 2022) or there is a complex line pattern (strong MWDs).
Current model atmospheres fit the spectra of MWDs with a weak or moderate field strength, B ≲ 50 MG, relatively well but have trouble achieving satisfactory precision for stronger ones (Külebi et al. 2009; Hardy et al. 2023; Vera-Rueda & Rohrmann 2024). Also, the description of the continuum polarization spectrum in highly MWDs remains incomplete. This is partly due to the lack of reliable opacity data for bound-free (bf) atomic transitions in strong magnetic fields. Detailed models of both flux and polarization spectra require magnetic-field-dependent photoionization opacities for atomic H. Accurate photoionization cross sections have been obtained since the 1990s (Delande et al. 1991; Wang & Greene 1991; Merani et al. 1995; Zhao & Stancil 2007; Zhao 2021). However, the available data are very limited and unsatisfactory for model atmosphere calculations, as they cover a small number of transitions over a sparse grid of field strengths and a usually incomplete polarization basis. Consequently, MWD models use an adapted form (Jordan 1992) of the bf absorption first formulated by Lamb & Sutherland (1972) for weak fields (B < 10 MG). This approach is justified for the linear Zeeman regime, where the field-perturbed, atomic Hamiltonian is diagonal and so the matrix elements associated with bf radiative transitions are unaltered. Such properties led to an ansatz to evaluating field-dependent cross sections in terms of their zero-field values, which are distributed among the various components according to the Wigner-Eckart theorem. Curiously, to the best of our knowledge, a direct comparison between the cross sections resulting from this procedure and those from full quantum mechanical calculations has never been made. Moreover, no definitive conclusions about bf cross sections can be reached from spectrum fits unless they are self-consistently applied with field-dependent atomic populations. The absolute value of the bf opacity depends on the occupation numbers of many sublevels, which must be calculated in a thermodynamically consistent way for partially ionized, magnetized atmospheres. The lack of chemical equilibrium models for hydrogen gas in the MWD field regime has been a significant limitation of all model atmospheres published until recently. Nowadays, an equation of state fulfills this requirement (Vera-Rueda & Rohrmann 2020).
This Letter presents the first self-consistent, total photoionization absorption of a partially ionized hydrogen gas in the magnetic field regime of MWDs. A comparison of the cross sections of the updated Lamb-Sutherland ansatz with full quantum values is given for the first time (Sect. 2). These calculations were combined with a detailed chemical model (Sect. 3) to obtain the total photoionization opacity as a sum of all relevant transitions arising from the discrete state spectrum of magnetized hydrogen atoms for the three basic polarizations of light (Sect. 4). The prediction of a previously unknown dichroic feature is reported.
2. Cross-section evaluation
First, we consider zero-field conditions. The photoionization cross section in dipole and single-particle approximation is proportional to the energy of the absorbed photon, ℰ, and the matrix element of the dipole moment (in the usual notation):
where |nlm⟩ is the bound initial state, |kl′m′⟩ the final state (the eigenstate for the ejected electron), r the electron position in the atom, and eq the unit polarization vector of the incident photon. Transitions with q = −1, 0, +1 are the so-called σ−, π, and σ+ transitions and are referred to as left-handed circular, linear, and right-handed circular polarizations, respectively. The main quantum number, n, and its analogous for the continuum, k, are related to the energies of the bound and continuum atomic states and to the energy of the absorbed photon, in Rydberg (Ry) units:
According to the Wigner-Eckart theorem, the matrix element in Eq. (1) can be separated into geometrical and physical parts,
where 𝒲 denotes the square of the Wigner (3-j) coefficient,
and the rest of terms in Eq. (3) are independent of projections of angular momenta (m, q, and m′). The cross section for transitions nl → kl′ follows the usual rule of summing over the final states and averaging over the initial states,
Combining Eqs. (3) and (5), one obtains
where a sum property of the Wigner coefficients was used, . On the other hand, σnl, kl′ can be written in terms of the mean cross section for transitions n → k ,
with Pnk and Qnl, kl′ polynomials on k whose explicit expressions can be found in the literature (e.g., Menzel & Pekeris 1935; Hatanaka 1946). With Eqs. (6) and (7), summing over all final sublevels, we obtain the cross section for a given photon polarization (q = m′−m) and for allowed dipolar transitions (l′=l ± 1) from a specific bound sublevel to the continuum,
Equation (8) is exact in the absence of a magnetic field. The proposal of Lamb & Sutherland (1972) to evaluate the cross section when a magnetic field is present consists in assuming that the initial and final wave functions remain essentially unchanged in a small region of overlap such that the matrix element in Eq. (1) remains unaltered. The cross section can thus be approximated as
where ℰ* and ℰ0 are the ionization thresholds with and without a magnetic field, and σ0 denotes the field-free cross section. This approach is asymptotically correct at B → 0 and, when combined with precise ionization threshold energies (Sect. 4), provides a surprisingly good approximation at mean field strengths as high as 100–200 MG.
Figure 1 compares results from Eqs. (8) and (9) with precise evaluations from Zhao & Stancil (2007) for the lowest field strengths published in that work. The cross sections based on full quantum mechanic calculations show rich resonance structures that arise due to quasi-bound Coulombic states embedded in the Landau continua. The evaluations based on Eqs. (8) and (9) are able to reproduce the mean values of the full theory. The comparison is even better if one takes into account the fact that much of the resonance structure can play a minor role in the formation of spectra. Since the magnetic field on the stellar surface varies significantly, by up to a factor of two for a dipole field, most narrow resonances will be smeared out or smoothed as they are averaged over even a fraction of the visible disk. These results support the use of Eqs. (8) and (9), supplemented with accurate energy data, as an empirical approach that provides reasonable cross sections at moderately high field strengths beyond the linear Zeeman regime.
![]() |
Fig. 1. Partial photoionization cross section in megabarns as a function of the light wavelength for different nlm sublevels ( |
3. Atomic energies and occupation numbers
The Hamiltonian of a non-relativistic electron in a fixed Coulomb potential and a uniform magnetic field, B (directed along the z-axis) can be written using cylindrical coordinates (ϕ, ρ, z) and atomic units (except energy in Ry) as
with , β = B/B0 a dimensionless field strength (B0 = 4.70103 × 109 G), and
and
the components of the orbital and spin angular momenta in the z-axis. Eigenstates of H can be identified via the magnetic and spin quantum numbers (m and
) complemented with the principal (n) and orbital (l) ones in the weak-field regime, and by the Landau number, N, and the longitudinal quantum number, ν (associated with excitations of the wave function along the field lines) in the strong-field regime. The correspondence between (n, l) and (N, ν) is based on the non-crossing rule (Simola & Virtamo 1978), and their explicit relations are given in Vera-Rueda & Rohrmann (2020).
Eigenvalues of H have contributions (ℰ′) from the reduced Hamiltonian, H′ (which contains the diamagnetic term, quadratic in β) and paramagnetic contributions (linear in β):
According to Eq. (11), states with m > 0 lie 4βm above those with −m,
Hereafter, we use a compact notation, ℰm ≡ ℰnlmms, where the magnetic quantum number is highlighted because of its relevance in the discussion. High-precision values for ℰm ≤ 0 were provided by Schimeczek & Wunner (2014) and fitted in Vera-Rueda & Rohrmann (2020) for a large number of bound states as a function of β. On the other hand, continuum energies that define the ionization thresholds are specified by Landau states of a free electron with linear momentum, pz, along the field
with nρ = 0, 1, …, the radial quantum number of a Landau state.
The ionization equilibrium and atomic populations were calculated as described in Vera-Rueda & Rohrmann (2020). Since the field strengths considered here are not too high, we used atomic internal partition in the form given by Pavlov & Meszaros (1993), with finite-velocity effects incorporated using effective atomic masses. Most of the MWDs have fields between 10 and 100 MG, a range in which the onset of profound changes in the eigenenergy spectrum and the distribution of atomic populations takes place. As an illustrative case, Fig. 2 shows occupation numbers and energies corresponding to sublevels (n, l)=(11, 6) for a temperature of 20 000 K. Energies of states with m > 0 are increased by multiples of cyclotron energy ℏωe = 4β with respect to levels ℰ−|m| (Eq. (12)), and their populations decrease accordingly due to the usual Boltzmann factor. The abundance of atoms in spin-up states is 48% at B = 14 MG and falls to 30% at 117 MG. Around 32 MG there is a peak of 3.5% of atoms in metastable bound states (m > 0) above the first Landau level.
![]() |
Fig. 2. Occupation numbers (ξ; top) and energies of bound energies (bottom) for sublevels (n, l)=(11, 6) for various B (gas temperature 20 000 K and density 10−8 g/cm3), shown as circles. The ionization edges (red squares) for σ+ transitions and the size of the cyclotron energy are also shown in the lower panel. |
4. Absolute photoionization opacity
The monochromatic extinction coefficient due to photoionization for the three polarizations, q = 0, ±1, is given by
where ξnlmms is the number density of particles in a bound state and the sum extends over all bound states, which we label with the quantum numbers of an atom in the zero-field limit.
First we consider transitions from spin-down levels. The minimum photoionization energy, ℰ*, of a bound state (given by Eqs. (12) and (13) with pz = 0, nρ = 0, ) depends on the photon polarization (q = m′−m), which characterizes the corresponding continuum edge (Jordan 1992; Merani et al. 1995),
In general, ℰ* equals the binding energy, −ℰ−|m|, except in two notable cases. First, for right-handed circular polarization (q = +1) and m ≥ 0, the continuum edge has an additional increment of 4β, which is the cyclotron energy, ℏωe, expressed in Ry. Second, for left-handed circular polarization (q = −1) and positive m, the continuum threshold is reduced by 4β, which can place the continuum edge near or even below the bound state for a strong enough field; in this case, photoionization takes place from zero photon energy. This scheme is the same for spin-up states because the spin-field coupling affects both discrete and continuum states equally and the spin does not change in the transition; the only difference is that the initial spin-down levels are statistically more populated and contribute more to the opacity.
Figure 3 shows the contributions to the opacity arising from sublevels at (n, l)=(11, 6) for different field strengths and polarizations. At relatively low field strengths, B < 200 MG (β < 0.05), relevant binding energies, −ℰ−|m|, are moderately modified around their zero-field values, Ry/n2, which are lower than the cyclotron energy for . This means that a great number of transitions from excited sublevels, especially π and σ− transitions, are shifted toward infrared and blue wavelengths as the field strength increases. However, the most significant changes occur for σ+ transitions from B ≈ 14 MG; this is when the cyclotron energy becomes appreciably greater than most −ℰ−|m| values (as shown in Fig. 2). In such a case, σ+ photoionization has two origins: transitions from states with negative m have thresholds distributed as in the previous cases, while those from m ≥ 0 move to energies higher by ℏωe = 4β. As shown in Fig. 3, there is a clear second jump in χ+ on the cyclotron fundamental, and it moves toward higher energies as the field size increases. This occurs with similar characteristics for sublevels originated from other n and l, n > (4β)−1/2; all those with m ≥ 0 reinforce the σ+ continuum with an edge along the cyclotron line.
![]() |
Fig. 3. Contributions to the extinction χq from sublevels (n, l)=(11, 6) as a function of the photon energy, for the conditions shown in Fig. 2. |
Figure 4 shows the total photoionization opacities at q = 0, ±1 produced by a hydrogen gas as a function of the photon wavelength (λ) and for various field strengths. They are compared with the zero-field opacity, which shows the usual discontinuities associated with initial levels n = 1 − 17. At low strengths (4.7 MG), the magnetic field mainly affects highly excited states, so the largest opacity changes occur at infrared and longer wavelengths. When B exceeds 14 MG, a bulge in χ+ forms at λ shorter than the cyclotron resonance. As the field strengthens, the bulge transforms into a well-defined jump followed by a deep trough at longer wavelengths. χ− opacities are very similar to χ0 opacities, with minor differences arising from ionization rules and Wigner coefficients; both are larger than χ+ at wavelengths longer than the cyclotron line. The opposite occurs at shorter wavelengths, where the bulge of σ+ transitions from m ≥ 0 sublevels takes place.
![]() |
Fig. 4. Monochromatic extinction coefficient due to photoionization for the three basic polarizations, q = 0, ±1, and for a hydrogen gas with T = 20 000 K and ρ = 10−8 g cm−3. The extinctions for B = 0 (thin gray lines) and cyclotron resonance (light blue line) are also shown. |
Figure 4 clearly shows how the photoionization continuum of a magnetized hydrogen gas becomes strongly dichroic (q-dependent), and the medium reveals handedness (circular chirality) through marked differences between χ− and χ+. While the total abundance of neutral atoms (mostly in the ground state) is slightly affected by the magnetic field for the studied strengths (the Lyman continuum remains almost unchanged at B ≲ 100 MG), the occupation numbers of the excited levels modeling the bf continuum in the visible and infrared spectra are strongly modified. Consequently, self-consistent calculations are needed to bring the continuum opacity into agreement with the chemical equilibrium for magnetized gases. Such calculations reveal the formation of a strong, right-handed circularly polarized absorption next to the cyclotron line, formed by the superposition of several thousand continua coming from Hilbert subspaces m ≥ 0. Much of the resonance pattern in the cross section of each isolated transition, as obtained in full quantum theory, is expected to be smoothed due to such a superposition. This provides additional support for the use of the ansatz based on Eqs. (8) and (9) for MWD atmospheres with low and moderate field strengths.
The cyclotron resonance, which is completely right-handed circularly polarized, is displayed in Fig. 4 as a comparison. It was calculated following the prescription given by Lamb & Sutherland (1974), which includes Doppler broadening and produces a very intense and very narrow absorption, as shown in the figure. Even with collision broadening included, this free-free opacity remains extremely sharp. However, Martin & Wickramasinghe (1979) showed that the cyclotron absorption cannot produce very deep features in the spectrum for realistic field patterns on the stellar disk; it can at most produce a smooth, extended depression. Similarly, the strong bf absorption reported here could also be smoothed by field spread. While field broadening may reduce the effects of the proposed opacity in the flux spectra, the expected strong continuous circular polarization signature (with a strong asymmetry around ℏωe) may be detectable in polarization measurements.
The implications of these results on the flux and polarization spectra observed in MWDs will be investigated in a future study.
References
- Amorim, L. L., Kepler, S. O., Külebi, B., Jordan, S., & Romero, A. D. 2023, ApJ, 944, 56 [NASA ADS] [CrossRef] [Google Scholar]
- Bagnulo, S., & Landstreet, J. D. 2021, MNRAS, 507, 5902 [NASA ADS] [CrossRef] [Google Scholar]
- Berdyugin, A. V., Piirola, V., Bagnulo, S., Landstreet, J. D., & Berdyugina, S. V. 2022, A&A, 657, A105 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Delande, D., Bommier, A., & Gay, J. C. 1991, Phys. Rev. Lett., 66, 141 [CrossRef] [Google Scholar]
- Ferrario, L., Wickramasinghe, D., & Kawka, A. 2020, Adv. Space Res., 66, 1025 [NASA ADS] [CrossRef] [Google Scholar]
- Hardy, F., Dufour, P., & Jordan, S. 2023, MNRAS, 520, 6111 [CrossRef] [Google Scholar]
- Hatanaka, T. 1946, Jpn. J. Astron. Geophys., 21, 1 [Google Scholar]
- Jordan, S. 1992, A&A, 265, 570 [NASA ADS] [Google Scholar]
- Külebi, B., Jordan, S., Euchner, F., Gänsicke, B. T., & Hirsch, H. 2009, A&A, 506, 1341 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Lamb, F. K., & Sutherland, P. G. 1972, Line Formation in the Presence of Magnetic Fields (Boulder, CO: National Center for Atmospheric Research), 183 [Google Scholar]
- Lamb, F. K., & Sutherland, P. G. 1974, in Physics of Dense Matter, ed. C. J. Hansen, 53, 265 [CrossRef] [Google Scholar]
- Martin, B., & Wickramasinghe, D. T. 1979, MNRAS, 189, 69 [NASA ADS] [CrossRef] [Google Scholar]
- Menzel, D. H., & Pekeris, C. L. 1935, MNRAS, 96, 77 [NASA ADS] [Google Scholar]
- Merani, N., Main, J., & Wunner, G. 1995, A&A, 298, 193 [NASA ADS] [Google Scholar]
- Pavlov, G. G., & Meszaros, P. 1993, ApJ, 416, 752 [NASA ADS] [CrossRef] [Google Scholar]
- Schimeczek, C., & Wunner, G. 2014, ApJS, 212, 26 [CrossRef] [Google Scholar]
- Simola, J., & Virtamo, J. 1978, J. Phys. B At. Mol. Phys., 11, 3309 [NASA ADS] [CrossRef] [Google Scholar]
- Vera-Rueda, M., & Rohrmann, R. D. 2020, A&A, 635, A180 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Vera-Rueda, M., & Rohrmann, R. D. 2024, A&A, 687, A141 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Wang, Q., & Greene, C. H. 1991, Phys. Rev. A, 44, 7448 [NASA ADS] [CrossRef] [Google Scholar]
- Zhao, L. B. 2021, ApJS, 254, 21 [NASA ADS] [CrossRef] [Google Scholar]
- Zhao, L. B., & Stancil, P. C. 2007, ApJ, 667, 1119 [NASA ADS] [CrossRef] [Google Scholar]
All Figures
![]() |
Fig. 1. Partial photoionization cross section in megabarns as a function of the light wavelength for different nlm sublevels ( |
In the text |
![]() |
Fig. 2. Occupation numbers (ξ; top) and energies of bound energies (bottom) for sublevels (n, l)=(11, 6) for various B (gas temperature 20 000 K and density 10−8 g/cm3), shown as circles. The ionization edges (red squares) for σ+ transitions and the size of the cyclotron energy are also shown in the lower panel. |
In the text |
![]() |
Fig. 3. Contributions to the extinction χq from sublevels (n, l)=(11, 6) as a function of the photon energy, for the conditions shown in Fig. 2. |
In the text |
![]() |
Fig. 4. Monochromatic extinction coefficient due to photoionization for the three basic polarizations, q = 0, ±1, and for a hydrogen gas with T = 20 000 K and ρ = 10−8 g cm−3. The extinctions for B = 0 (thin gray lines) and cyclotron resonance (light blue line) are also shown. |
In the text |
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.