Source author record

Matthias Bartelmann

Matthias Bartelmann appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

83works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

83 published item(s)

preprint2023arXiv

Baryon-photon interactions in Resummed Kinetic Field Theory

We explore how interactions between baryons and photons can be incorporated into Kinetic Field Theory (KFT), a description of cosmic structure formation based on classical Hamiltonian particle dynamics. In KFT, baryons are described as effective mesoscopic particles which represent fluid elements governed by the hydrodynamic equations. In this paper, we modify the mesoscopic particle model to include pressure effects exerted on baryonic matter through interactions with photons. As a proof of concept, we use this extended mesoscopic model to describe the tightly coupled baryon-photon fluid between matter-radiation equality and recombination. We show that this model can qualitatively reproduce the formation of baryon-acoustic oscillations in the cosmological power spectrum.

preprint2022arXiv

Kinetic Field Theory for Cosmic Structure Formation

We apply kinetic field theory to non-linear cosmic structure formation. Kinetic field theory decomposes the cosmic density field into particles and follows their trajectories through phase space. We assume that initial particle momenta are drawn from a Gaussian random field. We place particular emphasis on the late-time, asymptotic behaviour on small spatial scales of low-order statistical measures for the distribution of particles in configuration and velocity space. Our main result is that the power spectra for density and velocity fluctuations in ensembles of particles freely streaming along Zel'dovich trajectories asymptotically fall off with wave number $k$ like $k^{-3}$ for $k\to\infty$, irrespective of the cosmological model and the type of dark matter assumed, with the exponent set only by the number of spatial dimensions. This conclusion remains valid for density-fluctuation power spectra if particle interactions are taken into account in a mean-field approximation. We also show that the bispectrum of freely-streaming particles falls off asymptotically like $k^{-11/2}$ under the same general conditions.

preprint2022arXiv

Triaxiality in galaxy clusters: Mass versus Potential reconstructions

Accounting for the triaxial shapes of galaxy clusters will become important in the context of upcoming cosmological surveys. We show that, compared to the gas density distribution, the cluster gravitational potential can be better characterised by a simple 3D model and is more robust against fluctuations. Perturbations in the gas density distribution can have a substantial influence on the derived thermodynamic properties, while cluster potentials are smooth and well-approximated by a spheroidal model. We use a statistical sample of 85 galaxy clusters from a large cosmological hydrodynamical simulation to investigate cluster shapes as a function of radius. In particular, we examine the shape of isodensity and isopotential shells and analyze how it is affected by the choice of component (gas vs. potential), substructure removal (for the gas density) and the definition of the computation domain (interior vs. shells). We find that the orientation and axis ratios of gas isodensity contours are degenerate with the presence of substructures and unstable against fluctuations. We observe that, as the derived cluster shape depends on the method used for removing the substructures, thermodynamic properties extracted from e.g. X-ray emissivity profiles suffer from this additional, often underestimated bias. In contrast, the shape reconstruction of the potential is largely unaffected by these factors and converges towards simple geometric models for both relaxed and dynamically active clusters. The observation that cluster potentials are better represented by simple geometrical models and reconstructed with a low level of systematics for both dynamically active and relaxed clusters suggests that by characterising galaxy clusters by their potential rather than by their mass, dynamically active and relaxed clusters could be combined in cosmological studies, improving statistics and lowering scatter.

preprint2021arXiv

Gravitational wave induced baryon acoustic oscillations

We study the impact of gravitational waves originating from a first order phase transition on structure formation. To do so, we perform a second order perturbation analysis in the $1+3$ covariant framework and derive a wave equation in which second order, adiabatic density perturbations of the photon-baryon fluid are sourced by the gravitational wave energy density during radiation domination and on sub-horizon scales. The scale on which such waves affect the energy density perturbation spectrum is found to be proportional to the horizon size at the time of the phase transition times its inverse duration. Consequently, structure of the size of galaxies and bigger can only be affected in this way by relatively late phase transitions at $\ge 10^{6}\,\text{s}$. Using cosmic variance as a bound we derive limits on the strength $α$ and the relative duration $(β/H_*)^{-1}$ of phase transitions as functions of the time of their occurrence which results in a new exclusion region for the energy density in gravitational waves today. We find that the cosmic variance bound forbids only relative long lasting phase transitions, e.g. $β/H_*\lesssim 6.8$ for $t_*\approx 5\times10^{11}\,\text{s}$, which exhibit a substantial amount of supercooling $α>20$ to affect the matter power spectrum.

preprint2020arXiv

Analytic Derivation of the Halo Mass Function from the Non-Linear Cosmic Density Field

We estimate the halo mass function (HMF) by applying the excursion set approach to the non-linear cosmic density field. Thereby, we account for the non-Gaussianity of today's density distribution and constrain the HMF independent of the linear collapse threshold $δ_{\textrm{crit}}$. We consider a spherical region as a halo, if its density today exceeds the virial overdensity threshold $Δ$. We model the probability distribution of the non-linear density field by a superposition of a Gaussian and a lognormal distribution, which we constrain with the bispectrum of density fluctuations, predicted by the kinetic field theory description of cosmic structure formation. Two different excursion set approaches are compared. The first treats the density $δ$ as an uncorrelated random walk of the smoothing scale $R$. The second assumes $δ(R)$ to be correlated. We find that the resulting HMFs correspond well to the HMF found in numerical simulations if the correlation of $δ(R)$ is taken into account. Furthermore, the HMF depends only weakly on the choice of the density threshold $Δ$.

preprint2020arXiv

Model Independent Analysis of Supernova Data, Dark Energy, Trans-Planckian Censorship and the Swampland

In this Letter, we consider the model-independent reconstruction of the expansion and growth functions from the Pantheon supernova data. The method relies on developing the expansion function in terms of shifted Chebyshev polynomials and determining the coefficients of the polynomials by a maximum-likelihood fit to the data. Having obtained the expansion function in a model-independent way, we can then also determine the growth function without assuming a particular model. We then compare the results with the predictions of two classes of Dark Energy models, firstly a class of quintessence scalar field models consistent with the trans-Planckian censorship and swampland conjectures, and secondly a class of generalized Proca vector field models. We determine constraints on the parameters which appear in these models.

preprint2020arXiv

Resummed Kinetic Field Theory: a model of coupled baryonic and dark matter

We present a new analytical description of cosmic structure formation in a mixture of dark and baryonic matter, using the framework of Kinetic Field Theory (KFT) -- a statistical field theory for classical particle dynamics. So far, KFT has only been able to describe a single type of particles, sufficient to consider structure growth due to the gravitational interactions between dark matter. However, the influence of baryonic gas dynamics becomes increasingly relevant when describing smaller scales. In this paper, we thus demonstrate how to extend the KFT formalism as well as a previously presented resummation scheme towards describing such mixtures of two particle species. Thereby, the gas dynamics of baryons are accounted for using the recently developed model of Mesoscopic Particle Hydrodynamics. Assuming a flat $Λ$CDM Universe and a simplified model for the thermal gas evolution, we demonstrate the validity of this approach by computing the linear evolution of the individual and total matter power spectra between the epoch of recombination and today. Our results correctly reproduce the expected behaviour, showing a suppression of both baryonic and dark matter structure growth on scales smaller than the baryonic Jeans length, in good agreement with results from the numerical Boltzmann solver CLASS. Nonlinear corrections within this approach will be investigated in upcoming works.

preprint2019arXiv

Cosmic Structure Formation with Kinetic Field Theory

Kinetic Field Theory (KFT) is a statistical field theory for an ensemble of point-like classical particles in or out of equilibrium. We review its application to cosmological structure formation. Beginning with the construction of the generating functional of the theory, we describe in detail how the theory needs to be adapted to reflect the expanding spatial background and the homogeneous and isotropic, correlated initial conditions for cosmic structures. Based on the generating functional, we develop three main approaches to non-linear, late-time cosmic structures, which rest either on the Taylor expansion of an interaction operator, suitable averaging procedures for the interaction term, or a resummation of perturbation terms. We show how an analytic, parameter-free equation for the non-linear cosmic power spectrum can be derived. We explain how the theory can be used to derive the density profile of gravitationally bound structures and use it to derive power spectra of cosmic velocity densities. We further clarify how KFT relates to the BBGKY hierarchy. We then proceed to apply kinetic field theory to fluids, introduce a reformulation of KFT in terms of macroscopic quantities which leads to a resummation scheme, and use this to describe mixtures of gas and dark matter. We discuss how KFT can be applied to study cosmic structure formation with modified theories of gravity. As an example for an application to a non-cosmological particle ensemble, we show results on the spatial correlation function of cold Rydberg atoms derived from KFT.

preprint2016arXiv

Weak gravitational lensing

According to the theory of general relativity, masses deflect light in a way similar to convex glass lenses. This gravitational lensing effect is astigmatic, giving rise to image distortions. These distortions allow to quantify cosmic structures statistically on a broad range of scales, and to map the spatial distribution of dark and visible matter. We summarise the theory of weak gravitational lensing and review applications to galaxies, galaxy clusters and larger-scale structures in the Universe.

preprint2015arXiv

Constraints on $Ω_\mathrm{m}$ and $σ_8$ from the potential-based cluster temperature function

The abundance of galaxy clusters is in principle a powerful tool to constrain cosmological parameters, especially $Ω_\mathrm{m}$ and $σ_8$, due to the exponential dependence in the high-mass regime. While the best observables are the X-ray temperature and luminosity, the abundance of galaxy clusters, however, is conventionally predicted as a function of mass. Hence, the intrinsic scatter and the uncertainties in the scaling relations between mass and either temperature or luminosity lower the reliability of galaxy clusters to constrain cosmological parameters. In this article, we further refine the X-ray temperature function for galaxy clusters by Angrick et al., which is based on the statistics of perturbations in the cosmic gravitational potential and proposed to replace the classical mass-based temperature function, by including a refined analytic merger model and compare the theoretical prediction to results from a cosmological hydrodynamical simulation. Although we find already a good agreement if we compare with a cluster temperature function based on the mass-weighted temperature, including a redshift-dependent scaling between mass-based and spectroscopic temperature yields even better agreement between theoretical model and numerical results. As a proof of concept, incorporating this additional scaling in our model, we constrain the cosmological parameters $Ω_\mathrm{m}$ and $σ_8$ from an X-ray sample of galaxy clusters and tentatively find agreement with the recent cosmic microwave background based results from the Planck mission at 1$σ$-level.

preprint2015arXiv

Degeneracies of parametric lens model families near folds and cusps

We develop an approach to select families of lens models that can describe doubly and triply gravitationally lensed images near folds and cusps using the model-independent ratios of lensing-potential derivatives derived in Wagner & Bartelmann (2015). Models are selected by comparing these model-independent ratios of potential derivatives to (numerically determined) ratios of potential derivatives along critical curves for entire lens model families in a given range of parameter values. This comparison returns parameter ranges which lens model families can reproduce observation within, as well as sections of the critical curve where image sets of the observed type can appear. If the model-independent potential-derivative ratios inferred from the observation fall outside the range of these ratios derived for the lens model family, the entire family can be excluded as a feasible model in the given volume in parameter space. We employ this approach for the family of singular isothermal spheres with external shear to examples of lensing by a galaxy and two galaxy clusters (JVAS B1422+231, SDSS J2222+2745, and MACS J1149.5+2223) and show that the results obtained by our general method are in good agreement with results of previous model fits.

preprint2015arXiv

Evolution of linear perturbations in Lemaître-Tolman-Bondi void models

We study the evolution of linear perturbations in a Lemaître-Tolman-Bondi (LTB) void model with realistic cosmological initial conditions. Linear perturbation theory in LTB models is substantially more complicated than in standard Friedmann universes as the inhomogeneous background causes gauge-invariant perturbations to couple at first order. As shown by Clarkson et al. (2009), the evolution is constrained by a system of linear partial differential equations which need to be integrated numerically. We present a new numerical scheme using finite element methods to solve this equation system and generate scalar initial conditions based on Gaussian random fields with an underlying power spectrum for the Bardeen potential. After spherical harmonic decomposition, the initial fluctuations are propagated in time and estimates of angular power spectra of each gauge invariant variable are computed as functions of redshift. This allows to analyse the coupling strength in a statistical way. We find significant couplings up to $25\%$ for large and deep voids of Gpc scale as required to fit the distance redshift relations of SNe.

