Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “LSST”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Modelling the Milky Way – I. Method and first results fitting the thick disc and halo with DES-Y3 data

ABSTRACT We present a technique to fit the stellar components of the Galaxy by comparing Hess Diagrams (HDs) generated from trilegal models to real data. We apply this technique, which we call mwfitting, to photometric data from the first 3 yr of the Dark Energy Survey (DES). After removing regions containing known resolved stellar systems such as globular clusters, dwarf galaxies, nearby galaxies, the Large Magellanic Cloud, and the Sagittarius Stream, our main sample spans a total area of ∼2300 deg2. We further explore a smaller subset (∼1300 deg2) that excludes all regions with known stellar streams and stellar overdensities. Validation tests on synthetic data possessing similar properties to the DES data show that the method is able to recover input parameters with a precision better than 3 per cent. We fit the DES data with an exponential thick disc model and an oblate double power-law halo model. We find that the best-fitting thick disc model has radial and vertical scale heights of 2.67 ± 0.09 kpc and 925 ± 40 pc, respectively. The stellar halo is fit with a broken power-law density profile with an oblateness of 0.75 ± 0.01, an inner index of 1.82 ± 0.08, an outer index of 4.14 ± 0.05, and a break at 18.52 ± 0.27 kpc from the Galactic centre. Several previously discovered stellar overdensities are recovered in the residual stellar density map, showing the reliability of mwfitting in determining the Galactic components. Simulations made with the best-fitting parameters are a promising way to predict Milky Way star counts for surveys such as the LSST and Euclid.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Reliable photometric membership (RPM) of galaxies in clusters – I. A machine learning method and its performance in the local universe

ABSTRACT We introduce a new method to determine galaxy cluster membership based solely on photometric properties. We adopt a machine learning approach to recover a cluster membership probability from galaxy photometric parameters and finally derive a membership classification. After testing several machine learning techniques (such as stochastic gradient boosting, model averaged neural network and k-nearest neighbours), we found the support vector machine algorithm to perform better when applied to our data. Our training and validation data are from the Sloan Digital Sky Survey main sample. Hence, to be complete to $M_r^* + 3$, we limit our work to 30 clusters with $z$phot-cl ≤ 0.045. Masses (M200) are larger than $\sim 0.6\times 10^{14} \, \mathrm{M}_{\odot }$ (most above $3\times 10^{14} \, \mathrm{M}_{\odot }$). Our results are derived taking in account all galaxies in the line of sight of each cluster, with no photometric redshift cuts or background corrections. Our method is non-parametric, making no assumptions on the number density or luminosity profiles of galaxies in clusters. Our approach delivers extremely accurate results (completeness, C $\sim 92{\rm{ per\ cent}}$ and purity, P $\sim 87{\rm{ per\ cent}}$) within R200, so that we named our code reliable photometric membership. We discuss possible dependencies on magnitude, colour, and cluster mass. Finally, we present some applications of our method, stressing its impact to galaxy evolution and cosmological studies based on future large-scale surveys, such as eROSITA, EUCLID, and LSST.

Lopes, Paulo A. A.↗
DOE OSTI2020

Cross-correlating Planck with VST ATLAS LRGs: a new test for the ISW effect in the Southern hemisphere

ABSTRACT The integrated Sachs–Wolfe (ISW) effect probes the late-time expansion history of the Universe, offering direct constraints on dark energy. Here, we present our measurements of the ISW signal at redshifts of $\bar{z}=0.35$, 0.55, and 0.68, using the cross-correlation of the Planck cosmic microwave background temperature map with ∼0.5 million luminous red galaxies (LRGs) selected from the VST ATLAS survey. We then combine these with previous measurements based on WMAP and similar SDSS LRG samples, providing a total sample of ∼2.1 million LRGs covering ∼12 000 deg2 of sky. At $\bar{z}=0.35$ and $\bar{z}=0.55$, we detect the ISW signal at 1.2σ and 2.3σ (or 2.6σ combined), in agreement with the predictions of lambda cold dark matter (ΛCDM). We verify these results by repeating the measurements using the BOSS LOWZ and CMASS, spectroscopically confirmed LRG samples. We also detect the ISW effect in three magnitude limited ATLAS + SDSS galaxy samples extending to z ≈ 0.4 at ∼2σ per sample. However, we do not detect the ISW signal at $\bar{z}=0.68$ when combining the ATLAS and SDSS results. Further tests using spectroscopically confirmed eBOSS LRGs at this redshift remain inconclusive due to the current low sky coverage of the survey. If the ISW signal is shown to be redshift dependent in a manner inconsistent with the predictions of ΛCDM, it could open the door to alternative theories such as modified gravity. It is therefore important to repeat the high-redshift ISW measurement using the completed eBOSS sample, as well as deeper upcoming surveys such as DESI and LSST.

