Cosmic Shear Research Articles

B-modes from galaxy cluster alignments in future surveys

Intrinsic alignment (IA) of source galaxies represents an important contaminant for upcoming cosmic shear surveys. In particular, it is expected on general grounds that IA contains a B-mode while the weak lensing signal does not. Thus, a detection of B-modes offers the possibility to study directly the IA signal of the sources. Galaxy clusters exhibit strong IA and are therefore a natural candidate to look for a B-mode signal. We forecast the signal-to-noise ratio (SNR) for B-modes from IA of galaxy clusters in the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST). We use a perturbative model for the IA multipoles based on the Effective Field Theory of Intrinsic Alignments (EFT of IA), which has recently been validated against N-body simulations. We forecast SNR ≈12 and find that this detectability is not significantly impacted by different analysis choices. Lastly, we also apply our forecast to clusters in the redMaPPer SDSS and DESY1 samples. We find SNR ≈5 and SNR ≈3, respectively, suggesting a detection is within reach, provided accurate redshift information is available.

Accuracy requirements on intrinsic alignments for Stage-IV cosmic shear

In the context of cosmological weak lensing studies, intrinsic alignments (IAs) are one of the most In the context of cosmological weak lensing studies, intrinsic alignments (IAs) are one of the most complicated astrophysical systematic to model, given the poor understanding of the physical processes that cause them. A number of modelling frameworks for IAs have been proposed in the literature, both purely phenomenological or grounded on a perturbative treatment of symmetry-based arguments. However, the accuracy with which any of these approaches is able to describe the impact of IAs on cosmic shear data, particularly on the comparatively small scales () to which this observable is sensitive, is not clear. Here we quantify the level of disagreement between the true underlying intrinsic alignments and the theoretical model used to describe them that can be allowed in the context of cosmic shear analyses with future Stage-IV surveys. We consider various models describing this “IA residual’’, covering both physics-based approaches, as well as completely agnostic prescriptions. The same qualitative results are recovered in all cases explored: for a Stage-IV cosmic shear survey, a mis-modelling of the IA contribution at the ∼10% level produces shifts of ≲0.5σ on the final cosmological parameter constraints. Current and future IA models should therefore aim to achieve this level of accuracy, a prospect that is not unfeasible for models with sufficient flexibility.

Constraining hot dark matter sub-species with weak lensing and the cosmic microwave background radiation

