[Paper Review] Cyclotron resonant scattering feature simulations

Fritz-Walter Schwarm, Ralf Ballhausen+13 more|arXiv (Cornell University)|Jan 26, 2017

Astrophysical Phenomena and Observations46 references13 citations

TL;DR

This paper presents a new Monte Carlo simulation code for generating synthetic cyclotron resonant scattering features (CRSFs) in X-ray pulsars, enabling accurate modeling of complex, cylindrically symmetric accretion geometries. The method uses precomputed mean free path tables and Green's functions to efficiently simulate line profiles, which are then fitted to NuSTAR data using an XSPEC-compatible model (cyclofs), yielding a magnetic field strength of 3.6×10¹² G for Cep X-4—significantly different from Gaussian line fits due to the physically realistic CRSF shape.

ABSTRACT

Context. Cyclotron resonant scattering features (CRSFs) are formed by scattering of X-ray photons off quantized plasma electrons in the strong magnetic field (of the order 1012 G) close to the surface of an accreting X-ray pulsar. Due to the complex scattering cross-sections, the line profiles of CRSFs cannot be described by an analytic expression. Numerical methods, such as Monte Carlo (MC) simulations of the scattering processes, are required in order to predict precise line shapes for a given physical setup, which can be compared to observations to gain information about the underlying physics in these systems. Aims. A versatile simulation code is needed for the generation of synthetic cyclotron lines. Sophisticated geometries should be investigatable by making their simulation possible for the first time. Methods. The simulation utilizes the mean free path tables described in the first paper of this series for the fast interpolation of propagation lengths. The code is parallelized to make the very time-consuming simulations possible on convenient time scales. Furthermore, it can generate responses to monoenergetic photon injections, producing Green’s functions, which can be used later to generate spectra for arbitrary continua. Results. We develop a new simulation code to generate synthetic cyclotron lines for complex scenarios, allowing for unprecedented physical interpretation of the observed data. An associated XSPEC model implementation is used to fit synthetic line profiles to NuSTAR data of Cep X-4. The code has been developed with the main goal of overcoming previous geometrical constraints in MC simulations of CRSFs. By applying this code also to more simple, classic geometries used in previous works, we furthermore address issues of code verification and cross-comparison of various models. The XSPEC model and the Green’s function tables are available online (see link in footnote, page 1).

Motivation & Objective

• To overcome previous limitations in Monte Carlo simulations of cyclotron resonant scattering features (CRSFs), which were restricted to simple, predefined geometries.
• To develop a flexible, parallelized simulation code capable of modeling synthetic CRSFs in arbitrarily complex, cylindrically symmetric accretion column configurations.
• To generate Green’s functions for mono-energetic photon injections, enabling efficient spectral fitting with arbitrary continuum models.
• To provide an XSPEC-compatible model (cyclofs) that allows direct comparison of synthetic line profiles with observational data.
• To validate the code against classic geometries and previous works, ensuring consistency and reliability for future physical interpretation.

Proposed method

• The simulation uses mean free path (MFP) tables from the companion paper (Schwarm et al. 2017, Paper I) for fast interpolation of thermally averaged scattering cross sections.
• A parallelized Monte Carlo algorithm simulates photon propagation, scattering, and energy deposition in a 3D cylindrical geometry with user-defined magnetic field, temperature, and density gradients.
• The code generates Green’s functions by injecting mono-energetic photons and recording the emergent spectral response, decoupling the time-consuming simulation from continuum choice.
• The XSPEC model 'cyclofs' is implemented using precomputed Green’s function tables, allowing direct fitting of synthetic CRSF profiles to observational spectra.
• The simulation accounts for relativistic effects such as gravitational redshift and includes viewing angle dependence, which significantly affects line shape and width.
• The code supports averaging over multiple viewing angles via an approximate method to account for pulse-phase-averaged data.

Experimental results

Research questions

• RQ1 How do complex, non-uniform accretion geometries affect the shape and depth of cyclotron resonant scattering features in X-ray pulsars?
• RQ2 To what extent do viewing angle and temperature gradients influence the observed line width and centroid energy of CRSFs?
• RQ3 Can a physically grounded, non-Gaussian line profile model improve the accuracy of magnetic field inference compared to empirical Gaussian absorption models?
• RQ4 How does the inclusion of realistic scattering cross sections and propagation effects alter the derived physical parameters from observational fits?
• RQ5 What is the impact of continuum model choice on the derived CRSF parameters when using a physically consistent line shape model?

Key findings

• The new simulation code successfully generates synthetic cyclotron line profiles for complex, cylindrically symmetric geometries, overcoming prior geometric constraints in Monte Carlo simulations.
• The code produces line profiles with strong viewing-angle dependence, where the width and depth vary significantly with the angle to the magnetic field axis.
• When fitted to NuSTAR data of Cep X-4, the physical XSPEC model 'cyclofs' yields a magnetic field strength of 3.6×10¹² G, which differs significantly from values obtained with a Gaussian absorption line model.
• The best-fit physical model results in a reduced χ² of 1.17 for 875 degrees of freedom, outperforming the Gaussian model (χ² = 1.57) and showing clear residuals when using a fixed centroid energy.
• The fundamental cyclotron line shape is inherently complex and asymmetric, leading to systematic differences in magnetic field inference when non-Gaussian profiles are used, even for smooth, symmetric-appearing lines.
• The XSPEC model 'cyclofs' and associated Green’s function tables are publicly available, enabling widespread use and integration into standard X-ray spectral fitting workflows.