-
The Blending ToolKit: A simulation framework for evaluation of galaxy detection and deblending
Authors:
Ismael Mendoza,
Andrii Torchylo,
Thomas Sainrat,
Axel Guinot,
Alexandre Boucaud,
Maxime Paillasa,
Camille Avestruz,
Prakruth Adari,
Eric Aubourg,
Biswajit Biswas,
James Buchanan,
Patricia Burchat,
Cyrille Doux,
Remy Joseph,
Sowmya Kamath,
Alex I. Malz,
Grant Merz,
Hironao Miyatake,
Cécile Roucelle,
Tianqing Zhang,
the LSST Dark Energy Science Collaboration
Abstract:
We present an open source Python library for simulating overlapping (i.e., blended) images of galaxies and performing self-consistent comparisons of detection and deblending algorithms based on a suite of metrics. The package, named Blending Toolkit (BTK), serves as a modular, flexible, easy-to-install, and simple-to-use interface for exploring and analyzing systematic effects related to blended g…
▽ More
We present an open source Python library for simulating overlapping (i.e., blended) images of galaxies and performing self-consistent comparisons of detection and deblending algorithms based on a suite of metrics. The package, named Blending Toolkit (BTK), serves as a modular, flexible, easy-to-install, and simple-to-use interface for exploring and analyzing systematic effects related to blended galaxies in cosmological surveys such as the Vera Rubin Observatory Legacy Survey of Space and Time (LSST). BTK has three main components: (1) a set of modules that perform fast image simulations of blended galaxies, using the open source image simulation package GalSim; (2) a module that standardizes the inputs and outputs of existing deblending algorithms; (3) a library of deblending metrics commonly defined in the galaxy deblending literature. In combination, these modules allow researchers to explore the impacts of galaxy blending in cosmological surveys. Additionally, BTK provides researchers who are developing a new deblending algorithm a framework to evaluate algorithm performance and make principled comparisons with existing deblenders. BTK includes a suite of tutorials and comprehensive documentation. The source code is publicly available on GitHub at https://github.com/LSSTDESC/BlendingToolKit.
△ Less
Submitted 10 September, 2024;
originally announced September 2024.
-
YOLO-CL cluster detection in the Rubin/LSST DC2 simulation
Authors:
Kirill Grishin,
Simona Mei,
Stephane Ilic,
Michel Aguena,
Dominique Boutigny,
Marie Paturel,
the LSST Dark Energy Science Collaboration
Abstract:
LSST will provide galaxy cluster catalogs up to z$\sim$1 that can be used to constrain cosmological models once their selection function is well-understood. We have applied the deep convolutional network YOLO for CLuster detection (YOLO-CL) to LSST simulations from the Dark Energy Science Collaboration Data Challenge 2 (DC2), and characterized the LSST YOLO-CL cluster selection function. We have t…
▽ More
LSST will provide galaxy cluster catalogs up to z$\sim$1 that can be used to constrain cosmological models once their selection function is well-understood. We have applied the deep convolutional network YOLO for CLuster detection (YOLO-CL) to LSST simulations from the Dark Energy Science Collaboration Data Challenge 2 (DC2), and characterized the LSST YOLO-CL cluster selection function. We have trained and validated the network on images from a hybrid sample of (1) clusters observed in the Sloan Digital Sky Survey and detected with the red-sequence Matched-filter Probabilistic Percolation, and (2) simulated DC2 dark matter haloes with masses $M_{200c} > 10^{14} M_{\odot}$. We quantify the completeness and purity of the YOLO-CL cluster catalog with respect to DC2 haloes with $M_{200c} > 10^{14} M_{\odot}$. The YOLO-CL cluster catalog is 100% and 94% complete for halo mass $M_{200c} > 10^{14.6} M_{\odot}$ at $0.2<z<0.8$, and $M_{200c} > 10^{14} M_{\odot}$ and redshift $z \lesssim 1$, respectively, with only 6% false positive detections. All the false positive detections are dark matter haloes with $ 10^{13.4} M_{\odot} \lesssim M_{200c} \lesssim 10^{14} M_{\odot}$. The YOLO-CL selection function is almost flat with respect to the halo mass at $0.2 \lesssim z \lesssim 0.9$. The overall performance of YOLO-CL is comparable or better than other cluster detection methods used for current and future optical and infrared surveys. YOLO-CL shows better completeness for low mass clusters when compared to current detections in surveys using the Sunyaev Zel'dovich effect, and detects clusters at higher redshifts than X-ray-based catalogs. The strong advantage of YOLO-CL over traditional galaxy cluster detection techniques is that it works directly on images and does not require photometric and photometric redshift catalogs, nor does it need to mask stellar sources and artifacts.
△ Less
Submitted 5 September, 2024;
originally announced September 2024.
-
Impact of survey spatial variability on galaxy redshift distributions and the cosmological $3\times2$-point statistics for the Rubin Legacy Survey of Space and Time (LSST)
Authors:
Qianjun Hang,
Benjamin Joachimi,
Eric Charles,
John Franklin Crenshaw,
Patricia Larsen,
Alex I. Malz,
Sam Schmidt,
Ziang Yan,
Tianqing Zhang,
the LSST Dark Energy Science Collaboration
Abstract:
We investigate the impact of spatial survey non-uniformity on the galaxy redshift distributions for forthcoming data releases of the Rubin Observatory Legacy Survey of Space and Time (LSST). Specifically, we construct a mock photometry dataset degraded by the Rubin OpSim observing conditions, and estimate photometric redshifts of the sample using a template-fitting photo-$z$ estimator, BPZ, and a…
▽ More
We investigate the impact of spatial survey non-uniformity on the galaxy redshift distributions for forthcoming data releases of the Rubin Observatory Legacy Survey of Space and Time (LSST). Specifically, we construct a mock photometry dataset degraded by the Rubin OpSim observing conditions, and estimate photometric redshifts of the sample using a template-fitting photo-$z$ estimator, BPZ, and a machine learning method, FlexZBoost. We select the Gold sample, defined as $i<25.3$ for 10 year LSST data, with an adjusted magnitude cut for each year and divide it into five tomographic redshift bins for the weak lensing lens and source samples. We quantify the change in the number of objects, mean redshift, and width of each tomographic bin as a function of the coadd $i$-band depth for 1-year (Y1), 3-year (Y3), and 5-year (Y5) data. In particular, Y3 and Y5 have large non-uniformity due to the rolling cadence of LSST, hence provide a worst-case scenario of the impact from non-uniformity. We find that these quantities typically increase with depth, and the variation can be $10-40\%$ at extreme depth values. Based on these results and using Y3 as an example, we propagate the variable depth effect to the weak lensing $3\times2$pt data vector in harmonic space. We find that galaxy clustering is most susceptible to variable depth, causing significant deviations at large scales if not corrected for, due to the depth-dependent number density variations. For galaxy-shear and shear-shear power spectra, we find little impact given the expected LSST Y3 noise.
△ Less
Submitted 4 September, 2024;
originally announced September 2024.
-
Modeling the 3-point correlation function of projected scalar fields on the sphere
Authors:
Abraham Arvizu,
Alejandro Aviles,
Juan Carlos Hidalgo,
Eladio Moreno,
Gustavo Niz,
Mario A. Rodriguez-Meza,
Sofía Samario,
The LSST Dark Energy Science Collaboration
Abstract:
One of the main obstacles for the signal extraction of the three point correlation function using photometric surveys, such as the Rubin Observatory Legacy Survey of Space and Time (LSST), will be the prohibitive computation time required for dealing with a vast quantity of sources. Brute force algorithms, which naively scales as $\mathcal{O}(N^3)$ with the number of objects, can be further improv…
▽ More
One of the main obstacles for the signal extraction of the three point correlation function using photometric surveys, such as the Rubin Observatory Legacy Survey of Space and Time (LSST), will be the prohibitive computation time required for dealing with a vast quantity of sources. Brute force algorithms, which naively scales as $\mathcal{O}(N^3)$ with the number of objects, can be further improved with tree methods but not enough to deal with large scale correlations of Rubin's data. However, a harmonic basis decomposition of these higher order statistics reduces the time dramatically, to scale as a two-point correlation function with the number of objects, so that the signal can be extracted in a reasonable amount of time. In this work, we aim to develop the framework to use these expansions within the Limber approximation for scalar (or spin-0) fields, such as galaxy counts, weak lensing convergence or aperture masses. We develop an estimator to extract the signal from catalogs and different phenomenological and theoretical models for its description. The latter includes halo model and standard perturbation theory, to which we add a simple effective field theory prescription based on the short range of non-locality of cosmic fields, significantly improving the agreement with simulated data. In parallel to the modeling of the signal, we develop a code that can efficiently calculate three points correlations of more than 200 million data points (a full sky simulation with Nside=4096) in $\sim$40 minutes on a single high-performance computing node, enabling a feasible analysis for the upcoming LSST data.
△ Less
Submitted 29 August, 2024;
originally announced August 2024.
-
MADNESS Deblender: Maximum A posteriori with Deep NEural networks for Source Separation
Authors:
Biswajit Biswas,
Eric Aubourg,
Alexandre Boucaud,
Axel Guinot,
Junpeng Lao,
Cécile Roucelle,
the LSST Dark Energy Science Collaboration
Abstract:
Due to the unprecedented depth of the upcoming ground-based Legacy Survey of Space and Time (LSST) at the Vera C. Rubin Observatory, approximately two-thirds of the galaxies are likely to be affected by blending - the overlap of physically separated galaxies in images. Thus, extracting reliable shapes and photometry from individual objects will be limited by our ability to correct blending and con…
▽ More
Due to the unprecedented depth of the upcoming ground-based Legacy Survey of Space and Time (LSST) at the Vera C. Rubin Observatory, approximately two-thirds of the galaxies are likely to be affected by blending - the overlap of physically separated galaxies in images. Thus, extracting reliable shapes and photometry from individual objects will be limited by our ability to correct blending and control any residual systematic effect. Deblending algorithms tackle this issue by reconstructing the isolated components from a blended scene, but the most commonly used algorithms often fail to model complex realistic galaxy morphologies.
As part of an effort to address this major challenge, we present MADNESS, which takes a data-driven approach and combines pixel-level multi-band information to learn complex priors for obtaining the maximum a posteriori solution of deblending. MADNESS is based on deep neural network architectures such as variational auto-encoders and normalizing flows. The variational auto-encoder reduces the high-dimensional pixel space into a lower-dimensional space, while the normalizing flow models a data-driven prior in this latent space.
Using a simulated test dataset with galaxy models for a 10-year LSST survey and a galaxy density ranging from 48 to 80 galaxies per arcmin2 we characterize the aperture-photometry g-r color, structural similarity index, and pixel cosine similarity of the galaxies reconstructed by MADNESS. We compare our results against state-of-the-art deblenders including scarlet. With the r-band of LSST as an example, we show that MADNESS performs better than in all the metrics. For instance, the average absolute value of relative flux residual in the r-band for MADNESS is approximately 29% lower than that of scarlet. The code is publicly available on GitHub.
△ Less
Submitted 27 August, 2024;
originally announced August 2024.
-
Analytical Weak Lensing Shear Inference for Precision Cosmology
Authors:
Xiangchong Li,
Rachel Mandelbaum,
The LSST Dark Energy Science Collaboration
Abstract:
Noise bias is a significant source of systematic error in weak gravitational lensing measurements that must be corrected to satisfy the stringent standards of modern imaging surveys in the era of precision cosmology. This paper reviews the analytical noise bias correction method and provides analytical derivations demonstrating that we can recover shear to its second order using the 'renoising' no…
▽ More
Noise bias is a significant source of systematic error in weak gravitational lensing measurements that must be corrected to satisfy the stringent standards of modern imaging surveys in the era of precision cosmology. This paper reviews the analytical noise bias correction method and provides analytical derivations demonstrating that we can recover shear to its second order using the 'renoising' noise bias correction approach introduced by Metacalibration. We implement this analytical noise bias correction within the AnaCal shear estimation framework and propose several enhancements to the noise bias correction algorithm. We evaluate the improved AnaCal using simulations designed to replicate Rubin LSST imaging data. These simulations feature semi-realistic galaxies and stars, complete with representative distributions of magnitudes and Galactic spatial density. We conduct tests under various observational challenges, including cosmic rays, defective CCD columns, bright star saturation, bleed trails, and spatially variable point spread functions. Our results indicate a multiplicative bias in weak lensing shear recovery of less than a few tenths of a percent, meeting LSST DESC requirements without requiring calibration from external image simulations. Additionally, our algorithm achieves rapid processing, handling one galaxy in less than a millisecond.
△ Less
Submitted 12 August, 2024;
originally announced August 2024.
-
The little coadd that could: Estimating shear from coadded images
Authors:
Robert Armstrong,
Erin Sheldon,
Eric Huff,
Jim Bosch,
Eli Rykoff,
Rachel Mandelbaum,
Arun Kannawadi,
Peter Melchior,
Robert Lupton,
Matthew R. Becker,
Yusra Al-Sayyed,
The LSST Dark Energy Science Collaboration
Abstract:
Upcoming wide field surveys will have many overlapping epochs of the same region of sky. The conventional wisdom is that in order to reduce the errors sufficiently for systematics-limited measurements, like weak lensing, we must do simultaneous fitting of all the epochs. Using current algorithms this will require a significant amount of computing time and effort. In this paper, we revisit the pote…
▽ More
Upcoming wide field surveys will have many overlapping epochs of the same region of sky. The conventional wisdom is that in order to reduce the errors sufficiently for systematics-limited measurements, like weak lensing, we must do simultaneous fitting of all the epochs. Using current algorithms this will require a significant amount of computing time and effort. In this paper, we revisit the potential of using coadds for shear measurements. We show on a set of image simulations that the multiplicative shear bias can be constrained below the 0.1% level on coadds, which is sufficient for future lensing surveys. We see no significant differences between simultaneous fitting and coadded approaches for two independent shear codes: Metacalibration and BFD. One caveat of our approach is the assumption of a principled coadd, i.e. the PSF is mathematically well-defined for all the input images. This requires us to reject CCD images that do not fully cover the coadd region. We estimate that the number of epochs that must be rejected for a survey like LSST is on the order of 20%, resulting in a small loss in depth of less than 0.1 magnitudes. We also put forward a cell-based coaddition scheme that meets the above requirements for unbiased weak lensing shear estimation in the context of LSST.
△ Less
Submitted 1 July, 2024;
originally announced July 2024.
-
Accurate measurement of telescope filter bandpasses with a Collimated Beam Projector and impact on cosmological parameters
Authors:
Jérémy Neveu,
Dylan Kuhn,
Thierry Souverin,
LEMAITRE collaboration
Abstract:
The measurement of magnitudes with different filters in photometric surveys gives access to cosmological distances and parameters. However, for current and future large surveys like the ZTF, DES, HSC or LSST, the photometric calibration uncertainties are almost comparable to statistical uncertainties in the error budget of type Ia cosmology analysis, which limits our ability to use type Ia superno…
▽ More
The measurement of magnitudes with different filters in photometric surveys gives access to cosmological distances and parameters. However, for current and future large surveys like the ZTF, DES, HSC or LSST, the photometric calibration uncertainties are almost comparable to statistical uncertainties in the error budget of type Ia cosmology analysis, which limits our ability to use type Ia supernovae for precision cosmology. The knowledge of the bandpasses of the survey filters at the per-mill level can help reach the sub-percent precision for magnitudes. We show how a misknowledge of the bandpasses central wavelengths or of the presence of out-of-band leakages leads to biased cosmological measurements. Then, we present how to measure the filter throughputs at the required precision with a Collimated Beam Projector.
△ Less
Submitted 1 July, 2024;
originally announced July 2024.
-
On the importance of Earth's atmosphere for SNIa precision cosmology
Authors:
Jérémy Neveu,
LEMAITRE collaboration,
AuxTel collaboration
Abstract:
The measurement of colours in photometric surveys is a path to get access to cosmological distances. But for future large surveys like the Large Survey of Space and Time undertaken by the Vera Rubin Observatory in Chile, the large statistical power of the promised catalogues will make the photometric calibration uncertainties dominant in the error budget and will limit our ability to use it for pr…
▽ More
The measurement of colours in photometric surveys is a path to get access to cosmological distances. But for future large surveys like the Large Survey of Space and Time undertaken by the Vera Rubin Observatory in Chile, the large statistical power of the promised catalogues will make the photometric calibration uncertainties dominant in the error budget and will limit our ability to use it for precision cosmology. The knowledge of the in situ atmospheric transmission for each exposure, each season and on average for the entire survey can help reaching the subpercent precision for magnitudes. I will show the impact of precipitable water vapour during a cosmological survey on supernova cosmology. Then I will present how AuxTel at the Vera Rubin Observatory have the capability to measure the on-site atmospheric transmission in real time to improve LSST photometric calibration for precision cosmology.
△ Less
Submitted 1 July, 2024;
originally announced July 2024.
-
MNTES: Modeling Nonlinearity of TES detectors for Enhanced Cosmic Microwave Background measurements with LiteBIRD
Authors:
Tijmen de Haan,
LiteBIRD Collaboration
Abstract:
Traditional methods of converting electronic readout counts to optical power incident on Transition Edge Sensors (TES) for Cosmic Microwave Background (CMB) observations involve a linear approximation. For the upcoming LiteBIRD CMB satellite, strict nonlinearity requirements must be met to prevent contamination of the science band at 4f by the 2f signal, which arises from differential transmission…
▽ More
Traditional methods of converting electronic readout counts to optical power incident on Transition Edge Sensors (TES) for Cosmic Microwave Background (CMB) observations involve a linear approximation. For the upcoming LiteBIRD CMB satellite, strict nonlinearity requirements must be met to prevent contamination of the science band at 4f by the 2f signal, which arises from differential transmission or emissivity related to the half-wave plate's rotation rate fHWP. These constraints cannot be met using hardware solutions alone and therefore require a form of nonlinearity correction. We present MNTES, a novel physics-based, nonlinear calibration technique. This method leverages our physical understanding of the TES power balance equation, accounts for imperfect voltage bias by casting the bias network as its Thévenin equivalent, and can incorporate external information such as time-varying magnetic fields and focal plane temperature variations. The detector-specific parameters of MNTES will be measured during the ground calibration campaign prior to the LiteBIRD launch, yielding conversion functions that can take raw time-ordered data and output the reconstructed incident optical power. MNTES will allow us to achieve LiteBIRD's goal of measuring the primordial tensor fluctuation spectrum to δr < 0.001.
△ Less
Submitted 27 June, 2024;
originally announced June 2024.
-
Large Interferometer For Exoplanets (LIFE): XIII. The Value of Combining Thermal Emission and Reflected Light for the Characterization of Earth Twins
Authors:
E. Alei,
S. P. Quanz,
B. S. Konrad,
E. O. Garvin,
V. Kofman,
A. Mandell,
D. Angerhausen,
P. Mollière,
M. R. Meyer,
T. Robinson,
S. Rugheimer,
the LIFE Collaboration
Abstract:
Following the recommendations to NASA and ESA, the search for life on exoplanets will be a priority in the next decades. Two direct imaging space mission concepts are being developed: the Habitable Worlds Observatory (HWO) and the Large Interferometer for Exoplanets (LIFE). HWO focuses on reflected light spectra in the ultraviolet/visible/near-infrared (UV/VIS/NIR), while LIFE captures the mid-inf…
▽ More
Following the recommendations to NASA and ESA, the search for life on exoplanets will be a priority in the next decades. Two direct imaging space mission concepts are being developed: the Habitable Worlds Observatory (HWO) and the Large Interferometer for Exoplanets (LIFE). HWO focuses on reflected light spectra in the ultraviolet/visible/near-infrared (UV/VIS/NIR), while LIFE captures the mid-infrared (MIR) emission of temperate exoplanets. We assess the potential of HWO and LIFE in characterizing a cloud-free Earth twin orbiting a Sun-like star at 10 pc, both separately and synergistically, aiming to quantify the increase in information from joint atmospheric retrievals on a habitable planet. We perform Bayesian retrievals on simulated data from an HWO-like and a LIFE-like mission separately, then jointly, considering the baseline spectral resolutions currently assumed for these concepts and using two increasingly complex noise simulations. HWO would constrain H$_2$O, O$_2$, and O$_3$, in the atmosphere, with ~ 100 K uncertainty on the temperature profile. LIFE would constrain CO$_2$, H$_2$O, O$_3$ and provide constraints on the thermal atmospheric structure and surface temperature (~ 10 K uncertainty). Both missions would provide an upper limit on CH$_4$. Joint retrievals on HWO and LIFE data would accurately define the atmospheric thermal profile and planetary parameters, decisively constrain CO$_2$, H$_2$O, O$_2$, and O$_3$, and weakly constrain CO and CH$_4$. The detection significance is greater or equal to single-instrument retrievals. Both missions provide specific information to characterize a terrestrial habitable exoplanet, but the scientific yield is maximized with synergistic UV/VIS/NIR+MIR observations. Using HWO and LIFE together will provide stronger constraints on biosignatures and life indicators, potentially transforming the search for life in the universe.
△ Less
Submitted 18 June, 2024;
originally announced June 2024.
-
Joint Modelling of Astrophysical Systematics for Cosmology with LSST
Authors:
Nikolina Šarčević,
C. Danielle Leonard,
Markus M. Rau,
LSST Dark Energy Science Collaboration
Abstract:
We present a novel framework for jointly modelling the weak lensing source galaxy redshift distribution and the intrinsic alignment of galaxies via a shared luminosity function (LF). Considering this framework within the context of a Rubin Observatory's Legacy Survey of Space and Time (LSST) Year 1 and Year 10 cosmic shear analysis, we first demonstrate the substantial impact of the LF on both sou…
▽ More
We present a novel framework for jointly modelling the weak lensing source galaxy redshift distribution and the intrinsic alignment of galaxies via a shared luminosity function (LF). Considering this framework within the context of a Rubin Observatory's Legacy Survey of Space and Time (LSST) Year 1 and Year 10 cosmic shear analysis, we first demonstrate the substantial impact of the LF on both source galaxy redshift distributions and the intrinsic alignment contamination. We establish how the individual parameters of a Schechter LF model impact the redshift distribution of a magnitude-limited sample, and we demonstrate the effect of marginalising over the LF parameters as incorporated in the intrinsic alignment modelling of a standard cosmic shear analysis set-up. We forecast the impact of our joint modelling framework on cosmological parameter constraints. Our preliminary results are promising, indicating that this framework can yield cosmological constraints consistent with those expected from standard analyses, enhanced by the flexibility of not fixing LF parameters. We plan to further validate these findings with comprehensive Markov chain Monte Carlo simulations to robustly quantify bias avoidance and underscore the framework's efficacy. Taking advantage of our forecasting results and the parameter degeneracies, we identify the specific impact of the shape of the LF of source galaxies on the cosmic shear data vector. We also discuss the potential of this method in providing a way to model generic selection functions in redshift distribution estimation, as well as its possibilities for extension to a 3x2pt analysis, particularly with respect to incorporating galaxy bias in this luminosity-function-based framework. Although we consider the context of LSST cosmic shear in this work, the proposed joint modelling framework is generically applicable to weak lensing surveys.
△ Less
Submitted 12 June, 2024; v1 submitted 5 June, 2024;
originally announced June 2024.
-
The LiteBIRD mission to explore cosmic inflation
Authors:
T. Ghigna,
A. Adler,
K. Aizawa,
H. Akamatsu,
R. Akizawa,
E. Allys,
A. Anand,
J. Aumont,
J. Austermann,
S. Azzoni,
C. Baccigalupi,
M. Ballardini,
A. J. Banday,
R. B. Barreiro,
N. Bartolo,
S. Basak,
A. Basyrov,
S. Beckman,
M. Bersanelli,
M. Bortolami,
F. Bouchet,
T. Brinckmann,
P. Campeti,
E. Carinos,
A. Carones
, et al. (134 additional authors not shown)
Abstract:
LiteBIRD, the next-generation cosmic microwave background (CMB) experiment, aims for a launch in Japan's fiscal year 2032, marking a major advancement in the exploration of primordial cosmology and fundamental physics. Orbiting the Sun-Earth Lagrangian point L2, this JAXA-led strategic L-class mission will conduct a comprehensive mapping of the CMB polarization across the entire sky. During its 3-…
▽ More
LiteBIRD, the next-generation cosmic microwave background (CMB) experiment, aims for a launch in Japan's fiscal year 2032, marking a major advancement in the exploration of primordial cosmology and fundamental physics. Orbiting the Sun-Earth Lagrangian point L2, this JAXA-led strategic L-class mission will conduct a comprehensive mapping of the CMB polarization across the entire sky. During its 3-year mission, LiteBIRD will employ three telescopes within 15 unique frequency bands (ranging from 34 through 448 GHz), targeting a sensitivity of 2.2\,$μ$K-arcmin and a resolution of 0.5$^\circ$ at 100\,GHz. Its primary goal is to measure the tensor-to-scalar ratio $r$ with an uncertainty $δr = 0.001$, including systematic errors and margin. If $r \geq 0.01$, LiteBIRD expects to achieve a $>5σ$ detection in the $\ell=$2-10 and $\ell=$11-200 ranges separately, providing crucial insight into the early Universe. We describe LiteBIRD's scientific objectives, the application of systems engineering to mission requirements, the anticipated scientific impact, and the operations and scanning strategies vital to minimizing systematic effects. We will also highlight LiteBIRD's synergies with concurrent CMB projects.
△ Less
Submitted 4 June, 2024;
originally announced June 2024.
-
Probabilistic Forward Modeling of Galaxy Catalogs with Normalizing Flows
Authors:
John Franklin Crenshaw,
J. Bryce Kalmbach,
Alexander Gagliano,
Ziang Yan,
Andrew J. Connolly,
Alex I. Malz,
Samuel J. Schmidt,
The LSST Dark Energy Science Collaboration
Abstract:
Evaluating the accuracy and calibration of the redshift posteriors produced by photometric redshift (photo-z) estimators is vital for enabling precision cosmology and extragalactic astrophysics with modern wide-field photometric surveys. Evaluating photo-z posteriors on a per-galaxy basis is difficult, however, as real galaxies have a true redshift but not a true redshift posterior. We introduce P…
▽ More
Evaluating the accuracy and calibration of the redshift posteriors produced by photometric redshift (photo-z) estimators is vital for enabling precision cosmology and extragalactic astrophysics with modern wide-field photometric surveys. Evaluating photo-z posteriors on a per-galaxy basis is difficult, however, as real galaxies have a true redshift but not a true redshift posterior. We introduce PZFlow, a Python package for the probabilistic forward modeling of galaxy catalogs with normalizing flows. For catalogs simulated with PZFlow, there is a natural notion of "true" redshift posteriors that can be used for photo-z validation. We use PZFlow to simulate a photometric galaxy catalog where each galaxy has a redshift, noisy photometry, shape information, and a true redshift posterior. We also demonstrate the use of an ensemble of normalizing flows for photo-z estimation. We discuss how PZFlow will be used to validate the photo-z estimation pipeline of the Dark Energy Science Collaboration (DESC), and the wider applicability of PZFlow for statistical modeling of any tabular data.
△ Less
Submitted 7 May, 2024;
originally announced May 2024.
-
Astrometric Redshifts of Supernovae
Authors:
Jaemyoung Jason Lee,
Masao Sako,
Richard Kessler,
Alex I. Malz,
The LSST Dark Energy Science Collaboration
Abstract:
Differential Chromatic Refraction (DCR) is caused by the wavelength dependence of our atmosphere's refractive index, which shifts the apparent positions of stars and galaxies and distorts their shapes depending on their spectral energy distributions (SEDs). While this effect is typically mitigated and corrected for in imaging observations, we investigate how DCR can instead be used to our advantag…
▽ More
Differential Chromatic Refraction (DCR) is caused by the wavelength dependence of our atmosphere's refractive index, which shifts the apparent positions of stars and galaxies and distorts their shapes depending on their spectral energy distributions (SEDs). While this effect is typically mitigated and corrected for in imaging observations, we investigate how DCR can instead be used to our advantage to infer the redshifts of supernovae from multi-band, time-series imaging data. We simulate Type Ia supernovae (SNe Ia) in the proposed Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) Deep Drilling Field (DDF), and evaluate astrometric redshifts. We find that the redshift accuracy improves dramatically with the statistical quality of the astrometric measurements as well as with the accuracy of the astrometric solution. For a conservative choice of a 5-mas systematic uncertainty floor, we find that our redshift estimation is accurate at $z < 0.6$. We then combine our astrometric redshifts with both host galaxy photometric redshifts and supernovae photometric (light-curve) redshifts and show that this considerably improves the overall redshift estimates. These astrometric redshifts will be valuable especially since Rubin will discover a vast number of supernovae for which we will not be able to obtain spectroscopic redshifts.
△ Less
Submitted 7 May, 2024;
originally announced May 2024.
-
Broadband Multi-wavelength Properties of M87 during the 2018 EHT Campaign including a Very High Energy Flaring Episode
Authors:
J. C. Algaba,
M. Balokovic,
S. Chandra,
W. Y. Cheong,
Y. Z. Cui,
F. D'Ammando,
A. D. Falcone,
N. M. Ford,
M. Giroletti,
C. Goddi,
M. A. Gurwell,
K. Hada,
D. Haggard,
S. Jorstad,
A. Kaur,
T. Kawashima,
S. Kerby,
J. Y. Kim,
M. Kino,
E. V. Kravchenko,
S. S. Lee,
R. S. Lu,
S. Markoff,
J. Michail,
J. Neilsen
, et al. (14 additional authors not shown)
Abstract:
The nearby elliptical galaxy M87 contains one of the only two supermassive black holes whose emission surrounding the event horizon has been imaged by the Event Horizon Telescope (EHT). In 2018, more than two dozen multi-wavelength (MWL) facilities (from radio to gamma-ray energies) took part in the second M87 EHT campaign. The goal of this extensive MWL campaign was to better understand the physi…
▽ More
The nearby elliptical galaxy M87 contains one of the only two supermassive black holes whose emission surrounding the event horizon has been imaged by the Event Horizon Telescope (EHT). In 2018, more than two dozen multi-wavelength (MWL) facilities (from radio to gamma-ray energies) took part in the second M87 EHT campaign. The goal of this extensive MWL campaign was to better understand the physics of the accreting black hole M87*, the relationship between the inflow and inner jets, and the high-energy particle acceleration. Understanding the complex astrophysics is also a necessary first step towards performing further tests of general relativity. The MWL campaign took place in April 2018, overlapping with the EHT M87* observations. We present a new, contemporaneous spectral energy distribution (SED) ranging from radio to very high energy (VHE) gamma-rays, as well as details of the individual observations and light curves. We also conduct phenomenological modelling to investigate the basic source properties. We present the first VHE gamma-ray flare from M87 detected since 2010. The flux above 350 GeV has more than doubled within a period of about 36 hours. We find that the X-ray flux is enhanced by about a factor of two compared to 2017, while the radio and millimetre core fluxes are consistent between 2017 and 2018. We detect evidence for a monotonically increasing jet position angle that corresponds to variations in the bright spot of the EHT image. Our results show the value of continued MWL monitoring together with precision imaging for addressing the origins of high-energy particle acceleration. While we cannot currently pinpoint the precise location where such acceleration takes place, the new VHE gamma-ray flare already presents a challenge to simple one-zone leptonic emission model approaches, and emphasises the need for combined image and spectral modelling.
△ Less
Submitted 9 September, 2024; v1 submitted 24 April, 2024;
originally announced April 2024.
-
Machine Learning LSST 3x2pt analyses -- forecasting the impact of systematics on cosmological constraints using neural networks
Authors:
Supranta S. Boruah,
Tim Eifler,
Vivian Miranda,
Elyas Farah,
Jay Motka,
Elisabeth Krause,
Xiao Fang,
Paul Rogozenski,
The LSST Dark Energy Science Collaboration
Abstract:
Validating modeling choices through simulated analyses and quantifying the impact of different systematic effects will form a major computational bottleneck in the preparation for 3$\times$2 analysis with Stage-IV surveys such as Vera Rubin Observatory's Legacy Survey of Space and Time (LSST). We can significantly reduce the computational requirements by using machine learning based emulators, whi…
▽ More
Validating modeling choices through simulated analyses and quantifying the impact of different systematic effects will form a major computational bottleneck in the preparation for 3$\times$2 analysis with Stage-IV surveys such as Vera Rubin Observatory's Legacy Survey of Space and Time (LSST). We can significantly reduce the computational requirements by using machine learning based emulators, which allow us to run fast inference while maintaining the full realism of the data analysis pipeline. In this paper, we use such an emulator to run simulated 3$\times$2 (cosmic shear, galaxy-galaxy lensing, and galaxy clustering) analyses for mock LSST-Y1/Y3/Y6/Y10 surveys and study the impact of various systematic effects (galaxy bias, intrinsic alignment, baryonic physics, shear calibration and photo-$z$ uncertainties). Closely following the DESC Science Requirement Document (with several updates) our main findings are: {\it a)} The largest contribution to the `systematic error budget' of LSST 3$\times$2 analysis comes from galaxy bias uncertainties, while the contribution of baryonic and shear calibration uncertainties are significantly less important. {\it b)} Tighter constraints on intrinsic alignment and photo-$z$ parameters can improve cosmological constraints noticeably, which illustrates synergies of LSST and spectroscopic surveys. {\it c)} The scale cuts adopted in the DESC SRD may be too conservative and pushing to smaller scales can increase cosmological information significantly. {\it d)} We investigate the impact of photo-$z$ outliers on 3$\times$2 pt analysis and find that we need to determine the outlier fraction to within $5-10\%$ accuracy to ensure robust cosmological analysis. We caution that these findings depend on analysis choices (parameterizations, priors, scale cuts) and can change for different settings.
△ Less
Submitted 18 March, 2024;
originally announced March 2024.
-
Recovered SN Ia rate from simulated LSST images
Authors:
Vincenzo Petrecca,
Maria Teresa Botticella,
Enrico Cappellaro,
Laura Greggio,
Bruno Sánchez,
Anais Möller,
Masao Sako,
Melissa Graham,
Maurizio Paolillo,
Federica Bianco,
the LSST Dark Energy Science Collaboration
Abstract:
The Legacy Survey of Space and Time (LSST) will revolutionize Time Domain Astronomy by detecting millions of transients. In particular, it is expected to increment the number of type Ia supernovae (SNIa) of a factor of 100 compared to existing samples up to z~1.2. Such a high number of events will dramatically reduce statistical uncertainties in the analysis of SNIa properties and rates. However,…
▽ More
The Legacy Survey of Space and Time (LSST) will revolutionize Time Domain Astronomy by detecting millions of transients. In particular, it is expected to increment the number of type Ia supernovae (SNIa) of a factor of 100 compared to existing samples up to z~1.2. Such a high number of events will dramatically reduce statistical uncertainties in the analysis of SNIa properties and rates. However, the impact of all other sources of uncertainty on the measurement must still be evaluated. The comprehension and reduction of such uncertainties will be fundamental both for cosmology and stellar evolution studies, as measuring the SNIa rate can put constraints on the evolutionary scenarios of different SNIa progenitors. We use simulated data from the DESC Data Challenge 2 (DC2) and LSST Data Preview 0 (DP0) to measure the SNIa rate on a 15 deg2 region of the Wide-Fast-Deep area. We select a sample of SN candidates detected on difference images, associate them to the host galaxy, and retrieve their photometric redshifts (z-phot). Then, we test different light curves classification methods, with and without redshift priors. We discuss how the distribution in redshift measured for the SN candidates changes according to the selected host galaxy and redshift estimate. We measure the SNIa rate analyzing the impact of uncertainties due to z-phot, host galaxy association and classification on the distribution in redshift of the starting sample. We found a 17% average lost fraction of SNIa with respect to the simulated sample. As 10% of the bias is due to the uncertainty on the z-phot alone (which also affects classification when used as a prior), it results to be the major source of uncertainty. We discuss possible reduction of the errors in the measurement of the SNIa rate, including synergies with other surveys, which may help using the rate to discriminate different progenitor models.
△ Less
Submitted 27 February, 2024;
originally announced February 2024.
-
Improving Photometric Redshift Estimates with Training Sample Augmentation
Authors:
Irene Moskowitz,
Eric Gawiser,
John Franklin Crenshaw,
Brett H. Andrews,
Alex I. Malz,
Samuel Schmidt,
The LSST Dark Energy Science Collaboration
Abstract:
Large imaging surveys will rely on photometric redshifts (photo-z's), which are typically estimated through machine learning methods. Currently planned spectroscopic surveys will not be deep enough to produce a representative training sample for LSST, so we seek methods to improve the photo-z estimates that arise from non-representative training samples. Spectroscopic training samples for photo-z'…
▽ More
Large imaging surveys will rely on photometric redshifts (photo-z's), which are typically estimated through machine learning methods. Currently planned spectroscopic surveys will not be deep enough to produce a representative training sample for LSST, so we seek methods to improve the photo-z estimates that arise from non-representative training samples. Spectroscopic training samples for photo-z's are biased towards redder, brighter galaxies, which also tend to be at lower redshift than the typical galaxy observed by LSST, leading to poor photo-z estimates with outlier fractions nearly 4 times larger than for a representative training sample. In this paper, we apply the concept of training sample augmentation, where we augment simulated non-representative training samples with simulated galaxies possessing otherwise unrepresented features. When we select simulated galaxies with (g-z) color, i-band magnitude and redshift outside the range of the original training sample, we are able to reduce the outlier fraction of the photo-z estimates for simulated LSST data by nearly 50% and the normalized median absolute deviation (NMAD) by 56%. When compared to a fully representative training sample, augmentation can recover nearly 70% of the degradation in the outlier fraction and 80% of the degradation in NMAD. Training sample augmentation is a simple and effective way to improve training samples for photo-z's without requiring additional spectroscopic samples.
△ Less
Submitted 14 May, 2024; v1 submitted 23 February, 2024;
originally announced February 2024.
-
Detecting strongly-lensed type Ia supernovae with LSST
Authors:
Nikki Arendse,
Suhail Dhawan,
Ana Sagués Carracedo,
Hiranya V. Peiris,
Ariel Goobar,
Radek Wojtak,
Catarina Alves,
Rahul Biswas,
Simon Huber,
Simon Birrer,
The LSST Dark Energy Science Collaboration
Abstract:
Strongly-lensed supernovae are rare and valuable probes of cosmology and astrophysics. Upcoming wide-field time-domain surveys, such as the Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST), are expected to discover an order-of-magnitude more lensed supernovae than have previously been observed. In this work, we investigate the cosmological prospects of lensed type Ia supernovae (…
▽ More
Strongly-lensed supernovae are rare and valuable probes of cosmology and astrophysics. Upcoming wide-field time-domain surveys, such as the Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST), are expected to discover an order-of-magnitude more lensed supernovae than have previously been observed. In this work, we investigate the cosmological prospects of lensed type Ia supernovae (SNIa) in LSST by quantifying the expected annual number of detections, the impact of stellar microlensing, follow-up feasibility, and how to best separate lensed and unlensed SNIa. We simulate SNIa lensed by galaxies, using the current LSST baseline v3.0 cadence, and find an expected number of 44 lensed SNIa detections per year. Microlensing effects by stars in the lens galaxy are predicted to lower the lensed SNIa detections by $\sim 8 \%$. The lensed events can be separated from the unlensed ones by jointly considering their colours and peak magnitudes. We define a `gold sample' of $\sim 10$ lensed SNIa per year with time delay $> 10$ days, $> 5$ detections before light-curve peak, and sufficiently bright ($m_i < 22.5$ mag) for follow-up observations. In three years of LSST operations, such a sample is expected to yield a $1.5\%$ measurement of the Hubble constant.
△ Less
Submitted 16 May, 2024; v1 submitted 7 December, 2023;
originally announced December 2023.
-
The impact of environment on size: Galaxies are 50% smaller in the Fornax Cluster compared to the field
Authors:
Nushkia Chamba,
Matthew Hayes,
The LSST Dark Energy Science Collaboration
Abstract:
Size is a fundamental parameter for measuring the growth of galaxies and the role of the environment on their evolution. However, the conventional size definitions used for this purpose are often biased and miss the diffuse, outermost signatures of galaxy growth, including star formation and gas accretion. This issue is addressed by examining low surface brightness truncations or galaxy ``edges''…
▽ More
Size is a fundamental parameter for measuring the growth of galaxies and the role of the environment on their evolution. However, the conventional size definitions used for this purpose are often biased and miss the diffuse, outermost signatures of galaxy growth, including star formation and gas accretion. This issue is addressed by examining low surface brightness truncations or galaxy ``edges'' as a physically motivated tracer of size based on star formation thresholds. Our total sample consists of $\sim900$ galaxies with stellar masses ranging from $10^5 M_{\odot} < M_{\star} < 10^{11} M_{\odot}$. This sample of nearby cluster, group satellite and nearly isolated field galaxies was compiled using multi-band imaging from the Fornax Deep Survey, deep IAC Stripe 82 and Dark Energy Camera Legacy Surveys. We find that the edge radii scale as $R_{\rm edge} \propto M_{\star}^{0.42}$ with a very small intrinsic scatter ($\sim 0.07$ dex). The scatter is driven by the morphology and environment of galaxies. In both the cluster and field, early-type dwarfs are systematically smaller by $\sim20\%$ than the late-types. However, compared to the field galaxies in the Fornax cluster are the most impacted. At a fixed stellar mass, edges in the cluster can be found at $\sim$ 50\% smaller radii and the average stellar surface density at the edges is a factor of two higher $\sim 1\,M_{\odot}$/pc$^2$. Our findings support the rapid removal of loosely bound neutral hydrogen in hot, crowded environments which truncates galaxies outside-in earlier, preventing the formation of more extended sizes and lower density edges. Our results highlight the importance of deep imaging surveys to study the low surface brightness imprints of the large scale structure and environment on galaxy evolution.
△ Less
Submitted 30 May, 2024; v1 submitted 16 November, 2023;
originally announced November 2023.
-
An Empirical Model For Intrinsic Alignments: Insights From Cosmological Simulations
Authors:
Nicholas Van Alfen,
Duncan Campbell,
Jonathan Blazek,
C. Danielle Leonard,
Francois Lanusse,
Andrew Hearin,
Rachel Mandelbaum,
The LSST Dark Energy Science Collaboration
Abstract:
We extend current models of the halo occupation distribution (HOD) to include a flexible, empirical framework for the forward modeling of the intrinsic alignment (IA) of galaxies. A primary goal of this work is to produce mock galaxy catalogs for the purpose of validating existing models and methods for the mitigation of IA in weak lensing measurements. This technique can also be used to produce n…
▽ More
We extend current models of the halo occupation distribution (HOD) to include a flexible, empirical framework for the forward modeling of the intrinsic alignment (IA) of galaxies. A primary goal of this work is to produce mock galaxy catalogs for the purpose of validating existing models and methods for the mitigation of IA in weak lensing measurements. This technique can also be used to produce new, simulation-based predictions for IA and galaxy clustering. Our model is probabilistically formulated, and rests upon the assumption that the orientations of galaxies exhibit a correlation with their host dark matter (sub)halo orientation or with their position within the halo. We examine the necessary components and phenomenology of such a model by considering the alignments between (sub)halos in a cosmological dark matter only simulation. We then validate this model for a realistic galaxy population in a set of simulations in the IllustrisTNG suite. We create an HOD mock with TNG-like correlations using our method, constraining the associated IA model parameters, with the $χ^2_{\rm dof}$ between our model's correlations and those of Illustris matching as closely as 1.4 and 1.1 for orientation--position and orientation--orientation correlation functions, respectively. By modeling the misalignment between galaxies and their host halo, we show that the 3-dimensional two-point position and orientation correlation functions of simulated (sub)halos and galaxies can be accurately reproduced from quasi-linear scales down to $0.1~h^{-1}{\rm Mpc}$. We also find evidence for environmental influence on IA within a halo. Our publicly-available software provides a key component enabling efficient determination of Bayesian posteriors on IA model parameters using observational measurements of galaxy-orientation correlation functions in the highly nonlinear regime.
△ Less
Submitted 11 September, 2024; v1 submitted 13 November, 2023;
originally announced November 2023.
-
Generation of realistic input parameters for simulating atmospheric point-spread functions at astronomical observatories
Authors:
Claire-Alice Hébert,
Joshua E. Meyers,
My H. Do,
Patricia R. Burchat,
the LSST Dark Energy Science Collaboration
Abstract:
High-fidelity simulated astronomical images are an important tool in developing and measuring the performance of image-processing algorithms, particularly for high precision measurements of cosmic shear -- correlated distortions of images of distant galaxies due to weak gravitational lensing caused by the large-scale mass distribution in the Universe. For unbiased measurements of cosmic shear, all…
▽ More
High-fidelity simulated astronomical images are an important tool in developing and measuring the performance of image-processing algorithms, particularly for high precision measurements of cosmic shear -- correlated distortions of images of distant galaxies due to weak gravitational lensing caused by the large-scale mass distribution in the Universe. For unbiased measurements of cosmic shear, all other sources of correlated image distortions must be modeled or removed. One such source is the correlated blurring of images due to optical turbulence in the atmosphere, which dominates the point-spread function (PSF) for ground-based instruments. In this work, we leverage data from weather forecasting models to produce wind speeds and directions, and turbulence parameters, that are realistically correlated with altitude. To study the resulting correlations in the size and shape of the PSF, we generate simulated images of the PSF across a ~10 square-degree field of view -- the size of the camera focal plane for the Vera C. Rubin Observatory in Chile -- using weather data and historical seeing for a geographic location near the Observatory. We make quantitative predictions for two-point correlation functions (2PCF) that are used in analyses of cosmic shear. We observe a strong anisotropy in the two-dimensional 2PCF, which is expected based on observations in real images, and study the dependence of the orientation of the anisotropy on dominant wind directions near the ground and at higher altitudes.
The code repository for producing the correlated weather parameters for input to simulations (psf-weather-station) is public at https://github.com/LSSTDESC/psf-weather-station.
△ Less
Submitted 21 March, 2024; v1 submitted 16 October, 2023;
originally announced October 2023.
-
Cosmology with imaging galaxy surveys: the impact of evolving galaxy bias and magnification
Authors:
Shivam Pandey,
Carles Sánchez,
Bhuvnesh Jain,
The LSST Dark Energy Science Collaboration
Abstract:
The joint analysis of galaxy clustering and galaxy-shear cross-correlations (galaxy-galaxy lensing) in imaging surveys constitutes one of the main avenues to obtain cosmological information. Analyses from Stage III surveys have assumed constant galaxy bias and weak lensing magnification coefficients within redshift bins. We use simple models of bias and magnification evolution to show that plausib…
▽ More
The joint analysis of galaxy clustering and galaxy-shear cross-correlations (galaxy-galaxy lensing) in imaging surveys constitutes one of the main avenues to obtain cosmological information. Analyses from Stage III surveys have assumed constant galaxy bias and weak lensing magnification coefficients within redshift bins. We use simple models of bias and magnification evolution to show that plausible levels of evolution impact cosmological parameter inference significantly for Stage IV surveys and possibly for current datasets as well. The distance weighting factors in galaxy-galaxy lensing cause it to have a different effective mean redshift within a lens bin than galaxy clustering. If the evolution is not included carefully in the theory modeling, this produces an apparent de-correlation of clustering and lensing even on large scales and leads to biased cosmological constraints. The level of bias is sensitive to the redshift and the width of the lens redshift distributions: for magnitude-limited-like samples it is comparable to the statistical precision even for the full DES survey, and exceeds it for Rubin-LSST. We find that ignoring the evolution can lead to approximately $1σ$ biases in the cosmological constraints in future Rubin-LSST analyses. We describe how these effects can be mitigated from the modeling and data sides, in particular by using cross-correlations to constrain the evolution parameters. We also demonstrate that the combination of galaxy clustering and CMB lensing does not suffer from this problem, due to the different distance weighting.
△ Less
Submitted 2 October, 2023;
originally announced October 2023.
-
A Joint Fermi-GBM and Swift-BAT Analysis of Gravitational-Wave Candidates from the Third Gravitational-wave Observing Run
Authors:
C. Fletcher,
J. Wood,
R. Hamburg,
P. Veres,
C. M. Hui,
E. Bissaldi,
M. S. Briggs,
E. Burns,
W. H. Cleveland,
M. M. Giles,
A. Goldstein,
B. A. Hristov,
D. Kocevski,
S. Lesage,
B. Mailyan,
C. Malacaria,
S. Poolakkil,
A. von Kienlin,
C. A. Wilson-Hodge,
The Fermi Gamma-ray Burst Monitor Team,
M. Crnogorčević,
J. DeLaunay,
A. Tohuvavohu,
R. Caputo,
S. B. Cenko
, et al. (1674 additional authors not shown)
Abstract:
We present Fermi Gamma-ray Burst Monitor (Fermi-GBM) and Swift Burst Alert Telescope (Swift-BAT) searches for gamma-ray/X-ray counterparts to gravitational wave (GW) candidate events identified during the third observing run of the Advanced LIGO and Advanced Virgo detectors. Using Fermi-GBM on-board triggers and sub-threshold gamma-ray burst (GRB) candidates found in the Fermi-GBM ground analyses,…
▽ More
We present Fermi Gamma-ray Burst Monitor (Fermi-GBM) and Swift Burst Alert Telescope (Swift-BAT) searches for gamma-ray/X-ray counterparts to gravitational wave (GW) candidate events identified during the third observing run of the Advanced LIGO and Advanced Virgo detectors. Using Fermi-GBM on-board triggers and sub-threshold gamma-ray burst (GRB) candidates found in the Fermi-GBM ground analyses, the Targeted Search and the Untargeted Search, we investigate whether there are any coincident GRBs associated with the GWs. We also search the Swift-BAT rate data around the GW times to determine whether a GRB counterpart is present. No counterparts are found. Using both the Fermi-GBM Targeted Search and the Swift-BAT search, we calculate flux upper limits and present joint upper limits on the gamma-ray luminosity of each GW. Given these limits, we constrain theoretical models for the emission of gamma-rays from binary black hole mergers.
△ Less
Submitted 25 August, 2023;
originally announced August 2023.
-
Large Interferometer For Exoplanets (LIFE). X. Detectability of currently known exoplanets and synergies with future IR/O/UV reflected-starlight imaging missions
Authors:
Óscar Carrión-González,
Jens Kammerer,
Daniel Angerhausen,
Felix Dannert,
Antonio García Muñoz,
Sascha P. Quanz,
Olivier Absil,
Charles A. Beichman,
Julien H. Girard,
Bertrand Mennesson,
Michael R. Meyer,
Karl R. Stapelfeldt,
The LIFE Collaboration
Abstract:
The next generation of space-based observatories will characterize the atmospheres of low-mass, temperate exoplanets with the direct-imaging technique. This will be a major step forward in our understanding of exoplanet diversity and the prevalence of potentially habitable conditions beyond the Earth. We compute a list of currently known exoplanets detectable with the mid-infrared Large Interferom…
▽ More
The next generation of space-based observatories will characterize the atmospheres of low-mass, temperate exoplanets with the direct-imaging technique. This will be a major step forward in our understanding of exoplanet diversity and the prevalence of potentially habitable conditions beyond the Earth. We compute a list of currently known exoplanets detectable with the mid-infrared Large Interferometer For Exoplanets (LIFE) in thermal emission. We also compute the list of known exoplanets accessible to a notional design of the Habitable Worlds Observatory (HWO), observing in reflected starlight. With a pre-existing method, we processed the NASA Exoplanet Archive and computed orbital realizations for each known exoplanet. We derived their mass, radius, equilibrium temperature, and planet-star angular separation. We used the LIFEsim simulator to compute the integration time ($t_{int}$) required to detect each planet with LIFE. A planet is considered detectable if a broadband signal-to-noise ratio $S/N$=7 is achieved over the spectral range $4-18.5μ$m in $t_{int}\leq$100 hours. We tested whether the planet is accessible to HWO in reflected starlight based on its notional inner and outer working angles, and minimum planet-to-star contrast. LIFE's reference configuration (four 2-m telescopes with 5% throughput and a nulling baseline between 10-100 m) can detect 212 known planets within 20 pc. Of these, 55 are also accessible to HWO in reflected starlight, offering a unique opportunity for synergies in atmospheric characterization. LIFE can also detect 32 known transiting exoplanets. Furthermore, 38 LIFE-detectable planets orbit in the habitable zone, of which 13 with $M_p<5M_\oplus$ and 8 with $5M_\oplus<M_p<10M_\oplus$. LIFE already has enough targets to perform ground-breaking analyses of low-mass, habitable-zone exoplanets, a fraction of which will also be accessible to other instruments.
△ Less
Submitted 18 August, 2023;
originally announced August 2023.
-
Large Interferometer For Exoplanets (LIFE): XI. Phase-space synthesis decomposition for planet detection and characterization
Authors:
Taro Matsuo,
Felix Dannert,
Romain Laugier,
Sascha P. Quanz,
Andjelka B. Kovacevic,
LIFE collaboration
Abstract:
A mid-infrared nulling-space interferometer is a promising way to characterize thermal light from habitable planet candidates around Sun-like stars. However, one of the main challenges for achieving this ambitious goal is a high-precision stability of the optical path difference (OPD) and amplitude over a few days for planet detection and up to a few weeks for in-depth characterization. Here we pr…
▽ More
A mid-infrared nulling-space interferometer is a promising way to characterize thermal light from habitable planet candidates around Sun-like stars. However, one of the main challenges for achieving this ambitious goal is a high-precision stability of the optical path difference (OPD) and amplitude over a few days for planet detection and up to a few weeks for in-depth characterization. Here we propose a new method called phase-space synthesis decomposition (PSSD) to shorten the stability requirement to minutes, significantly relaxing the technological challenges of the mission. Focusing on what exactly modulates the planet signal in the presence of the stellar leak and systematic error, PSSD prioritizes the modulation of the signals along the wavelength domain rather than baseline rotation. Modulation along the wavelength domain allows us to extract source positions in parallel to the baseline vector for each exposure. The sum of the one-dimensional data converts into two-dimensional information. Based on the reconstructed image, we construct a continuous equation and extract the spectra through the singular value decomposition (SVD) while efficiently separating them from a long-term systematic stellar leak. We performed numerical simulations to investigate the feasibility of PSSD for the LIFE mission concept. We confirm that multiple terrestrial planets in the habitable zone around a Sun-like star at 10 pc can be detected and characterized despite high levels and long durations of systematic noise. We also find that PSSD is more robust against a sparse sampling of the array rotation compared to purely rotation-based signal extraction. Using PSSD as signal extraction method significantly relaxes the technical requirements on signal stability and further increases the feasibility of the LIFE mission.
△ Less
Submitted 2 August, 2023;
originally announced August 2023.
-
Slitless spectrophotometry with forward modelling: principles and application to atmospheric transmission measurement
Authors:
Jérémy Neveu,
Vincent Brémaud,
Pierre Antilogus,
Florent Barret,
Sébastien Bongard,
Yannick Copin,
Sylvie Dagoret-Campagne,
Claire Juramy,
Laurent Le-Guillou,
Marc Moniez,
Eduardo Sepulveda,
The LSST Dark Energy Science Collaboration
Abstract:
In the next decade, many optical surveys will aim to tackle the question of dark energy nature, measuring its equation of state parameter at the permil level. This requires trusting the photometric calibration of the survey with a precision never reached so far, controlling many sources of systematic uncertainties. The measurement of the on-site atmospheric transmission for each exposure, or on av…
▽ More
In the next decade, many optical surveys will aim to tackle the question of dark energy nature, measuring its equation of state parameter at the permil level. This requires trusting the photometric calibration of the survey with a precision never reached so far, controlling many sources of systematic uncertainties. The measurement of the on-site atmospheric transmission for each exposure, or on average for each season or for the full survey, can help reach the permil precision for magnitudes. This work aims at proving the ability to use slitless spectroscopy for standard star spectrophotometry and its use to monitor on-site atmospheric transmission as needed, for example, by the Vera C. Rubin Observatory Legacy Survey of Space and Time supernova cosmology program. We fully deal with the case of a disperser in the filter wheel, which is the configuration chosen in the Rubin Auxiliary Telescope. The theoretical basis of slitless spectrophotometry is at the heart of our forward model approach to extract spectroscopic information from slitless data. We developed a publicly available software called Spectractor (https://github.com/LSSTDESC/Spectractor) that implements each ingredient of the model and finally performs a fit of a spectrogram model directly on image data to get the spectrum. We show on simulations that our model allows us to understand the structure of spectrophotometric exposures. We also demonstrate its use on real data, solving specific issues and illustrating how our procedure allows the improvement of the model describing the data. Finally, we discuss how this approach can be used to directly extract atmospheric transmission parameters from data and thus provide the base for on-site atmosphere monitoring. We show the efficiency of the procedure on simulations and test it on the limited data set available.
△ Less
Submitted 10 June, 2024; v1 submitted 10 July, 2023;
originally announced July 2023.
-
Fat cosmic ray tracks in charge-coupled devices
Authors:
Theodore A. Grosson,
Andrei Nomerotski,
The LSST Dark Energy Science Collaboration
Abstract:
Cosmic rays are particles from the upper atmosphere which often leave bright spots and trails in images from telescope CCDs. We investigate so-called ``fat" cosmic rays seen in images from Vera C. Rubin Observatory and the Subaru Telescope. These tracks are much wider and brighter than typical cosmic ray tracks, and therefore are more capable of obscuring data in science images. By understanding t…
▽ More
Cosmic rays are particles from the upper atmosphere which often leave bright spots and trails in images from telescope CCDs. We investigate so-called ``fat" cosmic rays seen in images from Vera C. Rubin Observatory and the Subaru Telescope. These tracks are much wider and brighter than typical cosmic ray tracks, and therefore are more capable of obscuring data in science images. By understanding the origins of these tracks, we can better ensure that they do not interfere with on-sky data. We compare the properties of these tracks to simulated and theoretical models in order to identify both the particles causing these tracks as well as the reason for their excess spread. We propose that the origin of these tracks is cosmic ray protons, which deposit much greater charge in the CCDs than typical cosmic rays due to their lower velocities. The generated charges then repel each other while drifting through the detector, resulting in a track which is much wider than typical tracks.
△ Less
Submitted 19 June, 2023;
originally announced June 2023.
-
Are classification metrics good proxies for SN Ia cosmological constraining power?
Authors:
Alex I. Malz,
Mi Dai,
Kara A. Ponder,
Emille E. O. Ishida,
Santiago Gonzalez-Gaitain,
Rupesh Durgesh,
Alberto Krone-Martins,
Rafael S. de Souza,
Noble Kennamer,
Sreevarsha Sreejith,
Lluis Galbany,
The LSST Dark Energy Science Collaboration,
The Cosmostatistics Initiative
Abstract:
Context: When selecting a classifier to use for a supernova Ia (SN Ia) cosmological analysis, it is common to make decisions based on metrics of classification performance, i.e. contamination within the photometrically classified SN Ia sample, rather than a measure of cosmological constraining power. If the former is an appropriate proxy for the latter, this practice would save those designing an…
▽ More
Context: When selecting a classifier to use for a supernova Ia (SN Ia) cosmological analysis, it is common to make decisions based on metrics of classification performance, i.e. contamination within the photometrically classified SN Ia sample, rather than a measure of cosmological constraining power. If the former is an appropriate proxy for the latter, this practice would save those designing an analysis pipeline from the computational expense of a full cosmology forecast. Aims: This study tests the assumption that classification metrics are an appropriate proxy for cosmology metrics. Methods: We emulate photometric SN Ia cosmology samples with controlled contamination rates of individual contaminant classes and evaluate each of them under a set of classification metrics. We then derive cosmological parameter constraints from all samples under two common analysis approaches and quantify the impact of contamination by each contaminant class on the resulting cosmological parameter estimates. Results: We observe that cosmology metrics are sensitive to both the contamination rate and the class of the contaminating population, whereas the classification metrics are insensitive to the latter. Conclusions: We therefore discourage exclusive reliance on classification-based metrics for cosmological analysis design decisions, e.g. classifier choice, and instead recommend optimizing using a metric of cosmological parameter constraining power.
△ Less
Submitted 23 May, 2023;
originally announced May 2023.
-
Forecasting the power of Higher Order Weak Lensing Statistics with automatically differentiable simulations
Authors:
Denise Lanzieri,
François Lanusse,
Chirag Modi,
Benjamin Horowitz,
Joachim Harnois-Déraps,
Jean-Luc Starck,
The LSST Dark Energy Science Collaboration
Abstract:
We present the Differentiable Lensing Lightcone (DLL), a fully differentiable physical model designed for being used as a forward model in Bayesian inference algorithms requiring access to derivatives of lensing observables with respect to cosmological parameters. We extend the public FlowPM N-body code, a particle-mesh N-body solver, simulating lensing lightcones and implementing the Born approxi…
▽ More
We present the Differentiable Lensing Lightcone (DLL), a fully differentiable physical model designed for being used as a forward model in Bayesian inference algorithms requiring access to derivatives of lensing observables with respect to cosmological parameters. We extend the public FlowPM N-body code, a particle-mesh N-body solver, simulating lensing lightcones and implementing the Born approximation in the Tensorflow framework. Furthermore, DLL is aimed at achieving high accuracy with low computational costs. As such, it integrates a novel Hybrid Physical-Neural parameterisation able to compensate for the small-scale approximations resulting from particle-mesh schemes for cosmological N-body simulations. We validate our simulations in an LSST setting against high-resolution $κ$TNG simulations by comparing both the lensing angular power spectrum and multiscale peak counts. We demonstrate an ability to recover lensing $C_\ell$ up to a 10% accuracy at $\ell=1000$ for sources at redshift 1, with as few as $\sim 0.6$ particles per Mpc/h. As a first use case, we use this tool to investigate the relative constraining power of the angular power spectrum and peak counts statistic in an LSST setting. Such comparisons are typically very costly as they require a large number of simulations, and do not scale well with the increasing number of cosmological parameters. As opposed to forecasts based on finite differences, these statistics can be analytically differentiated with respect to cosmology, or any systematics included in the simulations at the same computational cost of the forward simulation. We find that the peak counts outperform the power spectrum on the cold dark matter parameter $Ω_c$, on the amplitude of density fluctuations $σ_8$, and on the amplitude of the intrinsic alignment signal $A_{IA}$.
△ Less
Submitted 12 May, 2023;
originally announced May 2023.
-
Stellar Parameters and Chemical Abundances Estimated from LAMOST-II DR8 MRS based on Cycle-StarNet
Authors:
Rui Wang,
A-Li Luo,
Shuo Zhang,
Yuan-Sen Ting,
Teaghan O'Briain,
LAMOST MRS Collaboration
Abstract:
Deriving stellar atmospheric parameters and chemical abundances from stellar spectra is crucial for understanding the evolution of the Milky Way. By performing a fitting with MARCS model atmospheric theoretical synthetic spectra combined with a domain-adaptation method, we estimate the fundamental stellar parameters (Teff, log g, [Fe/H], vmic, and vmac) and 11 chemical abundances for 1.38 million…
▽ More
Deriving stellar atmospheric parameters and chemical abundances from stellar spectra is crucial for understanding the evolution of the Milky Way. By performing a fitting with MARCS model atmospheric theoretical synthetic spectra combined with a domain-adaptation method, we estimate the fundamental stellar parameters (Teff, log g, [Fe/H], vmic, and vmac) and 11 chemical abundances for 1.38 million FGKM-type stars of the Medium-Resolution Spectroscopic Survey (MRS) from LAMOST-II DR8. The domain-adaptation method, Cycle-StarNet, is employed to reduce the gap between observed and synthetic spectra, and the L-BFGS algorithm is used to search for the best-fit synthetic spectra. By combining the 2MASS photometric survey data, Gaia EDR3 parallax, and MIST isochrones, the surface gravities of the stars are constrained after estimating their bolometric luminosities. The accuracy of Teff, log g, and [Fe/H] can reach 150 K, 0.11 dex, and 0.15 dex, evaluated by the PASTEL catalog, asteroseismic samples, and other spectroscopic surveys. The precision of these parameters and elemental abundances ([C/Fe], [Na/Fe], [Mg/Fe], [Si/Fe], [Ca/Fe], [Ti/Fe], [Cr/Fe], [Mn/Fe], [Co/Fe], [Ni/Fe], and [Cu/Fe]) is assessed by repeated observations and validated by cluster members. For spectra with signal-to-noise (S/N) ratios greater than 10, the precision of the three stellar parameters and elemental abundances can achieve 76 K, 0.014 dex, 0.096 dex, and 0.04-0.15 dex. For spectra with S/N ratios higher than 100, the precision stabilizes at 22 K, 0.006 dex, 0.043 dex, and 0.01-0.06 dex. The full LAMOST MRS stellar properties catalog is available online.
△ Less
Submitted 9 May, 2023;
originally announced May 2023.
-
Large Interferometer For Exoplanets (LIFE): IX. Assessing the Impact of Clouds on Atmospheric Retrievals at Mid-Infrared Wavelengths with a Venus-Twin Exoplanet
Authors:
B. S. Konrad,
E. Alei,
S. P. Quanz,
P. Mollière,
D. Angerhausen,
J. J. Fortney,
K. Hakim,
S. Jordan,
D. Kitzmann,
S. Rugheimer,
O. Shorttle,
R. Wordsworth,
the LIFE Collaboration
Abstract:
The Large Interferometer For Exoplanets (LIFE) initiative aims to develop a space based mid-infrared (MIR) nulling interferometer to measure the thermal emission spectra of temperate terrestrial exoplanets.
We investigate how well LIFE could characterize a cloudy Venus-twin exoplanet to: (1) test our retrieval routine on a realistic non-Earth-like MIR spectrum of a known planet, (2) investigate…
▽ More
The Large Interferometer For Exoplanets (LIFE) initiative aims to develop a space based mid-infrared (MIR) nulling interferometer to measure the thermal emission spectra of temperate terrestrial exoplanets.
We investigate how well LIFE could characterize a cloudy Venus-twin exoplanet to: (1) test our retrieval routine on a realistic non-Earth-like MIR spectrum of a known planet, (2) investigate how clouds impact retrievals, (3) refine the LIFE requirements derived in previous Earth-centered studies.
We run retrievals for simulated LIFE observations of a Venus-twin exoplanet orbiting a Sun-like star located 10 pc from the observer. By assuming different models (cloudy and cloud-free) we analyze the performance as a function of the quality of the LIFE observation. This allows us to determine how well atmosphere and clouds are characterizable depending on the quality of the spectrum.
Our study shows that the current minimal resolution ($R=50$) and signal-to-noise ($S/N=10$ at $11.2μ$m) requirements for LIFE suffice to characterize the structure and composition of a Venus-like atmosphere above the cloud deck if an adequate model is chosen. However, we cannot infer cloud properties. The accuracy of the retrieved planet radius ($R_{pl}$), equilibrium temperature ($T_{eq}$), and Bond albedo ($A_B$) depend on the choice of model. Generally, a cloud-free model performs best and thus the presence of clouds cannot be inferred. This model dependence of retrieval results emphasizes the importance of developing a community-wide best-practice for atmospheric retrieval studies. If we consider higher quality spectra (especially $S/N=20$), we can infer the presence of clouds and pose first constraints on their structure.
△ Less
Submitted 21 March, 2023; v1 submitted 8 March, 2023;
originally announced March 2023.
-
Metadetection Weak Lensing for the Vera C. Rubin Observatory
Authors:
Erin S. Sheldon,
Matthew R. Becker,
Michael Jarvis,
Robert Armstrong,
The LSST Dark Energy Science Collaboration
Abstract:
Forthcoming astronomical imaging surveys will use weak gravitational lensing shear as a primary probe to study dark energy, with accuracy requirements at the 0.1% level. We present an implementation of the Metadetection shear measurement algorithm for use with the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST). This new code works with the data products produced by the LSST Scien…
▽ More
Forthcoming astronomical imaging surveys will use weak gravitational lensing shear as a primary probe to study dark energy, with accuracy requirements at the 0.1% level. We present an implementation of the Metadetection shear measurement algorithm for use with the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST). This new code works with the data products produced by the LSST Science Pipelines, and uses the pipeline algorithms when possible. We tested the code using a new set of simulations designed to mimic LSST imaging data. The simulated images contained semi-realistic galaxies, stars with representative distributions of magnitudes and galactic spatial density, cosmic rays, bad CCD columns and spatially variable point spread functions. Bright stars were saturated and simulated ``bleed trails'' were drawn. Problem areas were interpolated, and the images were coadded into small cells, excluding images not fully covering the cell to guarantee a continuous point spread function. In all our tests the measured shear was accurate within the LSST requirements.
△ Less
Submitted 11 May, 2023; v1 submitted 6 March, 2023;
originally announced March 2023.
-
Measurement of telescope transmission using a Collimated Beam Projector
Authors:
Nicholas Mondrik,
Michael Coughlin,
Marc Betoule,
Sébastien Bongard,
Joseph P. Rice,
Ping-Shine Shaw,
Christopher W. Stubbs,
John T. Woodward,
LSST Dark Energy Science Collaboration
Abstract:
With the increasingly large number of type Ia supernova being detected by current-generation survey telescopes, and even more expected with the upcoming Rubin Observatory Legacy Survey of Space and Time, the precision of cosmological measurements will become limited by systematic uncertainties in flux calibration rather than statistical noise. One major source of systematic error in determining SN…
▽ More
With the increasingly large number of type Ia supernova being detected by current-generation survey telescopes, and even more expected with the upcoming Rubin Observatory Legacy Survey of Space and Time, the precision of cosmological measurements will become limited by systematic uncertainties in flux calibration rather than statistical noise. One major source of systematic error in determining SNe Ia color evolution (needed for distance estimation) is uncertainty in telescope transmission, both within and between surveys. We introduce here the Collimated Beam Projector (CBP), which is meant to measure a telescope transmission with collimated light. The collimated beam more closely mimics a stellar wavefront as compared to flat-field based instruments, allowing for more precise handling of systematic errors such as those from ghosting and filter angle-of-incidence dependence. As a proof of concept, we present CBP measurements of the StarDICE prototype telescope, achieving a standard (1 sigma) uncertainty of 3 % on average over the full wavelength range measured with a single beam illumination.
△ Less
Submitted 22 February, 2023;
originally announced February 2023.
-
The catalog-to-cosmology framework for weak lensing and galaxy clustering for LSST
Authors:
J. Prat,
J. Zuntz,
Y. Omori,
C. Chang,
T. Tröster,
E. Pedersen,
C. García-García,
E. Phillips-Longley,
J. Sanchez,
D. Alonso,
X. Fang,
E. Gawiser,
K. Heitmann,
M. Ishak,
M. Jarvis,
E. Kovacs,
P. Larsen,
Y. -Y. Mao,
L. Medina Varela,
M. Paterno,
S. D. Vitenti,
Z. Zhang,
The LSST Dark Energy Science Collaboration
Abstract:
We present TXPipe, a modular, automated and reproducible pipeline for ingesting catalog data and performing all the calculations required to obtain quality-assured two-point measurements of lensing and clustering, and their covariances, with the metadata necessary for parameter estimation. The pipeline is developed within the Rubin Observatory Legacy Survey of Space and Time (LSST) Dark Energy Sci…
▽ More
We present TXPipe, a modular, automated and reproducible pipeline for ingesting catalog data and performing all the calculations required to obtain quality-assured two-point measurements of lensing and clustering, and their covariances, with the metadata necessary for parameter estimation. The pipeline is developed within the Rubin Observatory Legacy Survey of Space and Time (LSST) Dark Energy Science Collaboration (DESC), and designed for cosmology analyses using LSST data. In this paper, we present the pipeline for the so-called 3x2pt analysis -- a combination of three two-point functions that measure the auto- and cross-correlation between galaxy density and shapes. We perform the analysis both in real and harmonic space using TXPipe and other LSST-DESC tools. We validate the pipeline using Gaussian simulations and show that it accurately measures data vectors and recovers the input cosmology to the accuracy level required for the first year of LSST data under this simplified scenario. We also apply the pipeline to a realistic mock galaxy sample extracted from the CosmoDC2 simulation suite (Korytov et al. 2019). TXPipe establishes a baseline framework that can be built upon as the LSST survey proceeds. Furthermore, the pipeline is designed to be easily extended to science probes beyond the 3x2pt analysis.
△ Less
Submitted 21 April, 2023; v1 submitted 19 December, 2022;
originally announced December 2022.
-
Improved Tomographic Binning of 3x2pt Lens Samples: Neural Network Classifiers and Optimal Bin Assignments
Authors:
Irene Moskowitz,
Eric Gawiser,
Abby Bault,
Adam Broussard,
Jeffrey A. Newman,
Joe Zuntz,
the LSST Dark Energy Science Collaboration
Abstract:
Large imaging surveys, such as the Legacy Survey of Space and Time, rely on photometric redshifts and tomographic binning for 3x2pt analyses that combine galaxy clustering and weak lensing. In this paper, we propose a method for optimizing the tomographic binning choice for the lens sample of galaxies. We divide the CosmoDC2 and Buzzard simulated galaxy catalogs into a training set and an applicat…
▽ More
Large imaging surveys, such as the Legacy Survey of Space and Time, rely on photometric redshifts and tomographic binning for 3x2pt analyses that combine galaxy clustering and weak lensing. In this paper, we propose a method for optimizing the tomographic binning choice for the lens sample of galaxies. We divide the CosmoDC2 and Buzzard simulated galaxy catalogs into a training set and an application set, where the training set is nonrepresentative in a realistic way, and then estimate photometric redshifts for the application sets. The galaxies are sorted into redshift bins covering equal intervals of redshift or comoving distance, or with an equal number of galaxies in each bin, and we consider a generalized extension of these approaches. We find that bins of equal comoving distance produce the highest dark energy figure of merit of the initial binning choices, but that the choice of bin edges can be further optimized. We then train a neural network classifier to identify galaxies that are either highly likely to have accurate photometric redshift estimates or highly likely to be sorted into the correct redshift bin. The neural network classifier is used to remove poor redshift estimates from the sample, and the results are compared to the case when none of the sample is removed. We find that the neural network classifiers are able to improve the figure of merit by ~13% and are able to recover ~25% of the loss in the figure of merit that occurs when a nonrepresentative training sample is used.
△ Less
Submitted 14 June, 2023; v1 submitted 13 December, 2022;
originally announced December 2022.
-
The N5K Challenge: Non-Limber Integration for LSST Cosmology
Authors:
C. D. Leonard,
T. Ferreira,
X. Fang,
R. Reischke,
N. Schoeneberg,
T. Tröster,
D. Alonso,
J. E. Campagne,
F. Lanusse,
A. Slosar,
M. Ishak,
the LSST Dark Energy Science Collaboration
Abstract:
The rapidly increasing statistical power of cosmological imaging surveys requires us to reassess the regime of validity for various approximations that accelerate the calculation of relevant theoretical predictions. In this paper, we present the results of the 'N5K non-Limber integration challenge', the goal of which was to quantify the performance of different approaches to calculating the angula…
▽ More
The rapidly increasing statistical power of cosmological imaging surveys requires us to reassess the regime of validity for various approximations that accelerate the calculation of relevant theoretical predictions. In this paper, we present the results of the 'N5K non-Limber integration challenge', the goal of which was to quantify the performance of different approaches to calculating the angular power spectrum of galaxy number counts and cosmic shear data without invoking the so-called 'Limber approximation', in the context of the Rubin Observatory Legacy Survey of Space and Time (LSST). We quantify the performance, in terms of accuracy and speed, of three non-Limber implementations: ${\tt FKEM (CosmoLike)}$, ${\tt Levin}$, and ${\tt matter}$, themselves based on different integration schemes and approximations. We find that in the challenge's fiducial 3x2pt LSST Year 10 scenario, ${\tt FKEM (CosmoLike)}$ produces the fastest run time within the required accuracy by a considerable margin, positioning it favourably for use in Bayesian parameter inference. This method, however, requires further development and testing to extend its use to certain analysis scenarios, particularly those involving a scale-dependent growth rate. For this and other reasons discussed herein, alternative approaches such as ${\tt matter}$ and ${\tt Levin}$ may be necessary for a full exploration of parameter space. We also find that the usual first-order Limber approximation is insufficiently accurate for LSST Year 10 3x2pt analysis on $\ell=200-1000$, whereas invoking the second-order Limber approximation on these scales (with a full non-Limber method at smaller $\ell$) does suffice.
△ Less
Submitted 14 February, 2023; v1 submitted 8 December, 2022;
originally announced December 2022.
-
Fringing Analysis and Simulation for the Vera C. Rubin Observatory's Legacy Survey of Space and Time
Authors:
Zhiyuan Guo,
Chris W. Walter,
Craig Lage,
Robert H. Lupton,
The LSST Dark Energy Science Collaboration
Abstract:
The presence of fringing in astronomical CCD images will impact photometric quality and measurements. Yet its impact on the Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST) has not been fully studied. We present a detailed study on fringing for Charge-Coupled Devices (CCDs) already implemented on the Rubin Observatory LSST Camera's focal plane. After making physical measurements…
▽ More
The presence of fringing in astronomical CCD images will impact photometric quality and measurements. Yet its impact on the Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST) has not been fully studied. We present a detailed study on fringing for Charge-Coupled Devices (CCDs) already implemented on the Rubin Observatory LSST Camera's focal plane. After making physical measurements and knowing the compositions, we have developed a model for the e2v CCDs. We present a method to fit for the internal height variation of the epoxy layer within the sensors based on fringing measurements in a laboratory setting. This method is generic enough that it can be easily modified to work for other CCDs. Using the derived fringing model, we successfully reproduce comparable fringing amplitudes that match the observed levels in images taken by existing telescopes with different optical designs. This model is then used to forecast the expected level of fringing in a single LSST y-band sky background exposure with Rubin telescope optics in the presence of a realistic time varying sky spectrum. The predicted fringing amplitude in LSST images ranges from $0.04\%$ to $0.2\%$ depending on the location of a CCD on the focal plane. We find that the predicted variation in surface brightness caused by fringing in LSST y-band skybackground images is about $0.6\ μ\rm{Jy}\ \rm{arcsec}^{-2}$, which is 40 times larger than the current measurement error. We conclude that it is necessary to include fringing correction in the Rubin's LSST image processing pipeline.
△ Less
Submitted 16 November, 2022;
originally announced November 2022.
-
Hierarchical Inference of the Lensing Convergence from Photometric Catalogs with Bayesian Graph Neural Networks
Authors:
Ji Won Park,
Simon Birrer,
Madison Ueland,
Miles Cranmer,
Adriano Agnello,
Sebastian Wagner-Carena,
Philip J. Marshall,
Aaron Roodman,
the LSST Dark Energy Science Collaboration
Abstract:
We present a Bayesian graph neural network (BGNN) that can estimate the weak lensing convergence ($κ$) from photometric measurements of galaxies along a given line of sight. The method is of particular interest in strong gravitational time delay cosmography (TDC), where characterizing the "external convergence" ($κ_{\rm ext}$) from the lens environment and line of sight is necessary for precise in…
▽ More
We present a Bayesian graph neural network (BGNN) that can estimate the weak lensing convergence ($κ$) from photometric measurements of galaxies along a given line of sight. The method is of particular interest in strong gravitational time delay cosmography (TDC), where characterizing the "external convergence" ($κ_{\rm ext}$) from the lens environment and line of sight is necessary for precise inference of the Hubble constant ($H_0$). Starting from a large-scale simulation with a $κ$ resolution of $\sim$1$'$, we introduce fluctuations on galaxy-galaxy lensing scales of $\sim$1$''$ and extract random sightlines to train our BGNN. We then evaluate the model on test sets with varying degrees of overlap with the training distribution. For each test set of 1,000 sightlines, the BGNN infers the individual $κ$ posteriors, which we combine in a hierarchical Bayesian model to yield constraints on the hyperparameters governing the population. For a test field well sampled by the training set, the BGNN recovers the population mean of $κ$ precisely and without bias, resulting in a contribution to the $H_0$ error budget well under 1\%. In the tails of the training set with sparse samples, the BGNN, which can ingest all available information about each sightline, extracts more $κ$ signal compared to a simplified version of the traditional method based on matching galaxy number counts, which is limited by sample variance. Our hierarchical inference pipeline using BGNNs promises to improve the $κ_{\rm ext}$ characterization for precision TDC. The implementation of our pipeline is available as a public Python package, Node to Joy.
△ Less
Submitted 14 November, 2022;
originally announced November 2022.
-
Large Interferometer For Exoplanets (LIFE): VIII. Where is the phosphine? Observing exoplanetary PH3 with a space based MIR nulling interferometer
Authors:
D. Angerhausen,
M. Ottiger,
F. Dannert,
Y. Miguel,
C. Sousa-Silva,
J. Kammerer,
F. Menti,
E. Alei,
B. S. Konrad,
H. S. Wang,
S. P. Quanz,
the LIFE collaboration
Abstract:
Phosphine could be a key molecule in the understanding of exotic chemistry happening in (exo)planetary atmospheres. While it has been detected in the Solar System's giant planets, it has not been observed in exoplanets yet. In the exoplanetary context however it has been theorized as a potential biosignature molecule. The goal of our study is to identify which illustrative science cases for PH3 ch…
▽ More
Phosphine could be a key molecule in the understanding of exotic chemistry happening in (exo)planetary atmospheres. While it has been detected in the Solar System's giant planets, it has not been observed in exoplanets yet. In the exoplanetary context however it has been theorized as a potential biosignature molecule. The goal of our study is to identify which illustrative science cases for PH3 chemistry are observable with a space-based mid-infrared nulling interferometric observatory like the LIFE (Large Interferometer For Exoplanets) concept. We identified a representative set of scenarios for PH3 detections in exoplanetary atmospheres varying over the whole dynamic range of the LIFE mission. We used chemical kinetics and radiative transfer calculations to produce forward models of these informative, prototypical observational cases for LIFEsim, our observation simulator software for LIFE. In a detailed, yet first order approximation it takes a mission like LIFE: (i) about 1h to find phosphine in a warm giant around a G star at 10 pc, (ii) about 10 h in H2 or CO2 dominated temperate super-Earths around M star hosts at 5 pc, (iii) and even in 100h it seems very unlikely that phosphine would be detectable in a Venus-Twin with extreme PH3 concentrations at 5 pc. Phosphine in concentrations previously discussed in the literature is detectable in 2 out of the 3 cases and about an order of magnitude faster than comparable cases with JWST. We show that there is a significant number of objects accessible for these classes of observations. These results will be used to prioritize the parameter range for the next steps with more detailed retrieval simulations. They will also inform timely questions in the early design phase of a mission like LIFE and guide the community by providing easy-to-scale first estimates for a large part of detection space of such a mission.
△ Less
Submitted 9 November, 2022;
originally announced November 2022.
-
Using Host Galaxy Photometric Redshifts to Improve Cosmological Constraints with Type Ia Supernova in the LSST Era
Authors:
Ayan Mitra,
Richard Kessler,
Surhud More,
Renee Hlozek,
The LSST Dark Energy Science Collaboration
Abstract:
We perform a rigorous cosmology analysis on simulated type Ia supernovae (SN~Ia) and evaluate the improvement from including photometric host-galaxy redshifts compared to using only the "zspec" subset with spectroscopic redshifts from the host or SN. We use the Deep Drilling Fields (~50 deg^2) from the Photometric LSST Astronomical Time-Series Classification Challenge (PLaSTiCC), in combination wi…
▽ More
We perform a rigorous cosmology analysis on simulated type Ia supernovae (SN~Ia) and evaluate the improvement from including photometric host-galaxy redshifts compared to using only the "zspec" subset with spectroscopic redshifts from the host or SN. We use the Deep Drilling Fields (~50 deg^2) from the Photometric LSST Astronomical Time-Series Classification Challenge (PLaSTiCC), in combination with a low-z sample based on Data Challenge2 (DC2). The analysis includes light curve fitting to standardize the SN brightness, a high-statistics simulation to obtain a bias-corrected Hubble diagram, a statistical+systematics covariance matrix including calibration and photo-z uncertainties, and cosmology fitting with a prior from the cosmic microwave background. Compared to using the zspec subset, including events with SN+host photo-z results in i) more precise distances for z>0.5, ii) a Hubble diagram that extends 0.3 further in redshift, and iii) a 50 % increase in the Dark Energy Task Force figure of merit (FoM) based on the w0-wa CDM model. Analyzing 25 simulated data samples, the average bias on w0 and wa is consistent with zero. The host photo-z systematic of 0.01 reduces FoM by only 2 % because i) most z<0.5 events are in the zspec subset, ii) the combined SN+host photo-z has X 2 smaller bias, and iii) the anti-correlation between fitted redshift and color self corrects distance errors. To prepare for analysing real data, the next SNIa-cosmology analysis with photo-z's should include non SN-Ia contamination and host galaxy mis-associations.
△ Less
Submitted 16 October, 2022; v1 submitted 14 October, 2022;
originally announced October 2022.
-
Large Interferometer For Exoplanets (LIFE): VI. Detecting rocky exoplanets in the habitable zones of Sun-like stars
Authors:
Jens Kammerer,
Sascha P. Quanz,
Felix Dannert,
the LIFE Collaboration
Abstract:
While previous studies have shown a strong preference for a future mid-infrared nulling interferometer space mission to detect planets within the HZ around M dwarfs, we here focus on a more conservative approach toward the concept of habitability and present yield estimates for two stellar samples consisting of nearby (d<20 pc) Sun-like stars (4800-6300 K) and nearby FGK-type stars (3940-7220 K) a…
▽ More
While previous studies have shown a strong preference for a future mid-infrared nulling interferometer space mission to detect planets within the HZ around M dwarfs, we here focus on a more conservative approach toward the concept of habitability and present yield estimates for two stellar samples consisting of nearby (d<20 pc) Sun-like stars (4800-6300 K) and nearby FGK-type stars (3940-7220 K) accessible to such a mission. Our yield estimates are based on recently derived occurrence rates of rocky planets from the Kepler mission and our LIFE exoplanet observation simulation tool LIFEsim, which includes all main astrophysical noise sources, but no instrumental noise sources as yet. Depending on a pessimistic or optimistic extrapolation of the Kepler results, we find that during a 2.5-year search phase, LIFE could detect between ~10-16 (average) or ~5-34 (including 1$σ$ uncertainties) rocky planets (0.5-1.5 R${}_\rm{Earth}$) within the optimistic HZ of Sun-like stars and between ~4-6 (average) or ~1-13 (including 1$σ$ uncertainties) exo-Earth candidates (EECs) assuming four collector spacecraft equipped with 2 m mirrors and a conservative instrument throughput of 5%. With D=3.5 m or 1 m mirrors, the yield $Y$ changes strongly, following approximately $Y \propto D^{3/2}$. With the larger sample of FGK-type stars, the yield increases to ~16-22 (average) rocky planets within the optimistic HZ and ~5-8 (average) EECs. Furthermore, we find that in addition to the mirror diameter, the yield depends strongly on the total throughput, but only weakly on the exozodiacal dust level and the accessible wavelength range of the mission. When the focus lies entirely on Sun-like stars, larger mirrors (~3 m with 5% total throughput) or a better total throughput (~20% with 2 m mirrors) are required to detect a statistically relevant sample of ~30 rocky planets within the optimistic HZ.
△ Less
Submitted 19 October, 2022; v1 submitted 4 October, 2022;
originally announced October 2022.
-
Atmospheric retrievals for LIFE and other future space missions: the importance of mitigating systematic effects
Authors:
Eleonora Alei,
Björn S. Konrad,
Paul Mollière,
Sascha P. Quanz,
Daniel Angerhausen,
Mohanakrishna Ranganathan,
the LIFE collaboration
Abstract:
Atmospheric retrieval studies are essential to determine the science requirements for future generation missions, such as the Large Interferometer for Exoplanets (LIFE). The use of heterogeneous absorption cross-sections might be the cause of systematic effects in retrievals, which could bias a correct characterization of the atmosphere. In this contribution we quantified the impact of differences…
▽ More
Atmospheric retrieval studies are essential to determine the science requirements for future generation missions, such as the Large Interferometer for Exoplanets (LIFE). The use of heterogeneous absorption cross-sections might be the cause of systematic effects in retrievals, which could bias a correct characterization of the atmosphere. In this contribution we quantified the impact of differences in line list provenance, broadening coefficients, and line wing cut-offs in the retrieval of an Earth twin exoplanet orbiting a Sun-like star at 10 pc from the observer, as it would be observed with LIFE. We ran four different retrievals on the same input spectrum, by varying the opacity tables that the Bayesian retrieval framework was allowed to use. We found that the systematics introduced by the opacity tables could bias the correct estimation of the atmospheric pressure at the surface level, as well as an accurate retrieval of the abundance of some species in the atmosphere (such as CO$_2$ and N$_2$O). We argue that differences in the line wing cut-off might be the major source of errors. We highlight the need for more laboratory and modeling efforts, as well as inter-model comparisons of the main radiative transfer models and Bayesian retrieval frameworks. This is especially relevant in the context of LIFE and future generation missions, to identify issues and critical points for the community to jointly work together to prepare for the analysis of the upcoming observations.
△ Less
Submitted 30 September, 2022;
originally announced September 2022.
-
PSFs of coadded images
Authors:
Rachel Mandelbaum,
Mike Jarvis,
Robert H. Lupton,
James Bosch,
Arun Kannawadi,
Michael D. Murphy,
Tianqing Zhang,
the LSST Dark Energy Science Collaboration
Abstract:
We provide a detailed exploration of the connection between choice of coaddition schemes and the point-spread function (PSF) of the resulting coadded images. In particular, we investigate what properties of the coaddition algorithm lead to the final coadded image having a well-defined PSF. The key elements of this discussion are as follows:
1. We provide an illustration of how linear coaddition…
▽ More
We provide a detailed exploration of the connection between choice of coaddition schemes and the point-spread function (PSF) of the resulting coadded images. In particular, we investigate what properties of the coaddition algorithm lead to the final coadded image having a well-defined PSF. The key elements of this discussion are as follows:
1. We provide an illustration of how linear coaddition schemes can produce a coadd that lacks a well-defined PSF even for relatively simple scenarios and choices of weight functions.
2. We provide a more formal demonstration of the fact that a linear coadd only has a well-defined PSF in the case that either (a) each input image has the same PSF or (b) the coadd is produced with weights that are independent of the signal.
3. We discuss some reasons that two plausible nonlinear coaddition algorithms (median and clipped-mean) fail to produce a consistent PSF profile for stars.
4. We demonstrate that all nonlinear coaddition procedures fail to produce a well-defined PSF for extended objects.
In the end, we conclude that, for any purpose where a well-defined PSF is desired, one should use a linear coaddition scheme with weights that do not correlate with the signal and are approximately uniform across typical objects of interest.
△ Less
Submitted 2 February, 2023; v1 submitted 19 September, 2022;
originally announced September 2022.
-
A Joint Roman Space Telescope and Rubin Observatory Synthetic Wide-Field Imaging Survey
Authors:
M. A. Troxel,
C. Lin,
A. Park,
C. Hirata,
R. Mandelbaum,
M. Jarvis,
A. Choi,
J. Givans,
M. Higgins,
B. Sanchez,
M. Yamamoto,
H. Awan,
J. Chiang,
O. Dore,
C. W. Walter,
T. Zhang,
J. Cohen-Tanugi,
E. Gawiser,
A. Hearin,
K. Heitmann,
M. Ishak,
E. Kovacs,
Y. -Y. Mao,
M. Wood-Vasey,
the LSST Dark Energy Science Collaboration
Abstract:
We present and validate 20 deg$^2$ of overlapping synthetic imaging surveys representing the full depth of the Nancy Grace Roman Space Telescope High-Latitude Imaging Survey (HLIS) and five years of observations of the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST). The two synthetic surveys are summarized, with reference to the existing 300 deg$^2$ of LSST simulated imaging prod…
▽ More
We present and validate 20 deg$^2$ of overlapping synthetic imaging surveys representing the full depth of the Nancy Grace Roman Space Telescope High-Latitude Imaging Survey (HLIS) and five years of observations of the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST). The two synthetic surveys are summarized, with reference to the existing 300 deg$^2$ of LSST simulated imaging produced as part of Dark Energy Science Collaboration (DESC) Data Challenge 2 (DC2). Both synthetic surveys observe the same simulated DESC DC2 universe. For the synthetic Roman survey, we simulate for the first time fully chromatic images along with the detailed physics of the Sensor Chip Assemblies derived from lab measurements using the flight detectors. The simulated imaging and resulting pixel-level measurements of photometric properties of objects span a wavelength range of $\sim$0.3 to 2.0 $μ$m. We also describe updates to the Roman simulation pipeline, changes in how astrophysical objects are simulated relative to the original DC2 simulations, and the resulting simulated Roman data products. We use these simulations to explore the relative fraction of unrecognized blends in LSST images, finding that 20-30% of objects identified in LSST images with $i$-band magnitudes brighter than 25 can be identified as multiple objects in Roman images. These simulations provide a unique testing ground for the development and validation of joint pixel-level analysis techniques of ground- and space-based imaging data sets in the second half of the 2020s -- in particular the case of joint Roman--LSST analyses.
△ Less
Submitted 14 September, 2022;
originally announced September 2022.
-
A Unified Catalog-level Reanalysis of Stage-III Cosmic Shear Surveys
Authors:
Emily P. Longley,
Chihway Chang,
Christopher W. Walter,
Joe Zuntz,
Mustapha Ishak,
Rachel Mandelbaum,
Hironao Miyatake,
Andrina Nicola,
Eske M. Pedersen,
Maria E. S. Pereira,
Judit Prat,
J. Sánchez,
Tilman Tröster,
Michael Troxel,
Angus Wright,
The LSST Dark Energy Science Collaboration
Abstract:
Cosmological parameter constraints from recent galaxy imaging surveys are reaching $2-3\%$-level accuracy. The upcoming Legacy Survey of Space and Time (LSST) of the Vera C. Rubin Observatory will produce sub-percent level measurements of cosmological parameters, providing a milestone test of the $Λ$CDM model. To supply guidance to the upcoming LSST analysis, it is important to understand thorough…
▽ More
Cosmological parameter constraints from recent galaxy imaging surveys are reaching $2-3\%$-level accuracy. The upcoming Legacy Survey of Space and Time (LSST) of the Vera C. Rubin Observatory will produce sub-percent level measurements of cosmological parameters, providing a milestone test of the $Λ$CDM model. To supply guidance to the upcoming LSST analysis, it is important to understand thoroughly the results from different recent galaxy imaging surveys and assess their consistencies. In this work we perform a unified catalog-level reanalysis of three cosmic shear datasets: the first year data from the Dark Energy Survey (DES-Y1), the 1,000 deg$^{2}$ dataset from the Kilo-Degree Survey (KiDS-1000), and the first year data from the Hyper Suprime-Cam Subaru Strategic Program (HSC-Y1). We utilize a pipeline developed and rigorously tested by the LSST Dark Energy Science Collaboration to perform the reanalysis and assess the robustness of the results to analysis choices. We find the $S_{8}$ constraint to be robust to two different small-scale modeling approaches, and varying choices of cosmological priors. Our unified analysis allows the consistency of the surveys to be rigorously tested and we find the three surveys to be statistically consistent. Due to the partially overlapping footprint, we model the cross-covariance between KiDS-1000 and HSC-Y1 approximately when combining all three datasets, resulting in a $1.6-1.9\%$ constraint on $S_8$ given different assumptions on the cross-covariance.
△ Less
Submitted 15 August, 2022;
originally announced August 2022.
-
Search for multi-messenger events during LIGO/Virgo era
Authors:
Tosta e Melo,
I.,
LVK collaboration
Abstract:
Multi-messenger astronomy is a vast and expanding field as electromagnetic observations (EM) are no longer the only way of exploring the Universe. Due to the new messengers, astrophysical events with both gravitational waves (GWs) and EM emission are no longer a dream of the astronomical community. A breakthrough for GW multi-messenger astronomy came when the LIGO-Virgo network detected a GW signa…
▽ More
Multi-messenger astronomy is a vast and expanding field as electromagnetic observations (EM) are no longer the only way of exploring the Universe. Due to the new messengers, astrophysical events with both gravitational waves (GWs) and EM emission are no longer a dream of the astronomical community. A breakthrough for GW multi-messenger astronomy came when the LIGO-Virgo network detected a GW signal of two low-mass compact objects consistent with a binary neutron star (BNS, GW170817) an event that generated a short gamma-ray burst (GRBs) and a kilonova. While GW170817 represents the testimony to BNS mergers being the progenitor of at least some GRBs, a wide range of highly energetic astrophysical phenomena is expected to be accompanied by the emission of GWs and photons. Here we present an unmodelled method to search for GWs having gamma and radio counterparts, using the LIGO-Virgo data and observations of partner telescopes. We also discuss the most recent results of the unmodelled coherent searches targeting astrophysical events during the first part of the LIGO-Virgo third observing run (O3a): 105 GRBs detected by the Fermi and Swift satellites.
△ Less
Submitted 28 June, 2022;
originally announced June 2022.
-
The Simulated Catalogue of Optical Transients and Correlated Hosts (SCOTCH)
Authors:
Martine Lokken,
Alexander Gagliano,
Gautham Narayan,
Renée Hložek,
Richard Kessler,
John Franklin Crenshaw,
Laura Salo,
Catarina S. Alves,
Deep Chatterjee,
Maria Vincenzi,
Alex I. Malz,
The LSST Dark Energy Science Collaboration
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-dr…
▽ More
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 modeling 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 ($\lesssim$2 day) cadence. In conjunction with SCOTCH, we also release an associated set of tutorials and the 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.
△ Less
Submitted 27 February, 2023; v1 submitted 6 June, 2022;
originally announced June 2022.
-
Optimising the shape of photometric redshift distributions with clustering cross-correlations
Authors:
Benjamin Stölzner,
Benjamin Joachimi,
Andreas Korn,
the LSST Dark Energy Science Collaboration
Abstract:
We present an optimisation method for the assignment of photometric galaxies into a chosen set of redshift bins. This is achieved by combining simulated annealing, an optimisation algorithm inspired by solid-state physics, with an unsupervised machine learning method, a self-organising map (SOM) of the observed colours of galaxies. Starting with a sample of galaxies that is divided into redshift b…
▽ More
We present an optimisation method for the assignment of photometric galaxies into a chosen set of redshift bins. This is achieved by combining simulated annealing, an optimisation algorithm inspired by solid-state physics, with an unsupervised machine learning method, a self-organising map (SOM) of the observed colours of galaxies. Starting with a sample of galaxies that is divided into redshift bins based on a photometric redshift point estimate, the simulated annealing algorithm repeatedly reassigns SOM-selected subsamples of galaxies, which are close in colour, to alternative redshift bins. We optimise the clustering cross-correlation signal between photometric galaxies and a reference sample of galaxies with well-calibrated redshifts. Depending on the effect on the clustering signal, the reassignment is either accepted or rejected. By dynamically increasing the resolution of the SOM, the algorithm eventually converges to a solution that minimises the number of mismatched galaxies in each tomographic redshift bin and thus improves the compactness of their corresponding redshift distribution. This method is demonstrated on the synthetic LSST cosmoDC2 catalogue. We find a significant decrease in the fraction of catastrophic outliers in the redshift distribution in all tomographic bins, most notably in the highest redshift bin with a decrease in the outlier fraction from 57 per cent to 16 per cent.
△ Less
Submitted 7 February, 2023; v1 submitted 26 May, 2022;
originally announced May 2022.