Searching for Narrow Emission Lines in X-ray Spectra: Computation and Methods

Abstract

The detection and quantification of narrow emission lines in X-ray spectra is a challenging statistical task. The Poisson nature of the photon counts leads to local random fluctuations in the observed spectrum that often results in excess emission in a narrow band of energy resembling a weak narrow line. From a formal statistical perspective, this leads to a (sometimes highly) multimodal likelihood. Many standard statistical procedures are based on (asymptotic) Gaussian approximations to the likelihood and simply cannot be used in such settings. Bayesian methods offer a more direct paradigm for accounting for such complicated likelihood functions but even here multimodal likelihoods pose significant computational challenges. The new Markov chain Monte Carlo (MCMC) methods developed in 2008 by van Dyk and Park, however, are able to fully explore the complex posterior distribution of the location of a narrow line, and thus provide valid statistical inference. Even with these computational tools, standard statistical quantities such as means and standard deviations cannot adequately summarize inference and standard testing procedures cannot be used to test for emission lines. In this paper, we use new efficient MCMC algorithms to fit the location of narrow emission lines, we develop new statistical strategies for summarizing highly multimodal distributions and quantifying valid statistical inference, and we extend the method of posterior predictive p-values proposed by Protassov et al. (2002) to test for the presence of narrow emission lines in X-ray spectra. We illustrate and validate our methods using simulation studies and apply them to the Chandra observations of the high redshift quasar PG1634+706.

43 pages, 11 figures; accepted for publication in ApJ

References in corpus (22)

• An Overview of the Performance and Scientific Results from the Chandra X-Ray Observatory (CXO)
• Statistics: Handle with Care, Detecting Multiple Model Components with the Likelihood Ratio Test
• Going with the flow: can the base of jets subsume the role of compact accretion disk coronae?
• Sherpa: a Mission-Independent Data Analysis Application
• Relativistic X-ray Lines from the Inner Accretion Disks Around Black Holes
• Chandra Detects Relativistic Broad Absorption Lines from APM 08279+5255
• X-ray Observations of Optically-Selected, Radio-Quiet Quasars I: The ASCA results
• X-ray Spectroscopy of QSOs with Broad Ultraviolet Absorption Lines
• An absorption origin for the X-ray spectral variability of MCG-6-30-15
• On the dependence of the iron K-line profiles with luminosity in Active Galactic Nuclei
• Physical Diagnostics from a Narrow Fe-Kalpha Emission Line Detected by Chandra in the Seyfert 1 Galaxy NGC 5548
• XMM-Newton spectroscopy of high-redshift QSOs
• Analysis of Energy Spectra with Low Photon Counts via Bayesian Posterior Simulation
• The XMM-Newton view of PG quasars: II. Properties of the Fe K-alpha line
• On the X-ray Baldwin effect for narrow Fe Kalpha emission line
• Statistics of relativistically broadened Fe K-alpha lines in AGN
• Variable iron-line emission near the black hole of Markarian 766
• Suzaku Observations of Iron Lines and Reflection in AGN
• Is the X-ray spectrum of the narrow emission line QSO PG1211+143 defined by its energetic outflow?
• High redshift radio quiet quasars -- exploring parameter space of accretion models. Part II: patchy corona model
• Energy distribution of individual quasars from far-UV to X-rays: I. Intrinsic UV hardness and dust opacities
• A Short Review of Relativistic Iron Lines from Stellar-Mass Black Holes

Cited by in corpus (13)

• A Bayesian test for periodic signals in red noise
• On Computing Upper Limits to Source Intensities
• The rapid X-ray variability of NGC 4051
• Quasi-Periodic Oscillations and broadband variability in short magnetar bursts
• Accounting for Calibration Uncertainties in X-ray Analysis: Effective Areas in Spectral Fitting
• Metropolis-Hastings within Partially Collapsed Gibbs Samplers
• Cross-Fertilizing Strategies for Better EM Mountain Climbing and DA Field Exploration: A Graphical Guide Book
• Disentangling Overlapping Astronomical Sources using Spatial and Spectral Information
• A Bayesian Analysis of the Correlations Among Sunspot Cycles
• A Hierarchical Model for the Ages of Galactic Halo White Dwarfs
• Data-Space Validation of High-Dimensional Models by Comparing Sample Quantiles
• A novel approach to detect line emission under high background in high-resolution X-ray spectra
• X-ray Properties of the Luminous Quasar PG 1634+706 at z = 1.337 from SRG and XMM-Newton Data