preprint2015arXiv

Hubble Space Telescope Combined Strong and Weak Lensing Analysis of the CLASH Sample: Mass and Magnification Models and Systematic Uncertainties

We present results from a comprehensive lensing analysis in HST data, of the complete CLASH cluster sample. We identify new multiple-images previously undiscovered allowing improved or first constraints on the cluster inner mass distributions and profiles. We combine these strong-lensing constraints with weak-lensing shape measurements within the HST FOV to jointly constrain the mass distributions. The analysis is performed in two different common parameterizations (one adopts light-traces-mass for both galaxies and dark matter while the other adopts an analytical, elliptical NFW form for the dark matter), to provide a better assessment of the underlying systematics - which is most important for deep, cluster-lensing surveys, especially when studying magnified high-redshift objects. We find that the typical (median), relative systematic differences throughout the central FOV are $\sim40\%$ in the (dimensionless) mass density, $κ$, and $\sim20\%$ in the magnification, $μ$. We show maps of these differences for each cluster, as well as the mass distributions, critical curves, and 2D integrated mass profiles. For the Einstein radii ($z_{s}=2$) we find that all typically agree within $10\%$ between the two models, and Einstein masses agree, typically, within $\sim15\%$. At larger radii, the total projected, 2D integrated mass profiles of the two models, within $r\sim2\arcmin$, differ by $\sim30\%$. Stacking the surface-density profiles of the sample from the two methods together, we obtain an average slope of $d\log (Σ)/d\log(r)\sim-0.64\pm0.1$, in the radial range [5,350] kpc. Lastly, we also characterize the behavior of the average magnification, surface density, and shear differences between the two models, as a function of both the radius from the center, and the best-fit values of these quantities.

preprint2015arXiv

Model-independent characterisation of strong gravitational lenses

We develop a new approach to extracting model-independent information from observations of strong gravitational lenses. The approach is based on the generic properties of images near the fold and cusp catastrophes in caustics and critical curves. Observables used are the relative image positions, the magnification ratios and ellipticities of extended images, and time delays between images with temporally varying intensity. We show how these observables constrain derivatives and ratios of derivatives of the lensing potential near a critical curve. Based on these measured properties of the lensing potential, classes of parametric lens models can then easily be restricted to such parameter values compatible with the measurements, thus allowing fast scans of large varieties of models. Applying our approach to a representative galaxy (JVAS B1422+231) and a galaxy-cluster lens (MACS J1149.5+2223), we show which model-independent information can be extracted in those cases and demonstrate that the parameters obtained by our approach for known parametric lens models agree well with those found by detailed model fitting.

preprint2015arXiv

Non-equilibrium statistical field theory for classical particles: Basic kinetic theory

Recently Mazenko and Das and Mazenko introduced a non-equilibrium field theoretical approach to describe the statistical properties of a classical particle ensemble starting from the microscopic equations of motion of each individual particle. We use this theory to investigate the transition from those microscopic degrees of freedom to the evolution equations of the macroscopic observables of the ensemble. For the free theory, we recover the continuity and Jeans equations of a collisionless gas. For a theory containing two-particle interactions in a canonical perturbation series, we find the macroscopic evolution equations to be described by the Born-Bogoliubov-Green-Kirkwood-Yvon hierarchy (BBGKY hierarchy) with a truncation criterion depending on the order in perturbation theory. This establishes a direct link between the classical and the field-theoretical approaches to kinetic theory that might serve as a starting point to investigate kinetic theory beyond the classical limits.

preprint2015arXiv

Reconstructing the projected gravitational potential of Abell 1689 from X-ray measurements

Context. Galaxy clusters can be used as cosmological probes, but to this end, they need to be thoroughly understood. Combining all cluster observables in a consistent way will help us to understand their global properties and their internal structure. Aims. We provide proof of the concept that the projected gravitational potential of galaxy clusters can directly be reconstructed from X-ray observations. We also show that this joint analysis can be used to locally test the validity of the equilibrium assumptions in galaxy clusters. Methods. We used a newly developed reconstruction method, based on Richardson-Lucy deprojection, that allows reconstructing projected gravitational potentials of galaxy clusters directly from X-ray observations. We applied this algorithm to the well-studied cluster Abell 1689 and compared the gravitational potential reconstructed from X-ray observables to the potential obtained from gravitational lensing measurements. [...] Results. Assuming spherical symmetry and hydrostatic equilibrium, the potentials recovered from gravitational lensing and from X-ray emission agree very well beyond 500 kpc. Owing to the fact that the Richardson-Lucy deprojection algorithm allows deprojecting each line of sight independently, this result may indicate that non-gravitational effects and/or asphericity are strong in the central regions of the clusters. Conclusions. We demonstrate the robustness of the potential reconstruction method based on the Richardson-Lucy deprojection algorithm and show that gravitational lensing and X-ray emission lead to consistent gravitational potentials. Our results illustrate the power of combining galaxy-cluster observables in a single, non-parametric, joint reconstruction of consistent cluster potentials that can be used to locally constrain the physical state of the gas.

preprint2015arXiv

The projected gravitational potential of the galaxy cluster MACS~J1206 derived from galaxy kinematics

We reconstruct the radial profile of the projected gravitational potential of the galaxy cluster MACS-J1206 from 592 spectroscopic measurements of velocities of cluster members. For doing so, we use a method we have developed recently based on the Richardson-Lucy deprojection algorithm and an inversion of the spherically-symmetric Jeans equation. We find that, within the uncertainties, our reconstruction agrees very well with a potential reconstruction from weak and strong gravitational lensing as well as with a potential obtained from X-ray measurements. In addition, our reconstruction is in good agreement with several common analytic profiles of the lensing potential. Varying the anisotropy parameter in the Jeans equation, we find that isotropy parameters which are either small, $β\lesssim0.2$, or decrease with radius yield potential profiles which strongly disagree with that obtained from gravitational lensing. We achieve the best agreement between our potential profile and the profile from gravitational lensing if the anisotropy parameter rises quite steeply to $β\approx0. 6$ within $\approx0.5\,\mathrm{Mpc}$ and stays constant further out.

preprint2014arXiv

CLASH-VLT: Constraints on the Dark Matter Equation of State from Accurate Measurements of Galaxy Cluster Mass Profiles

A pressureless scenario for the Dark Matter (DM) fluid is a widely adopted hypothesis, despite the absence of a direct observational evidence. According to general relativity, the total mass-energy content of a system shapes the gravitational potential well, but different test particles perceive this potential in different ways depending on their properties. Cluster galaxy velocities, being $\ll$c, depend solely on the gravitational potential, whereas photon trajectories reflect the contributions from the gravitational potential plus a relativistic-pressure term that depends on the cluster mass. We exploit this phenomenon to constrain the Equation of State (EoS) parameter of the fluid, primarily DM, contained in galaxy clusters. We use the complementary information provided by the kinematic and lensing mass profiles of the galaxy cluster MACS 1206.2-0847 at $z=0.44$, as obtained in an extensive imaging and spectroscopic campaign within the CLASH survey. The unprecedented high quality of our data-set and the properties of this cluster are well suited to determine the EoS parameter of the cluster fluid. Since baryons contribute at most $15\%$ to the total mass in clusters and their pressure is negligible, the EoS parameter we derive describes the behavior of the DM fluid. We obtain the most stringent constraint on the DM EoS parameter to date, $w=(p_r+2\,p_t)/(3\,c^2ρ)=0.00\pm0.15\mathrm{(stat)}\pm0.08\mathrm{(syst)}$, averaged over the radial range $0.5\,\mathrm{Mpc}\leq$$r$$\leq$$r_{200}$, where $p_r$ and $p_t$ are the radial and tangential pressure, and $ρ$ is the density. We plan to further improve our constraint by applying the same procedure to all clusters from the ongoing CLASH-VLT program.

preprint2014arXiv

CLASH: Extreme Emission Line Galaxies and Their Implication on Selection of High-Redshift Galaxies

We utilize the CLASH (Cluster Lensing And Supernova survey with Hubble) observations of 25 clusters to search for extreme emission-line galaxies (EELGs). The selections are carried out in two central bands: F105W (Y105) and F125W (J125), as the flux of the central bands could be enhanced by the presence of [O III] 4959, 5007 at redshift of about 0.93-1.14 and 1.57-1.79, respectively. The multi-band observations help to constrain the equivalent widths of emission lines. Thanks to cluster lensing, we are able to identify 52 candidates down to an intrinsic limiting magnitude of 28.5 and to a rest-frame [O III] 4959,5007 equivalent width of about 3737 angstrom. Our samples include a number of EELGs at lower luminosities that are missed in other surveys, and the extremely high equivalent width can be only found in such faint galaxies. These EELGs can mimic the dropout feature similar to that of high redshift galaxies and contaminate the color-color selection of high redshift galaxies when the S/N ratio is limited or the band coverage is incomplete. We predict that the fraction of EELGs in the future high redshift galaxy selections cannot be neglected.

preprint2014arXiv

Non-equilibrium statistical field theory for classical particles: Impact of correlated initial conditions on non-ideal gases

We use the non-equilibrium statistical field theory for classical particles recently developed by Mazenko and Das and Mazenko, together with the free generating functional for particles initially correlated in phase space derived in Bartelmann et al. to study the impact of initial correlations on the equation of state of real gases. We first show that we can reproduce the well known van der Waals equation of state for uncorrelated initial conditions using this approach. We then impose correlated initial conditions and study their qualitative and quantitative effect on the equation of state of a van der Waals gas. The correlations impose a significant correction to the pressure of an ideal gas which is an order of magnitude larger than the correction due to particle interactions.

preprint2014arXiv

Non-equilibrium statistical field theory for classical particles: Initially correlated grand canonical ensembles

It was recently shown by Bartelmann et al. how correlated initial conditions can be introduced into the statistical field theory for classical particles pioneered by Das and Mazenko. In this paper we extend this development from the canonical to the grand canonical ensemble for a system satisfying statistical homogeneity and isotropy. We do this by translating the probability distribution for the initial phase space coordinates of the particles into an easy diagrammatic representation and then using a variant of the Mayer cluster expansion to sum over particle numbers. The grand canonical generating functional is then used in a structured approach to the derivation of the non-interacting cumulants of the two core collective fields, the density $ρ$ and the response field $B$. As a side-product we find several theorems pertaining to these cumulants which will be useful when investigating the interacting regime of the theory in future work.

preprint2014arXiv

Non-equilibrium statistical field theory for classical particles: Linear and mildly non-linear evolution of cosmological density power spectra

We use the non-equlibrium statistical field theory for classical particles, recently developed by Mazenko and Das and Mazenko, together with the free generating functional we have previously derived for point sets initially correlated in phase space, to calculate the time evolution of power spectra in the free theory, i.e. neglecting particle interactions. We provide expressions taking linear and quadratic momentum correlations into account. Up to this point, the expressions are general with respect to the free propagator of the microscopic degrees of freedom. We then specialise the propagator to that expected for particles in cosmology treated within the Zel'dovich approximation and show that, to linear order in the momentum correlations, the linear growth of the cosmological power spectrum is reproduced. Quadratic momentum correlations return a first contribution to the non-linear evolution of the power spectrum, for which we derive a simple closed expression valid for arbitrary wave numbers. This expression is a convolution of the initial density power spectrum with itself, multiplied by a mode-coupling kernel. We also derive the bispectrum expected in this theory within these approximations and show that its connected part reproduces almost, but not quite, the bispectrum expected in Eulerian perturbation theory of the density contrast.

preprint2014arXiv

Non-equilibrium statistical field theory for classical particles: Non-linear structure evolution with first-order interaction