Ansarinejad, Behzad↗
DOE OSTI2020

First cosmology results using type Ia supernovae from the Dark Energy Survey: the effect of host galaxy properties on supernova luminosity

ABSTRACT We present improved photometric measurements for the host galaxies of 206 spectroscopically confirmed type Ia supernovae discovered by the Dark Energy Survey Supernova Program (DES-SN) and used in the first DES-SN cosmological analysis. For the DES-SN sample, when considering a 5D (z, x1, c, α, β) bias correction, we find evidence of a Hubble residual ‘mass step’, where SNe Ia in high-mass galaxies (>1010M⊙) are intrinsically more luminous (after correction) than their low-mass counterparts by $\gamma =0.040\pm 0.019$ mag. This value is larger by 0.031 mag than the value found in the first DES-SN cosmological analysis. This difference is due to a combination of updated photometric measurements and improved star formation histories and is not from host-galaxy misidentification. When using a 1D (redshift-only) bias correction the inferred mass step is larger, with $\gamma =0.066\pm 0.020$ mag. The 1D−5D γ difference for DES-SN is $0.026\pm 0.009$ mag. We show that this difference is due to a strong correlation between host galaxy stellar mass and the x1 component of the 5D distance-bias correction. Including an intrinsic correlation between the observed properties of SNe Ia, stretch and colour, and stellar mass in simulated SN Ia samples, we show that a 5D fit recovers γ with −9 mmag bias compared to a +2 mmag bias for a 1D fit. This difference can explain part of the discrepancy seen in the data. Improvements in modelling correlations between galaxy properties and SN is necessary to ensure unbiased precision estimates of the dark energy equation of state as we enter the era of LSST.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Predicting bulge to total luminosity ratio of galaxies using deep learning

ABSTRACT We present a deep learning model to predict the r-band bulge-to-total luminosity ratio (B/T) of nearby galaxies using their multiband JPEG images alone. Our Convolutional Neural Network (CNN) based regression model is trained on a large sample of galaxies with reliable decomposition into the bulge and disc components. The existing approaches to estimate the B/T ratio use galaxy light-profile modelling to find the best fit. This method is computationally expensive, prohibitively so for large samples of galaxies, and requires a significant amount of human intervention. Machine learning models have the potential to overcome these shortcomings. In our CNN model, for a test set of 20 000 galaxies, 85.7 per cent of the predicted B/T values have absolute error (AE) less than 0.1. We see further improvement to 87.5 per cent if, while testing, we only consider brighter galaxies (with r-band apparent magnitude <17) with no bright neighbours. Our model estimates the B/T ratio for the 20 000 test galaxies in less than a minute. This is a significant improvement in inference time from the conventional fitting pipelines, which manage around 2–3 estimates per minute. Thus, the proposed machine learning approach could potentially save a tremendous amount of time, effort, and computational resources while predicting B/T reliably, particularly in the era of next-generation sky surveys such as the Legacy Survey of Space and Time (LSST) and the Euclid sky survey which will produce extremely large samples of galaxies.

Grover, Harsh (ORCID:0000000321338142)↗
DOE OSTI2021

Optimizing a magnitude-limited spectroscopic training sample for photometric classification of supernovae

