Unbiased and precise mass calibration of galaxy clusters is crucial to fully exploit galaxy clusters as cosmological probes. Stacking of weak lensing signal allows us to measure observable-mass relations down to less massive halos halos without extrapolation. We propose a Bayesian inference method to constrain the intrinsic scatter of the mass proxy in stacked analyses. The scatter of the stacked data is rescaled with respect to the individual scatter based on the number of binned clusters. We apply this method to the galaxy clusters detected with the AMICO (Adaptive Matched Identifier of Clustered Objects) algorithm in the third data release of the Kilo-Degree Survey. The results confirm the optical richness as a low scatter mass proxy. Based on the optical richness and the calibrated weak lensing mass-richness relation, mass of individual objects down to ~10^13 solar masses can be estimated with a precision of ~20 per cent.

In this paper, we present resolved observations of the Sunyaev-Zel'dovich (SZ) effect, obtained with the NIKA2 camera, towards the cluster of galaxies XLSSC102, a relatively low-mass system ($M_{500} \sim 2 \times 10^{14}$ M$_{\odot}$) at $z = 0.97$ detected from the XXL survey. We combine NIKA2 SZ data, XMM-Newton X-ray data, and Megacam optical data to explore, respectively, the spatial distribution of the gas electron pressure, the gas density, and the galaxies themselves. We find significant offsets between the X-ray peak, the SZ peak, the brightest cluster galaxy, and the peak of galaxy density. Additionally, the galaxy distribution and the gas present elongated morphologies. This is interpreted as the sign of a recent major merging event, which induced a local boost of the gas pressure towards the north of XLSSC102 and stripped the gas out of the galaxy group. The NIKA2 data are also combined with XXL data to construct the thermodynamic profiles of XLSSC102, obtaining relatively tight constraints up to about $\sim r_{500}$, and revealing properties that are typical of disturbed systems. We also explore the impact of the cluster centre definition and the implication of local pressure substructure on the recovered profiles. Finally, we derive the global properties of XLSSC102 and compare them to those of high-mass-and-low-redshift systems, finding no strong evidence for non-standard evolution. We also use scaling relations to obtain alternative mass estimates from our profiles. The variation between these different mass estimates reflects the difficulty to accurately measure the mass of low-mass clusters at z$\sim$1, especially with low signal-to-noise ratio (S/N) data and for a disturbed system. However, it also highlights the strength of resolved SZ observations alone and in combination with survey-like X-ray data.

We present a weak-lensing analysis of X-ray galaxy groups and clusters selected from the XMM-XXL survey using the first-year data from the Hyper Suprime-Cam (HSC) Subaru Strategic Program. Our joint weak-lensing and X-ray analysis focuses on 136 spectroscopically confirmed X-ray-selected systems at 0.031 < z < 1.033 detected in the 25sqdeg XXL-N region. We characterize the mass distributions of individual clusters and establish the concentration-mass (c-M) relation for the XXL sample, by accounting for selection bias and statistical effects, and marginalizing over the remaining mass calibration uncertainty. We find the mass-trend parameter of the c-M relation to be \beta = -0.07 \pm 0.28 and the normalization to be c200 = 4.8 \pm 1.0 (stat) \pm 0.8 (syst) at M200=10^{14}Msun/h and z = 0.3. We find no statistical evidence for redshift evolution. Our weak-lensing results are in excellent agreement with dark-matter-only c-M relations calibrated for recent LCDM cosmologies. The level of intrinsic scatter in c200 is constrained as \sigma(\ln[c200]) < 24% (99.7% CL), which is smaller than predicted for the full population of LCDM halos. This is likely caused in part by the X-ray selection bias in terms of the relaxation state. We determine the temperature-mass (Tx-M500) relation for a subset of 105 XXL clusters that have both measured HSC lensing masses and X-ray temperatures. The resulting Tx-M500 relation is consistent with the self-similar prediction. Our Tx-M500 relation agrees with the XXL DR1 results at group scales, but has a slightly steeper mass trend, implying a smaller mass scale in the cluster regime. The overall offset in the Tx-M500 relation is at the $1.5\sigma$ level, corresponding to a mean mass offset of (34\pm 20)%. We also provide bias-corrected, weak-lensing-calibrated M200 and M500 mass estimates of individual XXL clusters based on their measured X-ray temperatures.

We report the scaling relations derived by fitting the X-ray parameters determined from analyzing the XMM-Newton observations of 120 galaxy clusters in the Planck Early Sunyaev-Zel'dovich sample spanning the redshift range of 0.059$<$$z$$<$0.546. We find that the slopes of all the investigated scaling relations significantly deviate from the self-similar predictions, if self-similar redshift evolution is assumed. When the redshift evolution is left free to vary, the derived slopes are more in agreement with the self-similar predictions. Relaxed clusters have on average $\sim$30$\%$ higher X-ray luminosity than disturbed clusters at a given mass, a difference that, depending on the relative fraction of relaxed and disturbed clusters in the samples (e.g. SZ vs X-ray selected), have a strong impact in the normalization obtained in different studies. Using the core-excised cluster luminosities reduces the scatter and brings into better agreement the $L$-$M_{tot}$ and $L$-$T$ relations determined for different samples. $M_{tot}$-$T$, $M_{tot}$-$Y_X$, and $M_{tot}$-$M_{gas}$ relations show little dependence on the dynamical state of the clusters, but the normalizations of these relations may depend on the mass range investigated. Although most of the clusters investigated in this work reside at relatively low redshift, the fits prefer values of $\gamma$, the parameter accounting for the redshift evolution, different from the self-similar predictions. This suggests an evolution ($<$2$\sigma$ level, with the exception of the $M_{tot}$-$T$ relation) of the scaling relations. For the first time, we find significant evolution ($>$3$\sigma$) of the $M_{tot}$-$T$ relation, pointing to an increase of the kinetic-to-thermal energy ratio with redshift. This is consistent with a scenario in which higher redshift clusters are on average more disturbed than their lower redshift counterparts.

The Beyond Ultra-deep Frontier Fields and Legacy Observations (BUFFALO) is a 101 orbit + 101 parallel Cycle 25 Hubble Space Telescope Treasury program taking data from 2018-2020. BUFFALO will expand existing coverage of the Hubble Frontier Fields (HFF) in WFC3/IR F105W, F125W, and F160W and ACS/WFC F606W and F814W around each of the six HFF clusters and flanking fields. This additional area has not been observed by HST but is already covered by deep multi-wavelength datasets, including Spitzer and Chandra. As with the original HFF program, BUFFALO is designed to take advantage of gravitational lensing from massive clusters to simultaneously find high-redshift galaxies which would otherwise lie below HST detection limits and model foreground clusters to study properties of dark matter and galaxy assembly. The expanded area will provide a first opportunity to study both cosmic variance at high redshift and galaxy assembly in the outskirts of the large HFF clusters. Five additional orbits are reserved for transient followup. BUFFALO data including mosaics, value-added catalogs and cluster mass distribution models will be released via MAST on a regular basis, as the observations and analysis are completed for the six individual clusters.

Scaling relations trace the formation and evolution of galaxy clusters. We exploited multi-wavelength surveys -- the XXL survey at \emph{XMM-Newton} in the X-ray band, and the Hyper Suprime-Cam Subaru Strategic Program for optical weak lensing -- to study an X-ray selected, complete sample of clusters and groups. The scalings of gas mass, temperature, and soft-band X-ray luminosity with the weak lensing mass show imprints of radiative cooling and AGN feedback in groups. From the multi-variate analysis, we found some evidence for steeper than self-similar slopes for gas mass ($\beta_{m_\text{g}|m}=1.73 \pm0.80$) and luminosity ($\beta_{l|m}=1.91\pm0.94$) and a nearly self-similar slope for the temperature ($\beta_{t|m}=0.78\pm0.43$). Intrinsic scatters of X-ray properties appear to be positively correlated at a fixed mass (median correlation factor $\rho_{X_1X_2|m}\sim0.34$) due to dynamical state and merger history of the halos. Positive correlations with the weak lensing mass (median correlation factor $\rho_{m_\text{wl}X|m}\sim0.35$) can be connected to triaxiality and orientation. Comparison of weak lensing and hydrostatic masses suggests a small role played by non-thermal pressure support ($9\pm17\%$).

Context. Scaling relations between cluster properties embody the formation and evolution of cosmic structure. Intrinsic scatters and correlations between X-ray properties are determined from merger history, baryonic processes, and dynamical state. Aims. We look for an unbiased measurement of the scatter covariance matrix between the three main X-ray observable quantities attainable in large X-ray surveys -- temperature, luminosity, and gas mass. This also gives us the cluster property with the lowest conditional intrinsic scatter at fixed mass. Methods. Intrinsic scatters and correlations can be measured under the assumption that the observable properties of the intra-cluster medium hosted in clusters are log-normally distributed around power-law scaling relations. The proposed method is self-consistent, based on minimal assumptions, and requires neither the external calibration by weak lensing, dynamical, or hydrostatic masses nor the knowledge of the mass completeness. Results. We analyzed the 100 brightest clusters detected in the XXL Survey and their X-ray properties measured within a fixed radius of 300 kpc. The gas mass is the less scattered proxy (~8%). The temperature (~20%) is intrinsically less scattered than the luminosity (~30%) but it is measured with a larger observational uncertainty. We found some evidence that gas mass, temperature and luminosity are positively correlated. Time-evolutions are in agreement with the self-similar scenario, but the luminosity-temperature and the gas mass-temperature relations are steeper. Conclusions. Positive correlations between X-ray properties can be determined by the dynamical state and the merger history of the halos. The slopes of the scaling relations are affected by radiative processes.

We present results from simultaneous modeling of high angular resolution GBT/MUSTANG-2 90 GHz Sunyaev-Zel'dovich effect (SZE) measurements and XMM-XXL X-ray images of three rich galaxy clusters selected from the HSC-SSP Survey. The combination of high angular resolution SZE and X-ray imaging enables a spatially resolved multi-component analysis, which is crucial to understand complex distributions of cluster gas properties. The targeted clusters have similar optical richnesses and redshifts, but exhibit different dynamical states in their member galaxy distributions; a single-peaked cluster, a double-peaked cluster, and a cluster belonging to a supercluster. A large-scale residual pattern in both regular Compton-parameter $y$ and X-ray surface brightness distributions is found in the single-peaked cluster, indicating a sloshing mode. The double-peaked cluster shows an X-ray remnant cool core between two SZE peaks associated with galaxy concentrations. The temperatures of the two peaks reach $\sim20-30$ keV in contrast to the cool core component of $\sim2$ keV, indicating a violent merger. The highest SZE signal for the supercluster is elongated along a direction perpendicular to the major axis of the X-ray core, suggesting a minor merger before core passage. The $S_X$ and $y$ distributions are thus perturbed at some level, regardless of the optical properties. We find that an integrated Compton $y$ parameter and a temperature for the major merger are significantly boosted from those expected by the weak-lensing mass and those for the other two clusters show no significant deviations, which is consistent with predictions of numerical simulations.

We present the on-going activity to characterize the geometrical properties of the gas and dark matter haloes using multi-wavelength observations of galaxy clusters. The role of the SZ signal in describing the gas distribution is discussed for the pilot case of the CLASH object MACS J1206.2-0847. Preliminary images of the NIKA2 and ALMA exposures are presented.

Nearby clusters of galaxies, z<0.1, are cosmic structures still under formation. Understanding the thermodynamic properties of merging clusters can provide crucial information on how they grow in the local universe. A detailed study of the intra-cluster medium (ICM) properties of un-relaxed systems is essential to understand the fate of in-falling structures and, more generally, the virialization process. We analyzed a mosaic of XMM-Newton observations (240 ks) of the Lyra system (z=0.067) that shows a complex dynamical state. We find the main cluster RXC J1825.3+3026 to be in a late merger phase, whereas its companion CIZA J1824.1+3029 is a relaxed cool-core cluster. We estimate a mass ratio of ~1:2 for the pair. No diffuse X-ray emission is found in the region between them, indicating that these clusters are in a pre-merger phase. We found evidence of a galaxy group infalling on RXC J1825.3+3026 in an advanced state of disruption. The Southern Galaxy, one of the brightest galaxies in the Lyra complex, was very likely at the center of the infalling group. This galaxy has a gaseous corona indicating that it was able to retain some of its gas after the ram-pressure stripping of the intra-group medium. In this scenario the diffuse emission excess observed southwest of RXC J1825.3+3026 could be due to gas once belonging to the group and/or to cluster ICM dislocated by the passage of the group. Finally, we identified three high-velocity galaxies aligned between RXC J1825.3+3026 and the SG, two of these showing evidence of gas stripped from them during infall. We estimate them to be currently falling onto the main cluster at an infall velocity of ~ 3000 km/s. Our study of the Lyra complex provides important clues about the processes presiding over the virialization of massive clusters in the local Universe.

Diffuse radio emission associated with the intra-cluster medium (ICM) is observed in a number of merging galaxy clusters. It is currently believed that in mergers a fraction of the kinetic energy is channeled into non-thermal components, such as turbulence, cosmic rays and magnetic fields, that may lead to the formation of giant synchrotron sources in the ICM. Studying merging galaxy clusters in different evolutionary phases is fundamental to understanding the origin of radio emission in the ICM. We observed the nearby galaxy cluster pair RXC J1825.3+3026 ($z\sim0.065$) and CIZA J1824.1+3029 ($z\sim0.071$) at 120-168 MHz with the LOw Frequency ARray (LOFAR) and made use of a deep (240 ks) XMM-Newton dataset to study the non-thermal and thermal properties of the system. RXC J1825.3+3026 is in a complex dynamical state, with a primary on-going merger in the E-W direction and a secondary later stage merger with a group of galaxies in the SW, while CIZA J1824.1+3029 is dynamically relaxed. These two clusters are in a pre-merger phase. We report the discovery of a Mpc-scale radio halo with a low surface brightness extension in RXC J1825.3+3026 that follows the X-ray emission from the cluster center to the remnant of a galaxy group in the SW. This is among the least massive systems and the faintest giant radio halo known to date. Contrary to this, no diffuse radio emission is observed in CIZA J1824.1+3029 nor in the region between the pre-merger cluster pair. The power spectra of the X-ray surface brightness fluctuations of RXC J1825.3+3026 and CIZA J1824.1+3029 are in agreement with the findings for clusters exhibiting a radio halo and the ones where no radio emission has been detected, respectively. We provide quantitative support to the idea that cluster mergers play a crucial role in the generation of non-thermal components in the ICM.