We calculate the power spectrum of density fluctuations in the statistical non-equilibrium field theory for classical, microscopic degrees of freedom to first order in the interaction potential. We specialise our result to cosmology by choosing appropriate initial conditions and propagators and show that the non-linear growth of the density power spectrum found in numerical simulations of cosmic structure evolution is reproduced well to redshift zero and for arbitrary wave numbers. The main difference of our approach to ordinary cosmological perturbation theory is that we do not perturb a dynamical equation for the density contrast. Rather, we transport the initial phase-space distribution of a canonical particle ensemble forward in time and extract any collective information from it at the time needed. Since even small perturbations of particle trajectories can lead to large fluctuations in density, our approach allows to reach high density contrast already at first order in the perturbations of the particle trajectories. We argue why the expected asymptotic behaviour of the non-linear power spectrum at large wave numbers can be reproduced in our approach at any order of the perturbation series.

preprint2014arXiv

Probing spatial homogeneity with LTB models: a detailed discussion

Do current observational data confirm the assumptions of the cosmological principle, or is there statistical evidence for deviations from spatial homogeneity on large scales? To address these questions, we developed a flexible framework based on spherically symmetric, but radially inhomogeneous Lemaitre-Tolman-Bondi (LTB) models with synchronous Big Bang. We expanded the (local) matter density profile in terms of flexible interpolation schemes and orthonormal polynomials. A Monte Carlo technique in combination with recent observational data was used to systematically vary the shape of these profiles. In the first part of this article, we reconsider giant LTB voids without dark energy to investigate whether extremely fine-tuned mass profiles can reconcile these models with current data. While the local Hubble rate and supernovae can easily be fitted without dark energy, however, model-independent constraints from the Planck 2013 data require an unrealistically low local Hubble rate, which is strongly inconsistent with the observed value; this result agrees well with previous studies. In the second part, we explain why it seems natural to extend our framework by a non-zero cosmological constant, which then allows us to perform general tests of the cosmological principle. Moreover, these extended models facilitate explorating whether fluctuations in the local matter density profile might potentially alleviate the tension between local and global measurements of the Hubble rate, as derived from Cepheid-calibrated type Ia supernovae and CMB experiments, respectively. We show that current data provide no evidence for deviations from spatial homogeneity on large scales. More accurate constraints are required to ultimately confirm the validity of the cosmological principle, however.

preprint2014arXiv

The strongest gravitational lenses: III. The order statistics of the largest Einstein radii

The Einstein radius (ER) of a gravitational lens encodes information about decisive quantities such as halo mass, concentration, triaxiality, and orientation with respect to the observer. Thus, the largest Einstein radii can potentially be utilised to test the predictions of the LCDM model. Hitherto, studies have focussed on the single largest observed ER. We extend those studies by employing order statistics to formulate exclusion criteria based on the n largest Einstein radii and apply these criteria to the strong lensing analysis of 12 MACS clusters at z>0.5. We obtain the order statistics of Einstein radii by a MC approach, based on the semi-analytic modelling of the halo population on the past lightcone. After sampling the order statistics, we fit a GEV distribution to the first-order distribution, which allows us to derive analytic relations for the order statistics of the Einstein radii. We find that the Einstein radii of the 12 MACS clusters are not in conflict with the LCDM expectations. Our exclusion criteria indicate that, in order to exhibit tension with the concordance model, one would need to observe approximately twenty Einstein radii >30", ten >35" or five >42" in the range of 0.5<z<1.0 on the full sky. Furthermore, we find that, with increasing order, the haloes with the largest Einstein radii are on average less aligned along the line-of-sight and less triaxial. In general, the cumulative distribution functions steepen for higher orders, giving them better constraining power. (abridged)

preprint2014arXiv

The strongest gravitational lenses: IV. The order statistics of the largest Einstein radii with cluster mergers

Based on techniques developed in the previous papers of this series, we investigate the impact of galaxy-cluster mergers on the order statistics of the largest Einstein radii. We show that the inclusion of mergers significantly shifts the extreme value distribution of the largest Einstein radius to higher values, typically increasing the expected value by ${\sim}10\%$. A comparison with current data reveals that the largest observed Einstein radius agrees excellently well with the theoretical predictions of the $Λ$CDM model at redshifts $z > 0.5$. At redshifts $z < 0.5$, our results are somewhat more controversial. Although cluster mergers also increase the expected values of the order statistics of the $n$ largest Einstein radii by ${\sim}10\%$, the theoretically expected values are notably lower (${\sim}3σ$ deviation for $n = 12$) than the largest Einstein radii of a selected sample of SDSS clusters in the redshift range $0.1 \leq z \leq 0.55$. The uncertainties of the observed Einstein radii are still large, however, and thus the measurements need to be carefully revised in future works. Therefore, given the premature state of current observational data, overall, there is still no reliable statistical evidence for observed Einstein radii to exceed the theoretical expectations of the standard cosmological model.

preprint2014arXiv

Three Gravitationally Lensed Supernovae Behind CLASH Galaxy Clusters

We report observations of three gravitationally lensed supernovae (SNe) in the Cluster Lensing And Supernova survey with Hubble (CLASH) Multi-Cycle Treasury program. These objects, SN CLO12Car (z = 1.28), SN CLN12Did (z = 0.85), and SN CLA11Tib (z = 1.14), are located behind three different clusters, MACSJ1720.2+3536 (z = 0.391), RXJ1532.9+3021 (z = 0.345), and Abell 383 (z = 0.187), respectively. Each SN was detected in Hubble Space Telescope (HST) optical and infrared images. Based on photometric classification, we find that SNe CLO12Car and CLN12Did are likely to be Type Ia supernovae (SNe Ia), while the classification of SN CLA11Tib is inconclusive. Using multi-color light-curve fits to determine a standardized SN Ia luminosity distance, we infer that SN CLO12Car was approximately 1.0 +/- 0.2 mag brighter than field SNe Ia at a similar redshift and ascribe this to gravitational lens magnification. Similarly, SN CLN12Did is approximately 0.2 +/- 0.2 mag brighter than field SNe Ia. We derive independent estimates of the predicted magnification from CLASH strong+weak lensing maps of the clusters: 0.83 +/- 0.16 mag for SN CLO12Car, 0.28 +/- 0.08 mag for SN CLN12Did, and 0.43 +/- 0.11 mag for SN CLA11Tib. The two SNe Ia provide a new test of the cluster lens model predictions: we find that the magnifications based on the SN Ia brightness and those predicted by the lens maps are consistent. Our results herald the promise of future observations of samples of cluster-lensed SNe Ia (from the ground or space) to help illuminate the dark-matter distribution in clusters of galaxies, through the direct determination of absolute magnifications.

preprint2014arXiv

Trajectories of point particles in cosmology and the Zel'dovich approximation

Using a Green's function approach, we compare the trajectories of classical Hamiltonian point particles in an expanding space-time to the effectively inertial trajectories in the Zel'dovich approximation. It is shown that the effective gravitational potential accelerating the particles relative to the Zel'dovich trajectories vanishes exactly initially as a consequence of the continuity equation, and acts only during a short, early period. The Green's function approach suggests an iterative scheme for improving the Zel'dovich trajectories, which can be analytically solved. We construct these trajectories explicitly and show how they interpolate between the Zel'dovich and the exact trajectories. The effective gravitational potential acting on the improved trajectories is substantially smaller at late times than the potential acting on the exact trajectories. The results may be useful for Lagrangian perturbation theory and for numerical simulations.

preprint2013arXiv

A comparison of structure formation in minimally and non-minimally coupled quintessence models

We study structure formation in non-minimally coupled dark energy models, where there is a coupling in the Lagrangian between a quintessence scalar field and gravity via the Ricci scalar. We consider models with a range of different non-minimal coupling strengths and compare these to minimally coupled quintessence models with time-dependent dark energy densities. The equations of state of the latter are tuned to either reproduce the equation of state of the non-minimally coupled models or their background history. Thereby they provide a reference to study the unique imprints of coupling on structure formation. We show that the coupling between gravity and the scalar field, which effectively results in a time-varying gravitational constant G, is not negligible and its effect can be distinguished from a minimally coupled model. We extend previous work on this subject by showing that major differences appear in the determination of the mass function at high masses, where we observe differences of the order of 40% at z=0. Our new results concern effects on the non-linear matter power spectrum and on the lensing signal (differences of ~10% for both quantities), where we find that non-minimally coupled models could be distinguished from minimally coupled ones.

preprint2013arXiv

A study of relative velocity statistics in Lagrangian perturbation theory with PINOCCHIO

Subject of this paper is a detailed analysis of the PINOCCHIO algorithm for studying the relative velocity statistics of merging haloes in Lagrangian perturbation theory. Given a cosmological background model, a power spectrum of fluctuations as well as a Gaussian linear density contrast field $δ_{\rm l}$ is generated on a cubic grid, which is then smoothed repeatedly with Gaussian filters. For each Lagrangian particle at position $\bmath{q}$ and each smoothing radius $R$, the collapse time, the velocities and ellipsoidal truncation are computed using Lagrangian Perturbation Theory. The collapsed medium is then fragmented into isolated objects by an algorithm designed to mimic the accretion and merger events of hierarchical collapse. Directly after the fragmentation process the mass function, merger histories of haloes and the statistics of the relative velocities at merging are evaluated. We reimplemented the algorithm in C++, recovered the mass function and optimised the construction of halo merging histories. Comparing our results with the output of the Millennium simulation suggests that PINOCCHIO is well suited for studying relative velocities of merging haloes and is able to reproduce the pairwise velocity distribution.

preprint2013arXiv

CLASH: Complete Lensing Analysis of the Largest Cosmic Lens MACS J0717.5+3745 and Surrounding Structures

The galaxy cluster MACS J0717.5+3745 (z=0.55) is the largest known cosmic lens, with complex internal structures seen in deep X-ray, Sunyaev-Zel'dovich effect and dynamical observations. We perform a combined weak and strong lensing analysis with wide-field BVRi'z' Subaru/Suprime-Cam observations and 16-band Hubble Space Telescope observations taken as part of the Cluster Lensing And Supernova survey with Hubble (CLASH). We find consistent weak distortion and magnification measurements of background galaxies, and combine these signals to construct an optimally estimated radial mass profile of the cluster and its surrounding large-scale structure out to 5 Mpc/h. We find consistency between strong-lensing and weak-lensing in the region where these independent data overlap, <500 kpc/h. The two-dimensional weak-lensing map reveals a clear filamentary structure traced by distinct mass halos. We model the lensing shear field with 9 halos, including the main cluster, corresponding to mass peaks detected above 2.5σ_κ. The total mass of the cluster as determined by the different methods is M_{vir}=(2.8\pm0.4) \times 10^15 M_sun. Although this is the most massive cluster known at z>0.5, in terms of extreme value statistics we conclude that the mass of MACS J0717.5+3745 by itself is not in serious tension with LambdaCDM, representing only a ~2σ departure above the maximum simulated halo mass at this redshift.

preprint2013arXiv

Galaxy halo truncation and Giant Arc Surface Brightness Reconstruction in the Cluster MACSJ1206.2-0847

In this work we analyze the mass distribution of MACSJ1206.2-0847, especially focusing on the halo properties of its cluster members. The cluster appears relaxed in its X-ray emission, but has significant amounts of intracluster light which is not centrally concentrated, suggesting that galaxy-scale interactions are still ongoing despite the overall relaxed state. The cluster lenses 12 background galaxies into multiple images and one galaxy at $z=1.033$ into a giant arc and its counterimage. The multiple image positions and the surface brightness distribution (SFB) of the arc which is bent around several cluster members are sensitive to the cluster galaxy halo properties. We model the cluster mass distribution with a NFW profile and the galaxy halos with two parameters for the mass normalization and extent of a reference halo assuming scalings with their observed NIR--light. We match the multiple image positions at an r.m.s. level of $0.85\arcsec$ and can reconstruct the SFB distribution of the arc in several filters to a remarkable accuracy based on this cluster model. The length scale where the enclosed galaxy halo mass is best constrained is about 5 effective radii -- a scale in between those accessible to dynamical and field strong lensing mass estimates on one hand and galaxy--galaxy weak lensing results on the other hand. The velocity dispersion and halo size of a galaxy with $m_{\rm 160W,AB}=19.2$ or $M_{\rm B,Vega}=-20.7$ are $σ=150 \rm kms^{-1}$ and $r\approx 26\pm 6 \rm kpc$, indicating that the halos of the cluster galaxies are tidally stripped. We also reconstruct the unlensed source (which is smaller by a factor of $\sim5.8$ in area), demonstrating the increase of morphological information due to lensing and conclude that this galaxy has likely star--forming spiral arms with a red (older) central component.