ABSTRACT In preparation for photometric classification of transients from the Legacy Survey of Space and Time (LSST) we run tests with different training data sets. Using estimates of the depth to which the 4-m Multi-Object Spectroscopic Telescope (4MOST) Time Domain Extragalactic Survey (TiDES) can classify transients, we simulate a magnitude-limited sample reaching rAB ≈ 22.5 mag. We run our simulations with the software snmachine, a photometric classification pipeline using machine learning. The machine-learning algorithms struggle to classify supernovae when the training sample is magnitude limited, in contrast to representative training samples. Classification performance noticeably improves when we combine the magnitude-limited training sample with a simulated realistic sample of faint high-redshift supernovae observed from larger spectroscopic facilities; the algorithms’ range of average area under receiver operator characteristic curve (AUC) scores over 10 runs increases from 0.547–0.628 to 0.946–0.969 and purity of the classified sample reaches 95 per cent in all runs for two of the four algorithms. By creating new, artificial light curves using the augmentation software avocado, we achieve a purity in our classified sample of 95 per cent in all 10 runs performed for all machine-learning algorithms considered. We also reach a highest average AUC score of 0.986 with the artificial neural network algorithm. Having ‘true’ faint supernovae to complement our magnitude-limited sample is a crucial requirement in optimization of a 4MOST spectroscopic sample. However, our results are a proof of concept that augmentation is also necessary to achieve the best classification results.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Mitigating baryonic effects with a theoretical error covariance

ABSTRACT One of the primary sources of uncertainties in modelling the cosmic-shear power spectrum on small scales is the effect of baryonic physics. Accurate cosmology for stage-IV surveys requires knowledge of the matter power spectrum deep in the non-linear regime at the per cent level. Therefore, it is important to develop reliable mitigation techniques to take into account baryonic uncertainties if information from small scales is to be considered in the cosmological analysis. In this work, we develop a new mitigation method for dealing with baryonic physics for the case of the shear angular power spectrum. The method is based on an augmented covariance matrix that incorporates baryonic uncertainties informed by hydrodynamical simulations. We use the results from 13 hydrodynamical simulations and the residual errors arising from a fit to a ΛCDM model using the extended halo model code HMCode to account for baryonic physics. These residual errors are used to model a so-called theoretical error covariance matrix that is added to the original covariance matrix. In order to assess the performance of the method, we use the 2D tomographic shear from four hydrodynamical simulations that have different extremes of baryonic parameters as mock data and run a likelihood analysis comparing the residual bias on Ωm and σ8 of our method and the HMCode for an LSST-like survey. We use different modelling of the theoretical error covariance matrix to test the robustness of the method. We show that it is possible to reduce the bias in the determination of the tested cosmological parameters at the price of a modest decrease in the precision.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Modelling type 1 quasar colours in the era of Rubin and Euclid

We construct a parametric SED model which is able to reproduce the average observed SDSS–UKIDSS–WISE quasar colours to within one-tenth of a magnitude across a wide range of redshift (0 < z < 5) and luminosity (–22 > M i > –29). This model is shown to provide accurate predictions for the colours of known quasars which are less luminous than those used to calibrate the model parameters, and also those at higher redshifts z > 5. Using a single parameter, the model encapsulates an up-to-date understanding of the intra-population variance in the rest-frame ultraviolet and optical emission lines of luminous quasars. At fixed redshift, there are systematic changes in the average quasar colours with apparent i-band magnitude, which we find to be well explained by the contribution from the host galaxy and our parametrization of the emission-line properties. By including redshift as an additional free parameter, the model could be used to provide photometric redshifts for individual objects. For the population as a whole we find that the average emission line and host-galaxy contributions can be well described by simple functions of luminosity which account for the observed changes in the average quasar colours across 18.1 < i AB < 21.5. Finally, we use these trends to provide predictions for quasar colours at the luminosities and redshifts which will be probed by the Rubin Observatory LSST and ESA-Euclid wide survey. The model code is applicable to a wide range of upcoming photometric and spectroscopic surveys, and is made publicly available.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Discovery of two bright high-redshift gravitationally lensed quasars revealed by Gaia

