Search NASA⌕ Search

Engineering topics

Peiris, Hiranya V.

Publications and source records attributed to Peiris, Hiranya V..

A robust estimator of mutual information for deep learning interpretability

Abstract We develop the use of mutual information (MI), a well-established metric in information theory, to interpret the inner workings of deep learning (DL) models. To accurately estimate MI from a finite number of samples, we present GMM-MI (pronounced ‘Jimmie’), an algorithm based on Gaussian mixture models that can be applied to both discrete and continuous settings. GMM-MI is computationally efficient, robust to the choice of hyperparameters and provides the uncertainty on the MI estimate due to the finite sample size. We extensively validate GMM-MI on toy data for which the ground truth MI is known, comparing its performance against established MI estimators. We then demonstrate the use of our MI estimator in the context of representation learning, working with synthetic data and physical datasets describing highly non-linear processes. We train DL models to encode high-dimensional data within a meaningful compressed (latent) representation, and use GMM-MI to quantify both the level of disentanglement between the latent variables, and their association with relevant physical quantities, thus unlocking the interpretability of the latent representation. We make GMM-MI publicly available in this GitHub repository.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Impact of Rubin Observatory Cadence Choices on Supernovae Photometric Classification

The Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST) will discover an unprecedented number of supernovae (SNe), making spectroscopic classification for all the events infeasible. LSST will thus rely on photometric classification, whose accuracy depends on the not-yet-finalized LSST observing strategy. In this work, we analyze the impact of cadence choices on classification performance using simulated multiband light curves. First, we simulate SNe with an LSST baseline cadence, a nonrolling cadence, and a presto-color cadence, which observes each sky location three times per night instead of twice. Each simulated data set includes a spectroscopically confirmed training set, which we augment to be representative of the test set as part of the classification pipeline. Then we use the photometric transient classification library snmachine to build classifiers. We find that the active region of the rolling cadence used in the baseline observing strategy yields a 25% improvement in classification performance relative to the background region. This improvement in performance in the actively rolling region is also associated with an increase of up to a factor of 2.7 in the number of cosmologically useful Type Ia SNe relative to the background region. However, adding a third visit per night as implemented in presto-color degrades classification performance due to more irregularly sampled light curves. Overall, our results establish desiderata on the observing cadence related to classification of full SNe light curves, which in turn impacts photometric SNe cosmology with LSST.

79 ASTRONOMY AND ASTROPHYSICS↗

PSTN-055: Survey Cadence Optimization Committee's Phase 2 Recommendations

This document describes the recommendations formulated by the Survey Cadence Optimization Committee (SCOC) of Vera C. Rubin Observatory for the observing strategy of the Legacy Survey of Space and Time (LSST). These "Phase 2" recommendations follow the initial "Phase 1" recommendations released in December 2021 as PSTN-053, and convey the current recommendations for the initial implementation of the LSST. These recommendations are based on the SCOC analysis of LSST simulations generated according to PSTN-053 and the feedback received from the scientific community on those simulations. These recommendations should be considered as the current guidance for implementing the LSST observing survey at the start of the survey, a strategy to be reviewed regularly over the survey lifetime, and are accompanied by a simulation, baseline_v3.0, which represents a specific implementation of these recommendations. While these recommendations dramatically narrow the possible LSST observing parameter choices, we emphasize that at this stage they do not yet represent a final or complete set of recommendations for the survey strategy, and several details are still to be refined. Furthermore, the strategy should be considered "living" and subject to evolution throughout the survey. The SCOC continues to engage with the community to finalize its recommendations. At the time of the release of this document, the final phases of construction of Rubin Observatory are ongoing with an expected first light in mid-2024. We expect the SCOC recommendations to be revised as needed, as knowledge of the system as built and the outcome of commissioning and Science Verification mature. This document describes the Phase 2 SCOC recommendations —highlighting elements that can be considered final and elements that remain to be finalized—, the process through which the SCOC has converged on these recommendations, and the working plan for the SCOC starting in January 2023 through the start of the survey to finalize the initial implementation of the Legacy Survey of Space and Time.

79 ASTRONOMY AND ASTROPHYSICS↗

Modelling populations of kilonovae

Abstract The 2017 detection of a kilonova coincident with gravitational-wave emission has identified neutron star mergers as the major source of the heaviest elements and dramatically constrained alternative theories of gravity. Observing a population of such sources has the potential to transform cosmology, nuclear physics, and astrophysics. However, with only one confident multi-messenger detection currently available, modelling the diversity of signals expected from such a population requires improved theoretical understanding. In particular, models that are quick to evaluate and are calibrated with more detailed multi-physics simulations are needed to design observational strategies for kilonovae detection and to obtain rapid-response interpretations of new observations. We use grey-opacity models to construct populations of kilonovae, spanning ejecta parameters predicted by numerical simulations. Our modelling focuses on wavelengths relevant for upcoming optical surveys, such as the Rubin Observatory Legacy Survey of Space and Time (LSST). In these simulations, we implement heating rates that are based on nuclear reaction network calculations. We create a Gaussian-process emulator for kilonova grey opacities, calibrated with detailed radiative transfer simulations. Using recent fits to numerical relativity simulations, we predict how the ejecta parameters from binary neutron star (BNS) mergers shape the population of kilonovae, accounting for the viewing-angle dependence. Our simulated population of BNS mergers produce peak i-band absolute magnitudes of −20 ≤ Mi ≤ −11. A comparison with detailed radiative transfer calculations indicates that further improvements are needed to accurately reproduce spectral shapes over the full light curve evolution.

