Logo
ScienceStack
Table of Contents
Fetching ...
Sign In

Modeling X-ray Bursting Neutron Star Atmospheres

Lourenzo Colleyn, Zach Medin, Alan Calder

TL;DR

This work verifies a LANL X-ray bursting neutron star atmosphere model (ZCODE) by showing convergence of radiative-transfer calculations to steady-state atmospheric snapshots under fixed composition and surface gravity. The approach uses an Implicit Monte Carlo method with a base-temperature proxy for luminosity, 100 atmospheric cells, and 300 frequency groups to capture frequency-dependent opacities. Results reveal a diluted blackbody-like emergent spectrum shaped by frequency-dependent opacity and Compton scattering, with a color-correction factor fc(l) that matches prior theory and models, including a local minimum near mid-luminosity and an exponential rise at high luminosity. The study demonstrates robust convergence across modeling choices and provides a framework for extending to different compositions and gravities, advancing the interpretation of X-ray burst observations for neutron star properties.

Abstract

We present a verification of a computational model, developed at the Los Alamos National Laboratory (LANL) for simulating radiation transfer in X-ray bursting neutron star atmospheres. We tested a baseline case and demonstrated strong agreement in the behavior of the outgoing spectrum's color-correction factor with earlier work and theoretical expectations. By analyzing the relationship between the simulation time and outgoing flux, we also demonstrated how the model calculates through a sequence of time-independent atmospheric snapshots, each iteratively refined, and uses them to progressively converge toward the correct atmospheric state (as would be observed during a burst). We examined the behavior of the outgoing flux across different optical depths and explored the physical explanations for deviations from a pure blackbody spectrum, attributed to frequency-dependent opacity sources. Additionally, we assessed the impact of Compton scattering, highlighting its role in redistributing photon energies.

Modeling X-ray Bursting Neutron Star Atmospheres

TL;DR

This work verifies a LANL X-ray bursting neutron star atmosphere model (ZCODE) by showing convergence of radiative-transfer calculations to steady-state atmospheric snapshots under fixed composition and surface gravity. The approach uses an Implicit Monte Carlo method with a base-temperature proxy for luminosity, 100 atmospheric cells, and 300 frequency groups to capture frequency-dependent opacities. Results reveal a diluted blackbody-like emergent spectrum shaped by frequency-dependent opacity and Compton scattering, with a color-correction factor fc(l) that matches prior theory and models, including a local minimum near mid-luminosity and an exponential rise at high luminosity. The study demonstrates robust convergence across modeling choices and provides a framework for extending to different compositions and gravities, advancing the interpretation of X-ray burst observations for neutron star properties.

Abstract

We present a verification of a computational model, developed at the Los Alamos National Laboratory (LANL) for simulating radiation transfer in X-ray bursting neutron star atmospheres. We tested a baseline case and demonstrated strong agreement in the behavior of the outgoing spectrum's color-correction factor with earlier work and theoretical expectations. By analyzing the relationship between the simulation time and outgoing flux, we also demonstrated how the model calculates through a sequence of time-independent atmospheric snapshots, each iteratively refined, and uses them to progressively converge toward the correct atmospheric state (as would be observed during a burst). We examined the behavior of the outgoing flux across different optical depths and explored the physical explanations for deviations from a pure blackbody spectrum, attributed to frequency-dependent opacity sources. Additionally, we assessed the impact of Compton scattering, highlighting its role in redistributing photon energies.
Paper Structure (14 sections, 6 equations, 5 figures)

This paper contains 14 sections, 6 equations, 5 figures.

Table of Contents

  1. Introduction
  2. Equations
  3. Simulation Methods
  4. Results
  5. Total Outgoing Flux
  6. Outgoing Spectrum
  7. Color-correction Factors
  8. A note on performance
  9. Discussion
  10. Behavior of the Total Outgoing Flux
  11. The Dilute Blackbody Distribution
  12. Outgoing Spectra at Different Luminosities
  13. Color-Correction Factors
  14. Conclusion

Figures (5)

  • Figure 1: Schematic depicting the atmosphere's spherical cell structure with material quantities defined within each discrete cell and at the base. Note that in the terminology of this section, $r_{\rm base} = R$, and $g_{\rm base} = g$ is assumed constant.
  • Figure 2: Total outgoing flux vs. cycle number at $l = 0.5$ with the fiducial value of $\tau_{\rm{min}}$ (left panel) and total outgoing flux vs. cycle number at $l = 0.5$ with varying values of $\tau_{\rm{min}}$.
  • Figure 3: Total outgoing flux vs. cycle number at $l = 0.5$ showing simulations with different numbers of photon packets. The panel on the right is detailed view of the results at large cycle numbers.
  • Figure 4: Outgoing radiation spectrum at $l = 0.5$ (blue curve) with the blackbody approximation at the atmosphere effective temperature (black curve) (left panel). Outgoing radiation spectra at different luminosities, ranging from $l = 0.003$ to $l = 1.02$ (right panel). Note that in the panel labels the luminosity ratio is referred to as $l_{\rm proj}$; the added subscript there refers to the luminosity "projected" on to the base of the atmosphere, as described in Medin_2016.
  • Figure 5: Plot showing outgoing luminosity (in units of the Eddington luminosity) vs. color-correction values, for both ZCODE and SPW12 Suleimanov_2012. Comparable line trends and color-correction values indicate the validity of ZCODE's methods.