In this Letter, we report the discovery of a radio halo in the high-redshift galaxy cluster PSZ2 G099.86+58.45 ($z=0.616$) with the LOw Frequency ARray (LOFAR) at 120-168 MHz. This is one of the most distant radio halos discovered so far. The diffuse emission extends over $\sim$ 1 Mpc and has a morphology similar to that of the X-ray emission as revealed by XMM-Newton data. The halo is very faint at higher frequencies and is barely detected by follow-up 1-2 GHz Karl G.~Jansky Very Large Array (JVLA) observations, which enable us to constrain the radio spectral index to be $\alpha\leq 1.5-1.6$, i.e.; with properties between canonical and ultra-steep spectrum radio halos. Radio halos are currently explained as synchrotron radiation from relativistic electrons that are re-accelerated in the intra-cluster medium (ICM) by turbulence driven by energetic mergers. We show that in such a framework radio halos are expected to be relatively common at $\sim150$ MHz ($\sim30-60\%$) in clusters with mass and redshift similar to PSZ2 G099.86+58.45; however, at least 2/3 of these radio halos should have steep spectrum and thus be very faint above $\sim 1$ GHz frequencies. Furthermore, since the luminosity of radio halos at high redshift depends strongly on the magnetic field strength in the hosting clusters, future LOFAR observations will also provide vital information on the origin and amplification of magnetic fields in galaxy clusters.

Galaxy cluster counts in bins of mass and redshift have been shown to be a competitive probe to test cosmological models. This method requires an efficient blind detection of clusters from surveys with a well-known selection function and robust mass estimates. The Euclid wide survey will cover 15000 deg$^2$ of the sky in the optical and near-infrared bands, down to magnitude 24 in the $H$-band. The resulting data will make it possible to detect a large number of galaxy clusters spanning a wide-range of masses up to redshift $\sim 2$. This paper presents the final results of the Euclid Cluster Finder Challenge (CFC). The objective of these challenges was to select the cluster detection algorithms that best meet the requirements of the Euclid mission. The final CFC included six independent detection algorithms, based on different techniques, such as photometric redshift tomography, optimal filtering, hierarchical approach, wavelet and friend-of-friends algorithms. These algorithms were blindly applied to a mock galaxy catalog with representative Euclid-like properties. The relative performance of the algorithms was assessed by matching the resulting detections to known clusters in the simulations. Several matching procedures were tested, thus making it possible to estimate the associated systematic effects on completeness to $<3$%. All the tested algorithms are very competitive in terms of performance, with three of them reaching $>80$% completeness for a mean purity of 80% down to masses of $10^{14}$ M$_{\odot}$ and up to redshift $z=2$. Based on these results, two algorithms were selected to be implemented in the Euclid pipeline, the AMICO code, based on matched filtering, and the PZWav code, based on an adaptive wavelet approach. [abridged]

The possibility of determining the value of the Hubble constant using observations of galaxy clusters in X-ray and microwave wavelengths through the Sunyaev Zel\'dovich (SZ) effect has long been known. Previous measurements have been plagued by relatively large errors in the observational data and severe biases induced, for example, by cluster triaxiality and clumpiness. The advent of \textit{Planck} allows us to map the Compton parameter y, that is, the amplitude of the SZ effect, with unprecedented accuracy at intermediate cluster-centric radii, which in turn allows performing a detailed spatially resolved comparison with X-ray measurements. Given such higher quality observational data, we developed a Bayesian approach that combines informed priors on the physics of the intracluster medium obtained from hydrodynamical simulations of massive clusters with measurement uncertainties. We apply our method to a sample of 61 galaxy clusters with redshifts up to z < 0.5 observed with Planck and XMM-Newton observations and find H_0=67 \pm 3 km s^{-1} Mpc^{-1}.

We present the mass calibration for galaxy clusters detected with the AMICO code in KiDS DR3 data. The cluster sample comprises $\sim$ 7000 objects and covers the redshift range 0.1 < $z$ < 0.6. We perform a weak lensing stacked analysis by binning the clusters according to redshift and two different mass proxies provided by AMICO, namely the amplitude $A$ (measure of galaxy abundance through an optimal filter) and the richness $\lambda^*$ (sum of membership probabilities in a consistent radial and magnitude range across redshift). For each bin, we model the data as a truncated NFW profile plus a 2-halo term, taking into account uncertainties related to concentration and miscentring. From the retrieved estimates of the mean halo masses, we construct the $A$-$M_{200}$ and the $\lambda^*$-$M_{200}$ relations. The relations extend over more than one order of magnitude in mass, down to $M_{200} \sim 2 (5) \times 10^{13} M_\odot/h$ at $z$ = 0.2 (0.5), with small evolution in redshift. The logarithmic slope is $\sim 2.0$ for the $A$-mass relation, and $\sim 1.7$ for the $\lambda^*$-mass relation, consistent with previous estimations on mock catalogues and coherent with the different nature of the two observables.

We present the reconstruction of hydrostatic mass profiles in 13 X-ray luminous galaxy clusters that have been mapped in their X-ray and SZ signal out to $R_{200}$ for the XMM-Newton Cluster Outskirts Project (X-COP). Using profiles of the gas temperature, density and pressure that have been spatially resolved out to (median value) 0.9 $R_{500}$, 1.8 $R_{500}$, and 2.3 $R_{500}$, respectively, we are able to recover the hydrostatic gravitating mass profile with several methods and using different mass models. The hydrostatic masses are recovered with a relative (statistical) median error of 3% at $R_{500}$ and 6% at $R_{200}$. By using several different methods to solve the equation of the hydrostatic equilibrium, we evaluate some of the systematic uncertainties to be of the order of 5% at both $R_{500}$ and $R_{200}$. A Navarro-Frenk-White profile provides the best-fit in nine cases out of 13, with the remaining four cases that do not show a statistically significant tension with it. The distribution of the mass concentration follows the correlations with the total mass predicted from numerical simulations with a scatter of 0.18 dex, with an intrinsic scatter on the hydrostatic masses of 0.15 dex. We compare them with the estimates of the total gravitational mass obtained through X-ray scaling relations applied to $Y_X$, gas fraction and $Y_{SZ}$, and from weak lensing and galaxy dynamics techniques, and measure a substantial agreement with the results from scaling laws, from WL at both $R_{500}$ and $R_{200}$ (with differences below 15%), from cluster velocity dispersions, but a significant tension with the caustic masses that tend to underestimate the hydrostatic masses by 40% at $R_{200}$. We also compare these measurements with predictions from alternative models to the Cold Dark Matter, like the "Emergent Gravity" and MOND scenarios.

We present the first catalogue of galaxy cluster candidates derived from the third data release of the Kilo Degree Survey (KiDS-DR3). The sample of clusters has been produced using the Adaptive Matched Identifier of Clustered Objects (AMICO) algorithm. In this analysis AMICO takes advantage of the luminosity and spatial distribution of galaxies only, not considering colours. In this way, we prevent any selection effect related to the presence or absence of the red-sequence in the clusters. The catalogue contains 7988 candidate galaxy clusters in the redshift range 0.1<z<0.8 down to S/N>3.5 with a purity approaching 95% over the entire redshift range. In addition to the catalogue of galaxy clusters we also provide a catalogue of galaxies with their probabilistic association to galaxy clusters. We quantify the sample purity, completeness and the uncertainties of the detection properties, such as richness, redshift, and position, by means of mock galaxy catalogues derived directly from the data. This preserves their statistical properties including photo-z uncertainties, unknown absorption across the survey, missing data, spatial correlation of galaxies and galaxy clusters. Being based on the real data, such mock catalogues do not have to rely on the assumptions on which numerical simulations and semi-analytic models are based on. This paper is the first of a series of papers in which we discuss the details and physical properties of the sample presented in this work.

Context. We present an estimation of cosmological parameters with clusters of galaxies. Aims. We constrain the $\Omega_m$, $\sigma_8$, and $w$ parameters from a stand-alone sample of X-ray clusters detected in the 50 deg$^2$ XMM-XXL survey with a well-defined selection function. Methods. We analyse the redshift distribution of a sample comprising 178 high S/N clusters out to a redshift of unity. The cluster sample scaling relations are determined in a self-consistent manner. Results. In a lambda cold dark matter ($\Lambda$CDM) model, the cosmology favoured by the XXL clusters compares well with results derived from the Planck S-Z clusters for a totally different sample (mass/redshift range, selection biases, and scaling relations). However, with this preliminary sample and current mass calibration uncertainty, we find no inconsistency with the Planck CMB cosmology. If we relax the $w$ parameter, the Planck CMB uncertainties increase by a factor of $\sim$10 and become comparable with those from XXL clusters. Combining the two probes allows us to put constraints on $\Omega_m$=0.316$\pm$0.060, $\sigma_8$=0.814$\pm$0.054, and $w$=-1.02$\pm$0.20. Conclusions. This first self-consistent cosmological analysis of a sample of serendipitous XMM clusters already provides interesting insights into the constraining power of the XXL survey. Subsequent analysis will use a larger sample extending to lower confidence detections and include additional observable information, potentially improving posterior uncertainties by roughly a factor of 3.

Galaxy clusters trace the highest density peaks in the large-scale structure of the Universe. Their clustering provides a powerful probe that can be exploited in combination with cluster mass measurements to strengthen the cosmological constraints provided by cluster number counts. We investigate the spatial properties of a homogeneous sample of X-ray selected galaxy clusters from the XXL survey, the largest programme carried out by the XMM-Newton satellite. The measurements are compared to $\Lambda$-cold dark matter predictions, and used in combination with self-calibrated mass scaling relations to constrain the effective bias of the sample, $b_{eff}$, and the matter density contrast, $\Omega_{\rm M}$. We measured the angle-averaged two-point correlation function of the XXL cluster sample. The analysed catalogue consists of $182$ X-ray selected clusters from the XXL second data release, with median redshift $\langle z \rangle=0.317$ and median mass $\langle M_{500} \rangle\simeq1.3\cdot10^{14} M_\odot$. A Markov chain Monte Carlo analysis is performed to extract cosmological constraints using a likelihood function constructed to be independent of the cluster selection function. Modelling the redshift-space clustering in the scale range $10<r\,[$\Mpch$]<40$, we obtain $\Omega_{\rm M}=0.27_{-0.04}^{+0.06}$ and $b_{eff}=2.73_{-0.20}^{+0.18}$. This is the first time the two-point correlation function of an X-ray selected cluster catalogue at such relatively high redshifts and low masses has been measured. The XXL cluster clustering appears fully consistent with standard cosmological predictions. The analysis presented in this work demonstrates the feasibility of a cosmological exploitation of the XXL cluster clustering, paving the way for a combined analysis of XXL cluster number counts and clustering.

Galaxy clusters form at the highest density nodes of the cosmic web. The clustering of massive halos is enhanced relative to the general mass distribution and matter beyond the virial region is strongly correlated to the halo mass (halo bias). Clustering can be further enhanced depending on halo properties other than mass (secondary bias). The questions of how much and why the regions surrounding rich clusters are over-dense are still unanswered. Here, we report the analysis of the environment bias in a sample of very massive clusters, selected through the Sunyaev-Zel'dovich effect by the Planck mission. We present the first detection of the correlated dark matter associated to a single cluster, PSZ2 G099.86+58.45. The system is extremely rare in the current paradigm of structure formation. The gravitational lensing signal was traced up to 30 megaparsecs with high signal-to-noise ratio ~3.4. The measured shear is very large and points at environment matter density in notable excess of the cosmological mean. The boosting of the correlated dark matter density around high mass halos can be very effective. Together with ensemble studies of the large scale structure, lensing surveys can picture the surroundings of single haloes.

The main goal of this study is to investigate the LF of a sample of 142 X-ray selected clusters, with spectroscopic redshift confirmation and a well defined selection function, spanning a wide redshift and mass range, and to test the LF dependence on cluster global properties, in a homogeneous and unbiased way. Our study is based on the Canada-France-Hawaii Telescope Legacy Survey (CFHTLS) photometric galaxy catalogue,associated with photometric redshifts. We constructed LFs inside a scaled radius using a selection in photometric redshift around the cluster spectroscopic redshift in order to reduce projection effects. The width of the photometric redshift selection was carefully determined to avoid biasing the LF and depended on both the cluster redshift and the galaxy magnitudes. The purity was then enhanced by applying a precise background subtraction. We constructed composite luminosity functions (CLFs) by stacking the individual LFs and studied their evolution with redshift and richness, analysing separately the brightest cluster galaxy (BCG) and non-BCG members. We fitted the dependences of the CLFs and BCG distributions parameters with redshift and richness conjointly in order to distinguish between these two effects. We find that the usual photometric redshift selection methods can bias the LF estimate if the redshift and magnitude dependence of the photometric redshift quality is not taken into account. Our main findings concerning the evolution of the galaxy luminosity distribution with redshift and richness are that, in the inner region of clusters and in the redshift-mass range we probe (about $0<z<1$ and $10^{13} M_{\odot}<M_{500}<5\times10^{14}M_{\odot}$), the bright part of the LF (BCG excluded) does not depend much on mass or redshift except for its amplitude, whereas the BCG luminosity increases both with redshift and richness, and its scatter decreases with redshift.

