Projected Matter Density Research Articles

AMICO galaxy clusters in KiDS-DR3: Measuring the splashback radius from weak gravitational lensing

Context. Weak gravitational lensing offers a powerful method to investigate the projected matter density distribution within galaxy clusters, granting crucial insights into the broader landscape of dark matter on cluster scales. Aims. In this study, we make use of the large photometric galaxy cluster data set derived from the publicly available Third Data Release of the Kilo-Degree Survey, along with the associated shear signal. Our primary objective is to model the peculiar sharp transition in the cluster profile slope, that is what is commonly referred to as the splashback radius. The data set under scrutiny includes 6962 galaxy clusters, selected by AMICO (an optimised detection algorithm of galaxy clusters) on the KiDS-DR3 data, in the redshift range of 0.1 ≤ z ≤ 0.6, all observed at a signal-to-noise ratio greater than 3.5. Methods. Employing a comprehensive Bayesian analysis, we model the stacked excess surface mass density distribution of the clusters. We adopt a model from recent results on numerical simulations that capture the dynamics of both orbiting and infalling materials, separated by the region where the density profile slope undergoes a pronounced deepening. Results. We find that the adopted profile successfully characterizes the cluster masses, consistent with previous works, and models the deepening of the slope of the density profiles measured with weak-lensing data up to the outskirts. Moreover, we measure the splashback radius of galaxy clusters and show that its value is close to the radius within which the enclosed overdensity is 200 times the mean matter density of the Universe, while theoretical models predict a larger value consistent with a low accretion rate. This points to a potential bias of optically selected clusters preferentially characterized by a high density at small scales compared to a pure mass-selected cluster sample.

Where shadows lie: reconstruction of anisotropies in the neutrino sky

The Cosmic Neutrino Background (CNB) encodes a wealth of information, but has not yet been observed directly. To determine the prospects of detection and to study its information content, we reconstruct the phase-space distribution of local relic neutrinos from the three-dimensional distribution of matter within 200 h -1 Mpc of the Milky Way. Our analysis relies on constrained realization simulations and forward modelling of the 2M++ galaxy catalogue. We find that the angular distribution of neutrinos is anti-correlated with the projected matter density, due to the capture and deflection of neutrinos by massive structures along the line of sight. Of relevance to tritium capture experiments, we find that the gravitational clustering effect of the large-scale structure on the local number density of neutrinos is more important than that of the Milky Way for neutrino masses less than 0.1 eV. Nevertheless, we predict that the density of relic neutrinos is close to the cosmic average, with a suppression or enhancement over the mean of (-0.3%, +7%, +27%) for masses of (0.01, 0.05, 0.1) eV. This implies no more than a marginal increase in the event rate for tritium capture experiments like PTOLEMY. We also predict that the CNB and CMB rest frames coincide for 0.01 eV neutrinos, but that neutrino velocities are significantly perturbed for masses larger than 0.05 eV. Regardless of mass, we find that the angle between the neutrino dipole and the ecliptic plane is small, implying a near-maximal annual modulation in the bulk velocity. Along with this paper, we publicly release our simulation data, comprising more than 100 simulations for six different neutrino masses.

Calibration of bias and scatter involved in cluster mass measurements using optical weak gravitational lensing

ABSTRACT Cosmological inference from cluster number counts is systematically limited by the accuracy of the mass calibration, i.e. the empirical determination of the mapping between cluster selection observables and halo mass. In this work we demonstrate a method to quantitatively determine the bias and uncertainties in weak-lensing (WL) mass calibration. To this end, we extract a library of projected matter density profiles from hydrodynamical simulations. Accounting for shear bias and noise, photometric redshift uncertainties, mis-centreing, cluster member contamination, cluster morphological diversity, and line-of-sight projections, we produce a library of shear profiles. Fitting a one-parameter model to these profiles, we extract the so-called WL mass MWL. Relating the WL mass to the halo mass from gravity-only simulations with the same initial conditions as the hydrodynamical simulations allows us to estimate the impact of hydrodynamical effects on cluster number counts experiments. Creating new shear libraries for ∼1000 different realizations of the systematics provides a distribution of the parameters of the WL to halo mass relation, reflecting their systematic uncertainty. This result can be used as a prior for cosmological inference. We also discuss the impact of the inner fitting radius on the accuracy, and determine the outer fitting radius necessary to exclude the signal from neighbouring structures. Our method is currently being applied to different Stage III lensing surveys, and can easily be extended to Stage IV lensing surveys.

