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.

We present a novel proposal strategy for the Metropolis-Hastings algorithm designed to efficiently sample general convex polytopes in 100 or more dimensions. This improves upon previous sampling strategies used for free-form reconstruction of gravitational lenses, but is general enough to be applied to other fields. We have written a parallel implementation within the lens modeling framework GLASS. Testing shows that we are able to produce uniform uncorrelated random samples which are necessary for exploring the degeneracies inherent in lens reconstruction.

For a general class of analytic f(R)-gravity theories, we discuss the weak field limit in view of gravitational lensing. Though an additional Yukawa term in the gravitational potential modifies dynamics with respect to the standard Newtonian limit of General Relativity, the motion of massless particles results unaffected thanks to suitable cancellations in the post-Newtonian limit. Thus, all the lensing observables are equal to the ones known from General Relativity. Since f(R)-gravity is claimed, among other things, to be a possible solution to overcome for the need of dark matter in virialized systems, we discuss the impact of our results on the dynamical and gravitational lensing analyses. In this framework, dynamics could, in principle, be able to reproduce the astrophysical observations without recurring to dark matter, but in the case of gravitational lensing we find that dark matter is an unavoidable ingredient. Another important implication is that gravitational lensing, in the post-Newtonian limit, is not able to constrain these extended theories, since their predictions do not differ from General Relativity.

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.

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.