Weak Lensing Observations Research Articles

Inferring high-redshift large-scale structure dynamics from the Lyman-αforest

One of the major science goals over the coming decade is to test fundamental physics with probes of the cosmic large-scale structure out to high redshift. Here we present a fully Bayesian approach to infer the three-dimensional cosmic matter distribution and its dynamics atz > 2 from observations of the Lyman-αforest. We demonstrate that the method recovers the unbiased mass distribution and the correct matter power spectrum at all scales. Our method infers the three-dimensional density field from a set of one-dimensional spectra, interpolating the information between the lines of sight. We show that our algorithm provides unbiased mass profiles of clusters, becoming an alternative for estimating cluster masses complementary to weak lensing or X-ray observations. The algorithm employs a Hamiltonian Monte Carlo method to generate realizations of initial and evolved density fields and the three-dimensional large-scale flow, revealing the cosmic dynamics at high redshift. The method correctly handles multi-modal parameter distributions, which allow constraining the physics of the intergalactic medium with high accuracy. We performed several tests using realistic simulated quasar spectra to test and validate our method. Our results show that detailed and physically plausible inference of three-dimensional large-scale structures at high redshift has become feasible.

Covariance matrices for galaxy cluster weak lensing: from virial regime to uncorrelated large-scale structure

ABSTRACT Next-generation optical imaging surveys will revolutionize the observations of weak gravitational lensing by galaxy clusters and provide stringent constraints on growth of structure and cosmic acceleration. In these experiments, accurate modelling of covariance matrices of cluster weak lensing plays the key role in obtaining robust measurements of the mean mass of clusters and cosmological parameters. We use a combination of analytical calculations and high-resolution N-body simulations to derive accurate covariance matrices that span from the virial regime to linear scales of the cluster-matter cross-correlation. We validate this calculation using a public ray-tracing lensing simulation and provide a software package for calculating covariance matrices for a wide range of cluster and source sample choices. We discuss the relative importance of shape noise and density fluctuations, the impact of radial bin size, and the impact of off-diagonal elements. For a weak lensing source density ns = 10 arcmin−2, shape noise typically dominates the variance on comoving scales $r_{\rm p}\lesssim 5\ h^{-1} \, \rm Mpc$. However, for ns = 60 arcmin−2, potentially achievable with future weak lensing experiments, density fluctuations typically dominate the variance at $r_{\rm p}\gtrsim 1\ h^{-1} \, \rm Mpc$ and remain comparable to shape noise on smaller scales.

Weak lensing cosmology with convolutional neural networks on noisy data

ABSTRACT Weak gravitational lensing is one of the most promising cosmological probes of the late universe. Several large ongoing (DES, KiDS, HSC) and planned (LSST, Euclid, WFIRST) astronomical surveys attempt to collect even deeper and larger scale data on weak lensing. Due to gravitational collapse, the distribution of dark matter is non-Gaussian on small scales. However, observations are typically evaluated through the two-point correlation function of galaxy shear, which does not capture non-Gaussian features of the lensing maps. Previous studies attempted to extract non-Gaussian information from weak lensing observations through several higher order statistics such as the three-point correlation function, peak counts, or Minkowski functionals. Deep convolutional neural networks (CNN) emerged in the field of computer vision with tremendous success, and they offer a new and very promising framework to extract information from 2D or 3D astronomical data sets, confirmed by recent studies on weak lensing. We show that a CNN is able to yield significantly stricter constraints of (σ8, Ωm) cosmological parameters than the power spectrum using convergence maps generated by full N-body simulations and ray-tracing, at angular scales and shape noise levels relevant for future observations. In a scenario mimicking LSST or Euclid, the CNN yields 2.4–2.8 times smaller credible contours than the power spectrum, and 3.5–4.2 times smaller at noise levels corresponding to a deep space survey such as WFIRST. We also show that at shape noise levels achievable in future space surveys the CNN yields 1.4–2.1 times smaller contours than peak counts, a higher order statistic capable of extracting non-Gaussian information from weak lensing maps.

The large-scale general-relativistic correction for Newtonian mocks

