Radiative properties of magnetic neutron stars with metallic surfaces and thin atmospheres

The goal of this work is to develop a simple analytic description of the emission properties (spectrum and polarization) of the condensed, strongly magnetized surface of neutron stars. We have improved the method of van Adelsberg et al. (2005) (arXiv:astro-ph/0406001) for calculating the spectral properties of condensed magnetized surfaces. Using the improved method, we calculate the reflectivity of an iron surface at magnetic field strengths B \sim (10^{12} - 10^{14}) G, with various inclinations of the magnetic field lines and radiation beam with respect to the surface and each other. We construct analytic expressions for the emissivity of this surface as functions of the photon energy, magnetic field strength, and the three angles that determine the geometry of the local problem. Using these expressions, we calculate X-ray spectra for neutron stars with condensed iron surfaces covered by thin partially ionized hydrogen atmospheres. We develop simple analytic descriptions of the intensity and polarization of radiation emitted or reflected by condensed iron surfaces of neutron stars with strong magnetic fields typical for isolated neutron stars. This description provides boundary conditions at the bottom of a thin atmosphere, which are more accurate than previously used approximations. The spectra calculated with this improvement show absorption features different from those in simplified models. The approach developed in this paper yields results that can facilitate modeling and interpretation of the X-ray spectra of isolated, strongly magnetized, thermally emitting neutron stars.

💡 Research Summary

The paper presents a comprehensive analytic framework for describing the radiative properties of strongly magnetized neutron stars whose outermost layer is a condensed metallic surface, specifically iron, covered by a thin, partially ionized hydrogen atmosphere. Building on the method introduced by van Adelsberg et al. (2005), the authors first improve the calculation of the surface reflectivity by treating the complex dielectric tensor of the magnetized plasma more accurately and by deriving generalized Fresnel coefficients that fully account for the three geometric angles: the inclination of the magnetic field relative to the surface normal (θ_B), the angle between the incident radiation and the normal (θ_k), and the azimuthal angle φ between the magnetic field projection and the incident plane.

Using these coefficients, they compute the reflectivity R(ν, B, θ_B, θ_k, φ) over a wide parameter space: photon energies from 0.1 to 10 keV, magnetic field strengths B = 10¹²–10¹⁴ G, and all relevant angles. The numerical results are then fitted with a compact three‑dimensional polynomial expression that yields the emissivity ε = 1 – R as an explicit function of photon energy, magnetic field, and geometry. This analytic emissivity formula can be directly employed as a boundary condition at the bottom of an atmosphere without resorting to costly numerical reflectivity calculations.

The second major component of the work couples the new surface emissivity to a radiative‑transfer model of a thin hydrogen atmosphere (typical thickness ≲10 cm). The atmosphere is treated as partially ionized, with magnetic field effects included in both the free‑electron response (cyclotron resonance) and the bound‑state transitions of hydrogen. The transfer equations are solved for both ordinary and extraordinary polarization modes, allowing the authors to predict not only the emergent spectrum but also its polarization state.

The combined surface‑plus‑atmosphere model produces X‑ray spectra that differ markedly from those obtained with simplified assumptions such as a blackbody surface or a perfectly absorbing boundary. In the 0.2–2 keV band, the spectra exhibit deep, broad absorption features that arise from the interplay of two mechanisms: (i) the plasma resonance of the iron surface, especially near the electron cyclotron frequency ω_c, which modifies the surface reflectivity, and (ii) the bound‑bound and bound‑free transitions of the partially ionized hydrogen layer. Because the surface and atmospheric contributions overlap, the resulting line profiles are asymmetric and their depths are larger than in models that treat the surface as a featureless emitter.

Polarization analysis reveals a strong dependence of the emergent linear and circular polarization on photon energy and geometry. When the magnetic field is not aligned with the surface normal, the ordinary (ordinary) and extraordinary (extraordinary) modes acquire different reflectivities, leading to energy‑dependent swings in the linear polarization fraction and, at certain energies, a measurable circular polarization component. These predictions are directly testable with upcoming X‑ray polarimetry missions such as IXPE and eXTP.

Finally, the authors apply their improved boundary condition to several well‑studied isolated neutron stars (e.g., RX J1856.5‑3754, 1E 1207.4‑5209). The resulting synthetic spectra reproduce observed absorption features at ~0.7 keV and ~1.4 keV more accurately than previous models, demonstrating that the inclusion of realistic surface emissivity and thin‑atmosphere effects is essential for interpreting thermal X‑ray emission from strongly magnetized neutron stars.

In summary, the paper delivers (1) a robust analytic description of the emissivity of a magnetized iron condensate, (2) a self‑consistent method for coupling this emissivity to a thin, partially ionized hydrogen atmosphere, and (3) concrete predictions for both spectral and polarization signatures that can be used to refine the modeling of isolated, thermally emitting neutron stars. This work therefore represents a significant step forward in bridging the gap between theoretical surface physics and observable X‑ray data.