Identification of high-energy astrophysical point sources via hierarchical Bayesian nonparametric clustering

The light we receive from distant astrophysical objects carries information about their origins and the physical mechanisms that power them. The study of these signals, however, is complicated by the fact that observations are often a mixture of the light emitted by multiple localized sources situated in a spatially-varying background. A general algorithm to achieve robust and accurate source identification in this case remains an open question in astrophysics. This paper focuses on high-energy light (such as X-rays and γ-rays), for which observatories can detect individual photons (quanta of light), measuring their incoming direction, arrival time, and energy. Our proposed Bayesian methodology uses both the spatial and energy information to identify point sources, that is, separate them from the spatiallyvarying background, to estimate their number, and to compute the posterior probabilities that each photon originated from each identified source. This is accomplished via a Dirichlet process mixture while the background is simultaneously reconstructed via a flexible Bayesian nonparametric model based on B-splines. Our proposed method is validated with a suite of simulation studies and illustrated with an application to a complex region of the sky observed by the Fermi Gamma-ray Space Telescope. 2 Sottosanti et al. (2021)

[1]  V. Dose,et al.  Background-Source separation in astronomical images with Bayesian Probability Theory , 2009, 0903.2342.

[2]  R Core Team,et al.  R: A language and environment for statistical computing. , 2014 .

[3]  C. Kouveliotou,et al.  A Chandra Search for Coronal X-Rays from the Cool White Dwarf GD 356 , 2006, astro-ph/0609585.

[4]  P. Müller,et al.  A Bayesian mixture model for differential gene expression , 2005 .

[5]  V. Kashyap,et al.  Determining X-Ray Source Intensity and Confidence Bounds in Crowded Fields , 2014, 1410.2564.

[6]  T. Ensslin,et al.  Denoising, deconvolving, and decomposing photon observations - Derivation of the D3PO algorithm , 2015 .

[7]  Vinay L. Kashyap,et al.  Bayesian Estimation of Hardness Ratios: Modeling and Computations , 2006 .

[8]  C. Biller.  Adaptive Bayesian Regression Splines in Semiparametric Generalized Linear Models , 2000 .

[9]  Ryan E. Keeley,et al.  Strong constraints on thermal relic dark matter from Fermi-LAT observations of the Galactic Center , 2020, 2003.10416.

[10]  Radford M. Neal.  Markov Chain Sampling Methods for Dirichlet Process Mixture Models , 2000 .

[11]  Statistics, Handle with Care: Detecting Multiple Model Components with the Likelihood Ratio Test , 2002, astro-ph/0201547.

[12]  Sergey E. Koposov,et al.  Snake in the Clouds: a new nearby dwarf galaxy in the Magellanic bridge* , 2018, Monthly Notices of the Royal Astronomical Society.

[13]  D. Blackwell.  Discreteness of Ferguson Selections , 1973 .

[14]  The Fermi-LAT Collaboration.  Fermi Large Area Telescope Third Source Catalog , 2015, 1501.02003.

[15]  Markus Ackermann,et al.  Fermi Large Area Telescopeによる宇宙線電子非等方性の探査 , 2010 .

[16]  Sylvia Frühwirth-Schnatter,et al.  Dealing with Label Switching under Model Uncertainty , 2011 .

[17]  David E. Jones,et al.  DISENTANGLING OVERLAPPING ASTRONOMICAL SOURCES USING SPATIAL AND SPECTRAL INFORMATION , 2014, 1411.7447.

[18]  R. Kass,et al.  Bayesian curve-fitting with free-knot splines , 2001 .

[19]  D. Thompson,et al.  DEVELOPMENT OF THE MODEL OF GALACTIC INTERSTELLAR EMISSION FOR STANDARD POINT-SOURCE ANALYSIS OF FERMI LARGE AREA TELESCOPE DATA , 2016, 1602.07246.

[20]  B -splines, hypergeometric functions, and Dirichlet averages , 1991 .

[21]  Sergey E. Koposov,et al.  The hidden giant: discovery of an enormous Galactic dwarf satellite in Gaia DR2 , 2018, Monthly Notices of the Royal Astronomical Society.

[22]  P. Green,et al.  Corrigendum: On Bayesian analysis of mixtures with an unknown number of components , 1997 .

[23]  Fernando A. Quintana,et al.  Bayesian Nonparametric Data Analysis , 2015 .

[24]  M. Escobar.  Estimating Normal Means with a Dirichlet Process Prior , 1994 .

[25]  D. Burrows,et al.  Determination of Confidence Limits for Experiments with Low Numbers of Counts , 1991 .

[26]  Göran Kauermann,et al.  Density estimation and comparison with a penalized mixture approach , 2012, Comput. Stat..

[27]  Adrian F. M. Smith,et al.  Automatic Bayesian curve fitting , 1998 .

[28]  C. R. Deboor,et al.  A practical guide to splines , 1978 .

[29]  M. Knoetig.  SIGNAL DISCOVERY, LIMITS, AND UNCERTAINTIES WITH SPARSE ON/OFF MEASUREMENTS: AN OBJECTIVE BAYESIAN ANALYSIS , 2014, 1406.2922.

[30]  Analysis of energy spectra with low photon counts via Bayesian posterior simulation , 2001, astro-ph/0008170.

[31]  S. Oliver,et al.  Bayesian Methods of Astronomical Source Extraction , 2005, astro-ph/0512597.

[32]  David J. Spiegelhalter,et al.  Local computations with probabilities on graphical structures and their application to expert systems , 1990 .

[33]  F. Feroz,et al.  Multimodal nested sampling: an efficient and robust alternative to Markov Chain Monte Carlo methods for astronomical data analyses , 2007, 0704.3704.

[34]  S. MacEachern,et al.  A semiparametric Bayesian model for randomised block designs , 1996 .

[35]  A. Brazzale,et al.  Bayesian mixture modelling of the high-energy photon counts collected by the Fermi Large Area Telescope , 2020, Statistical Modelling.

[36]  D. V. Dyk,et al.  DETECTING UNSPECIFIED STRUCTURE IN LOW-COUNT IMAGES , 2015, 1510.04662.

[37]  Jun S. Liu,et al.  The Collapsed Gibbs Sampler in Bayesian Computations with Applications to a Gene Regulation Problem , 1994 .

[38]  S. MacEachern,et al.  Estimating mixture of dirichlet process models , 1998 .

[39]  D. L. Bertsch,et al.  The Likelihood Analysis of EGRET Data , 1996 .

[40]  R. Trotta.  Bayesian Methods in Cosmology , 2017, 1701.01467.

[41]  T. Ensslin,et al.  The denoised, deconvolved, and decomposed Fermi γ-ray sky - An application of the D3PO algorithm , 2014, 1410.4562.

[42]  T. Ferguson.  A Bayesian Analysis of Some Nonparametric Problems , 1973 .

[43]  M. Escobar,et al.  Bayesian Density Estimation and Inference Using Mixtures , 1995 .

[44]  D. Thompson,et al.  PRECISE γ-RAY TIMING AND RADIO OBSERVATIONS OF 17 FERMI γ-RAY PULSARS , 2010, 1011.2468.

[45]  D. Ruppert,et al.  Bayesian adaptive B-spline estimation in proportional hazards frailty models , 2010 .

[46]  D. Finkbeiner,et al.  Inference of Unresolved Point Sources at High Galactic Latitudes Using Probabilistic Catalogs , 2016, 1607.04637.

[47]  F. Collaboration,et al.  The Fermi Large Area Telescope On Orbit: Event Classification, Instrument Response Functions, and Calibration , 2012, 1206.1896.