Engineering topics

Mandelbaum, Rachel (ORCID:0000000322711527)

Publications and source records attributed to Mandelbaum, Rachel (ORCID:0000000322711527).

15 records

DOE OSTI2025

Accurate shear estimation with fourth-order moments

ABSTRACT As imaging surveys progress in exploring the large-scale structure of the Universe through the use of weak gravitational lensing, achieving sub-per cent accuracy in estimating shape distortions caused by lensing, or shear, is imperative for precision cosmology. In this paper, we extend the Fourier power function shapelets (FPFS) shear estimator using fourth-order shapelet moments and combine it with the original second-order shear estimator to reduce galaxy shape noise. We calibrate this novel shear estimator analytically to a sub-per cent level-accuracy using the AnaCal framework. This higher order shear estimator is tested with realistic image simulations, and after analytical correction for the detection/selection bias and noise bias, the multiplicative shear bias $|m|$ is below $3\times 10^{-3}$ (99.7 per cent confidence interval) for both isolated and blended galaxies. Once combined with the second-order FPFS shear estimator, the shape noise is reduced by $\sim 35~{{\ \rm per\ cent}}$ for isolated galaxies in simulations with Hyper Suprime-Cam and Vera C. Rubin Observatory Legacy Survey of Space and Time observational conditions. However, for blended galaxies, the effective number density does not significantly improve with the combination of the two estimators. Based on these results, we recommend exploration of how this framework can further reduce the systematic uncertainties in shear due to point spread function leakage and modelling error, and potentially provide improved precision in shear inference in high-resolution space-based images.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Photometric redshifts and intrinsic alignments: Degeneracies and biases in the 3 × 2 pt analysis

We present a systematic study of cosmological parameter bias in weak lensing and large-scale structure analyses for upcoming imaging surveys induced by the interplay of intrinsic alignments (IA) and photometric redshift (photo-z) model misspecification error. We first examine the degeneracies between the parameters of the tidal alignment–tidal torquing (TATT) model for IA and of a photo-z model including a mean shift ( Δ z ¯ ) and variance ( σ z ) for each tomographic bin of lenses and sources, under a variety of underlying true IA behaviors. We identify strong degeneracies between: (1) the redshift scaling of the tidal alignment amplitude and the mean shift and variances of source bins, (2) the redshift scaling of the tidal torquing amplitude and the variance of the lowest- z source bin, and (3) the IA source density weighting and the mean shift and variance of several source bins. We then use this information to guide our exploration of the level of cosmological parameter bias which can be induced given incorrect modeling of IA, photo-z, or both. We find that marginalizing over all the parameters of TATT is generally sufficient to preclude cosmological parameter bias in the scenarios we consider. However, this does not necessarily mean that IA and photo-z parameters are themselves unbiased, nor does it mean that the best-fit model is a good fit to the data. We also find scenarios where the inferred parameters produce χ d . o . f . 2 values indicative of a good fit but cosmological parameter bias is significant, particularly when the IA source density weighting parameter is not marginalized over. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2024

Increasing the power of weak lensing data with multipole-based intrinsic alignment estimators

ABSTRACT It has long been known that galaxy shapes align coherently with the large-scale density field. Characterizing this effect is essential to interpreting measurements of weak gravitational lensing, the deflection of light from distant galaxies by matter overdensities along the line of sight, as it also produces coherent galaxy alignments that we wish to interpret in terms of a cosmological model. Existing direct measurements of intrinsic alignments using galaxy samples with high-quality shape and redshift measurements typically use well-understood but sub-optimal projected estimators, which do not make good use of the information in the data when comparing those estimators to theoretical models. We demonstrate a more optimal estimator, based on a multipole expansion of the correlation functions or power spectra, for direct measurements of galaxy intrinsic alignments. We show that even using the lowest order multipole alone increases the significance of inferred model parameters using simulated and real data, without any additional modelling complexity. We apply this estimator to measurements of parameters of the non-linear alignment model using data from the Sloan Digital Sky survey, demonstrating consistent results with a factor of ∼2 greater precision in parameter fits to intrinsic alignments models. This result is functionally equivalent to quadrupling the survey area, but without the attendant costs – thereby demonstrating the value in using this new estimator in current and future intrinsic alignments measurements using spectroscopic galaxy samples.

Astronomy & Astrophysics↗

DOE OSTI2023