We present the discovery and preliminary characterisation of two high-redshift gravitationally lensed quasar systems in Gaia Data Release 2 (DR2). Candidates with multiple close-separation Gaia detections and quasar-like colours in WISE, Pan-STARRS, and DES are selected for follow-up spectroscopy with the New Technology Telescope. We confirm DES J215028.71-465251.3 as a $z$ = 4.130 ± 0.006 asymmetric, doubly imaged lensed quasar system and model the lensing mass distribution as a singular isothermal sphere. The system has an Einstein radius of 1.202 ± 0.005 arcsec and a predicted time delay of ~122.0 d between the quasar images, assuming a lensing galaxy redshift of $z$ = 0.5, making this a priority system for future optical monitoring. We confirm PS J042913.17+142840.9 as a $z$ = 3.866 ± 0.003 four-image quasar system in a cusp configuration, lensed by two foreground galaxies. The system is well modelled using a singular isothermal ellipsoid for the primary lens and a singular isothermal sphere for the secondary lens with Einstein radii 0.704 ± 0.006 and 0.241 ± 0.030 arcsec, respectively. A maximum predicted time delay of 9.6 d is calculated, assuming lensing galaxy redshifts of $z$ = 1.0. Furthermore, PS J042913.17+142840.9 exhibits a large flux ratio anomaly, up to a factor of 2.66 ± 0.37 in i band, that varies across optical and near-infrared wavelengths. We discuss LSST and its implications for future high-redshift lens searches and outline an extension to the search using supervised machine learning techniques.

Astronomy & Astrophysics↗
DOE OSTI2021

Probing gravity with the DES-CMASS sample and BOSS spectroscopy

ABSTRACT The DES-CMASS sample (DMASS) is designed to optimally combine the weak lensing measurements from the Dark Energy Survey (DES) and redshift-space distortions (RSD) probed by the CMASS galaxy sample from the Baryonic Oscillation Spectroscopic Survey. In this paper, we demonstrate the feasibility of adopting DMASS as the equivalent of CMASS for a joint analysis of DES and BOSS in the framework of modified gravity. We utilize the angular clustering of the DMASS galaxies, cosmic shear of the DES metacalibration sources, and cross-correlation of the two as data vectors. By jointly fitting the combination of the data with the RSD measurements from the CMASS sample and Planck data, we obtain the constraints on modified gravity parameters $\mu _0=-0.37^{+0.47}_{-0.45}$ and $\Sigma _0=0.078^{+0.078}_{-0.082}$. Our constraints of modified gravity with DMASS are tighter than those with the DES Year 1 redMaGiC sample with the same external data sets by 29 per cent for μ0 and 21 per cent for Σ0, and comparable to the published results of the DES Year 1 modified gravity analysis despite this work using fewer external data sets. This improvement is mainly because the galaxy bias parameter is shared and more tightly constrained by both CMASS and DMASS, effectively breaking the degeneracy between the galaxy bias and other cosmological parameters. Such an approach to optimally combine photometric and spectroscopic surveys using a photometric sample equivalent to a spectroscopic sample can be applied to combining future surveys having a limited overlap such as DESI and LSST.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Map-based cosmology inference with lognormal cosmic shear maps

ABSTRACT Most cosmic shear analyses to date have relied on summary statistics (e.g. ξ+ and ξ−). These types of analyses are necessarily suboptimal, as the use of summary statistics is lossy. In this paper, we forward-model the convergence field of the Universe as a lognormal random field conditioned on the observed shear data. This new map-based inference framework enables us to recover the joint posterior of the cosmological parameters and the convergence field of the Universe. Our analysis properly accounts for the covariance in the mass maps across tomographic bins, which significantly improves the fidelity of the maps relative to single-bin reconstructions. We verify that applying our inference pipeline to Gaussian random fields recovers posteriors that are in excellent agreement with their analytical counterparts. At the resolution of our maps – and to the extent that the convergence field can be described by the lognormal model – our map posteriors allow us to reconstruct all summary statistics (including non-Gaussian statistics). We forecast that a map-based inference analysis of LSST-Y10 data can improve cosmological constraints in the σ8–Ωm plane by $\approx\!{30}{{\ \rm per\ cent}}$ relative to the currently standard cosmic shear analysis. This improvement happens almost entirely along the $S_8=\sigma _8\Omega _{\rm m}^{1/2}$ directions, meaning map-based inference fails to significantly improve constraints on S8.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Photometric calibration in u -band using blue halo stars

