Accurate X-ray Timing in the Presence of Systematic Biases With Simulation-Based Inference

Abstract

Because many of our X-ray telescopes are optimized towards observing faint sources, observations of bright sources like X-ray binaries in outburst are often affected by instrumental biases. These effects include dead time and photon pile-up, which can dramatically change the statistical inference of physical parameters from these observations. While dead time is difficult to take into account in a statistically consistent manner, simulating dead time-affected data is often straightforward. This structure makes the issue of inferring physical properties from dead time-affected observations fall into a class of problems common across many scientific disciplines. There is a growing number of methods to address them under the name of Simulation-Based Inference (SBI), aided by new developments in density estimation and statistical machine learning. In this paper, we introduce SBI as a principled way to infer variability properties from dead time-affected light curves. We use Sequential Neural Posterior Estimation to estimate the posterior probability for variability properties. We show that this method can recover variability parameters on simulated data even when dead time is variable, and present results of an application of this approach to NuSTAR observations of the galactic black hole X-ray binary GRS 1915+105.

21 pages, 24 figures; accepted for publication in MNRAS; code available at https://github.com/dhuppenkothen/DeadTimeSBI and data products at http://doi.org/10.5281/zenodo.4670161

References in corpus (25)

• scikit-image: Image processing in Python
• The frontier of simulation-based inference
• Fast likelihood-free cosmology with neural density estimators and active learning
• Searching the Entirety of Kepler Data. II. Occurrence Rate Estimates for FGK Stars
• No Time for Dead Time: Timing analysis of bright black hole binaries with NuSTAR
• Gaia DR2 reveals a star formation burst in the disc 2-3 Gyr ago
• Occurrence Rates of Planets Orbiting M Stars: Applying ABC to Kepler DR25, Gaia DR2, and 2MASS Data
• Flexible statistical inference for mechanistic models of neural dynamics
• Rapid Variability of Sgr A* across the Electromagnetic Spectrum
• astroABC: An Approximate Bayesian Computation Sequential Monte Carlo sampler for cosmological parameter estimation
• A Relationship Between Stellar Age and Spot Coverage
• Low frequency X-ray timing with Gaussian processes and reverberation in the radio-loud AGN 3C 120
• Analysis of coronal mass ejection flux rope signatures using 3DCORE and approximate Bayesian Computation
• Measurement of the B-band Galaxy Luminosity Function with Approximate Bayesian Computation
• Towards constraining warm dark matter with stellar streams through neural simulation-based inference
• Error formulae for the energy-dependent cross-spectrum
• Exploring the radial disc ionization profile of the black hole X-ray binary GRS 1915+105
• Systematic errors in strong gravitational lensing reconstructions, a numerical simulation perspective
• Constraining the recent star formation history of galaxies : an Approximate Bayesian Computation approach
• A unified framework for 21cm tomography sample generation and parameter inference with Progressively Growing GANs
• Reliability Correction is Key for Robust Kepler Occurrence Rates
• Testing the chemical homogeneity of chemically tagged dissolved birth clusters
• A Preferential Attachment Model for the Stellar Initial Mass Function
• Combining Transit and RV: A Synthesized Population Model
• HNPE: Leveraging Global Parameters for Neural Posterior Estimation

Cited by in corpus (12)

• Accelerated Bayesian SED Modeling using Amortized Neural Posterior Estimation
• $\mathrm{S}\mathrm{IM}\mathrm{BIG}$: Mock Challenge for a Forward Modeling Approach to Galaxy Clustering
• GLASS: Generator for Large Scale Structure
• Fast and robust Bayesian Inference using Gaussian Processes with GPry
• Searching for quasi-periodic oscillations in astrophysical transients using Gaussian processes
• Calibrating cosmological simulations with implicit likelihood inference using galaxy growth observables
• DIGS: Deep Inference of Galaxy Spectra with Neural Posterior Estimation
• Isolated pulsar population synthesis with simulation-based inference
• Fourier Domain
• Searching for quasi-periodicities in short transients: the curious case of GRB 230307A
• Deep inference of simulated strong lenses in ground-based surveys
• The Effects of Instrumental Deadtime on NICER Timing Products