Weak Lensing Constraints Research Articles

Weak lensing constraints on dark matter-baryon interactions with 𝖭-body simulations and machine learning

We investigate the elastic scattering cross section between dark matter and protons using the DES Year 3 weak lensing data. This scattering induces a dark acoustic oscillation structure in the matter power spectra. To address non-linear effects at low redshift, we utilize principal component analysis alongside a limited set of N-body simulations, improving the reliability of our matter power spectrum prediction. We further perform a robust Markov Chain Monte Carlo analysis to derive the upper bounds on the DM-proton elastic scattering cross-section, assuming different velocity dependencies. Our results, presented as the first Frequentist upper limits, are compared with the ones obtained by Bayesian approach. Compared with the upper limits derived from the Planck cosmic microwave background data, our findings from DES Year 3 data exhibit improvements of up to a factor of five. In addition, we forecast the future sensitivities of the China Space Station Telescope, the upcoming capabilities of this telescope could improve the current limits by approximately one order of magnitude.

Inferring the impact of feedback on the matter distribution using the Sunyaev Zel’dovich effect: insights from CAMELS simulations and ACT + DES data

ABSTRACT Feedback from active galactic nuclei and stellar processes changes the matter distribution on small scales, leading to significant systematic uncertainty in weak lensing constraints on cosmology. We investigate how the observable properties of group-scale haloes can constrain feedback’s impact on the matter distribution using Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS). Extending the results of previous work to smaller halo masses and higher wavenumber, k, we find that the baryon fraction in haloes contains significant information about the impact of feedback on the matter power spectrum. We explore how the thermal Sunyaev Zel’dovich (tSZ) signal from group-scale haloes contains similar information. Using recent Dark Energy Survey weak lensing and Atacama Cosmology Telescope tSZ cross-correlation measurements and models trained on CAMELS, we obtain 10 per cent constraints on feedback effects on the power spectrum at $k \sim 5\, h\, {\rm Mpc}^{-1}$. We show that with future surveys, it will be possible to constrain baryonic effects on the power spectrum to $\mathcal {O}(\lt 1~{{\ \rm per\ cent}})$ at $k = 1\, h\, {\rm Mpc}^{-1}$ and $\mathcal {O}(3~{{\ \rm per\ cent}})$ at $k = 5\, h\, {\rm Mpc}^{-1}$ using the methods that we introduce here. Finally, we investigate the impact of feedback on the matter bispectrum, finding that tSZ observables are highly informative in this case.

MAximum-entropy ReconStruction (MARS): A New Strong-lensing Reconstruction Algorithm for the JWST Era

AbstractThe MAximum-entropy ReconStruction (MARS) method is a free-form strong-lensing (SL) reconstruction algorithm, which adopts the maximum cross-entropy as a regularization. MARS shows remarkable convergence of multiple images in both source (∼0.”02) and image planes (∼0.”05 – 0.”1) while suppressing spurious fluctuations. Although the reconstruction requires a large number of free parameters exceeding ∼19,000, our implementation through PyTorch can obtain the reconstruction within hours. From our test using the publicly available synthetic clusters, we have verified that the reconstructed radial mass profiles are consistent with the truth within 1 percent. This makes MARS one of the best-performing SL reconstruction methods. We apply MARS to the six Hubble Frontier Fields clusters and present new mass reconstruction results. We also reconstruct a mass model of Abell 2744 using both weak-lensing (WL) and SL data from the JWST observations, with the largest dataset of Abell 2744, including 286 SL multiple images and ∼350 arcmin−2 WL constraints.

Decaying warm dark matter revisited

Decaying dark matter models provide a physically motivated way of channeling energy between the matter and radiation sectors. In principle, this could affect the predicted value of the Hubble constant in such a way as to accommodate the discrepancies between CMB inferences and local measurements of the same. Here, we revisit the model of warm dark matter decaying non-relativistically to invisible radiation. In particular, we rederive the background and perturbation equations starting from a decaying neutrino model and describe a new, computationally efficient method of computing the decay product perturbations up to large multipoles. We conduct MCMC analyses to constrain all three model parameters, for the first time including the mass of the decaying species, and assess the ability of the model to alleviate the Hubble and σ 8 tensions, the latter being the discrepancy between the CMB and weak gravitational lensing constraints on the amplitude of matter fluctuations on an 8 h -1 Mpc-1 scale. We find that the model reduces the H 0 tension from ∼ 4σ to ∼ 3σ and neither alleviates nor worsens the S 8 ≡ σ 8 (Ω m /0.3)0.5 tension, ultimately showing only mild improvements with respect to ΛCDM. However, the values of the model-specific parameters favoured by data is found to be well within the regime of relativistic decays where inverse processes are important, rendering a conclusive evaluation of the decaying warm dark matter model open to future work.