ABSTRACT We develop a method to calibrate u-band photometry based on the observed colour of blue Galactic halo stars. The Galactic halo stars belong to an old stellar population of the Milky Way and have relatively low metallicity. The ‘blue tip’ of the halo population – the main sequence turn-off (MSTO) stars – is known to have a relatively uniform intrinsic edge u-g colour with only slow spatial variation. In SDSS data, the observed variation is correlated with Galactic Latitude, which we attribute to contamination by higher metallicity disc stars and fit with an empirical curve. This curve can then be used to calibrate u-band imaging if g-band imaging of matching depth is available. Our approach can be applied to single-field observations at |b| > 30°, and removes the need for standard star observations or overlap with calibrated u-band imaging. We include in our method the calibration of g-band data with ATLAS-Refcat2. We test our approach on stars in KiDS DR 4, ATLAS DR 4, and DECam imaging from the NOIRLab Source Catalog (NSC DR2), and compare our calibration with SDSS. For this process, we use synthetic magnitudes to derive the colour equations between these data sets, in order to improve zero-point accuracy. We find an improvement for all data sets, reaching a zero-point precision of 0.016 mag for KiDS (compared to the original 0.033 mag), 0.020 mag for ATLAS (originally 0.027 mag), and 0.016 mag for DECam (originally 0.041 mag). Thus, this method alone reaches the goal of 0.02 mag photometric precision in u-band for the Rubin Observatory’s Legacy Survey of Space and Time (LSST).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Deep learning methods for obtaining photometric redshift estimations from images

ABSTRACT Knowing the redshift of galaxies is one of the first requirements of many cosmological experiments, and as it is impossible to perform spectroscopy for every galaxy being observed, photometric redshift (photo-z) estimations are still of particular interest. Here, we investigate different deep learning methods for obtaining photo-z estimates directly from images, comparing these with ‘traditional’ machine learning algorithms which make use of magnitudes retrieved through photometry. As well as testing a convolutional neural network (CNN) and inception-module CNN, we introduce a novel mixed-input model that allows for both images and magnitude data to be used in the same model as a way of further improving the estimated redshifts. We also perform benchmarking as a way of demonstrating the performance and scalability of the different algorithms. The data used in the study comes entirely from the Sloan Digital Sky Survey (SDSS) from which 1 million galaxies were used, each having 5-filtre (ugriz) images with complete photometry and a spectroscopic redshift which was taken as the ground truth. The mixed-input inception CNN achieved a mean squared error (MSE) =0.009, which was a significant improvement ($30{{\ \rm per\ cent}}$) over the traditional random forest (RF), and the model performed even better at lower redshifts achieving a MSE = 0.0007 (a $50{{\ \rm per\ cent}}$ improvement over the RF) in the range of z < 0.3. This method could be hugely beneficial to upcoming surveys, such as Euclid and the Vera C. Rubin Observatory’s Legacy Survey of Space and Time (LSST), which will require vast numbers of photo-z estimates produced as quickly and accurately as possible.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

OzDES Reverberation Mapping Program: Hβ lags from the 6-yr survey

ABSTRACT Reverberation mapping measurements have been used to constrain the relationship between the size of the broad-line region and luminosity of active galactic nuclei (AGN). This R–L relation is used to estimate single-epoch virial black hole masses, and has been proposed to use to standardize AGN to determine cosmological distances. We present reverberation measurements made with Hβ from the 6-yr Australian Dark Energy Survey (OzDES) Reverberation Mapping Program. We successfully recover reverberation lags for eight AGN at 0.12 < z < 0.71, probing higher redshifts than the bulk of Hβ measurements made to date. Our fit to the R–L relation has a slope of α = 0.41 ± 0.03 and an intrinsic scatter of σ = 0.23 ± 0.02 dex. The results from our multi-object spectroscopic survey are consistent with previous measurements made by dedicated source-by-source campaigns, and with the observed dependence on accretion rate. Future surveys, including LSST, TiDES, and SDSS-V, which will be revisiting some of our observed fields, will be able to build on the results of our first-generation multi-object reverberation mapping survey.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Sky subtraction in an era of low surface brightness astronomy

