Selection Functions in Astronomical Data Modeling, with the Space Density of White Dwarfs as Worked Example

Abstract

Statistical studies of astronomical data sets, in particular of cataloged properties for discrete objects, are central to astrophysics. One cannot model those objects' population properties or incidences without a quantitative understanding of the conditions under which these objects ended up in a catalog or sample, the sample's selection function. As systematic and didactic introductions to this topic are scarce in the astrophysical literature, we aim to provide one, addressing generically the following questions: What is a selection function? What arguments $\vec{q}$ should a selection function depend on? Over what domain must a selection function be defined? What approximations and simplifications can be made? And, how is a selection function used in `modelling'? We argue that volume-complete samples, with the volume drastically curtailed by the faintest objects, reflect a highly sub-optimal selection function that needlessly reduces the number of bright and usually rare objects in the sample. We illustrate these points by a worked example, deriving the space density of white dwarfs (WD) in the Galactic neighbourhood as a function of their luminosity and Gaia color, ${\Phi }_0(M_G,B-R)$ in [mag${}^{-2}$pc${}^{-3}$]. We construct a sample of $10^5$ presumed WDs through straightforward selection cuts on the Gaia EDR3 catalog, in magnitude, color, parallax, and astrometric fidelity $\vec{q}=(m_G,B-R,\varpi ,p_{af})$. We then combine a simple model for ${\Phi }_0$ with the effective survey volume derived from this selection function $S_C(\vec{q})$ to derive a detailed and robust estimate of ${\Phi }_0(M_G,B-R)$. This resulting white dwarf luminosity-color function ${\Phi }_0(M_G,B-R)$ differs dramatically from the initial number density distribution in the luminosity-color plane: by orders of magnitude in density and by four magnitudes in density peak location.

References in corpus (8)

• The Gaia mission
• Gaia-2MASS 3D maps of Galactic interstellar dust within 3 kpc
• Exoplanet population inference and the abundance of Earth analogs from noisy, incomplete catalogs
• $Gaia$ white dwarfs within 40 pc II: the volume-limited northern hemisphere sample
• Completeness of the Gaia-verse II: what are the odds that a star is missing from Gaia DR2?
• Completeness of the Gaia-verse I: when and where were Gaia's eyes on the sky during DR2?
• Completeness of the Gaia-verse IV: The Astrometry Spread Function of Gaia DR2
• Completeness of the Gaia-verse III: using hidden states to infer gaps, detection efficiencies and the scanning law from the DR2 light curves