Tightening weak lensing constraints on the ellipticity of galaxy-scale dark matter haloes

Cosmological simulations predict that galaxies are embedded into triaxial dark matter haloes, which appear approximately elliptical in projection. Weak gravitational lensing allows us to constrain these halo shapes and thereby test the nature of dark matter. Weak lensing has already provided robust detections of the signature of halo flattening at the mass scales of groups and clusters, whereas results for galaxies have been somewhat inconclusive. Here we combine data from five weak lensing surveys (NGVSLenS, KiDS/KV450, CFHTLenS, CS82, and RCSLenS, listed in order of most to least constraining) in order to tighten observational constraints on galaxy-scale halo ellipticity for photometrically selected lens samples. We constrain fh, the average ratio between the aligned component of the halo ellipticity and the ellipticity of the light distribution, finding fh = 0.303−0.079+0.080 for red lens galaxies and fh = 0.217−0.159+0.160 for blue lens galaxies when assuming elliptical Navarro-Frenk-White density profiles and a linear scaling between halo ellipticity and galaxy ellipticity. Our constraints for red galaxies constitute the currently most significant (3.8σ) systematics-corrected detection of the signature of halo flattening at the mass scale of galaxies. Our results are in good agreement with expectations from the Millennium Simulation that apply the same analysis scheme and incorporate models for galaxy–halo misalignment. Assuming these misalignment models and the analysis assumptions stated above are correct, our measurements imply an average dark matter halo ellipticity for the studied red galaxy samples of ⟨|ϵh|⟩ = 0.174 ± 0.046, where |ϵh| = (1 − q)/(1 + q) relates to the ratio q = b/a of the minor and major axes of the projected mass distribution. Similar measurements based on larger upcoming weak lensing data sets can help to calibrate models for intrinsic galaxy alignments, which constitute an important source of systematic uncertainty in cosmological weak lensing studies.

KiDS-1000 methodology: Modelling and inference for joint weak gravitational lensing and spectroscopic galaxy clustering analysis

We present the methodology for a joint cosmological analysis of weak gravitational lensing from the fourth data release of the ESO Kilo-Degree Survey (KiDS-1000) and galaxy clustering from the partially overlapping Baryon Oscillation Spectroscopic Survey (BOSS) and the 2-degree Field Lensing Survey (2dFLenS). Cross-correlations between BOSS and 2dFLenS galaxy positions and source galaxy ellipticities have been incorporated into the analysis, necessitating the development of a hybrid model of non-linear scales that blends perturbative and non-perturbative approaches, and an assessment of signal contributions by astrophysical effects. All weak lensing signals were measured consistently via Fourier-space statistics that are insensitive to the survey mask and display low levels of mode mixing. The calibration of photometric redshift distributions and multiplicative gravitational shear bias has been updated, and a more complete tally of residual calibration uncertainties was propagated into the likelihood. A dedicated suite of more than 20 000 mocks was used to assess the performance of covariance models and to quantify the impact of survey geometry and spatial variations of survey depth on signals and their errors. The sampling distributions for the likelihood and the χ2 goodness-of-fit statistic have been validated, with proposed changes for calculating the effective number of degrees of freedom. The prior volume was explicitly mapped, and a more conservative, wide top-hat prior on the key structure growth parameter S8 = σ8 (Ωm/0.3)1/2 was introduced. The prevalent custom of reporting S8 weak lensing constraints via point estimates derived from its marginal posterior is highlighted to be easily misinterpreted as yielding systematically low values of S8, and an alternative estimator and associated credible interval are proposed. Known systematic effects pertaining to weak lensing modelling and inference are shown to bias S8 by no more than 0.1 standard deviations, with the caveat that no conclusive validation data exist for models of intrinsic galaxy alignments. Compared to the previous KiDS analyses, S8 constraints are expected to improve by 20% for weak lensing alone and by 29% for the joint analysis.