We reconstruct the two-dimensional (2D) matter distributions in 20 high-mass galaxy clusters selected from the CLASH survey by using the new approach of performing a joint weak lensing analysis of 2D shear and azimuthally averaged magnification measurements. This combination allows for a complete analysis of the field, effectively breaking the mass-sheet degeneracy. In a Bayesian framework, we simultaneously constrain the mass profile and morphology of each individual cluster assuming an elliptical Navarro-Frenk-White halo characterized by the mass, concentration, projected axis ratio, and position angle of the projected major axis.. We find that spherical mass estimates of the clusters from azimuthally averaged weak-lensing measurements in previous work are in excellent agreement with our results from a full 2D analysis. Combining all 20 clusters in our sample, we detect the elliptical shape of weak-lensing halos at the $5\sigma$ significance level within a scale of 2Mpc$/h$. The median projected axis ratio is $0.67\pm 0.07$ at a virial mass of $M_\mathrm{vir}=(15.2\pm 2.8) \times 10^{14} M_\odot$, which is in agreement with theoretical predictions of the standard collisionless cold dark matter model. We also study misalignment statistics of the brightest cluster galaxy, X-ray, thermal Sunyaev-Zel'dovich effect, and strong-lensing morphologies with respect to the weak-lensing signal. Among the three baryonic tracers studied here, we find that the X-ray morphology is best aligned with the weak-lensing mass distribution, with a median misalignment angle of $21\pm 7$ degrees. We also conduct a stacked quadrupole shear analysis assuming that the X-ray major axis is aligned with that of the projected mass distribution. This yields a consistent axis ratio of $0.67\pm 0.10$, suggesting again a tight alignment between the intracluster gas and dark matter.

We present a new cosmological probe for galaxy clusters, the halo sparsity. This characterises halos in terms of the ratio of halo masses measured at two different radii and carries cosmological information encoded in the halo mass profile. Building upon the work of Balmes et al. (2014) we test the properties of the sparsity using halo catalogs from a numerical N-body simulation of ($2.6$ Gpc/h)$^3$ volume with $4096^3$ particles. We show that at a given redshift the average sparsity can be predicted from prior knowledge of the halo mass function. This provides a quantitative framework to infer cosmological parameter constraints using measurements of the sparsity of galaxy clusters. We show this point by performing a likelihood analysis of synthetic datasets with no systematics, from which we recover the input fiducial cosmology. We also perform a preliminary analysis of potential systematic errors and provide an estimate of the impact of baryonic effects on sparsity measurements. We evaluate the sparsity for a sample of 104 clusters with hydrostatic masses from X-ray observations and derive constraints on the cosmic matter density $\Omega_m$ and the normalisation amplitude of density fluctuations at the $8$ Mpc h$^{-1}$ scale, $\sigma_8$. Assuming no systematics, we find $\Omega_m=0.42\pm 0.17$ and $\sigma_8=0.80\pm 0.31$ at $1\sigma$, corresponding to $S_8\equiv \sigma_8\sqrt{\Omega_m}=0.48\pm 0.11$. Future cluster surveys may provide opportunities for precise measurements of the sparsity. A sample of a few hundreds clusters with mass estimate errors at a few percent level can provide competitive cosmological parameter constraints complementary to those inferred from other cosmic probes.

The $\Lambda$CDM model of structure formation makes strong predictions on concentration and shape of DM (dark matter) halos, which are determined by mass accretion processes. Comparison between predicted shapes and observations provides a geometric test of the $\Lambda$CDM model. Accurate and precise measurements needs a full three-dimensional analysis of the cluster mass distribution. We accomplish this with a multi-probe 3D analysis of the X-ray regular CLASH (Cluster Lensing And Supernova survey with Hubble) clusters combining strong and weak lensing, X-ray photometry and spectroscopy, and the Sunyaev-Zel'dovich effect. The cluster shapes and concentrations are consistent with $\Lambda$CDM predictions. The CLASH clusters are randomly oriented, as expected given the sample selection criteria. Shapes agree with numerical results for DM-only halos, which hints at baryonic physics being not so effective in making halos rounder.

We perform a three-dimensional triaxial analysis of 16 X-ray regular and 4 high-magnification galaxy clusters selected from the CLASH survey by combining two-dimensional weak-lensing and central strong-lensing constraints. In a Bayesian framework, we constrain the intrinsic structure and geometry of each individual cluster assuming a triaxial Navarro-Frenk-White halo with arbitrary orientations, characterized by the mass $M_{200\mathrm{c}}$, halo concentration $C_{200\mathrm{c}}$, and triaxial axis ratios ($q_{\mathrm{a}} \le q_{\mathrm{b}}$), and investigate scaling relations between these halo structural parameters. From triaxial modeling of the X-ray-selected subsample, we find that the halo concentration decreases with increasing cluster mass, with a mean concentration of $C_{200\mathrm{c}} = 4.82\pm0.30$ at the pivot mass $M_{200\mathrm{c}}=10^{15}M_{\odot}h^{-1}$. This is consistent with the result from spherical modeling, $C_{200\mathrm{c}}=4.51\pm 0.14$. Independently of the priors, the minor-to-major axis ratio $q_{\mathrm{a}}$ of our full sample exhibits a clear deviation from the spherical configuration ($q_{\mathrm{a}}=0.52 \pm 0.04$ at $10^{15}M_{\odot}h^{-1}$ with uniform priors), with a weak dependence on the cluster mass. Combining all 20 clusters, we obtain a joint ensemble constraint on the minor-to-major axis ratio of $q_{\mathrm{a}}=0.652^{+0.162}_{-0.078}$ and a lower bound on the intermediate-to-major axis ratio of $q_{\mathrm{b}}>0.63$ at the $2\sigma$ level from an analysis with uniform priors. Assuming priors on the axis ratios derived from numerical simulations, we constrain the degree of triaxiality for the full sample to be $\mathcal{T}=0.79 \pm 0.03$ at $10^{15}M_{\odot}h^{-1}$, indicating a preference for a prolate geometry of cluster halos. We find no statistical evidence for an orientation bias ($f_{\mathrm{geo}}=0.93 \pm 0.07$) (abridged)

In their "How proper are Bayesian models in the astronomical literature?" [arXiv:1712.03549], Hyungsuk Tak, Sujit K. Ghosh and Justin A. Ellis criticised my work with false statements. This is an infamous case of straw man fallacy. They give the impression of refuting an opponent's argument, while they refute an argument that was not presented.

An X-ray survey with the XMM-Newton telescope, XMM-XXL, has identified hundreds of galaxy groups and clusters in two 25 deg$^2$ fields. Combining spectroscopic and X-ray observations in one field, we determine how the kinetic energy of galaxies scales with hot gas temperature and also, by imposing prior constraints on the relative energies of galaxies and dark matter, infer a power-law scaling of total mass with temperature. Our goals are: i) to determine parameters of the scaling between galaxy velocity dispersion and X-ray temperature, $T_{\rm 300kpc}$, for the halos hosting XXL-selected clusters, and; ii) to infer the log-mean scaling of total halo mass with temperature, $\langle \ln M_{200} \, | \, T, z \rangle$. We apply an ensemble velocity likelihood to a sample of $> 1500$ spectroscopic redshifts within $132$ spectroscopically confirmed clusters with redshifts $z < 0.6$ to model, $\langle \ln \sigma_{\rm gal}\,|\,T,z\rangle$, where $\sigma_{\rm gal}$ is the velocity dispersion of XXL cluster member galaxies and $T$ is a 300 kpc aperture temperature. To infer total halo mass we use a precise virial relation for massive halos calibrated by N-body simulations along with a single degree of freedom summarizing galaxy velocity bias with respect to dark matter. For the XXL-N cluster sample, we find $\sigma_{\rm gal} \propto T^{0.63\pm0.05}$, a slope significantly steeper than the self-similar expectation of $0.5$. Assuming scale-independent galaxy velocity bias, we infer a mean logarithmic mass at a given X-ray temperature and redshift, $\langle\ln (E(z) M_{200}/10^{14}\,{\rm M}_{\odot})|T,z\rangle=\pi+\alpha \ln(T/T_p )+\beta\ln (E(z)/E(z_p) )$ using pivot values ${\rm k}T_{p}=2.2\,{\rm keV}$ and $z_p=0.25$, with normalization $\pi=0.45\pm0.24$ and slope $\alpha=1.89\pm0.15$. We obtain only weak constraints on redshift evolution, $\beta=-1.29\pm1.14$.

The fraction of galaxies bound in groups in the nearby Universe is high (50% at z~0). Systematic studies of galaxy properties in groups are important in order to improve our understanding of the evolution of galaxies and of the physical phenomena occurring within this environment. We have built a complete spectrophotometric sample of galaxies within X-ray detected, optically spectroscopically confirmed groups and clusters (G&C), covering a wide range of halo masses at z<= 0.6. In the context of the XXL survey, we analyse a sample of 164 G&C in the XXL-North region (XXL-N), at z <= 0.6, with a wide range of virial masses (1.24 x 10^13 <=M_500 M_sun <= 6.63 x 10^14) and X-ray luminosities ( 2.27 x 10^41 <= L^XXL_500 (erg/s)<= 2.15 x10^44). The G&C are X-ray selected and spectroscopically confirmed. We describe the membership assignment and the spectroscopic completeness analysis, and compute stellar masses. As a first scientific exploitation of the sample, we study the dependence of the galaxy stellar mass function (GSMF) on global environment. We present a spectrophotometric characterisation of the G&C and their galaxies. The final sample contains 132 G&C, 22111 field galaxies and 2225 G&C galaxies with r-band magnitude <20. Of the G&C, 95% have at least three spectroscopic members, and 70% at least ten. The shape of the GSMF seems not to depend on environment (field versus G&C) or X-ray luminosity ( used as a proxy for the virial mass of the system). These results are confirmed by the study of the correlation between mean stellar mass of G&C members and L^XXL_500.We release the spectrophotometric catalogue of galaxies with all the quantities computed in this work. As a first homogeneous census of galaxies within X-ray spectroscopically confirmed G&C at these redshifts, this sample will allow environmental studies of the evolution of galaxy properties.

The possibly unbiased selection process in surveys of the Sunyaev Zel'dovich effect can unveil new populations of galaxy clusters. We performed a weak lensing analysis of the PSZ2LenS sample, i.e. the PSZ2 galaxy clusters detected by the Planck mission in the sky portion covered by the lensing surveys CFHTLenS and RCSLenS. PSZ2LenS consists of 35 clusters and it is a statistically complete and homogeneous subsample of the PSZ2 catalogue. The Planck selected clusters appear to be unbiased tracers of the massive end of the cosmological haloes. The mass concentration relation of the sample is in excellent agreement with predictions from the Lambda cold dark matter model. The stacked lensing signal is detected at 14 sigma significance over the radial range 0.1<R<3.2 Mpc/h, and is well described by the cuspy dark halo models predicted by numerical simulations. We confirmed that Planck estimated masses are biased low by b_SZ= 27+-11(stat)+-8(sys) per cent with respect to weak lensing masses. The bias is higher for the cosmological subsample, b_SZ= 40+-14+-(stat)+-8(sys) per cent.

The Large Synoptic Survey Telescope (LSST) will enable revolutionary studies of galaxies, dark matter, and black holes over cosmic time. The LSST Galaxies Science Collaboration has identified a host of preparatory research tasks required to leverage fully the LSST dataset for extragalactic science beyond the study of dark energy. This Galaxies Science Roadmap provides a brief introduction to critical extragalactic science to be conducted ahead of LSST operations, and a detailed list of preparatory science tasks including the motivation, activities, and deliverables associated with each. The Galaxies Science Roadmap will serve as a guiding document for researchers interested in conducting extragalactic science in anticipation of the forthcoming LSST era.

The XXL survey currently covers two 25 sq. deg. patches with XMM observations of ~10ks. We summarise the scientific results associated with the first release of the XXL data set, that occurred mid 2016. We review several arguments for increasing the survey depth to 40 ks during the next decade of XMM operations. X-ray (z<2) cluster, (z<4) AGN and cosmic background survey science will then benefit from an extraordinary data reservoir. This, combined with deep multi-$\lambda$ observations, will lead to solid standalone cosmological constraints and provide a wealth of information on the formation and evolution of AGN, clusters and the X-ray background. In particular, it will offer a unique opportunity to pinpoint the z>1 cluster density. It will eventually constitute a reference study and an ideal calibration field for the upcoming eROSITA and Euclid missions.

Galaxy clusters are the most recent products of hierarchical accretion over cosmological scales. The gas accreted from the cosmic field is thermalized inside the cluster halo. Gas entropy and pressure are expected to have a self-similar behaviour with their radial distribution following a power law and a generalized Navarro-Frenk-White profile, respectively. This has been shown also in many different hydrodynamical simulations. We derive the spatially-resolved thermodynamical properties of 47 X-ray galaxy clusters observed with Chandra in the redshift range 0.4 < z < 1.2, the largest sample investigated so far in this redshift range with X-rays spectroscopy, with a particular care in reconstructing the gas entropy and pressure radial profiles. We search for deviation from the self-similar behaviour and look for possible evolution with redshift. The entropy and pressure profiles lie very close to the baseline prediction from gravitational structure formation. We show that these profiles deviate from the baseline prediction as function of redshift, in particular at z > 0.75, where, in the central regions, we observe higher values of the entropy (by a factor of 2.2) and systematically lower estimates (by a factor of 2.5) of the pressure. The effective polytropic index, which retains informations about the thermal distribution of the gas, shows a slight linear positive evolution with the redshift and the concentration of the dark matter distribution. A prevalence of non-cool-core, disturbed systems, as we observe at higher redshifts, can explain such behaviours.

Mass measurements of astronomical objects are most wanted but still elusive. We need them to trace the formation and evolution of cosmic structure but we can get direct measurements only for a minority. This lack can be circumvented with a proxy and a scaling relation. The twofold goal of estimating the unbiased relation and finding the right proxy value to plug in can be hampered by systematics, selection effects, Eddington/Malmquist biases and time evolution. We present a Bayesian hierarchical method which deals with these issues. Masses to be predicted are treated as missing data in the regression and are estimated together with the scaling parameters. The calibration subsample with measured masses does not need to be representative of the full sample as far as it follows the same scaling relation. We apply the method to forecast weak lensing calibrated masses of the Planck, redMaPPer and MCXC clusters. Planck masses are biased low with respect to weak lensing calibrated masses, with a bias more pronounced for high redshift clusters. MCXC masses are under-estimated by ~20 per cent, which may be ascribed to hydrostatic bias. Packages and catalogs are made available with the paper.

Multi-wavelength techniques can probe the distribution and the physical properties of baryons and dark matter in galaxy clusters from the inner regions out to the peripheries. We present a full three-dimensional analysis combining strong and weak lensing, X-ray surface brightness and temperature, and the Sunyaev-Zel'dovich effect. The method is applied to MACS J1206.2-0847, a remarkably regular, face-on, massive, M_{200}=(1.1+-0.2)*10^{15}M_Sun/h, cluster at z=0.44. The measured concentration, c_{200}=6.3+-1.2, and the triaxial shape are common to halos formed in a LCDM scenario. The gas has settled in and follows the shape of the gravitational potential, which is evidence of pressure equilibrium via the shape theorem. There is no evidence for significant non-thermal pressure and the equilibrium is hydrostatic.