We clarify the subtle issue of finding the correct mapping of Newtonian simulations to light-cone observables at very large distance scales. A faithful general-relativistic interpretation specifies a gauge, i.e. a chart that relates the simulation data to points of the space-time manifold. It has already been pointed out that the implicit gauge choice of Newtonian simulations is indeed different from the Poisson gauge that is commonly adopted for relativistic calculations, the difference being most significant at large scales. It is therefore inconsistent, for example, to predict weak-lensing observables from simulations unless this gauge issue is properly accounted for. Using perturbation theory as well as fully relativistic N-body simulations we quantify the systematic error introduced this way, and we discuss several solutions that would render the calculations relativistically self-consistent.

The impact of baryonic physics and massive neutrinos on weak lensing peak statistics

ABSTRACT We study the impact of baryonic processes and massive neutrinos on weak lensing peak statistics that can be used to constrain cosmological parameters. We use the BAHAMAS suite of cosmological simulations, which self-consistently include baryonic processes and the effect of massive neutrino free-streaming on the evolution of structure formation. We construct synthetic weak lensing catalogues by ray tracing through light-cones, and use the aperture mass statistic for the analysis. The peaks detected on the maps reflect the cumulative signal from massive bound objects and general large-scale structure. We present the first study of weak lensing peaks in simulations that include both baryonic physics and massive neutrinos (summed neutrino mass Mν = 0.06, 0.12, 0.24, and 0.48 eV assuming normal hierarchy), so that the uncertainty due to physics beyond the gravity of dark matter can be factored into constraints on cosmological models. Assuming a fiducial model of baryonic physics, we also investigate the correlation between peaks and massive haloes, over a range of summed neutrino mass values. As higher neutrino mass tends to suppress the formation of massive structures in the Universe, the halo mass function and lensing peak counts are therefore modified as a function of Mν. Over most of the S/N range, the impact of fiducial baryonic physics is greater (less) than neutrinos for 0.06 and 0.12 (0.24 and 0.48) eV models. Both baryonic physics and massive neutrinos should be accounted for when deriving cosmological parameters from weak lensing observations.

Distinguishing standard and modified gravity cosmologies with machine learning

We present a convolutional neural network to classify distinct cosmological scenarios based on the statistically similar weak-lensing maps they generate. Modified gravity (MG) models that include massive neutrinos can mimic the standard concordance model [Lambda cold dark matter ($\mathrm{\ensuremath{\Lambda}}\mathrm{CDM}$)] in terms of Gaussian weak-lensing observables. An inability to distinguish viable models that are based on different physics potentially limits a deeper understanding of the fundamental nature of cosmic acceleration. For a fixed redshift of sources, we demonstrate that a machine learning network trained on simulated convergence maps can discriminate between such models better than conventional higher-order statistics. Results improve further when multiple source redshifts are combined. To accelerate training, we implement a novel data compression strategy that incorporates our prior knowledge of the morphology of typical convergence map features. Our method fully distinguishes $\mathrm{\ensuremath{\Lambda}}\mathrm{CDM}$ from its most similar MG model on noise-free data, and it correctly identifies among the MG models with at least 80% accuracy when using the full redshift information. Adding noise lowers the correct classification rate of all models, but the neural network still significantly outperforms the peak statistics used in a previous analysis.

Physical Bayesian modelling of the non-linear matter distribution: New insights into the nearby universe

Accurate analyses of present and next-generation cosmological galaxy surveys require new ways to handle effects of non-linear gravitational structure formation processes in data. To address these needs we present an extension of our previously developed algorithm for Bayesian Origin Reconstruction from Galaxies (BORG) to analyse matter clustering at non-linear scales in observations. This is achieved by incorporating a numerical particle mesh model of gravitational structure formation into our Bayesian inference framework. The algorithm simultaneously infers the three-dimensional primordial matter fluctuations from which present non-linear observations formed and provides reconstructions of velocity fields and structure formation histories. The physical forward modelling approach automatically accounts for the non-Gaussian features in gravitationally evolved matter density fields and addresses the redshift space distortion problem associated with peculiar motions of observed galaxies. Our algorithm employs a hierarchical Bayes approach to jointly account for various observational effects, such as unknown galaxy biases, selection effects, and observational noise. Corresponding parameters of the data model are marginalized out via a sophisticated Markov chain Monte Carlo approach relying on a combination of a multiple block sampling framework and an efficient implementation of a Hamiltonian Monte Carlo sampler. We demonstrate the performance of the method by applying it to the 2M++ galaxy compilation, tracing the matter distribution of the nearby universe. We show accurate and detailed inferences of the three-dimensional non-linear dark matter distribution of the nearby universe. As exemplified in the case of the Coma cluster, our method provides complementary mass estimates that are compatible with those obtained from weak lensing and X-ray observations. For the first time, we also present a reconstruction of the vorticity of the non-linear velocity field from observations. In summary, our method provides plausible and very detailed inferences of the dark matter and velocity fields of our cosmic neighbourhood.

