Search

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

The dynamics of correlated systems is relevant in many fields ranging from cosmology to plasma physics. However, they are challenging to predict and understand even for classical systems due to the typically large numbers of particles involved. Here, we study the evolution of an ultracold, correlated many-body system with repulsive interactions and initial correlations set by the Rydberg blockade using the analytical framework of Kinetic Field Theory (KFT). The KFT formalism is based on the path-integral formulation for classical mechanics and was first developed and successfully used in cosmology to describe structure formation in Dark Matter. The theoretical framework offers a high flexibility regarding the initial configuration and interactions between particles and, in addition, is computationally cheap. More importantly, the analytic approach allows us to gain better insight into the processes which dominate the dynamics. In this work we show that KFT can be applied in a much more general context and study the evolution of a correlated ion plasma. We find good agreement between the analytical KFT results for the evolution of the correlation function and results obtained from numerical simulations. We use the correlation functions obtained with KFT to compute the temperature increase in the ionic system due to disorder-induced heating. For certain choices of parameters we observe that the effect can be reversed, leading to correlation cooling. Due to its numerical efficiency as compared to numerical simulations, a detailed study using KFT can help to constrain parameter spaces where disorder-induced heating is minimal in order to reach the regime of strong coupling.

We present recent improvements in the perturbative treatment of particle interactions in Kinetic Field Theory (KFT) for inertial Zel'dovich trajectories. KFT has been developed for the systematic analytical calculation of non-linear cosmic structure formation on the basis of microscopic phase-space dynamics. We improve upon the existing treatment of the interaction operator by deriving a more rigorous treatment of phase-space trajectories of particles in an expanding universe. We then show how these results can be applied to KFT perturbation theory by calculating corrections to the late-time dark matter power spectrum at second order in the interaction operator. We find that the modified treatment of interactions w.r.t. inertial Zel'dovich trajectories improves the agreement of KFT with simulation results on intermediate scales compared to earlier results. Additionally, we illustrate that including particle interactions up to second order leads to a systematic improvement of the non-linear power spectrum compared to the first-order result., 26 pages, 4 figures

We study the asymptotic behaviour of integrals of the Laplace-Fourier type $P(k) = \int_\Omega\mathrm{e}^{-|k|^sf(x)}\mathrm{e}^{\mathrm{i} kx}\mathrm{d} x\;, $ with $k\in\mathbb{R}^d$ in $d\ge1$ dimensions, with $\Omega\subset\mathbb{R}^d$ and sufficiently well-behaved functions $f:\Omega\to\mathbb{R}$. Our main result is $ P(k) \sim \frac{\mathrm{e}^{-|k|^sf(0)}}{|k|^{sd/2}}\sqrt{\frac{(2\pi)^d}{\det A}} \exp\left(-\frac{k^\top A^{-1}k}{2|k|^s}\right) $ for $|k|\to\infty$, where $A$ is the Hessian matrix of the function $f$ at its critical point, assumed to be at $x_0 = 0$. In one dimension, the Hessian is replaced by the second derivative, $A = f''(0)$. We also show that the integration domain $\Omega$ can be extended to $\mathbb{R}^d$ without changing the asymptotic behaviour., Comment: 13 pages, to be submitted to SciPost Physics

We analyse the functional derivative of the cosmic-shear power spectrum $C_\ell^\gamma$ with respect to the cosmic expansion function. Our interest in doing so is two-fold: (i) In view of attempts to detect minor changes of the cosmic expansion function which may be due to a possibly time-dependent dark-energy density, we wish to know how sensitive the weak-lensing power spectrum is to changes in the expansion function. (ii) In view of recent empirical determinations of the cosmic expansion function from distance measurements, independent of specific cosmological models, we wish to find out how uncertainties in the expansion function translate to uncertainties in the cosmic-shear power spectrum. We find the following answers: Relative changes of the expansion function are amplified by the cosmic-shear power spectrum by a factor $\approx 4-6$, weakly depending on the scale factor where the change is applied, and the current uncertainty of one example for an empirically determined expansion function translates to a relative uncertainty of the cosmic-shear power spectrum of $\approx15\,\%$., Comment: 15 pages, 12 figures, to be submitted to SciPost Astronomy

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., Comment: 8 pages, 5 figues, published version

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, 5 pages, 3 figures

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., Comment: 10 pages, 8 figures

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$., Comment: Withdrawn due to issues with the used probability distribution function and the application of the excursion-set formalism