preprint2013arXiv

Internal Cluster Structure

The core structure of galaxy clusters is fundamentally important. Even though self-gravitating systems have no stable equilibrium state due to their negative heat capacity, numerical simulations find density profiles which are universal in the sense that they are fairly flat within a scale radius and gradually steepen farther outward, asymptotically approaching a logarithmic slope of $\approx-3$ near the virial radius. We argue that the reason for the formation of this profile is not satisfactorily understood. The ratio between the virial radius and the scale radius, the so-called concentration, is found in simulations to be closely related to the mass and the redshift and low for cluster-sized haloes, but observed to be substantially higher at least in a subset of observed clusters. Haloes formed from cold dark matter should furthermore be richly substructured. We review theoretical and observational aspects of cluster cores here, discuss modifications by baryonic physics and observables that can provide better insight into the internal structure of clusters.

preprint2013arXiv

Joint reconstruction of galaxy clusters from gravitational lensing and thermal gas I. Outline of a non-parametric method

We present a method to estimate the lensing potential from massive galaxy clusters for given observational X-ray data. The concepts developed and applied in this work can easily be combined with other techniques to infer the lensing potential, e.g. weak gravitational lensing or galaxy kinematics, to obtain an overall best fit model for the lensing potential. After elaborating on the physical details and assumptions the method is based on, we explain how the numerical algorithm itself is implemented with a Richardson-Lucy algorithm as a central part. Our reconstruction method is tested on simulated galaxy clusters with a spherically symmetric NFW density profile filled with gas in hydrostatic equilibrium. We describe in detail how these simulated observational data sets are created and how they need to be fed into our algorithm. We test the robustness of the algorithm against small parameter changes and estimate the quality of the reconstructed lensing potentials. As it turns out we achieve a very high degree of accuracy in reconstructing the lensing potential. The statistical errors remain below 2.0% whereas the systematical error does not exceed 1.0%.

preprint2013arXiv

Joint reconstruction of galaxy clusters from gravitational lensing and thermal gas. II. Inversion of the thermal Sunyaev-Zel'dovich effect

This paper continues a series in which we intend to show how all observables of galaxy clusters can be combined to recover the two-dimensional, projected gravitational potential of individual clusters. Our goal is to develop a non-parametric algorithm for joint cluster reconstruction taking all cluster observables into account. In this paper, we begin with the relation between the Compton-y parameter and the Newtonian gravitational potential, assuming hydrostatic equilibrium and a polytropic stratification of the intracluster gas. We show how Richardson-Lucy deconvolution can be used to convert the intensity change of the CMB due to the thermal Sunyaev-Zel'dovich effect into an estimate for the two-dimensional gravitational potential. Synthetic data simulated with characteristics of the ALMA telescope show that the two-dimensional potential of a cluster with mass 5*10^14 M_sun/h at redshift 0.2 is possible with an error of < 5% between the cluster centre and a radius r < 0.9 Mpc/h.

preprint2013arXiv

Masses of galaxy clusters from gravitational lensing

Despite consistent progress in numerical simulations, the observable properties of galaxy clusters are difficult to predict ab initio. It is therefore important to compare both theoretical and observational results to a direct measure of the cluster mass. This can be done by measuring the gravitational lensing effects caused by the bending of light by the cluster mass distribution. In this review we discuss how this phenomenon can be used to determine cluster masses and study the mass distribution itself. As sample sizes increase, the accuracy of the weak lensing mass estimates needs to improve accordingly. We discuss the main practical aspects of these measurements. We review a number of applications and highlight some recent results.

preprint2013arXiv

Reconstructing the projected gravitational potential of galaxy clusters from galaxy kinematics

We develop a method for reconstructing the two-dimensional, projected gravitational potential of galaxy clusters from observed line-of-sight velocity dispersions of cluster galaxies. It is the third of an intended series of papers aiming at a unique reconstruction method for cluster potentials combining lensing, X-ray, Sunyaev-Zel'dovich and kinematic data. The observed galaxy velocity dispersions are deprojected using the Richardson-Lucy algorithm. The obtained radial velocity dispersions are then related to the gravitational potential by using the tested assumption of a polytropic relation between the effective galaxy pressure and the density. Once the gravitational potential is obtained in three dimensions, projection along the line-of-sight yields the two-dimensional potential. For simplicity we adopt spherical symmetry and a known profile for the anisotropy parameter of the galaxy velocity dispersions. We test the method with a numerically simulated galaxy cluster and galaxies identified therein. We extract a projected velocity-dispersion profile from the simulated cluster and pass it through our algorithm, showing that the deviation between the true and the reconstructed gravitational potential is less then 10% within approximately 1.2 Mpc/h from the cluster centre.

preprint2013arXiv

The Contribution of Halos with Different Mass Ratios to the Overall Growth of Cluster-Sized Halos

We provide a new observational test for a key prediction of the ΛCDM cosmological model: the contributions of mergers with different halo-to-main-cluster mass ratios to cluster-sized halo growth. We perform this test by dynamically analyzing seven galaxy clusters, spanning the redshift range $0.13 < z_c < 0.45$ and caustic mass range $0.4-1.5$ $10^{15} h_{0.73}^{-1}$ M$_{\odot}$, with an average of 293 spectroscopically-confirmed bound galaxies to each cluster. The large radial coverage (a few virial radii), which covers the whole infall region, with a high number of spectroscopically identified galaxies enables this new study. For each cluster, we identify bound galaxies. Out of these galaxies, we identify infalling and accreted halos and estimate their masses and their dynamical states. Using the estimated masses, we derive the contribution of different mass ratios to cluster-sized halo growth. For mass ratios between ~0.2 and ~0.7, we find a ~1 $σ$ agreement with ΛCDM expectations based on the Millennium simulations I and II. At low mass ratios, $\lesssim 0.2$, our derived contribution is underestimated since the detection efficiency decreases at low masses, $\sim 2 \times 10^{14}$ $h_{0.73}^{-1}$ M$_{\odot}$. At large mass ratios, $\gtrsim 0.7$, we do not detect halos probably because our sample, which was chosen to be quite X-ray relaxed, is biased against large mass ratios. Therefore, at large mass ratios, the derived contribution is also underestimated.

preprint2012arXiv

A highly magnified candidate for a young galaxy seen when the Universe was 500 Myrs old

The early Universe at redshift z\sim6-11 marks the reionization of the intergalactic medium, following the formation of the first generation of stars. However, those young galaxies at a cosmic age of \lesssim 500 million years (Myr, at z \gtrsim 10) remain largely unexplored as they are at or beyond the sensitivity limits of current large telescopes. Gravitational lensing by galaxy clusters enables the detection of high-redshift galaxies that are fainter than what otherwise could be found in the deepest images of the sky. We report the discovery of an object found in the multi-band observations of the cluster MACS1149+22 that has a high probability of being a gravitationally magnified object from the early universe. The object is firmly detected (12 sigma) in the two reddest bands of HST/WFC3, and not detected below 1.2 μm, matching the characteristics of z\sim9 objects. We derive a robust photometric redshift of z = 9.6 \pm 0.2, corresponding to a cosmic age of 490 \pm 15Myr (i.e., 3.6% of the age of the Universe). The large number of bands used to derive the redshift estimate make it one of the most accurate estimates ever obtained for such a distant object. The significant magnification by cluster lensing (a factor of \sim15) allows us to analyze the object's ultra-violet and optical luminosity in its rest-frame, thus enabling us to constrain on its stellar mass, star-formation rate and age. If the galaxy is indeed at such a large redshift, then its age is less than 200 Myr (at the 95% confidence level), implying a formation redshift of zf \lesssim 14. The object is the first z>9 candidate that is bright enough for detailed spectroscopic studies with JWST, demonstrating the unique potential of galaxy cluster fields for finding highly magnified, intrinsically faint galaxies at the highest redshifts.

preprint2012arXiv

CLASH: Mass Distribution in and around MACS J1206.2-0847 from a Full Cluster Lensing Analysis

We derive an accurate mass distribution of the galaxy cluster MACS J1206.2-0847 (z=0.439) from a combined weak-lensing distortion, magnification, and strong-lensing analysis of wide-field Subaru BVRIz' imaging and our recent 16-band Hubble Space Telescope observations taken as part of the Cluster Lensing And Supernova survey with Hubble (CLASH) program. We find good agreement in the regions of overlap between several weak and strong lensing mass reconstructions using a wide variety of modeling methods, ensuring consistency. The Subaru data reveal the presence of a surrounding large scale structure with the major axis running approximately north-west south-east (NW-SE), aligned with the cluster and its brightest galaxy shapes, showing elongation with a \sim 2:1 axis ratio in the plane of the sky. Our full-lensing mass profile exhibits a shallow profile slope dlnΣ/dlnR\sim -1 at cluster outskirts (R>1Mpc/h), whereas the mass distribution excluding the NW-SE excess regions steepens further out, well described by the Navarro-Frenk-White form. Assuming a spherical halo, we obtain a virial mass M_{vir}=(1.1\pm 0.2\pm 0.1)\times 10^{15} M_{sun}/h and a halo concentration c_{vir} = 6.9\pm 1.0\pm 1.2 (\sim 5.7 when the central 50kpc/h is excluded), which falls in the range 4< <c> <7 of average c(M,z) predictions for relaxed clusters from recent Lambda cold dark matter simulations. Our full lensing results are found to be in agreement with X-ray mass measurements where the data overlap, and when combined with Chandra gas mass measurements, yield a cumulative gas mass fraction of 13.7^{+4.5}_{-3.0}% at 0.7Mpc/h (\approx 1.7r_{2500}), a typical value observed for high mass clusters.

preprint2012arXiv

CLASH: Precise New Constraints on the Mass Profile of Abell 2261

We precisely constrain the inner mass profile of Abell 2261 (z=0.225) for the first time and determine this cluster is not "over-concentrated" as found previously, implying a formation time in agreement with ΛCDM expectations. These results are based on strong lensing analyses of new 16-band HST imaging obtained as part of the Cluster Lensing and Supernova survey with Hubble (CLASH). Combining this with revised weak lensing analyses of Subaru wide field imaging with 5-band Subaru + KPNO photometry, we place tight new constraints on the halo virial mass M_vir = 2.2\pm0.2\times10^15 M\odot/h70 (within r \approx 3 Mpc/h70) and concentration c = 6.2 \pm 0.3 when assuming a spherical halo. This agrees broadly with average c(M,z) predictions from recent ΛCDM simulations which span 5 <~ <c> <~ 8. Our most significant systematic uncertainty is halo elongation along the line of sight. To estimate this, we also derive a mass profile based on archival Chandra X-ray observations and find it to be ~35% lower than our lensing-derived profile at r2500 ~ 600 kpc. Agreement can be achieved by a halo elongated with a ~2:1 axis ratio along our line of sight. For this elongated halo model, we find M_vir = 1.7\pm0.2\times10^15 M\odot/h70 and c_vir = 4.6\pm0.2, placing rough lower limits on these values. The need for halo elongation can be partially obviated by non-thermal pressure support and, perhaps entirely, by systematic errors in the X-ray mass measurements. We estimate the effect of background structures based on MMT/Hectospec spectroscopic redshifts and find these tend to lower Mvir further by ~7% and increase cvir by ~5%.

preprint2012arXiv

Cluster Lensing And Supernova survey with Hubble (CLASH): An Overview