Mass is a fundamental property of galaxy groups and clusters. In theory weak gravitational lensing will enable an approximately unbiased measurement of mass, but parametric methods for extracting cluster masses from data require the additional knowledge of concentration. Measurements of both mass and concentration are limited by the degeneracy between the two parameters, particularly in low mass, high redshift systems where the signal-to-noise is low. In this paper we develop a hierarchical model of mass and concentration for mass inference we test our method on toy data and then apply it to a sample of galaxy groups and poor clusters down to masses of $\sim$ 1e13 M$_\odot$. Our fit and model gives a relationship among masses, concentrations and redshift that allow prediction of these parameters from incomplete and noisy future measurements. Additionally the underlying population can be used to infer an observationally based concentration-mass relation. Our method is equivalent to a quasi- stacking approach with the degree of stacking set by the data. We also demonstrate that mass and concentration derived from pure stacking can be offset from the population mean with differing values depending on the method of stacking.

Galaxy clusters are the most recent, gravitationally-bound products of the hierarchical mass accretion over cosmological scales. How the mass is concentrated is predicted to correlate with the total mass in the cluster's halo, with systems at higher mass being less concentrated at given redshift and for any given mass, systems with lower concentration are found at higher redshifts. Through a spatial and spectral X-ray analysis, we reconstruct the total mass profile of 47 galaxy clusters observed with Chandra in the redshift range $0.4<z<1.2$, selected to have no major mergers, to investigate the relation between the mass and the dark matter concentration, and the evolution of this relation with redshift. The sample in exam is the largest one investigated so far at $z>0.4$, and is well suited to provide the first constraint on the concentration--mass relation at $z>0.7$ from X-ray analysis. Under the assumptions that the distribution of the X-ray emitting gas is spherically symmetric and in hydrostatic equilibrium, we combine the deprojected gas density and spectral temperature profiles through the hydrostatic equilibrium equation to recover the parameters that describe a NFW total mass distribution. The comparison with results from weak lensing analysis reveals a very good agreement both for masses and concentrations. Uncertainties are however too large to make any robust conclusion on the hydrostatic bias of these systems. The relation is well described by the form $c \propto M^B (1+z)^C$, with $B=-0.50 \pm 0.20$, $C=0.12 \pm 0.61$ (at 68.3\% confidence), it is slightly steeper than the one predicted by numerical simulations ($B\sim-0.1$) and does not show any evident redshift evolution. We obtain the first constraints on the properties of the concentration--mass relation at $z > 0.7$ from X-ray data, showing a reasonable good agreement with recent numerical predictions.

Recent N-body simulations have shown that Einasto radial profiles provide the most accurate description of dark matter halos. Predictions based on the traditional NFW functional form may fail to describe the structural properties of cosmic objects at the percent level required by precision cosmology. We computed the systematic errors expected for weak lensing analyses of clusters of galaxies if one wrongly models the lens density profile. Even though the NFW fits of observed tangential shear profiles can be excellent, viral masses and concentrations of very massive halos (>~ 10^{15}M_Sun/h) can be over- and underestimated by ~10 per cent, respectively. Misfitting effects also steepen the observed mass-concentration relation, as observed in multi-wavelength observations of galaxy groups and clusters. Based on shear analyses, Einasto and NFW halos can be set apart either with deep observations of exceptionally massive structures (>~ 2\times10^{15}M_Sun/h) or by stacking the shear profiles of thousands of group-sized lenses (>~ 10^{14}M_Sun/h).

We present the XXL Survey, the largest XMM programme totaling some 6.9 Ms to date and involving an international consortium of roughly 100 members. The XXL Survey covers two extragalactic areas of 25 deg2 each at a point-source sensitivity of ~ 5E-15 erg/sec/cm2 in the [0.5-2] keV band (completeness limit). The survey's main goals are to provide constraints on the dark energy equation of state from the space-time distribution of clusters of galaxies and to serve as a pathfinder for future, wide-area X-ray missions. We review science objectives, including cluster studies, AGN evolution, and large-scale structure, that are being conducted with the support of approximately 30 follow-up programmes. We describe the 542 XMM observations along with the associated multi-lambda and numerical simulation programmes. We give a detailed account of the X-ray processing steps and describe innovative tools being developed for the cosmological analysis. The paper provides a thorough evaluation of the X-ray data, including quality controls, photon statistics, exposure and background maps, and sky coverage. Source catalogue construction and multi-lambda associations are briefly described. This material will be the basis for the calculation of the cluster and AGN selection functions, critical elements of the cosmological and science analyses. The XXL multi-lambda data set will have a unique lasting legacy value for cosmological and extragalactic studies and will serve as a calibration resource for future dark energy studies with clusters and other X-ray selected sources. With the present article, we release the XMM XXL photon and smoothed images along with the corresponding exposure maps. The XMM XXL observation list (Table B.1) is available in electronic form at the CDS. The present paper is the first in a series reporting results of the XXL-XMM survey.

Scaling relations of clusters have made them particularly important cosmological probes of structure formation. In this work, we present a comprehensive study of the relation between two profile observables, concentration ($\mathrm{c_{vir}}$) and mass ($\mathrm{M_{vir}}$). We have collected the largest known sample of measurements from the literature which make use of one or more of the following reconstruction techniques: Weak gravitational lensing (WL), strong gravitational lensing (SL), Weak+Strong Lensing (WL+SL), the Caustic Method (CM), Line-of-sight Velocity Dispersion (LOSVD), and X-ray. We find that the concentration-mass (c-M) relation is highly variable depending upon the reconstruction technique used. We also find concentrations derived from dark matter only simulations (at approximately $\mathrm{M_{vir} \sim 10^{14} M_{\odot}}$) to be inconsistent with the WL and WL+SL relations at the $\mathrm{1\sigma}$ level, even after the projection of triaxial halos is taken into account. However, to fully determine consistency between simulations and observations, a volume-limited sample of clusters is required, as selection effects become increasingly more important in answering this. Interestingly, we also find evidence for a steeper WL+SL relation as compared to WL alone, a result which could perhaps be caused by the varying shape of cluster isodensities, though most likely reflects differences in selection effects caused by these two techniques. Lastly, we compare concentration and mass measurements of individual clusters made using more than one technique, highlighting the magnitude of the potential bias which could exist in such observational samples.

Linear regression is common in astronomical analyses. I discuss a Bayesian hierarchical modeling of data with heteroscedastic and possibly correlated measurement errors and intrinsic scatter. The method fully accounts for time evolution. The slope, the normalization, and the intrinsic scatter of the relation can evolve with the redshift. The intrinsic distribution of the independent variable is approximated using a mixture of Gaussian distributions whose means and standard deviations depend on time. The method can address scatter in the measured independent variable (a kind of Eddington bias), selection effects in the response variable (Malmquist bias), and departure from linearity in form of a knee. I tested the method with toy models and simulations and quantified the effect of biases and inefficient modeling. The R-package LIRA (LInear Regression in Astronomy) is made available to perform the regression.

We perform a 3D multi-probe analysis of the rich galaxy cluster A1689 by combining improved weak-lensing data from new BVRi'z' Subaru/Suprime-Cam observations with strong-lensing, X-ray, and Sunyaev-Zel'dovich effect (SZE) data sets. We reconstruct the projected matter distribution from a joint weak-lensing analysis of 2D shear and azimuthally integrated magnification constraints, the combination of which allows us to break the mass-sheet degeneracy. The resulting mass distribution reveals elongation with axis ratio ~0.7 in projection. When assuming a spherical halo, our full weak-lensing analysis yields a projected concentration of $c_{200c}^{2D}=8.9\pm 1.1$ ($c_{vir}^{2D}\sim 11$), consistent with and improved from earlier weak-lensing work. We find excellent consistency between weak and strong lensing in the region of overlap. In a parametric triaxial framework, we constrain the intrinsic structure and geometry of the matter and gas distributions, by combining weak/strong lensing and X-ray/SZE data with minimal geometric assumptions. We show that the data favor a triaxial geometry with minor-major axis ratio 0.39+/-0.15 and major axis closely aligned with the line of sight (22+/-10 deg). We obtain $M_{200c}=(1.2\pm 0.2)\times 10^{15} M_{\odot}/h$ and $c_{200c}=8.4\pm 1.3$, which overlaps with the $>1\sigma$ tail of the predicted distribution. The shape of the gas is rounder than the underlying matter but quite elongated with minor-major axis ratio 0.60+/-0.14. The gas mass fraction within 0.9Mpc is 10^{+3}_{-2}%. The thermal gas pressure contributes to ~60% of the equilibrium pressure, indicating a significant level of non-thermal pressure support. When compared to Planck's hydrostatic mass estimate, our lensing measurements yield a spherical mass ratio of $M_{Planck}/M_{GL}=0.70\pm 0.15$ and $0.58\pm 0.10$ with and without corrections for lensing projection effects, respectively.

The measurement of the mass of clusters of galaxies is crucial for their use in cosmology and astrophysics. Masses can be efficiently determined with weak lensing (WL) analyses. I compiled Literature Catalogs of weak Lensing Clusters (LC$^2$). Cluster identifiers, coordinates, and redshifts have been standardised. WL masses were reported to over-densities of 2500, 500, 200, and to the virial one in the reference $\Lambda$CDM model. Duplicate entries were carefully handled. I produced three catalogs: LC$^2$-single, with 485 unique groups and clusters analysed with the single-halo model; LC$^2$-substructure, listing substructures in complex systems; LC$^2$-all, listing all the 822 WL masses found in literature. The catalogs and future updates are publicly available at http://pico.bo.astro.it/~sereno/CoMaLit/LC2/.

We discuss the scaling relation between mass and integrated Compton parameter of a sample of galaxy clusters from the all-sky {\it Planck} Sunyaev-Zel'dovich catalogue. Masses were measured with either weak lensing, caustics techniques, or assuming hydrostatic equilibrium. The retrieved $Y_{500}$-$M_{500}$ relation does not strongly depend on the calibration sample. We found a slope of 1.4-1.9, in agreement with self-similar predictions, with an intrinsic scatter of $20\pm10$ per cent. The absolute calibration of the relation can not be ascertained due to systematic differences of $\sim$20-40 per cent in mass estimates reported by distinct groups. Due to the scatter, the slope of the conditional scaling relation, to be used in cosmological studies of number counts, is shallower, $\sim$1.1-1.6. The regression methods employed account for intrinsic scatter in the mass measurements too. We found that Planck mass estimates suffer from a mass dependent bias.

