From Blackbirds to Black Holes: Investigating Capture-Recapture Methods for Time Domain Astronomy

Abstract

In time domain astronomy, recurrent transients present a special problem: how to infer total populations from limited observations. Monitoring observations may give a biassed view of the underlying population due to limitations on observing time, visibility and instrumental sensitivity. A similar problem exists in the life sciences, where animal populations (such as migratory birds) or disease prevalence, must be estimated from sparse and incomplete data. The class of methods termed Capture-Recapture is used to reconstruct population estimates from time-series records of encounters with the study population. This paper investigates the performance of Capture-Recapture methods in astronomy via a series of numerical simulations. The Blackbirds code simulates monitoring of populations of transients, in this case accreting binary stars (neutron star or black hole accreting from a stellar companion) under a range of observing strategies. We first generate realistic light-curves for populations of binaries with contrasting orbital period distributions. These models are then randomly sampled at observing cadences typical of existing and planned monitoring surveys. The classical capture-recapture methods, Lincoln-Peterson, Schnabel estimators, related techniques, and newer methods implemented in the Rcapture package are compared. A general exponential model based on the radioactive decay law is introduced, and demonstrated to recover (at 95% confidence) the underlying population abundance and duty cycle, in a fraction of the observing visits (10-50%) required to discover all the sources in the simulation. Capture-Recapture is a promising addition to the toolbox of time domain astronomy, and methods implemented in R by the biostats community can be readily called from within Python.

arXiv commentAccepted to New Astronomy. 11 pages, 8 figures (refereed version prior to editorial process)

References in corpus (18)

• Be/X-ray binaries
• The MACHO Project LMC Microlensing Results from the First Two Years and the Nature of the Galactic Dark Halo
• Chandra Orion Ultradeep Project: Observations and Source Lists
• Master Robotic Net
• An Overabundance of Transient X-ray Binaries within 1 pc of the Galactic Center
• The XMM-Newton/Chandra monitoring campaign of the Galactic center region: Description of the program and preliminary results
• The Spectra and Variability of X-ray Sources in a Deep Chandra Observation of the Galactic Center
• Stellar coronal astronomy - a review
• A Long Look at the Be/X-Ray Binaries of the Small Magellanic Cloud
• The INTEGRAL Galactic bulge monitoring program: the first 1.5 years
• Chandra Multiwavelength Plane (ChaMPlane) Survey: an Introduction
• A Catalog of Transient X-ray Sources in M31
• The X-ray source content of the XMM-Newton Galactic Plane Survey
• Discovery of a Significant Magnetic CV Population in the Limiting Window
• Luminosity Functions and Point Source Properties from Multiple Chandra Observations of M81
• INTEGRAL deep observations of the Small Magellanic Cloud
• The Chandra ACIS Survey of M33 (ChASeM33): Transient X-ray Sources Discovered in M33
• Properties of X-ray Binaries in the Magellanic Clouds from RXTE and Chandra Observations