A Probabilistic Framework for Population Studies of the Solar Neighborhood: Application to SDSS-V and Gaia
Ilija Medan, Keivan G. Stassun, Zachary Way, Guy S. Stringfellow, Alexandre Roman-Lopes, Jiadong Li, Madeline Lucey, Andrew R. Casey, Bárbara Rojas-Ayala, Ricardo López-Valdivia, José G. Fernández-Trincado
astro-ph.SR, astro-ph.GA, astro-ph.IM
Submitted: 2026-07-21
Comments: 23 pages, 8 figures
Code: https://github.com/imedan/snc_sf
License: http://creativecommons.org/licenses/by-nc-nd/4.0/
The gist: Studies of the Solar Neighborhood require spectroscopic follow-up of stars identified in astrometric surveys to fully characterize their physical properties.
Terminology
Abstract
Studies of the Solar Neighborhood require spectroscopic follow-up of stars identified in astrometric surveys to fully characterize their physical properties. The SDSS-V Solar Neighborhood Census (SNC) is a dedicated program to observe stars within 100 pc. However, due to competing observing programs and fiber assignment constraints, the resulting sample carries severe and complex selection effects. A framework is presented for characterizing the selection function of the SDSS-V SNC relative to the Gaia Catalog of Nearby Stars (GCNS), along with a forward modeling method to infer the properties of stellar subpopulations across the GCNS-defined 100 pc sample. The selection function is based on a method that models the selection probability as a function of sky position, Gaia G magnitude, and BP-RP. The resulting detection probabilities faithfully reproduce the known survey planning logic. This work further introduces the concept of a "subpopulation probability" -- a grid of posterior estimates across the Hertzsprung-Russell (HR) diagram representing the likelihood that a GCNS member belongs to a given SDSS-V defined subpopulation. The framework is validated with a mock dataset and its scientific utility is demonstrated through two applications using data from the Data Release 19: mapping H alpha emission across the HR diagram and measuring the variation of stellar density with mass and metallicity. These results illustrate how statistically robust population studies can be conducted with an incomplete spectroscopic survey when the selection function is well characterized. The code is made publicly available with this work, which will serve as an important tool for future studies.
Sources
- robostrategy: Field and Target Assignment Optimization in the Sloan Digital Sky Survey V
- The No-U-Turn Sampler: Adaptively Setting Path Lengths in Hamiltonian Monte Carlo
- Procedures for Constraining Robotic Fiber Positioning for Highly Multiplexed Spectroscopic Surveys: The Case of FPS for SDSS-V
- Composable Effects for Flexible and Accelerated Probabilistic Programming in NumPyro
- Stellar Parameters of BOSS M dwarfs in SDSS-V DR19
- The Nineteenth Data Release of the Sloan Digital Sky Survey
Related papers
- HXI-DLA2: A Physics-Constrained Deep Learning Algorithm for the ASO-S Hard X-ray Imager
- Effect of Neutron Star Jets on Common Envelope Evolution
- Constraining the origin of magnetic white dwarfs
- JW-FD: A 15-Year Multimodal Dataset for Solar Flare Forecasting
- Phlegethon: a fully compressible magnetohydrodynamic code for simulations in stellar astrophysics
- Can MHD Oscillations Modulate Quasi-Periodic Plasma Release from Coronal Streamers?