High Energy Astrophysics 3 JAN, 2015

Magnetized neutron star atmospheres: beyond cold plasma approximation

Magnetized neutron star atmospheres: beyond cold plasma approximation

All the neutron star (NS) atmosphere models published so far have been calculated in the “cold plasma approximation”, which neglects the relativistic effects in the radiative processes, such as cyclotron emission/absorption at harmonics of cyclotron frequency. Here we present new NS atmosphere models which include such effects. We calculate a set of models for effective temperatures T_eff =1-3 MK and magnetic fields B \sim 10^{10}-10^{11} G, typical for the so-called central compact objects (CCOs) in supernova remnants, for which the electron cyclotron energy E_{c,e} and its first harmonics are in the observable soft X-ray range. Although the relativistic parameters, such as kT_eff /(m_e c^2) and E_{c,e} /(m_e c^2), are very small for CCOs, the relativistic effects substantially change the emergent spectra at the cyclotron resonances, E \approx sE_{c,e} (s=1, 2,…). Although the cyclotron absorption features can form in a cold plasma due to the quantum oscillations of the free-free opacity, the shape and depth of these features change substantially if the relativistic effects are included. In particular, the features acquire deep Doppler cores, in which the angular distribution of the emergent intensity is quite different from that in the cold plasma approximation. The relative contributions of the Doppler cores to the equivalent widths of the features grow with increasing the quantization parameter b_eff = E_{c,e}/kT_eff and harmonic number s. The total equivalent widths of the features can reach \sim 150-250 eV; they increase with growing b_eff and are smaller for higher harmonics.

💡 Research Summary

The paper addresses a long‑standing limitation in neutron‑star (NS) atmosphere modeling: the widespread use of the “cold plasma approximation,” which neglects relativistic effects in radiative processes such as cyclotron emission and absorption at harmonic frequencies. While this approximation is justified when both the electron temperature (kT_eff) and the electron cyclotron energy (E_c,e) are much smaller than the electron rest‑mass energy (m_ec²), it fails for a class of objects known as central compact objects (CCOs) in supernova remnants. CCOs typically have surface magnetic fields B ≈ 10¹⁰–10¹¹ G and effective temperatures T_eff ≈ 1–3 MK, placing the fundamental electron cyclotron energy and its first few harmonics squarely in the soft X‑ray band (∼0.1–2 keV). Consequently, even though the relativistic parameters kT_eff/(m_ec²) and E_c,e/(m_ec²) are of order 10⁻³, relativistic Doppler broadening and harmonic emission become spectroscopically significant.

To overcome this, the authors construct a new set of NS atmosphere models that incorporate full relativistic treatment of cyclotron processes. The key ingredients of the model are:

  1. Quantum‑oscillatory free‑free opacity – retained from the cold‑plasma framework, representing the periodic modulation of the free‑free absorption coefficient due to Landau quantization of electron motion.
  2. Relativistic cyclotron transition rates – derived from the full quantum‑electrodynamic (QED) cross‑sections, including both absorption and spontaneous emission, and explicitly accounting for the Doppler shift of electrons moving along the magnetic field.
  3. Quantization parameter b_eff = E_c,e / kT_eff – introduced as a dimensionless measure of how many cyclotron quanta fit within the thermal energy spread. Larger b_eff values correspond to stronger relativistic effects.
  4. Angular dependence of emergent intensity – the model solves the radiative transfer equation for different propagation angles relative to the magnetic field, thereby capturing the anisotropic beaming that arises from relativistic motion.

The authors compute a grid of atmosphere models for B = 10¹⁰–10¹¹ G, T_eff = 1–3 MK, and various viewing angles. The results reveal several crucial departures from the cold‑plasma predictions:

  • Deep Doppler cores appear at each cyclotron resonance (E ≈ s E_c,e, with harmonic number s = 1, 2, …). These cores are much narrower (tens of eV) and deeper (optical depth up to several) than the broad, shallow features produced solely by quantum free‑free oscillations.
  • Equivalent widths (EWs) of the combined features (core + surrounding quantum oscillation) can reach 150–250 eV for the fundamental (s = 1) and remain sizable (80–150 eV) for the second and third harmonics, especially when b_eff > 8. The EW grows with increasing b_eff and with harmonic number, but the relative contribution of the Doppler core becomes more pronounced for higher s.
  • Angular redistribution is dramatic: within the Doppler core the emergent intensity is strongly beamed perpendicular to the magnetic field, whereas outside the core the intensity follows the broader pattern dictated by the quantum free‑free opacity. This anisotropy implies that the observed depth of a cyclotron line can vary dramatically with the observer’s line‑of‑sight, offering a potential diagnostic of magnetic geometry.
  • Spectral shape modification – the inclusion of relativistic effects not only deepens the lines but also shifts their centroids slightly due to the average Doppler shift of the thermal electron distribution. The shift is modest (a few percent of E_c,e) but measurable with high‑resolution spectrometers.

These findings have immediate observational relevance. Current and upcoming X‑ray missions (e.g., XRISM Resolve, Athena X‑IFU) will achieve energy resolutions of ∼2–5 eV in the soft X‑ray band, sufficient to resolve the narrow Doppler cores predicted here. Detecting such cores would provide a direct measurement of b_eff, and thus of the ratio E_c,e/kT_eff, allowing independent constraints on both surface magnetic field strength and temperature. Moreover, the angular beaming pattern encoded in the line profiles could be exploited to infer the magnetic inclination and viewing geometry of CCOs, complementing pulse‑profile modeling.

The paper also discusses limitations and future directions. The present models assume a plane‑parallel, uniform‑magnetic‑field atmosphere and ignore possible magnetic field gradients, surface temperature inhomogeneities, and higher‑order QED effects (e.g., vacuum polarization) that become important at B ≳ 10¹³ G. Extending the framework to multi‑dimensional radiative transfer and incorporating realistic magnetic topology will be essential for applying the method to a broader class of neutron stars, including magnetars and rotation‑powered pulsars.

In summary, by moving beyond the cold‑plasma approximation and incorporating relativistic cyclotron physics, the authors demonstrate that even modest magnetic fields and temperatures can produce pronounced, Doppler‑shaped cyclotron absorption features in NS atmospheres. This work opens a new avenue for high‑precision X‑ray spectroscopy of CCOs and provides a robust theoretical foundation for interpreting forthcoming high‑resolution observations.