Probability distribution function of the aperture mass field with large deviation theory

ABSTRACT In the context of tomographic cosmic shear surveys, a theoretical model for the one-point statistics of the aperture mass (Map) is developed. This formalism is based on the application of the large deviation principle to the projected matter density field and more specifically to the angular aperture masses. The latter holds the advantage of being an observable that can be directly extracted from the observed shear field and to be, by construction, independent from the long wave modes. Furthermore, we show that, with the help of a nulling procedure based on the so-called BNT transform, it is possible to build observables that depend only on a finite range of redshifts making them also independent from the small-scale modes. This procedure makes predictions for the shape of the one-point probability distribution function of such an observable very accurate, comparable to what had been previously obtained for 3D observables. Comparisons with specific simulations reveal however inconsistent results showing that synthetic lensing maps were not accurate enough for such refined observables. It points to the need for more precise dedicated numerical developments whose performances could be benchmarked with such observables. We furthermore review the possible systematics that could affect such a formalism in future weak-lensing surveys like Euclid, notably the impact of shape noise as well as leading corrections coming from lend–lens couplings, geodesic deviation, reduced shear and magnification bias.

Testing the lognormality of the galaxy and weak lensing convergence distributions from Dark Energy Survey maps

It is well known that the probability distribution function (PDF) of galaxy density contrast is approximately lognormal; whether the PDF of mass fluctuations derived from weak lensing convergence (kappa_WL) is lognormal is less well established. We derive PDFs of the galaxy and projected matter density distributions via the Counts in Cells (CiC) method. We use maps of galaxies and weak lensing convergence produced from the Dark Energy Survey (DES) Science Verification data over 139 deg^2. We test whether the underlying density contrast is well described by a lognormal distribution for the galaxies, the convergence and their joint PDF. We confirm that the galaxy density contrast distribution is well modeled by a lognormal PDF convolved with Poisson noise at angular scales from 10-40 arcmin (corresponding to physical scales of 3-10 Mpc). We note that as kappa_WL is a weighted sum of the mass fluctuations along the line of sight, its PDF is expected to be only approximately lognormal. We find that the kappa_WL distribution is well modeled by a lognormal PDF convolved with Gaussian shape noise at scales between 10 and 20 arcmin, with a best-fit chi^2/DOF of 1.11 compared to 1.84 for a Gaussian model, corresponding to p-values 0.35 and 0.07 respectively, at a scale of 10 arcmin. Above 20 arcmin a simple Gaussian model is sufficient. The joint PDF is also reasonably fitted by a bivariate lognormal. As a consistency check we compare the variances derived from the lognormal modelling with those directly measured via CiC. Our methods are validated against maps from the MICE Grand Challenge N-body simulation.

THE SIGNATURE OF THE WARM-HOT INTERGALACTIC MEDIUM INWMAPAND THE FORTHCOMINGPLANCKDATA

