Models of magnetized neutron star atmospheres: thin atmospheres and partially ionized hydrogen atmospheres with vacuum polarization

Observed X-ray spectra of some isolated magnetized neutron stars display absorption features, sometimes interpreted as ion cyclotron lines. Modeling the observed spectra is necessary to check this hypothesis and to evaluate neutron star parameters.We develop a computer code for modeling magnetized neutron star atmospheres in a wide range of magnetic fields (10^{12} - 10^{15} G) and effective temperatures (3 \times 10^5 - 10^7 K). Using this code, we study the possibilities to explain the soft X-ray spectra of isolated neutron stars by different atmosphere models. The atmosphere is assumed to consist either of fully ionized electron-ion plasmas or of partially ionized hydrogen. Vacuum resonance and partial mode conversion are taken into account. Any inclination of the magnetic field relative to the stellar surface is allowed. We use modern opacities of fully or partially ionized plasmas in strong magnetic fields and solve the coupled radiative transfer equations for the normal electromagnetic modes in the plasma. Spectra of outgoing radiation are calculated for various atmosphere models: fully ionized semi-infinite atmosphere, thin atmosphere, partially ionized hydrogen atmosphere, or novel “sandwich” atmosphere (thin atmosphere with a hydrogen layer above a helium layer. Possibilities of applications of these results are discussed. In particular, the outgoing spectrum using the “sandwich” model is constructed. Thin partially ionized hydrogen atmospheres with vacuum polarization are shown to be able to improve the fit to the observed spectrum of the nearby isolated neutron star RBS 1223 (RX J1308.8+2127).

💡 Research Summary

The paper presents a comprehensive numerical framework for modeling the atmospheres of isolated, strongly magnetized neutron stars (NSs) over a broad range of magnetic field strengths (10¹²–10¹⁵ G) and effective temperatures (3 × 10⁵–10⁷ K). The authors aim to explain the soft X‑ray spectra of several isolated NSs, many of which exhibit absorption features that have been interpreted as ion cyclotron lines. To test this hypothesis, they construct a versatile atmosphere code that can handle both fully ionized electron‑ion plasmas and partially ionized hydrogen, while explicitly incorporating vacuum polarization (vacuum resonance) and partial mode conversion between the two normal electromagnetic modes (the ordinary “O‑mode” and extraordinary “X‑mode”) that propagate in a magnetized plasma.

Key physical ingredients:

1. Opacities in strong B‑fields – The code uses up‑to‑date atomic data for hydrogen (including bound‑bound, bound‑free, and free‑free processes) and for fully ionized plasmas, taking into account the Landau quantization of electron motion and the modification of atomic energy levels by the magnetic field.
2. Vacuum resonance – In fields ≳10¹⁴ G the dielectric tensor of the vacuum becomes comparable to that of the plasma, producing a resonance where the O‑ and X‑mode refractive indices cross. The authors compute the energy‑ and angle‑dependent conversion probability rather than assuming either complete conversion or none, which allows a realistic treatment of the “partial mode conversion” that strongly influences line depths.
3. Geometry – The magnetic field can be inclined at any angle θ relative to the local surface normal, and the radiative transfer equations are solved in two dimensions (depth and angle) for each mode. This flexibility enables the generation of phase‑resolved spectra for rotating NSs.

Atmosphere configurations explored:

• Semi‑infinite fully ionized atmosphere (the traditional benchmark).
• Thin atmosphere where the column depth is smaller than the photon mean free path, leading to a spectrum that deviates markedly from a blackbody, especially at higher energies.
• Partially ionized hydrogen atmosphere with vacuum polarization, which produces narrower, deeper absorption features due to bound‑bound transitions of hydrogen in a strong field.
• “Sandwich” atmosphere – a novel construct consisting of a thin hydrogen layer over a thin helium layer. Because hydrogen and helium have different transition energies in a magnetic field, the combined spectrum can display multiple absorption lines that mimic the observed complex features.

Numerical method: The coupled radiative‑transfer equations for the two normal modes are solved iteratively using an accelerated Lambda‑iteration scheme. Boundary conditions enforce radiative equilibrium deep in the atmosphere and free‑streaming at the surface. The code outputs emergent spectra for any chosen set of (B, T_eff, θ) and for any of the atmosphere structures above.

Results and astrophysical implications:

• The inclusion of vacuum resonance and partial mode conversion significantly modifies the shape and depth of cyclotron‑like features, often reducing their equivalent width compared with models that neglect these effects.
• Thin partially ionized hydrogen atmospheres with vacuum polarization provide a markedly better fit to the X‑ray spectrum of the nearby isolated NS RBS 1223 (RX J1308.8+2127) than the standard fully ionized semi‑infinite models.
• The sandwich atmosphere can reproduce spectra that show two or more absorption lines at distinct energies, offering a natural explanation for sources where multiple lines are detected.
• By varying the magnetic inclination θ, the authors demonstrate how the observed pulse‑phase dependence of line strength can arise from geometric effects alone, without invoking surface temperature anisotropies.

In summary, this work delivers the first unified atmosphere modeling tool that simultaneously treats strong‑field atomic physics, vacuum polarization, and a variety of realistic atmospheric thicknesses and compositions. The resulting spectral library is directly applicable to current and upcoming X‑ray missions (e.g., NICER, Athena) and provides a robust framework for extracting neutron‑star parameters—magnetic field strength, surface temperature distribution, and composition—from high‑resolution X‑ray observations.