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.

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 $\delta_{\textrm{crit}}$. We consider a spherical region as a halo, if its density today exceeds the virial overdensity threshold $\Delta$. 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 $\delta$ as an uncorrelated random walk of the smoothing scale $R$. The second assumes $\delta(R)$ to be correlated. We find that the resulting HMFs correspond well to the HMF found in numerical simulations if the correlation of $\delta(R)$ is taken into account. Furthermore, the HMF depends only weakly on the choice of the density threshold $\Delta$.

Since the discovery of the accelerated cosmic expansion, one of the most important tasks in observational cosmology is to determine the nature of the dark energy. We should build our understanding on a minimum of assumptions in order to avoid biases from assumed cosmological models. The two most important functions describing the evolution of the universe and its structures are the expansion function E(a) and the linear growth factor D_+(a). The expansion function has been determined in previous papers in a model-independent way using distance moduli to type-Ia supernovae and assuming only a metric theory of gravity, spatial isotropy and homogeneity. Here, we extend this analysis in three ways: (1) We extend the data sample by combining the Pantheon measurements of type-Ia supernovae with measurements of baryonic acoustic oscillations; (2) we substantially simplify and generalise our method for reconstructing the expansion function; and (3) we use the reconstructed expansion function to determine the linear growth factor of cosmic structures, equally independent of specific assumptions on an underlying cosmological model other than the usual spatial symmetries. We show that the result is quite insensitive to the initial conditions for solving the growth equation, leaving the present-day matter-density parameter {\Omega}_m0 as the only relevant parameter for an otherwise purely empirical and accurate determination of the growth factor.

The formation of cosmic structures is an important diagnostic for both the dynamics of the cosmological model and the underlying theory of gravity. At the linear level of these structures, certain degeneracies remain between different cosmological models and alternative gravity theories. It is thus indispensable to study the non-linear, late-time evolution of cosmic structures to try and disentangle their fundamental properties caused by the cosmological model or gravity theory itself. Conventionally, non-linear cosmic structure formation is studied by means of computationally expensive numerical simulations. Since these inevitably suffer from shot noise and are too time consuming to systematically scrutinize large parameter spaces of cosmological models or fundamental theories, analytical methods are needed to overcome the limitations of numerical simulations. Recently, a new analytic approach to non-linear cosmic structure formation has been proposed based on kinetic field theory for classical particle ensembles. Within this theory, a closed, analytic, non-perturbative and parameter-free equation could be derived for the non-linear power spectrum of cosmic density perturbations which agrees very well with numerically simulated results to wave numbers $k\lesssim10\,h\,\mathrm{Mpc}^{-1}$ at redshift $z = 0$. In this Letter, we study for the first time the implications of alternative gravity theories for non-linear cosmic structure formation applying this promising new analytic framework. As an illustrative example, we consider vector-tensor theories, which support very interesting isotropic cosmological solutions.

We investigate the implications of string Swampland criteria for alternative theories of gravity. Even though this has not rigorously been proven, there is some evidence that exact de-Sitter solutions with a positive cosmological constant cannot be successfully embedded into string theory, and that the low energy effective field theories containing a scalar field $\pi$ with a potential $V$ in the habitable landscape should satisfy the Swampland criteria $|V'|/V\ge c\sim\mathcal{O}(1)$. As a paradigmatic class of modified gravity theories for inflation and dark energy, we consider the extensively studied family of Horndeski Lagrangians in view of cosmological observations. Apart from a possible potential term, they contain derivative self-interactions as the Galileon and non-minimal couplings to the gravity sector. Hence, our conclusions on the Galileon sector can be also applied to many other alternative theories with scalar helicity modes containing derivative interactions such as massive gravity and generalized Proca. In the presence of such derivative terms, the dynamics of the scalar mode is substantially modified, and imposing the cosmological evolution constrained by observations places tight constraints on $c$ within the Swampland conjecture.

Galaxy clusters provide a multitude of observational data across wavelengths and their structure and morphology are of considerable interest in cosmology as well as astrophysics. We develop a framework that allows the combination of lensing and non-lensing observations in a free-form and mesh-free approach to infer the projected mass distribution of individual galaxy clusters. This method can be used to test common assumptions on the morphology of clusters in parametric models. We make use of the lensing reconstruction code SaWLens2 and expand its capabilities by incorporating an estimate of the projected gravitational potential based on X-ray data that are deprojected using the local Richardson-Lucy method and used to infer the Newtonian potential of the cluster and we discuss how potentially arising numerical artefacts can be treated. We demonstrate the feasibility of our method on a simplified mock NFW halo and on a cluster from a realistic hydrodynamical simulation and show how the combination of X-ray and weak lensing data can affect a free-form reconstruction, improving the accuracy in the central region in some cases by a factor of two.

Recently, Bartelmann et al. presented a 'Kinetic Field Theory' (KFT) formalism to tackle the difficulties of large scale structure formation. In this approach, the dynamics of a non-equilibrium ensemble of classical particles are examined based on methods of statistical field theory. So far, only contributions coming from dark matter were considered, which is assumed to pose an accurate description of our universe on very large scales. Nevertheless, going to smaller scales, also baryonic contributions have to be taken into account. Building on the ideas of Viermann et al. we present an effective particle model of hydrodynamics to describe baryonic matter in a cosmological framework. Using this model, the baryonic density contrast power spectrum is computed to lowest perturbative order within the resummed KFT framework of Lilow et al. We discuss the qualitative differences to the dark matter case and perform a quantitative comparison to the baryonic spectrum obtained from Eulerian perturbation theory. A subsequent paper will resolve the problem of coupling both theories describing dark and baryonic matter, respectively, to gain a full model of cosmic matter. Even though our focus is on cosmological systems only, we want to emphasize that all methods presented here are of a quite general fashion, making it applicable also to other fields.

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.

In earlier work, we have developed a nonequilibrium statistical field theory description of cosmic structure formation, dubbed Kinetic Field Theory (KFT), which is based on the Hamiltonian phase-space dynamics of classical particles and thus remains valid beyond shell-crossing. Here, we present an exact reformulation of the KFT framework that allows to resum an infinite subset of terms appearing in the original perturbative expansion of KFT. We develop the general formalism of this resummed KFT, including a diagrammatic language for the resummed perturbation theory, and compute the lowest-order results for the power spectra of the dark matter density contrast and momentum density. This allows us to derive analytically how the linear growth of the largest structures emerges from Newtonian particle dynamics alone, which, to our knowledge, is the first time this has been achieved.

In this paper we study long range modifications of gravity in the consistent framework of bigravity, which introduces a second massive spin-2 field and allows to continuously interpolate between the regime of General Relativity (mediated by a massless spin-2 field) and massive gravity (mediated by a massive spin-2 field). In particular we derive for the first time the equations for light deflection in this framework and study the effect on the lensing potential of galaxy clusters. By comparison of kinematic and lensing mass reconstructions, stringent bounds can be set on the parameter space of the new spin-2 fields. Furthermore, we investigate galactic rotation curves and the effect on the observable dark matter abundance within this framework.

A recent comparison of the massive galaxy cluster Abell 2744 with the Millennium XXL (MXXL) N-body simulation has hinted at a tension between the observed substructure distribution and the predictions of LambdaCDM. Follow-up investigations indicated that this could be due to the contribution from the host halo and the subhalo finding algorithm used. To be independent of any subhalo finding algorithm, we therefore investigate the particle data of the MXXL simulation directly. We propose a new method to find substructures in 2D mass maps using a wavelet transform, which treats the simulation and observations equally. Using the same criteria to define a subhalo in observations and simulated data, we find three Abell 2744 analogues in the MXXL simulation. Thus the observations in Abell 2744 are in agreement with the predictions of LambdaCDM. We investigate the reasons for the discrepancy between the results obtained from the SUBFIND and full particle data analyses. We find that this is due to incompatible substructure definitions in observations and SUBFIND.

In this recent Letter [1], we studied the implications of string Swampland conjectures for quintessence models of dark energy. A couple of days ago, a paper appeared [2] which refers to our work but misrepresents what we have done. Here, we set the record straight.

In this Letter, we study the implications of string Swampland criteria for dark energy in view of ongoing and future cosmological observations. If string theory should be the ultimate quantum gravity theory, there is evidence that exact de Sitter solutions with a positive cosmological constant cannot describe the fate of the late-time universe. Even though cosmological models with dark energy given by a scalar field $\pi$ evolving in time are not in direct tension with string theory, they have to satisfy the Swampland criteria $|\Delta\pi|<d\sim\mathcal{O}(1)$ and $|V'|/V>c\sim\mathcal{O}(1)$, where $V$ is the scalar field potential. In view of the restrictive implications that the Swampland criteria have on dark energy, we investigate the accuracy needed for future observations to tightly constrain standard dark-energy models. We find that current 3-$\sigma$ constraints with $c \lesssim 1.35$ are still well in agreement with the string Swampland criteria. However, Stage-4 surveys such as Euclid, LSST and DESI, tightly constraining the equation of state $w(z)$, will start putting surviving quintessence models into tensions with the string Swampland criteria by demanding $c<0.4$. We further investigate whether any idealised futuristic survey will ever be able to give a decisive answer to the question whether the cosmological constant would be preferred over a time-evolving dark-energy model within the Swampland criteria. Hypothetical surveys with a reduction in the uncertainties by a factor of $\sim20$ compared to Euclid would be necessary to reveal strong tension between quintessence models obeying the string Swampland criteria and observations by pushing the allowed values down to $c<0.1$. In view of such perspectives, there will be fundamental observational limitations with future surveys.

We report the results of $EasyCritics$, a fully automated algorithm for the efficient search of strong-lensing (SL) regions in wide-field surveys, applied to the Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS). By using only the photometric information of the brightest elliptical galaxies distributed over a wide redshift range ($\smash{0.2 \lesssim z \lesssim 0.9}$) and without requiring the identification of arcs, our algorithm produces lensing potential models and catalogs of critical curves of the entire survey area. We explore several parameter set configurations in order to test the efficiency of our approach. In a specific configuration, $EasyCritics$ generates only $\sim1200$ possibly super-critical regions in the CFHTLS area, drastically reducing the effective area for inspection from $154$ sq. deg to $\sim0.623$ sq. deg, $i.e.$ by more than two orders of magnitude. Among the pre-selected SL regions, we identify 32 of the 44 previously known lenses on the group and cluster scale, and discover 9 new promising lens candidates. The detection rate can be easily improved to $\sim82\%$ by a simple modification in the parameter set, but at the expense of increasing the total number of possible SL candidates. Note that $EasyCritics$ is fully complementary to other arc-finders since we characterize lenses instead of directly identifying arcs. Although future comparisons against numerical simulations are required for fully assessing the efficiency of $EasyCritics$, the algorithm seems very promising for upcoming surveys covering $\smash{10^{4}}$ sq. deg, such as the $Euclid$ mission and $LSST$, where the pre-selection of candidates for any kind of SL analysis will be indispensable due to the expected enormous data volume.

Building upon the recently developed formalism of Kinetic Field Theory (KFT) for cosmic structure formation by Bartelmann et al., we investigate a kinematic relationship between diffusion and gravitational interactions in cosmic structure formation. In the first part of this work we explain how the process of structure formation in KFT can be separated into three processes, particle diffusion, the accumulation of structure due to initial momentum correlations and interactions relative to the inertial motion of particles. We study these processes by examining the time derivative of the non-linear density power spectrum in the Born approximation. We observe that diffusion and accumulation are delicately balanced because of the Gaussian form of the initial conditions, and that the net diffusion, resulting from adding these two counteracting contributions, approaches the contributions from the interactions in amplitude over time. This hints at a kinematic relation between diffusion and interactions in KFT. Indeed, in the second part, we show that the response of the system to arbitrary gradient forces is directly related to the evolution of particle diffusion in the form of kinematic fluctuation-dissipation relations (FDRs). This result is independent of the interaction potential. We show that this relationship roots in a time-reversal symmetry of the underlying generating functional. Furthermore, our studies demonstrate how FDRs originating from purely kinematic arguments can be used in theories far from equilibrium.

The mass of galaxy clusters is not a direct observable, nonetheless it is commonly used to probe cosmological models. Based on the combination of all main cluster observables, that is, the X-ray emission, the thermal Sunyaev-Zel'dovich (SZ) signal, the velocity dispersion of the cluster galaxies, and gravitational lensing, the gravitational potential of galaxy clusters can be jointly reconstructed. We derive the two main ingredients required for this joint reconstruction: the potentials individually reconstructed from the observables and their covariance matrices, which act as a weight in the joint reconstruction. We show here the method to derive these quantities. The result of the joint reconstruction applied to a real cluster will be discussed in a forthcoming paper. We apply the Richardson-Lucy deprojection algorithm to data on a two-dimensional (2D) grid. We first test the 2D deprojection algorithm on a $\beta$-profile. Assuming hydrostatic equilibrium, we further reconstruct the gravitational potential of a simulated galaxy cluster based on synthetic SZ and X-ray data. We then reconstruct the projected gravitational potential of the massive and dynamically active cluster Abell 2142, based on the X-ray observations collected with XMM-Newton and the SZ observations from the Planck satellite. Finally, we compute the covariance matrix of the projected reconstructed potential of the cluster Abell 2142 based on the X-ray measurements collected with XMM-Newton. The gravitational potentials of the simulated cluster recovered from synthetic X-ray and SZ data are consistent, even though the potential reconstructed from X-rays shows larger deviations from the true potential. Regarding Abell 2142, the projected gravitational cluster potentials recovered from SZ and X-ray data reproduce well the projected potential inferred from gravitational-lensing observations. (abridged)

In this work we review the theory of the spherical collapse model and critically analyse the aspects of the numerical implementation of its fundamental equations. By extending a recent work by Herrera et al. 2017, we show how different aspects, such as the initial integration time, the definition of constant infinity and the criterion for the extrapolation method (how close the inverse of the overdensity has to be to zero at the collapse time) can lead to an erroneous estimation (a few per mill error which translates to a few percent in the mass function) of the key quantity in the spherical collapse model: the linear critical overdensity $\delta_{\rm c}$, which plays a crucial role for the mass function of halos. We provide a better recipe to adopt in designing a code suitable to a generic smooth dark energy model and we compare our numerical results with analytic predictions for the EdS and the $\Lambda$CDM models. We further discuss the evolution of $\delta_{\rm c}$ for selected classes of dark energy models as a general test of the robustness of our implementation. We finally outline which modifications need to be taken into account to extend the code to more general classes of models, such as clustering dark energy models and non-minimally coupled models.

We derive a non-perturbative, closed, analytic equation for the non-linear power spectrum of cosmic density fluctuations. This result is based on our kinetic field theory (KFT) of cosmic structure formation and on evaluating particle interactions in the Born approximation. At redshift zero, relative deviations between our analytic result and typical numerical results are ~ 15 % on average up to wave numbers of k <= 10 h/Mpc. The theory underlying our analytic equation is fully specified once the statistical properties of the initial density and momentum fluctuations are set. It has no further adjustable parameters. Apart from this equation, our main result is that the characteristic non-linear deformations of the power spectrum at late cosmic times are determined by the initial momentum correlations and a partial compensation between diffusion damping and particle interactions.

In recent work we developed a description of cosmic large scale structure formation in terms of non-equilibrium ensembles of classical particles, with time evolution obtained in the framework of a statistical field theory. In these works, the initial Gaussian correlations between particles have so far been treated perturbatively or restricted to pure momentum correlations. Here we treat the correlations between all phase-space coordinates exactly by adopting a diagrammatic language for the different forms of correlations, directly inspired by the Mayer cluster expansion. We will demonstrate that explicit expressions for phase-space density cumulants of arbitrary $n$-point order, which fully capture the non-linear coupling of free streaming kinematics due to initial correlations, can be obtained from a simple set of Feynman rules. These cumulants will be the foundation for further investigations of interacting perturbation theory.

We present a new strong lensing analysis of the galaxy cluster MACS J1206.2-0847 (MACS 1206), at z=0.44, using deep spectroscopy from CLASH-VLT and VLT/MUSE archival data in combination with imaging from the Cluster Lensing and Supernova survey with Hubble. MUSE observations enable the spectroscopic identification of 23 new multiply imaged sources, extending the previous compilations by a factor of approximately five. In total, we use the positional measurements of 82 spectroscopic multiple images belonging to 27 families at z=1.0-6.1 to reconstruct the projected total mass distribution of MACS 1206. Remarkably, 11 multiple images are found within 50 kpc of the brightest cluster galaxy, making this an unprecedented set of constraints for the innermost projected mass distribution of a galaxy cluster. We thus find that, although dynamically relaxed, the smooth matter component (dark matter plus hot gas) of MACS 1206 shows a significant asymmetry, which closely follows the asymmetric distribution of the stellar component (galaxy members and intracluster light). We determine the value of the innermost logarithmic slope of the projected total mass density profile and find it to be close to the canonical Navarro-Frenk-White value. We demonstrate that this quantity is very robust against different parametrizations of the diffuse mass component; however, this is not the case when only one central image is used in the mass reconstruction. We also show that the mass density profile from our new strong lensing model is in very good agreement with dynamical and X-ray measurements at larger radii, where they overlap.

We study the dynamical behaviour of gauge-invariant linear perturbations in spherically symmetric dust cosmologies including a cosmological constant. In contrast to spatially homogeneous FLRW models, the reduced degree of spatial symmetry causes a non-trivial dynamical coupling of gauge-invariant quantities already at first order perturbation theory and the strength and influence of this coupling on the spacetime evolution is investigated here. We present results on the underlying dynamical equations augmented by a cosmological constant and integrate them numerically. We also present a method to derive cosmologically relevant initial variables for this setup. Estimates of angular power spectra for each metric variable are computed and evaluated on the central observer's past null cone. By comparing the full evolution to the freely evolved initial profiles, the coupling strength will be determined for a best fit radially inhomogeneous patch obtained in previous works (see Redlich et. al. (2014)). We find that coupling effects are not noticeable within the cosmic variance limit and can therefore safely be neglected for a relevant cosmological scenario. On the contrary, we find very strong coupling effects in a best fit spherical void model matching the distance redshift relation of SNe which is in accordance with previous findings using parametric void models.

In earlier work, we have developed a Kinetic Field Theory (KFT) for cosmological structure formation and showed that the non-linear density-fluctuation power spectrum known from numerical simulations can be reproduced quite well even if particle interactions are taken into account to first order only. Besides approximating gravitational interactions, we had to truncate the initial correlation hierarchy of particle momenta at the second order. Here, we substantially simplify KFT. We show that its central object, the free generating functional, can be factorized, taking the full hierarchy of momentum correlations into account. The factors appearing in the generating functional, which we identify as non-linearly evolved density-fluctuation power spectra, have a universal form and can thus be tabulated for fast access in perturbation schemes. In this paper, we focus on a complete evaluation of the free generating functional of KFT, not including particle interactions yet. This implies that the non-linearly evolved power spectra contain a damping term which reflects that structures are being wiped out at late times by free streaming. Once particle interactions will be taken into account, they will compensate this damping. If we suppress this damping in a way suggested by the fluctuation-dissipation relations of KFT, our results show that the complete hierarchy of initial momentum correlations is responsible for a large part of the characteristic non-linear deformation and the mode transport in the density-fluctuation power spectrum. Without any adjustable parameters, KFT accurately reproduces the scale at which non-linear evolution sets in. Finally, we further develop perturbation theory based on the factorization of the generating functional and propose a diagrammatic scheme for the perturbation terms.