Extended lens reconstructions with grale: exploiting time-domain, substructural, and weak lensing information

ABSTRACT The information about the mass density of galaxy clusters provided by the gravitational lens effect has inspired many inversion techniques. In this article, updates to the previously introduced method in grale are described, and explored in a number of examples. The first looks into a different way of incorporating time delay information, not requiring the unknown source position. It is found that this avoids a possible bias that leads to ‘overfocusing’ the images, i.e. providing source position estimates that lie in a considerably smaller region than the true positions. The second is inspired by previous reconstructions of the cluster of galaxies MACS J1149.6+2223, where a multiply imaged background galaxy contained a supernova, SN Refsdal, of which four additional images were produced by the presence of a smaller cluster galaxy. The inversion for the cluster as a whole was not able to recover sufficient detail interior to this quad. We show how constraints on such different scales, from the entire cluster to a single member galaxy, can now be used, allowing such small-scale substructures to be resolved. Finally, the addition of weak lensing information to this method is investigated. While this clearly helps recover the environment around the strong lensing region, the mass sheet degeneracy may make a full strong and weak inversion difficult, depending on the quality of the ellipticity information at hand. We encounter ring-like structure at the boundary of the two regimes, argued to be the result of combining strong and weak lensing constraints, possibly affected by degeneracies.

Hints, neutrino bounds, and WDM constraints from SDSS DR14 Lyman-α and Planck full-survey data

The Lyman-α forest 1D flux power spectrum is a powerful probe of several cosmological parameters. Assuming a ΛCDM cosmology including massive neutrinos, we find that the latest SDSS DR14 BOSS and eBOSS Lyman-α forest data is in very good agreement with current weak lensing constraints on (Ωm, σ8) and has the same small level of tension with Planck. We did not identify a systematic effect in the data analysis that could explain this small tension, but we show that it can be reduced in extended cosmological models where the spectral index is not the same on the very different times and scales probed by CMB and Lyman-α data. A particular case is that of a ΛCDM model including a running of the spectral index on top of massive neutrinos. With combined Lyman-α and Planck data, we find a slight (3σ) preference for negative running, αs= −0.010 ± 0.004 (68%CL). Neutrino mass bounds are found to be robust against different assumptions. In the ΛCDM model with running, we find ∑ mν <0.11 eV at the 95% confidence level for combined Lyman-α and Planck (temperature and polarisation) data, or ∑ mν < 0.09 eV when adding CMB lensing and BAO data. We further provide strong and nearly model-independent bounds on the mass of thermal warm dark matter. For a conservative configuration consisting of SDSS data restricted to z<4.5 combined with XQ-100 Lyman-α data, we find mX > 5.3 keV (95%CL).

Hybrid-lenstool: a self-consistent algorithm to model galaxy clusters with strong- and weak-lensing simultaneously

ABSTRACT We present a new galaxy cluster lens modelling approach, hybrid-lenstool, that is implemented in the publicly available modelling software lenstool. hybrid-lenstool combines a parametric approach to model the core of the cluster, and a non-parametric (free-form) approach to model the outskirts. hybrid-lenstool optimizes both strong- and weak-lensing constraints simultaneously (Joint-Fit), providing a self-consistent reconstruction of the cluster mass distribution on all scales. In order to demonstrate the capabilities of the new algorithm, we tested it on a simulated cluster. hybrid-lenstool yields more accurate reconstructed mass distributions than the former Sequential-Fit approach where the parametric and the non-parametric models are optimized successively. Indeed, we show with the simulated cluster that the mass density profile reconstructed with a Sequential-Fit deviates from the input by 2–3σ at all scales while the Joint-Fit gives a profile that is within 1–1.5σ of the true value. This gain in accuracy is consequential for recovering mass distributions exploiting cluster lensing and therefore for all applications of clusters as cosmological probes. Finally we found that the Joint-Fit approach yields shallower slope of the inner density profile than the Sequential-Fit approach, thus revealing possible biases in previous lensing studies.

Weak lensing constraints on splashback around massive clusters

