popclass: a python package for classifying microlensing events
Authors:
Greg Sallaberry,
Zofia Kaczmarek,
Peter McGill,
Scott E. Perkins,
William A. Dawson,
Caitlin G. Begbie
Abstract:
popclass is a python package that provides a flexible, probabilistic framework for classifying the lens of a gravitational microlensing event. popclass allows a user to match characteristics of a microlensing signal to a simulation of the Galaxy to calculate lens type probabilities for an event. Constraints on any microlensing signal characteristics and any Galactic model can be used. popclass com…
▽ More
popclass is a python package that provides a flexible, probabilistic framework for classifying the lens of a gravitational microlensing event. popclass allows a user to match characteristics of a microlensing signal to a simulation of the Galaxy to calculate lens type probabilities for an event. Constraints on any microlensing signal characteristics and any Galactic model can be used. popclass comes with an interface to common inference libraries for microlensing signal constraints, pre-loaded Galactic models, plotting functionality, and classification uncertainty quantification methods.
△ Less
Submitted 17 October, 2024;
originally announced October 2024.
A Scalable Gaussian Process Approach to Shear Mapping with MuyGPs
Authors:
Gregory Sallaberry,
Benjamin W. Priest,
Robert Armstrong,
Michael D. Schneider,
Amanda Muyskens,
Trevor Steil,
Keita Iwabuchi
Abstract:
Analysis of cosmic shear is an integral part of understanding structure growth across cosmic time, which in-turn provides us with information about the nature of dark energy. Conventional methods generate \emph{shear maps} from which we can infer the matter distribution in the universe. Current methods (e.g., Kaiser-Squires inversion) for generating these maps, however, are tricky to implement and…
▽ More
Analysis of cosmic shear is an integral part of understanding structure growth across cosmic time, which in-turn provides us with information about the nature of dark energy. Conventional methods generate \emph{shear maps} from which we can infer the matter distribution in the universe. Current methods (e.g., Kaiser-Squires inversion) for generating these maps, however, are tricky to implement and can introduce bias. Recent alternatives construct a spatial process prior for the lensing potential, which allows for inference of the convergence and shear parameters given lensing shear measurements. Realizing these spatial processes, however, scales cubically in the number of observations - an unacceptable expense as near-term surveys expect billions of correlated measurements. Therefore, we present a linearly-scaling shear map construction alternative using a scalable Gaussian Process (GP) prior called MuyGPs. MuyGPs avoids cubic scaling by conditioning interpolation on only nearest-neighbors and fits hyperparameters using batched leave-one-out cross validation. We use a suite of ray-tracing results from N-body simulations to demonstrate that our method can accurately interpolate shear maps, as well as recover the two-point and higher order correlations. We also show that we can perform these operations at the scale of billions of galaxies on high performance computing platforms.
△ Less
Submitted 30 September, 2024;
originally announced October 2024.