We present accurate photometric redshifts for galaxies observed by the Cluster Lensing and Supernova survey with Hubble (CLASH). CLASH observed 25 massive galaxy cluster cores with the Hubble Space Telescope in 16 filters spanning 0.2 - 1.7 $\mu$m. Photometry in such crowded fields is challenging. Compared to our previously released catalogs, we make several improvements to the photometry, including smaller apertures, ICL subtraction, PSF matching, and empirically measured uncertainties. We further improve the Bayesian Photometric Redshift (BPZ) estimates by adding a redder elliptical template and by inflating the photometric uncertainties of the brightest galaxies. The resulting photometric redshift accuracies are dz/(1+z) $\sim$ 0.8\%, 1.0\%, and 2.0\% for galaxies with I-band F814W AB magnitudes $<$ 18, 20, and 23, respectively. These results are consistent with our expectations. They improve on our previously reported accuracies by a factor of 4 at the bright end and a factor of 2 at the faint end. Our new catalog includes 1257 spectroscopic redshifts, including 382 confirmed cluster members. We also provide stellar mass estimates. Finally, we include lensing magnification estimates of background galaxies based on our public lens models. Our new catalog of all 25 CLASH clusters is available via MAST. The analysis techniques developed here will be useful in other surveys of crowded fields, including the Frontier Fields and surveys carried out with J-PAS and JWST.

The massive substructures found in Abell 2744 by Jauzac et al. (2016) present a challenge to the cold dark matter paradigm due to their number and proximity to the cluster centre. We use one of the biggest N-body simulations, the Millennium XXL, to investigate the substructure in a large sample of massive dark matter haloes. A range of effects which influence the comparison with the observations is considered, extending the preliminary evaluation carried out by Jauzac et al. (2016). There are many tens of haloes in the simulation with a total mass comparable with or larger than that of Abell 2744. However, we find no haloes with a number and distribution of massive substructures (> 5 10^13 Msun) that is close to that inferred from the observations of Abell 2744. The application of extreme value statistics suggests that we would need a simulation of at least ten times the volume of the Millennium XXL to find a single dark matter halo with a similar internal structure to Abell 2744. Explaining the distribution of massive substructures in clusters is a new hurdle for hierarchical models to negotiate, which is not weakened by appeals to baryonic physics or uncertainty over the nature of the dark matter particle.

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.

Building upon the recent pioneering work by Mazenko and by Das and Mazenko, we develop a microscopic, non-equilibrium, statistical field theory for initially correlated canonical ensembles of classical microscopic particles obeying Hamiltonian dynamics. Our primary target is cosmic structure formation, where initial Gaussian correlations in phase space are believed to be set by inflation. We give an exact expression for the generating functional of this theory and work out suitable approximations. We specify the initial correlations by a power spectrum and derive general expressions for the correlators of the density and the response field. We derive simple closed expressions for the lowest-order contributions to the non-linear cosmological power spectrum, valid for arbitrary wave numbers. We further calculate the bispectrum expected in this theory within these approximations and the power spectrum of cosmic density fluctuations to first order in the gravitational interaction, using a recent improvement of the Zel'dovich approximation. We show that, with a modification motivated by the adhesion approximation, 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 even within first-order perturbation theory. Our results present the first fully analytic calculation of the non-linear power spectrum of cosmic structures.

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.

The statistics of peaks in weak gravitational lensing maps is a promising technique to constrain cosmological parameters in present and future surveys. Here we investigate its power when using general extreme value statistics which is very sensitive to the exponential tail of the halo mass function. To this end, we use an analytic method to quantify the number of weak lensing peaks caused by galaxy clusters, large-scale structures and observational noise. Doing so, we further improve the method in the regime of high signal-to-noise ratios dominated by non-linear structures by accounting for the embedding of those counts into the surrounding shear caused by large scale structures. We derive the extreme value and order statistics for both over-densities (positive peaks) and under-densities (negative peaks) and provide an optimized criterion to split a wide field survey into sub-fields in order to sample the distribution of extreme values such that the expected objects causing the largest signals are mostly due to galaxy clusters. We find good agreement of our model predictions with a ray-tracing $N$-body simulation. For a Euclid-like survey, we find tight constraints on $\sigma_8$ and $\Omega_\text{m}$ with relative uncertainties of $\sim 10^{-3}$. In contrast, the equation of state parameter $w_0$ can be constrained only with a $10\%$ level, and $w_\text{a}$ is out of reach even if we include redshift information.

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.

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, $\beta\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 $\beta\approx0. 6$ within $\approx0.5\,\mathrm{Mpc}$ and stays constant further out.

The abundance of galaxy clusters is in principle a powerful tool to constrain cosmological parameters, especially $\Omega_\mathrm{m}$ and $\sigma_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 $\Omega_\mathrm{m}$ and $\sigma_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$\sigma$-level.

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.

We present profiles of temperature (Tx), gas mass, and hydrostatic mass estimated from new and archival X-ray observations of CLASH clusters. We compare measurements derived from XMM and Chandra observations with one another and compare both to gravitational lensing mass profiles derived with CLASH HST and ground-based lensing data. Radial profiles of Chandra and XMM electron density and enclosed gas mass are nearly identical, indicating that differences in hydrostatic masses inferred from X-ray observations arise from differences in Tx measurements. Encouragingly, cluster Txs are consistent with one another at ~100-200 kpc radii but XMM Tx systematically decline relative to Chandra Tx at larger radii. The angular dependence of the discrepancy suggests additional investigation on systematics such as the XMM point spread function correction, vignetting and off-axis responses. We present the CLASH-X mass-profile comparisons in the form of cosmology-independent and redshift-independent circular-velocity profiles. Ratios of Chandra HSE mass profiles to CLASH lensing profiles show no obvious radial dependence in the 0.3-0.8 Mpc range. However, the mean mass biases inferred from the WL and SaWLens data are different. e.g., the weighted-mean value at 0.5 Mpc is <b> = 0.12 for the WL comparison and <b> = -0.11 for the SaWLens comparison. The ratios of XMM HSE mass profiles to CLASH lensing profiles show a pronounced radial dependence in the 0.3-1.0 Mpc range, with a weighted mean mass bias of value rising to <b>~0.3 at ~1 Mpc for the WL comparison and <b> of 0.25 for SaWLens comparison. The enclosed gas mass profiles from both Chandra and XMM rise to a value 1/8 times the total-mass profiles inferred from lensing at 0.5 Mpc and remain constant outside of that radius, suggesting that [8xMgas] profiles may be an excellent proxy for total-mass profiles at >0.5 Mpc in massive galaxy clusters.

We study the evolution of linear perturbations in a Lema\^itre-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.

We update the all-sky Planck catalogue of 1227 clusters and cluster candidates (PSZ1) published in March 2013, derived from Sunyaev-Zeldovich (SZ) effect detections using the first 15.5 months of Planck satellite observations. Addendum. We deliver an updated version of the PSZ1 catalogue, reporting the further confirmation of 86 Planck-discovered clusters. In total, the PSZ1 now contains 947 confirmed clusters, of which 214 were confirmed as newly discovered clusters through follow-up observations undertaken by the Planck Collaboration. The updated PSZ1 contains redshifts for 913 systems, of which 736 (~80.6%) are spectroscopic, and associated mass estimates derived from the Y_z mass proxy. We also provide a new SZ quality flag, derived from a novel artificial neural network classification of the SZ signal, for the remaining 280 candidates. Based on this assessment, the purity of the updated PSZ1 catalogue is estimated to be 94%. In this release, we provide the full updated catalogue and an additional readme file with further information on the Planck SZ detections.

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, $\kappa$, and $\sim20\%$ in the magnification, $\mu$. 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 (\Sigma)/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.

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.

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.

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 $\rho$ 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.

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.

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.

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.

We utilise a two-color Lyman-Break selection criterion to search for z~9-10 galaxies over the first 19 clusters in the CLASH program. A systematic search yields three z~9-10 candidates. While we have already reported the most robust of these candidates, MACS1149-JD, two additional z~9 candidates are also found and have H_{160}-band magnitudes of ~26.2-26.9. A careful assessment of various sources of contamination suggests <~1 contaminants for our z~9-10 selection. To determine the implications of these search results for the LF and SFR density at z~9, we introduce a new differential approach to deriving these quantities in lensing fields. Our procedure is to derive the evolution by comparing the number of z~9-10 galaxy candidates found in CLASH with the number of galaxies in a slightly lower redshift sample (after correcting for the differences in selection volumes), here taken to be z~8. This procedure takes advantage of the fact that the relative volumes available for the z~8 and z~9-10 selections behind lensing clusters are not greatly dependent on the details of the lensing models. We find that the normalization of the UV LF at z~9 is just 0.28_{-0.20}^{+0.39}\times that at z~8, ~1.4_{-0.8}^{+3.0}x lower than extrapolating z~4-8 LF results. While consistent with the evolution in the UV LF seen at z~4-8, these results marginally favor a more rapid evolution at z>8. Compared to similar evolutionary findings from the HUDF, our result is less insensitive to large-scale structure uncertainties, given our many independent sightlines on the high-redshift universe.

We utilize 16 band Hubble Space Telescope (HST) observations of 18 lensing clusters obtained as part of the Cluster Lensing And Supernova survey with Hubble (CLASH) Multi-Cycle Treasury program to search for $z\sim6-8$ galaxies. We report the discovery of 204, 45, and 13 Lyman-break galaxy candidates at $z\sim6$, $z\sim7$, and $z\sim8$, respectively, identified from purely photometric redshift selections. This large sample, representing nearly an order of magnitude increase in the number of magnified star-forming galaxies at $z\sim 6-8$ presented to date, is unique in that we have observations in four WFC3/UVIS UV, seven ACS/WFC optical, and all five WFC3/IR broadband filters, which enable very accurate photometric redshift selections. We construct detailed lensing models for 17 of the 18 clusters to estimate object magnifications and to identify two new multiply lensed $z \gtrsim 6$ candidates. The median magnifications over the 17 clusters are 4, 4, and 5 for the $z\sim6$, $z\sim7$, and $z\sim8$ samples, respectively, over an average area of 4.5 arcmin$^2$ per cluster. We compare our observed number counts with expectations based on convolving "blank" field UV luminosity functions through our cluster lens models and find rough agreement down to $\sim27$ mag, where we begin to suffer significant incompleteness. In all three redshift bins, we find a higher number density at brighter observed magnitudes than the field predictions, empirically demonstrating for the first time the enhanced efficiency of lensing clusters over field surveys. Our number counts also are in general agreement with the lensed expectations from the cluster models, especially at $z\sim6$, where we have the best statistics.

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 $\Lambda$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\sigma$ 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.

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.

The ESA's Planck satellite, dedicated to studying the early Universe and its subsequent evolution, was launched 14 May 2009 and has been scanning the microwave and submillimetre sky continuously since 12 August 2009. This paper gives an overview of the mission and its performance, the processing, analysis, and characteristics of the data, the scientific results, and the science data products and papers in the release. The science products include maps of the CMB and diffuse extragalactic foregrounds, a catalogue of compact Galactic and extragalactic sources, and a list of sources detected through the SZ effect. The likelihood code used to assess cosmological models against the Planck data and a lensing likelihood are described. Scientific results include robust support for the standard six-parameter LCDM model of cosmology and improved measurements of its parameters, including a highly significant deviation from scale invariance of the primordial power spectrum. The Planck values for these parameters and others derived from them are significantly different from those previously determined. Several large-scale anomalies in the temperature distribution of the CMB, first detected by WMAP, are confirmed with higher confidence. Planck sets new limits on the number and mass of neutrinos, and has measured gravitational lensing of CMB anisotropies at greater than 25 sigma. Planck finds no evidence for non-Gaussianity in the CMB. Planck's results agree well with results from the measurements of baryon acoustic oscillations. Planck finds a lower Hubble constant than found in some more local measures. Some tension is also present between the amplitude of matter fluctuations derived from CMB data and that derived from SZ data. The Planck and WMAP power spectra are offset from each other by an average level of about 2% around the first acoustic peak.

We present the results of a numerical study based on the analysis of the MUSIC-2 simulations, aimed at estimating the expected concentration-mass relation for the CLASH cluster sample. We study nearly 1400 halos simulated at high spatial and mass resolution, which were projected along many lines-of-sight each. We study the shape of both their density and surface-density profiles and fit them with a variety of radial functions, including the Navarro-Frenk-White, the generalised Navarro-Frenk-White, and the Einasto density profiles. We derive concentrations and masses from these fits and investigate their distributions as a function of redshift and halo relaxation. We use the X-ray image simulator X-MAS to produce simulated Chandra observations of the halos and we use them to identify objects resembling the X-ray morphologies and masses of the clusters in the CLASH X-ray selected sample. We also derive a concentration-mass relation for strong-lensing clusters. We find that the sample of simulated halos which resemble the X-ray morphology of the CLASH clusters is composed mainly by relaxed halos, but it also contains a significant fraction of un-relaxed systems. For such a sample we measure an average 2D concentration which is ~11% higher than found for the full sample of simulated halos. After accounting for projection and selection effects, the average NFW concentrations of CLASH clusters are expected to be intermediate between those predicted in 3D for relaxed and super-relaxed halos. Matching the simulations to the individual CLASH clusters on the basis of the X-ray morphology, we expect that the NFW concentrations recovered from the lensing analysis of the CLASH clusters are in the range [3-6], with an average value of 3.87 and a standard deviation of 0.61. Simulated halos with X-ray morphologies similar to those of the CLASH clusters are affected by a modest orientation bias.

We present a new determination of the concentration-mass relation for galaxy clusters based on our comprehensive lensing analysis of 19 X-ray selected galaxy clusters from the Cluster Lensing and Supernova Survey with Hubble (CLASH). Our sample spans a redshift range between 0.19 and 0.89. We combine weak lensing constraints from the Hubble Space Telescope (HST) and from ground-based wide field data with strong lensing constraints from HST. The result are reconstructions of the surface-mass density for all CLASH clusters on multi-scale grids. Our derivation of NFW parameters yields virial masses between 0.53 x 10^15 and 1.76 x 10^15 M_sol/h and the halo concentrations are distributed around c_200c ~ 3.7 with a 1-sigma significant negative trend with cluster mass. We find an excellent 4% agreement between our measured concentrations and the expectation from numerical simulations after accounting for the CLASH selection function based on X-ray morphology. The simulations are analyzed in 2D to account for possible biases in the lensing reconstructions due to projection effects. The theoretical concentration-mass (c-M) relation from our X-ray selected set of simulated clusters and the c-M relation derived directly from the CLASH data agree at the 90% confidence level.

We describe the all-sky Planck catalogue of clusters and cluster candidates derived from Sunyaev--Zeldovich (SZ) effect detections using the first 15.5 months of Planck satellite observations. The catalogue contains 1227 entries, making it over six times the size of the Planck Early SZ (ESZ) sample and the largest SZ-selected catalogue to date. It contains 861 confirmed clusters, of which 178 have been confirmed as clusters, mostly through follow-up observations, and a further 683 are previously-known clusters. The remaining 366 have the status of cluster candidates, and we divide them into three classes according to the quality of evidence that they are likely to be true clusters. The Planck SZ catalogue is the deepest all-sky cluster catalogue, with redshifts up to about one, and spans the broadest cluster mass range from (0.1 to 1.6) 10^{15}Msun. Confirmation of cluster candidates through comparison with existing surveys or cluster catalogues is extensively described, as is the statistical characterization of the catalogue in terms of completeness and statistical reliability. The outputs of the validation process are provided as additional information. This gives, in particular, an ensemble of 813 cluster redshifts, and for all these Planck clusters we also include a mass estimated from a newly-proposed SZ-mass proxy. A refined measure of the SZ Compton parameter for the clusters with X-ray counter-parts is provided, as is an X-ray flux for all the Planck clusters not previously detected in X-ray surveys.

We aim at constraining the assembly history of clusters by studying the intra cluster light (ICL) properties, estimating its contribution to the fraction of baryons in stars, f*, and understanding possible systematics/bias using different ICL detection techniques. We developed an automated method, GALtoICL, based on the software GALAPAGOS to obtain a refined version of typical BCG+ICL maps. We applied this method to our test case MACS J1206.2-0847, a massive cluster located at z=0.44, that is part of the CLASH sample. Using deep multi-band SUBARU images, we extracted the surface brightness (SB) profile of the BCG+ICL and we studied the ICL morphology, color, and contribution to f* out to R500. We repeated the same analysis using a different definition of the ICL, SBlimit method, i.e. a SB cut-off level, to compare the results. The most peculiar feature of the ICL in MACS1206 is its asymmetric radial distribution, with an excess in the SE direction and extending towards the 2nd brightest cluster galaxy which is a Post Starburst galaxy. This suggests an interaction between the BCG and this galaxy that dates back to t <= 1.5 Gyr. The BCG+ICL stellar content is 8% of M_(*,500) and the (de-) projected baryon fraction in stars is f*=0.0177 (0.0116), in excellent agreement with recent results. The SBlimit method provides systematically higher ICL fractions and this effect is larger at lower SB limits. This is due to the light from the outer envelopes of member galaxies that contaminate the ICL. Though more time consuming, the GALtoICL method provides safer ICL detections that are almost free of this contamination. This is one of the few ICL study at redshift z > 0.3. At completion, the CLASH/VLT program will allow us to extend this analysis to a statistically significant cluster sample spanning a wide redshift range: 0.2<z<0.6.

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)

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.

We present a strong lensing system in which a double source is imaged 5 times by 2 early-type galaxies. We take advantage in this target of the multi-band photometry obtained as part of the CLASH program, complemented by the spectroscopic data of the VLT/VIMOS and FORS2 follow-up campaign. We use a photometric redshift of 3.7 for the source and confirm spectroscopically the membership of the 2 lenses to the galaxy cluster MACS J1206.2-0847 at redshift 0.44. We exploit the excellent angular resolution of the HST/ACS images to model the 2 lenses in terms of singular isothermal sphere profiles and derive robust effective velocity dispersions of (97 +/- 3) and (240 +/- 6) km/s. The total mass distribution of the cluster is also well characterized by using only the local information contained in this lensing system, that is located at a projected distance of more than 300 kpc from the cluster luminosity center. According to our best-fitting lensing and composite stellar population models, the source is magnified by a total factor of 50 and has a luminous mass of about (1.0 +/- 0.5) x 10^{9} M_{Sun}. By combining the total and luminous mass estimates of the 2 lenses, we measure luminous over total mass fractions projected within the effective radii of 0.51 +/- 0.21 and 0.80 +/- 0.32. With these lenses we can extend the analysis of the mass properties of lens early-type galaxies by factors that are about 2 and 3 times smaller than previously done with regard to, respectively, velocity dispersion and luminous mass. The comparison of the total and luminous quantities of our lenses with those of astrophysical objects with different physical scales reveals the potential of studies of this kind for investigating the internal structure of galaxies. These studies, made possible thanks to the CLASH survey, will allow us to go beyond the current limits posed by the available lens samples in the field.

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\rho)=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 $\rho$ is the density. We plan to further improve our constraint by applying the same procedure to all clusters from the ongoing CLASH-VLT program.

The Cluster Lensing And Supernovae survey with Hubble (CLASH) is an Hubble Space Telescope (HST) Multi-Cycle Treasury program observing 25 massive galaxy clusters. CLASH observations are carried out in 16 bands from UV to NIR to derive accurate and reliable estimates of photometric redshifts. We present the CLASH photometric redshifts and study the photometric redshift accuracy of the arcs in more detail for the case of MACS1206.2-0847. We use the publicly available Le Phare and BPZ photometric redshift codes on 17 CLASH galaxy clusters. Using Le Phare code for objects with StoN>=10, we reach a precision of 3%(1+z) for the strong lensing arcs, which is reduced to 2.4%(1+z) after removing outliers. For galaxies in the cluster field the corresponding values are 4%(1+z) and 3%(1+z). Using mock galaxy catalogues, we show that 3%(1+z) precision is what one would expect from the CLASH photometry when taking into account extinction from dust, emission lines and the finite range of SEDs included in the photo-z template library. We study photo-z results for different aperture photometry and find that the SExtractor isophotal photometry works best. Le Phare and BPZ give similar photo-z results for the strong lensing arcs as well as galaxies of the cluster field. Results are improved when optimizing the photometric aperture shape showing an optimal aperture size around 1" radius giving results which are equivalent to isophotal photometry. Tailored photometry of the arcs improve the photo-z results.