The Cluster Lensing And Supernova survey with Hubble (CLASH) is a 524-orbit multi-cycle treasury program to use the gravitational lensing properties of 25 galaxy clusters to accurately constrain their mass distributions. The survey, described in detail in this paper, will definitively establish the degree of concentration of dark matter in the cluster cores, a key prediction of CDM. The CLASH cluster sample is larger and less biased than current samples of space-based imaging studies of clusters to similar depth, as we have minimized lensing-based selection that favors systems with overly dense cores. Specifically, twenty CLASH clusters are solely X-ray selected. The X-ray selected clusters are massive (kT > 5 keV; 5 - 30 x 10^14 M_solar) and, in most cases, dynamically relaxed. Five additional clusters are included for their lensing strength (Einstein radii > 35 arcsec at z_source = 2) to further quantify the lensing bias on concentration, to yield high resolution dark matter maps, and to optimize the likelihood of finding highly magnified high-redshift (z > 7) galaxies. The high magnification, in some cases, provides angular resolutions unobtainable with any current UVOIR facility and can yield z > 7 candidates bright enough for spectroscopic follow-up. A total of 16 broadband filters, spanning the near-UV to near-IR, are employed for each 20-orbit campaign on each cluster. These data are used to measure precise (sigma_phz < 0.02(1+z)) photometric redshifts for dozens of newly discovered multiply-lensed images per cluster. Observations of each cluster are spread over 8 epochs to enable a search, primarily in the parallel fields, for Type Ia supernovae at z > 1 to improve constraints on the time dependence of the dark energy equation of state and the evolution of such supernovae in an epoch when the universe is matter dominated.

preprint2012arXiv

Constraints on Cosmological Parameters from Strong Gravitational Lensing by Galaxy Clusters

We investigate how observations of strong lensing can be used to infer cosmological parameters, in particular the equation of state of dark energy. We focus on the growth of the critical lines of lensing clusters with the source redshift as this behaviour depends on the distance-redshift relation and is therefore cosmologically sensitive. Purely analytical approaches are generally insufficient because they rely on axisymmetric mass distributions and thus cannot take irregular critical curves into account. We devise a numerical method based on the Metropolis-Hastings algorithm: an elliptical generalization of the NFW density profile is used to fit a lens model to an observed configuration of giant luminous arcs while simultaneously optimizing the geometry. A semi-analytic method, which derives geometric parameters from critical points, is discussed as a faster alternative. We test the approaches on mock observations of gravitational lensing by a numerically simulated cluster. We find that no constraints can be derived from observations of individual clusters if no knowledge of the underlying mass distribution is assumed. Uncertainties are improved if a fixed lens model is used for a purely geometrical optimization, but the choice of a parametric model may produce strong biases.

preprint2012arXiv

Estimation of halo ellipticity using spin-3 flexion

Estimating the ellipticity of dark matter haloes at the galaxy or galaxy-cluster scale can provide important constraints on the formation of galaxies or clusters, as well as on the nature of dark matter. We show in this paper that the spin-3 gravitational flexion can add useful information on the ellipticity of lensing haloes. We introduce a general formalism to decompose fields with arbitrary spin into radial and tangential components. The ratio of the tangential and radial flexion components directly estimates the lens ellipticity. We point out that any centroid offset will significantly bias our estimate, which on the other hand can be used to determine the centre of the lens halo.

preprint2012arXiv

Miscentring in Galaxy Clusters: Dark Matter to Brightest Cluster Galaxy Offsets in 10,000 SDSS Clusters

We characterise the typical offset between the Dark Matter (DM) projected centre and the Brightest Cluster Galaxy (BCG) in 10,000 SDSS clusters. To place constraints on the centre of DM, we use an automated strong-lensing analysis, mass-modelling technique which is based on the well-tested assumption that light traces mass. The cluster galaxies are modelled with a steep power-law, and the DM component is obtained by smoothing the galaxy distribution fitting a low-order 2D polynomial (via spline interpolation), while probing a whole range of polynomial degrees and galaxy power laws. We find that the offsets between the BCG and the peak of the smoothed light map representing the DM, Δ, are distributed equally around zero with no preferred direction, and are well described by a log-normal distribution with <log_{10}(Δ[h^{-1} Mpc])>=-1.895^{+0.003}_{-0.004}, and σ=0.501\pm0.004 (95% confidence levels), or <log_{10}(Δ[\arcsec])>=0.564\pm0.005, and σ=0.475\pm0.007. Some of the offsets originate in prior misidentifications of the BCG or other bright cluster members by the cluster finding algorithm, whose level we make an additional effort to assess, finding that ~10% of the clusters in the probed catalogue are likely to be misidentified, contributing to higher-end offsets in general agreement with previous studies. Our results constitute the first statistically-significant high-resolution distributions of DM-to-BCG offsets obtained in an observational analysis, and importantly show that there exists such a typical non-zero offset in the probed catalogue. The offsets show a weak positive correlation with redshift, so that higher separations are generally found for higher-z clusters in agreement with the hierarchical growth of structure, which in turn could help characterise the merger, relaxation and evolution history of clusters, in future studies. [ABRIDGED]

preprint2012arXiv

MOKA: a new tool for Strong Lensing Studies

Strong gravitational lensing is a powerful tool for probing the matter distribution in the cores of massive dark matter haloes. Recent and ongoing analyses of galaxy cluster surveys (MACS, CFHTLS, SDSS, SGAS, CLASH, LoCuSS) have adressed the question of the nature of the dark matter distribution in clusters. N-body simulations of cold dark-matter haloes consistently find that haloes should be characterized by a concentration-mass relation that decreases monotonically with halo mass, and populated by a large amount of substructures, representing the cores of accreted progenitor halos. It is important for our understanding of dark matter to test these predictions. We present MOKA, a new algorithm for simulating the gravitational lensing signal from cluster-sized haloes. It implements the most recent results from numerical simulations to create realistic cluster-scale lenses with properties independent of numerical resolution. We perform systematic studies of the strong lensing cross section as a function of halo structures. We find that the strong lensing cross sections depend most strongly on the concentration and on the inner slope of the density profile of a halo, followed in order of importance by halo triaxiality and the presence of a bright central galaxy.

preprint2012arXiv

Order statistics applied to the most massive and most distant galaxy clusters

In this work we present for the first time an analytic framework for calculating the individual and joint distributions of the n-th most massive or n-th highest redshift galaxy cluster for a given survey characteristic allowing to formulate LCDM exclusion criteria. We show that the cumulative distribution functions steepen with increasing order, giving them a higher constraining power with respect to the extreme value statistics. Additionally, we find that the order statistics in mass (being dominated by clusters at lower redshifts) is sensitive to the matter density and the normalisation of the matter fluctuations, whereas the order statistics in redshift is particularly sensitive to the geometric evolution of the Universe. For a fixed cosmology, both order statistics are efficient probes of the functional shape of the mass function at the high mass end. To allow a quick assessment of both order statistics, we provide fits as a function of the survey area that allow percentile estimation with an accuracy better than two per cent. Furthermore, we discuss the joint distributions in the two-dimensional case for different combinations of order. Having introduced the theory, we apply the order statistical analysis to the SPT massive cluster sample and MCXC catalogue and find that the ten most massive clusters in the sample are consistent with LCDM and the Tinker mass function. In turn, by assuming the LCDM reference cosmology, order statistics can also be utilised for consistency checks of the completeness of the observed sample and of the modelling of the survey selection function. [abridged]

preprint2012arXiv

Perturbation Theory Trispectrum in the Time Renormalisation Approach

An accurate theoretical description of structure formation at least in the mildly non-linear regime is essential for comparison with data from next generation galaxy surveys. In a recent approach one follows the time evolution of correlators directly and finds a hierarchy of evolution equations with increasing order (Pietroni 2008). So far, in this so called time renormalisation group method the trispectrum was neglected in order to obtain a closed set of equations. In this work we study the influence of the trispectrum on the evolution of the power spectrum. In order to keep the numerical cost at a manageable level we use the tree-level trispectrum from Eulerian perturbation theory. In comparison to numerical simulations we find improvement in the mildly non-linear regime up to k = 0.25 h/Mpc. Beyond k = 0.25 h/Mpc the perturbative description of the trispectrum fails and the method performs worse than without the trispectrum included. Our results reinforce the conceptual advantage of the time renormalisation group method with respect to perturbation theory.

preprint2012arXiv

Probing ionizing radiation of L<~0.1L* star-forming galaxies at z>~3 with strong lensing

We show the effectiveness of strong lensing in the characterisation of Lyman continuum emission from faint L<~0.1L* star-forming galaxies at redshift >~ 3. Past observations of L>~L* galaxies at redshift >~3 have provided upper limits of the average escape fraction of ionising radiation of fesc~5%. Galaxies with relatively high fesc (>10%) seem to be particularly rare at these luminosities, there is therefore the need to explore fainter limits. Before the advent of giant ground based telescopes, one viable way to probe fesc down to 0.05-0.15L* is to exploit strong lensing magnification. This is investigated with Monte Carlo simulations that take into account the current observational capabilities. Adopting a lensing cross-section of 10 arcmin^2 within which the magnification is higher than 1 (achievable with about 4-5 galaxy clusters), with a U-band survey depth of 30(30.5) (AB, 1-sigma), it is possible to constrain fesc for z~3 star-forming galaxies down to 15(10)% at 3-sigma for L<0.15L* luminosities. This is particularly interesting if fesc increases at fainter luminosities, as predicted from various HI reionization scenarios and radiation transfer modelling. Ongoing observational programs on galaxy clusters are discussed and offer positive prospects for the future, even though from space the HST/WFC3 instrument represents the only option we have to investigate details of the spatial distribution of the Lyman continuum emission arising from z~2-4 galaxies.

preprint2012arXiv

Relativistic virialization in the Spherical Collapse model for Einstein-de Sitter and ΛCDM cosmologies

Spherical collapse has turned out to be a successful semi-analytic model to study structure formation in different DE models and theories of gravity. Nevertheless, the process of virialization is commonly studied on the basis of the virial theorem of classical mechanics. In the present paper, a fully generally-relativistic virial theorem based on the Tolman-Oppenheimer-Volkoff (TOV) solution for homogeneous, perfect-fluid spheres is constructed for the Einstein-de Sitter and ΛCDM cosmologies. We investigate the accuracy of classical virialization studies on cosmological scales and consider virialization from a more fundamental point of view. Throughout, we remain within general relativity and the class of FLRW models. The virialization equation is set up and solved numerically for the virial radius, y_{vir}, from which the virial overdensity Δ_{V} is directly obtained. Leading order corrections in the post-Newtonian framework are derived and quantified. In addition, problems in the application of this formalism to dynamical DE models are pointed out and discussed explicitly. We show that, in the weak field limit, the relative contribution of the leading order terms of the post-Newtonian expansion are of the order of 10^{-3}% and the solution of Wang & Steinhardt (1998) is precisely reproduced. Apart from the small corrections, the method could provide insight into the process of virialization from a more fundamental point of view.

preprint2012arXiv

The influence of mergers on the cluster temperature function and cosmological parameters derived from it

We develop a parameter-free analytic model to include the effects of mergers into the theoretical modelling of the X-ray temperature function of galaxy clusters. We include this description into our model for the cluster population based on fluctuations of the gravitational potential, which avoids any reference to mass. Comparisons with a numerical simulation reveal that the theoretical model is in good agreement with the simulation results. We show that building the model on the dynamics of spherical rather than ellipsoidal collapse yields better results if emission-weighted temperatures are used, while ellipsoidal collapse yields good agreement between model and simulation for mass-weighted temperatures. Analysing two different samples of X-ray clusters, we quantify the influence of mergers and a conversion between different temperature definitions on the joint determination of Omega_m0 and sigma_8. If effects of mergers are included, temperature functions based on cluster masses and on the gravitational potential built on spherical collapse are in good agreement with other cosmological probes without any conversion of temperatures.

preprint2012arXiv

The strongest gravitational lenses: I. The statistical impact of cluster mergers

