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.

arXiv comment21 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 (41)

• emcee: The MCMC Hammer
• scikit-image: Image processing in Python
• The Nuclear Spectroscopic Telescope Array (NuSTAR) Mission
• The Fermi Gamma-Ray Burst Monitor
• Towards a unified model for black hole X-ray binary jets
• The frontier of simulation-based inference
• The flux-dependent amplitude of broadband noise variability in X-ray binaries and active galaxies
• Stingray: A Modern Python Library For Spectral Timing
• Fast likelihood-free cosmology with neural density estimators and active learning
• Fast $\epsilon$-free Inference of Simulation Models with Bayesian Conditional Density Estimation
• 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
• Living on a Flare: Relativistic Reflection in V404 Cyg Observed by NuSTAR During its Summer 2015 Outburst
• 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
• Approximate Bayesian Computation in Large Scale Structure: constraining the galaxy-halo connection
• The multiplicity distribution of Kepler's exoplanets
• ABC-CDE: Towards Approximate Bayesian Computation with Complex High-Dimensional Data and Limited Simulations
• No Time for Dead Time: Use the Fourier Amplitude Differences to Normalize Dead-time-affected Periodograms
• Low frequency X-ray timing with Gaussian processes and reverberation in the radio-loud AGN 3C 120
• The redshift distribution of cosmological samples: a forward modeling approach
• Analysis of coronal mass ejection flux rope signatures using 3DCORE and approximate Bayesian Computation
• Towards constraining warm dark matter with stellar streams through neural simulation-based inference
• Measurement of the B-band Galaxy Luminosity Function with Approximate Bayesian Computation
• Error formulae for the energy-dependent cross-spectrum
• Monte Carlo Control Loops for cosmic shear cosmology with DES Year 1
• 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
• Primordial power spectrum and cosmology from black-box galaxy surveys
• On the Statistical Properties of Cospectra
• 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