X-ray bursting neutron star atmosphere models using an exact relativistic kinetic equation for Compton scattering

Theoretical spectra of X-ray bursting neutron star (NS) model atmospheres are widely used to determine the basic NS parameters such as their masses and radii. We construct accurate NS atmosphere models using for the first time an exact treatment of Compton scattering via the integral relativistic kinetic equation. We also compare the results with the previous calculations based on the Kompaneets operator. We solve the radiation transfer equation together with the hydrostatic equilibrium equation accounting exactly for the radiation pressure by electron scattering. We thus construct a new set of plane-parallel atmosphere models in LTE for hot NSs. The models were computed for six chemical compositions (pure H, pure He, solar H/He mix with various heavy elements abundances Z = 1, 0.3, 0.1, and 0.01 Z_sun, and three log g = 14.0, 14.3, and 14.6. For each chemical composition and log g, we compute more than 26 model atmospheres with various luminosities relative to the Eddington luminosity L_Edd computed for the Thomson cross-section. The maximum relative luminosities L/L_Edd reach values of up to 1.1 for high gravity models. The emergent spectra of all models are redshifted and fitted by diluted blackbody spectra in the 3–20 keV energy range appropriate for the RXTE/PCA. We also compute the color correction factors f_c. The radiative acceleration g_rad in our luminous, hot-atmosphere models is significantly smaller than in corresponding models based on the Kompaneets operator, because of the Klein-Nishina reduction of the electron scattering cross-section, and therefore formally “super-Eddington” model atmospheres do exist. The differences between the new and old model atmospheres are small for L / L_Edd < 0.8. For the same g_rad / g, the new f_c are slightly larger (by approximately 1%) than the old values.

💡 Research Summary

The paper presents a new generation of neutron‑star (NS) atmosphere models that treat Compton scattering with an exact relativistic kinetic equation rather than the approximate Kompaneets operator traditionally used in such calculations. By solving the radiation‑transfer equation together with hydrostatic equilibrium while accounting for the full Klein‑Nishina reduction of the electron‑scattering cross‑section, the authors obtain a self‑consistent description of the radiative acceleration (g_rad) and emergent spectra for hot, luminous NS atmospheres.

A comprehensive grid of plane‑parallel, LTE models is constructed. Six chemical compositions are considered: pure hydrogen, pure helium, and a solar H/He mixture with heavy‑element abundances Z = 1, 0.3, 0.1, 0.01 Z⊙. For each composition three surface gravities (log g = 14.0, 14.3, 14.6) are used, and for each (composition, log g) pair more than 26 luminosity points are calculated, spanning relative luminosities L/L_Edd from ≈0.1 up to ≈1.1 (the Eddington limit is defined using the Thomson cross‑section). The maximum L/L_Edd reaches 1.1 for the highest‑gravity models, indicating that “super‑Eddington” atmospheres can exist when the Klein‑Nishina reduction is taken into account.

The emergent spectra are redshifted and fitted in the 3–20 keV band (the range of the RXTE/PCA instrument) with diluted blackbody functions of the form w Bν(T_c), where the dilution factor w = f_c⁻⁴ and f_c is the color‑correction factor. The authors find that, because the exact treatment lowers the effective scattering opacity, the radiative acceleration g_rad is significantly smaller than in Kompaneets‑based models. Consequently, models with L/L_Edd > 1 are dynamically stable, a situation that cannot be reproduced with the Kompaneets approximation. For luminosities below about 0.8 L_Edd, the differences between the two approaches are modest; however, at higher luminosities the new models yield color‑correction factors that are roughly 1 % larger than the older values for the same g_rad/g ratio.

The study highlights several key physical insights. First, the Klein‑Nishina effect is crucial for accurately modeling the radiation pressure in hot NS atmospheres, especially near the Eddington limit. Second, the dependence of f_c on chemical composition, surface gravity, and luminosity is quantified, showing that lower heavy‑element abundances and higher gravities tend to produce harder spectra (larger f_c). Third, the precise calculation of g_rad enables the identification of genuinely super‑Eddington atmospheres, which has implications for interpreting photospheric radius‑expansion (PRE) bursts and for constraining NS mass–radius relations.

By providing a more accurate set of color‑correction factors and radiative‑acceleration profiles, the paper supplies essential input for the analysis of X‑ray burst observations from instruments such as RXTE, NICER, and the upcoming eXTP mission. The authors suggest that future extensions could incorporate non‑LTE effects, magnetic fields, rapid rotation, and multi‑dimensional geometry, thereby further refining the connection between observed burst spectra and the underlying neutron‑star equation of state.