We compute the cross-correlation between the Warm-Hot Intergalactic Medium and maps of Cosmic Microwave Background temperature anisotropies using a log-normal probability density function to describe the weakly non-linear matter density field. We search for this contribution in the data measured by the Wilkinson Microwave Anisotropy Probe. We use a template of projected matter density reconstructed from the Two-Micron All-Sky Redshift Survey as a tracer of the electron distribution. The spatial distribution of filaments is modeled using the recently developed Augmented Lagrangian Perturbation Theory. On the scales considered here, the reconstructed density field is very well described by the assumed log-normal distribution function. We predict that the cross-correlation will have an amplitude of $0.03-0.3\mu$K. The measured value is close to $1.5\mu$K, compatible with random alignments between structure in the template and in the temperature anisotropy data. Using the W1 Differencing Assembly to remove this systematic gives a residual correlation dominated by Galactic foregrounds. Planck could detect the Warm-Hot Medium if it is well traced by the density field reconstructed from galaxy surveys. The 217GHz channel will allow to eliminate spurious contributions and its large frequency coverage can show the sign change from the Rayleigh-Jeans to the Wien part of the spectrum, characteristic of the thermal Sunyaev-Zeldovich effect.

Towards an understanding of third-order galaxy-galaxy lensing

Third-order galaxy-galaxy lensing (G3L) is a next generation galaxy-galaxy lensing technique that either measures the excess shear about lens pairs or the excess shear-shear correlations about lenses. It is clear that these statistics assess the three-point correlations between galaxy positions and projected matter density. For future applications of these novel statistics, we aim at a more intuitive understanding of G3L to isolate the main features that possibly can be measured. We construct a toy model ("isolated lens model"; ILM) for the distribution of galaxies and associated matter to determine the measured quantities of the two G3L correlation functions and traditional galaxy-galaxy lensing (GGL) in a simplified context. The ILM presumes single lens galaxies to be embedded inside arbitrary matter haloes that, however, are statistically independent ("isolated") from any other halo or lens position. In the ILM, the average mass-to-galaxy number ratio of clusters of any size cannot change. GGL and galaxy clustering alone cannot distinguish an ILM from any more complex scenario. The lens-lens-shear correlator in combination with second-order statistics enables us to detect deviations from a ILM, though. This can be quantified by a difference signal defined in the paper. We demonstrate with the ILM that this correlator picks up the excess matter distribution about galaxy pairs inside clusters. The lens-shear-shear correlator is sensitive to variations among matter haloes. In principle, it could be devised to constrain the ellipticities of haloes, without the need for luminous tracers, or maybe even random halo substructure. [Abridged]

Is Gravitational Lensing by Intercluster Filaments Always Negligible?

Intercluster filaments negligibly contribute to the weak lensing signal in general relativity (GR), $\gamma_{N}\sim 10^{-4}-10^{-3}$. In the context of relativistic modified Newtonian dynamics (MOND) introduced by Bekenstein, however, a single filament inclined by $\approx 45^\circ$ from the line of sight can cause substantial distortion of background sources pointing towards the filament's axis ($\kappa=\gamma=(1-A^{-1})/2\sim 0.01$); this is rigorous for infinitely long uniform filaments, but also qualitatively true for short filaments ($\sim 30$Mpc), and even in regions where the projected matter density of the filament is equal to zero. Since galaxies and galaxy are generally embedded in filaments or are projected on such structures, this contribution complicates the interpretation of the weak lensing shear map in the context of MOND. While our analysis is of mainly theoretical interest providing order-of-magnitude estimates only, it seems safe to conclude that when modeling systems with anomalous weak lensing signals, e.g. the bullet cluster of Clowe et al., the cosmic train wreck of Abell 520 from Mahdavi et al., and the dark clusters of Erben et al., filamentary structures might contribute in a significant and likely complex fashion. On the other hand, our predictions of a (conceptual) difference in the weak lensing signal could, in principle, be used to falsify MOND/TeVeS and its variations.

Poisson equation for weak gravitational lensing

Using the Newman and Penrose [E. T. Newman and R. Penrose, J. Math. Phys. (N.Y.) 3, 566 (1962).] spin-coefficient formalism, we examine the full Bianchi identities of general relativity in the context of gravitational lensing, where the matter and space-time curvature are projected into a lens plane perpendicular to the line of sight. From one component of the Bianchi identity, we provide a rigorous, new derivation of a Poisson equation for the projected matter density where the source term involves second derivatives of the observed weak gravitational lensing shear. We also show that the other components of the Bianchi identity reveal no new results. Numerical integration of the Poisson equation in test cases shows an accurate mass map can be constructed from the combination of a ground-based, wide-field image and a Hubble Space Telescope image of the same system.