We present a quintuply lensed z ~ 6 candidate discovered in the field of the galaxy cluster RXC J2248.7-4431 (z ~ 0.348) targeted within the Cluster Lensing and Supernova survey with Hubble (CLASH) and selected in the deep HST Frontier Fields survey. Thanks to the CLASH 16-band HST imaging, we identify the quintuply lensed z ~ 6 candidate as an optical dropout in the inner region of the cluster, the brightest image having magAB=24.81+-0.02 in the f105w filter. We perform a detailed photometric analysis to verify its high-z and lensed nature. We get as photometric redshift z_phot ~ 5.9, and given the extended nature and NIR colours of the lensed images, we rule out low-z early type and galactic star contaminants. We perform a strong lensing analysis of the cluster, using 13 families of multiple lensed images identified in the HST images. Our final best model predicts the high-z quintuply lensed system with a position accuracy of 0.8''. The magnifications of the five images are between 2.2 and 8.3, which leads to a delensed UV luminosity of L_1600 ~ 0.5L*_1600 at z=6. We also estimate the UV slope from the observed NIR colours, finding a steep beta=-2.89+-0.38. We use singular and composite stellar population SEDs to fit the photometry of the hiz candidate, and we conclude that it is a young (age <300 Myr) galaxy with mass of M ~ 10^8Msol, subsolar metallicity (Z<0.2Zsol) and low dust content (AV ~ 0.2-0.4).

We present VIsible Multi-Object Spectrograph (VIMOS) observations of a z 6 galaxy quintuply imaged by the Frontier Fields galaxy cluster RXC J2248.7-4431 (z=0.348). This sub-L^*, high-z galaxy has been recently discovered by Monna et al. (2013) using dropout techniques with the 16-band HST photometry acquired as part of the Cluster Lensing And Supernova survey with Hubble (CLASH). Obtained as part of the CLASH-VLT survey, the VIMOS medium-resolution spectra of this source show a very faint continuum between ~8700A and ~9300A and a prominent emission line at 8643A, which can be readily identified with Lyman-alpha at z=6.110. The emission line exhibits an asymmetric profile, with a more pronounced red wing. The rest-frame equivalent width of the line is EW=79+-10A. After correcting for magnification, the star-formation rate (SFR) estimated from the Lya line is SFR(Lya)=11 M_{sol}/yr and that estimated from the UV data is SFR(UV)=3 M_{sol}/yr. We estimate that the effective radius of the source is R_e<~0.4 kpc, which implies a star formation surface mass density Sigma_{SFR}>6 M_{sol}/yr/kpc^2 and, using the Kennicutt-Schmidt relation, a gas surface mass density Sigma_{gas}>10^3 M_{sol}/pc^2. Our results support the idea that this magnified, distant galaxy is a young and compact object with 0.4 L^* at z=6, with comparable amount of mass in gas and stars. Future follow-up observations with ALMA will provide valuable insight into the SFR and molecular gas content of this source. In the spirit of the Frontier Fields initiative, we also publish the redshifts of several multiply imaged sources and other background objects which will help improving the strong lensing model of this galaxy cluster.

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.

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\sigma_\kappa. 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{\sigma} departure above the maximum simulated halo mass at this redshift.

We perform a model independent reconstruction of the cosmic expansion rate based on type Ia supernova data. Using the Union 2.1 data set, we show that the Hubble parameter behaviour allowed by the data without making any hypothesis about cosmological model or underlying gravity theory is consistent with a flat LCDM universe having H_0 = 70.43 +- 0.33 and Omega_m=0.297 +- 0.020, weakly dependent on the choice of initial scatter matrix. This is in closer agreement with the recently released Planck results (H_0 = 67.3 +- 1.2, Omega_m = 0.314 +- 0.020) than other standard analyses based on type Ia supernova data. We argue this might be an indication that, in order to tackle subtle deviations from the standard cosmological model present in type Ia supernova data, it is mandatory to go beyond parametrized approaches.

We provide a new observational test for a key prediction of the \Lambda 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 $\sigma$ agreement with \Lambda 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.

We use an unprecedented data-set of about 600 redshifts for cluster members, obtained as part of a VLT/VIMOS large programme, to constrain the mass profile of the z=0.44 cluster MACS J1206.2-0847 over the radial range 0-5 Mpc (0-2.5 virial radii) using the MAMPOSSt and Caustic methods. We then add external constraints from our previous gravitational lensing analysis. We invert the Jeans equation to obtain the velocity-anisotropy profiles of cluster members. With the mass-density and velocity-anisotropy profiles we then obtain the first determination of a cluster pseudo-phase-space density profile. The kinematics and lensing determinations of the cluster mass profile are in excellent agreement. This is very well fitted by a NFW model with mass M200=(1.4 +- 0.2) 10^15 Msun and concentration c200=6 +- 1, only slightly higher than theoretical expectations. Other mass profile models also provide acceptable fits to our data, of (slightly) lower (Burkert, Hernquist, and Softened Isothermal Sphere) or comparable (Einasto) quality than NFW. The velocity anisotropy profiles of the passive and star-forming cluster members are similar, close to isotropic near the center and increasingly radial outside. Passive cluster members follow extremely well the theoretical expectations for the pseudo-phase-space density profile and the relation between the slope of the mass-density profile and the velocity anisotropy. Star-forming cluster members show marginal deviations from theoretical expectations. This is the most accurate determination of a cluster mass profile out to a radius of 5 Mpc, and the only determination of the velocity-anisotropy and pseudo-phase-space density profiles of both passive and star-forming galaxies for an individual cluster [abridged]

Growing observational evidence now indicates that nebular line emission has a significant impact on the rest-frame optical fluxes of z~5-7 galaxies observed with Spitzer. This line emission makes z~5-7 galaxies appear more massive, with lower specific star formation rates. However, corrections for this line emission have been very difficult to perform reliably due to huge uncertainties on the overall strength of such emission at z>~5.5. Here, we present the most direct observational evidence yet for ubiquitous high-EW [OIII]+Hbeta line emission in Lyman-break galaxies at z~7, while also presenting a strategy for an improved measurement of the sSFR at z~7. We accomplish this through the selection of bright galaxies in the narrow redshift window z~6.6-7.0 where the IRAC 4.5 micron flux provides a clean measurement of the stellar continuum light. Observed 4.5 micron fluxes in this window contrast with the 3.6 micron fluxes which are contaminated by the prominent [OIII]+Hbeta lines. To ensure a high S/N for our IRAC flux measurements, we consider only the brightest (H_{160}<26 mag) magnified galaxies we have identified in CLASH and other programs targeting galaxy clusters. Remarkably, the mean rest-frame optical color for our bright seven-source sample is very blue, [3.6]-[4.5]=-0.9+/-0.3. Such blue colors cannot be explained by the stellar continuum light and require that the rest-frame EW of [OIII]+Hbeta be greater than 637 Angstroms for the average source. The bluest four sources from our seven-source sample require an even more extreme EW of 1582 Angstroms. Our derived lower limit for the mean [OIII]+Hbeta EW could underestimate the true EW by ~2x based on a simple modeling of the redshift distribution of our sources. We can also set a robust lower limit of >~4 Gyr^-1 on the specific star formation rates based on the mean SED for our seven-source sample. (abridged)

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 $\sigma=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.

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.

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.

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%.

The existence of an arc statistics problem was at the center of a strong debate in the last fifteen years. With the aim to clarify if the optical depth for giant gravitational arcs by galaxy clusters in the so called concordance model is compatible with observations, several studies were carried out which helped to significantly improve our knowledge of strong lensing clusters, unveiling their extremely complex internal structure. In particular, the abundance and the frequency of strong lensing events like gravitational arcs turned out to be a potentially very powerful tool to trace the structure formation. However, given the limited size of observational and theoretical data-sets, the power of arc statistics as a cosmological tool has been only minimally exploited so far. On the other hand, the last years were characterized by significant advancements in the field, and several cluster surveys that are ongoing or planned for the near future seem to have the potential to make arc statistics a competitive cosmological probe. Additionally, recent observations of anomalously large Einstein radii and concentrations in galaxy clusters have reinvigorated the debate on the arc statistics problem. In this paper, we review the work done so far on arc statistics, focussing on what is the lesson we learned and what is likely to improve in the next years.

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.

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.

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 $\delta_{\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.

{abridged} We study the imprints on the formation and evolution of cosmic structures of dynamical dark energy models, characterized by an oscillating equation of state. The redshift evolution of the equation of state parameter w(z) for dark energy is characterized by two parameters, describing the amplitude and the frequency of the oscillations. We consider six different oscillating dark energy models, each characterized by a different set of parameter values. Under the common assumption that dark energy is not clustering on the scales of interest, we study different aspects of cosmic structure formation. In particular, we self-consistently solve the spherical collapse problem. We then estimate the behavior of several cosmological observables, such as the linear growth factor, the Integrated Sachs-Wolfe (ISW) effect, the number counts of massive structures, and the matter and cosmic shear power spectra. We show that, independently of the amplitude and the frequency of the dark energy oscillations, none of the aforementioned observables show an oscillating behavior as a function of redshift. This is a consequence of the said observables' being integrals over some functions of the expansion rate over cosmic history. We also notice that deviations with respect to the expectations for a fiducial LambdaCDM cosmology are generically small, and in the majority of the cases distinguishing an oscillating dark energy model would be difficult. Exceptions to this conclusion are provided by the cosmic shear power spectrum, which for some of the models shows a difference at the level of \sim 10% over a wide range of angular scales, and the abundance of galaxy clusters, which is modified at the $\sim 10-20%$ level at $z \gtrsim 0.6$ for future wide weak lensing surveys.

While usually cosmological initial conditions are assumed to be Gaussian, inflationary theories can predict a certain amount of primordial non-Gaussianity which can have an impact on the statistical properties of the lensing observables. In order to evaluate this effect, we build a large set of realistic maps of different lensing quantities starting from light-cones extracted from large dark-matter only N-body simulations with initial conditions corresponding to different levels of primordial local non-Gaussianity strength $f_{\rm NL}$. Considering various statistical quantities (PDF, power spectrum, shear in aperture, skewness and bispectrum) we find that the effect produced by the presence of primordial non-Gaussianity is relatively small, being of the order of few per cent for values of $|f_{\rm NL}|$ compatible with the present CMB constraints and reaching at most 10-15 per cent for the most extreme cases with $|f_{\rm NL}|=1000$. We also discuss the degeneracy of this effect with the uncertainties due to the power spectrum normalization $\sigma_8$ and matter density parameter $\Omega_{\rm m}$, finding that an error in the determination of $\sigma_8$ ($\Omega_{\rm m}$) of about 3 (10) per cent gives differences comparable with non-Gaussian models having $f_{\rm NL}=\pm 1000$. These results suggest that the possible presence of an amount of primordial non-Gaussianity corresponding to $|f_{\rm NL}|=100$ is not hampering a robust determination of the main cosmological parameters in present and future weak lensing surveys, while a positive detection of deviations from the Gaussian hypothesis is possible only breaking the degeneracy with other cosmological parameters and using data from deep surveys covering a large fraction of the sky.

We perform a strong-lensing analysis of the merging galaxy cluster MACS J0416.1-2403 (M0416; z=0.42) in recent CLASH/HST observations. We identify 70 new multiple images and candidates of 23 background sources in the range 0.7<z_{phot}<6.14 including two probable high-redshift dropouts, revealing a highly elongated lens with axis ratio ~5:1, and a major axis of ~100\arcsec (z_{s}~2). Compared to other well-studied clusters, M0416 shows an enhanced lensing efficiency. Although the critical area is not particularly large (~0.6 \square\arcmin; z_{s}~2), the number of multiple images, per critical area, is anomalously high. We calculate that the observed elongation boosts the number of multiple images, \emph{per critical area}, by a factor of ~2.5\times, due to the increased ratio of the caustic area relative to the critical area. Additionally, we find that the observed separation between the two main mass components enlarges the critical area by a factor of ~2. These geometrical effects can account for the high number (density) of multiple images observed. We find in numerical simulations, that only ~4% of the clusters (with M_{vir}>6 x 10^{14} h^{-1}M_{\odot}) exhibit as elongated critical curves as M0416.

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]

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.

We evaluate the construction methodology of an all-sky catalogue of galaxy clusters detected through the Sunyaev-Zel'dovich (SZ) effect. We perform an extensive comparison of twelve algorithms applied to the same detailed simulations of the millimeter and submillimeter sky based on a Planck-like case. We present the results of this "SZ Challenge" in terms of catalogue completeness, purity, astrometric and photometric reconstruction. Our results provide a comparison of a representative sample of SZ detection algorithms and highlight important issues in their application. In our study case, we show that the exact expected number of clusters remains uncertain (about a thousand cluster candidates at |b|> 20 deg with 90% purity) and that it depends on the SZ model and on the detailed sky simulations, and on algorithmic implementation of the detection methods. We also estimate the astrometric precision of the cluster candidates which is found of the order of ~2 arcmins on average, and the photometric uncertainty of order ~30%, depending on flux.

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.

We examine the cosmological constraining power of future large-scale weak lensing surveys on the model of \emph{Euclid}, with particular reference to primordial non-Gaussianity. Our analysis considers several different estimators of the projected matter power spectrum, based on both shear and flexion, for which we review the covariances and Fisher matrices. The bounds provided by cosmic shear alone for the local bispectrum shape, marginalized over $\sigma_8$, are at the level of $\Delta f_\mathrm{NL} \sim 100$. We consider three additional bispectrum shapes, for which the cosmic shear constraints range from $\Delta f_\mathrm{NL}\sim 340$ (equilateral shape) up to $\Delta f_\mathrm{NL}\sim 500$ (orthogonal shape). The competitiveness of cosmic flexion constraints against cosmic shear ones depends on the galaxy intrinsic flexion noise, that is still virtually unconstrained. Adopting the very high value that has been occasionally used in the literature results in the flexion contribution being basically negligible with respect to the shear one, and for realistic configurations the former does not improve significantly the constraining power of the latter. Since the flexion noise decreases with decreasing scale, by extending the analysis up to $\ell_\mathrm{max} = 20,000$ cosmic flexion, while being still subdominant, improves the shear constraints by $\sim 10%$ when added. However on such small scales the highly non-linear clustering of matter and the impact of baryonic physics make any error estimation uncertain. By considering lower, and possibly more realistic, values of the flexion intrinsic shape noise results in flexion constraining power being a factor of $\sim 2$ better than that of shear, and the bounds on $\sigma_8$ and $f_\mathrm{NL}$ being improved by a factor of $\sim 3$ upon their combination. (abridged)

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.

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, \Delta, are distributed equally around zero with no preferred direction, and are well described by a log-normal distribution with <log_{10}(\Delta [h^{-1} Mpc])>=-1.895^{+0.003}_{-0.004}, and \sigma=0.501\pm0.004 (95% confidence levels), or <log_{10}(\Delta [\arcsec])>=0.564\pm0.005, and \sigma=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]

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.