For more than a decade now, it has been controversial whether or not the high rate of giant gravitational arcs and the largest observed Einstein radii are consistent with the standard cosmological model. Recent studies indicate that mergers provide an efficient mechanism to substantially increase the strong-lensing efficiency of individual clusters. Based on purely semi-analytic methods, we investigated the statistical impact of cluster mergers on the distribution of the largest Einstein radii and the optical depth for giant gravitational arcs of selected cluster samples. Analysing representative all-sky realizations of clusters at redshifts z < 1 and assuming a constant source redshift of z_s = 2.0, we find that mergers increase the number of Einstein radii above 10 arcsec (20 arcsec) by ~ 35 % (~ 55 %). Exploiting the tight correlation between Einstein radii and lensing cross sections, we infer that the optical depth for giant gravitational arcs with a length-to-width ratio > 7.5 of those clusters with Einstein radii above 10 arcsec (20 arcsec) increases by ~ 45 % (85 %). Our findings suggest that cluster mergers significantly influence in particular the statistical lensing properties of the strongest gravitational lenses. We conclude that semi-analytic studies must inevitably take these events into account before questioning the standard cosmological model on the basis of the largest observed Einstein radii and the statistics of giant gravitational arcs.

preprint2012arXiv

The strongest gravitational lenses: II. Is the large Einstein radius of MACS J0717.5+3745 in conflict with LCDM?

Can the standard cosmological model be questioned on the basis of a single observed extreme galaxy cluster? Usually, the word extreme refers directly to cluster mass, which is not a direct observable and thus subject to substantial uncertainty. Hence, it is desirable to extend studies of extreme clusters to direct observables, such as the Einstein radius (ER). We aim to evaluate the occurrence probability of the large observed ER of MACS J0717.5 within the standard LCDM cosmology. In particular, we want to model the distribution function of the single largest ER in a given cosmological volume and to study which underlying assumptions and effects have the strongest impact on the results. We obtain this distribution by a Monte Carlo approach, based on the semi-analytic modelling of the halo population on the past lightcone. After sampling the distribution, we fit the results with the general extreme value (GEV) distribution which we use for the subsequent analysis. We find that the distribution of the maximum ER is particularly sensitive to the precise choice of the halo mass function, lens triaxiality, the inner slope of the halo density profile and the mass-concentration relation. Using the distributions so obtained,we study the occurrence probability of the large ER of MACS J0717.5, finding that this system is not in tension with LCDM. We also find that the GEV distribution can be used to fit very accurately the sampled distributions and that all of them can be described by a Frechet distribution. With a multitude of effects that strongly influence the distribution of the single largest ER, it is more than doubtful that the standard LCDM cosmology can be ruled out on the basis of a single observation. If, despite the large uncertainties in the underlying assumptions, one wanted to do so, a much larger ER (> 100 arcsec) than that of MACS J0717.5 would have to be observed.

preprint2012arXiv

The Three-Dimensional Shapes of Galaxy Clusters

While clusters of galaxies are considered one of the most important cosmological probes, the standard spherical modelling of the dark matter and the intracluster medium is only a rough approximation. Indeed, it is well established both theoretically and observationally that galaxy clusters are much better approximated as triaxial objects. However, investigating the asphericity of galaxy clusters is still in its infancy. We review here this topic which is currently gathering a growing interest from the cluster community. We begin by introducing the triaxial geometry. Then we discuss the topic of deprojection and demonstrate the need for combining different probes of the cluster's potential. We discuss the different works that have been addressing these issues. We present a general parametric framework intended to simultaneously fit complementary data sets (X-ray, Sunyaev Zel'dovich and lensing data). We discuss in details the case of Abell 1689 to show how different models/data sets lead to different haloe parameters. We present the results obtained from fitting a 3D NFW model to X-ray, SZ, and lensing data for 4 strong lensing clusters. We argue that a triaxial model generally allows to lower the inferred value of the concentration parameter compared to a spherical analysis. This may alleviate tensions regarding, e.g. the over-concentration problem. However, we stress that predictions from numerical simulations rely on a spherical analysis of triaxial halos. Given that triaxial analysis will have a growing importance in the observational side, we advocate the need for simulations to be analysed in the very same way, allowing reliable and meaningful comparisons. Besides, methods intended to derive the three dimensional shape of galaxy clusters should be extensively tested on simulated multi-wavelength observations.

preprint2012arXiv

The Universal Einstein Radius Distribution from 10,000 SDSS Clusters

We present results from strong-lens modelling of 10,000 SDSS clusters, to establish the universal distribution of Einstein radii. Detailed lensing analyses have shown that the inner mass distribution of clusters can be accurately modelled by assuming light traces mass, successfully uncovering large numbers of multiple-images. Approximate critical curves and the effective Einstein radius of each cluster can therefore be readily calculated, from the distribution of member galaxies and scaled by their luminosities. We use a subsample of 10 well-studied clusters covered by both SDSS and HST to calibrate and test this method, and show that an accurate determination of the Einstein radius and mass can be achieved by this approach "blindly", in an automated way, and without requiring multiple images as input. We present the results of the first 10,000 clusters analysed in the range $0.1<z<0.55$, and compare them to theoretical expectations. We find that for this all-sky representative sample the Einstein radius distribution is log-normal in shape, with $< Log(θ_{e}\arcsec)>=0.73^{+0.02}_{-0.03}$, $σ=0.316^{+0.004}_{-0.002}$, and with higher abundance of large $θ_{e}$ clusters than predicted by $Λ$CDM. We visually inspect each of the clusters with $θ_{e}>40 \arcsec$ ($z_{s}=2$) and find that $\sim20%$ are boosted by various projection effects detailed here, remaining with $\sim40$ real giant-lens candidates, with a maximum of $θ_{e}=69\pm12 \arcsec$ ($z_{s}=2$) for the most massive candidate, in agreement with semi-analytic calculations. The results of this work should be verified further when an extended calibration sample is available.

preprint2012arXiv

Weak gravitational lensing with DEIMOS

We introduce a novel method for weak-lensing measurements, which is based on a mathematically exact deconvolution of the moments of the apparent brightness distribution of galaxies from the telescope's PSF. No assumptions on the shape of the galaxy or the PSF are made. The (de)convolution equations are exact for unweighted moments only, while in practice a compact weight function needs to be applied to the noisy images to ensure that the moment measurement yields significant results. We employ a Gaussian weight function, whose centroid and ellipticity are iteratively adjusted to match the corresponding quantities of the source. The change of the moments caused by the application of the weight function can then be corrected by considering higher-order weighted moments of the same source. Because of the form of the deconvolution equations, even an incomplete weighting correction leads to an excellent shear estimation if galaxies and PSF are measured with a weight function of identical size. We demonstrate the accuracy and capabilities of this new method in the context of weak gravitational lensing measurements with a set of specialized tests and show its competitive performance on the GREAT08 challenge data. A complete C++ implementation of the method can be requested from the authors.

preprint2011arXiv

Calibration biases in measurements of weak lensing

As recently shown by Viola et al., the common (KSB) method for measuring weak gravitational shear creates a non-linear relation between the measured and the true shear of objects. We investigate here what effect such a non-linear calibration relation may have on cosmological parameter estimates from weak lensing if a simpler, linear calibration relation is assumed. We show that the non-linear relation introduces a bias in the shear-correlation amplitude and thus a bias in the cosmological parameters Omega_matter and sigma_8. Its direction and magnitude depends on whether the point-spread function is narrow or wide compared to the galaxy images from which the shear is measured. Substantial over- or underestimates of the cosmological parameters are equally possible, depending also on the variant of the KSB method. Our results show that for trustable cosmological-parameter estimates from measurements of weak lensing, one must verify that the method employed is free from ellipticity-dependent biases or monitor that the calibration relation inferred from simulations is applicable to the survey at hand.

preprint2011arXiv

Lensed arc statistics: comparison of Millennium-simulation galaxy clusters to Hubble Space Telescope observations of an X-ray selected sample

It has been debated for a decade whether there is a large overabundance of strongly lensed arcs in galaxy clusters, compared to expectations from LambdaCDM cosmology. We perform ray tracing through the most massive halos of the Millennium simulation at several redshifts in their evolution, using the Hubble Ultra Deep Field as a source image, to produce realistic simulated lensed images. We compare the lensed arc statistics measured from the simulations to those of a sample of 45 X-ray selected clusters, observed with the Hubble Space Telescope, that we have analysed in Horesh et al. (2010). The observations and the simulations are matched in cluster masses, redshifts, observational effects, and the algorithmic arc detection and selection. At z=0.6 there are too few massive-enough clusters in the Millennium volume for a proper statistical comparison with the observations. At redshifts 0.3<z<0.5, however, we have large numbers of simulated and observed clusters, and the latter are an unbiased selection from a complete sample. For these redshifts, we find excellent agreement between the observed and simulated arc statistics, in terms of the mean number of arcs per cluster, the distribution of number of arcs per cluster, and the angular separation distribution. At z ~ 0.2 some conflict remains, with real clusters being ~3 times more efficient arc producers than their simulated counterparts. This may arise due to selection biases in the observed subsample at this redshift, to some mismatch in masses between the observed and simulated clusters, or to physical effects that arise at low redshift and enhance the lensing efficiency, but which are not represented by the simulations.

preprint2011arXiv

Shear-flexion cross-talk in weak-lensing measurements

Gravitational flexion, caused by derivatives of the gravitational tidal field, is potentially important for the analysis of the dark-matter distribution in gravitational lenses, such as galaxy clusters or the dark-matter haloes of galaxies. Flexion estimates rely on measurements of galaxy-shape distortions with spin-1 and spin-3 symmetry. We show in this paper that and how such distortions are generally caused not only by the flexion itself, but also by coupling terms of the form (shear $\times$ flexion), which have hitherto been neglected. Similar coupling terms occur between intrinsic galaxy ellipticities and the flexion. We show, by means of numerical tests, that neglecting these terms can introduce biases of up to 85% on the $F$ flexion and 150% on the $G$ flexion for galaxies with an intrinsic ellipticity dispersion of $σ_ε=0.3$. In general, this bias depends on the strength of the lensing fields, the ellipticity dispersion, and the concentration of the lensed galaxies. We derive a new set of equations relating the measured spin-1 and spin-3 distortions to the lensing fields up to first order in the shear, the flexion, the product of shear and flexion, and the morphological properties of the galaxy sample. We show that this new description is accurate with a bias $\leq 7%$ (spin-1 distortion) and $\leq 3%$ (spin-3 distortion) even close to points where the flexion approach breaks down due to merging of multiple images. We propose an explanation why a spin-3 signal could not be measured yet and comment on the difficulties in using a model-fitting approach to measure the flexion signal.

preprint2010arXiv

Gravitational Lensing

Gravitational lensing has developed into one of the most powerful tools for the analysis of the dark universe. This review summarises the theory of gravitational lensing, its main current applications and representative results achieved so far. It has two parts. In the first, starting from the equation of geodesic deviation, the equations of thin and extended gravitational lensing are derived. In the second, gravitational lensing by stars and planets, galaxies, galaxy clusters and large-scale structures is discussed and summarised.

preprint2010arXiv

Halo model description of the non-linear dark matter power spectrum at $k \gg 1$ Mpc$^{-1}$

Accurate knowledge of the non-linear dark-matter power spectrum is important for understanding the large-scale structure of the Universe, the statistics of dark-matter haloes and their evolution, and cosmological gravitational lensing. We analytically model the dark-matter power spectrum and its cross-power spectrum with dark-matter haloes. Our model extends the halo-model formalism, including realistic substructure population within individual dark-matter haloes and the scatter of the concentration parameter at fixed halo mass. We consider three prescriptions for the mass-concentration relation and two for the substructure distribution in dark-matter haloes. We show that this extension of the halo model mainly increases the predicted power on the small scales, and is crucial for proper modeling the cosmological weak-lensing signal due to low-mass haloes. Our extended formalism shows how the halo model approach can be improved in accuracy as one increases the number of ingredients that are calibrated from n-body simulations.

preprint2010arXiv

Soft clustering analysis of galaxy morphologies: A worked example with SDSS

