-
Closely-Spaced Object Classification Using MuyGPyS
Authors:
Kerianne Pruett,
Nathan McNaughton,
Michael Schneider
Abstract:
Accurately detecting rendezvous and proximity operations (RPO) is crucial for understanding how objects are behaving in the space domain. However, detecting closely-spaced objects (CSO) is challenging for ground-based optical space domain awareness (SDA) algorithms as two objects close together along the line-of-sight can appear blended as a single object within the point-spread function (PSF) of…
▽ More
Accurately detecting rendezvous and proximity operations (RPO) is crucial for understanding how objects are behaving in the space domain. However, detecting closely-spaced objects (CSO) is challenging for ground-based optical space domain awareness (SDA) algorithms as two objects close together along the line-of-sight can appear blended as a single object within the point-spread function (PSF) of the optical system. Traditional machine learning methods can be useful for differentiating between singular objects and closely-spaced objects, but many methods require large training sample sizes or high signal-to-noise conditions. The quality and quantity of realistic data make probabilistic classification methods a superior approach, as they are better suited to handle these data inadequacies. We present CSO classification results using the Gaussian process python package, MuyGPyS, and examine classification accuracy as a function of angular separation and magnitude difference between the simulated satellites. This orbit-independent analysis is done on highly accurate simulated SDA images that emulate realistic ground-based commercial-of-the-shelf (COTS) optical sensor observations of CSOs. We find that MuyGPyS outperforms traditional machine learning methods, especially under more challenging circumstances.
△ Less
Submitted 17 November, 2023;
originally announced November 2023.
-
Disentangling the Black Hole Mass Spectrum with Photometric Microlensing Surveys
Authors:
Scott Ellis Perkins,
Peter McGill,
William Dawson,
Natasha S. Abrams,
Casey Y. Lam,
Ming-Feng Ho,
Jessica R. Lu,
Simeon Bird,
Kerianne Pruett,
Nathan Golovich,
George Chapline
Abstract:
From the formation mechanisms of stars and compact objects to nuclear physics, modern astronomy frequently leverages surveys to understand populations of objects to answer fundamental questions. The population of dark and isolated compact objects in the Galaxy contains critical information related to many of these topics, but is only practically accessible via gravitational microlensing. However,…
▽ More
From the formation mechanisms of stars and compact objects to nuclear physics, modern astronomy frequently leverages surveys to understand populations of objects to answer fundamental questions. The population of dark and isolated compact objects in the Galaxy contains critical information related to many of these topics, but is only practically accessible via gravitational microlensing. However, photometric microlensing observables are degenerate for different types of lenses, and one can seldom classify an event as involving either a compact object or stellar lens on its own. To address this difficulty, we apply a Bayesian framework that treats lens type probabilistically and jointly with a lens population model. This method allows lens population characteristics to be inferred despite intrinsic uncertainty in the lens-class of any single event. We investigate this method's effectiveness on a simulated ground-based photometric survey in the context of characterizing a hypothetical population of primordial black holes (PBHs) with an average mass of $30 M_{\odot}$. On simulated data, our method outperforms current black hole (BH) lens identification pipelines and characterizes different subpopulations of lenses while jointly constraining the PBH contribution to dark matter to ${\approx}25$\%. Key to robust inference, our method can marginalize over population model uncertainty. We find the lower mass cutoff for stellar origin BHs, a key observable in understanding the BH mass gap, particularly difficult to infer in our simulations. This work lays the foundation for cutting-edge PBH abundance constraints to be extracted from current photometric microlensing surveys.
△ Less
Submitted 5 October, 2023;
originally announced October 2023.
-
Markov Chain Monte Carlo for Bayesian Parametric Galaxy Modeling in LSST
Authors:
James J. Buchanan,
Michael D. Schneider,
Kerianne Pruett,
Robert E. Armstrong
Abstract:
We apply Markov Chain Monte Carlo (MCMC) to the problem of parametric galaxy modeling, estimating posterior distributions of galaxy properties such as ellipticity and brightness for more than 100,000 images of galaxies taken from DC2, a simulated telescope survey resembling the upcoming Rubin Observatory Legacy Survey of Space and Time (LSST). We use a physically informed prior and apply selection…
▽ More
We apply Markov Chain Monte Carlo (MCMC) to the problem of parametric galaxy modeling, estimating posterior distributions of galaxy properties such as ellipticity and brightness for more than 100,000 images of galaxies taken from DC2, a simulated telescope survey resembling the upcoming Rubin Observatory Legacy Survey of Space and Time (LSST). We use a physically informed prior and apply selection corrections to the likelihood. The resulting posterior samples enable rigorous probabilistic inference of galaxy model parameters and their uncertainties. These posteriors are one key ingredient in a fully probabilistic description of galaxy catalogs, which can ultimately enable a refined Bayesian estimate of cosmological parameters. We systematically examine the reliability of the posterior mean as a point estimator of galaxy parameters, and of the posterior width as a measure of uncertainty, under some common modeling approximations. We implement the probabilistic modeling and MCMC inference using the JIF (Joint Image Framework) tool, which we make freely available online.
△ Less
Submitted 19 September, 2023;
originally announced September 2023.
-
MEGASIM: Distribution and Detection of Earth Trojan Asteroids
Authors:
Travis Yeager,
Nathan Golovich,
Kerianne Pruett
Abstract:
Using N-body simulation results from the MEGASIM dataset, we present spatial distributions of Earth Trojan Asteroids and assess the detectability of the population in current and next-generation ground-based astronomical surveys (Yeager & Golovich 2022). Our high-fidelity Earth Trojan Asteroid (ETA) distribution maps show never-before-seen high-resolution spatial features that evolve over timescal…
▽ More
Using N-body simulation results from the MEGASIM dataset, we present spatial distributions of Earth Trojan Asteroids and assess the detectability of the population in current and next-generation ground-based astronomical surveys (Yeager & Golovich 2022). Our high-fidelity Earth Trojan Asteroid (ETA) distribution maps show never-before-seen high-resolution spatial features that evolve over timescales up to 1 Gyr. The simulation was synchronized to start times and timelines of two observational astronomy surveys, 1) the Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST) and 2) the Zwicky Transient Facility (ZTF). We calculate upper limits for the number of ETAs potentially observable with both the ZTF and LSST surveys. Due to the Yarkovsky Effect, we find no stable ETAs on billion year timescales likely to be detected by any ETA survey, as no C-type or S-type ETAs (with $H<22$ and $H<24$, respectively) are likely to be stable on billion year timescales, and ETAs large enough to remain stable on billion year timescales are very rare relative to the rest of the ETA population. We find that a twilight ETA survey will not drastically increase the likelihood of individual ETA detection, but would provide orders of magnitude more observations of select ETA populations. The null detection to date from ZTF restricts the potential ETA population to hundreds of objects larger than 100 meters (at $H\approx22$, Lifset et al. 2021), while a null detection by LSST will further restrict the ETA population to tens of objects larger than 100 meters.
△ Less
Submitted 7 November, 2023; v1 submitted 21 February, 2023;
originally announced February 2023.
-
Primordial Black Hole Dark Matter Simulations Using PopSyCLE
Authors:
K. Pruett,
W. Dawson,
M. S. Medford,
C. Lam,
J. R. Lu,
S. Perkins,
P. McGill,
N. Golovich,
G. Chapline
Abstract:
Primordial black holes (PBHs), theorized to have originated in the early universe, are speculated to be a viable form of dark matter. If they exist, they should be detectable through photometric and astrometric signals resulting from gravitational microlensing of stars in the Milky Way. Population Synthesis for Compact-object Lensing Events, or PopSyCLE, is a simulation code that enables users to…
▽ More
Primordial black holes (PBHs), theorized to have originated in the early universe, are speculated to be a viable form of dark matter. If they exist, they should be detectable through photometric and astrometric signals resulting from gravitational microlensing of stars in the Milky Way. Population Synthesis for Compact-object Lensing Events, or PopSyCLE, is a simulation code that enables users to simulate microlensing surveys, and is the first of its kind to include both photometric and astrometric microlensing effects, which are important for potential PBH detection and characterization. To estimate the number of observable PBH microlensing events we modify PopSyCLE to include a dark matter halo consisting of PBHs. We detail our PBH population model, and demonstrate our PopSyCLE + PBH results through simulations of the OGLE-IV and Roman microlensing surveys. We provide a proof-of-concept analysis for adding PBHs into PopSyCLE, and thus include many simplifying assumptions, such as $f_{\text{DM}}$, the fraction of dark matter composed of PBHs, and $\bar{m}_{\text{PBH}}$, mean PBH mass. Assuming $\bar{m}_{\text{PBH}}=30$ $M_{\odot}$, we find $\sim$ 3.6$f_{\text{DM}}$ times as many PBH microlensing events than stellar evolved black hole events, a PBH average peak Einstein crossing time of $\sim$ 91.5 days, estimate on order of $10^2f_{\text{DM}}$ PBH events within the 8 year OGLE-IV results, and estimate Roman to detect $\sim$ 1,000$f_{\text{DM}}$ PBH microlensing events throughout its planned microlensing survey.
△ Less
Submitted 8 April, 2024; v1 submitted 12 November, 2022;
originally announced November 2022.
-
A Reanalysis of Public Galactic Bulge Gravitational Microlensing Events from OGLE-III and IV
Authors:
Nathan Golovich,
William A. Dawson,
Fran Bartolić,
Casey Y. Lam,
Jessica R. Lu,
Michael S. Medford,
Michael D. Schneider,
George Chapline,
Edward F. Schlafly,
Alex Drlica-Wagner,
Kerianne Pruett
Abstract:
Modern surveys of gravitational microlensing events have progressed to detecting thousands per year. Surveys are capable of probing Galactic structure, stellar evolution, lens populations, black hole physics, and the nature of dark matter. One of the key avenues for doing this is studying the microlensing Einstein radius crossing time distribution ($t_E$). However, systematics in individual light…
▽ More
Modern surveys of gravitational microlensing events have progressed to detecting thousands per year. Surveys are capable of probing Galactic structure, stellar evolution, lens populations, black hole physics, and the nature of dark matter. One of the key avenues for doing this is studying the microlensing Einstein radius crossing time distribution ($t_E$). However, systematics in individual light curves as well as over-simplistic modeling can lead to biased results. To address this, we developed a model to simultaneously handle the microlensing parallax due to Earth's motion, systematic instrumental effects, and unlensed stellar variability with a Gaussian Process model. We used light curves for nearly 10,000 OGLE-III and IV Milky Way bulge microlensing events and fit each with our model. We also developed a forward model approach to infer the timescale distribution by forward modeling from the data rather than using point estimates from individual events. We find that modeling the variability in the baseline removes a source of significant bias in individual events, and previous analyses over-estimated the number of long timescale ($t_E>100$ days) events due to their over simplistic models ignoring parallax effects and stellar variability. We use our fits to identify hundreds of events that are likely black holes.
△ Less
Submitted 16 September, 2020;
originally announced September 2020.
-
The Complete Calibration of the Color-Redshift Relation (C3R2) Survey: Analysis and Data Release 2
Authors:
Daniel C. Masters,
Daniel K. Stern,
Judith G. Cohen,
Peter L. Capak,
S. Adam Stanford,
Nina Hernitschek,
Audrey Galametz,
Iary Davidzon,
Jason D. Rhodes,
Dave Sanders,
Bahram Mobasher,
Francisco Castander,
Kerianne Pruett,
Sotiria Fotopoulou
Abstract:
The Complete Calibration of the Color-Redshift Relation (C3R2) survey is a multi-institution, multi-instrument survey that aims to map the empirical relation of galaxy color to redshift to i~24.5 (AB), thereby providing a firm foundation for weak lensing cosmology with the Stage IV dark energy missions Euclid and WFIRST. Here we present 3171 new spectroscopic redshifts obtained in the 2016B and 20…
▽ More
The Complete Calibration of the Color-Redshift Relation (C3R2) survey is a multi-institution, multi-instrument survey that aims to map the empirical relation of galaxy color to redshift to i~24.5 (AB), thereby providing a firm foundation for weak lensing cosmology with the Stage IV dark energy missions Euclid and WFIRST. Here we present 3171 new spectroscopic redshifts obtained in the 2016B and 2017A semesters with a combination of DEIMOS, LRIS, and MOSFIRE on the Keck telescopes. The observations come from all of the Keck partners: Caltech, NASA, the University of Hawaii, and the University of California. Combined with the 1283 redshifts published in DR1, the C3R2 survey has now obtained and published 4454 high quality galaxy redshifts. We discuss updates to the survey design and provide a catalog of photometric and spectroscopic data. Initial tests of the calibration method performance are given, indicating that the sample, once completed and combined with extensive data collected by other spectroscopic surveys, should allow us to meet the cosmology requirements for Euclid, and make significant headway toward solving the problem for WFIRST. We use the full spectroscopic sample to demonstrate that galaxy brightness is weakly correlated with redshift once a galaxy is localized in the Euclid or WFIRST color space, with potentially important implications for the spectroscopy needed to calibrate redshifts for faint WFIRST and LSST sources.
△ Less
Submitted 12 April, 2019;
originally announced April 2019.