Full-sky weak lensing: a nonlinear post-Friedmann treatment

We present a full-sky derivation of weak lensing observables in the Post-Friedmann (PF) formalism. Weak lensing has the characteristic of mixing small scales and large scales since it is affected by inhomogeneities integrated along the photon trajectory. With the PF formalism, we develop a modelling of lensing observables which encompasses both leading order relativistic effects and effects that are due to the fully non-linear matter distribution at small scales. We derive the reduced shear, convergence and rotation up to order 1/c4 in the PF approximation, accounting for scalar, vector and tensor perturbations, as well as galaxies' peculiar velocities. We discuss the various contributions that break the Kaiser-Squires relation between the shear and the convergence at different orders. We pay particular attention to the impact of the frame-dragging vector potential on lensing observables and we discuss potential ways to measure this effect in future lensing surveys.

Weak and strong deflection gravitational lensings by a charged Horndeski black hole

A charged black hole was predicted by the Einstein-Horndeski-Maxwell theory. In order to provide its observational signatures, we investigate its weak and strong deflection gravitational lensings. We find its weak deflection lensing observables, including the positions, magnifications and differential time delay of the lensed images. We also obtain its strong deflection lensing observables, including the apparent radius of the photon sphere as well as the angular separation, brightness difference and differential time delay between the relativistic images. Taking the supermassive black hole in the Galactic Center as the lens, we evaluate these observables and compare these signatures with those of the Schwarzschild, Reissner-Nordström, tidal Reissner-Nordström and charged Galileon black holes. After a detailed analysis of the feasibility of measuring these lensing observables, we conclude that although it is possible to detect some leading effects of the weak and strong deflection lensings by the charged Horndeski and other black holes with current technology, it would be unlikely to distinguish one kind of these black holes from the others based on these detections in the near future due to lack of enough highly angular resolution in astronomical observations to tell their differences.

Quantifying baryon effects on the matter power spectrum and the weak lensing shear correlation

Feedback processes from baryons are expected to strongly affect weak-lensing observables of current and future cosmological surveys. In this paper we present a new parametrisation of halo profiles based on gas, stellar, and dark matter density components. This parametrisation is used to modify outputs of gravity-only N-body simulations (following the prescription of [1]) in order to mimic baryonic effects on the matter density field. The resulting baryonic correction model relies on a few well motivated physical parameters and is able to reproduce the redshift zero clustering signal of hydrodynamical simulations at two percent accuracy below k∼10 h/Mpc. A detailed study of the baryon suppression effects on the matter power spectrum and the weak lensing shear correlation reveals that the signal is dominated by two parameters describing the slope of the gas profile in haloes and the maximum radius of gas ejection. We show that these parameters can be constrained with the observed gas fraction of galaxy groups and clusters from X-ray data. Based on these observations we predict a beyond percent effect on the power spectrum above k=0.2–1.0 h/Mpc with a maximum suppression of 15–25 percent around k∼ 10 h/Mpc. As a result, the weak lensing angular shear power spectrum is suppressed by 15–25 percent at scales beyond ℓ∼ 100–600 and the shear correlations ξ+ and ξ− are affected at the 10–25 percent level below 5 and 50 arc-minutes, respectively. The relatively large uncertainties of these predictions are a result of the poorly known hydrostatic mass bias of current X-ray observations as well as the generic difficulty to observe the low density gas outside of haloes.

Cosmological Tests of Gravity with the Latest Observations