ABSTRACT The Vera C. Rubin Observatory Wide-Fast Deep sky survey will reach unprecedented surface brightness depths over tens of thousands of square degrees. Surface brightness photometry has traditionally been a challenge. Current algorithms which combine object detection with sky estimation systematically oversubtract the sky, biasing surface brightness measurements at the faint end and destroying or severely compromising low surface brightness light. While it has recently been shown that properly accounting for undetected faint galaxies and the wings of brighter objects can in principle recover a more accurate sky estimate, this has not yet been demonstrated in practice. Obtaining a consistent spatially smooth underlying sky estimate is particularly challenging in the presence of representative distributions of bright and faint objects. In this paper, we use simulations of crowded and uncrowded fields designed to mimic Hyper Suprime-Cam data to perform a series of tests on the accuracy of the recovered sky. Dependence on field density, galaxy type, and limiting flux for detection are all considered. Several photometry packages are utilized: source extractor, gnuastro, and the LSST science pipelines. Each is configured in various modes, and their performance at extreme low surface brightness analysed. We find that the combination of the source extractor software package with novel source model masking techniques consistently produce extremely faint output sky estimates, by up to an order of magnitude, as well as returning high fidelity output science catalogues.

Astronomy & Astrophysics↗
DOE OSTI2023

Nucleosynthesis in outflows of compact objects and detection prospects of associated kilonovae

ABSTRACT We perform a comparative analysis of nucleosynthesis yields from binary neutron star (BNS) mergers, black hole-neutron star (BHNS) mergers, and core-collapse supernovae (CCSNe) with the goal of determining which are the most dominant sources of r-process enrichment observed in stars. We find that BNS and BHNS binaries may eject similar mass distributions of robust r-process nuclei post-merger (up to third peak and actinides, A ∼ 200−240), after accounting for the volumetric event rates. Magnetorotational (MR) CCSNe likely undergo a weak r-process (up to A ∼ 140) and contribute to the production of light element primary process (LEPP) nuclei, whereas typical thermal, neutrino-driven CCSNe only synthesize up to first r-process peak nuclei (A ∼ 80−90). We also find that the upper limit to the rate of MR CCSNe is $\lesssim 1~{{\ \rm per\ cent}}$ the rate of typical thermal CCSNe; if the rate was higher, then weak r-process nuclei would be overproduced. Although the largest uncertainty is from the volumetric event rate, the prospects are encouraging for confirming these rates in the next few years with upcoming surveys. Using a simple model to estimate the resulting kilonova light curve from mergers and our set of fiducial merger parameters, we predict that ∼7 BNS and ∼2 BHNS events will be detectable per year by the Vera C. Rubin Observatory (LSST), with prior gravitational wave (GW) triggers.

Astronomy & Astrophysics↗
DOE OSTI2023

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↗
DOE OSTI2023

The simulated catalogue of optical transients and correlated hosts (SCOTCH)

ABSTRACT As we observe a rapidly growing number of astrophysical transients, we learn more about the diverse host galaxy environments in which they occur. Host galaxy information can be used to purify samples of cosmological Type Ia supernovae, uncover the progenitor systems of individual classes, and facilitate low-latency follow-up of rare and peculiar explosions. In this work, we develop a novel data-driven methodology to simulate the time-domain sky that includes detailed modelling of the probability density function for multiple transient classes conditioned on host galaxy magnitudes, colours, star formation rates, and masses. We have designed these simulations to optimize photometric classification and analysis in upcoming large synoptic surveys. We integrate host galaxy information into the snana simulation framework to construct the simulated catalogue of optical transients and correlated hosts (SCOTCH, a publicly available catalogue of 5-million idealized transient light curves in LSST passbands and their host galaxy properties over the redshift range 0 < z < 3. This catalogue includes supernovae, tidal disruption events, kilonovae, and active galactic nuclei. Each light curve consists of true top-of-the-galaxy magnitudes sampled with high (≲2 d) cadence. In conjunction with SCOTCH, we also release an associated set of tutorials and transient-specific libraries to enable simulations of arbitrary space- and ground-based surveys. Our methodology is being used to test critical science infrastructure in advance of surveys by the Vera C. Rubin Observatory and the Nancy G. Roman Space Telescope.

Astronomy & Astrophysics↗