Mapping Dark Matter with Cosmic Magnification

We develop a new tool to generate statistically precise dark matter maps from the cosmic magnification of galaxies with distance estimates. We show how to overcome the intrinsic clustering problem using the slope of the luminosity function, because magnification changes strongly over the luminosity function, while intrinsic clustering only changes weakly. This may allow precision cosmology beyond most current systematic limitations. The Square Kilometre Array is able to reconstruct the projected matter density map at smoothing scale approximately 10' with S/N > or = 1, at the rate of 200-4000 deg2 per year, depending on the abundance and evolution of 21 cm emitting galaxies. This power of mapping dark matter is comparable to, or even better than, that of cosmic shear from deep optical surveys or 21 cm surveys.

Cosmological parameters from combined second- and third-order aperture mass statistics of cosmic shear

We present predictions for cosmological parameter constraints from combined measurements of second- and third-order statistics of cosmic shear. We define the generalized third-order aperture mass statistics and show that it contains much more information about the bispectrum of the projected matter density than the skewness of the aperture mass. From theoretical models as well as from ΛCDM ray-tracing simulations, we calculate and and their dependence on cosmological parameters. The covariances including shot noise and cosmic variance of , and their cross-correlation are calculated using ray-tracing simulations. We perform an extensive Fisher matrix analysis, and for various combinations of cosmological parameters, we predict 1-σ-errors corresponding to measurements from a deep 29 square degree cosmic shear survey. Although the parameter degeneracies can not be lifted completely, the (linear) combination of second- and third-order aperture mass statistics reduces the errors significantly. The strong degeneracy between and , present for all second-order cosmic shear measures, is diminished substantially, whereas less improvement is found for the near-degenerate pair consisting of the shape parameter Γ and the spectral index ns. Uncertainties in the source galaxy redshift z0 increase the errors of all other parameters.

Lensing effect on polarization in the microwave background: Extracting the convergence power spectrum

Matter inhomogeneities along the line of sight deflect the cosmic microwave background (CMB) photons originating at the last scattering surface at a redshift $z\ensuremath{\sim}1100.$ These distortions modify the pattern of CMB polarization. We identify specific combinations of Stokes Q and U parameters that correspond to spin 0, \ifmmode\pm\else\textpm\fi{}2 variables and can be used to reconstruct the projected matter density. We compute the expected signal-to-noise ratio as a function of detector sensitivity and angular resolution. With the Planck satellite the detection would be at a few $\ensuremath{\sigma}$ level. Several times better detector sensitivity would be needed to measure the projected dark matter power spectrum over a wider range of scales, which could provide an independent confirmation of the projected matter power spectrum as measured from other methods.

Reconstructing projected matter density power spectrum from cosmic microwave background

Gravitational lensing distorts the cosmic microwave background (CMB) anisotropies and imprints a characteristic pattern onto it. The distortions depend on the projected matter density between today and redshift $z \sim 1100$. In this paper we develop a method for a direct reconstruction of the projected matter density from the CMB anisotropies. This reconstruction is obtained by averaging over quadratic combinations of the derivatives of CMB field. We test the method using simulations and show that it can successfully recover projected density profile of a cluster of galaxies if there are measurable anisotropies on scales smaller than the characteristic cluster size. In the absence of sufficient small scale power the reconstructed maps have low signal to noise on individual structures, but can give a positive detection of the power spectrum or when cross correlated with other maps of large scale structure. We develop an analytic method to reconstruct the power spectrum including the effects of noise and beam smoothing. Tests with Monte Carlo simulations show that we can recover the input power spectrum both on large and small scales, provided that we use maps with sufficiently low noise and high angular resolution.