The splashback radius $r_\text{sp}$ separates the physical regimes of collapsed and infalling material around massive dark matter haloes. In cosmological simulations, this location is associated with a steepening of the spherically averaged density profile $\rho(r)$. In this work, we measure the splashback feature in the stacked weak gravitational lensing signal of $27$ massive clusters from the Cluster Canadian Comparison Project with careful control of residual systematics effects. We find that the shear introduced by the presence of additional structure along the line of sight significantly affects the noise at large clustercentric distances. Although we do not detect a significant steepening, the use of a simple parametric model enables us to measure both $r_\text{sp}=3.5^{+1.1}_{-0.7}$ comoving Mpc and the value of the logarithmic slope $\gamma = \log \rho / \log r$ at this point, $\gamma(r_\text{sp}) = -4.3^{+1.0}_{-1.5}$.

Mass and Light of Abell 370: A Strong and Weak Lensing Analysis

Abstract We present a new gravitational lens model of the Hubble Frontier Fields cluster Abell 370 (z = 0.375) using imaging and spectroscopy from Hubble Space Telescope and ground-based spectroscopy. We combine constraints from a catalog of 909 weakly lensed galaxies and 39 multiply imaged sources comprised of 114 multiple images, including a system of multiply imaged candidates at z = 7.84 ± 0.02, to obtain a best-fit mass distribution using the cluster lens modeling code Strong and Weak Lensing United. As the only analysis of A370 using strong and weak lensing constraints from Hubble Frontier Fields data, our method provides an independent check of assumptions on the mass distribution used in other methods. Convergence, shear, and magnification maps are made publicly available through the HFF website (http://www.stsci.edu/hst/campaigns/frontier-fields). We find that the model we produce is similar to models produced by other groups, with some exceptions due to the differences in lensing code methodology. In an effort to study how our total projected mass distribution traces light, we measure the stellar mass density distribution using Spitzer/Infrared Array Camera imaging. Comparing our total mass density to our stellar mass density in a radius of 0.3 Mpc, we find a mean projected stellar to total mass ratio of (stat.) using the diet Salpeter initial mass function. This value is in general agreement with independent measurements of in clusters of similar total mass and redshift.

Precise weak lensing constraints from deep high-resolution Ks images: VLT/HAWK-I analysis of the super-massive galaxy cluster RCS2 J 232727.7−020437 at z = 0.70

We demonstrate that deep good-seeing VLT/HAWK-I Ks images complemented with g + z-band photometry can yield a sensitivity for weak lensing studies of massive galaxy clusters at redshifts 0.7 ≲ z ≲ 1.1, which is almost identical to the sensitivity of HST/ACS mosaics of single-orbit depth. Key reasons for this good performance are the excellent image quality frequently achievable for Ks imaging from the ground, a highly effective photometric selection of background galaxies, and a galaxy ellipticity dispersion that is noticeably lower than for optically observed high-redshift galaxy samples. Incorporating results from the 3D-HST and UltraVISTA surveys we also obtained a more accurate calibration of the source redshift distribution than previously achieved for similar optical weak lensing data sets. Here we studied the extremely massive galaxy cluster RCS2 J232727.7−020437 (z = 0.699), combining deep VLT/HAWK-I Ks images (point spread function with a 0.′′35 full width at half maximum) with LBT/LBC photometry. The resulting weak lensing mass reconstruction suggests that the cluster consists of a single overdensity, which is detected with a peak significance of 10.1σ. We constrained the cluster mass to M200c/(1015 M⊙) = 2.06−0.26+0.28(stat.) ± 0.12(sys.) assuming a spherical Navarro, Frenk &amp; White model and simulation-based priors on the concentration, making it one of the most massive galaxy clusters known in the z ≳ 0.7 Universe. We also cross-checked the HAWK-I measurements through an analysis of overlapping HST/ACS images, yielding fully consistent estimates of the lensing signal.

THE GRISM LENS-AMPLIFIED SURVEY FROM SPACE (GLASS). VI. COMPARING THE MASS AND LIGHT IN MACS J0416.1-2403 USING FRONTIER FIELD IMAGING AND GLASS SPECTROSCOPY

ABSTRACT We present a model using both strong and weak gravitational lensing of the galaxy cluster MACS J0416.1-2403, constrained using spectroscopy from the Grism Lens-Amplified Survey from Space (GLASS) and Hubble Frontier Fields (HFF) imaging data. We search for emission lines in known multiply imaged sources in the GLASS spectra, obtaining secure spectroscopic redshifts of 30 multiple images belonging to 15 distinct source galaxies. The GLASS spectra provide the first spectroscopic measurements for five of the source galaxies. The weak lensing signal is acquired from 884 galaxies in the F606W HFF image. By combining the weak lensing constraints with 15 multiple image systems with spectroscopic redshifts and nine multiple image systems with photometric redshifts, we reconstruct the gravitational potential of the cluster on an adaptive grid. The resulting map of total mass density is compared with a map of stellar mass density obtained from the deep Spitzer Frontier Fields imaging data to study the relative distribution of stellar and total mass in the cluster. We find that the projected stellar mass to total mass ratio, f ⋆, varies considerably with the stellar surface mass density. The mean projected stellar mass to total mass ratio is (stat.), but with a systematic error as large as 0.004–0.005, dominated by the choice of the initial mass function. We find agreement with several recent measurements of f ⋆ in massive cluster environments. The lensing maps of convergence, shear, and magnification are made available to the broader community in the standard HFF format.

Cosmological hints of modified gravity?

The recent measurements of Cosmic Microwave Background temperature and polarization anisotropies made by the Planck satellite have provided impressive confirmation of the $\Lambda$CDM cosmological model. However interesting hints of slight deviations from $\Lambda$CDM have been found, including a $95 \%$ c.l. preference for a "modified gravity" structure formation scenario. In this paper we confirm the preference for a modified gravity scenario from Planck 2015 data, find that modified gravity solves the so-called $A_{lens}$ anomaly in the CMB angular spectrum, and constrains the amplitude of matter density fluctuations to $\sigma_8=0.815_{-0.048}^{+0.032}$, in better agreement with weak lensing constraints. Moreover, we find a lower value for the reionization optical depth of $\tau=0.059\pm0.020$ (to be compared with the value of $\tau= 0.079 \pm 0.017$ obtained in the standard scenario), more consistent with recent optical and UV data. We check the stability of this result by considering possible degeneracies with other parameters, including the neutrino effective number, the running of the spectral index and the amount of primordial helium. The indication for modified gravity is still present at about $95\%$ c.l., and could become more significant if lower values of $\tau$ were to be further confirmed by future cosmological and astrophysical data.

CFHTLenS: weak lensing constraints on the ellipticity of galaxy-scale matter haloes and the galaxy-halo misalignment

We present weak lensing constraints on the ellipticity of galaxy-scale matter haloes and the galaxy-halo misalignment. Using data from the Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS), we measure the weighted-average ratio of the aligned projected ellipticity components of galaxy matter haloes and their embedded galaxies, $f_\mathrm{h}$, split by galaxy type. We then compare our observations to measurements taken from the Millennium Simulation, assuming different models of galaxy-halo misalignment. Using the Millennium Simulation we verify that the statistical estimator used removes contamination from cosmic shear. We also detect an additional signal in the simulation, which we interpret as the impact of intrinsic shape-shear alignments between the lenses and their large-scale structure environment. These alignments are likely to have caused some of the previous observational constraints on $f_\mathrm{h}$ to be biased high. From CFHTLenS we find $f_\mathrm{h}=-0.04 \pm 0.25$ for early-type galaxies, which is consistent with current models for the galaxy-halo misalignment predicting $f_\mathrm{h}\simeq 0.20$. For late-type galaxies we measure $f_\mathrm{h}=0.69_{-0.36}^{+0.37}$ from CFHTLenS. This can be compared to the simulated results which yield $f_\mathrm{h}\simeq 0.02$ for misaligned late-type models.