The first building block to use galaxy clusters in astrophysics and cosmology is the accurate determination of their mass. Two of the most well regarded direct mass estimators are based on weak lensing (WL) determinations or X-ray analyses assuming hydrostatic equilibrium (HE). By comparing these two mass measurements in samples of rich clusters, we determined the intrinsic scatters, $\sigma_\mathrm{WL}\sim$15 per cent for WL masses and $\sigma_\mathrm{HE}\sim$25 per cent for HE masses. The certain assessment of the bias is hampered by differences as large as $\sim$40 per cent in either WL or HE mass estimates reported by different groups. If the intrinsic scatter in the mass estimate is not considered, the slope of any scaling relation `observable--mass' is biased towards shallower values, whereas the intrinsic scatter of the scaling is over-estimated.

The scaling of observable properties of galaxy clusters with mass evolves with time. Assessing the role of the evolution is crucial to study the formation and evolution of massive halos and to avoid biases in the calibration. We present a general method to infer the mass and the redshift dependence, and the time-evolving intrinsic scatter of the mass-observable relations. The procedure self-calibrates the redshift dependent completeness function of the sample. The intrinsic scatter in the mass estimates used to calibrate the relation is considered too. We apply the method to the scaling of mass M_Delta versus line of sight galaxy velocity dispersion sigma_v, optical richness, X-ray luminosity, L_X, and Sunyaev-Zel'dovich signal. Masses were calibrated with weak lensing measurements. The measured relations are in good agreement with time and mass dependencies predicted in the self-similar scenario of structure formation. The lone exception is the L_X-M_Delta relation whose time evolution is negative in agreement with formation scenarios with additional radiative cooling and uniform preheating at high redshift. The intrinsic scatter in the sigma_v-M_Delta relation is notably small, of order of 14 per cent. Robust predictions on the observed properties of the galaxy clusters in the CLASH sample are provided as cases of study. Catalogs and scripts are publicly available at http://pico.bo.astro.it/~sereno/CoMaLit/.

The joint analysis of clustering and stacked gravitational lensing of galaxy clusters in large surveys can constrain the formation and evolution of structures and the cosmological parameters. On scales outside a few virial radii, the halo bias, $b$, is linear and the lensing signal is dominated by the correlated distribution of matter around galaxy clusters. We discuss a method to measure the power spectrum amplitude $\sigma_8$ and $b$ based on a minimal modelling. We considered a sample of $\sim 120000$ clusters photometrically selected from the Sloan Digital Sky Survey in the redshift range $0.1<z<0.6$. The auto-correlation was studied through the two-point function of a subsample of $\sim 70000$ clusters; the matter-halo correlation was derived from the weak lensing signal of the subsample of $\sim 1200$ clusters with Canada-France-Hawaii Lensing Survey data. We obtained a direct measurement of $b$, which increases with mass in agreement with predictions of the $\Lambda$CDM paradigm. Assuming $\Omega_\mathrm{M}=0.3$, we found $\sigma_8=0.78\pm0.16$. We used the same clusters for measuring both lensing and clustering and the estimate of $\sigma_8$ did require neither the mass-richness relation, nor the knowledge of the selection function, nor the modelling of $b$. With an additional theoretical prior on the bias, we obtained $\sigma_8=0.75\pm0.08$.

The relation between mass and concentration of galaxy clusters traces their formation and evolution. Massive lensing clusters were observed to be over-concentrated and following a steep scaling in tension with predictions from the concordance $\Lambda$CDM paradigm. We critically revise the relation in the CLASH, the SGAS, the LOCUSS, and the high-redshift samples of weak lensing clusters. Measurements of mass and concentration are anti-correlated, which can bias the observed relation towards steeper values. We corrected for this bias and compared the measured relation to theoretical predictions accounting for halo triaxiality, adiabatic contraction of the halo, presence of a dominant BCG and, mostly, selection effects in the observed sample. The normalisation, the slope and the scatter of the expected relation are strongly sample-dependent. For the considered samples, the predicted slope is much steeper than that of the underlying relation characterising dark-matter only clusters. We found that the correction for statistical and selection biases in observed relations mostly solve the tension with the $\Lambda$CDM model.

The statistics of peaks in weak lensing convergence maps is a promising tool to investigate both the properties of dark matter haloes and constrain the cosmological parameters. We study how the number of detectable peaks and its scaling with redshift depend upon the cluster dark matter halo profiles and use peak statistics to constrain the parameters of the mass - concentration (MC) relation. We investigate which constraints the Euclid mission can set on the MC coefficients also taking into account degeneracies with the cosmological parameters. To this end, we first estimate the number of peaks and its redshift distribution for different MC relations. We find that the steeper the mass dependence and the larger the normalisation, the higher is the number of detectable clusters, with the total number of peaks changing up to $40\%$ depending on the MC relation. We then perform a Fisher matrix forecast of the errors on the MC relation parameters as well as cosmological parameters. We find that peak number counts detected by Euclid can determine the normalization $A_v$, the mass $B_v$ and redshift $C_v$ slopes and intrinsic scatter $\sigma_v$ of the MC relation to an unprecedented accuracy being $\sigma(A_v)/A_v = 1\%$, $\sigma(B_v)/B_v = 4\%$, $\sigma(C_v)/C_v = 9\%$, $\sigma(\sigma_v)/\sigma_v = 1\%$ if all cosmological parameters are assumed to be known. Should we relax this severe assumption, constraints are degraded, but remarkably good results can be restored setting only some of the parameters or combining peak counts with Planck data. This precision can give insight on competing scenarios of structure formation and evolution and on the role of baryons in cluster assembling. Alternatively, for a fixed MC relation, future peaks counts can perform as well as current BAO and SNeIa when combined with Planck.

Hubble Space Telescope imaging of the galaxy cluster Abell 1689 has revealed an exceptional number of strongly lensed multiply-imaged galaxies, including high-redshift candidates. Previous studies have used this data to obtain the most detailed dark matter reconstructions of any galaxy cluster to date, resolving substructures ~25 kpc across. We examine Abell 1689 (hereafter, A1689) non-parametrically, combining strongly lensed images and weak distortions from wider field Subaru imaging, and we incorporate member galaxies to improve the lens solution. Strongly lensed galaxies are often locally affected by member galaxies, however, these perturbations cannot be recovered in grid based reconstructions because the lensing information is too sparse to resolve member galaxies. By adding luminosity-scaled member galaxy deflections to our smooth grid we can derive meaningful solutions with sufficient accuracy to permit the identification of our own strongly lensed images, so our model becomes self consistent. We identify 11 new multiply lensed system candidates and clarify previously ambiguous cases, in the deepest optical and NIR data to date from Hubble and Subaru. Our improved spatial resolution brings up new features not seen when the weak and strong lensing effects are used separately, including clumps and filamentary dark matter around the main halo. Our treatment means we can obtain an objective mass ratio between the cluster and galaxy components, for examining the extent of tidal stripping of the luminous member galaxies. We find a typical mass-to-light ratios of M/L_B = 21 inside the r<1 arcminute region that drops to M/L_B = 17 inside the r<40 arcsecond region. Our model independence means we can objectively evaluate the competitiveness of stacking cluster lenses for defining the geometric lensing-distance-redshift relation in a model independent way.

We present the observational evidence of the 2-halo term in the stacked shear profile of a sample of about 1200 optically selected galaxy clusters based on imaging data and the public shear catalog from the CFHTLenS. We find that the halo bias, a measure of the correlated distribution of matter around galaxy clusters, has amplitude and correlation with galaxy cluster mass in very good agreement with the predictions based on the LCDM standard cosmological model. The mass-concentration relation is flat but higher than theoretical predictions. We also confirm the close scaling relation between the optical richness of galaxy clusters and their mass.

We present the final analysis of the observational campaign carried out by the PLAN (Pixel Lensing Andromeda) collaboration to detect a dark matter signal in form of MACHOs through the microlensing effect. The campaign consists of about 1 month/year observations carried out during 4 years (2007-2010) at the 1.5m Cassini telescope in Loiano ("Astronomical Observatory of BOLOGNA", OAB) plus 10 days of data taken in 2010 at the 2m Himalayan Chandra Telescope (HCT) monitoring the central part of M31 (two fields of about 13'x12.6'). We establish a fully automated pipeline for the search and the characterization of microlensing flux variations: as a result we detect 3 microlensing candidates. We evaluate the expected signal through a full Monte Carlo simulation of the experiment completed by an analysis of the detection efficiency of our pipeline. We consider both "self lensing" and "MACHO lensing" lens populations, given by M31 stars and dark matter halo MACHOs, in the M31 and the Milky Way (MW), respectively. The total number of events is compatible with the expected self-lensing rate. Specifically, we evaluate an expected signal of about 2 self-lensing events. As for MACHO lensing, for full 0.5 (0.01) solar mass MACHO halos, our prediction is for about 4 (7) events. The comparatively small number of expected MACHO versus self lensing events, together with the small number statistics at disposal, do not enable us to put strong constraints on that population. Rather, the hypothesis, suggested by a previous analysis, on the MACHO nature of OAB-07-N2, one of the microlensing candidates, translates into a sizeable lower limit for the halo mass fraction in form of the would be MACHO population, f, of about 15% for 0.5 solar mass MACHOs.

We develop a novel statistical strong lensing approach to probe the cosmological parameters by exploiting multiple redshift image systems behind galaxies or galaxy clusters. The method relies on free-form mass inversion of strong lenses and does not need any additional information other than gravitational lensing. Since in free-form lensing the solution space is a high-dimensional convex polytope, we consider Bayesian model comparison analysis to infer the cosmological parameters. The volume of the solution space is taken as a tracer of the probability of the underlying cosmological assumption. In contrast to parametric mass inversions, our method accounts for the mass-sheet degeneracy, which implies a degeneracy between the steepness of the profile and the cosmological parameters. Parametric models typically break this degeneracy, introducing hidden priors to the analysis that contaminate the inference of the parameters. We test our method with synthetic lenses, showing that it is able to infer the assumed cosmological parameters. Applied to the CLASH clusters, the method might be competitive with other probes.

Time delays between multiple images of lensed sources can probe the geometry of the universe. We propose a novel method based on free-form modelling of gravitational lenses to estimate time-delay distances and, in turn, cosmological parameters. This approach does not suffer from the degeneracy between the steepness of the profile and the cosmological parameters. We apply the method to 18 systems having time delay measurements and find H_0=69+-6(stat.)+-4(syst.) km s^{-1}Mpc^{-1}. In combination with WMAP9, the constraints on dark energy are Omega_w=0.68+-0.05 and w=-0.86+-0.17 in a flat model with constant equation-of-state.

Mass and concentration of clusters of galaxies are related and evolving with redshift. We study the properties of a sample of 31 massive galaxy clusters at high redshift, 0.8 < z < 1.5, using weak and strong lensing observations. Concentration is a steep function of mass, c_{200} ~M_{200}^{-0.83 +-0.39}, with higher-redshift clusters being less concentrated. Mass and concentration from the stacked analysis, M_{200}=(4.1+-0.4)x10^{14}M_Sun/h and c_{200}=2.3+-0.2, are in line with theoretical results extrapolated from the local universe. Clusters with signs of dynamical activity preferentially feature high concentrations. We discuss the possibility that the whole sample is a mix of two different kinds of haloes. Over-concentrated clusters might be accreting haloes out of equilibrium in a transient phase of compression, whereas less concentrated ones might be more relaxed.

The origin of the accelerated expansion of the universe is still unclear and new physics is needed on cosmological scales. We propose and test a novel interpretation of dark energy as originated by an elastic strain due to a cosmic defect in an otherwise Euclidean space-time. The strain modifies the expansion history of the universe. This new effective contribution tracks radiation at early times and mimics a cosmological constant at late times. The theory is tested against observations, from nucleosynthesis to the cosmic microwave background and formation and evolution of large scale structure to supernovae. Data are very well reproduced with Lam\'e parameters of the order of 10^{-52} m^{-2}.

Knowledge of mass and concentration of galaxy clusters is crucial to understand their formation and evolution. Unbiased estimates require the understanding of the shape and orientation of the halo as well as its equilibrium status. We propose a novel method to determine the intrinsic properties of galaxy clusters from a multi-wavelength data set spanning from X-ray spectroscopic and photometric data to gravitational lensing to the Sunyaev-Zel'dovich effect (SZe). The method relies on two quite non informative geometrical assumptions: the distributions of total matter or gas are approximately ellipsoidal and co-aligned; they have different, constant axial ratios but share the same degree of triaxiality. Weak and strong lensing probe the features of the total mass distribution in the plane of the sky. X-ray data measure size and orientation of the gas in the plane of the sky. Comparison with the SZ amplitude fixes the elongation of the gas along the line of sight. These constraints are deprojected thanks to Bayesian inference. The mass distribution is described as a Navarro-Frenk-White halo with arbitrary orientation, gas density and temperature are modelled with parametric profiles. We applied the method to Abell 1689. Independently of the priors, the cluster is massive, M_{200}=(1.3+-0.2)*10^{15}M_sun, and over-concentrated, c_{200}=8+-1, but still consistent with theoretical predictions. The total matter is triaxial (minor to major axis ratio ~0.5+-0.1 exploiting priors from N-body simulations) with the major axis nearly orientated along the line of sight. The gas is rounder (minor to major axis ratio ~0.6+-0.1) and deviates from hydrostatic equilibrium. The contribution of non-thermal pressure is ~20-50 per cent in inner regions, <~ 300 kpc, and ~25+-5 per cent at ~1.5 Mpc.

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.

The properties of a spherically symmetric static space-time permeated of dark energy are worked out. Dark energy is viewed as the strain energy of an elastically deformable four dimensional manifold. The metric is worked out in the vacuum region around a central spherical mass/defect in the linear approximation. We discuss analogies and differences with the analogue in the de Sitter space time and how these competing scenarios could be differentiated on an observational ground. The comparison with the tests at the solar system scale puts upper limits to the parameters of the theory, consistent with the values obtained applying the classical cosmological tests.

The core collapse supernova (CCSN) rate provides a strong lower limit for the star formation rate (SFR). Progress in using it as a cosmic SFR tracer requires some confidence that it is consistent with more conventional SFR diagnostics in the nearby Universe. This paper compares standard SFR measurements based on Halpha, FUV and TIR galaxy luminosities with the observed CCSN rate in the same galaxy sample. The comparison can be viewed from two perspectives. Firstly, by adopting an estimate of the minimum stellar mass to produce a CCSN one can determine a SFR from SN numbers. Secondly, the radiative SFRs can be assumed to be robust and then the SN statistics provide a constrain on the minimum stellar mass for CCSN progenitors. The novel aspect of this study is that Halpha, FUV and TIR luminosities are now available for a complete galaxy sample within the local 11Mpc volume and the number of discovered SNe in this sample within the last 13 years is high enough to perform a meaningful statistical comparison. We exploit the multi-wavelength dataset from 11HUGS, a volume-limited survey designed to provide a census of SFR in the local Volume. Assuming a lower limit for CCSN progenitors of 8 Msun, the CCSN rate matches the SFR from the FUV luminosity. However the SFR based on Halpha luminosity is lower than these two estimates by a factor of nearly 2. If we assume that the FUV or Halpha based luminosities are a true reflection of the SFR, we find that the minimum mass for CCSN progenitors is 8 +/- 1 Msun, and 6 +/- 1 Msun, respectively. The estimate of the minimum mass for CCSN progenitors obtained exploiting FUV data is in good agreement with that from the direct detection of CCSN progenitors. The concordant results by these independent methods point toward a constraint of 8 +/- 1 Msun on the lower mass limit for progenitor stars of CCSNe.

The high concentrations derived for several strong-lensing clusters present a major inconsistency between theoretical LambdaCDM expectations and measurements. Triaxiality and orientation biases might be at the origin of this disagreement, as clusters elongated along the line-of-sight would have a relatively higher projected mass density, boosting the resulting lensing properties. Analyses of statistical samples can probe further these effects and crucially reduce biases. In this work we perform a fully triaxial strong-lensing analysis of the 12 MACS clusters at z > 0.5, a complete X-ray selected sample, and fully account for the impact of the intrinsic 3D shapes on their strong lensing properties. We first construct strong-lensing mass models for each cluster based on multiple-images, and fit projected ellipsoidal Navarro-Frenk-White halos with arbitrary orientations to each mass distribution. We then invert the measured surface mass densities using Bayesian statistics. Although the Einstein radii of this sample are significantly larger than predicted by LambdaCDM, here we find that the mass-concentration relation is in full agreement with results from N-body simulations. The z > 0.5 MACS clusters suffer from a moderate form of orientation bias as may be expected for X-ray selected samples. Being mostly unrelaxed, at a relatively high redshift, with high X-ray luminosity and noticeable substructures, these clusters may lie outside the standard concentration-Einstein radius relation. Our results remark the importance of triaxiality and properly selected samples for understanding galaxy clusters properties, and suggest that higher-z, unrelaxed low-concentration clusters form a different class of prominent strong gravitational lenses. Arc redshift confirmation and weak lensing data in the outer region are needed to further refine our analysis.

Knowledge of intrinsic shape and orientation of galaxy clusters is crucial to understand their formation and evolution. We propose a novel model which uses Bayesian inference to determine the intrinsic form of the hot intracluster medium of galaxy clusters. The method exploits X-ray spectroscopic and photometric data plus measurements of the Sunyaev-Zel'dovich effect (SZe). The gas distribution is modelled with an ellipsoidal parametric profile who can fit observed X-ray surface-brightness and temperature. Comparison with the SZ amplitude fixes the elongation along the line of sight. Finally, Bayesian inference allows us to deproject the measured elongation and the projected ellipticity and constrain the intrinsic shape and orientation of the cluster. We apply the method to the rich cluster Abell 1689, which was targeted by the Chandra and XMM satellites as well as by several SZe observatories. Observations cover in detail a region <~ 1 Mpc. Our analysis favours a mildly triaxial cluster with a minor to major axis ratio of 0.70+-0.15, preferentially elongated along the line of sight, as expected for massive lensing clusters. The triaxial structure together with the orientation bias can reconcile X-ray with lensing analyses and supports the view of A1689 as a just slightly over-concentrated massive cluster not so far from hydrostatic equilibrium.

Halos formed in the standard Lambda cold dark matter framework should follow an universal mass density profile and fit a well defined mass-concentration relation. Lensing analyses of clusters with a large Einstein radius seem to contradict this scenario, with the massive cluster Abell 1689 being often claimed as a notable example of a highly over-concentrated halo. Shape and orientation biases in lensing studies might be at the basis of this disagreement between theory and observations. We developed a method for a full three-dimensional analysis of strong and weak lensing data. Surface density maps estimated from lensing are de-projected to infer the actual triaxial structure of the cluster, whose mass distribution is approximated as an ellipsoidal Navarro-Frenk-White halo with arbitrary orientation. Inversion is performed under competing a priori assumptions, integrated in the method thanks to Bayesian statistics. We applied the method to A1689. Whatever the considered priors on shape and orientation, both weak and strong lensing analyses found the halo to be slightly over-concentrated but still consistent with theoretical predictions. We found some evidence for a mildly triaxial lens (minor to major axis ratio ~ 0.5 +- 0.2) with the major axis orientated along the line of sight. Exploiting priors from N-body simulations, we found mass M_{200} = (1.3 +- 0.4)x 10^{15}M_Sun and concentration c_{200} =10+-3 for the weak lensing analysis of Subaru data, M_{200} = (1.7 +- 0.3)x 10^{15}M_Sun and c_{200} =6.1+-0.9 for the strong lensing analysis of multiple image systems, and M_{200} = (1.3 +- 0.2)\times 10^{15}M_Sun and c_{200} =7.3+-0.8 for the combined weak plus strong analysis.

LISA might detect gravitational waves from mergers of massive black hole binaries strongly lensed by intervening galaxies (Sereno et al. 2010). The detection of multiple gravitational lensing events would provide a new tool for cosmography. Constraints on cosmological parameters could be placed by exploiting either lensing statistics of strongly lensed sources or time delay measurements of lensed gravitational wave signals. These lensing methods do not need the measurement of the redshifts of the sources and the identification of their electromagnetic counterparts. They would extend cosmological probes to redshift z <= 10 and are then complementary to other lower or higher redshift tests, such as type Ia supernovae or cosmic microwave background. The accuracy of lensing tests strongly depends on the formation history of the merging binaries, and the related number of total detectable multiple images. Lensing amplification might also help to find the host galaxies. Any measurement of the source redshifts would allow to exploit the distance-redshift test in combination with lensing methods. Time-delay analyses might measure the Hubble parameter H_0 with accuracy of >= 10 km s^{-1}Mpc^{-1}. With prior knowledge of H_0, lensing statistics and time delays might constrain the dark matter density (delta Omega_M >= 0.08, due to parameter degeneracy). Inclusion of our methods with other available orthogonal techniques might significantly reduce the uncertainty contours for Omega_M and the dark energy equation of state.

The Cosmic Defect theory has been confronted with four observational constraints: primordial nuclear species abundances emerging from the big bang nucleosynthesis; large scale structure formation in the universe; cosmic microwave background acoustic scale; luminosity distances of type Ia supernovae. The test has been based on a statistical analysis of the a posteriori probabilities for three parameters of the theory. The result has been quite satisfactory and such that the performance of the theory is not distinguishable from the one of the Lambda-CDM theory. The use of the optimal values of the parameters for the calculation of the Hubble constant and the age of the universe confirms the compatibility of the Cosmic Defect approach with observations.

We discuss strong gravitational lensing of gravitational waves from merging of massive black hole binaries in the context of the LISA mission. Detection of multiple events would provide invaluable information on competing theories of gravity, evolution and formation of structures and, with complementary observations, constraints on H_0 and other cosmological parameters. Most of the optical depth for lensing is provided by intervening massive galactic halos, for which wave optics effects are negligible. Probabilities to observe multiple events are sizable for a broad range of formation histories. For the most optimistic models, up to 4 multiple events with a signal to noise ratio >= 8 are expected in a 5-year mission. Chances are significant even for conservative models with either light (<= 60%) or heavy (<= 40%) seeds. Due to lensing amplification, some intrinsically too faint signals are brought over threshold (<= 2 per year).

We present an updated analysis of the M31 pixel lensing candidate event OAB-N2 previously reported in Calchi Novati et al. (2009). Here we take advantage of new data both astrometric and photometric. Astrometry: using archival 4m-KPNO and HST/WFPC2 data we perform a detailed analysis on the event source whose result, although not fully conclusive on the source magnitude determination, is confirmed by the following light curve photometry analysis. Photometry: first, unpublished WeCAPP data allows us to confirm OAB-N2, previously reported only as a viable candidate, as a well constrained pixel lensing event. Second, this photometry enables a detailed analysis in the event parameter space including the effects due to finite source size. The combined results of these analyses allow us to put a strong lower limit on the lens proper motion. This outcome favors the MACHO lensing hypothesis over self lensing for this individual event and points the way toward distinguishing between the MACHO and self-lensing hypotheses from larger data sets.

Strong lensing studies can provide detailed mass maps of the inner regions even in dynamically active galaxy clusters. It is shown that proper modelling of the intracluster medium, i.e. the main baryonic component, can play an important role. In fact, the addition of a new contribution accounting for the gas can increase the statistical significance of the lensing model. We propose a parametric method for strong lensing analyses which exploits multi-wavelength observations. The mass model accounts for cluster-sized dark matter halos, galaxies (whose stellar mass can be obtained from optical analyses) and the intracluster medium. The gas distribution is fitted to lensing data exploiting prior knowledge from X-ray observations. This gives an unbiased look at each matter component and allows us to study the dynamical status of a cluster.}The method has been applied to AC 114, an irregular X-ray cluster. We find positive evidence for dynamical activity, with the dark matter distribution shifted and rotated with respect to the gas. On the other hand, the dark matter follows the galaxy density both for shape and orientation, which hints at its collisionless nature. The inner region (< 250kpc) is under-luminous in optical bands whereas the gas fraction (~20+- 5%) slightly exceeds typical values. Evidence from lensing and X-ray suggests that the cluster develops in the plane of the sky and is not affected by the lensing over-concentration bias. Despite the dynamical activity, the matter distribution seems to be in agreement with predictions from N-body simulations. An universal cusped profile provides a good description of either the overall or the dark matter distribution whereas theoretical scaling relations seem to be nicely fitted.

The analysis of non-radiative sources of static or time-dependent gravitational fields in the Solar System is crucial to accurately estimate the free-fall orbits of the LISA space mission. In particular, we take into account the gravitational effects of Interplanetary Dust (ID) on the spacecraft trajectories. The perturbing gravitational field has been calculated for some ID density distributions that fit the observed zodiacal light. Then we integrated the Gauss planetary equations to get the deviations from the LISA keplerian orbits around the Sun. This analysis can be eventually extended to Local Dark Matter (LDM), as gravitational fields are expected to be similar for ID and LDM distributions. Under some strong assumptions on the displacement noise at very low frequency, the Doppler data collected during the whole LISA mission could provide upper limits on ID and LDM densities.

Lambda cold dark matter paradigm predicts that galaxy clusters follow an universal mass density profile and fit a well defined mass-concentration relation, with lensing clusters being preferentially triaxial haloes elongated along the line of sight. Oddly, recent strong and weak lensing analyses of clusters with a large Einstein radius suggested those haloes to be highly over-concentrated. Here, we investigate what intrinsic shape and orientation an halo should have to account for both theoretical predictions and observations. We considered a sample of 10 strong lensing clusters. We first measured their elongation assuming a given mass-concentration relation. Then, for each cluster we found the intrinsic shape and orientation which are compatible with the inferred elongation and the measured projected ellipticity. We distinguished two groups. The first one (nearly one half) seems to be composed of outliers of the mass-concentration relation, which they would fit only if they were characterised by a filamentary structure extremely elongated along the line of sight, that is not plausible considering standard scenarios of structure formations. The second sample supports expectations of N-body simulations which prefer mildly triaxial lensing clusters with a strong orientation bias.

With one exception, previous analyses of the measurement accuracy of gravitational wave experiments for comparable-mass binary systems have neglected either spin-precession effects or subdominant harmonics and amplitude modulations. Here we give the first explicit description of how these effects combine to improve parameter estimation. We consider supermassive black hole binaries as expected to be observed with the planned space-based interferometer LISA, and study the measurement accuracy for several astrophysically interesting parameters obtainable taking into account the full 2PN waveform for spinning bodies, as well as spin-precession effects. We find that for binaries with a total mass in the range 10^5 M_Sun < M < 10^7 M_Sun at a redshift of 1, a factor ~1.5 is in general gained in accuracy, with the notable exception of the determination of the individual masses in equal-mass systems, for which a factor ~5 can be gained. We also find, as could be expected, that using the full waveform helps increasing the upper mass limit for detection, which can be as high as M = 10^8 M_Sun at a redshift of 1, as well as the redshift limit where some information can be extracted from a system, which is roughly z = 10 for M < 10^7 M_Sun, 1.5-5 times higher than with the restricted waveform. We computed that the full waveform allows to use supermassive black hole binaries as standard sirens up to a redshift of z = 1.6, about 0.4 larger than what previous studies allowed. We found that for lower unequal-mass binary systems, the measurement accuracy is not as drastically improved as for other systems. This suggests that for these systems, adding parameters such as eccentricity or alternative gravity parameters could be achieved without much loss in the accuracy.

We present comprehensive photometric and spectroscopic observations of the faint transient SN 2008S discovered in NGC 6946. SN 2008S exhibited slow photometric evolution and almost no spectral variability during the first nine months, implying a high density CS medium. The light curve is similar in shape to that of SN 1998S and SN 1979C, although significantly fainter at maximum light. Our quasi-bolometric lightcurve extends to 300 days and shows a tail phase decay rate consistent with that of ^{56}Co. We propose that this is evidence for an explosion and formation of ^{56}Ni (0.0015 +/- 0.0004 M_Sun). The large MIR flux detected shortly after explosion can be explained by a light echo from pre-exisiting dust. The late NIR flux excess is plausibly due to a combination of warm newly-formed ejecta dust together with shock-heated dust in the CS environment. We reassess the progenitor object detected previously in Spitzer archive images, supplementing this discussion with a model of the MIR spectral energy distribution. This supports the idea of a dusty, optically thick shell around SN 2008S with an inner radius of nearly 90AU and outer radius of 450AU, and an inferred heating source of 3000 K and luminosity of L ~ 10^{4.6} L_Sun. The combination of our monitoring data and the evidence from the progenitor analysis leads us to support the scenario of a weak electron capture supernova explosion in a super-AGB progenitor star (of initial mass 6-8 M_sun) embedded within a thick CS gaseous envelope. We suggest that all of main properties of the electron capture SN phenomenon are observed in SN 2008S and future observations may allow a definitive answer.

Microlensing observations towards M31 are a powerful tool for the study of the dark matter population in the form of MACHOs both in the Galaxy and the M31 halos, a still unresolved issue, as well as for the analysis of the characteristics of the M31 luminous populations. In this work we present the second year results of our pixel lensing campaign carried out towards M31 using the 152 cm Cassini telescope in Loiano. We have established an automatic pipeline for the detection and the characterisation of microlensing variations. We have carried out a complete simulation of the experiment and evaluated the expected signal, including an analysis of the efficiency of our pipeline. As a result, we select 1-2 candidate microlensing events (according to different selection criteria). This output is in agreement with the expected rate of M31 self-lensing events. However, the statistics are still too low to draw definitive conclusions on MACHO lensing.

The drag-free satellites of LISA will maintain the test masses in geodesic motion over many years with residual accelerations at unprecedented small levels and time delay interferometry (TDI) will keep track of their differential positions at level of picometers. This may allow investigations of fine details of the gravitational field in the Solar System previously inaccessible. In this spirit, we present the concept of a method to measure directly the gravitational effect of the density of diffuse Local Dark Matter (LDM) with a constellation of a few drag-free satellites, by exploiting how peculiarly it would affect their relative motion. Using as test bed an idealized LISA with rigid arms, we find that the separation in time between the test masses is uniquely perturbed by the LDM, so that they acquire a differential breathing mode. Such a LDM signal is related to the LDM density within the orbits and has characteristic spectral components, with amplitudes increasing in time, at various frequencies of the dynamics of the constellation. This is the relevant result, in that the LDM signal is brought to non-zero frequencies.

We study the effect of a cosmological constant $\Lambda$ on the propagation and detection of gravitational waves. To this purpose we investigate the linearised Einstein's equations with terms up to linear order in $\Lambda$ in a de Sitter and an anti-de Sitter background spacetime. In this framework the cosmological term does not induce changes in the polarization states of the waves, whereas the amplitude gets modified with terms depending on $\Lambda$. Moreover, if a source emits a periodic waveform, its periodicity as measured by a distant observer gets modified. These effects are, however, extremely tiny and thus well below the detectability by some twenty orders of magnitude within present gravitational wave detectors such as LIGO or future planned ones such as LISA.

The cosmological constant Lambda affects cosmological gravitational lensing. Effects due to Lambda can be studied in the framework of the Schwarzschild-de Sitter spacetime. Two novel contributions, which can not be accounted for by a proper use of angular diameter distances, are derived. First, a term 2m b Lambda/3 has to be added to the bending angle, where "m" is the lens mass and "b" the impact parameter. Second, Lambda brings about a difference in the redshifts of multiple images. Both effects are quite small for real astrophysical systems (contribution to the bending < 0.1 microarcsec and difference in redshift < 10^{-7}).

It is known that a relative translational motion between the deflector and the observer affects gravitational lensing. In this paper, a lens equation is obtained to describe such effects on actual lensing observables. Results can be easily interpreted in terms of aberration of light-rays. Both radial and transverse motions with relativistic velocities are considered. The lens equation is derived by first considering geodesic motion of photons in the rest-frame Schwarzschild spacetime of the lens, and, then, light-ray detection in the moving observer's frame. Due to the transverse motion images are displaced and distorted in the observer's celestial sphere, whereas the radial velocity along the line of sight causes an effective re-scaling of the lens mass. The Einstein ring is distorted to an ellipse whereas the caustics in the source plane are still point-like. Either for null transverse motion or up to linear order in velocities, the critical curve is still a circle with its radius corrected by a factor (1+z_d) with respect to the static case, z_d being the relativistic Doppler shift of the deflector. From the observational point of view, the orbital motion of the Earth can cause potentially observable corrections of the order of the microarcsec in lensing towards the super-massive black hole at the Galactic center. On a cosmological scale, tangential peculiar velocities of cluster of galaxies bring about a typical flexion in images of background galaxies in the weak lensing regime but future measurements seem to be too much challenging.

The primary optical caustic surface behind a Kerr black hole is a four-cusped tube displaced from the line of sight. We derive the caustic surface in the nearly asymptotic region far from the black hole through a Taylor expansion of the lightlike geodesics up to and including fourth-order terms in m/b and a/b, where $m$ is the black hole mass, a the spin and b the impact parameter. The corresponding critical locus in the observer's sky is elliptical and a point-like source inside the caustics will be imaged as an Einstein cross. With regard to lensing near critical points, a Kerr lens is analogous to a circular lens perturbed by a dipole and a quadrupole potential. The caustic structure of the supermassive black hole in the Galactic center could be probed by lensing of low mass X-ray binaries in the Galactic inner regions or by hot spots in the accretion disk.

The self-accelerating braneworld model (DGP) seems to provide a simple alternative to the the standard $\Lambda$CDM cosmology to explain the current cosmic acceleration, which is strongly indicated by measurements of Type Ia supernovae, as well as other concordant observations. In this work, we investigate observational constraints on this scenario from gravitational lensing statistics using the Cosmic Lens All-Sky Survey (CLASS) lensing sample. We show that a large parameter space of the DGP model is in good agreement with this radio source gravitational lensing sample. In the flat case, $\Omega_\mathrm{K}=0$, the likelihood is maximized, ${\cal L}={\cal L_\mathrm{max}}$, for $\1 = 0.30_{-0.11}^{+0.19}$. If we relax the prior on $\Omega_\mathrm{K}$, the likelihood peaks at $\{\1,\2 \} \simeq \{0.29, 0.12\}$, just slightly in the region of open models. However the confidence contours are pretty elongated so that we can not discard either close or flat or open models

The cosmological constant Lambda affects gravitational lensing phenomena. The contribution of Lambda to the observable angular positions of multiple images and to their amplification and time delay is here computed through a study in the weak deflection limit of the equations of motion in the Schwarzschild-de Sitter metric. Due to Lambda the unresolved images are slightly demagnified, the radius of the Einstein ring decreases and the time delay increases. The effect is however negligible for near lenses. In the case of null cosmological constant, we provide some updated results on lensing by a Schwarzschild black hole.

We perform a combined X-ray and strong lensing analysis of RX J1347.5-1145, one of the most luminous galaxy clusters at X-ray wavelengths. We show that evidence from strong lensing alone, based on published VLT and new HST data, strongly argues in favor of a complex structure. The analysis takes into account arc positions, shapes and orientations and is done thoroughly in the image plane. The cluster inner regions are well fitted by a bimodal mass distribution, with a total projected mass of $M_{tot} = (9.9 \pm 0.3)\times 10^{14} M_\odot/h$ within a radius of $360 \mathrm{kpc}/h$ ($1.5'$). Such a complex structure could be a signature of a recent major merger as further supported by X-ray data. A temperature map of the cluster, based on deep Chandra observations, reveals a hot front located between the first main component and an X-ray emitting South Eastern sub-clump. The map also unveils a filament of cold gas in the innermost regions of the cluster, most probably a cooling wake caused by the motion of the cD inside the cool core region. A merger scenario in the plane of the sky between two dark matter sub-clumps is consistent with both our lensing and X-ray analyses, and can explain previous discrepancies with mass estimates based on the virial theorem.

We discuss the influence of the cosmological constant $\Lambda$ on the gravitational equations of motion of bodies with arbitrary masses and eventually solve the two-body problem. Observational constraints are derived from measurements of the periastron advance in stellar systems, in particular binary pulsars and the solar system. For the latter we consider also the change in the mean motion due to $\Lambda$. Up to now, Earth and Mars data give the best constraint, $\Lambda \sim 10^{-36} \mathrm{km}^{-2}$. If properly accounting for the gravito-magnetic effect, this upper limit on $\Lambda$ could greatly improve in the near future thanks to new data from planned or already operating space-missions. Dark matter or modifications of the Newtonian inverse-square law in the solar system are discussed as well. Variations in the $1/r^2$ behavior are considered in the form of either a possible Yukawa-like interaction or a modification of gravity of MOND type.

Knowledge of the intrinsic shape of galaxy clusters is very important in investigating cosmic structure formation and astrophysical processes. The reconstruction of the 3-dimensional structure usually relies on deprojecting 2-dimensional X-ray, Sunyaev-Zeldovich (SZ) and/or gravitational lensing observations. As known, a joint analysis of these data sets can provide the elongation of the cluster along the line of sight together with its length and width in the plane of the sky. An unbiased measurement of the Hubble constant can be also inferred. Due to some intrinsic degeneracies, the observational constraints obtained from such projected data-sets are not enough to allow an unique inversion. In general, the projected maps can be at the same time compatible with prolate, oblate and with many triaxial configurations. Even a prolate cluster might be interpreted as an oblate system and vice versa. Assuming that the cluster is axially symmetric is likely to overestimate the intrinsic ellipticity, whereas the system always looks rounder performing the inversion under the hypothesis of a triaxial cluster aligned with the line of sight. In general, analysing triaxial clusters under the prolate or oblate assumption may introduce strong biases even when the clusters are actually near to axial symmetry whereas the systematics introduced assuming the cluster to be aligned with the line of sight are more under control.

Gravitational lensing provides an efficient tool for the investigation of matter structures, independent of the dynamical or hydrostatic equilibrium properties of the deflecting system. However, it depends on the kinematic status. In fact, either a translational motion or a coherent rotation of the mass distribution can affect the lensing properties. Here, light deflection by galaxy clusters in motion is considered. Even if gravitational lensing mass measurements of galaxy clusters are regarded as very reliable estimates, the kinematic effect should be considered. A typical peculiar motion with respect to the Hubble flow brings about a systematic error < 0.3%, independent of the mass of the cluster. On the other hand, the effect of the spin increases with the total mass. For cluster masses ~ 10^{15}M_{sun}, the effect of the gravitomagnetic term is < 0.04% on strong lensing estimates and < 0.5% in the weak lensing analyses. The total kinematic effect on the mass estimate is then < 1%, which is negligible in current statistical studies. In the weak lensing regime, the rotation imprints a typical angular modulation in the tangential shear distortion. This would allow in principle a detection of the gravitomagnetic field and a direct measurement of the angular velocity of the cluster but the required background source densities are well beyond current tecnological capabilities.

The gravitational action of the smooth energy-matter components filling in the universe can affect the orbit of a planetary system. Changes are related to the acceleration of the cosmological scale size R. In a universe with significant dark matter, a gravitational system expands or contracts according to the amount and equation of state of the dark energy. At present time, the Solar system, according to the LambdaCDM scenario emerging from observational cosmology, should be expanding if we consider only the effect of the cosmological background. Its fate is determined by the equation of state of the dark energy alone. The mean motion and periastron precession of a planet are directly sensitive to (d^2 R/d t^2)/R, whereas variations with time in the semi-major axis and eccentricity are related to its time variation. Actual bounds on the cosmological deceleration parameters q_0 from accurate astrometric data of perihelion precession and changes in the third Kepler's law in the Solar system fall short of ten orders of magnitude with respect to estimates from observational cosmology. Future radio-ranging measurements of outer planets could improve actual bounds by five orders of magnitude.

The shear triplet statistics is a geometric method to measure cosmological parameters with observations in the weak gravitational lensing regime towards massive haloes. Here, this proposal is considered to probe the dark energy equation of state and its time derivative in view of future wide-field galaxy surveys. A survey with a median redshift of nearly 0.7 and a total area of nearly 10000 square degrees would be pretty effective in determining the dark matter cosmological density and in putting useful constraints on the dark energy properties.

We present an analytical treatment of gravitational lensing by a Kerr black hole in the weak deflection limit. Lightlike geodesics are expanded as a Taylor series up to and including third-order terms in m/b and a/b, where m is the black hole mass, a the angular momentum and b the impact parameter of the light ray. Positions and magnifications of individual images are computed with a perturbative analysis. At this order, the degeneracy with the translated Schwarzschild lens is broken. The critical curve is still a circle displaced from the black hole position in the equatorial direction and the corresponding caustic is point-like. The degeneracy between the black hole spin and its inclination relative to the observer is broken through the angular coordinates of the perturbed images.

Dark matter or modifications of the Newtonian inverse-square law in the solar-system are studied with accurate planetary astrometric data. From extra-perihelion precession and possible changes in the third Kepler's law, we get an upper limit on the local dark matter density, rho_{DM} < 3*10^{-16} kg/m^3 at the 2-sigma confidence level. Variations in the 1/r^2 behavior are considered in the form of either a possible Yukawa-like interaction or a modification of gravity of MOND type. Up to scales of 10^{11} m, scale-dependent deviations in the gravitational acceleration are really small. We examined the MOND interpolating function mu in the regime of strong gravity. Gradually varying mu suggested by fits of rotation curves are excluded, whereas the standard form mu(x)= x/(1+x^2)^{1/2} is still compatible with data. In combination with constraints from galactic rotation curves and theoretical considerations on the external field effect, the absence of any significant deviation from inverse square attraction in the solar system makes the range of acceptable interpolating functions significantly narrow. Future radio ranging observations of outer planets with an accuracy of few tenths of a meter could either give positive evidence of dark matter or disprove modifications of gravity.

Deviations from the gravitational inverse-square law would imprint scale-dependent features on the power spectrum of mass density fluctuations. We model such deviations as a Yukawa-like contribution to the gravitational potential and discuss the growth function in a mixed dark matter model with adiabatic initial conditions. Evolution of perturbations is considered in general non-flat cosmological models with a cosmological constant, and an analytical approximation for the growth function is provided. The coupling between baryons and cold dark matter across recombination is negligibly affected by modified gravity physics if the proper cutoff length of the long-range Yukawa-like force is > 10 h^{-1} Mpc. Enhancement of gravity affects the subsequent evolution, boosting large-scale power in a way that resembles the effect of a lower matter density. This phenomenon is almost perfectly degenerate in power-spectrum shape with the effect of a background of massive neutrinos. Back-reaction on density growth from a modified cosmic expansion rate should however also affect the normalization of the power spectrum, with a shape distortion similar to the case of a non-modified background.

We investigate the strong deflection limit of gravitational lensing by a Schwarzschild black hole embedded in an external gravitational field. The study of this model, analogous to the Chang & Refsdal lens in the weak deflection limit, is important to evaluate the gravitational perturbations on the relativistic images that appear in proximity of supermassive black holes hosted in galactic centers. By a simple dimensional argument, we prove that the tidal effect on the light ray propagation mainly occurs in the weak field region far away from the black hole and that the external perturbation can be treated as a weak field quadrupole term. We provide a description of relativistic critical curves and caustics and discuss the inversion of the lens mapping. Relativistic caustics are shifted and acquire a finite diamond shape. Sources inside the caustics produce four sequences of relativistic images. On the other hand, retro-lensing caustics are only shifted while remaining point-like to the lowest order.

Some tests of gravity theories - periastron shift, geodetic precession, change in mean motion and gravitational redshift - are applied in solar and stellar systems to constrain the cosmological constant. We thus consider a length scale range from 10^8 to 10^{15} km. Best bounds from the solar system come from perihelion advance and change in mean motion of Earth and Mars, Lambda < 10^{-36} km^{-2}. Such a limit falls very short to estimates from observational cosmology analyses but a future experiment performing radio ranging observations of outer planets could improve it by four orders of magnitude. Beyond the solar system, together with future measurements of periastron advance in wide binary pulsars, gravitational redshift of white dwarfs can provide bounds competitive with Mars data.

The intrinsic shape of galaxy clusters can be obtained through a combination of X-ray and Sunyaev-Zeldovich effect observations once cosmological parameters are assumed to be known. In this paper we discuss the feasibility of modelling galaxy clusters as either prolate or oblate ellipsoids. We analyze the intra-cluster medium distribution for a sample of 25 X-ray selected clusters, with measured Sunyaev-Zeldovich temperature decrements. A mixed population of prolate and oblate ellipsoids of revolution fits the data well, with prolate shapes preferred on a 60-76% basis. We observe an excess of clusters nearly aligned along the line of sight, with respect to what is expected from a randomly oriented cluster population, which might imply the presence of a selection bias in our sample. We also find signs that a more general triaxial morphology might better describe the morphology of galaxy clusters. Additional constraints from gravitational lensing could disentangle the degeneracy between an ellipsoidal and a triaxial morphology, and could also allow an unbiased determination of the Hubble constant.

We discuss the influence of the cosmological constant on the gravitational equations of motion of bodies with arbitrary masses and eventually solve the two-body problem. Observational constraints are derived from measurements of the periastron advance in stellar systems, in particular binary pulsars and the solar system. Up to now, Earth and Mars data give the best constraint, Lambda < 10^{-36} km^{-2}; bounds from binary pulsars are potentially competitive with limits from interplanetary measurements. If properly accounting for the gravito-magnetic effect, this upper limit on $\Lambda$ could greatly improve in the near future thanks to new data from planned or already operating space-missions.

In this paper we present an analytical treatment of gravitational lensing by Kerr black holes in the limit of very large deflection angles, restricting to observers in the equatorial plane. We accomplish our objective starting from the Schwarzschild black hole and adding corrections up to second order in the black hole spin. This is sufficient to provide a full description of all caustics and the inversion of lens mapping for sources near them. On the basis of these formulae we argue that relativistic images of Low Mass X-ray Binaries around Sgr A* are very likely to be seen by future X-ray interferometry missions.

We here combine Sunyaev-Zel'dovich effect, X-ray observations and spectroscopic redshifts of member galaxies, to constrain the intrinsic three-dimensional shape of the galaxy cluster: Abell 370. The cluster turns out to be strongly elongated along the l.o.s., with two (or more) substructures in the process of merging. Spectroscopy further suggests that the process must be taking place at a small angle respect to the l.o.s.

We discuss a method to constrain the intrinsic three-dimensionale shapes of galaxy clusters by combining X-Ray and Sunyaev-Zeldovich observations. The method is applied to a sample of 25 X-Ray selected clusters, with measured Sunyaev-Zeldovich temperature decrements. The sample turns out to be slightly biased, with strongly elongated clusters preferentially aligned along the line of sight. This result demonstrates that X-Ray selected cluster samples may be affected by morphological and orientation effects even if a relatively high threshold signal-to-noise ratio is used to select the sample. A large majority of the clusters in our sample exhibit a marked triaxial structure, with prolate-like shapes being slightly more likely than oblate-like ones; the spherical hypothesis is strongly rejected for most sample members. Cooling flow clusters do not show preferentially regular morphologies.

Gravitational lensing by a spinning deflector in translational motion relative to the observer is discussed in the weak field, slow motion approximation. The effect of rotation, which generates an intrinsic gravito-magnetic field, separates from that due to radial motion. Corrections to the lens equation, deflection angle and time delay are derived.

We use the adiabatic compression theory to build a physically well - motivated Milky Way mass model in agreement with the observational data. The visible mass of the Galaxy is distributed in a spheroidal bulge and a multi - components disc parametrized by three galactic parameters, the Sun distance to the galactic centre, R_0, the total bulge mass, M_{bulge}, and the local disc surface density, \Sigma_{\odot}. To model the dark matter component, we adiabatically compress a Navarro, Frenk and White (NFW) halo (with concentration $c$ and total mass M_{vir}) for fixed values of the spin parameter, \lambda, the fraction of the mass in baryons, $m_b$, and the thin disc contribution to total angular momentum, j_d. An iterative selection procedure is used to explore in very detail the wide space of parameters only selecting those combinations of {R_0, M_{bulge}, \Sigma_{\odot}, \lambda, m_b, j_b, c, M_{vir}} that give rise to a Milky Way model in agreement with the observational constraints. This analysis leads us to conclude that only models with R_0 = 8.5 kpc, $0.8 x 10^{10} M_{\odot} < M_{bulge} < 1.6 x 10^{10} M_{\odot}$ and $49 M_{\odot} pc^{-2} \le \Sigma_{\odot} \le 56 M_{\odot} pc^{-2}$ can be reconciled with the set of observational constraints. As regard the parameters entering the adiabatic compression, we find $0.03 \le \lambda \le 0.10$ and $0.04 \le m_b \le 0.10$, while the final estimates of the parameters describing the initial halo profile turn out to be $5 \stackrel{<}{\sim} c \stackrel{<}{\sim} 12$ and $7 x 10^{11} M_{\odot} \stackrel{<}{\sim} M_\mathrm{vir} \stackrel{<}{\sim} 17 x 10^{11} M_{\odot}$ (all at 95.7% CL).

We evaluate the effect of small scale inhomogeneities on large scale observations within the statistics of gravitationally lensed quasars. At this aim, we consider a cosmological model whose large scale properties (dynamics, matter distribution) are the same as in Friedmann-Lemaitre models, but whose matter distribution is locally inhomogeneous. We use the well known Dyer-Roder distances to allow a simple analytical expression of the optical depth, and pay particular attention on the different role played by the different notions of distances (filled beam angular diameter distance and Dyer-Roder distances) when calculating this quantity, following the prescription from Ehlers & Schneider for a coherent formalism. We find that the expected number of gravitationally lensed quasars is a decreasing function of the clumpiness parameter.

Gravitational lensing by rotating extended deflectors is discussed. Due to spin, corrections on image positions, caustics and critical curves can be significant. In order to obtain realistic quantitative estimates, the lens is modeled as a singular isothermal sphere. Gravito-magnetic time delays of nearly 0.2 days between different images of background sources can occur.

Constraints on an exact quintessence scalar-field model with an exponential potential are derived from gravitational lens statistics. An exponential potential can account for data from both optical quasar surveys and radio selected sources. Based on the Cosmic Lens All-Sky Survey (CLASS) sample, lensing statistics provides, for the pressureless matter density parameter, an estimate of Omega_{M0} = 0.31(-0.14)(+0.12).

We discuss the effects of an isolated gravitational lens on the rotation of the plane of polarization of linearly polarized light rays, the so called gravitational Faraday rotation, in metric theories of gravity. By applying the thin lens approximation, we derive simple expressions for extended, rotating astrophysical systems deflecting electromagnetic radiation from background sources. Higher order corrections and the case of a Kerr black hole are also considered to show how the rotation of polarization is a purely gravito-magnetic effect. Prospects of a future detection of a gravitational Faraday rotation are examined.

We discuss a new method which is potentially capable of constraining cosmological parameters using observations of giant luminous arcs in rich X-ray clusters of galaxies. The mass profile and the mass normalization of the lenses are determined from X-ray measurements. The method also allows to probe the amount and equation of state of the dark energy in the universe. The analysis of a preliminary sample of 6 luminous, relatively relaxed clusters of galaxies strongly favours an accelerating expansion of the universe. Under the assumption that the dark energy is in the form of a cosmological constant, the data provide an estimate of $ \Omega_\Lambda = 1.1$ with a statistical error of $\pm 0.2$. Including the constraint of a flat universe and an equation of state for the dark energy $w_X \geq -1$, we obtain $\Omega_M = 0.10 \pm 0.10$ and $w_X = -0.84 \pm 0.14$. Relaxing the prior on $w_X$, we find that the null energy condition is satisfied at the 3-$\sigma$ confidence level.

We present Chandra observations of the galaxy cluster AC114, which shows a strongly irregular morphology, with signs of multiple merging activity. We report the discovery of a soft X-ray filament originating close to the core of the cluster. We confirm that X-ray emission is associated with two of three mass concentrations identified in previous gravitational lensing studies of this object. These two mass concentrations are located at opposite ends of the soft filament, evidence for interaction between them. In the northern part, the cluster shows two sharp discontinuities, both in surface brightness and in temperature, evincing another, more recent merger event which took place in, or close, to the cluster core. In spite of the merger activity, a combined mass and lensing analysis shows remarkably good agreement between lensing and X-ray masses. We therefore advocate for the robustness of the X-ray mass estimates, and conclude that the assumption of hydrostatic equilibrium can yield accurate mass estimates even in clusters as dynamically active as AC 114, once the gas density distribution is properly mapped.

We examine the rotation of the plane of polarization for linearly polarized light rays by the weak gravitational field of an isolated physical system. Based on the rotation of inertial frames, we review the general integral expression for the net rotation. We apply this formula, analogue to the usual electromagnetic Faraday effect, to some interesting astrophysical systems: uniformly shifting mass monopoles and a spinning external shell.

We study gravitational lensing by a Reissner-Nordstrom (RN) black hole in the weak field limit. We obtain the basic equations for the deflection angle and time delay and find analytical expressions for the positions and amplifications of the primary and secondary images. Due to a net positive charge, the separation between images increases, but no change in the total magnification occurs.

Clusters of galaxies contain a fair sample of the universal baryonic mass fraction. A combined analysis of the intracluster medium (ICM) within their hydrostatic regions, as derived from both Sunyaev-Zeldovich effect (SZE) measurements and X-ray images, makes it possible to constrain the cosmological parameters. We consider both gas fraction estimates and angular diameter distance measurements. Adopting median statistics, we find, at the 2-$\sigma$ level, the pressureless matter density, $\Omega_{\rm M0}$, to be between 0.30 and 0.40 and the Hubble constant, $H_0$, between 44 and 66 Km s$^{-1}$ Mpc$^{-1}$.

Gravitational lensing by spinning stars, approximated as homogeneous spheres, is discussed in the weak field limit. Dragging of inertial frames, induced by angular momentum of the deflector, breaks spherical symmetry. I examine how the gravito-magnetic field affects image positions, caustics and critical curves. Distortion in microlensing-induced light curves is also considered.

The effect of currents of mass on bending of light rays is considered in the weak field regime. Following Fermat's principle and the standard theory of gravitational lensing, we derive the gravitomagnetic correction to time delay function and deflection angle caused by a geometrically-thin lens. The cases of both rotating and shifting deflectors are discussed.

We investigate the effects of the gravitomagnetic corrections to the usual gravitational lens quantities for a specific lensing mass distribution modelled after spiral galaxies. An exponential disk is embedded into two different spherical halo models where disk and haloes parameters are fixed according to the observed mass to light ratios, galaxy magnitudes and rotation curves. The general expressions for the lensing deflection angle are given also taking into account the orientation of the galaxy disk plane with respect to the lens plane. It is found that the gravitomagnetic term changes the deflection angle by a typical amount of the order of ten microarcseconds.

Gravitational lensing in metric theories of gravity is discussed. I introduce a generalized approximate metric element, inclusive of both post-post-Newtonian (ppN) contributions and gravito-magnetic field. Following Fermat's principle and standard hyphoteses, I derive the time delay function and deflection angle caused by an isolated mass distribution. Several astrophysical systems are considered. In most of the cases, the gravito-magnetic correction offers the best perspectives for an observational detection. Actual measurements distinguish only marginally different metric theories one from another.

Taking into account a torsion field gives rise to a negative pressure contribution in cosmological dynamics and then to an accelerated behaviour of Hubble fluid. The presence of torsion has the same effect of a Lambda - term. We obtain a general exact solution which well fits the data coming from high - redshift supernovae and Sunyaev - Zeldovich/X - ray method for the determination of the cosmological parameters. On the other hand, it is possible to obtain observational constraints on the amount of torsion density. A dust dominated Friedmann behaviour is recovered as soon as torsion effects are not relevant.

Following Sereno (2002a), we discuss the bending of light rays by spherically symmetric lenses with angular momentum. For several astrophysical systems, such as white dwarfs and galaxies, gravitomagnetism induces a correction on the deflection angle as large as 0.1%.

Observations of clusters of galaxies that gravitationally lens faint background galaxies can probe the amount and the equation of state, $\5$, of the dark energy (quintessence) in the universe. Provided that the mass profile and the mass normalization of the cluster are determined, it is possible to constrain the cosmological parameters that enter the lensing equations by means of the angular diameter distances, by locating (either by observations of giant arcs and magnification bias effect) the critical lines corresponding to known redshift source populations of galaxies. This method can help to distinguish between accelerating and decelerating models of the universe. Furthermore, since the position of critical lines is affected, especially in low-matter density universes, by the properties of quintessence, the observations of a suitable number of lensing clusters at intermediate redshifts can determine the equation of state. A very preliminary application of the method to the cluster CL 0024+1654 seems to support a flat accelerating universe dominated by dark energy.

We discuss the amplification dispersion in the observed luminosity of standard candles, like supernovae (SNe) of type Ia, induced by gravitational lensing in a Universe with dark energy (quintessence). We derive the main features of the magnification probability distribution function (pdf) of SNe in the framework of on average Friedmann-Lemaitre-Robertson-Walker (FLRW) models for both lensing by large-scale structures and compact objects. The magnification pdf is strongly dependent on the equation of state, $w_Q$, of the quintessence. The dispersion increases with the redshift of the source and is maximum for dark energy with very large negative pressure; the effects of gravitational lensing on the magnification pdf, i.e. the mode biased towards de-amplified values and the long tail towards large magnifications, are reduced for both microscopic DM and quintessence with an intermediate $w_Q$. Different equations of state of the dark energy can deeply change the dispersion in amplification for the projected observed samples of SNe Ia by future space-born missions. The "noise" in the Hubble diagram due to gravitational lensing strongly affects the determination of the cosmological parameters from SNe data. The errors on the pressureless matter density parameter, $\Omega_M$, and on $w_Q$ are maximum for quintessence with not very negative pressure. The effect of the gravitational lensing is of the same order of the other systematics affecting observations of SNe Ia. Due to the lensing by large-scale structures, in a flat Universe with $\Omega_M =0.4$, at $z=1$ a cosmological constant ($w_Q=-1$) can be interpreted as dark energy with $w_Q <-0.84$ (at 2-$\sigma$ confidence limit).

We investigate the evolution of matter density perturbations and some properties of the peculiar velocity field for a special class of exponential potentials in a scalar field model for quintessence, for which a general exact solution is known. The data from the 2-degree Field Galaxy Redshift Survey (2dFGRS) suggest a value of the today pressureless matter density Omega_M0 = 0.18 +- 0.05.

We present the first results of the analysis of data collected during the 1998-99 observational campaign at the 1.3 meter McGraw-Hill Telescope, towards the Andromeda galaxy (M31), aimed to the detection of gravitational microlensing effects as a probe of the presence of dark matter in our and in M31 halo. The analysis is performed using the pixel lensing technique, which consists in the study of flux variations of unresolved sources and has been proposed and implemented by the AGAPE collaboration. We carry out a shape analysis by demanding that the detected flux variations be achromatic and compatible with a Paczynski light curve. We apply the Durbin-Watson hypothesis test to the residuals. Furthermore, we consider the background of variables sources. Finally five candidate microlensing events emerge from our selection. Comparing with the predictions of a Monte Carlo simulation, assuming a standard spherical model for the M31 and Galactic haloes, and typical values for the MACHO mass, we find that our events are only marginally consistent with the distribution of observable parameters predicted by the simulation.

We investigate the properties of cosmological distances in locally inhomogeneous universes with pressureless matter and dark energy (quintessence), with constant equation of state. We give exact solutions for angular diameter distances in theempty beam approximation. In this hypothesis, the distance-redshift equation is derived fron the multiple lens-plane theory. The case of a flat universe is considered with particular attention. We show how this general scheme makes distances degenerate with respect to w_X and the smoothness parameters, alpha, accounting for the homogeneously distributed fraction of energy of the i-components. We analyse how this degeneracy influences the critical redshift where the angular diameter distance takes its maximum, and put in evidence future prospects for measuring the smoothness parameter of the pressureless matter, alpha_M.