We examine the inner mass distribution of the relaxed galaxy cluster Abell 383 in deep 16-band HST/ACS+WFC3 imaging taken as part of the CLASH multi-cycle treasury program. Our program is designed to study the dark matter distribution in 25 massive clusters, and balances depth with a wide wavelength coverage to better identify lensed systems and generate precise photometric redshifts. This information together with the predictive strength of our strong-lensing analysis method identifies 13 new multiply-lensed images and candidates, so that a total of 27 multiple-images of 9 systems are used to tightly constrain the inner mass profile, $d\log \Sigma/d\log r\simeq -0.6\pm 0.1$ (r<160 kpc). We find consistency with the standard distance-redshift relation for the full range spanned by the lensed images, 1.01<z<6.03, with the higher redshift sources deflected through larger angles as expected. The inner mass profile derived here is consistent with the results of our independent weak-lensing analysis of wide-field Subaru images, with good agreement in the region of overlap. The overall mass profile is well fitted by an NFW profile with M_{vir}=(5.37^{+0.70}_{-0.63}\pm 0.26) x 10^{14}M_{\odot}/h and a relatively high concentration, c_{vir}=8.77^{+0.44}_{-0.42}\pm 0.23, which lies above the standard c-M relation similar to other well-studied clusters. The critical radius of Abell 383 is modest by the standards of other lensing clusters, r_{E}\simeq16\pm2\arcsec (for z_s=2.55), so the relatively large number of lensed images uncovered here with precise photometric redshifts validates our imaging strategy for the CLASH survey. In total we aim to provide similarly high-quality lensing data for 25 clusters, 20 of which are X-ray selected relaxed clusters, enabling a precise determination of the representative mass profile free from lensing bias. (ABRIDGED)

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\Sigma/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.

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 \Lambda 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 \Delta_{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.

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.

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.

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.

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 {\mu}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.

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.

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(\theta_{e}\arcsec)>=0.73^{+0.02}_{-0.03}$, $\sigma=0.316^{+0.004}_{-0.002}$, and with higher abundance of large $\theta_{e}$ clusters than predicted by $\Lambda$CDM. We visually inspect each of the clusters with $\theta_{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 $\theta_{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.

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.

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.

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 {\Lambda}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 {\Lambda}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%.

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.

We describe the processing of the 336 billion raw data samples from the High Frequency Instrument (HFI) which we performed to produce six temperature maps from the first 295 days of Planck-HFI survey data. These maps provide an accurate rendition of the sky emission at 100, 143, 217, 353, 545 and 857 GHz with an angular resolution ranging from 9.9 to 4.4^2. The white noise level is around 1.5 {\mu}K degree or less in the 3 main CMB channels (100--217GHz). The photometric accuracy is better than 2% at frequencies between 100 and 353 GHz and around 7% at the two highest frequencies. The maps created by the HFI Data Processing Centre reach our goals in terms of sensitivity, resolution, and photometric accuracy. They are already sufficiently accurate and well-characterised to allow scientific analyses which are presented in an accompanying series of early papers. At this stage, HFI data appears to be of high quality and we expect that with further refinements of the data processing we should be able to achieve, or exceed, the science goals of the Planck project.

We present a strong-lensing analysis of the galaxy cluster MACS J1206.2-0847 ($z$=0.44) using UV, Optical, and IR, HST/ACS/WFC3 data taken as part of the CLASH multi-cycle treasury program, with VLT/VIMOS spectroscopy for some of the multiply-lensed arcs. The CLASH observations, combined with our mass-model, allow us to identify 47 new multiply-lensed images of 12 distant sources. These images, along with the previously known arc, span the redshift range $1\la z\la5.5$, and thus enable us to derive a detailed mass distribution and to accurately constrain, for the first time, the inner mass-profile of this cluster. We find an inner profile slope of $d\log \Sigma/d\log \theta\simeq -0.55\pm 0.1$ (in the range [1\arcsec, 53\arcsec], or $5\la r \la300$ kpc), as commonly found for relaxed and well-concentrated clusters. Using the many systems uncovered here we derive credible critical curves and Einstein radii for different source redshifts. For a source at $z_{s}\simeq2.5$, the critical curve encloses a large area with an effective Einstein radius of $\theta_{E}=28\pm3\arcsec$, and a projected mass of $1.34\pm0.15\times10^{14} M_{\odot}$. From the current understanding of structure formation in concordance cosmology, these values are relatively high for clusters at $z\sim0.5$, so that detailed studies of the inner mass distribution of clusters such as MACS J1206.2-0847 can provide stringent tests of the $\Lambda$CDM paradigm.

All-sky data from the Planck survey and the Meta-Catalogue of X-ray detected Clusters of galaxies (MCXC) are combined to investigate the relationship between the thermal Sunyaev-Zeldovich (SZ) signal and X-ray luminosity. The sample comprises ~ 1600 X-ray clusters with redshifts up to ~ 1 and spans a wide range in X-ray luminosity. The SZ signal is extracted for each object individually, and the statistical significance of the measurement is maximised by averaging the SZ signal in bins of X-ray luminosity, total mass, or redshift. The SZ signal is detected at very high significance over more than two decades in X-ray luminosity (10^43 erg/s < L_500 E(z)^-7/3 < 2 X 10^45 erg/s). The relation between intrinsic SZ signal and X-ray luminosity is investigated and the measured SZ signal is compared to values predicted from X-ray data. Planck measurements and X-ray based predictions are found to be in excellent agreement over the whole explored luminosity range. No significant deviation from standard evolution of the scaling relations is detected. For the first time the intrinsic scatter in the scaling relation between SZ signal and X-ray luminosity is measured and found to be consistent with the one in the luminosity -- mass relation from X-ray studies. There is no evidence of any deficit in SZ signal strength in Planck data relative to expectations from the X-ray properties of clusters, underlining the robustness and consistency of our overall view of intra-cluster medium properties.

We present the XMM-Newton follow-up for confirmation of Planck cluster candidates. Twenty-five candidates have been observed to date using snapshot (~10 ksec) exposures, ten as part of a pilot programme to sample a low range of signal-to-noise ratios (4<S/N<6), and a further 15 in a programme to observe a sample of S/N>5 candidates. The sensitivity and spatial resolution of XMM-Newton allows unambiguous discrimination between clusters and false candidates. The 4 false candidates have S/N <= 4.1. A total of 21 candidates are confirmed as extended X-ray sources. Seventeen are single clusters, the majority of which are found to have highly irregular and disturbed morphologies (about ~70%). The remaining four sources are multiple systems, including the unexpected discovery of a supercluster at z=0.45. For 20 sources we are able to derive a redshift estimate from the X-ray Fe K line (albeit of variable quality). The new clusters span the redshift range 0.09 <= z <= 0.54, with a median redshift of z~0.37. A first determination is made of their X-ray properties including the characteristic size, which is used to improve the estimate of the SZ Compton parameter, Y_SZ. The follow-up validation programme has helped to optimise the Planck candidate selection process. It has also provided a preview of the X-ray properties of these newly-discovered clusters, allowing comparison with their SZ properties, and to the X-ray and SZ properties of known clusters observed in the Planck survey. Our results suggest that Planck may have started to reveal a non-negligible population of massive dynamically perturbed objects that is under-represented in X-ray surveys. However, despite their particular properties, these new clusters appear to follow the Y_SZ-Y_X relation established for X-ray selected objects, where Y_X is the product of the gas mass and temperature.

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 $\sigma_{\epsilon}=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.

We present precise Sunyaev-Zeldovich (SZ) effect measurements in the direction of 62 nearby galaxy clusters (z <0.5) detected at high signal-to-noise in the first Planck all-sky dataset. The sample spans approximately a decade in total mass, 10^14 < M_500 < 10^15, where M_500 is the mass corresponding to a total density contrast of 500. Combining these high quality Planck measurements with deep XMM-Newton X-ray data, we investigate the relations between D_A^2 Y_500, the integrated Compton parameter due to the SZ effect, and the X-ray-derived gas mass M_g,500, temperature T_X, luminosity L_X, SZ signal analogue Y_X,500 = M_g,500 * T_X, and total mass M_500. After correction for the effect of selection bias on the scaling relations, we find results that are in excellent agreement with both X-ray predictions and recently-published ground-based data derived from smaller samples. The present data yield an exceptionally robust, high-quality local reference, and illustrate Planck's unique capabilities for all-sky statistical studies of galaxy clusters.

Based on the largest homogeneously reduced set of Type Ia supernova luminosity data currently available -- the Union2 sample -- we reconstruct the expansion history of the Universe in a model-independent approach. Our method tests the geometry of the Universe directly without reverting to any assumptions made on its energy content. This allows us to constrain Dark Energy models and non-standard cosmologies in a straightforward way. The applicability of the presented method is not restricted to testing cosmological models. It can be a valuable tool for pointing out systematic errors hidden in the supernova data and planning future Type Ia supernova cosmology campaigns.

We present the first all-sky sample of galaxy clusters detected blindly by the Planck satellite through the Sunyaev-Zeldovich (SZ) effect from its six highest frequencies. This early SZ (ESZ) sample is comprised of 189 candidates, which have a high signal-to-noise ratio ranging from 6 to 29. Its high reliability (purity above 95%) is further ensured by an extensive validation process based on Planck internal quality assessments and by external cross-identification and follow-up observations. Planck provides the first measured SZ signal for about 80% of the 169 previously-known ESZ clusters. Planck furthermore releases 30 new cluster candidates, amongst which 20 meet the ESZ signal-to-noise selection criterion. At the submission date, twelve of the 20 ESZ candidates were confirmed as new clusters, with eleven confirmed using XMM-Newton snapshot observations, most of them with disturbed morphologies and low luminosities. The ESZ clusters are mostly at moderate redshifts (86% with z below 0.3) and span more than a decade in mass, up to the rarest and most massive clusters with masses above 10^15 Msol.

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.

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.

A long-standing problem of strong lensing by galaxy clusters regards the observed high rate of giant gravitational arcs as compared to the predictions in the framework of the "standard" cosmological model. Recently, few other inconsistencies between theoretical expectations and observations have been claimed which regard the large size of the Einstein rings and the high concentrations of few clusters with strong lensing features. All of these problems consistently indicate that observed galaxy clusters may be gravitational lenses stronger than expected. We use clusters extracted from the MareNostrum Universe to build up mock catalogs of galaxy clusters selected through their X-ray flux. We use these objects to estimate the probability distributions of lensing cross sections, Einstein rings, and concentrations for the sample of 12 MACS clusters at $z>0.5$ presented in Ebeling et al. (2007) and discussed in Zitrin et al. (2010). We find that simulated clusters produce $\sim 50%$ less arcs than observed clusters do. The medians of the distributions of the Einstein ring sizes differ by $\sim 25%$ between simulations and observations. We estimate that, due to cluster triaxiality and orientation biases affecting the lenses with the largest cross sections, the concentrations of the individual MACS clusters inferred from the lensing analysis should be up to a factor of $\sim 2$ larger than expected from the $\Lambda$CDM model. The arc statistics, the Einstein ring, and the concentration problems in strong lensing clusters are mitigated but not solved on the basis of our analysis. Nevertheless, due to the lack of redshifts for most of the multiple image systems used for modeling the MACS clusters, the results of this work will need to be verified with additional data. The upcoming CLASH program will provide an ideal sample for extending our comparison (abridged).

We present the Sunyaev-Zeldovich (SZ) signal-to-richness scaling relation (Y500-N200) for the MaxBCG cluster catalogue. Employing a multi-frequency matched filter on the Planck sky maps, we measure the SZ signal for each cluster by adapting the filter according to weak-lensing calibrated mass-richness relations (N200-M500). We bin our individual measurements and detect the SZ signal down to the lowest richness systems (N200=10) with high significance, achieving a detection of the SZ signal in systems with mass as low as M500~5e13 Msolar. The observed Y500-N200 relation is well modeled by a power law over the full richness range. It has a lower normalisation at given N200 than predicted based on X-ray models and published mass-richness relations. An X-ray subsample, however, does conform to the predicted scaling, and model predictions do reproduce the relation between our measured bin-average SZ signal and measured bin-average X-ray luminosities. At fixed richness, we find an intrinsic dispersion in the Y500-N200 relation of 60% rising to of order 100% at low richness. Thanks to its all-sky coverage, Planck provides observations for more than 13,000 MaxBCG clusters and an unprecedented SZ/optical data set, extending the list of known cluster scaling laws to include SZ-optical properties. The data set offers essential clues for models of galaxy formation. Moreover, the lower normalisation of the SZ-mass relation implied by the observed SZ-richness scaling has important consequences for cluster physics and cosmological studies with SZ clusters.

The Planck High Frequency Instrument (HFI) is designed to measure the temperature and polarization anisotropies of the Cosmic Microwave Background and galactic foregrounds in six wide bands centered at 100, 143, 217, 353, 545 and 857 GHz at an angular resolution of 10' (100 GHz), 7' (143 GHz), and 5' (217 GHz and higher). HFI has been operating flawlessly since launch on 14 May 2009. The bolometers cooled to 100 mK as planned. The settings of the readout electronics, such as the bolometer bias current, that optimize HFI's noise performance on orbit are nearly the same as the ones chosen during ground testing. Observations of Mars, Jupiter, and Saturn verified both the optical system and the time response of the detection chains. The optical beams are close to predictions from physical optics modeling. The time response of the detection chains is close to pre-launch measurements. The detectors suffer from an unexpected high flux of cosmic rays related to low solar activity. Due to the redundancy of Planck's observations strategy, the removal of a few percent of data contaminated by glitches does not affect significantly the sensitivity. The cosmic rays heat up significantly the bolometer plate and the modulation on periods of days to months of the heat load creates a common drift of all bolometer signals which do not affect the scientific capabilities. Only the high energy cosmic rays showers induce inhomogeneous heating which is a probable source of low frequency noise.

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.

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.

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.

We analyse the KSB method to estimate gravitational shear from surface-brightness moments of small and noisy galaxy images. We identify three potentially problematic assumptions. These are: (1) While gravitational shear must be estimated from averaged galaxy images, KSB derives a shear estimate from each individual image and then takes the average. Since the two operations do not commute, KSB gives biased results. (2) KSB implicitly assumes that galaxy ellipticities are small, while weak gravitational lensing assures only that the change in ellipticity due to the shear is small. (3) KSB does not invert the convolution with the point-spread function, but gives an approximate PSF correction which - even for a circular PSF - holds only in the limit of circular sources. The effects of assumptions (2) and (3) partially counter-act in a way dependent on the width of the weight function and of the PSF. We quantitatively demonstrate the biases due to all assumptions, extend the KSB approach consistently to third order in the shear and ellipticity and show that this extension lowers the biases substantially. The issue of proper PSF deconvolution will be addressed in a forthcoming paper.

We study the spherical collapse model for several dark energy scenarios using the fully nonlinear differential equation for the evolution of the density contrast within homogeneous spherical overdensities derived from Newtonian hydrodynamics. While mathematically equivalent to the more common approach based on the differential equation for the radius of the perturbation, this approach has substantial conceptual as well as numerical advantages. Among the most important are that no singularities at early times appear, which avoids numerical problems in particular in applications to cosmologies with dynamical and early dark energy, and that the assumption of time-reversal symmetry can easily be dropped where it is not strictly satisfied. We use this approach to derive the two parameters characterising the spherical-collapse model, i.e.~the linear density threshold for collapse $\delta_\mathrm{c}$ and the virial overdensity $\Delta_\mathrm{V}$, for a broad variety of dark-energy models and to reconsider these parameters in cosmologies with early dark energy. We find that, independently of the model under investigation, $\delta_\mathrm{c}$ and $\Delta_\mathrm{V}$ are always very close to the values obtained for the standard $\Lambda$CDM model, arguing that the abundance of and the mean density within non-linear structures are quite insensitive to the differences between dark-energy cosmologies. Regarding early dark energy, we thus arrive at a different conclusion than some earlier papers, including one from our group, and we explain why.

Aims. We calibrate the number density, completeness, reliability and the lower mass limit of galaxy-cluster detections through their thermal SZ signal, and compare them to X-ray cluster detections. Methods. We simulate maps of the thermal SZ effect and the X-ray emission from light cones constructed in a large, hydrodynamical, cosmological simulation volume, including realistic noise contributions. The maps are convolved with linear, optimised, single- and multi-band filters to identify local peaks and their signal-to-noise ratios. The resulting peak catalogues are then compared to the halo population in the simulation volume to identify true and spurious detections. Results. Multi-band filtering improves the statistics of SZ cluster detections considerably compared to single-band filtering. Observations with the characteristics of ACT detect clusters with masses M>6-9e13 M_o/h, quite independent of redshift, reach 50% completeness at ~1e14 M_o/h and 100% completeness at ~2e14 M_o/h. Samples are contaminated by only a few per cent spurious detections. This is broadly comparable to X-ray cluster detections with XMM-Newton with 100 ks exposure time in the soft band, except that the mass limit for X-ray detections increases much more steeply with redshift than for SZ detections. A comparison of true and filtered signals in the SZ and X-ray maps confirms that the filters introduce at most a negligible bias.

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.

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.

We construct two linear filtering techniques based on weak gravitational lensing to constrain the inner slope $\alpha$ of the density profile of dark-matter halos. Both methods combine all available information into an estimate of this single number. Under idealised assumptions, $\alpha$ is constrained to ~15% if the halo concentration c is known, and to < 30% if not. We argue that the inevitable degeneracy between density-profile slope and halo concentration cannot be lifted under realistic conditions, and show by means of Fisher-matrix methods which linear combination of $\alpha$ and c is best constrained by our filtering of lensing data. This defines a new parameter, called P1, which is then constrained to ~15% for a single massive halo. If the signals of many halos can be stacked, their density profiles should thus be well constrained by the linear filters proposed here with the advantagebe insensitive to the cluster substructures.

This paper provides an overview of the Low Frequency Instrument (LFI) programme within the ESA Planck mission. The LFI instrument has been developed to produce high precision maps of the microwave sky at frequencies in the range 27-77 GHz, below the peak of the cosmic microwave background (CMB) radiation spectrum. The scientific goals are described, ranging from fundamental cosmology to Galactic and extragalactic astrophysics. The instrument design and development are outlined, together with the model philosophy and testing strategy. The instrument is presented in the context of the Planck mission. The LFI approach to ground and inflight calibration is described. We also describe the LFI ground segment. We present the results of a number of tests demonstrating the capability of the LFI data processing centre (DPC) to properly reduce and analyse LFI flight data, from telemetry information to calibrated and cleaned time ordered data, sky maps at each frequency (in temperature and polarization), component emission maps (CMB and diffuse foregrounds), catalogs for various classes of sources (the Early Release Compact Source Catalogue and the Final Compact Source Catalogue). The organization of the LFI consortium is briefly presented as well as the role of the core team in data analysis and scientific exploitation. All tests carried out on the LFI flight model demonstrate the excellent performance of the instrument and its various subunits. The data analysis pipeline has been tested and its main steps verified. In the first three months after launch, the commissioning, calibration, performance, and verification phases will be completed, after which Planck will begin its operational life, in which LFI will have an integral part.

We develop and apply an analytic method to predict peak counts in weak-lensing surveys. It is based on the theory of Gaussian random fields and suitable to quantify the level of spurious detections caused by chance projections of large-scale structures as well as the shape and shot noise contributed by the background galaxies. We compare our method to peak counts obtained from numerical ray-tracing simulations and find good agreement at the expected level. The number of peak detections depends substantially on the shape and size of the filter applied to the gravitational shear field. Our main results are that weak-lensing peak counts are dominated by spurious detections up to signal-to-noise ratios of 3--5 and that most filters yield only a few detections per square degree above this level, while a filter optimised for suppressing large-scale structure noise returns up to an order of magnitude more.

We seek to understand the impact on shape estimators obtained from circular and elliptical shapelet models under two realistic conditions: (a) only a limited number of shapelet modes is available for the model, and (b) the intrinsic galactic shapes are not restricted to shapelet models. We create a set of simplistic simulations, in which the galactic shapes follow a Sersic profile. By varying the Sersic index and applied shear, we quantify the amount of bias on shear estimates which arises from insufficient modeling. Additional complications due to PSF convolution, pixelation and pixel noise are also discussed. Steep and highly elliptical galaxy shapes cannot be accurately modeled within the circular shapelet basis system and are biased towards shallower and less elongated shapes. This problem can be cured partially by allowing elliptical basis functions, but for steep profiles elliptical shapelet models still depend critically on accurate ellipticity priors. As a result, shear estimates are typically biased low. Independently of the particular form of the estimator, the bias depends on the true intrinsic galaxy morphology, but also on the size and shape of the PSF. As long as the issues discussed here are not solved, the shapelet method cannot provide weak-lensing measurements with an accuracy demanded by upcoming missions and surveys, unless one can provide an accurate and reliable calibration, specific for the dataset under investigation.

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.

AIMS. While weak lensing cannot resolve cluster cores and strong lensing is almost insensitive to density profiles outside the scale radius, combinations of both effects promise to constrain density profiles of galaxy clusters well, and thus to allow testing of the CDM expectation on dark-matter halo density profiles. METHODS. We develop an algorithm further that we had recently proposed for this purpose. It recovers a lensing potential optimally reproducing observations of both strong and weak-lensing effects by combining high resolution in cluster cores with the larger-scale information from weak lensing. The main extensions concern the accommodation of mild non-linearity in inner iterations, the progressive increase in resolution in outer iterations, and the introduction of a suitable regularisation term. The linearity of the method is essentially preserved. RESULTS. We demonstrate the success of the algorithm with both idealised and realistic simulated data, showing that the simulated lensing mass distribution and its density profile are well reproduced. We then apply it to weak and strong lensing data of the cluster MS 2137 and obtain a parameter-free solution which is in good qualitative agreement with earlier parametric studies.

We study the clustering properties of galaxy clusters expected to be observed by various forthcoming surveys both in the X-ray and sub-mm regimes by the thermal Sunyaev-Zel'dovich effect. Several different background cosmological models are assumed, including the concordance $\Lambda$CDM and various cosmologies with dynamical evolution of the dark energy. Particular attention is paid to models with a significant contribution of dark energy at early times which affects the process of structure formation. Past light cone and selection effects in cluster catalogs are carefully modeled by realistic scaling relations between cluster mass and observables and by properly taking into account the selection functions of the different instruments. The results show that early dark-energy models are expected to produce significantly lower values of effective bias and both spatial and angular correlation amplitudes with respect to the standard $\Lambda$CDM model. Among the cluster catalogues studied in this work, it turns out that those based on \emph{eRosita}, \emph{Planck}, and South Pole Telescope observations are the most promising for distinguishing between various dark-energy models.

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 $\psi_E$ and $\psi_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.

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.

Context. One science goal of the upcoming Planck mission is to perform a full-sky cluster survey based on the Sunyaev-Zel'dovich (SZ) effect, which leads to the question of how such a survey would be affected by cosmological models with a different history of structure formation than LCDM. One class of these models are early dark energy (EDE) cosmologies, where the dark energy contribution does not vanish at early times. Aims. Since structures grow slower in the presence of EDE, one expects an increase in the number of galaxy clusters compared to LCDM at intermediate and high redshifts, which could explain the reported excess of the angular CMB power spectrum on cluster scales via an enhanced SZ contribution. We study the impact of EDE on Planck's expected cluster sample. Methods. To obtain realistic simulations, we constructed full-sky SZ maps for EDE and LCDM cosmologies, taking angular cluster correlation into account. Using these maps, we simulated Planck observations with and without Galactic foregrounds and fed the results into our filter pipeline based on the spherical multi-frequency matched filters. Results. For the case of EDE cosmologies, we clearly find an increase in the detected number of clusters compared to the fiducial LCDM case. This shows that the spherical multi-frequency matched filter is sensitive enough to find deviations from the LCDM sample, being caused by EDE. In addition we find an interesting effect of EDE on the completeness of the cluster sample, such that EDE helps to obtain cleaner samples.

Previous works have demonstrated that the generation of secondary CMB anisotropies due to the molecular optical depth is likely too small to be observed. In this paper, we examine additional ways in which primordial chemistry and the dark ages might influence the CMB. We present a detailed and updated chemical network and give an overview of the interactions of molecules with the CMB. We consider the optical depth due to line absorption, photoionization, photodissociation and free-free processes, and estimate the resulting changes in the CMB temperature and its power spectrum. The most promising results are obtained for the negative hydrogen ion \HM and the \HeHII molecule. The free-free process of \HM yields a relative change in the CMB temperature of up to $2\times10^{-11}$, and leads to a frequency-dependent change in the power spectrum of the order $10^{-7}$ at 30 GHz. With a change of the order $10^{-10}$ in the power spectrum, our result for the bound-free process of \HM is significantly below a previous suggestion. \HeHII efficiently scatters CMB photons and smears out primordial fluctuations, leading to a change in the power spectrum of the order $10^{-8}$. Improvements in the accuracy of future CMB experiments may thus help to constrain and finally detect these interesting signals from the dark ages of the universe.

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.

We use the large cosmological Millennium Simulation (MS) to construct the first all-sky maps of the lensing potential and the deflection angle, aiming at gravitational lensing of the CMB, with the goal of properly including small-scale non-linearities and non-Gaussianity. Exploiting the Born approximation, we implement a map-making procedure based on direct ray-tracing through the gravitational potential of the MS. We stack the simulation box in redshift shells up to $z\sim 11$, producing continuous all-sky maps with arcminute angular resolution. A randomization scheme avoids repetition of structures along the line of sight and structures larger than the MS box size are added to supply the missing contribution of large-scale (LS) structures to the lensing signal. The angular power spectra of the projected lensing potential and the deflection-angle modulus agree quite well with semi-analytic estimates on scales down to a few arcminutes, while we find a slight excess of power on small scales, which we interpret as being due to non-linear clustering in the MS. Our map-making procedure, combined with the LS adding technique, is ideally suited for studying lensing of CMB anisotropies, for analyzing cross-correlations with foreground structures, or other secondary CMB anisotropies such as the Rees-Sciama effect.

We use semi-analytic modelling of the galaxy-cluster population and its strong lensing efficiency to explore how the expected abundance of large gravitational arcs on the sky depends on $\sigma_8$. Our models take all effects into account that have been shown to affect strong cluster lensing substantially, in particular cluster asymmetry, substructure, merging, and variations in the central density concentrations. We show that the optical depth for long and thin arcs increases by approximately one order of magnitude when $\sigma_8$ increases from 0.7 to 0.9, owing to a constructive combination of several effects. Models with high $\sigma_8$ are also several orders of magnitude more efficient in producing arcs at intermediate and high redshifts. Finally, we use realistic source number counts to quantitatively predict the total number of arcs brighter than several magnitude limits in the R and I bands. We confirm that, while $\sigma_8\sim0.9$ may come close to the known abundance of arcs, even $\sigma_8\sim0.8$ falls short by almost an order of magnitude in reproducing known counts. We conclude that, should $\sigma_8\sim0.8$ be confirmed, we would fail to understand the strong-lensing efficiency of the galaxy cluster population, and in particular the abundance of arcs in high-redshift clusters. We argue that early-dark energy or non-Gaussian density fluctuations may indicate one way out of this problem.

While both isothermal and NFW-based mass models for galaxy clusters are widely adopted in strong-lensing studies, they cannot easily be distinguished based solely on observed positions of arcs and arclets. We compare the magnifications predicted for giant arcs obtained from isothermal and NFW profiles, taking axially-symmetric and asymmetric mass distributions into account. We find that arc magnifications can differ strongly between the two types of density profiles even if the image morphology is well reproduced. Magnifications by lenses with NFW density profiles are usually larger than those for lenses with singular or nearly singular isothermal density profiles, unless the latter have large cores. Asymmetries play an important role. We illustrate our results with the two well-studied clusters MS 2137 and A~370. We confirm earlier results showing that both isothermal and NFW mass models can very well reproduce the observed arcs, radial arcs and other arclets. While the mass model for MS 2137 is not very well constrained, the two types of mass models produce strongly differing critical curves and caustics for A~370. We find that the NFW mass profile is preferred for A~370. We identify new arclet candidates in the field of A~370. Redshift estimates allowed by the lens model are consistent with results in the literature, where available. Three newly found counter-images are suggested to arise from an active, dust-enshrouded star-forming galaxy at $z\approx1.1$.

Aims: In view of the substantial uncertainties regarding the possible dynamics of the dark energy, we aim at constraining the expansion rate of the universe without reference to a specific Friedmann model and its parameters. Methods: We show that cosmological observables integrating over the cosmic expansion rate can be converted into a Volterra integral equation which is known to have a unique solution in terms of a Neumann series. Expanding observables such as the luminosity distances to type-Ia supernovae into a series of orthonormal functions, the integral equation can be solved and the cosmic expansion rate recovered within the limits allowed by the accuracy of the data. Results: We demonstrate the performance of the method applying it to synthetic data sets of increasing complexity, and to the first-year SNLS data. In particular, we show that the method is capable of reproducing a hypothetical expansion function containing a sudden transition.

In this paper we investigate the effects that a dynamic dark energy component dominant in the universe at late epochs has on reionization. We follow the evolution of HII regions with the analytic approach of Furlanetto and Oh (2005) in two different universes for which we assume the Peebles and Ratra (2003) and Brax and Martin (2000) quintessence models and we compare our results to the LCDM scenario. We show that, for a fixed ionization efficiency, at the same cosmological epoch the topology of bubbles is dominated by high-mass objects and the characteristic size of the ionized regions is slightly smaller than in the LCDM model, especially at the latest stages of reionization, due to the higher recombination efficiency. As a consequence, the bubbles' `epoch of overlap' happens earlier than in LCDM. Finally, we show how the different evolution of the HII regions affects the transmission of the high-z QSO spectra, reducing the Lyman flux absorption at small optical depths.

We present a newly developed code that allows simulations of optical observations of galaxy fields with a variety of instruments. The code incorporates gravitational lensing effects and is targetted at simulating lensing by galaxy clusters. Our goal is to create the tools required for comparing theoretical expectations with observations to obtain a better understanding of how observational noise affects lensing applications such as mass estimates, studies on the internal properties of galaxy clusters and arc statistics. Starting from a set of input parameters, characterizing both the instruments and the observational conditions, the simulator provides a virtual observation of a patch of the sky. It includes several sources of noise such as photon-noise, sky background, seeing, and instrumental noise. Ray-tracing through simulated mass distributions accounts for gravitational lensing. Source morphologies are realistically simulated based on shapelet decompositions of galaxy images retrieved from the GOODS-ACS archive. According to their morphological class, spectral-energy-distributions are assigned to the source galaxies in order to reproduce observations of each galaxy in arbitrary photometric bands. We illustrate our techniques showing virtual observations of a galaxy-cluster core as it would be observed with the space telescope DUNE, which was recently proposed to ESA within its "Cosmic vision" programme. (Abridged)

Aims: We explore the lensing properties of asymmetric matter density distributions in Bekenstein's Tensor-Vector-Scalar theory (TeVeS). Methods: Using an iterative Fourier-based solver for the resulting non-linear scalar field equation, we numerically calculate the total gravitational potential and derive the corresponding TeVeS lensing maps. Results: Considering variations on rather small scales, we show that the lensing properties significantly depend on the lens's extent along the line of sight. Furthermore, all simulated TeVeS convergence maps strongly track the dominant baryonic components, non-linear effects, being capable of counteracting this trend, turn out to be very small. Setting up a toy model for the cluster merger 1E0657-558, we infer that TeVeS cannot explain observations without assuming an additional dark mass component in both cluster centers, which is in accordance with previous work.

We study the reliability of dark-matter halo detections with three different linear filters applied to weak-lensing data. We use ray-tracing in the multiple lens-plane approximation through a large cosmological simulation to construct realizations of cosmic lensing by large-scale structures between redshifts zero and two. We apply the filters mentioned above to detect peaks in the weak-lensing signal and compare them with the true population of dark matter halos present in the simulation. We confirm the stability and performance of a filter optimized for suppressing the contamination by large-scale structure. It allows the reliable detection of dark-matter halos with masses above a few times 1e13 M_sun/h with a fraction of spurious detections below ~10%. For sources at redshift two, 50% of the halos more massive than ~7e13 M_sun/h are detected, and completeness is reached at ~2e14 M_sun/h.

Galaxy clusters are hotter and more X-ray luminous than in quiescence while they undergo major mergers, which also transiently increase their strong-lensing efficiency. We use semi-analytic models for both effects to study how cluster dynamics in different dark-energy models affects the X-ray selected cluster population and its strong-lensing optical depth. We find that mergers increase the number of observable X-ray clusters by factors of a few and considerably broaden their redshift distribution. Strong-lensing optical depths are increased by a very similar amount. Quite independent of cosmology, X-ray bright clusters above a flux limit of $10^{-13.5} \mathrm{erg s^{-1} cm^{-2}}$ produce $\sim60%$ of the strong-lensing optical depth, and only $\sim1%$ above a flux limit of $10^{-11.5} \mathrm{erg s^{-1} cm^{-2}}$ if mergers are taken into account.

We use simulated merger trees of galaxy-cluster halos to study the effect of the halo concentration distribution on strong lensing and X-ray emission. Its log-normal shape typically found in simulations favors outliers with high concentration. Since, at fixed mass, more concentrated halos tend to be more efficient lenses, the scatter in the concentration increases the strong-lensing optical depth by $\lesssim50%$. Within cluster samples, mass and concentration have counteracting effects on strong lensing and X-ray emission because the concentration decreases for increasing mass. Selecting clusters by concentration thus has no effect on the lensing cross section. The most efficiently lensing and hottest clusters are typically the \textit{least} concentrated in samples with a broad mass range. Among cluster samples with a narrow mass range, however, the most strongly lensing and X-ray brightest clusters are typically 10% to 25% more concentrated.

We address the question of how well the density profile of galaxy clusters can be determined by combining strong lensing and velocity dispersion data. We use cosmological dark matter simulations of clusters to test the reliability of the method, producing mock catalogues of tangential and radial gravitational arcs and simulating the radial velocity dispersion profile of the cluster brightest central galaxy. The density profiles of the simulated clusters closely follow the NFW form, but we find that the recovered values of the inner slope are systematically underestimated, by about 0.4 in the mean, if the lens is assumed to be axially symmetric. However, if the ellipticity and orientation of the iso-contours of the cluster lensing potential are taken into account, then the inner slopes can be recovered quite accurately for a significant subset of the clusters whose central surface density profiles appear the most regular. These have lensing potentials with ellipticities in the range 0.15-0.4. Further simulations projecting one cluster along many random lines-of-sight show that, even for lower ellipticities, the central slopes are underestimated by $\sim 10-35%$. These simulations closely mimic past observations (see e.g. Sand et al., 2004), suggesting that existing estimates of the central slopes may be biased towards low values. For the remaining clusters, where the lensing potential is strongly perturbed by active merging or by substructure, the correct determination of the inner slope requires a more accurate model for the lens. When the halo profile is modelled by a generalised NFW profile, we find that the inferred scale radius and characteristic density, unlike the inner slope, are generally poorly constrained, since there is a strong degeneracy between these two parameters.

We show how hydrostatic equilibrium in galaxy clusters can be quantitatively probed combining X-ray, SZ, and gravitational-lensing data. Our previously published method for recovering three-dimensional cluster gas distributions avoids the assumption of hydrostatic equilibrium. Independent reconstructions of cumulative total-mass profiles can then be obtained from the gas distribution, assuming hydrostatic equilibrium, and from gravitational lensing, neglecting it. Hydrostatic equilibrium can then be quantified comparing the two. We describe this procedure in detail and show that it performs well on progressively realistic synthetic data. An application to a cluster merger demonstrates how hydrostatic equilibrium is violated and restored as the merger proceeds.

Aims: We discuss the applicability and reliability of the shapelet technique for scientific image analysis. Methods: We quantify the effects of non-orthogonality of sampled shapelet basis functions and misestimation of shapelet parameters. We perform the shapelet decomposition on artificial galaxy images with underlying shapelet models and galaxy images from the GOODS survey, comparing the publicly available IDL implementation with our new C++ implementation. Results: Non-orthogonality of the sampled basis functions and misestimation of the shapelet parameters can cause substantial misinterpretation of the physical properties of the decomposed objects. Additional constraints, image preprocessing and enhanced precision have to be incorporated in order to achieve reliable decomposition results.

We use the semi-analytic method developed by Fedeli et al. for computing strong-lensing optical depths to study the statistics of gravitational arcs in four dark-energy cosmologies. Specifically, we focus on models with early dark energy and compare them to more conventional models. Merger trees are constructed for the cluster population because strong cluster lensing is amplified by factors of two to three during mergers. We find that the optical depth for gravitational arcs in the early dark-energy models is increased by up to a factor of about 3 compared to the other models because of the modified dynamics of cluster formation. In particular, the probability for gravitational arcs in high-redshift clusters is considerably increased, which may offer an explanation for the unexpectedly high lensing efficiency of distant clusters.

We present the first hydrodynamic N-body simulations of primordial gas clouds responsible for the reionisation process in dark energy cosmologies. We compare the cosmological constant scenario with a SUGRA quintessence model with marked dynamics in order to highlight effects due to the different acceleration histories imposed by the dark energy. We show that both the number density of gas clouds and their clumpiness keep a record of the expansion rate during evolution, similar to the non-linear dark matter profile at virialisation, as was recently demonstrated by Dolag et al. (2004). Varying the shape of the primordial power spectrum, we show how this effect is mitigated by a running spectral index decreasing the power at small scales. Our results demonstrate that, in order to constrain the dark energy from large scale structures, one must track its effects down to the distribution of luminous matter.

We propose a method for recovering the two-dimensional gravitational potential of galaxy clusters which combines data from weak and strong gravitational lensing. A first estimate of the potential from weak lensing is improved at the approximate locations of critical curves. The method can be fully linearised and does not rely on the existence and identification of multiple images. We use simulations to show that it recovers the surface-mass density profiles and distributions very accurately, even if critical curves are only partially known and if their location is realistically uncertain. We further describe how arcs at different redshifts can be combined, and how deviations from weak lensing can be included.

We propose a new smoothing method for obtaining surface densities from discrete particle positions from numerical simulations. This is an essential step for many applications in gravitational lensing. This method is based on the ``scatter'' interpretation of the discrete density field in the Smoothed Particle Hydrodynamics. We use Monte Carlo simulations of uniform density fields and one isothermal ellipsoid to empirically derive the noise properties, and best smoothing parameters (such as the number of nearest neighbors used). A cluster from high-resolution simulations is then used to assess the reality of high-order singularities such as swallowtails and butterflies in caustics, which are important for the interpretation of substructures in gravitational lenses. We also compare our method with the Delaunay tesselation field estimator using the galaxy studied by Bradac et al. (2004), and find good agreements. We show that higher order singularities are not only connected with bound subhaloes but also with the satellite streams. However, the presence of high-order singularities are sensitive to not only the fluctuation amplitude of the surface density, but also the detailed form of the underlying smooth lensing potential (such as ellipticity and external shear).

We investigate how ellipticity, asymmetries and substructures separately affect the ability of galaxy clusters to produce strong lensing events, i.e. gravitational arcs, and how they influence the arc morphologies and fluxes. This is important for those studies aiming, for example, at constraining cosmological parameters from statistical lensing, or at determining the inner structure of galaxy clusters through gravitational arcs. We do so by creating two-dimensional gradually smoothed, differently elliptical and asymmetric versions of some numerical models. On average, we find that the contributions of ellipticity, asymmetries and substructures amount to ~40%, ~10% and ~30% of the total strong lensing cross section, respectively. However, our analysis shows that substructures play a more important role in less elliptical and asymmetric clusters, even if located at large distances from the cluster centers (~1Mpc/h). Conversely, their effect is less important in highly asymmetric lenses. The morphology, position and flux of individual arcs are strongly affected by the presence of substructures in the clusters. Removing substructures on spatial scales <~50kpc/h, roughly corresponding to mass scales <~5 10^{10}M_\odot/h, alters the image multiplicity of ~35% of the sources used in the simulations and causes position shifts larger than 5'' for ~40% of the arcs longer than 5''. We conclude that any model for cluster lens cannot neglect the effects of ellipticity, asymmetries and substructures. On the other hand, the high sensitivity of gravitational arcs to deviations from regular, smooth and symmetric mass distributions suggests that strong gravitational lensing is potentially a powerfull tool to measure the level of substructures and asymmetries in clusters.

We present an efficient algorithm designed for and capable of detecting elongated, thin features such as lines and curves in astronomical images, and its application to the automatic detection of gravitational arcs. The algorithm is sufficiently robust to detect such features even if their surface brightness is near the pixel noise in the image, yet the amount of spurious detections is low. The algorithm subdivides the image into a grid of overlapping cells which are iteratively shifted towards a local centre of brightness in their immediate neighbourhood. It then computes the ellipticity for each cell, and combines cells with correlated ellipticities into objects. These are combined to graphs in a next step, which are then further processed to determine properties of the detected objects. We demonstrate the operation and the efficiency of the algorithm applying it to HST images of galaxy clusters known to contain gravitational arcs. The algorithm completes the analysis of an image with 3000x3000 pixels in about 4 seconds on an ordinary desktop PC. We discuss further applications, the method's remaining problems and possible approaches to their solution.

We apply the linear filter for the weak-lensing signal of dark-matter halos developed in Maturi et al. (2005) to the cosmic-shear data extracted from the Garching-Bonn-Deep-Survey (GaBoDS). We wish to search for dark-matter halos through weak-lensing signatures which are significantly above the random and systematic noise level caused by intervening large-scale structures. We employ a linear matched filter which maximises the signal-to-noise ratio by minimising the number of spurious detections caused by the superposition of large-scale structures (LSS). This is achieved by suppressing those spatial frequencies dominated by the LSS contamination. We confirm the improved stability and reliability of the detections achieved with our new filter compared to the commonly-used aperture mass (Schneider, 1996; Schneider et al., 1998) and to the aperture mass based on the shear profile expected for NFW haloes (see e.g. Schirmer et al., 2004; Hennawi & Spergel, 2005). Schirmer et al.~(2006) achieved results comparable to our filter, but probably only because of the low average redshift of the background sources in GaBoDS, which keeps the LSS contamination low. For deeper data, the difference will be more important, as shown by Maturi et al. (2005). We detect fourteen halos on about eighteen square degrees selected from the survey. Five are known clusters, two are associated with over-densities of galaxies visible in the GaBoDS image, and seven have no known optical or X-ray counterparts.

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.

We propose and test a new method based on Richardson-Lucy deconvolution to reconstruct three-dimensional gas density and temperature distributions in galaxy clusters from combined X-ray and thermal Sunyaev-Zel'dovich observations. Clusters are assumed to be axially symmetric and arbitrarily inclined with respect to the line-of-sight. No equilibrium assumption other than local thermal equilibrium is needed. We test the algorithm with synthetic observations of analytically modeled and numerically simulated galaxy clusters and discuss the quality of the density and temperature reconstructions in idealised situations and in presence of observational noise, deviations from axial symmetry and cluster substructure. We find that analytic and numerical gas density and temperature distributions can be accurately reconstructed in three dimensions, even if observational noise is present. We also discuss methods for determining the inclination angle from data and show that it can be constrained using X-ray temperature maps. For a realistic cluster and including observational noise the three-dimensional reconstructions reach a level of accuracy of about 15%.

Previous studies of strong gravitational lensing by galaxy clusters neglected the potential impact of the intracluster gas. Here, we compare simulations of strong cluster lensing including gas physics at increasing levels of complexity, i.e. with adiabatic, cooling, star-forming, feedback-receiving, and thermally conducting gas, and with different implementations of the artificial viscosity in the SPH simulations. Each cluster was simulated starting from the same initial conditions such as to allow directly comparing the simulated clusters. We compare the clusters' shapes, dynamics and density profiles and study their strong-lensing cross sections computed by means of ray-tracing simulations. We find that the impact of adiabatic gas depends on the amount of turbulence that builds up, which means that the artificial viscosity plays an important role. With the common viscosity implementation, adiabatic gas has little effect on strong cluster lensing, while lower viscosity allows stronger turbulence, thus higher non-thermal pressure and a generally broader gas distribution which tends to lower lensing cross sections. Conversely, cooling and star formation steepen the core density profiles and can thus increase the strong-lensing efficiency considerably.

Although N-body simulations of cosmic structure formation suggest that dark matter halos have density profiles shallower than isothermal at small radii and steeper at large radii, whether observed galaxy clusters follow this profile is still ambiguous. We use one such density profile, the asymmetric NFW profile, to model the mass distributions of 11 galaxy clusters with gravitational arcs observed by HST. We characterize the galaxy lenses in each cluster as NFW ellipsoids, each defined by an unknown scale convergence, scale radius, ellipticity, and position angle. For a given set of values of these parameters, we compute the arcs that would be produced by such a lens system. To define the goodness of fit to the observed arc system, we define a chi^2 function encompassing the overlap between the observed and reproduced arcs as well as the agreement between the predicted arc sources and the observational constraints on the source system. We minimize this chi^2 to find the values of the lens parameters that best reproduce the observed arc system in a given cluster. Here we report our best-fit lens parameters and corresponding mass estimates for each of the 11 lensing clusters. We find that cluster mass models based on lensing galaxies defined as NFW ellipsoids can accurately reproduce the observed arcs, and that the best-fit parameters to such a model fall within the reasonable ranges defined by simulations. These results assert NFW profiles as an effective model for the mass distributions of observed clusters.

Strong gravitational lensing by irregular mass distributions, such as galaxy clusters, is generally not well quantified by cross sections of analytic mass models. Computationally expensive ray-tracing methods have so far been necessary for accurate cross-section calculations. We describe a fast, semi-analytic method here which is based on surface integrals over high-magnification regions in the lens plane and demonstrate that it yields reliable cross sections even for complex, asymmetric mass distributions. The method is faster than ray-tracing simulations by factors of $\sim30$ and thus suitable for large cosmological simulations, saving large amounts of computing time. We apply this method to a sample of galaxy cluster-sized dark matter haloes with simulated merger trees and show that cluster mergers approximately double the strong-lensing optical depth for lens redshifts $z_\mathrm{l}\gtrsim0.5$ and sources near $z_\mathrm{s} = 2$. We believe that this result hints at one possibility for understanding the recently detected high arcs abundance in clusters at moderate and high redshifts, and is thus worth further studies.

We describe an assembly of numerical tools to model the output data of the Planck satellite. These start with the generation of a CMB sky in a chosen cosmology, add in various foreground sources, convolve the sky signal with arbitrary, even non-symmetric and polarised beam patterns, derive the time ordered data streams measured by the detectors depending on the chosen satellite-scanning strategy, and include noise signals for the individual detectors and electronic systems. The simulation products are needed to develop, verify, optimise, and characterise the accuracy and performance of all data processing and scientific analysis steps of the Planck mission, including data handling, data integrity checking, calibration, map making, physical component separation, and power spectrum estimation. In addition, the simulations allow detailed studies of the impact of many stochastic and systematic effects on the scientific results. The efficient implementation of the simulation allows the build-up of extended statistics of signal variances and co-variances. Although being developed specifically for the Planck mission, it is expected that the employed framework as well as most of the simulation tools will be of use for other experiments and CMB-related science in general.

We compare the statistical properties of giant gravitationally lensed arcs produced in matched simulated and observed cluster samples. The observed sample consists of 10 X-ray selected clusters at redshifts z ~ 0.2 imaged with HST by Smith et al. The simulated dataset is produced by lensing the Hubble Deep Field, which serves as a background source image, with 150 realizations (different projections and shifts) of five simulated z = 0.2 clusters from a LambdaCDM N-body simulation. The real and simulated clusters have similar masses, the real photometric redshift is used for each background source, and all the observational effects influencing arc detection in the real dataset, including light from cluster galaxies, are simulated in the artificial dataset. We develop, and apply to both datasets, an objective automatic arc-finding algorithm. We find consistent arc statistics in the real and in the simulated sample, with an average of ~ 1 detected giant (length to width ratio >= 10) arc per cluster and ~ 0.2 giant luminous (R<22.3 mag) arc per cluster. Thus, taking into account a realistic source population and observational effects, the clusters predicted by LambdaCDM have the same arc-production efficiency as the observed clusters. If, as suggested by other studies, there is a discrepancy between the predicted and the observed total number of arcs on the sky, it must be the result of differences between the redshift dependent cluster mass functions, and not due to differences in the lensing efficiency of the most massive clusters.

We use high-resolution N-body simulations to study the galaxy-cluster cross-sections and the abundance of giant arcs in the $\Lambda$CDM model. Clusters are selected from the simulations using the friends-of-friends method, and their cross-sections for forming giant arcs are analyzed. The background sources are assumed to follow a uniform ellipticity distribution from 0 to 0.5 and to have an area identical to a circular source with diameter $1\arcsec$. We find that the optical depth scales as the source redshift approximately as $\tau_{1''} = 2.25 \times 10^{-6}/[1+(\zs/3.14)^{-3.42}]$ ($0.6<\zs<7$). The amplitude is about 50% higher for an effective source diameter of $0.5\arcsec$. The optimal lens redshift for giant arcs with the length-to-width ratio ($L/W$) larger than 10 increases from 0.3 for $\zs=1$, to 0.5 for $\zs=2$, and to 0.7-0.8 for $\zs>3$. The optical depth is sensitive to the source redshift, in qualitative agreement with Wambsganss et al. (2004). However, our overall optical depth appears to be only $\sim$ 10% to 70% of those from previous studies. The differences can be mostly explained by different power spectrum normalizations ($\sigma_8$) used and different ways of determining the $L/W$ ratio. Finite source size and ellipticity have modest effects on the optical depth. We also found that the number of highly magnified (with magnification $|\mu|>10$) and ``undistorted'' images (with $L/W<3$) is comparable to the number of giant arcs with $|\mu|>10$ and $L/W>10$. We conclude that our predicted rate of giant arcs may be lower than the observed rate, although the precise `discrepancy' is still unclear due to uncertainties both in theory and observations.

We argue that a few per cent of "Early Dark Energy" can be detected by the statistics of nonlinear structures. The presence of Dark Energy during linear structure formation is natural in models where the present tiny Dark-Energy density is related to the age of the Universe rather than a new fundamental small parameter. Generalisation of the spherical collapse model shows that the linear collapse parameter delta_c is lowered. The corresponding relative enhancement of weak gravitational lensing on arc-minute scales lowers the value of sigma_8 inferred from a given lensing amplitude as compared to Lambda-CDM. In presence of Early Dark Energy, structures grow slower, such that for given sigma_8 the number of galaxies and galaxy clusters is substantially increased at moderate and high redshift. For realistic models, the number of clusters detectable through their thermal Sunyaev-Zel'dovich effect at redshift unity and above, e.g. with the Planck satellite, can be an order of magnitude larger than for Lambda-CDM.

We present an 8 sigma detection of cosmic magnification measured by the variation of quasar density due to gravitational lensing by foreground large scale structure. To make this measurement we used 3800 square degrees of photometric observations from the Sloan Digital Sky Survey (SDSS) containing \~200,000 quasars and 13 million galaxies. Our measurement of the galaxy-quasar cross-correlation function exhibits the amplitude, angular dependence and change in sign as a function of the slope of the observed quasar number counts that is expected from magnification bias due to weak gravitational lensing. We show that observational uncertainties (stellar contamination, Galactic dust extinction, seeing variations and errors in the photometric redshifts) are well controlled and do not significantly affect the lensing signal. By weighting the quasars with the number count slope, we combine the cross-correlation of quasars for our full magnitude range and detect the lensing signal at >4 sigma in all five SDSS filters. Our measurements of cosmic magnification probe scales ranging from 60 kpc/h to 10 Mpc/h and are in good agreement with theoretical predictions based on the WMAP concordance cosmology. As with galaxy-galaxy lensing, future measurements of cosmic magnification will provide useful constraints on the galaxy-mass power spectrum.

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.

We construct a linear filter optimised for detecting dark-matter halos in weak-lensing data. The filter assumes a mean radial profile of the halo shear pattern and modifies that shape by the noise power spectrum. Aiming at separating dark-matter halos from spurious peaks caused by large-scale structure lensing, we model the noise as being composed of weak lensing by large-scale structures and Poisson noise from random galaxy positions and intrinsic ellipticities. Optimal filtering against the noise requires the optimal filter scale to be smaller than typical halo sizes. Although a perfect separation of halos from spurious large-scale structure peaks is strictly impossible, we use numerical simulations to demonstrate that our filter produces substantially more sensitive, reliable and stable results than the conventionally used aperture-mass statistic.

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.

We adapt a non-linear filter proposed by Hu 2001 for detecting lensing of the CMB by large-scale structures to recover surface-density profiles of galaxy clusters from their localised, weak gravitational lensing effect on CMB fields. Shifting the band-pass of the filter to smaller scales, and normalising it such as to reproduce the convergence rather than the deflection angle, we find that the mean density profile of a sample of 100 clusters can be recovered to better than 10% from well within the scale radius to almost the virial radius. The kinetic Sunyaev-Zel'dovich effect is shown to be a negligible source of error. We test the filter applying it to data simulated using the characteristics of the Atacama Cosmology Telescope (ACT), showing that it will be possible to recover mean cluster profiles outside a radius of 1' corresponding to ACT's angular resolution.

We investigate statistical distributions of differences in gravitational-lensing deflections between two light rays, the so-called lensing excursion angles. A probability distribution function of the lensing excursion angles, which plays a key role in estimates of lensing effects on angular clustering of objects (such as galaxies, QSOs and also the CMB temperature map), is known to consist of two components; a Gaussian core and an exponential tail. We use numerical gravitational-lensing experiments in a Lambda-CDM cosmology for quantifying these two components. We especially focus on the physical processes responsible for generating those two components. We develop a simple empirical model for the exponential tail which allows us to explore its origin. We find that the tail is generated by the coherent lensing scatter by massive halos with M>10^14h^-1Msun at z<1 and that its exponential shape arises due to the exponential cut-off of the halo mass function at that mass range. On scales larger than 1 arc minute, the tail does not have a practical influence on the lensing effects on the angular clustering. Our model predicts that the coherent scatter may have non-negligible effects on angular clustering at sub-arcminute scales.

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.

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.

We study the efficiency of numerically simulated galaxy clusters for producing strong gravitational lensing events in different dark-energy cosmologies with constant and time-variable equation of state and we compare it with ``standard'' $\Lambda$CDM and OCDM models. Our main results are that: (1) the expected abundance of gravitational arcs with large length-to-width ratio depends on the equation of state of dark energy at the epoch of formation of the halo; (2) as found in previous studies, the strong lensing cross section of galaxy clusters is very sensitive to dynamical processes like mergers, whose impact however is different for different cosmologies, being stronger for models in which halos are less concentrated. As expected, the largest differences in the lensing optical depth occur at moderate and high redshift.

High-resolution, numerical simulations of 17 cluster-sized dark-matter haloes in eight different cosmologies with and without dynamical dark energy confirm the picture that core halo densities are imprinted early during their formation by the mean cosmological density. Quite independent of cosmology, halo concentrations have a log-normal distribution with a scatter of ~0.2 about the mean. We propose a simple scaling relation for halo concentrations in dark-energy cosmologies.

Compact groups (CGs) of galaxies -- relatively poor groups of galaxies in which the typical separations between members is of the order of a galaxy diameter -- offer an exceptional laboratory for the study of dense galaxian environments with short (<1Gyr) dynamical time-scales. In this paper, we present an objectively defined catalog of CGs in 153 sq deg of the Sloan Digital Sky Survey Early Data Release (SDSS EDR). To identify CGs, we applied a modified version of Hickson's (1982) criteria aimed at finding the highest density CGs and thus reducing the number of chance alignments. Our catalog contains 175 CGs down to a limiting galaxy magnitude of r* = 21. The resulting catalog has a median depth of approximately z = 0.13, substantially deeper than previous CG catalogs. Since the SDSS will eventually image up to one quarter of the celestial sphere, we expect our final catalog, based upon the completed SDSS, will contain on the order of 5,000 - 10,000 CGs. This catalog will be useful for conducting studies of the general characteristics of CGs, their environments, and their component galaxies.

It was recently claimed that several galaxy clusters containing radial and tangential gravitational arcs and having a measured velocity-dispersion profile for the brightest cluster galaxy had to have central density profiles considerably flatter than those found in CDM cluster simulations. Using a simple analytic mass model, we confirm this result_for axially symmetric_ mass distributions, but show that steep density profiles are well in agreement with the cluster requiring the flattest axially symmetric profile once even small deviations from axial symmetry are introduced.

We study the impact of merger events on the strong lensing properties of galaxy clusters. Previous lensing simulations were not able to resolve dynamical time scales of cluster lenses, which arise on time scales which are of order a Gyr. In this case study, we first describe qualitatively with an analytic model how some of the lensing properties of clusters are expected to change during merging events. We then analyse a numerically simulated lens model for the variation in its efficiency for producing both tangential and radial arcs while a massive substructure falls onto the main cluster body. We find that: (1) during the merger, the shape of the critical lines and caustics changes substantially; (2) the lensing cross sections for long and thin arcs can grow by one order of magnitude and reach their maxima when the extent of the critical curves is largest; (3) the cross section for radial arcs also grows, but the cluster can efficiently produce this kind of arcs only while the merging substructure crosses the main cluster centre; (4) while the arc cross sections pass through their maxima as the merger proceeds, the cluster's X-ray emission increases by a factor of $\sim5$. Thus, we conclude that accounting for these dynamical processes is very important for arc statistics studies. In particular, they may provide a possible explanation for the arc statistics problem.

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}).

We perform ray-tracing simulations evaluating the effect of a cD galaxy on the strong lensing properties of five galaxy cluster halos obtained from N-body simulations. The cD galaxy is modelled using both axially symmetric and elliptical models and assuming several masses for its dark matter halo. The effect of the cD orientation with respect to the mass distribution of the host galaxy cluster is also investigated. The simulations are carried out in an open and a flat model universe with cosmological constant. We find that the enhancement of the cluster lensing cross sections for long and thin arcs due to the presence of a massive cD at the cluster centre is typically less than 100%, depending on the model used for the cD galaxy and its orientation. The impact of the cD on the cluster efficiency for producing radially magnified images is larger only for those clusters whose lensing cross section for radial arcs is very small in absence of the central galaxy. We conclude that the presence of a cD galaxy at the cluster centre can only moderately influence the cluster efficiency for strong lensing and in particular fails to explain the discrepancy between the observed number of giant arcs in galaxy clusters and their abundance predicted from lensing simulations in the currently most favoured $\Lambda$CDM model.

We investigate how the probability of the formation of giant arcs in galaxy clusters is expected to change in cosmological models dominated by dark energy with an equation of state p=w rho c^2 compared to cosmological-constant or open models. To do so, we use a simple analytic model for arc cross sections based on the Navarro-Frenk-White density profile which we demonstrate reproduces essential features of numerically determined arc cross sections. Since analytic lens models are known to be inadequate for accurate absolute quantifications of arc probabilities, we use them only for studying changes relative to cosmological-constant models. Our main results are (1) the order of magnitude difference between the arc probabilities in low density, spatially flat and open CDM models found numerically is reproduced by our analytic model, and (2) dark-energy cosmologies with w>-1 increase the arc optical depth by at most a factor of two and are thus unlikely to reconcile arc statistics with spatially flat cosmological models with low matter density.

Arc statistics is known to be a powerful cosmological tool. Numerical lensing simulations show that orders of magnitude differences in the number of {\em giant} arcs on the whole sky are expected in different cosmological models. In this paper, we discuss the analytic and numerical methods in arc statistics and show that analytic models fail to reproduce the efficiency for strong lensing of more realistic numerical cluster models. Then, we discuss two recent extensions of the lensing simulations, i.e. the effects of cD galaxies in the lensing clusters and the impact of cluster mergers on arc statistics. We show that cD galaxies can increase the lensing cross sections for long and thin arcs by perhaps up to $\sim50%$, while major mergers can change the cluster efficiency for producing such arcs by up to one order of magnitude.

Most problems in gravitational lensing require numerical solutions. The most frequent types of problems are (1) finding multiple images of a single source and classifying the images according to their properties like magnification or distortion; (2) propagating light rays through large cosmological simulations; and (3) reconstructing mass distributions from their tidal field. This lecture describes methods for solving such problems. Emphasis is put on using adaptive-grid methods for finding images, issues of spatial resolution and reliability of statistics for weak lensing by large-scale structures, and methodical questions related to shear-inversion techniques.

The statistics of gravitationally lensed arcs was recognised earlier as a potentially powerful cosmological probe. However, while fully numerical models find orders of magnitude difference between the arc probabilities in different cosmological models, analytic models tend to find markedly different results. We introduce in this paper an analytic cluster lens model which improves upon existing analytic models in four ways. (1) We use the more realistic Navarro-Frenk-White profile instead of singular isothermal spheres, (2) we include the effect of cosmology on the compactness of the lenses, (3) we use elliptical instead of axially symmetric lenses, and (4) we take the intrinsic ellipticity of sources into account. While these improvements to the analytic model lead to a pronounced increase of the arc probability, comparisons with numerical models of the same virial mass demonstrate that the analytic models still fall short by a substantial margin of reproducing the results obtained with numerical models. Using multipole expansions of cluster mass distributions, we show that the remaining discrepancy can be attributed to substructure inside clusters and tidal fields contributed by the cluster surroundings, effects that cannot reasonably and reliably be mimicked in analytic models.

Intracluster magnetic fields with $ \sim \mu 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.

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.

Ongoing and planned wide-area surveys at optical and infrared wavelengths should detect a few times 10^5 galaxy clusters, roughly 10% of which are expected to be at redshifts >~ 0.8. We investigate what can be learned about the X-ray emission of these clusters from the ROSAT All-Sky Survey. While individual clusters at redshifts >~ 0.5 contribute at most a few photons to the survey, a significant measurement of the mean flux of cluster subsamples can be obtained by stacking cluster fields. We show that the mean X-ray luminosity of clusters with mass M >~ 2 x 10^{14} M_sun/h selected from the Sloan Digital Sky Survey should be measurable out to redshift unity with signal-to-noise >~ 10, even if clusters are binned with Delta z=0.1 and Delta ln M ~ 0.3. For such bins, suitably chosen hardness ratio allows the mean temperature of clusters to be determined out to z ~ 0.7 with a relative accuracy of Delta T/T <~ 0.15 for M > 10^{14} M_sun/h.

Via the magnification bias, gravitational lensing by large-scale structures causes angular cross-correlations between distant quasars and foreground galaxies on angular scales of arc minutes and above. We investigate the three-point cross-correlation between quasars and galaxy pairs measurable via the second moment of the galaxy counts around quasars and show that it reaches the level of a few per cent on angular scales near one arc minute. Combining two- and three-point correlations, a skewness parameter can be defined which is shown to be virtually independent on the shape and normalisation of the dark-matter power spectrum. If the galaxy bias is linear and deterministic, the skewness depends on the cosmic matter density parameter Omega_0 only; otherwise, it can be used to probe the linearity and stochasticity of the bias. We finally estimate the signal-to-noise ratio of a skewness determination and find that around twenty thousand distant quasars e.g. from the Sloan Digital Sky Survey should suffice for a direct measurement of Omega_0.

I present an overview of strong and weak gravitational lensing by galaxy clusters. After briefly introducing the principles of gravitational lensing, I discuss the main lessons learned from lensing on the mass distribution in clusters and their relation to cosmology.

We study the effects of a dark energy component with equation of state p=w rho with constant w >= -1 on the formation of Cold Dark Matter (CDM) haloes. We find two main effects: first, haloes form earlier as w increases, and second, the amplitude of the dark-matter power spectrum gets reduced in order to remain compatible with the large scale Cosmic Microwave Background (CMB) anisotropies. These effects counteract. Using recipes derived from numerical simulations, we show that haloes are expected to be up to ~ 50% more concentrated in CDM models with quintessence compared to Lambda-CDM models, the maximum increase being reached for w ~ -0.6. For larger w, the amplitude of the power spectrum decreases rapidly and makes expected halo concentrations drop. Halo detections through weak gravitational lensing are highly sensitive to halo concentrations. We show that weak-lensing halo counts with the aperture-mass technique can increase by a factor of ~ 2 as w is increased from -1 to -0.6, offering a new method for constraining the nature of dark energy.

We study the clustering properties of the galaxy clusters detectable for the Planck satellite due to their thermal Sunyaev-Zel'dovich effect. We take the past light-cone effect and the redshift evolution of both the underlying dark matter correlation function and the cluster bias factor into account. A theoretical mass-temperature relation allows us to convert the sensitivity limit of a catalogue into a minimum mass for the dark matter haloes hosting the clusters. We confirm that the correlation length is an increasing function of the sensitivity limits defining the survey. Using the expected characteristics of the Planck cluster catalogue, which will be a quite large and unbiased sample, we predict the two-point correlation function and power spectrum for different cosmological models. We show that the wide redshift distribution of the Planck survey, will allow to constrain the cluster clustering properties up to z=1. The dependence of our results on the main cosmological parameters (the matter density parameter, the cosmological constant and the normalisation of the density power-spectrum) is extensively discussed. We find that the future Planck clustering data place only mild constraints on the cosmological parameters, because the results depend on the physical characteristics of the intracluster medium, like the baryon fraction and the mass-temperature relation. Once the cosmological model and the Hubble constant are determined, the clustering data will allow a determination of the baryon fraction with an accuracy of few per cent.

The magnification bias of large-scale structures, combined with galaxy biasing, leads to a cross-correlation of distant quasars with foreground galaxies on angular scales of the order of arc minutes and larger. The amplitude and angular shape of the cross-correlation function w_QG contain information on cosmological parameters and the galaxy bias factor. While the existence of this cross-correlation has firmly been established, existing data did not allow an accurate measurement of w_QG yet, but wide area surveys like the Sloan Digital Sky Survey now provide an ideal database for measuring it. However, w_QG depends on several cosmological parameters and the galaxy bias factor. We study in detail the sensitivity of w_QG to these parameters and develop a strategy for using the data. We show that the parameter space can be reduced to the bias factor b, Omega_0 and sigma_8, and compute the accuracy with which these parameters can be deduced from SDSS data. Under reasonable assumptions, it should be possible to reach relative accuracies of the order of 5%-15% for b, Omega_0, and sigma_8. This method is complementary to other weak-lensing analyses based on cosmic shear.

We use cosmological magneto-hydrodynamic simulations to study the evolution of magnetic fields in galaxy clusters in two different cosmological models, a standard-CDM and a $\Lambda$-CDM model. We show that the magnetic field strength profiles closely follow the cluster density profiles outside a core region of radius $\sim200\kpc$. The magnetic field has a correlation length of order $50\kpc$ and reverses on scales of $\sim100\kpc$ along typical lines-of-sight. The power spectrum of the magnetic field can well be approximated by a steep power law with an exponent of $\sim-2.7$. The mean magnetic field in the cluster cores grows roughly exponentially with decreasing redshift, $B\sim10^{-2.5z} \muG$. Merger events have a pronounced effect on magnetic field evolution, which is strongly reflected in measurable quantities like the Faraday rotation. The field evolution in the two different cosmologies proceeds virtually identically. All our cluster models very well reproduce observed Faraday rotation measurements when starting with nG seed fields.

We show that the predictions of the recent physical model by Granato et al. (2001) are fully consistent with the available estimates of the clustering properties of galaxies at (sub)-mm wavelengths. The model is then exploited to compute the expected power spectrum of fluctuations due to clustering at the frequencies of the High Frequency Instrument (HFI) of ESA's {\sc Planck} satellite. The clustering signal is found to be detectable in regions of low interstellar dust emission. A further distinctive prediction of the adopted model is a remarkably high fraction of gravitationally lensed sources at bright mm/sub-mm fluxes. In fact, according to the adopted model, only thanks to gravitational lensing, a significant number of high-$z$ forming spheroidal galaxies will be detectable by large area, shallow surveys at mm/sub-mm wavelengths, such as those carried out by P{\sc lanck}/HFI. Allowing for other source populations, we find that the fraction of gravitationally lensed sources at bright mm/sub-mm fluxes is expected to be up to $\simeq 40%$.

Galaxy clusters can be detected as surface brightness enhancements in smoothed optical surveys. This method does not require individual galaxies to be identifiable, and enables clusters to be detected out to surprisingly high redshifts, as recently demonstrated by the Las Campanas Distant Cluster Survey (LCDCS). Here, we investigate redshift limits for cluster detection in the Sloan Digital Sky Survey (SDSS). Calibrating assumptions about the surface brightness profile, the mass-to-light ratio, and the spectral energy distribution of galaxy clusters using available observational data, we show that it should be possible to detect galaxy groups out to redshifts of ~0.5, and massive galaxy clusters out to redshifts of ~1.2 in summed r'+i'+z' SDSS data. Redshift estimates can be derived from the SDSS magnitudes of brightest cluster members out to redshifts near unity. Over the area of sky it covers, SDSS should find >~98% of the clusters detectable by the Planck satellite through the thermal Sunyaev-Zel'dovich effect. The few Planck clusters not detected in SDSS will almost all be at z>~1.2.

Cosmological simulations of magnetic fields in galaxy clusters show that remarkable agreement between simulations and observations of Faraday rotation and radio haloes can be achieved assuming that seed fields of 10^{-9} G were present at redshifts 15-20. The structure of the seed field is irrelevant for the final intracluster field. On average, the field grows exponentially with decreasing cluster redshift, but merger events cause steep transient increases in the field strength. Typical field-reversal scales are of order 50 kpc/h. In most cases, the intracluster fields are dynamically unimportant. Assuming secondary electron models, the average structure of cluster radio haloes can naturally be reproduced.

Upcoming wide-area surveys in the submillimetre regime will allow the construction of complete galaxy cluster samples through their thermal Sunyaev-Zel'dovich effect. We propose an analytic method to predict the number of gravitationally lensed giant arcs produced by the cluster sample expected to be detectable for the Planck satellite. The statistics of lensed arcs can be used to constrain cosmological parameters. This computation implies the choice of an analytic cluster model which has to reproduce the lensing properties of realistic galaxy clusters. Comparing the results of analytic computations and numerical ray-tracing simulations, we show that spherical models with the Navarro-Frenk-White density profile fail by far to produce as many giant arcs as realistic galaxy clusters. This confirms the great importance of asymmetries and substructures in the lensing matter distribution. By elliptically distorting the lensing potential of spherical NFW density profiles, we show that the discrepancy between analytic and numerical computations can almost completely be removed.

Integral measures of weak gravitational lensing by dark-matter haloes, like the aperture mass, are sensitive to different physical halo properties dependent on the halo mass density profile. For isothermal profiles, the relation between aperture mass and virial mass is steeper than for haloes with the universal NFW profile. Consequently, the halo mass range probed by the aperture mass is much wider for NFW than for isothermal haloes. We use recent modifications to the Press-Schechter mass function in CDM models normalised to the local abundance of rich clusters, to predict the properties of the halo sample expected to be accessible with the aperture mass technique. While ~10 haloes should be detected per square degree if the haloes have NFW profiles, their number density is lower by approximately an order of magnitude if they have isothermal profiles. These results depend only very mildly on the cosmological background model. We conclude that counts of haloes with a significant weak-lensing signal are a powerful discriminator between different dark-matter profiles.

We carry out ray tracing through five high resolution simulations of a galaxy cluster to study how its ability to produce giant gravitationally lensed arcs is influenced by the collision cross-section of its dark matter. In three cases typical dark matter particles in the cluster core undergo between 1 and 100 collisions per Hubble time; two more explore the long (``collisionless'') and short (``fluid'') mean free path limits. We study the size and shape distributions of arcs and compute the cross-section for producing ``extreme'' arcs of various sizes. Even a few collisions per particle modify the core structure enough to destroy the cluster's ability to produce long, thin arcs. For larger collision frequencies the cluster must be scaled up to unrealistically large masses before it regains the ability to produce giant arcs. None of our models with self-interacting dark matter (except the ``fluid'' limit) is able to produce radial arcs; even the case with the smallest scattering cross-section must be scaled to the upper limit of observed cluster masses before it produces radial arcs. Apparently the elastic collision cross-section of dark matter in clusters must be very small, below 0.1 cm^2 g^-1, to be compatible with the observed ability of clusters to produce both radial arcs and giant arcs.

We present detailed calculations of the magnification distribution, including both weak and strong lensing, using very recent solutions of the Dyer-Roeder (1973) equation for light propagation in a inhomogeneous universe with a cosmological constant and up-to-date models for the evolving cosmological distribution of dark matter halos. We address in particular the effect of a non-zero cosmological constant, of different density profiles of lenses, and of finite sizes of lensed sources. We show that, if dark matter fluctuations are normalized to the local cluster abundance, in the presence of a cosmological constant the optical depth for lensing {\em decreases} compared to the case of an Einstein-de Sitter universe, because halos in the relevant mass range are less abundant over a large redshift interval. We also discuss the differences in the magnification probability distributions produced by Navarro, Frenk & White (NFW) and by Singular Isothermal Sphere (SIS) density profiles of lenses. We find that the NFW lens is more efficient for moderate magnifications ($2\lsim A\lsim 4$), and less efficient for larger magnifications. Moreover, we discuss quantitatively the maximum magnification, $A_{\rm max}$, that can be achieved in the case of extended sources (galaxies) with realistic luminosity profiles, taking into account the possible ellipticity of the lens potential. Finally, we apply our results to a class of sources following the luminosity evolution typical for a unified scheme of QSO formation. We find that the lensed source counts at $850 \mu$m can be larger than the unlensed ones by several orders of magnitude at flux densities $\gta100 $mJy.

We use cosmological, magneto-hydrodynamic simulations of galaxy clusters to quantify the dynamical importance of magnetic fields in these clusters. The set-up of initial magnetic field strengths at high redshifts is chosen such that observed Faraday-rotation measurements in low-redshift clusters are well reproduced in the simulations. We compute the radial profiles of the intracluster gas temperature and of the thermal and magnetic pressure in a set of clusters simulated in the framework of an Einstein-de Sitter and a low-density, spatially-flat CDM cosmological model. We find that, for a realistic range of initial magnetic field strengths, the temperature of the intracluster gas changes by less than $\approx5%$.

We present a new method to determine the total mass-to-neutral gas ratio in Lyman-limits systems. The method exploits the relation between the neutral hydrogen column density and the magnification of background sources due to the weak gravitational lensing that these systems induce. Because weak lensing does not provide a direct measure of mass, one must use this relation in a statistical sense to solve for the average mass-to-gas ratio and its distribution. We use a detailed mock catalog of quasars (sources) and Lyman-limit systems (lenses) to demonstrate the applicability of this approach through our ability to recover the parameter. This mock catalog also allows us to check for systematics in the method and to sketch its limitations. For a universal constant mass-to-gas ratio and a sample of N quasars, we obtain an unbiased estimate of its value with 95% confidence limits (independent of its actual value) of +/- 140 {10^5/N)^0.5.

Full-sky microwave surveys like the upcoming Planck satellite mission will detect of order 10^4 galaxy clusters through their thermal Sunyaev-Zel'dovich effect. I investigate the properties of the gravitationally lensing subsample of these clusters. The main results are: (1) The combined sample comprises >~70% of the complete sample. (2) It is confined to redshifts 0.2+-0.1, and to masses (5+-3) x 10^14 solar masses. (3) Using a particular measure for the weak lensing effect, viz. the aperture mass, cluster masses can be determined with a relative accuracy of ~20% if their density profile is known. Consequently, the mass function of the combined sample can accurately be measured. (4) For low-density universes, I predict a sharp peak in the measured (aperture) mass function near 5 x 10^14 solar masses and explain its origin, showing that the peak will be absent in high-density universes. (5) The location of the peak and the exponential decrease of the mass function on its high-mass side will allow the determination of the amplitude of the dark-matter power spectrum on the cluster scale and the baryon fraction in clusters, and constrain the thermal history of the intracluster gas.

We present the results of a set of numerical simulations aiming at evaluating the effects of cluster galaxies on the arc statistics. At this goal we use nine different galaxy clusters obtained from N-body simulations. We mimic the presence of a population of galaxies inside each cluster, trying to reproduce the observed luminosity function and the spatial distribution. We compare the statistical distributions of the properties of the gravitational arcs produced by our clusters with and without galaxies. We find that the cluster galaxies do not introduce perturbations strong enough to significantly change the number of arcs and the distributions of the properties of long arcs.

We present the results of a set of numerical simulations evaluating the effect of cluster galaxies on arc statistics. We perform a first set of gravitational lensing simulations using three independent projections for each of nine different galaxy clusters obtained from N-body simulations. The simulated clusters consist of dark matter only. We add a population of galaxies to each cluster, mimicking the observed luminosity function and the spatial galaxy distribution, and repeat the lensing simulations including the effects of cluster galaxies, which themselves act as individual lenses. Each galaxy is represented by a spherical Navarro, Frenk & White (1997) density profile. We consider the statistical distributions of the properties of the gravitational arcs produced by our clusters with and without galaxies. We find that the cluster galaxies do not introduce perturbations strong enough to significantly change the number of arcs and the distributions of lengths, widths, curvature radii and length-to-width ratios of long arcs. We find some changes to the distribution of short-arc properties in presence of cluster galaxies. The differences appear in the distribution of curvature radii for arc lengths smaller than 12'', while the distributions of lengths, widths and length-to-width ratios are significantly changed only for arcs shorter than 4''.

We review theory and applications of weak gravitational lensing. After summarising Friedmann-Lemaitre cosmological models, we present the formalism of gravitational lensing and light propagation in arbitrary space-times. We discuss how weak-lensing effects can be measured. The formalism is then applied to reconstructions of galaxy-cluster mass distributions, gravitational lensing by large-scale matter distributions, QSO-galaxy correlations induced by weak lensing, lensing of galaxies by galaxies, and weak lensing of the cosmic microwave background. Contents: Introduction - Cosmological Background - Gravitational Light Deflection - Principles of Weak Gravitational Lensing - Weak Lensing by Galaxy Clusters - Weak Cosmological Lensing - QSO Magnification Bias and Large-Scale Structure - Galaxy-Galaxy Lensing - The Impact of Weak Gravitational Lensing on the Microwave Background Radiation - Summary and Outlook

We propose a new cluster deprojection algorithm for recovering cluster structure along the line of sight (los) based on the Richardson-Lucy (RL) algorithm, which recovers a nonnegative theoretical distribution from a given projection. To optimize our reconstruction we combine X-ray and Sunyaev-Zeldovich (SZ) data with gravitational lensing maps. The combination of these three data types constrains the structure of rich clusters of galaxies. To validate our new approach, we apply the algorithm to gas-dynamical cluster simulations. A reconstruction using the lensing potential is shown.

Given its extraordinary spatial resolution and sensitivity, the projected Next Generation Space Telescope (NGST) is likely to detect a large number of high-redshift QSOs lensed by spiral galaxies. Using realistic models for the QSO and spiral populations, we calculate the expected number density of detectable QSOs multiply imaged by spirals, and investigate the influence of various evolution effects on that number density. It is shown that NGST will probably find of order ten lensed QSOs per square degree at 26th magnitude in the V and L bands, and that various observable quantities like the total number density of lensed QSOs in these two bands, the ratio between the number densities of lensed QSOs in the V and L bands, the fraction of QSOs with more than two images and so forth can be used to constrain the evolution of the QSOs, the spirals, the dust in spirals, and the masses of spiral disks.

We perform cosmological, hydrodynamic simulations of magnetic fields in galaxy clusters. The computational code combines the special-purpose hardware Grape for calculating gravitational interaction, and smooth-particle hydrodynamics for the gas component. We employ the usual MHD equations for the evolution of the magnetic field in an ideally conducting plasma. As a first application, we focus on the question what kind of initial magnetic fields yield final field configurations within clusters which are compatible with Faraday-rotation measurements. Our main results can be summarised as follows: (i) Initial magnetic field strengths are amplified by approximately three orders of magnitude in cluster cores, one order of magnitude above the expectation from spherical collapse. (ii) Vastly different initial field configurations (homogeneous or chaotic) yield results that cannot significantly be distinguished. (iii) Micro-Gauss fields and Faraday-rotation observations are well reproduced in our simulations starting from initial magnetic fields of \~10^-9 G strength. Our results show that (i) shear flows in clusters are crucial for amplifying magnetic fields beyond simple compression, (ii) final field configurations in clusters are dominated by the cluster collapse rather than by the initial configuration, and (iii) initial magnetic fields of order 10^-9 G are required to match Faraday-rotation observations in real clusters.

We perform cosmological, hydrodynamic simulations of magnetic fields in galaxy clusters. The computational code combines the special-purpose hardware Grape for calculating gravitational interaction, and smooth-particle hydrodynamics for the gas component. We employ the usual MHD equations for the evolution of the magnetic field in an ideally conducting plasma. As a first application, we focus on the question what kind of initial magnetic fields yield final field configurations within clusters which are compatible with Faraday-rotation measurements. Our main results can be summarised as follows: (i) Initial magnetic field strengths are amplified by approximately three orders of magnitude in cluster cores, one order of magnitude above the expectation from spherical collapse. (ii) Vastly different initial field configurations (homogeneous or chaotic) yield results that cannot significantly be distinguished. (iii) Micro-Gauss fields and Faraday-rotation observations are well reproduced in our simulations starting from initial magnetic fields of \~ 10^-9 G strength at redshift 15. Our results show that (i) shear flows in clusters are crucial for amplifying magnetic fields beyond simple compression, (ii) final field configurations in clusters are dominated by the cluster collapse rather than by the initial configuration, and (iii) initial magnetic fields of order 10^-9 G are required to match Faraday-rotation observations in real clusters.

We examine gravitational lensing constraints on the structure of galaxy clusters and compare them with the results of cosmological N-body simulations of cluster formation in cold dark matter (CDM) dominated universes. We find that cluster core masses, as measured by the observed location of giant tangential arcs, generally exceed those of dark matter halos of similar velocity dispersion. The magnitude of the discrepancy is a strong function of cluster mass. Arc properties in the most massive clusters in the sample (i.e. those with velocity dispersion, \sigma \sim 1500-2000 km/s) are essentially consistent with the N-body predictions. On the other hand, giant arcs in \sigma \sim 1000 km/s clusters can only be reconciled with CDM cluster halos if their lensing power has been increased substantially by the presence of a massive (\sim 3 \times 10^{12} h^{-1} M_{\odot}) central galaxy and of significant substructure. Best agreement is found if the mass of the central galaxy and the effects of substructure are approximately independent of cluster mass. Massive central galaxies with steep inner density profiles are also needed to explain a clear trend observed in our dataset between the radial thickness of giant tangential arcs and the velocity dispersion of the cluster lens. The position and redshift of radial arcs may be used as independent tests of these results, but at present the dataset available is too limited to have a significant impact on these conclusions. Our results depend only weakly on the cosmological model adopted, and suggest that structural parameters of clusters derived from strong lensing studies cannot usefully constrain the values of cosmological parameters.

We demonstrate that the aperture mass as a measure for cosmic shear closely approximates (to better than ~5%) the scaled and shifted power spectrum of the projected mass density. This cosmological weak-lensing information can thus be used to directly infer the projected matter power spectrum with high accuracy. As an application, we show that aperture-mass observations can be used to constrain the cosmic density parameter and the power-spectrum amplitude. We show that, for a particular example, it should be possible to constrain Omega_0 to within ~ +/- 27%, and sigma_8 to within ~ +/- 8% using weak-shear data on a square-shaped field of 8 degrees side length.

Selection of galaxy clusters by mass is now possible due to weak gravitational lensing effects. It is an important question then whether this type of selection reduces the projection effects prevalent in optically selected cluster samples. We address this question using simulated data, from which we construct synthetic cluster catalogues both with Abell's criterion and an aperture-mass estimator sensitive to gravitational tidal effects. The signal-to-noise ratio of the latter allows to some degree to control the properties of the cluster sample. For the first time, we apply the cluster-detection algorithm proposed by Schneider to large-scale structure simulations. We find that selection of clusters through weak gravitational lensing is more reliable in terms of completeness and spurious detections. Choosing the signal-to-noise threshold appropriately, the completeness can be increased up to 100%, and the fraction of spurious detections can significantly be reduced compared to Abell-selected cluster samples. We also investigate the accuracy of mass estimates in cluster samples selected by both luminosity and weak-lensing effects. We find that mass estimates from gravitational lensing, for which we employ the zeta-statistics by Kaiser et al., are significantly more accurate than those obtained from galaxy kinematics via the virial theorem.

We present a new method for reconstructing two-dimensional mass maps of galaxy clusters from the image distortion of background galaxies. In contrast to most previous approaches, which directly convert locally averaged image ellipticities to mass maps (direct methods), our entropy-regularized maximum-likelihood method is an inverse approach. Albeit somewhat more expensive computationally, our method allows high spatial resolution in those parts of the cluster where the lensing signal is strong enough. Furthermore, it allows to straightforwardly incorporate additional constraints, such as magnification information or strong-lensing features. Using synthetic data, we compare our new approach to direct methods and find indeed a substantial improvement especially in the reconstruction of mass peaks. The main differences to previously published inverse methods are discussed.

Gravitational lensing of a quasar by a spiral galaxy should often be accompanied by damped Lyman-alpha absorption and dust extinction due to the intervening gaseous disk. In nearly edge-on configurations, the surface mass density of the gas and stars in the disk could by itself split the quasar image and contribute significantly to the overall lensing cross section. We calculate the lensing probability of a disk+halo mass model for spiral galaxies, including cosmic evolution of the lens parameters. A considerable fraction of the lens systems contains two images with sub-arcsecond separation, straddling a nearly edge-on disk. Because of that, extinction by dust together with observational selection effects (involving a minimum separation and a maximum flux ratio for the lensed images), suppress the detection efficiency of spiral lenses in optical wavebands by at least an order of magnitude. The missing lenses could be recovered in radio surveys. In modifying the statistics of damped Lyman-alpha absorbers, the effect of extinction dominates over the magnification bias due to lensing.

We propose a method to remove the mass sheet degeneracy that arises when the mass of galaxy clusters is inferred from gravitational shear. The method utilizes high-redshift standard candles that undergo weak lensing. Natural candidates for such standard candles are type Ia supernovae (SN Ia). When corrected with the light-curve shape (LCS), the peak magnitude of SN Ia provides a standard candle with an uncertainty in apparent magnitude of $\Delta m\simeq 0.1-0.2$. Gravitational magnification of a background SN Ia by an intervening cluster would cause a mismatch between the observed SN Ia peak magnitude compared to that expected from its LCS and redshift. The average detection rate for SN Ia with a significant mismatch of $\ge2\Delta m$ behind a cluster at $z\simeq0.05-0.15$ is about $1-2$ supernovae per cluster per year at $J,I,R\lesssim25-26$. Since SNe are point-like sources for a limited period, they can experience significant microlensing by MACHOs in the intracluster medium. Microlensing events caused by MACHOs of $\sim10^{-4} M_\odot$ are expected to have time scales similar to that of the SN light curve. Both the magnification curve by a MACHO and the light curve of a SN Ia have characteristic shapes that allow to separate them. Microlensing events due to MACHOs of smaller mass can unambiguously be identified in the SN light curve if the latter is continuously monitored. The average number of identifiable microlensing events per nearby cluster ($z\lesssim0.05$) per year is $\sim 0.02 (f/0.01)$, where $f$ is the fraction of the cluster mass in MACHOs of masses $10^{-7} < M_{macho}/M_\odot < 10^{-4}$.

These lectures give an introduction to Gravitational Lensing. We discuss lensing by point masses, lensing by galaxies, and lensing by clusters and larger-scale structures in the Universe. The relevant theory is developed and applications to astrophysical problems are discussed.

The gravitational lens effect of galaxy clusters can produce large arcs from source galaxies in their background. Typical source redshifts of ~ 1 require clusters at z ~ 0.3 for arcs to form efficiently. Given the cluster abundance at the present epoch, the fewer clusters exist at z ~ 0.3 the higher Omega_0 is, because the formation epoch of galaxy clusters strongly depends on Omega_0. In addition, at fixed Omega_0, clusters are less concentrated, and hence less efficient lenses, when the cosmological constant is positive, Omega_Lambda > 0. Numerical cluster simulations show that the expected number of arcs on the sky is indeed a sensitive function of Omega_0 and Omega_Lambda. The numerical results are compatible with the statistics of observed arcs only in a universe with low matter density, Omega_0 ~ 0.3, and zero cosmological constant. Other models fail by one or two orders of magnitude, rendering arc statistics a sensitive probe for cosmological parameters.

We use numerical simulations of galaxy clusters in different cosmologies to study their ability to form large arcs. The cosmological models are: Standard CDM (SCDM; Omega_0=1, Omega_Lambda=0); tauCDM with reduced small-scale power (parameters as SCDM, but with a smaller shape parameter of the power spectrum); open CDM (OCDM; Omega_0=0.3, Omega_Lambda=0); and spatially flat, low-density CDM (LambdaCDM; Omega_0=0.3, Omega_Lambda=0.7). All models are normalised to the local number density of rich clusters. Simulating gravitational lensing by these clusters, we compute optical depths for the formation of large arcs. For large arcs with length-to-width ratio >= 10, the optical depth is largest for OCDM. Relative to OCDM, the optical depth is lower by about an order of magnitude for LambdaCDM, and by about two orders of magnitude for S/tauCDM. These differences originate from the different epochs of cluster formation across the cosmological models, and from the non-linearity of the strong lensing effect. We conclude that only the OCDM model can reproduce the observed arc abundance well, while the other models fail to do so by orders of magnitude.

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.

The projected gravitational potential of galaxy clusters is reflected in both their X-ray emission and their imprint on the images of background sources due to their gravitational lensing effects. Since these projections of the potential are weighted differently along the line-of-sight, we propose a method to combine them and remove the degeneracy between two cases: (i) a cluster consisting of a single potential well, or (ii) an apparent cluster composed of several potential wells projected onto each other. We demonstrate with simulated data of potential models that this method indeed allows to significantly distinguish multiple from single clusters. The confidence limit for this distinction depends on the mass ratio between the clusters. It ranges from ~ 15 sigma for mass ratio 1:1 to ~ 4 sigma for mass ratio 1:6. Furthermore, the method reconstructs the correct cluster mass, the correct mass ratio of the two clusters, and the correct scale radii with typical fractional accuracies of a few percent at 3 sigma confidence. As an aside, our method allows to accurately determine gas fractions in clusters, also with 3 sigma fractional accuracies of order a few percent. We argue that our method provides an alternative to the commonly used beta-fit technique, and yields more reliable results in a broader range of cases.

Gravitational lensing of an optical QSO by a spiral galaxy is often counteracted by dust obscuration, since the line-of-sight to the QSO passes close to the center of the galactic disk. The dust in the lens is likely to be correlated with neutral hydrogen, which in turn should leave a Lyman-alpha absorption signature on the QSO spectrum. We use the estimated dust-to-gas ratio of the Milky-Way galaxy as a mean and allow a spread in its values to calculate the effects of dust on lensing by low redshift spiral galaxies. Using a no-evolution model for spirals at z<1 we find (in Lambda=0 cosmologies) that the magnification bias due to lensing is stronger than dust obscuration for QSO samples with a magnitude limit B<16. The density parameter of neutral hydrogen, Omega_HI, is overestimated in such samples and is underestimated for fainter QSOs.

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.

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

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 $\rho\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.

Precise determinations of the image positions in quad gravitational lenses using VLBI can be used to measure the transverse velocity of the lens galaxy and the observer. The typical proper motions are $\mu$as yr$^{-1}$, so the time scale to measure the motion is ten years. By measuring the dipole of the proper motions in an ensemble of lenses we can set limits on the deviation of the inertial frame defined by the lenses from that defined by the CMB dipole and estimate the Hubble constant. The residual proper motions after subtracting the dipole probe the evolution of peculiar velocities with redshift and can be used to estimate the density parameter $\Omega_0$. For $N$ lenses, VLBI measurement accuracies of $\sigma_\theta$, and a baseline of $T$ years, we estimate that the 2$\sigma$ limit on the rms peculiar velocity of the lens galaxies is $3100 (\sigma_\theta/10\mu\{as})({yrs}/T)/N^{1/2} \kms$, and that the time required for the 2--$\sigma$ limit to reach the level of the local rms peculiar velocity $v_{0,rms}$ is approximately $10 N^{-1/2} (v_{0,rms}/600\kms)(\sigma_\theta/10\mu as)$ years. For a ten year baseline and $N=10$ lenses we expect the 1$\sigma$ limit on the misalignment with the CMB dipole to be $\Delta \theta=20^{\circ}$ or equivalently to obtain an upper limit of $\Delta H_0 /H_0 < 0.34 $.

We develop a new statistical method to reanalyse angular correlations between background QSOs and foreground galaxies that are supposed to be a consequence of dark matter inhomogeneities acting as weak gravitational lenses. The method is based on a weighted average over the galaxy positions and is optimized to distinguish between a random distribution of galaxies around QSOs and a distribution which follows an assumed QSO-galaxy two-point correlation function, by choosing an appropriate weight function. With simulations we demonstrate that this weighted average is slightly more significant than Spearman's rank-order test which was used in previous investigations. In particular, the advantages of the weighted average show up if the two-point correlation function is weak. We then reanalyze the correlation between high-redshift 1-Jansky QSOs and IRAS galaxies, taken from the IRAS Faint Source Catalog; these samples were analyzed previously using Spearman's rank-order test. In agreement with the previous work, we find moderate to strong correlations between these two samples; considering the angular two-point correlation function of these samples, we find a typical scale of order $5'$ from which most of the correlation signal derives. However, the statistical significance of the correlation changes with the redshift slices of the QSO sample one considers. Comparing with simple theoretical estimates of the expected correlation, we find that the signal we derive is considerably stronger than expected. On the other hand, recent direct verifications of the overdensity of matter in the line-of-sight to high-redshift radio QSOs obtained from the shear field around these sources, indicates that the observed association can be attributed to a gravitational lens effect.

We present a novel method to recontruct the mass distribution of galaxy clusters from their gravitational lens effect on background galaxies. The method is based on a least-chisquare fit of the two-dimensional gravitational cluster potential. The method combines information from shear and magnification by the cluster lens and is designed to easily incorporate possible additional information. We describe the technique and demonstrate its feasibility with simulated data. Both the cluster morphology and the total cluster mass are well reproduced.

Damped Lyman-alpha absorbers are believed to be associated with galactic disks. We show that gravitational lensing can therefore affect the statistics of these systems. First, the magnification bias due to lensing raises faint QSOs above a given magnitude threshold and thereby enhances the probability for observing damped absorption systems. Second, the bending of light rays from the source effectively limits the minimum impact parameter of the line-of-sight relative to the center of the absorber, thus providing an upper cut-off to the observed neutral hydrogen (HI) column density. The combination of these effects yields a pronounced peak in the observed abundance of absorbers with high column densities (>2*10^{21} cm^{-2}) and low redshifts (z<1). The inferred value of the cosmological density parameter of neutral hydrogen, Omega_{HI}, increases with increasing redshift and luminosity of the sources even if the true HI density remains constant. This trend resembles the observed evolution of Omega_{HI}(z). Damped Lyman-alpha absorbers with column densities >10^{21} cm^{-2} and redshifts 0.5<z<1 are reliable flags for lensed QSOs with a close pair of images separated by 0.3 arcsec. Detection of these gravitational lensing signatures with the Hubble Space Telescope can be used to constrain the depth of the absorber potential-wells and the cosmological constant.

We propose a new technique, which we call the lens parallax method, to determine simultaneously the redshift distribution of the faint blue galaxies and the mass distributions of foreground clusters of galaxies. The method is based on gravitational lensing and makes use of the following: (1) the amplitude of lensing-induced distortions of background galaxies increases with redshift; (2) the surface brightnesses of galaxies decrease steeply with redshift. The distortions of galaxy images due to lensing are thus expected to be inversely correlated with surface brightness, allowing us to obtain relative distances to galaxies as a function of surface brightness. If the redshifts of the brightest galaxies are measured, then the relative distance scale can be converted to mean galaxy redshifts as a function of surface brightness. Further, by comparing the angular sizes of lensed galaxies with those of similar galaxies in empty control fields, it is possible to break the so-called mass sheet degeneracy inherent to cluster mass reconstruction techniques which are based purely on image ellipticities. This allows an unambiguous determination of the surface density of a lensing cluster. We describe an iterative algorithm based on these ideas and present numerical simulations which show that the proposed techniques are feasible with a sample of ~ 10 rich clusters at moderate redshifts ~ 0.3-0.4 and an equal number of control fields. The numerical tests show that the method can be used to determine the redshifts of galaxies with an accuracy of dz ~ 0.1-0.2 at z ~ 1-1.7, and to measure the masses of lensing clusters to about 5% accuracy.

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%.

We propose a new technique, which we call the lens parallax method, to determine simultaneously the redshift distribution of the faint blue galaxies and the mass distribution of foreground clusters of galaxies. The method is based on gravitational lensing and uses the fact that the mean angular size of the faint blue galaxies is a well-determined function of surface brightness. The method requires first a calibration of the angular-size vs. surface brightness relation in unlensed blank fields and a determination of the mean redshift of the brightest of the faint blue galaxies. By combining this information with observations of the distorted images of background galaxies in the fields of foreground clusters, we show that it is possible to obtain the mean redshift of the galaxies as a function of their surface brightness. With a sample of about ten moderate redshift rich clusters and using ten bins in surface brightness, a bin-to-bin signal-to-noise ratio of $\sim3.5$ can be achieved. The method simultaneously allows a determination of the convergence and shear of the lensing cluster as a function of position, and through this the mass distribution of the lens can be obtained. The lens parallax method can be used in conjunction with, and will improve the accuracy of, other previously proposed techniques for mapping clusters.

We review and discuss information on the mass distribution of galaxy clusters obtained from gravitational lensing.

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$.

We construct a sample of numerical models for clusters of galaxies and employ these to investigate their capability of imaging background sources into long arcs. Emphasis is laid on the statistics of these arcs. We study cross sections for arc length and length-to-width ratio and optical depths for these arc properties, we examine the distribution of arc widths and curvature radii among long arcs, and we compare these results to predictions based on simplified (radially symmetric) cluster models. We find that the capability of the numerically modeled clusters to produce long arcs is larger by about two orders of magnitude than that of spherically symmetric cluster models with the same observable parameters (core radii and velocity dispersions), and that they are similarly efficient as singular isothermal spheres with the same velocity-dispersion distribution. The influence of source ellipticity is also investigated; we find that the optical depth for arcs with a length-to-width ratio $\ga10$ is significantly larger for elliptical than for circular sources. Given these results, we conclude that spherically symmetric lens models for galaxy clusters, adapted to the observable parameters of these clusters, grossly underestimate the frequency of long arcs. We attribute this difference between numerically constructed and simplified analytical lens models to the abundance and the extent of intrinsic asymmetry and to substructure in galaxy clusters.

The angular two-point correlation function between background QSOs and foreground galaxies induced by gravitational lensing is derived. It is shown that the shape of this correlation function depends sensitively on the spectrum of the density fluctuations in the Universe, thus providing a possibility to distinguish between different models for the spectrum. Using numerical large-scale structure simulations, I estimate that the QSO-galaxy correlation function can be measured from galaxy counts down to $\simeq21^{\rm st}\ldots22^{\rm nd}$ magnitude in fields with radius $\la25'$ around $50\ldots100$ QSOs with redshift $z\ge1$. Since the QSO-galaxy correlation function is proportional to $(-a-1)b$, where $b$ is the biasing factor of galaxy formation and $a$ is the effective slope of the QSO number counts, steep number counts are favorable for this kind of analysis. I show that $a=-4\ldots-5$ can be achieved with the 1-Jansky sample of radio-loud QSOs when the double-waveband magnification bias is employed. Moreover, the cross-correlation analysis allows to determine the galaxy-formation bias factor $b$.

We announce the detection of correlations on angular scales of $\ga10'$ between optically bright, high-redshift, radio-loud QSOs with diffuse X-ray emission seen by ROSAT in the {\it All-Sky Survey}. These correlations reach significance levels of up to $99.8\%$. A comparison of the results with a sample of control fields, bootstrapping analyses, and Kolmogorov-Smirnov tests provide unambiguous support for the statistical significance of the correlations found. We argue that the detected enhanced diffuse X-ray emission is in the foreground of the QSOs, and that it is probably due to galaxy clusters which magnify the QSOs by their gravitational lensing effect, thereby giving rise to a magnification bias in the background source sample. A comparison of the results presented below with correlations previously found between the same QSO sample and either Lick or IRAS galaxies provides further evidence for this interpretation, and identifies positions in the sky where weak gravitational lensing may be detected by searching for coherent distortions of background galaxy images.

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 $\lambda\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 $\Upsilon$ from which the dark halo forms. Within this approximation, an expression for the spin parameter can be derived, which depends on the geometry of $\Upsilon$, the cosmological density parameter $\Omega_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 $\lambda\simeq0.07^{+0.04}_{-0.05}$, consistent with numerical investigations. This probability distribution steeply rises with increasing spin parameter, reaching its maximum at $\lambda\simeq0.025$. The 10 (50,90) percentile values are $\lambda=0.02$ (0.05,0.11, respectively). There is a weak anticorrelation of the spin parameter with the peak height $\nu$ of the density fluctuation field $\lambda\propto \nu^{-0.29}$. The dependence on $\Omega_0$ and the variance $\sigma$ 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 $\sigma$, our results should be quite generally applicable and independent on

{}From the lensing properties of a numerically simulated cluster, from the statistical properties of the giant arcs produced by the cluster, and from comparing the cluster properties derived both from lensing and from the X-ray properties of the intracluster gas, we conclude (1) that clusters may be significantly more efficient for lensing than estimated from their velocity dispersion on the basis of iso- thermal-sphere models, (2) that the discrepancy between cluster core radii derived from lensing and from their X-ray surface brightness may be understood assuming that the intracluster gas was expelled from the cluster galaxies when the dark `body' of the cluster had already formed, and (3) that clusters may be efficient in producing giant arcs even if they have extended cores and a fairly shallow sur- face-density profile. This may mean that arcs may be much easier to produce, and therefore more abundant, than previously estimated.

We demonstrate that high-redshift 1-Jansky radio sources are correlated with foreground IRAS galaxies on angular scales of several 10 arc minutes on a statistical significance level of up to 99.8\%. Based on a theoretical analysis published earlier, we propose to interpret these correlations in terms of gravitational lensing by large- scale structures.

We construct a simple elliptical gravitational lens model for the quadruple lens B1422+231 and show that the details of the configuration cannot be easily understood in terms of this model; in particular, the flux ratios of the images are hard to reproduce. This qualitatively verifies the results from a different lens model constructed for the same object by Hogg and Blandford.