Abstract We perform observational tests of modified gravity on cosmological scales following model-dependent and model-independent approaches using the latest astronomical observations, including measurements of the local Hubble constant, cosmic microwave background, the baryonic acoustic oscillations and redshift space distortions derived from galaxy surveys including the Sloan Digital Sky Survey BOSS and eBOSS, as well as the weak lensing observations performed by the CFHTLenS team. Combining all these data sets, we find a deviation from the prediction of general relativity in both the effective Newton’s constant, μ(a, k), and in the gravitational slip, η(a, k). The deviation is at a 3.1σ level in the joint {μ(a, k), η(a, k)} space using a two-parameter phenomenological model for μ and η, and reaches a 3.7σ level if a general parameterization is used. This signal, which may be subject to unknown observational systematics, or a sign of new physics, is worth further investigating with forthcoming observations.

Weak lensing distortions beyond shear

When a luminous source is extended, its distortions by weak gravitational lensing are richer than a mere combination of magnification and shear. In a recent work, we proposed an elegant formalism based on complex analysis to describe and calculate such distortions. The present article further elaborates this finite-beam approach, and applies it to a realistic cosmological model. In particular, the cosmic correlations of image distortions beyond shear are predicted for the first time. These constitute new weak-lensing observables, sensitive to very-small-scale features of the distribution of matter in the Universe. While the major part of the analysis is performed in the approximation of circular sources, a general method for extending it to noncircular sources is presented and applied to the astrophysically relevant case of elliptic sources.

The Massively Accreting Cluster A2029

Abstract We explore the structure of galaxy cluster A2029 and its surroundings based on intensive spectroscopy along with X-ray and weak lensing observations. The redshift survey includes 4376 galaxies (1215 spectroscopic cluster members) within 40′ of the cluster center; the redshifts are included here. Two subsystems, A2033 and a southern infalling group (SIG), appear in the infall region based on the spectroscopy, as well as on the weak lensing and X-ray maps. The complete redshift survey of A2029 also identifies at least 12 foreground and background systems (10 are extended X-ray sources) in the A2029 field; we include a census of their properties. The X-ray luminosities (L X )–velocity dispersions ( ) scaling relations for A2029, A2033, SIG, and the foreground/background systems are consistent with the known cluster scaling relations. The combined spectroscopy, weak lensing, and X-ray observations provide a robust measure of the masses of A2029, A2033, and SIG. The total mass of the infalling groups (A2033 and SIG) is ∼60% of the M 200 of the primary cluster, A2029. Simple dynamical considerations suggest that A2029 will accrete these subsystems in the next few Gyr. In agreement with simulations and other clusters observed in a similar redshift range, the total mass in the A2029 infall region is comparable to the A2029 M 200 and will mostly be accreted in the long-term future.

On the degeneracy between baryon feedback and massive neutrinos as probed by matter clustering and weak lensing

Massive neutrinos, due to their free streaming, produce a suppression in the matter power spectrum at intermediate and small scales which could be probed by galaxy clustering and/or weak lensing observables. This effect happens at scales that are also influenced by baryon feedback, i.e. galactic winds or Active Galactic Nuclei (AGN) feedback, which in realistic hydrodynamic simulations has also been shown to produce a suppression of power. Leaving aside, for the moment, the complex issue of galaxy bias, we focus here on matter clustering and tomographic weak lensing, we investigate the possible degeneracy between baryon feedback and neutrinos showing that it is not likely to degrade significantly the measurement of neutrino mass in future surveys. To do so, we generate mock data sets and fit them using the Markov Chain Monte Carlo (MCMC) technique and explore degeneracies between feedback parameters and neutrino mass. We model baryon feedback through fitting functions, while massive neutrinos are accounted for, also in the non-linear regime, using HALOFIT calibrated against accurate N-body neutrino simulations. In the error budget, we include the uncertainty in the modelling of non-linearities. For both matter clustering and weak lensing, we always recover the input neutrino mass within ∼0.25σ confidence level. Finally, we also take into account the intrinsic alignment effect in the weak lensing mock data. Even in this case, we are able to recover the right parameters: in particular, we find a significant degeneracy pattern between Mν and the intrinsic alignment parameter AIA.

Sunyaev–Zel’dovich effect and X-ray scaling relations from weak lensing mass calibration of 32 South Pole Telescope selected galaxy clusters