Although it is well known that the bulk of dark matter (DM) has to be cold, the existence of an additional sub-dominant, hot species remains a valid possibility. In this paper we investigate the potential of the cosmic shear power spectrum to constrain such a mixed (hot plus cold) DM scenario with two additional free parameters, the hot-to-total DM fraction ($f_ hdm $) and the thermal mass of the hot component ($m_ hdm $). Running a Bayesian inference analysis for both the Kilo-Degree Survey cosmic shear data ( as well as the cosmic microwave background (CMB) temperature and polarisation data from Planck we derive new constraints for the mixed DM scenario. We find a 95 per cent confidence limit of $f_ hdm <0.08$ for a very hot species of $m_ hdm eV. This constraint is weakened to hdm <0.25$ for $m_ hdm eV. Scenarios with masses above hdm eV remain unconstrained by the data. Next to providing limits, we investigate the potential of mixed DM to address the clustering (or $S_8$) tension between lensing and the CMB. We find a reduction of the 2D ($ - S_8$) tension from 2.9sigma to 1.6sigma when going from a pure cold DM to a mixed DM scenario. When computing the 1D Gaussian tension on $S_8$ the improvement is milder, from 2.4sigma to 2.0sigma

Generation of realistic input parameters for simulating atmospheric point-spread functions at astronomical observatories

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 (https://github.com/LSSTDESC/psf-weather-station)

Sensitivity of Cosmological Parameter Estimation to Nonlinear Prescription from Cosmic Shear

Several ongoing and upcoming large-scale structure surveys aim to explore the nonlinear regime of structure formation with high precision. Making reliable cosmological inferences from these observations necessitates precise theoretical modeling of the mildly nonlinear regime. In this work we explore how the choice of nonlinear prescription would impact parameter estimation from cosmic shear measurements for a Euclid-like survey. Specifically, we employ two different nonlinear prescriptions of halofit and the Effective Field Theory of the Large Scale Structure and compare their measurements for the three different cosmological scenarios of ΛCDM, wCDM, and (w 0, w a ) CDM. We also investigate the impact of different nonlinear cutoff schemes on parameter estimation. We find that the predicted errors on most parameters shrink considerably as smaller scales are included in the analysis, with the amount depending on the nonlinear prescription and the cutoff scheme used. We use predictions from the halofit model to analyze the mock data from DarkSky N-body simulations and quantify the parameter bias introduced in the measurements due to the choice of nonlinear prescription. We observe that σ 8 and n s have the largest measurement bias induced by inaccuracies of the halofit model.

Observational bounds on extended minimal theories of massive gravity: new limits on the graviton mass

In this work, we derive for the first time observational constraints on the extended Minimal Theory of Massive Gravity (eMTMG) framework in light of Planck-CMB data, geometrical measurements from Baryon Acoustic Oscillation (BAO), Type Ia supernovae from the recent Pantheon+ samples, and also using the auto and cross-correlations cosmic shear measurements from KIDS-1000 survey. Given the great freedom of dynamics choice for the theory, we consider an observationally motivated subclass in which the background evolution of the Universe goes through a transition from a (positive or negative) value of the effective cosmological constant to another value. From the statistical point of view, we did not find evidence of such a transition, i.e. deviation from the standard ΛCDM behavior, and from the joint analysis using Planck + BAO + Pantheon+ data, we constrain the graviton mass to < 6.6 × 10-34 eV at 95% CL. We use KIDS-1000 survey data to constrain the evolution of the scalar perturbations of the model and its limits for the growth of structure predicted by the eMTMG scenario. In this case, we find small evidence at 95% CL for a non-zero graviton mass. We interpret and discuss these results in light of the current tension on the S 8 parameter. We conclude that, within the subclass considered, the current data are only able to impose upper bounds on the eMTMG dynamics. Given its potentialities beyond the subclass, eMTMG can be classified as a good candidate for modified gravity, serving as a framework in which observational data can effectively constrain (or confirm) the graviton mass and deviations from the standard ΛCDM behavior.

Euclidpreparation

This work considers which higher order modeling effects on the cosmic shear angular power spectra must be taken into account forEuclid. We identified the relevant terms and quantified their individual and cumulative impact on the cosmological parameter inferences fromEuclid. We computed the values of these higher order effects using analytic expressions and calculated the impact on cosmological parameter estimations using the Fisher matrix formalism. We reviewed 24 effects and determined the ones that potentially need to be accounted for, namely: the reduced shear approximation, magnification bias, source-lens clustering, source obscuration, local Universe effects, and the flat Universe assumption. After computing these effects explicitly and calculating their cosmological parameter biases, using a maximum multipole ofℓ = 5000, we find that the magnification bias, source-lens clustering, source obscuration, and local Universe terms individually produce significant (&gt; 0.25σ) cosmological biases in one or more parameters; accordingly, these effects must be accounted for and warrant further investigation. In total, we find biases in Ωm, Ωb,h, andσ8of 0.73σ, 0.28σ, 0.25σ, and −0.79σ, respectively, for the flat ΛCDM. For thew0waCDM case, we found biases in Ωm, Ωb,h,ns,σ8, andwaof 1.49σ, 0.35σ, −1.36σ, 1.31σ, −0.84σ, and −0.35σ, respectively. These are increased relative to the ΛCDM due to additional degeneracies as a function of redshift and scale.

Dark dimension and decaying dark matter gravitons

We explore the cosmology of the dark-dimension scenario taking into account perturbations in the linear regime. In the context of the dark-dimension scenario, a natural candidate for dark matter in our Universe is the excitations of a tower of massive spin-2 Kaluza-Klein (KK) gravitons. These dark gravitons are produced in the early Universe and decay to lighter KK gravitons during the course of cosmological evolution. The decay causes the average dark matter mass to decrease as the Universe evolves. In addition, the kinetic energy liberated in each decay leads to a kick velocity for the dark matter particles, leading to a suppression of structure formation. Using current cosmic microwave background (), baryon acoustic oscillation, and cosmic shear (KiDS-1000) data, we put a bound on the dark matter kick velocity today vtoday≤2.2×10−4c at 95% CL. This leads to rather specific regions of parameter space for the dark-dimension scenario. The combination of the experimental bounds from cosmology, astrophysics, and table-top experiments leads to the range l5∼1–10 μm for the size of the dark dimension. The dark-dimension scenario is found to be remarkably consistent with current observations and provides signatures that are within reach of near-future experiments. Published by the American Physical Society 2024

Making the leap. Part I. Modelling the reconstructed lensing convergence PDF from cosmic shear with survey masks and systematics

The last few years have seen the development of a promising theoretical framework for statistics of the cosmic large-scale structure — the theory of large deviations (LDT) for modelling weak-lensing one-point statistics in the mildly nonlinear regime. The goal of this series of papers is to make the leap and lay out the steps to perform an actual data analysis with this theoretical tool. Building upon the LDT framework, in this work (paper I) we demonstrate how to accurately model the Probability Distribution Function (PDF) of a reconstructed Kaiser-Squires convergence field under a realistic mask, that of the third data release of the Dark Energy Survey (DES). We also present how weak lensing systematics and higher-order lensing corrections due to intrinsic alignments, shear biases, photo-z errors and baryonic feedback can be incorporated in the modelling of the reconstructed convergence PDF. In an upcoming work (paper II) we will then demonstrate the robustness of our modelling through simulated likelihood analyses, the final step required before applying our method to actual data.

Probing the two-body decaying dark matter scenario with weak lensing and the cosmic microwave background

Decaying dark matter (DDM) scenarios have recently regained attention due to their potential ability to resolve the well-known clustering (or S8) tension between weak lensing (WL) and cosmic microwave background (CMB) measurements. In this paper, we investigate a well-established model where the original dark matter particle decays into a massless particle and a massive daughter particle. The latter obtains a velocity kick during the decay process that results in the suppression of the matter power spectrum at scales that are observable with WL shear observations. We perform the first fully non-linear WL analysis of this two-body decaying dark matter (ΛDDM) scenario, including intrinsic alignment and baryonic feedback processes. We used the cosmic shear band power spectra from KiDS-1000 data and combined it with temperature and polarisation data from Planck in order to constrain the ΛDDM model. We report new limits on the decay rate and mass splitting parameters that are significantly stronger than previous results, especially in the case of low-mass splittings. Regarding the S8 tension, we found a reduction from about 3 to 2σ, depending on which statistical measure is applied. We therefore conclude that the two-body ΛDDM model is able to reduce the S8 tension without convincingly solving it.

Constraining ultralight axions with CSST weak gravitational lensing and galaxy clustering photometric surveys

ABSTRACT Ultralight axion (ULA) can be one of the potential candidates for dark matter. The extremely low mass of the ULA can lead to a de Broglie wavelength the size of galaxies which results in a suppression of the growth of structure on small scales. In this work, we forecast the constraint on the ULA particle mass ma and relative fraction to dark matter fa = Ωa/Ωd for the forthcoming Stage IV space-based optical survey equipment CSST (China Space Station Telescope). We focus on the CSST cosmic shear and galaxy clustering photometric surveys, and forecast the measurements of shear, galaxy, and galaxy–galaxy lensing power spectra (i.e. 3 × 2 pt). The effects of neutrino, baryonic feedback, and uncertainties of intrinsic alignment, shear calibration, galaxy bias, and photometric redshift are also included in the analysis. After performing a joint constraint on all the cosmological and systematical parameters based on the simulated data from the theoretical prediction, we obtain a lower limit of the ULA particle mass log10 (ma/eV) ≥ −22.5 and an upper limit of the ULA fraction fa ≤ 0.83 at 95 per cent confidence level, and log10 (ma/eV) ≥ −21.9 with fa ≤ 0.77 when ignoring the baryonic feedback. We find that the CSST photometric surveys can improve the constraint on the ULA mass by about one order of magnitude, compared to the current constraints using the same kind of observational data.

Galaxy bias in the era of LSST: perturbative bias expansions

Upcoming imaging surveys will allow for high signal-to-noise measurements of galaxy clustering at small scales. In this work, we present the results of the Rubin Observatory Legacy Survey of Space and Time (LSST) bias challenge, the goal of which is to compare the performance of different nonlinear galaxy bias models in the context of LSST Year 10 (Y10) data. Specifically, we compare two perturbative approaches, Lagrangian perturbation theory (LPT) and Eulerian perturbation theory (EPT) to two variants of Hybrid Effective Field Theory (HEFT), with our fiducial implementation of these models including terms up to second order in the bias expansion as well as nonlocal bias and deviations from Poissonian stochasticity. We consider a variety of different simulated galaxy samples and test the performance of the bias models in a tomographic joint analysis of LSST-Y10-like galaxy clustering, galaxy-galaxy-lensing and cosmic shear. We find both HEFT methods as well as LPT and EPT combined with non-perturbative predictions for the matter power spectrum to yield unbiased constraints on cosmological parameters up to at least a maximal scale of k max = 0.4 Mpc-1 for all samples considered, even in the presence of assembly bias. While we find that we can reduce the complexity of the bias model for HEFT without compromising fit accuracy, this is not generally the case for the perturbative models. We find significant detections of non-Poissonian stochasticity in all cases considered, and our analysis shows evidence that small-scale galaxy clustering predominantly improves constraints on galaxy bias rather than cosmological parameters. These results therefore suggest that the systematic uncertainties associated with current nonlinear bias models are likely to be subdominant compared to other sources of error for tomographic analyses of upcoming photometric surveys, which bodes well for future galaxy clustering analyses using these high signal-to-noise data.

Euclid: Validation of the MontePython forecasting tools

Context. The Euclid mission of the European Space Agency will perform a survey of weak lensing cosmic shear and galaxy clustering in order to constrain cosmological models and fundamental physics. Aims. We expand and adjust the mock Euclid likelihoods of the MontePython software in order to match the exact recipes used in previous Euclid Fisher matrix forecasts for several probes: weak lensing cosmic shear, photometric galaxy clustering, the cross-correlation between the latter observables, and spectroscopic galaxy clustering. We also establish which precision settings are required when running the Einstein–Boltzmann solvers CLASS and CAMB in the context of Euclid. Methods. For the minimal cosmological model, extended to include dynamical dark energy, we perform Fisher matrix forecasts based directly on a numerical evaluation of second derivatives of the likelihood with respect to model parameters. We compare our results with those of previously validated Fisher codes using an independent method based on first derivatives of the Euclid observables. Results. We show that such MontePython forecasts agree very well with previous Fisher forecasts published by the Euclid Collab oration, and also, with new forecasts produced by the CosmicFish code, now interfaced directly with the two Einstein–Boltzmann solvers CAMB and CLASS. Moreover, to establish the validity of the Gaussian approximation, we show that the Fisher matrix marginal error contours coincide with the credible regions obtained when running Monte Carlo Markov chains with MontePython while using the exact same mock likelihoods. Conclusions. The new Euclid forecast pipelines presented here are ready for use with additional cosmological parameters, in order to explore extended cosmological models.

Determining the baryon impact on the matter power spectrum with galaxy clusters

ABSTRACT The redistribution of baryonic matter in massive haloes through processes like active galactic nuclei feedback and star formation leads to a suppression of the matter power spectrum on small scales. This redistribution can be measured empirically via the gas and stellar mass fractions in galaxy clusters, and leaves imprints on their electron density profiles. We constrain two semi-analytical baryon correction models with a compilation of recent Bayesian population studies of galaxy groups and clusters sampling a mass range above ∼3 × 1013 M⊙, and with cluster gas density profiles derived from deep, high-resolution X-ray observations. We are able to fit all the considered observational data, but highlight some anomalies in the observations. The constraints allow us to place precise, physically informed priors on the matter power spectrum suppression. At a scale of k = 1 h Mpc−1 we find a suppression of $0.042^{+0.012}_{-0.014}$ ($0.049^{+0.016}_{-0.012}$), while at k = 3 h Mpc−1 we find $0.184^{+0.026}_{-0.031}$ ($0.179^{+0.018}_{-0.020}$), depending on the model used. In our fiducial setting, we also predict at 97.5 per cent credibility, that at scales k &amp;lt; 0.37 h Mpc−1 baryon feedback impacts the matter power less than 1 per cent. This puts into question if baryon feedback is the driving factor for the discrepancy between cosmic shear and primary CMB results. We independently confirm results on this suppression from small-scale cosmic shear studies, while we exclude some hydro-dynamical simulations with too strong and too weak baryonic feedback. Our empirical prediction of the power spectrum suppression shows that studies of galaxy groups and clusters will be instrumental in unlocking the cosmological constraining power of future cosmic shear experiments like Euclid and Rubin-LSST, and invites further investigation of the baryon correction models.

Cosmology from cross-correlation of ACT-DR4 CMB lensing and DES-Y3 cosmic shear

ABSTRACT Cross-correlation between weak lensing of the Cosmic Microwave Background (CMB) and weak lensing of galaxies offers a way to place robust constraints on cosmological and astrophysical parameters with reduced sensitivity to certain systematic effects affecting individual surveys. We measure the angular cross-power spectrum between the Atacama Cosmology Telescope (ACT) DR4 CMB lensing and the galaxy weak lensing measured by the Dark Energy Survey (DES) Y3 data. Our baseline analysis uses the CMB convergence map derived from ACT-DR4 and Planck data, where most of the contamination due to the thermal Sunyaev Zel’dovich effect is removed, thus avoiding important systematics in the cross-correlation. In our modelling, we consider the nuisance parameters of the photometric uncertainty, multiplicative shear bias and intrinsic alignment of galaxies. The resulting cross-power spectrum has a signal-to-noise ratio = 7.1 and passes a set of null tests. We use it to infer the amplitude of the fluctuations in the matter distribution (S8 ≡ σ8(Ωm/0.3)0.5 = 0.782 ± 0.059) with informative but well-motivated priors on the nuisance parameters. We also investigate the validity of these priors by significantly relaxing them and checking the consistency of the resulting posteriors, finding them consistent, albeit only with relatively weak constraints. This cross-correlation measurement will improve significantly with the new ACT-DR6 lensing map and form a key component of the joint 6×2pt analysis between DES and ACT.

A Lagrangian theory for galaxy shape statistics

We formulate the Lagrangian perturbation theory of galaxy intrinsic alignments and compute the resulting auto and cross power spectra of galaxy shapes, densities and matter to 1-loop order. Our model represents a consistent effective-theory description of galaxy shape including the resummation of long-wavelength displacements which damp baryon acoustic oscillations, and includes one linear, three quadratic and two cubic dimensionless bias coefficients at this order, along with counterterms and stochastic contributions whose structure we derive. We compare this Lagrangian model against the three-dimensional helicity spectra of halo shapes measured in N-body simulations by ref. [1] and find excellent agreement on perturbative scales while testing a number of more restrictive bias parametrizations. The calculations presented are immediately relevant to analyses of both cosmic shear surveys and spectroscopic shape measurements, and we make a fast FFTLog-based code spinosaurus publicly available at https://github.com/sfschen/spinosaurus.

A differentiable perturbation-based weak lensing shear estimator

ABSTRACT Upcoming imaging surveys will use weak gravitational lensing to study the large-scale structure of the Universe, demanding sub-per cent accuracy for precise cosmic shear measurements. We present a new differentiable implementation of our perturbation-based shear estimator (fpfs), using jax, which is publicly available as part of a new suite of analytic shear algorithms called anacal. This code can analytically calibrate the shear response of any non-linear observable constructed with the fpfs shapelets and detection modes utilizing autodifferentiation (ad), generalizing the formalism to include a family of shear estimators with corrections for detection and selection biases. Using the ad capability of jax, it calculates the full Hessian matrix of the non-linear observables, which improves the previously presented second-order noise bias correction in the shear estimation. As an illustration of the power of the new anacal framework, we optimize the effective galaxy number density in the space of the generalized shear estimators using an LSST-like galaxy image simulation for the 10 yr LSST. For the generic shear estimator, the magnitude of the multiplicative bias |m| is below 3 × 10−3 (99.7 per cent confidence interval), and the effective galaxy number density is improved by 5 per cent . We also discuss some planned future additions to the anacal software suite to extend its applicability beyond the fpfs measurements.

Euclid preparation. XXXIV. The effect of linear redshift-space distortions in photometric galaxy clustering and its cross-correlation with cosmic shear

The cosmological surveys that are planned for the current decade will provide us with unparalleled observations of the distribution of galaxies on cosmic scales, by means of which we can probe the underlying large-scale structure (LSS) of the Universe. This will allow us to test the concordance cosmological model and its extensions. However, precision pushes us to high levels of accuracy in the theoretical modelling of the LSS observables so that no biases are introduced into the estimation of the cosmological parameters. In particular, effects such as redshift-space distortions (RSD) can become relevant in the computation of harmonic-space power spectra even for the clustering of the photometrically selected galaxies, as has previously been shown in literature. In this work, we investigate the contribution of linear RSD, as formulated in the Limber approximation by a previous work, in forecast cosmological analyses with the photometric galaxy sample of the survey. We aim to assess their impact and to quantify the bias on the measurement of cosmological parameters that would be caused if this effect were neglected. We performed this task by producing mock power spectra for photometric galaxy clustering and weak lensing as is expected to be obtained from the survey. We then used a Markov chain Monte Carlo approach to obtain the posterior distributions of cosmological parameters from these simulated observations. When the linear RSD is neglected, significant biases are caused when galaxy correlations are used alone and when they are combined with cosmic shear in the so-called 3times 2pt approach. These biases can be equivalent to as much as $5\ when an underlying Lambda CDM cosmology is assumed When the cosmological model is extended to include the equation-of-state parameters of dark energy, the extension parameters can be shifted by more than $1\ sigma$.

Propagating photo-z uncertainties: a functional derivative approach

ABSTRACT Photometric redshifts are a key ingredient in the analysis and interpretation of large-scale structure (LSS) surveys. The accuracy and precision of these redshift estimates are directly linked to the constraining power of photometric surveys. It is hence necessary to define precision and accuracy requirements for the redshift calibration to not infer biased results in the final analysis. For weak gravitational lensing of the LSS, the photometry culminates in the estimation of the source-redshift distribution (SRD) in each of the tomographic bins used in the analysis. The focus has been on shifts of the mean of the SRDs and how well the calibration must be able to recover those. Since the estimated SRDs are usually given as a normalized histogram with corresponding errors, it would be advantageous to propagate these uncertainties accordingly to see whether the requirements of the given survey are indeed fulfilled. Here, we propose the use of functional derivatives to calculate the sensitivity of the final observables, for examples, the lensing angular power spectrum, with respect to the SRD at a specific redshift. This allows the propagation of arbitrarily shaped small perturbations to the SRD, without having to run the whole analysis pipeline for each realization again. We apply our method to a survey with Euclid-like specifications and demonstrate it with SRDs of the KV450 data set, recovering previous results. Lastly, we note that the moments of the SRD of order larger than two will probably not be relevant when propagating redshift uncertainties in cosmic shear analysis.

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

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