Context: The huge and still rapidly growing amount of galaxies in modern sky surveys raises the need of an automated and objective classification method. Unsupervised learning algorithms are of particular interest, since they discover classes automatically. Aims: We briefly discuss the pitfalls of oversimplified classification methods and outline an alternative approach called "clustering analysis". Methods: We categorise different classification methods according to their capabilities. Based on this categorisation, we present a probabilistic classification algorithm that automatically detects the optimal classes preferred by the data. We explore the reliability of this algorithm in systematic tests. Using a small sample of bright galaxies from the SDSS, we demonstrate the performance of this algorithm in practice. We are able to disentangle the problems of classification and parametrisation of galaxy morphologies in this case. Results: We give physical arguments that a probabilistic classification scheme is necessary. The algorithm we present produces reasonable morphological classes and object-to-class assignments without any prior assumptions. Conclusions: There are sophisticated automated classification algorithms that meet all necessary requirements, but a lot of work is still needed on the interpretation of the results.

preprint2010arXiv

The lensing efficiencies of MACS X-ray selected versus RCS optically selected galaxy clusters

The statistics of strongly lensed arcs in samples of galaxy clusters provide information on cluster structure that is complementary to that from individual clusters. However, samples of clusters that have been analyzed to date have been either small, heterogeneous, or observed with limited angular resolution. We measure the lensed-arc statistics of 97 clusters imaged at high angular resolution with the Hubble Space Telescope, identifying lensed arcs using two automated arc detection algorithms. The sample includes similar numbers of X-ray selected (MACS) and optically selected (RCS) clusters, and spans cluster redshifts in the range 0.2 < z < 1. We compile a catalogue of 42 arcs in the X-ray selected subsample and 7 arcs in the optical subsample. All but five of these arcs are reported here for the first time. At 0.3 < z < 0.7, the X-ray selected clusters have a significantly higher mean frequency of arcs, 1.2+/-0.2 per cluster, versus 0.2+/-0.1 in the optical sample. The strikingly different lensing efficiencies indicate that X-ray clusters trace much larger mass concentrations, despite the similar optical luminosities of the X-ray and optical clusters. The mass difference is supported also by the lower space density of the X-ray clusters, and by the small Einstein radii of the few arcs in the optical sample. Higher-order effects, such as differences in concentration or substructure, may also contribute.

preprint2010arXiv

Triaxial collapse and virialisation of dark-matter haloes

We reconsider the ellipsoidal-collapse model and extend it in two ways: We modify the treatment of the external gravitational shear field, introducing a hybrid model in between linear and non-linear evolution, and we introduce a virialisation criterion derived from the tensor virial theorem to replace the ad-hoc criterion employed so far. We compute the collapse parameters delta_c and Delta_v and find that they increase with ellipticity e and decrease with prolaticity p. We marginalise them over the appropriate distribution of e and p and show the marginalised results as functions of halo mass and virialisation redshift. While the hybrid model for the external shear gives results very similar to those obtained from the non-linear model, ellipsoidal collapse changes the collapse parameters typically by (20...50)%, in a way increasing with decreasing halo mass and decreasing virialisation redshift. We qualitatively confirm the dependence on mass and virialisation redshift of a fitting formula for delta_c, but find noticeable quantitative differences in particular at low mass and high redshift. The derived mass function is in good agreement with mass functions recently proposed in the literature.

preprint2009arXiv

Lensed CMB temperature and polarization maps from the Millennium Simulation

We have constructed the first all-sky CMB temperature and polarization lensed maps based on a high-resolution cosmological N-body simulation, the Millennium Simulation (MS). We have exploited the lensing potential map obtained using a map-making procedure (Carbone et al. 2008) which integrates along the line-of-sight the MS dark matter distribution by stacking and randomizing the simulation boxes up to $z = 127$, and which semi-analytically supplies the large-scale power in the angular lensing potential that is not correctly sampled by the N-body simulation. The lensed sky has been obtained by properly modifying the latest version of the LensPix code (Lewis 2005) to account for the MS structures. We have also produced all-sky lensed maps of the so-called $ψ_E$ and $ψ_B$ potentials, which are directly related to the electric and magnetic types of polarization. The angular power spectra of the simulated lensed temperature and polarization maps agree well with semi-analytic estimates up to $l < 2500$, while on smaller scales we find a slight excess of power which we interpret as being due to non-linear clustering in the MS. We also observe how non-linear lensing power in the polarised CMB is transferred to large angular scales by suitably misaligned modes in the CMB and the lensing potential. This work is relevant in view of the future CMB probes, as a way to analyse the lensed sky and disentangle the contribution from primordial gravitational waves.

preprint2009arXiv

Statistics of gravitational potential perturbations: A novel approach to deriving the X-ray temperature function

Context. While the halo mass function is theoretically a very sensitive measure of cosmological models, masses of dark-matter halos are poorly defined, global, and unobservable quantities. Aims. We argue that local, observable quantities such as the X-ray temperatures of galaxy clusters can be directly compared to theoretical predictions without invoking masses. We derive the X-ray temperature function directly from the statistics of Gaussian random fluctuations in the gravitational potential. Methods. We derive the abundance of potential minima constrained by the requirement that they belong to linearly collapsed structures. We then use the spherical-collapse model to relate linear to non-linear perturbations, and the virial theorem to convert potential depths to temperatures. No reference is made to mass or other global quantities in the derivation. Results. Applying a proper high-pass filter that removes large enough modes from the gravitational potential, we derive an X-ray temperature function that agrees very well with the classical Press-Schechter approach on relevant temperature scales, but avoids the necessity of measuring masses. Conclusions. TThis first study shows that and how an X-ray temperature function of galaxy clusters can be analytically derived, avoiding the introduction of poorly defined global quantities such as halo masses. This approach will be useful for reducing scatter in observed cluster distributions and thus in cosmological conclusions drawn from them.

preprint2009arXiv

The Dark Universe

For a few years now, cosmology has a standard model. By this term, we mean a consistent theoretical background which is at the same time simple and broad enough to offer coherent explanations for the vast majority of cosmological phenomena. This review will briefly summarise the cosmological model, then proceed to discuss what we know from observations about the evolution of the Universe and its contents, and what we conclude about the origin and the future of the Universe and the structures it contains.

preprint2008arXiv

Deconvolution with Shapelets

We seek to find a shapelet-based scheme for deconvolving galaxy images from the PSF which leads to unbiased shear measurements. Based on the analytic formulation of convolution in shapelet space, we construct a procedure to recover the unconvolved shapelet coefficients under the assumption that the PSF is perfectly known. Using specific simulations, we test this approach and compare it to other published approaches. We show that convolution in shapelet space leads to a shapelet model of order $n_{max}^h = n_{max}^g + n_{max}^f$ with $n_{max}^f$ and $n_{max}^g$ being the maximum orders of the intrinsic galaxy and the PSF models, respectively. Deconvolution is hence a transformation which maps a certain number of convolved coefficients onto a generally smaller number of deconvolved coefficients. By inferring the latter number from data, we construct the maximum-likelihood solution for this transformation and obtain unbiased shear estimates with a remarkable amount of noise reduction compared to established approaches. This finding is particularly valid for complicated PSF models and low $S/N$ images, which renders our approach suitable for typical weak-lensing conditions.

preprint2006arXiv

Detecting Sunyaev-Zel'dovich clusters with PLANCK: III. Properties of the expected SZ-cluster sample

The PLANCK-mission is the most sensitive all-sky submillimetric mission currently being planned and prepared. Special emphasis is given to the observation of clusters of galaxies by their thermal Sunyaev-Zel'dovich (SZ) effect. In this work, the results of a simulation are presented that combines all-sky maps of the thermal and kinetic SZ-effect with cosmic microwave background (CMB) fluctuations, Galactic foregrounds (synchrotron emission, thermal emission from dust, free-free emission and rotational transitions of carbon monoxide molecules) and sub-millimetric emission from planets and asteroids of the Solar System. Observational issues, such as PLANCKs beam shapes, frequency response and spatially non-uniform instrumental noise have been incorporated. Matched and scale-adaptive multi-frequency filtering schemes have been extended to spherical coordinates and are now applied to the data sets in order to isolate and amplify the weak thermal SZ-signal. The properties of the resulting SZ-cluster sample are characterised in detail: Apart from the number of clusters as a function of cluster parameters such as redshift z and total mass M, the distribution n(sigma)d sigma of the detection significance sigma, the number of detectable clusters in relation to the model cluster parameters entering the filter construction, the position accuracy of an SZ-detection and the cluster number density as a function of ecliptic latitude beta is examined.

preprint2005arXiv

Weak lensing in the second post-Newtonian approximation: Gravitomagnetic potentials and the integrated Sachs-Wolfe effect

Dark matter currents in the large-scale structure give rise to gravitomagnetic terms in the metric, which affect the light propagation. Corrections to the weak lensing power spectrum due to these gravitomagnetic potentials are evaluated by perturbation theory. A connection between gravitomagnetic lensing and the integrated Sachs-Wolfe (iSW) effect is drawn, which can be described by a line-of-sight integration over the divergence of the gravitomagnetic vector potential. This allows the power spectrum of the iSW-effect to be derived within the framework of the same formalism as derived for gravitomagnetic lensing and reduces the iSW-effect to a second order lensing phenomenon. The three-dimensional power spectra are projected by means of a generalised Limber-equation to yield the angular power spectra. While gravitomagnetic corrections to the weak lensing spectrum are negligible at observationally accessible scales, the angular power spectrum of the iSW-effect should be detectable as a correction to the CMB spectrum up to multipoles of l \~ 100 with the PLANCK-satellite.

preprint2004arXiv

Constraints on Dark Energy Models from Galaxy Clusters with Multiple Arcs

We make an exploratory study of how well dark energy models can be constrained using lensed arcs at different redshifts behind cluster lenses. Arcs trace the critical curves of clusters, and the growth of critical curves with source redshift is sensitive to the dark-energy equation of state. Using analytical models and numerically simulated clusters, we explore the key factors involved in using cluster arcs as a probe of dark energy. We quantify the sensitivity to lens mass, concentration and ellipticity with analytical models that include the effects of dark energy on halo structure. We show with simple examples how degeneracies between mass models and cosmography may be broken using arcs at multiple redshifts or additional constraints on the lens density profile. However we conclude that the requirements on the data are so stringent that it is very unlikely that robust constraints can be obtained from individual clusters. We argue that surveys of clusters, analyzed in conjunction with numerical simulations, are a more promising prospect for arc-cosmography. We use such numerically simulated clusters to estimate how large a sample of clusters/arcs could provide interesting constraints on dark energy models. We focus on the scatter produced by differences in the mass distribution of individual clusters. We find from our sample of simulated clusters that at least 1000 pairs of arcs are needed to obtain constraints if the mass distribution of individual clusters is taken to be undetermined. We discuss several unsolved problems that need study to fully develop this method for precision studies with future surveys.

preprint2004arXiv

Detecting Sunyaev-Zel'dovich clusters with PLANCK: I. Construction of all-sky thermal and kinetic SZ-maps

All-sky thermal and kinetic Sunyaev-Zel'dovich (SZ) maps are presented for assessing how well the PLANCK-mission can find and characterise clusters of galaxies, especially in the presence of primary anisotropies of the cosmic microwave background (CMB) and various galactic and ecliptic foregrounds. The maps have been constructed from numerical simulations of structure formation in a standard LCDM cosmology and contain all clusters out to redshifts of z = 1.46 with masses exceeding 5e13 M_solar/h. By construction, the maps properly account for the evolution of cosmic structure, the halo-halo correlation function, the evolving mass function, halo substructure and adiabatic gas physics. The velocities in the kinetic map correspond to the actual density environment at the cluster positions. We characterise the SZ-cluster sample by measuring the distribution of angular sizes, the integrated thermal and kinetic Comptonisations, the source counts in the three relevant PLANCK-channels, and give the angular power-spectra of the SZ-sky. While our results are broadly consistent with simple estimates based on scaling relations and spherically symmetric cluster models, some significant differences are seen which may affect the number of cluster detectable by PLANCK.

preprint2004arXiv

Detecting Sunyaev-Zel'dovich clusters with PLANCK: II. Foreground components and optimised filtering schemes