Gravitational lensing properties of an isothermal universal halo profile

N-body simulations predict that dark matter halos with different mass scales are described by a universal model, the Navarro-Frenk-White (NFW) density profiles. As a consequence of baryonic cooling effects, these halos will become more concentrated, and similar to an isothermal sphere over a large range in radii (∼ 300 h−1 kpc). The singular isothermal sphere (SIS) model however has to be truncated artificially at large radii since it extends to infinity. We model a massive galaxy halo as a combination of an isothermal sphere and an NFW density profile. We give an approximation for the mass concentration at different baryon fractions and present exact expressions for the weak lensing shear and flexion for such a halo. We compare the lensing properties with the SIS and NFW profiles. We find that the combined profile can generate higher order lensing signals at small radii and is more efficient in generating strong lensing events. In order to distinguish such a halo profile from the SIS or NFW profiles, one needs to combine strong and weak lensing constraints for small and large radii.

THE ATACAMA COSMOLOGY TELESCOPE: A MEASUREMENT OF THE PRIMORDIAL POWER SPECTRUM

We present constraints on the primordial power spectrum of adiabatic fluctuations using data from the 2008 Southern Survey of the Atacama Cosmology Telescope (ACT). The angular resolution of ACT provides sensitivity to scales beyond \ell = 1000 for resolution of multiple peaks in the primordial temperature power spectrum, which enables us to probe the primordial power spectrum of adiabatic scalar perturbations with wavenumbers up to k \simeq 0.2 Mpc^{-1}. We find no evidence for deviation from power-law fluctuations over two decades in scale. Matter fluctuations inferred from the primordial temperature power spectrum evolve over cosmic time and can be used to predict the matter power spectrum at late times; we illustrate the overlap of the matter power inferred from CMB measurements (which probe the power spectrum in the linear regime) with existing probes of galaxy clustering, cluster abundances and weak lensing constraints on the primordial power. This highlights the range of scales probed by current measurements of the matter power spectrum.

