[Paper Review] Cyclotron resonant scattering feature simulations. I. Thermally averaged cyclotron scattering cross sections, mean free photon-path tables, and electron momentum sampling

Fritz-Walter Schwarm, G. Schönherr+12 more|arXiv (Cornell University)|Sep 16, 2016

Pulsars and Gravitational Waves Research18 references14 citations

TL;DR

This paper presents a computationally efficient tabular interpolation scheme for simulating cyclotron resonant scattering features (CRSFs) in X-ray pulsars. By precomputing thermally averaged mean free paths and electron momentum distributions using adaptive numerical integration, the method accelerates Monte Carlo simulations—reducing time-intensive cross-section calculations—while maintaining a relative error of ≤1/15 across a range of magnetic fields (0.01–0.12 Bcrit) and temperatures (3–15 keV). The tables are publicly available and enable high-fidelity synthetic CRSF spectrum generation.

ABSTRACT

Context. Electron cyclotron resonant scattering features (CRSFs) are observed as absorption-like lines in the spectra of X-ray pulsars. A significant fraction of the computing time for Monte Carlo simulations of these quantum mechanical features is spent on the calculation of the mean free path for each individual photon before scattering, since it involves a complex numerical integration over the scattering cross section and the (thermal) velocity distribution of the scattering electrons. Aims. We aim to numerically calculate interpolation tables which can be used in CRSF simulations to sample the mean free path of the scattering photon and the momentum of the scattering electron. The tables also contain all the information required for sampling the scattering electron’s final spin. Methods. The tables were calculated using an adaptive Simpson integration scheme. The energy and angle grids were refined until a prescribed accuracy is reached. The tables are used by our simulation code to produce artificial CRSF spectra. The electron momenta sampled during these simulations were analyzed and justified using theoretically determined boundaries. Results. We present a complete set of tables suited for mean free path calculations of Monte Carlo simulations of the cyclotron scattering process for conditions expected in typical X-ray pulsar accretion columns (0.01 ≤ B/Bcᵣᵢₜ ≤ 0.12, where Bcᵣᵢₜ = 4.413 × 10¹³ G, and 3 keV ≤ kBT ≤ 15 keV). The sampling of the tables is chosen such that the results have an estimated relative error of at most 1/15 for all points in the grid. The tables are available online (see link in footnote, page 1).

Motivation & Objective

• To reduce the computational burden of Monte Carlo simulations of cyclotron resonant scattering features (CRSFs) in X-ray pulsars.
• To develop accurate, precomputed interpolation tables for photon mean free paths and electron momentum sampling.
• To enable faster, more scalable simulations of complex X-ray pulsar geometries and physical conditions.
• To validate the physical consistency of electron momentum sampling and scattering resonance conditions.
• To provide publicly accessible, high-accuracy tables for use in future CRSF modeling and spectral fitting.

Proposed method

• Employed adaptive Simpson integration to compute thermally averaged cyclotron scattering cross sections over energy, angle, and electron momentum grids.
• Generated mean free path (MFP) tables by numerically integrating over electron velocity distributions and scattering cross sections for all possible final states and polarizations.
• Calculated electron momentum sampling boundaries using relativistic Maxwellian distributions and Doppler broadening estimates.
• Incorporated spin-dependent scattering outcomes and polarization-averaged cross sections into the tables.
• Refined energy and angle grids iteratively until a target accuracy (relative error ≤1/15) was achieved.
• Provided online access to the complete set of interpolation tables for use in Monte Carlo simulation codes.

Experimental results

Research questions

• RQ1 How can the time-consuming calculation of photon mean free paths in CRSF simulations be accelerated without sacrificing accuracy?
• RQ2 What are the physically motivated bounds on electron parallel momentum that ensure accurate sampling in resonant scattering simulations?
• RQ3 How do the resonance conditions and Doppler broadening affect the spatial distribution of scattering events in energy-momentum space?
• RQ4 To what extent do the simulated scattering events align with zero-line-width analytical solutions of the resonance condition?
• RQ5 Can the interpolation tables accurately represent the thermally averaged cross sections across a wide range of magnetic fields and plasma temperatures?

Key findings

• The interpolation tables achieve a relative error of at most 1/15 across all grid points, ensuring high numerical fidelity.
• The tables cover a physically relevant range: 0.01 ≤ B/Bcrit ≤ 0.12 and 3 keV ≤ kBT ≤ 15 keV, suitable for typical X-ray pulsar accretion columns.
• Electron momentum sampling is physically bounded by the 99%-containment region of the relativistic Maxwellian distribution, validating the simulation limits.
• Resonant scattering events are confined to regions near the zero-line-width solutions of the resonance condition, confirming the physical consistency of the method.
• The Doppler broadening estimate (ΔE/E = √(8 ln 2 kBT / mec² cos ϑ)) provides a reliable approximation for CRSF line widths in the simulation framework.
• The tables successfully enable efficient Monte Carlo simulation of synthetic CRSF spectra, allowing exploration of complex physical scenarios beyond previous capabilities.