The PLANCK mission is the most sensitive all-sky CMB experiment currently planned. The High Frequency Instrument (HFI) will be especially suited for observing clusters of galaxies by their thermal Sunyaev-Zel'dovich (SZ) effect. In order to assess PLANCKs SZ-capabilities in the presence of spurious signals, a simulation is presented that combines maps of the thermal and kinetic SZ-effects with a realisation of the cosmic microwave background (CMB), in addition to Galactic foregrounds (synchrotron emission, free-free emission, thermal emission from dust, CO-line radiation) as well as the sub-millimetric emission from celestial bodies of our Solar system. Additionally, observational issues such as the finite angular resolution and spatially non-uniform instrumental noise of PLANCKs sky maps are taken into account, yielding a set of all-sky flux maps, the auto-correlation and cross-correlation properties of which are examined in detail. In the second part of the paper, filtering schemes based on scale-adaptive and matched filtering are extended to spherical data sets, that enable the amplification of the weak SZ-signal in the presence of all contaminations stated above. The theory of scale-adaptive and matched filtering in the framework of spherical maps is developed, the resulting filter kernel shapes are discussed and their functionality is verified.

preprint2003arXiv

Numerical study of halo concentrations in dark-energy cosmologies

We study the concentration parameters, their mass dependence and redshift evolution, of dark-matter halos in different dark-energy cosmologies with constant and time-variable equation of state, and compare them with "standard'' Lambda-CDM and OCDM models. We find that previously proposed algorithms for predicting halo concentrations can be well adapted to dark-energy models. When centred on the analytically expected values, halo concentrations show a log-normal distribution with a uniform standard deviation of ~0.2. The dependence of averaged halo concentrations on mass and redshift permits a simple fit of the form (1+z) c=c0 (M/M0)^a, with a~-0.1 throughout. We find that the cluster concentration depends on the dark energy equation of state at the cluster formation redshift z_{coll} through the linear growth factor D_+(z_{coll}). As a simple correction accounting for dark-energy cosmologies, we propose scaling c0 from Lambda-CDM with the ratio of linear growth factors, c0 -> c0 D_+(z_{coll})/D_{+,Lambda-CDM}(z_{coll}).

preprint2002arXiv

Improving the Accuracy of Cosmic Magnification Statistics

The systematic magnification of background sources by the weak gravitational-lensing effects of foreground matter, also called cosmic magnification, is becoming an efficient tool both for measuring cosmological parameters and for exploring the distribution of galaxies relative to the dark matter. We extend here the formalism of magnification statistics by estimating the contribution of second-order terms in the Taylor expansion of the magnification and show that the effect of these terms was previously underestimated. We test our analytical predictions against numerical simulations and demonstrate that including second-order terms allows the accuracy of magnification-related statistics to be substantially improved. We also show, however, that both numerical and analytical estimates can provide only lower bounds to real correlation functions, even in the weak lensing regime. We propose to use count-in-cells estimators rather than correlation functions for measuring cosmic magnification since they can more easily be related to correlations measured in numerical simulations.

preprint2002arXiv

Probing Intracluster Magnetic Fields with Cosmic Microwave Background Polarization

Intracluster magnetic fields with $ \sim μG $ strength induce Faraday rotation on the cosmic microwave background (CMB) polarization. Measurements of this effect can potentially probe the detailed structure of intracluster magnetic fields across clusters, since the CMB polarization is a continuously varying field on the sky, in contrast to the conventional method restricted by the limited number of radio sources behind or inside a cluster. We here construct a method for extracting information on magnetic fields from measurements of the effect, combined with possible observations of the Sunyaev-Zel'dovich effect and X-ray emission for the same cluster which are needed to reconstruct the electron density fields. Employing the high-resolution magneto-hydrodynamic simulations performed by Dolag, Bartelmann & Lesch (1999) as a realistic model of magnetized intracluster gas distribution, we demonstrate how our reconstruction technique can well reproduce the magnetic fields, i.e., the spherically average radial profiles of the field strength and the coherence length.

preprint1997arXiv

QSO-galaxy correlations due to weak lensing in arbitrary Friedmann-Lemaitre cosmologies

We calculate the angular cross-correlation function between background QSOs and foreground galaxies induced by the weak lensing effect of large-scale structures. Results are given for arbitrary Friedmann-Lemaitre cosmologies. The non-linear growth of density perturbations is included. Compared to the linear growth, the non-linear growth increases the correlation amplitude by about an order of magnitude in an Einstein-de Sitter universe, and by even more for lower Omega_0. The dependence of the correlation amplitude on the cosmological parameters strongly depends on the normalization of the power spectrum. The QSO-galaxy cross-correlation function is most sensitive to density structures on scales in the range (1-10) Mpc/h, where the normalization of the power spectrum to the observed cluster abundance appears most appropriate. In that case, the correlation strength changes by less than a factor of <~ 2 when Omega_0 varies between 0.3 and 1, quite independent of the value of Omega_Lambda. For Omega_0 <~ 0.3, the correlation strength increases with decreasing Omega_0, and it scales approximately linearly with the Hubble constant h.

preprint1996arXiv

A Comparison of X-ray and Strong Lensing Properties of Simulated X-ray Clusters

We use gas-dynamical simulations of galaxy clusters to compare their X-ray and strong lensing properties. Special emphasis is laid on mass estimates. The cluster masses range between 6 x 10^14 solar masses and 4 x 10^15 solar masses, and they are examined at redshifts between 1 and 0. We compute the X-ray emission of the intracluster gas by thermal bremsstrahlung, add background contamination, and mimic imaging and spectral observations with current X-ray telescopes. Although the beta model routinely provides excellent fits to the X-ray emission profiles, the derived masses are typically biased low because of the restricted range of radii within which the fit can be done. For beta values of ~ 2/3, which is the average in our numerically simulated sample, the mass is typically underestimated by ~ 40 per cent. The masses of clusters which exhibit pronounced substructure are often substantially underestimated. We suggest that the ratio between peak temperature and emission-weighted average cluster temperature may provide a good indicator for ongoing merging and, therefore, for unreliable mass estimates. X-ray mass estimates are substantially improved if we fit a King density profile rather than the beta model to the X-ray emission, thereby dropping the degree of freedom associated with beta. Clusters selected for their strong lensing properties are typically dynamically more active than typical clusters. Bulk flows in the intracluster gas contain a larger than average fraction of the internal energy of the gas in such objects, hence the measured gas temperatures are biased low. The bulk of the optical depth for arc formation is contributed by clusters with intermediate rather than high X-ray luminosity. Arcs occur predominantly in clusters which exhibit substructure and are not in an equilibrium state. Finally we explain why the

preprint1996arXiv

Aperture Multipole Moments from Weak Gravitational Lensing

The projected mass of a gravitational lens inside (circular) apertures can be derived from the measured shear inside an annulus which is caused by the tidal field of the deflecting mass distribution. Here we show that also the multipoles of the two-dimensional mass distribution can be derived from the shear in annuli. We derive several expressions for these mass multipole moments in terms of the shear, which allow large flexibility in the choice of a radial weight function. In contrast to determining multipole moments from weak-lensing mass reconstructions, this approach allows to quantify the signal-to-noise ratio of the multipole moments directly from the observed galaxy ellipticities, and thus to estimate the significance of the multipole detection. Radial weight functions can therefore be chosen such as to optimize the significance of the detection given an assumed radial mass profile. Application of our formulae to numerically simulated clusters demonstrates that the quadrupole moment of realistic cluster models can be detected with high signal-to-noise ratio S/N; in about 85 per cent of the simulated cluster fields S/N >~ 3. We also show that the shear inside a circular annulus determines multipole moments inside and outside the annulus. This is relevant for clusters whose central region is too bright to allow the observation of the shear of background galaxies, or which extend beyond the CCD. We also generalize the aperture mass equation to the case of `radial' weight functions which are constant on arbitrarily-shaped curves which are not necessarily self-similar.

preprint1996arXiv

Arcs from a Universal Dark-Matter Halo Profile

Navarro, Frenk, \& White have recently found numerically that the density profile of dark-matter halos can be described by a universal two-parameter function over a broad range of halo masses. The profile is singular, approaching the halo center with $ρ\propto r^{-1}$. It had been argued previously that radially distorted, gravitationally lensed images of background sources in galaxy clusters, so-called radial arcs, required a flat core in the cluster density profile. Such radial arcs have so far been detected in two galaxy clusters, in apparent contradiction with a singular density profile. I show here that the profile suggested by Navarro et al. can produce radial arcs despite its central singularity, and describe how the two parameters of the profile can be determined in clusters where radial and tangential arcs are observed. I then apply this analysis to the two clusters where radial arcs were detected. In both cases, the redshifts of the radial arcs are yet unknown, hence definitive conclusions on the profile parameters cannot yet be drawn. Numerically determined values for the parameters of cluster-sized halos can, however, be used to predict the range of the unknown arc redshifts, thus providing a direct observational test for the proposed density profile. A potential difficulty with the profile is that the radial magnification of tangential arcs is large, hence tangential arcs should be thick or their sources should be very thin in the radial direction.

preprint1994arXiv

Arc Statistics with Realistic Cluster Potentials. III. a Systematic Effect on Cluster Mass Estimates

It is shown in this paper that deviations of galaxy cluster lenses from spherical symmetry can render mass estimates for galaxy clusters based on the formation of large arcs systematically to high. Numerical models show that the mass needed for producing large arcs in clusters can be notably smaller than expected from simple spherically symmetric lens models. The reason is that the enhanced tidal effect in asymmetric and substructured lenses can compensate for part of the convergence necessary for strong lensing effects. An analytic argument is given to explain why deviations from radial symmetry will in general decrease the required lens mass. Mass estimates assuming radially symmetric lenses are on average too high by a factor of $\simeq1.6$, and with a probability of $\simeq20\%$ by a factor of $\simeq2$.

preprint1994arXiv

Cluster mass estimates from weak lensing

The reliability of cluster lens reconstruction techniques based on weak lensing is studied in terms of the accuracy of their reproduction of the total cluster mass as a function of distance from the cluster center. To do so, a variety of reconstruction algorithms is applied to synthetic lensing data created using a sample of 60 numerically modeled clusters, and the mass reconstruction is compared to the known deflector mass. The results can be summarized as follows: (1) Reconstruction algorithms which require integrations extending over the entire real plane yield unreliable results, because they give rise to boundary effects which are hard to control; mass overestimates are more likely and more substantial in this case than underestimates. (2) Reconstruction techniques which avoid these boundary effects yield reliable lower bounds to the cluster mass. The tightness of such bounds depends on the size of the field, which can be extended synthetically to improve the results considerably. For the sample of numerical cluster models, the best lower bounds, achieved by combining synthetic field extension with non-linear, finite-field reconstruction, decline from 100% to 80% of the true cluster mass going from the cluster center to an angular distance of 2.5'. The 80% error bars of the lower mass bounds are +/- 10% to 15%.

preprint1994arXiv

On the spin parameter of dark-matter haloes

The study by White (1984) on the growth of angular momentum in dark haloes is extended towards a more detailed investigation of the spin parameter $λ\equiv L\sqrt{E}/{G M^{2.5}}$. Starting from the Zel'dovich approximation to structure formation, a dark halo is approximated by a homogeneous ellipsoid with the inertial tensor of the (highly irregular) Lagrangian region $Υ$ from which the dark halo forms. Within this approximation, an expression for the spin parameter can be derived, which depends on the geometry of $Υ$, the cosmological density parameter $Ω_0$, the overdensity of the dark halo, and the tidal torque exerted on it. For Gaussian random fields, this expression can be evaluated statistically. As a result, we derive a probability distribution of the spin parameter which gives $λ\simeq0.07^{+0.04}_{-0.05}$, consistent with numerical investigations. This probability distribution steeply rises with increasing spin parameter, reaching its maximum at $λ\simeq0.025$. The 10 (50,90) percentile values are $λ=0.02$ (0.05,0.11, respectively). There is a weak anticorrelation of the spin parameter with the peak height $ν$ of the density fluctuation field $λ\propto ν^{-0.29}$. The dependence on $Ω_0$ and the variance $σ$ of the density-contrast field is very weak; there is only a marginal tendency for the spin parameter to be slightly larger for late-forming objects in an open universe. Due to the weak dependence on $σ$, our results should be quite generally applicable and independent on