Uncertainty in the mass-observable scaling relations is currently the limiting factor for galaxy cluster based cosmology. Weak gravitational lensing can provide a direct mass calibration and reduce the mass uncertainty. We present new ground-based weak lensing observations of 19 South Pole Telescope (SPT) selected clusters at redshifts $0.29 \leq z \leq 0.61$ and combine them with previously reported space-based observations of 13 galaxy clusters at redshifts $0.576 \leq z \leq 1.132$ to constrain the cluster mass scaling relations with the Sunyaev-Zel'dovich effect (SZE), the cluster gas mass \mgas, and \yx, the product of \mgas\ and X-ray temperature. We extend a previously used framework for the analysis of scaling relations and cosmological constraints obtained from SPT-selected clusters to make use of weak lensing information. We introduce a new approach to estimate the effective average redshift distribution of background galaxies and quantify a number of systematic errors affecting the weak lensing modelling. These errors include a calibration of the bias incurred by fitting a Navarro-Frenk-White profile to the reduced shear using $N$-body simulations. We blind the analysis to avoid confirmation bias. We are able to limit the systematic uncertainties to 5.6% in cluster mass (68% confidence). Our constraints on the mass--X-ray observable scaling relations parameters are consistent with those obtained by earlier studies, and our constraints for the mass--SZE scaling relation are consistent with the simulation-based prior used in the most recent SPT-SZ cosmology analysis. We can now replace the external mass calibration priors used in previous SPT-SZ cosmology studies with a direct, internal calibration obtained on the same clusters.

Modified gravity black hole lensing observables in weak and strong field of gravity

We extend a recent work on weak field first order light deflection in the MOdified Gravity (MOG) by comprehensively analyzing the actual observables in gravitational lensing both in the weak and strong field regime. The static spherically symmetric black hole (BH) obtained by Moffat is what we call here the Schwarzschild-MOG (abbreviated as SMOG) containing repulsive Yukawa-like force characterized by the MOG parameter $\alpha>0$ diminishing gravitational attraction. We point out a remarkable feature of SMOG, viz., it resembles a regular \textit{brane-world} BH in the range $-1<\alpha <0$ giving rise to a negative "tidal charge" $Q$ ($=\frac{1}{4}\frac{\alpha }{1+\alpha}$) interpreted as an imprint from the $5D$ bulk with an imaginary source charge $q$ in the brane. The Yukawa-like force of MOG is attractive in the brane-world range enhancing gravitational attraction. For $-\infty <\alpha <-1$, the SMOG represents a naked singularity. Specifically, we shall investigate the effect of $\alpha $ or Yukawa-type forces on the weak (up to third PPN order) and strong field lensing observables. For illustration, we consider the supermassive BH SgrA* with $\alpha =0.055$ for the weak field to quantify the deviation of observables from GR but in general we leave $\alpha$ unrestricted both in sign and magnitude so that future accurate lensing measurements, which are quite challenging, may constrain $\alpha$.

Weak lensing light-cones in modified gravity simulations with and without massive neutrinos

We present a novel suite of cosmological N-body simulations called the DUSTGRAIN-pathfinder, implementing simultaneously the effects of an extension to General Relativity in the form of $f(R)$ gravity and of a non-negligible fraction of massive neutrinos. We describe the generation of simulated weak lensing and cluster counts observables within a past light-cone extracted from these simulations. The simulations have been performed by means of a combination of the MG-GADGET code and a particle-based implementation of massive neutrinos, while the light-cones have been generated using the MapSim pipeline allowing us to compute weak lensing maps through a ray-tracing algorithm for different values of the source plane redshift. The mock observables extracted from our simulations will be employed for a series of papers focussed on understanding and possibly breaking the well-known observational degeneracy between $f(R)$ gravity and massive neutrinos, i.e. the fact that some specific combinations of the characteristic parameters for these two phenomena (the $f_{R0}$ scalar amplitude and the total neutrino mass $\Sigma m_{\nu}$) may result indistinguishable from the standard $\mathrm{\Lambda CDM}$ cosmology through several standard observational probes. In particular, in the present work we show how a tomographic approach to weak lensing statistics could allow - especially for the next generation of wide-field surveys - to disentangle some of the models that appear statistically indistinguishable through standard single-redshift weak lensing probe.