79 ASTRONOMY AND ASTROPHYSICS↗

PSTN-053: Survey Cadence Optimization Committee's Phase 1 Recommendations

This document is authored by the Survey Cadence Optimization Committee (SCOC), guided by the community input, and it will serve as a roadmap for the concluding round of survey cadence simulations. These simulations, to be delivered during early calendar year 2022, will inform the phase 2 survey strategy recommendation which will define the baseline strategy for starting the Rubin Observatory's Legacy Survey of Space and Time (LSST). The phase 2 SCOC recommendation will be delivered to the Rubin Observatory Operations Director before the end of calendar year 2022.

79 ASTRONOMY AND ASTROPHYSICS↗

Considerations for Optimizing the Photometric Classification of Supernovae from the Rubin Observatory

The Vera C. Rubin Observatory will increase the number of observed supernovae (SNe) by an order of magnitude; however, it is impossible to spectroscopically confirm the class for all SNe discovered. Thus, photometric classification is crucial, but its accuracy depends on the not-yet-finalized observing strategy of Rubin Observatory's Legacy Survey of Space and Time (LSST). We quantitatively analyze the impact of the LSST observing strategy on SNe classification using simulated multiband light curves from the Photometric LSST Astronomical Time-Series Classification Challenge (PLAsTiCC). First, we augment the simulated training set to be representative of the photometric redshift distribution per SNe class, the cadence of observations, and the flux uncertainty distribution of the test set. Then we build a classifier using the photometric transient classification library snmachine, based on wavelet features obtained from Gaussian process fits, yielding a similar performance to the winning PLAsTiCC entry. We study the classification performance for SNe with different properties within a single simulated observing strategy. We find that season length is important, with light curves of 150 days yielding the highest performance. Cadence also has an important impact on SNe classification; events with median inter-night gap <3.5 days yield higher classification performance. Interestingly, we find that large gaps (>10 days) in light-curve observations do not impact performance if sufficient observations are available on either side, due to the effectiveness of the Gaussian process interpolation. This analysis is the first exploration of the impact of observing strategy on photometric SN classification with LSST.

79 ASTRONOMY AND ASTROPHYSICS↗

The large-scale environment of thermonuclear and core-collapse supernovae

ABSTRACT The new generation of wide-field time-domain surveys has made it feasible to study the clustering of supernova (SN) host galaxies in the large-scale structure (LSS) for the first time. We investigate the LSS environment of SN populations, using 106 dark matter density realisations with a resolution of ∼3.8 Mpc, constrained by the 2M+ + galaxy survey. We limit our analysis to redshift z < 0.036, using samples of 498 thermonuclear and 782 core-collapse SNe from the Zwicky Transient Facility’s Bright Transient Survey and Census of the Local Universe catalogues. We detect clustering of SNe with high significance; the observed clustering of the two SNe populations is consistent with each other. Further, the clustering of SN hosts is consistent with that of the Sloan Digital Sky Survey (SDSS) Baryon Oscillation Spectroscopic Survey DR12 spectroscopic galaxy sample in the same redshift range. Using a tidal shear classifier, we classify the LSS into voids, sheets, filaments, and knots. We find that both SNe and SDSS galaxies are predominantly found in sheets and filaments. SNe are significantly under-represented in voids and over-represented in knots compared to the volume fraction in these structures. This work opens the potential for using forthcoming wide-field deep SN surveys as a complementary LSS probe.

Astronomy & Astrophysics↗

An emulator for the Lyman-α forest in beyond-ΛCDM cosmologies

Interpreting observations of the Lyman-α forest flux power spectrum requires interpolation between a small number of expensive simulations. We present a Gaussian process emulator modelling the 1D flux power spectrum as a function of the amplitude and slope of the small-scale linear matter power spectrum, and the state of the intergalactic medium at the epoch of interest (2 < z < 4). This parameterisation enables the prediction of the flux power spectrum in extended cosmological models that are not explicitly included in the training set, eliminating the need to construct bespoke emulators for a number of extensions to ΛCDM. Our emulator is appropriate for cosmologies in which the linear matter power spectrum is described to percent level accuracy by just an amplitude and slope across the epoch of interest, and in the regime probed by eBOSS/DESI data. We demonstrate this for massive neutrino cosmologies, where the emulator is able to predict the flux power spectrum in a Σm ν = 0.3 eV neutrino cosmology to sub-percent accuracy, without including massive neutrinos in the training simulations. Further parameters would be required to describe models with sharp features in the linear power, such as warm or light axion dark matter. This work will facilitate the combination of upcoming DESI data with observations of the cosmic microwave background, to obtain constraints on neutrino mass and other extensions to ΛCDM cosmology.

79 ASTRONOMY AND ASTROPHYSICS↗