Hyper Suprime-Cam Year 3 results: Cosmology from galaxy clustering and weak lensing with HSC and SDSS using the emulator based halo model

We present cosmology results from a blinded joint analysis of cosmic shear, galaxy-galaxy weak lensing, and projected galaxy clustering, measured from the Hyper Suprime-Cam three-year (HSC-Y3) shape catalog and the Sloan Digital Sky Survey (SDSS) DR11 spectroscopic galaxy catalog—a cosmology analysis. We define luminosity-cut, and therefore nearly volume-limited, samples of SDSS galaxies to serve as the tracers of w p and as the lens samples for Δ⁢Σ in three spectroscopic redshift bins spanning the range 0.15± and Δ⁢Σ measurements, we use a single sample of about seven million source galaxies over 416 deg 2 , selected from HSC-Y3 based on having photometric redshifts (photo-z) greater than 0.75. The deep, high-quality HSC-Y3 data enable significant detections of the Δ⁢Σ signals, with integrated signal-to-noise ratio S/N~24 in the range 3≤R/[ℎ –1 Mpc]≤30 over the three lens samples. ξ ± has S/N ~19 in the range 8' ≤$\theta$ ≤5⁢0' and 30' ≤$\theta$ ≤15⁢0' for ξ + and ξ – , respectively. For cosmological parameter inference, we use the dark emulator package, combined with a halo occupation distribution prescription for the relation between galaxies and halos, to model w p and Δ⁢Σ down to quasinonlinear scales, and we estimate cosmological parameters after marginalizing over nuisance parameters. In our baseline analysis we employ an uninformative flat prior of the residual photo-z error, given by Π⁡(Δz ph )=$U$⁡(–1,1), to model a residual bias in the mean redshift of HSC source galaxies. Comparing the relative lensing amplitudes for Δ⁢Σ in the three redshift bins and for ξ ± with the single HSC source galaxy sample allows us to calibrate the photo-z parameter Δ⁢z ph to the precision of σ⁡(Δ⁢z ph )≃0.09. With these methods, we obtain a robust constraint on the cosmological parameters for the flat Λ⁢CDM model: S 8 =σ 8 ⁢(Ω m /0.3) 0.5 =$0.763^{+0.040}_{–0.036}$, or the best-constrained parameter given by $S^{'}_{8}$ = σ 8 ⁢(Ω m /0.3) 0.22 =0.721±0.028, determined with about 4% fractional precision. Based on multidimensional tension metrics, HSC-Y3 data exhibits about 2.5σ tension with the cosmological constraint inferred by Planck for the Λ⁢CDM model, and hints at a nonzero residual photo-z bias implying that the true mean redshift of the HSC galaxies at z ≳0.75 is higher than that implied by the original photo-z estimates.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

Hyper Suprime-Cam Year 3 results: Measurements of clustering of SDSS-BOSS galaxies, galaxy-galaxy lensing, and cosmic shear

Here we utilize the Sloan Digital Sky Survey Baryon Oscillation Spectroscopic Survey (SDSS-BOSS) galaxies and its overlap with approximately 416 sq degrees of deep grizy-band imaging from the Subaru Hyper Suprime-Cam Survey (HSC). We perform measurements of three two-point correlations which form the basis of the cosmological inference presented in our companion papers, Miyatake et al. and Sugiyama et al. We use three approximately volume limited subsamples of spectroscopic galaxies by their i-band magnitude from the SDSS-BOSS: LOWZ (0.1 < z < 0.35), CMASS1 (0.43 < z < 0.55) and CMASS2 (0.55 < z < 0.7), respectively. We present high signal-to-noise ratio measurements of the projected correlation functions of these galaxies, which is expected to be proportional to the projected matter correlation function on large scales with a proportionality constant dependent on the bias of galaxies. In order to help break the degeneracy between the amplitude of the matter correlation and the bias of these spectroscopic galaxies, we use the distortions of the shapes of fainter galaxies in HSC due to weak gravitational lensing, to measure the galaxy-galaxy lensing signal, which probes the projected galaxy-matter cross-correlation function of the SDSS-BOSS galaxies. We also measure the cosmic shear correlation functions from HSC galaxies which is related to the projected matter correlation function. We demonstrate the robustness of our measurements by subjecting each of them to a variety of systematic tests. Our use of a single sample of HSC source galaxies is crucial to calibrate any residual systematic biases in the inferred redshifts of our galaxies. We also describe the construction of a suite of mocks: (i) spectroscopic galaxy catalogs which obey the clustering and abundance of each of the three SDSS-BOSS subsamples, and (ii) galaxy shape catalogs which obey the footprint of the HSC survey and have been appropriately sheared by the large-scale structure expected in a Lambda Cold Dark Matter model. Finally, we use these mock catalogs to compute the covariance of each of our observables.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

Hyper Suprime-Cam Year 3 results: Cosmology from galaxy clustering and weak lensing with HSC and SDSS using the minimal bias model

Here, we present cosmological parameter constraints from a blind joint analysis of three two-point correlation functions measured from the Year 3 Hyper Suprime-Cam (HSC-Y3) imaging data, covering about 416 deg 2 , and the SDSS DR11 spectroscopic galaxies spanning the redshift range [0.15, 0.70]. We subdivide the SDSS galaxies into three luminosity-cut, and therefore nearly volume-limited samples separated in redshift, each of which acts as a large-scale structure tracer characterized by the measurement of the projected correlation function, w p ⁡(R). We also use the measurements of the galaxy-galaxy weak-lensing signal Δ⁢Σ⁡(R) for each of these SDSS samples which act as lenses for a secure sample of source galaxies selected from the HSC-Y3 shape catalog based on their photometric redshifts. We combine these measurements with the cosmic shear correlation functions, ξ ± ⁡($\vartheta$) measured for our HSC source sample. We model these observables with the minimal bias model of the galaxy clustering observables in the context of a flat Λ⁢CDM cosmology. We use conservative scale cuts, R >12 and 8⁢ℎ -1 Mpc for Δ⁢Σ and w p , respectively, where the minimal bias model is valid, in addition to conservative prior on the residual bias in the mean redshift of the HSC photometric source galaxies. We present various validation tests of our model as well as analysis methods. Our baseline analysis yields S 8 = $0.77⁢5^{+0.043}_{-0.038}$ (68% C.I.) for the Λ⁢CDM model, after marginalizing over uncertainties in other parameters. Our value of S 8 is consistent with that from the Planck 2018 data, but the credible interval of our result is still relatively large. We show that various internal consistency tests based on different splits of the data are passed. Our results are statistically consistent with those of a companion paper, which extends this analysis to smaller scales with an emulator-based halo model, using Δ⁢Σ⁡(R) and w p ⁡(R) down to R >3 and 2⁢ℎ -1 Mpc, respectively.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

A general framework for removing point-spread function additive systematics in cosmological weak lensing analysis

ABSTRACT Cosmological weak lensing measurements rely on a precise measurement of the shear two-point correlation function (2PCF) along with a deep understanding of systematics that affect it. In this work, we demonstrate a general framework for detecting and modelling the impact of PSF systematics on the cosmic shear 2PCF and mitigating its impact on cosmological analysis. Our framework can detect PSF leakage and modelling error from all spin-2 quantities contributed by the PSF second and higher moments, rather than just the second moments, using the cross-correlations between galaxy shapes and PSF moments. We interpret null tests using the HSC Year 3 (Y3) catalogs with this formalism and find that leakage from the spin-2 combination of PSF fourth moments is the leading contributor to additive shear systematics, with total contamination that is an order-of-magnitude higher than that contributed by PSF second moments alone. We conducted a mock cosmic shear analysis for HSC Y3 and find that, if uncorrected, PSF systematics can bias the cosmological parameters Ωm and S8 by ∼0.3σ. The traditional second moment-based model can only correct for a 0.1σ bias, leaving the contamination largely uncorrected. We conclude it is necessary to model both PSF second and fourth moment contaminations for HSC Y3 cosmic shear analysis. We also reanalyse the HSC Y1 cosmic shear analysis with our updated systematics model and identify a 0.07σ bias on Ωm when using the more restricted second moment model from the original analysis. We demonstrate how to self-consistently use the method in both real space and Fourier space, assess shear systematics in tomographic bins, and test for PSF model overfitting.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2023

A unified catalogue-level reanalysis of stage-III cosmic shear surveys

ABSTRACT Cosmological parameter constraints from recent galaxy imaging surveys are reaching percent-level accuracy on the effective amplitude of the lensing signal, S8. The upcoming Legacy Survey of Space and Time (LSST) of the Vera C. Rubin Observatory will produce subpercent 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 catalogue-level reanalysis of three cosmic shear data sets: the first year data from the Dark Energy Survey (DES-Y1), the 1000 deg2 data set 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 S8 constraint to be robust to two different small-scale modelling 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 data sets, resulting in a 1.6–1.9 per cent constraint on S8 given different assumptions on the cross-covariance.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Impact of point spread function higher moments error on weak gravitational lensing – II. A comprehensive study

ABSTRACT Weak lensing is one of the most powerful probes for dark matter and dark energy science, although it faces increasing challenges in controlling systematic uncertainties as the statistical errors become smaller. The point spread function (PSF) needs to be precisely modeled to avoid systematic error on the weak lensing measurements. The weak lensing biases induced by errors in the PSF model second moments, i.e. its size and shape, are well-studied. However, Zhang et al. showed that errors in the higher moments of the PSF may also be a significant source of systematics for upcoming weak lensing surveys. Therefore, this work comprehensively investigate the modelling quality of PSF moments from the 3rd to 6th order, and propagate the PSFEx higher moments modelling error in the HSC survey data set to the weak lensing shear–shear correlation functions and their cosmological analyses. The overall multiplicative shear bias associated with errors in PSF higher moments can cause an ∼0.1σ shift on the cosmological parameters for LSST Y10, while the associated additive biases can induce 1σ uncertainties in cosmology parameter inference for LSST Y10, if not accounted. We compare the PSFEx model with PSF in Full FOV, and find similar performance in modelling the PSF higher moments. We conclude that PSF higher moment errors of the future PSF models should be reduced from those in current methods, otherwise needed to be explicitly modeled in the weak lensing analysis.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Photometric redshift uncertainties in weak gravitational lensing shear analysis: models and marginalization

ABSTRACT Recovering credible cosmological parameter constraints in a weak lensing shear analysis requires an accurate model that can be used to marginalize over nuisance parameters describing potential sources of systematic uncertainty, such as the uncertainties on the sample redshift distribution n(z). Due to the challenge of running Markov chain Monte Carlo (MCMC) in the high-dimensional parameter spaces in which the n(z) uncertainties may be parametrized, it is common practice to simplify the n(z) parametrization or combine MCMC chains that each have a fixed n(z) resampled from the n(z) uncertainties. In this work, we propose a statistically principled Bayesian resampling approach for marginalizing over the n(z) uncertainty using multiple MCMC chains. We self-consistently compare the new method to existing ones from the literature in the context of a forecasted cosmic shear analysis for the HSC three-year shape catalogue, and find that these methods recover statistically consistent error bars for the cosmological parameter constraints for predicted HSC three-year analysis, implying that using the most computationally efficient of the approaches is appropriate. However, we find that for data sets with the constraining power of the full HSC survey data set (and, by implication, those upcoming surveys with even tighter constraints), the choice of method for marginalizing over n(z) uncertainty among the several methods from the literature may modify the 1σ uncertainties on Ωm–S8 constraints by ∼4 per cent, and a careful model selection is needed to ensure credible parameter intervals.

Zhang, Tianqing (ORCID:000000025596198X)↗

DOE OSTI2022

The Impact of Observing Strategy on Cosmological Constraints with LSST

The generation-defining Vera C. Rubin Observatory will make state-of-the-art measurements of both the static and transient universe through its Legacy Survey for Space and Time (LSST). With such capabilities, it is immensely challenging to optimize the LSST observing strategy across the survey's wide range of science drivers. Many aspects of the LSST observing strategy relevant to the LSST Dark Energy Science Collaboration, such as survey footprint definition, single-visit exposure time, and the cadence of repeat visits in different filters, are yet to be finalized. Here, we present metrics used to assess the impact of observing strategy on the cosmological probes considered most sensitive to survey design; these are large-scale structure, weak lensing, type Ia supernovae, kilonovae, and strong lens systems (as well as photometric redshifts, which enable many of these probes). We evaluate these metrics for over 100 different simulated potential survey designs. Our results show that multiple observing strategy decisions can profoundly impact cosmological constraints with LSST; these include adjusting the survey footprint, ensuring repeat nightly visits are taken in different filters, and enforcing regular cadence. We provide public code for our metrics, which makes them readily available for evaluating further modifications to the survey design. We conclude with a set of recommendations and highlight observing strategy factors that require further research.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Optimization of the Observing Cadence for the Rubin Observatory Legacy Survey of Space and Time: A Pioneering Process of Community-focused Experimental Design

Vera C. Rubin Observatory is a ground-based astronomical facility under construction, a joint project of the National Science Foundation and the U.S. Department of Energy, designed to conduct a multipurpose 10 yr optical survey of the Southern Hemisphere sky: the Legacy Survey of Space and Time. Significant flexibility in survey strategy remains within the constraints imposed by the core science goals of probing dark energy and dark matter, cataloging the solar system, exploring the transient optical sky, and mapping the Milky Way. The survey’s massive data throughput will be transformational for many other astrophysics domains and Rubin’s data access policy sets the stage for a huge community of potential users. To ensure that the survey science potential is maximized while serving as broad a community as possible, Rubin Observatory has involved the scientific community at large in the process of setting and refining the details of the observing strategy. The motivation, history, and decision-making process of this strategy optimization are detailed in this paper, giving context to the science-driven proposals and recommendations for the survey strategy included in this Focus Issue.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2021

Impact of point spread function higher moments error on weak gravitational lensing

ABSTRACT Weak gravitational lensing is one of the most powerful tools for cosmology, while subject to challenges in quantifying subtle systematic biases. The point spread function (PSF) can cause biases in weak lensing shear inference when the PSF model does not match the true PSF that is convolved with the galaxy light profile. Although the effect of PSF size and shape errors – i.e. errors in second moments – is well studied, weak lensing systematics associated with errors in higher moments of the PSF model require further investigation. The goal of our study is to estimate their potential impact for LSST weak lensing analysis. We go beyond second moments of the PSF by using image simulations to relate multiplicative bias in shear to errors in the higher moments of the PSF model. We find that the current level of errors in higher moments of the PSF model in data from the Hyper Suprime-Cam survey can induce a ∼0.05 per cent shear bias, making this effect unimportant for ongoing surveys but relevant at the precision of upcoming surveys such as LSST.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

The LSST DESC DC2 Simulated Sky Survey

Here, we describe the simulated sky survey underlying the second data challenge (DC2) carried out in preparation for analysis of the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) by the LSST Dark Energy Science Collaboration (LSST DESC). Significant connections across multiple science domains will be a hallmark of LSST; the DC2 program represents a unique modeling effort that stresses this interconnectivity in a way that has not been attempted before. This effort encompasses a full end-to-end approach: starting from a large N-body simulation, through setting up LSST-like observations including realistic cadences, through image simulations, and finally processing with Rubin's LSST Science Pipelines. This last step ensures that we generate data products resembling those to be delivered by the Rubin Observatory as closely as is currently possible. The simulated DC2 sky survey covers six optical bands in a wide-fast-deep area of approximately 300 deg 2 , as well as a deep drilling field of approximately 1 deg 2 . We simulate 5 yr of the planned 10 yr survey. The DC2 sky survey has multiple purposes. First, the LSST DESC working groups can use the data set to develop a range of DESC analysis pipelines to prepare for the advent of actual data. Second, it serves as a realistic test bed for the image processing software under development for LSST by the Rubin Observatory. In particular, simulated data provide a controlled way to investigate certain image-level systematic effects. Finally, the DC2 sky survey enables the exploration of new scientific ideas in both static and time domain cosmology.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2019