We apply a generic formalism of light propagation to linearly perturbed spherically symmetric dust models including a cosmological constant. For a comoving observer on the central worldline, we derive the equation of geodesic deviation and perform a suitable spherical harmonic decomposition. This allows to map the abstract gauge-invariant perturbation variables to well-known quantities from weak gravitational lensing like convergence or cosmic shear. The resulting set of differential equations can effectively be solved by a Green's function approach leading to line-of-sight integrals sourced by the perturbation variables on the backward lightcone. The resulting spherical harmonic coefficients of the lensing observables are presented and the shear field is decomposed into its E- and B-modes. Results of this work are an essential tool to add information from linear structure formation to the analysis of spherically symmetric dust models with the purpose of testing the Copernican Principle with multiple cosmological probes., 18 pages, 1 figure, prepared for submission to JCAP

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., 8 pages, 7 figures, submitted for publication to Astronomy & Astrophysics

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

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

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., Comment: 9 pages, 2 figures

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., Comment: 9 pages, 1 figure

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., Comment: 16 pages, 15 figures

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., 7 pages, 4 figures. Submitted

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 =-1.895^{+0.003}_{-0.004}, and \sigma=0.501\pm0.004 (95% confidence levels), or =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], Comment: 15 pages, 11 figures; MNRAS in press; V3 includes minor text updates

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., 9 pages, 7 figures; submitted to A&A

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, Comment: Submitted to the Astrophysical Journal. 19 pages, 14 figures

We show the effectiveness of strong lensing in the characterisation of Lyman continuum emission from faint L~ 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, 5 pages, 2 figures, letter accepted for publication in MNRAS

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., Comment: Submitted to the Nature Journal. 39 Pages, 13 figures

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., Comment: 10 pages, 4 figures

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., Comment: 11 pages, MNRAS accepted

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=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., Comment: 18 pages, 15 figures, 1 table; V2 accepted to MNRAS, includes a significant revision, in particular a new discussion of the results

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, Comment: 10 pages, 8 figures, Accepted by MNRAS

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., Accepted for publication in the Astrophysical Journal Supplements, 22 pages, 16 figures. Updated Tables 3,4,8 and figures 6 and 8 to reflect replacement of Abell 963 with Abell 1423 in CLASH survey. A963 cannot be observed with WFC3 due to the lack of usable guide stars

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., Accepted for publication in MNRAS - minor changes

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., 20 pages, 9 figures, Accepted by MNRAS, A full resolution version can be found at ftp://wise-ftp.tau.ac.il/pub/assafh/arc_stat_paper.pdf

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., Comment: LaTeX file, 10 pages, MNRAS in press, scales larger than the Millennium Simulation box size semi-analytically added, maps changed, references added, typos corrected

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., 14 pages, 16 figures, 3 tables, submitted to MNRAS, 05.Jul.2004

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., Comment: 10 pages, 8 figures, submitted to MNRAS

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., gravitational lensing winter school, Aussois, la vieille Europe, 2003

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., Comment: 14 pages, 6 figures, to appear in "Gravitational Lensing: A Unique Tool For Cosmology", Aussois 2003, eds. D. Valls-Gabaud & J.-P. Kneib

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., Invited Review at "Matter and Energy in Clusters of Galaxies", Taipei 2002 (7 pages, no figures)

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., 5 pages, Latex using mrs2001 style, to appear in "Where's the Matter? Tracing Dark and Bright Matter with the New Generation of Large Scale Surveys", Eds. Treyer & Tresse, 2001, Frontier Group

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 ��$m can be larger than the unlensed ones by several orders of magnitude at flux densities $\gta100 $mJy., 13 pages including 14 figures, MNRAS submitted

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., 2 pages including 2 figures, uses PASP conference proceedings style. Poster contribution to "Gravitational Lensing: Recent Progress and Future Goals", Boston University, July 1999, ed. T.G. Brainerd and C.S. Kochanek

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., submitted to Astronomy & Astrophysics

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., accepted by Astronomy & Astrophysics

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., 4 pages including 4 figures. Accepted for publication in Astronomy & Astrophysics

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., Comment: 14 pages incuding 17 figures. Accepted for publication in Astronomy & Astrophysics

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., Comment: 4 pages, 5 figures, o appear in proceedings of the Ringberg workshop on "Diffuse Thermal and Relativistic Plasma in Galaxy Clusters", April 1999, H. Boehringer, L. Feretti, P. Schuecker (eds.)

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., Comment: 19 pages including 7 postscript figures; submitted to Astronomy and Astrophysics

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., submitted to MNRAS

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., 15 pages including figures; submitted to Astronomy & Astrophysics

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., Comment: 29 pages, LaTex, including 9 figures, submitted to Astronomy and Astrophysics

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., Comment: revised version: references updated, some new results included; 53 pages without any figures; complete versions can be found at http://www.mpa-garching.mpg.de/Lenses/Preprints/JeruLect.html

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

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 (z10^{21} cm^{-2} and redshifts 0.5, Comment: 20 pages, uuencoded compressed postscript file, 6 figures included