Future weak lensing constraints in a dark coupled universe

Coupled cosmologies can predict values for the cosmological parameters at low redshifts which may differ substantially from the parameters values within non-interacting cosmologies. Therefore, low redshift probes, as the growth of structure and the dark matter distribution via galaxy and weak lensing surveys constitute a unique tool to constrain interacting dark sector models. We focus here on weak lensing forecasts from future Euclid and LSST-like surveys combined with the ongoing Planck cosmic microwave background experiment. We find that these future data could constrain the dimensionless coupling to be smaller than a few $\times 10^{-2}$. The coupling parameter $\xi$ is strongly degenerate with the cold dark matter energy density $\Omega_{c}h^2$ and the Hubble constant $H_0$.These degeneracies may cause important biases in the cosmological parameter values if in the universe there exists an interaction among the dark matter and dark energy sectors.

COSMOS weak-lensing constraints on modified gravity

The observed acceleration of the universe, explained through dark energy, could alternatively be explained through a modification of gravity that would also induce modifications in the evolution of cosmological perturbations. We use new weak lensing data from the COSMOS survey to test for deviations from General Relativity. The departure from GR is parametrized in a model-independent way that consistently parametrizes the two-point cosmic shear amplitude and growth. Using CMB priors, we perform a likelihood analysis. We find constraints on the amplitude of the signal that do not indicate a deviation from General Relativity.

THE DARK MATTER DISTRIBUTION IN A383: EVIDENCE FOR A SHALLOW DENSITY CUSP FROM IMPROVED LENSING, STELLAR KINEMATIC, AND X-RAY DATA

We extend our analyses of the dark matter (DM) distribution in relaxed clusters to the case of Abell 383, a luminous X-ray cluster at z=0.189 with a dominant central galaxy and numerous strongly-lensed features. Following our earlier papers, we combine strong and weak lensing constraints secured with Hubble Space Telescope and Subaru imaging with the radial profile of the stellar velocity dispersion of the central galaxy, essential for separating the baryonic mass distribution in the cluster core. Hydrostatic mass estimates from Chandra X-ray observations further constrain the solution. These combined datasets provide nearly continuous constraints extending from 2 kpc to 1.5 Mpc in radius, allowing stringent tests of results from recent numerical simulations. Two key improvements in our data and its analysis make this the most robust case yet for a shallow slope \beta of the DM density profile \rho_DM ~ r^-\beta on small scales. First, following deep Keck spectroscopy, we have secured the stellar velocity dispersion profile to a radius of 26 kpc for the first time in a lensing cluster. Secondly, we improve our previous analysis by adopting a triaxial DM distribution and axisymmetric dynamical models. We demonstate that in this remarkably well-constrained system, the logarithmic slope of the DM density at small radii is \beta < 1.0 (95% confidence). An improved treatment of baryonic physics is necessary, but possibly insufficient, to reconcile our observations with the recent results of high-resolution simulations.