CosmoDC2: A Synthetic Sky Catalog for Dark Energy Science with LSST

This paper introduces cosmoDC2, a large synthetic galaxy catalog designed to support precision dark energy science with the Large Synoptic Survey Telescope (LSST). CosmoDC2 is the starting point for the second data challenge (DC2) carried out by the LSST Dark Energy Science Collaboration (LSST DESC). The catalog is based on a trillion-particle, (4.225 Gpc) 3 box cosmological N-body simulation, the Outer Rim run. It covers 440 deg 2 of sky area to a redshift of z = 3 and matches expected number densities from contemporary surveys to a magnitude depth of 28 in the r band. Each galaxy is characterized by a multitude of galaxy properties including stellar mass, morphology, spectral energy distributions, broadband filter magnitudes, host halo information, and weak lensing shear. The size and complexity of cosmoDC2 requires an efficient catalog generation methodology; our approach is based on a new hybrid technique that combines data-based empirical approaches with semianalytic galaxy modeling. A wide range of observation-based validation tests has been implemented to ensure that cosmoDC2 enables the science goals of the planned LSST DESC DC2 analyses. This paper also represents the official release of the cosmoDC2 data set, including an efficient reader that facilitates interaction with the data.

79 ASTRONOMY AND ASTROPHYSICS↗