Multiwavelength scaling relations in galaxy groups: a detailed comparison of GAMA and KiDS observations to BAHAMAS simulations

We study the scaling relations between the baryonic content and total mass of groups of galaxies, as these systems provide a unique way to examine the role of non-gravitational processes in structure formation. Using Planck and ROSAT data, we conduct detailed comparisons of the stacked thermal Sunyaev-Zel'dovich (tSZ) and X-ray scaling relations of galaxy groups found in the Galaxy And Mass Assembly (GAMA) survey and the BAHAMAS hydrodynamical simulation. We use weak gravitational lensing data from the Kilo Degree Survey (KiDS) to determine the average halo mass of the studied systems. We analyse the simulation in the same way, using realistic weak lensing, X-ray, and tSZ synthetic observations. Furthermore, to keep selection biases under control, we employ exactly the same galaxy selection and group identification procedures to the observations and simulation. Applying this comparison, we find that the simulations reproduce the richness, size, and stellar mass functions of GAMA groups, as well as the stacked weak lensing and tSZ signals in bins of group stellar mass. However, the simulations predict X-ray luminosities that are higher than observed for this optically-selected group sample. As the same simulations were previously shown to match the luminosities of X-ray-selected groups, this suggests that X-ray-selected systems may form a biased subset. Finally, we demonstrate that our observational processing of the X-ray and tSZ signals is free of significant biases. We find that our optical group selection procedure has, however, some room for improvement.

Jacobi mapping approach for a precise cosmological weak lensing formalism

Cosmological weak lensing has been a highly successful and rapidly developing research field since the first detection of cosmic shear in 2000. However, it has recently been pointed out in Yoo et al. that the standard weak lensing formalism yields gauge-dependent results and, hence, does not meet the level of accuracy demanded by the next generation of weak lensing surveys. Here, we show that the Jacobi mapping formalism provides a solid alternative to the standard formalism, as it accurately describes all the relativistic effects contributing to the weak lensing observables. We calculate gauge-invariant expressions for the distortion in the luminosity distance, the cosmic shear components and the lensing rotation to linear order including scalar, vector and tensor perturbations. In particular, the Jacobi mapping formalism proves that the rotation is fully vanishing to linear order. Furthermore, the cosmic shear components contain an additional term in tensor modes which is absent in the results obtained with the standard formalism. Our work provides further support and confirmation of the gauge-invariant lensing formalism needed in the era of precision cosmology.

Halo Profiles and the Concentration–Mass Relation for a ΛCDM Universe

Abstract Profiles of dark matter-dominated halos at the group and cluster scales play an important role in modern cosmology. Using results from two very large cosmological N-body simulations, which increase the available volume at their mass resolution by roughly two orders of magnitude, we robustly determine the halo concentration–mass (c−M) relation over a wide range of masses, employing multiple methods of concentration measurement. We characterize individual halo profiles, as well as stacked profiles, relevant for galaxy–galaxy lensing and next-generation cluster surveys; the redshift range covered is 0 ≤ z ≤ 4, with a minimum halo mass of M 200c ∼ 2 × 1011 M ⊙. Despite the complexity of a proper description of a halo (environmental effects, merger history, nonsphericity, relaxation state), when the mass is scaled by the nonlinear mass scale M ⋆(z), we find that a simple non-power-law form for the c–M/M ⋆ relation provides an excellent description of our simulation results across eight decades in M/M ⋆ and for 0 ≤ z ≤ 4. Over the mass range covered, the c–M relation has two asymptotic forms: an approximate power law below a mass threshold M/M ⋆ ∼ 500–1000, transitioning to a constant value, c 0 ∼ 3 at higher masses. The relaxed halo fraction decreases with mass, transitioning to a constant value of ∼0.5 above the same mass threshold. We compare Navarro–Frenk–White (NFW) and Einasto fits to stacked profiles in narrow mass bins at different redshifts; as expected, the Einasto profile provides a better description of the simulation results. At cluster scales at low redshift, however, both NFW and Einasto profiles are in very good agreement with the simulation results, consistent with recent weak lensing observations.