We discuss axion dark matter detection via two mechanisms: spontaneous decays and resonant conversion in neutron star magnetospheres. For decays, we show that the brightness temperature signal, rather than flux, is a less ambiguous measure for selecting candidate objects. This is owing principally to the finite beam width of telescopes which prevents one from being sensitive to the total flux from the object. With this in mind, we argue that the large surface-mass-density of the galactic centre or the Virgo cluster centre offers the best chance of improving current constraints on the axion-photon coupling via spontaneous decays. For the neutron star case, we first carry out a detailed study of mixing in magnetised plasmas. We derive transport equations for the axion-photon system via a controlled gradient expansion, allowing us to address inhomogeneous mass-shell constraints for arbitrary momenta. We then derive a non-perturbative Landau-Zener formula for the conversion probability valid across the range of relativistic and non-relativistic axions and show that the standard perturbative resonant conversion amplitude is a truncation of this result in the non-adiabatic limit. Our treatment reveals that that infalling dark matter axions typically convert non-adiabatically in magnetospheres. We describe the limitations of one-dimensional mixing equations and explain how three-dimensional effects activate new photon polarisations, including longitudinal modes and illustrate these arguments with numerical simulations in higher dimensions. We find that the bandwidth of the radio signal is dominated by Doppler broadening from the relative motion of the neutron star with respect to the observer. Therefore, we conclude that the radio signal from the resonant decay is weaker than previously thought, which means one relies on local density peaks to probe weaker axion-photon couplings.

We present a detailed study of the collapse of a spherical matter overdensity and the non-linear growth of large scale structures in the Galileon ghost condensate (GGC) model. This model is an extension of the cubic covariant Galileon (G3) which includes a field derivative of type $(\nabla_\mu\phi\nabla^\mu\phi)^2$ in the Lagrangian. We find that the cubic term activates the modifications in the main physical quantities and the additional term has an extended impact in the evolution of the latter. The GGC model shows largely mitigated effects in the linearised critical density contrast, non-linear effective gravitational coupling and the virial overdensity with respect to G3 but still preserves peculiar features with respect to the standard $\Lambda$CDM cosmological model, e.g., both the linear critical density contrast and the virial overdensity are larger than the ones in $\Lambda$CDM. The results of the spherical collapse model are then used to predict the evolution of the halo mass function, non-linear matter and lensing power spectra. While at low masses the GGC model presents about 10% fewer objects with respect to $\Lambda$CDM, at higher masses for $z>0$ it predicts 10% ($z=0.5$)-20% ($z=1$) more objects per comoving volume. Using a phenomenological approach to include the screening effect in the matter power spectrum, we show that the difference induced by the modifications of gravity are strongly dependent on the screening scale and that differences can be up to 60% with respect to $\Lambda$CDM. These differences translate to the lensing power spectrum where qualitatively the largest differences with respect to the standard cosmological model are for $\ell<10^3$. Depending on the screening scale, they can reach 60% on larger angular scales and then decrease for growing $\ell$. The results are obtained for the best fit parameters from cosmological data for each model.

The influence of considering a generalized dark matter (GDM) model, which allows for a non-pressure-less dark matter and a non-vanishing sound speed in the non-linear spherical collapse model is discussed for the Einstein-de Sitter-like (EdSGDM) and $\Lambda$GDM models. By assuming that the vacuum component responsible for the accelerated expansion of the Universe is not clustering and therefore behaving similarly to the cosmological constant $\Lambda$, we show how the change in the GDM characteristic parameters affects the linear density threshold for collapse of the non-relativistic component ($\delta_{\rm c}$) and its virial overdensity ($\Delta_{\rm V}$). We found that the generalized dark matter equation of state parameter $w_{\rm gdm}$ is responsible for lower values of the linear overdensity parameter as compared to the standard spherical collapse model and that this effect is much stronger than the one induced by a change in the generalized dark matter sound speed $c^2_{\rm s, gdm}$. We also found that the virial overdensity is only slightly affected and mostly sensitive to the generalized dark matter equation of state parameter $w_{\rm gdm}$. These effects could be relatively enhanced for lower values of the matter density. Finally, we found that the effects of the additional physics on $\delta_{\rm c}$ and $\Delta_{\rm V}$, when translated to non-linear observables such as the halo mass function, induce an overall deviation of about 40\% with respect to the standard $\Lambda$CDM model at late times for high mass objects. However, within the current linear constraints for $c^2_{\rm s, gdm}$ and $w_{\rm gdm}$, we found that these changes are the consequence of properly taking into account the correct linear matter power spectrum for the GDM model while the effects coming from modifications in the spherical collapse model remain negligible.

The evolution of the virial overdensity $\Delta_{\rm vir}$ for $\Lambda$CDM and seven dynamical dark-energy models is investigated in the extended spherical collapse model (SCM). Here the virialization process is naturally achieved by introducing shear and rotation instead of using the virial theorem. We generalise two approaches proposed in the literature and show that, regardless of the dark-energy model, the new virialization term can be calibrated on the peculiar velocity of the shell as measured from Einstein-de Sitter simulations. The two virialization recipes qualitatively reproduce the features of the ordinary SCM, i.e., a constant $\Delta_{\rm vir}$ for the EdS model and time-variation for dark-energy models, but without any mass dependence. Depending on the actual description of virialization and on the dark-energy model, the value of $\Delta_{\rm vir}$ varies between 10 and 40 percent. We use the new recipes to predict the surface-mass-density profile of dark matter haloes and the number of convergence density peaks for LSST- and Euclid-like weak lensing surveys.

We use the Equation of State (EoS) approach to study the evolution of the dark sector in Horndeski models, the most general scalar-tensor theories with second order equations of motion. By including the effects of the dark sector into our code EoS\_class, we demonstrate the numerical stability of the formalism and excellent agreement with results from other publicly available codes for a range of parameters describing the evolution of the function characterising the perturbations for Horndeski models, $\alpha_{\rm x}$, with ${\rm x}=\{{\rm K}, {\rm B}, {\rm M}, {\rm T}\}$. After demonstrating that on sub-horizon scales ($k\gtrsim 10^{-3}~{\rm Mpc}^{-1}$ at $z=0$) velocity perturbations in both the matter and the dark sector are typically subdominant with respect to density perturbations in the equation of state for perturbations, we find an attractor solution for the dark sector gauge-invariant density perturbation $\Delta_{\rm ds}$ by neglecting its time derivatives in the equation describing its time evolution, as commonly done in the well-known quasi-static approximation. Using this result, we provide simplified expressions for the equation-of-state functions: the dark sector entropy perturbations $w_{\rm ds}\Gamma_{\rm ds}$ and anisotropic stress $w_{\rm ds}\Pi_{\rm ds}$. From this we derive a growth factor-like equation for both matter and dark sector and are able to capture the relevant physics for several observables with great accuracy. We finally present new analytical expressions for the well-known modified gravity phenomenological functions $\mu$, $\eta$ and $\Sigma$ for a generic Horndeski model as functions of $\alpha_{\rm x}$. We show that on small scales they reproduce expressions presented in previous works, but on large scales, we find differences with respect to other works.

We derive the mass-temperature relation using an improved top-hat model and a continuous formation model which takes into account the effects of the ordered angular momentum acquired through tidal-torque interaction between clusters, random angular momentum, dynamical friction, and modifications of the virial theorem to include an external pressure term usually neglected. We show that the mass-temperature relation differs from the classical self-similar behavior, $M \propto T^{3/2}$, and shows a break at $3--4$ keV, and a steepening with a decreasing cluster temperature. We then compare our mass-temperature relation with those obtained in the literature with $N$-body simulations for $f(R)$ and symmetron models. We find that the mass-temperature relation is not a good probe to test gravity theories beyond Einstein's general relativity, because the mass-temperature relation of the $\Lambda$CDM model is similar to that of the modified gravity theories.

We investigate generalized Einstein-Aether theories that are compatible with the Planck Cosmic Microwave Background (CMB) temperature anisotropy, polarisation, and lensing data. For a given dark energy equation of state, $w_{\rm de}$, we formulate a designer approach and we investigate their impact on the CMB temperature anisotropy and matter power spectra. We use the Equation of State approach to parametrize the perturbations and find that this approach is particularly useful in identifying the most suitable and numerically efficient parameters to explore in a Markov chain Monte Carlo (MCMC) analysis. We find the data constrains models with $w_{\rm de} = -1$ to be compatible with $\Lambda$CDM. For $w_{\rm de} \not= -1$ models, which avoid the gravitational waves constraint through the entropy perturbation, we constrain $w_{\rm de}$ to be $w_{\rm de } = -1.06^{+0.08}_{-0.03}$ (CMB) and $w_{\rm de } =-1.04^{+0.05}_{-0.02}$ (CMB+Lensing) at $68\%$C.L., and find that these models can be different from $\Lambda$CDM and still be compatible with the data. We also find that these models can ameliorate some anomalies in $\Lambda$CDM when confronted with data, such as the low-$\ell$ and high-$k$ power in the CMB temperature anisotropy and matter power spectra respectively, but not simultaneously. We also investigate the anomalous lensing amplitude, quantified by $A_{\rm lens}$, and find that that for $w_{\rm de} = -1$ models, $A_{\rm lens} = 1.15^{+0.07}_{-0.08}$ (CMB) and $A_{\rm lens} = 1.12\pm0.05$ (CMB+Lensing) at $68\%$C.L. $\sim$ 2$\sigma$ larger than expected, similar to previous analyses of $\Lambda$CDM.

Using a semianalytic code, we show how baryon physics in a $\Lambda$CDM cosmology could solve the discrepancy between numerical predictions of dark matter haloes and observations, ranging from dwarf galaxies to clusters, without the need of nonstandard dark matter models as advocated, for example, by [Kaplinghat et al., Phys. Rev. Lett. 116, 041302, (2016)]. Combining well established results, we show, for the first time, how accounting for baryon physics, in particular dynamical friction mechanisms, leads to flat galaxy-cluster profiles and correlations in several of their properties, solves the so-called `diversity problem' and reproduces very well the challenging, extremely low-rising rotation curve of IC2574. We therefore suggest treating baryonic physics properly before introducing new exotic features, albeit legitimate, in the standard cosmological model.

We explore the constraints on dark sector models imposed by the recent observation of coincident gravitational waves and gamma rays from a binary neutron star merger, GW170817. Rather than focusing on specific models as has been considered by other authors, we explore this in the context of the equation of state approach of which the specific models are special cases. After confirming the strong constraints found by others for Horndeski, Einstein-Aether and massive gravity models, we discuss how it is possible to construct models which might evade the constraints from GW170817 but still leading to cosmologically interesting modifications to gravity. Possible examples are ``miracle cancellations" such as in $f(R)$ models, nonlocal models and higher-order derivatives. The latter two rely on the dimensionless ratio of the wave number of the observed gravitational waves to the Hubble expansion rate being very large ($\sim10^{19}$) which is used to suppress modifications to the speed of gravitational waves.

We review the Equation of State (EoS) approach to dark sector perturbations and apply it to $f(\mathcal{R})$ gravity models of dark energy. We show that the EoS approach is numerically stable and use it to set observational constraints on designer models. Within the EoS approach we build an analytical understanding of the dynamics of cosmological perturbations for the designer class of $f(\mathcal{R})$ gravity models, characterised by the parameter $B_0$ and the background equation of state of dark energy $w$. When we use the Planck Cosmic Microwave Background (CMB) temperature anisotropy, polarisation and lensing data as well as the Baryonic Acoustic Oscillation (BAO) data from SDSS and WiggleZ, we find $B_0<0.006$ (95\%CL) for the designer models with $w=-1$. Furthermore, we find $B_0<0.0045$ and $|w+1|<0.002$ (95\%CL) for the designer models with $w\neq -1$. Previous analyses found similar results for designer and Hu-Sawicki $f(\mathcal{R})$ gravity models using the Effective Field Theory (EFT) approach (Raveri et al. 2014, Hu et al. 2016), therefore this hints for the fact that generic $f(\mathcal{R})$ models with $w\neq-1$ can be tightly constrained by current cosmological data, complementary to solar system tests (Brax et al. 2008, Faulkner et al. 2007). When compared to a $w$CDM fluid with the same sound speed, we find that the equation of state for $f(\mathcal{R})$ models is better constrained to be close to -1 by about an order of magnitude, due to the strong dependence of the perturbations on $w$.

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

We compare Einstein-Boltzmann solvers that include modifications to General Relativity and find that, for a wide range of models and parameters, they agree to a high level of precision. We look at three general purpose codes that primarily model general scalar-tensor theories, three codes that model Jordan-Brans-Dicke (JBD) gravity, a code that models f(R) gravity, a code that models covariant Galileons, a code that models Ho\v{r}ava-Lifschitz gravity and two codes that model non-local models of gravity. Comparing predictions of the angular power spectrum of the cosmic microwave background and the power spectrum of dark matter for a suite of different models, we find agreement at the sub-percent level. This means that this suite of Einstein-Boltzmann solvers is now sufficiently accurate for precision constraints on cosmological and gravitational parameters.

Traditionally the spherical collapse of objects is studied with respect to a uniform background density, yielding the critical over-density $\delta_\mathrm{c}$ as key ingredient to the mass function of virialized objects. Here we investigate the shear and rotation acting on a peak in a Gaussian random field. By assuming that collapsing objects mainly form at those peaks, we use this shear and rotation as external effects changing the dynamics of the spherical collapse, which is described by the Raychaudhuri equation. We therefore assume that the shear and rotation have no additional dynamics on top of their cosmological evolution and thus only appear as inhomogeneities in the differential equation.

We investigate the evolution of cosmological perturbations in models of dark energy described by a time-like unit normalized vector field specified by a general function $\mathcal{F}(\mathcal{K})$, so-called Generalized Einstein-Aether models. First we study the background dynamics of such models via a designer approach in an attempt to model this theory as dark energy. We find that only one specific form of this designer approach matches $\Lambda$CDM at background order and we also obtain a differential equation which $\mathcal{F}(\mathcal{K})$ must satisfy for general $w$CDM cosmologies. We also present the equations of state for perturbations in Generalized Einstein-Aether models, which completely parametrize these models at the level of linear perturbations. A generic feature of modified gravity models is that they introduce new degrees of freedom. By fully eliminating these we are able to express the gauge invariant entropy perturbation and the scalar, vector, and tensor anisotropic stresses in terms of the perturbed fluid variables and metric perturbations only. These can then be used to study the evolution of perturbations in the scalar, vector, and tensor sectors and we use these to evolve the Newtonian gravitational potentials.

In this paper, extending past works of Del Popolo, we show how a high precision mass function (MF) can be obtained using the excursion set approach and an improved barrier taking implicitly into account a non-zero cosmological constant, the angular momentum acquired by tidal interaction of proto-structures and dynamical friction. In the case of the $\Lambda$CDM paradigm, we find that our MF is in agreement at the 3\% level to Klypin's Bolshoi simulation, in the mass range $M_{\rm vir} = 5 \times 10^9 h^{-1} M_{\odot} -- 5 \times 10^{14} h^{-1} M_{\odot}$ and redshift range $0 \lesssim z \lesssim 10$. For $z=0$ we also compared our MF to several fitting formulae, and found in particular agreement with Bhattacharya's within 3\% in the mass range $10^{12}-10^{16} h^{-1} M_{\odot}$. Moreover, we discuss our MF validity for different cosmologies.

We extend a previous work by Reischke et al., 2016 by studying the effects of tidal shear on clustering dark energy models within the framework of the extended spherical collapse model and using the Zel'dovich approximation. As in previous works on clustering dark energy, we assumed a vanishing effective sound speed describing the perturbations in dark energy models. To be self-consistent, our treatment is valid only on linear scales since we do not intend to introduce any heuristic models. This approach makes the linear overdensity $\delta_{\rm c}$ mass dependent and similarly to the case of smooth dark energy, its effects are predominant at small masses and redshifts. Tidal shear has effects of the order of percent or less, regardless of the model and preserves a well known feature of clustering dark energy: When dark energy perturbations are included, the models resemble better the $\Lambda$CDM evolution of perturbations. We also showed that effects on the comoving number density of halos are small and qualitatively and quantitatively in agreement with what previously found for smooth dark energy models.

A precise determination of the mass function is an important tool to verify cosmological predictions of the $\Lambda$CDM model and to infer more precisely the better model describing the evolution of the Universe. Galaxy clusters have been currently used to infer cosmological parameters, in particular the matter density parameter $\Omega_{\rm m}$, the matter power spectrum normalization $\sigma_8$ and the equation of state parameter $w_{\rm de}$ of the dark energy fluid. In this work, using data on massive galaxy clusters ($M>8\times 10^{14}~h^{-1}~M_{\odot}$) in the redshift range $0.05\lesssim z\lesssim 0.83$ we put constraints on the parameter $\alpha$ introduced within the formalism of the extended spherical collapse model to quantify deviations from sphericity due to shear and rotation. Since at the moment there is no physical model describing its functional shape, we assume it to be a logarithmic function of the cluster mass. By holding $\sigma_8$ fixed and restricting our analysis to a $\Lambda$CDM model, we find, at $1-\sigma$ confidence level, $\Omega_{\rm m}=0.284\pm0.0064$, $h=0.678\pm0.017$ and $\beta=0.0019^{+0.0008}_{-0.0015}$, where $\beta$ represents the slope of the parameter $\alpha$. This results translates into a $9\%$ decrement of the number of massive clusters with respect to a standard $\Lambda$CDM mass function, but better data are required to better constrain this quantity, since at the $2-\sigma$ and $3-\sigma$ confidence level we are only able to infer upper limits.

We study the nature of potentials in scalar field based models for dark energy - with both canonical and noncanonical kinetic terms. We calculate numerically, and using an analytic approximation around $a\approx 1$, potentials for models with constant equation-of-state parameter, $w_{\phi}$. We find that for a wide range of models with canonical and noncanonical kinetic terms there is a simple approximation for the potential that holds when the scale factor is in the range $0.6\lesssim a\lesssim 1.4$. We discuss how this form of the potential can also be used to represent models with non-constant $w_{\phi}$ and, hence, how it could be used in reconstruction from cosmological data.

We study the spherical collapse model in the presence of external gravitational tidal shear fields for different dark energy scenarios and investigate the impact on the mass function and cluster number counts. While previous studies of the influence of shear and rotation on $\delta_\mathrm{c}$ have been performed with heuristically motivated models, we try to avoid this model dependence and sample the external tidal shear values directly from the statistics of the underlying linearly evolved density field based on first order Lagrangian perturbation theory. Within this self-consistent approach, in the sense that we restrict our treatment to scales where linear theory is still applicable, only fluctuations larger than the scale of the considered objects are included into the sampling process which naturally introduces a mass dependence of $\delta_\mathrm{c}$. We find that shear effects are predominant for smaller objects and at lower redshifts, i. e. the effect on $\delta_\mathrm{c}$ is at or below the percent level for the $\Lambda$CDM model. For dark energy models we also find small but noticeable differences, similar to $\Lambda$CDM. The virial overdensity $\Delta_\mathrm{V}$ is nearly unaffected by the external shear. The now mass dependent $\delta_c$ is used to evaluate the mass function for different dark energy scenarios and afterwards to predict cluster number counts, which indicate that ignoring the shear contribution can lead to biases of the order of $1\sigma$ in the estimation of cosmological parameters like $\Omega_\mathrm{m}$, $\sigma_8$ or $w$.

In the present paper, we compare the predictions of two well known mechanisms considered able to solve the cusp/core problem (a. supernova feedback; b. baryonic clumps-DM interaction) by comparing their theoretical predictions to recent observations of the inner slopes of galaxies with masses ranging from dSphs to normal spirals. We compare the $\alpha$-$V_{\rm rot}$ and the $\alpha$-$M_{\ast}$ relationships, predicted by the two models with high resolution data coming from \citep{Adams2014,Simon2005}, LITTLE THINGS \citep{Oh2015}, THINGS dwarves \citep{Oh2011a,Oh2011b}, THINGS spirals \citep{Oh2015}, Sculptor, Fornax and the Milky Way. The comparison of the theoretical predictions with the complete set of data shows that the two models perform similarly, while when we restrict the analysis to a smaller subsample of higher quality, we show that the method presented in this paper (baryonic clumps-DM interaction) performs better than the one based on supernova feedback. We also show that, contrarily to the first model prediction, dSphs of small mass could have cored profiles. This means that observations of cored inner profiles in dSphs having a stellar mass $<10^6 M_{\odot}$ not necessarily imply problems for the $\Lambda$CDM model.

The accelerated expansion of the universe is a rather established fact in cosmology and many different models have been proposed as a viable explanation. Many of these models are based on the standard general relativistic framework of non-interacting fluids or more recently of coupled (interacting) dark energy models, where dark energy (the scalar field) is coupled to the dark matter component giving rise to a fifth-force. An interesting alternative is to couple the scalar field directly to the gravity sector via the Ricci scalar. These models are dubbed non-minimally coupled models and give rise to a time-dependent gravitational constant. In this work we study few models falling into this category and describe how observables depend on the strength of the coupling. We extend recent work on the subject by taking into account also the effects of the perturbations of the scalar field and showing their relative importance on the evolution of the mass function. By working in the framework of the spherical collapse model, we show that perturbations of the scalar field have a limited impact on the growth factor (for small coupling constant) and on the mass function with respect to the case where perturbations are neglected.

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

While in the standard cosmological model the accelerated expansion of the Universe is explained by invoking the presence of the cosmological constant term, it is still unclear the true origin of this stunning observational fact. It is therefore interesting to explore alternatives to the simplest scenario, in particular by assuming a more general framework where the fluid responsible of the accelerated expansion is characterised by a time-dependant equation of state. Usually these models, dubbed dark energy models, are purely phenomenological, but in this work we concentrate on a theoretically justified model, the ghost dark energy model. Within the framework of the spherical collapse model, we evaluate effects of dark energy perturbations both at the linear and non-linear level and transfer these results into an observable quantity, the mass function, by speculatively taking into account contributions of dark energy to the mass of the halos. We showed that the growth rate is higher in ghost models and that perturbations enhance the number of structures with respect to the $\Lambda$CDM model, with stronger effects when the total mass takes into account dark energy clumps.

The rate of structure formation in the Universe is different in homogeneous and clustered dark energy models. The degree of dark energy clustering depends on the magnitude of its effective sound speed $c^{2}_{\rm eff}$ and for $c_{\rm eff}=0$ dark energy clusters in a similar fashion to dark matter while for $c_{\rm eff}=1$ it stays (approximately) homogeneous. In this paper we consider two distinct equations of state for the dark energy component, $w_{\rm d}=const$ and $w_{\rm d}=w_0+w_1\left(\frac{z}{1+z}\right)$ with $c_{\rm eff}$ as a free parameter and we try to constrain the dark energy effective sound speed using current available data including SnIa, Baryon Acoustic Oscillation, CMB shift parameter ({\em Planck} and {\em WMAP}), Hubble parameter, Big Bang Nucleosynthesis and the growth rate of structures $f\sigma_{8}(z)$. At first we derive the most general form of the equations governing dark matter and dark energy clustering under the assumption that $c_{\rm eff}=const$. Finally, performing an overall likelihood analysis we find that the likelihood function peaks at $c_{\rm eff}=0$, however the dark energy sound speed is degenerate with respect to the cosmological parameters, namely $\Omega_{\rm m}$ and $w_{\rm d}$.

In this work we study the relevance of the cosmic web and substructures on the matter and lensing power spectra measured from halo mock catalogues extracted from the N-body simulations. Since N-body simulations are computationally expensive, it is common to use faster methods that approximate the dark matter field as a set of halos. In this approximation, we replace mass concentrations in N-body simulations by a spherically symmetric Navarro-Frenk-White halo density profile. We also consider the full mass field as the sum of two distinct fields: dark matter halos ($M>9\times 10^{12}~M_{\odot}$/h) and particles not included into halos. Mock halos reproduce well the matter power spectrum, but underestimate the lensing power spectrum on large and small scales. For sources at $z_{\rm s}=1$ the lensing power spectrum is underestimated by up to 40% at $\ell\approx 10^4$ with respect to the simulated halos. The large scale effect can be alleviated by combining the mock catalogue with the dark matter distribution outside the halos. In addition, to evaluate the contribution of substructures we have smeared out the intra-halo substructures in a N-body simulation while keeping the halo density profiles unchanged. For the matter power spectrum the effect of this smoothing is only of the order of 5%, but for lensing substructures are much more important: for $\ell\approx 10^4$ the internal structures contribute 30% of the total spectrum. These findings have important implications in the way mock catalogues have to be created, suggesting that some approximate methods currently used for galaxy surveys will be inadequate for future weak lensing surveys.

We present a simple model for describing intrinsic correlations for galaxy sizes based on the halo model. Studying these correlations is important both to improve our understanding of galaxy properties and because it is an important potential systematic for weak lensing size magnification measurements. Our model assumes that the density field drives these intrinsic correlations and we also model the distribution of satellite galaxies. We calculate the possible contamination to measurements of lensing convergence power spectrum from galaxy sizes, and show that the cross-correlation of intrinsic sizes with convergence is potentially an important systematic. We also explore how these intrinsic size correlations may affect surveys with different redshift depth. We find that, in this simple approach, intrinsic size correlations cannot be neglected in order to estimate lensing convergence power spectrum for constraining cosmological parameters.

In this work we investigate the spherical collapse model in flat FRW dark energy universes. We consider the Holographic Dark Energy (HDE) model as a dynamical dark energy scenario with a slowly time-varying equation-of-state (EoS) parameter $w_{\rm de}$ in order to evaluate the effects of the dark energy component on structure formation in the universe. We first calculate the evolution of density perturbations in the linear regime for both phantom and quintessence behavior of the HDE model and compare the results with standard Einstein-de Sitter (EdS) and $\Lambda$CDM models. We then calculate the evolution of two characterizing parameters in the spherical collapse model, i.e., the linear density threshold $\delta_{\rm c}$ and the virial overdensity parameter $\Delta_{\rm vir}$. We show that in HDE cosmologies the growth factor $g(a)$ and the linear overdensity parameter $\delta_{\rm c}$ fall behind the values for a $\Lambda$CDM universe while the virial overdensity $\Delta_{\rm vir}$ is larger in HDE models than in the $\Lambda$CDM model. We also show that the ratio between the radius of the spherical perturbations at the virialization and turn-around time is smaller in HDE cosmologies than that predicted in a $\Lambda$CDM universe. Hence the growth of structures starts earlier in HDE models than in $\Lambda$CDM cosmologies and more concentrated objects can form in this case. It has been shown that the non-vanishing surface pressure leads to smaller virial radius and larger virial overdensity $\Delta_{\rm vir}$. We compare the predicted number of halos in HDE cosmologies and find out that in general this value is smaller than for $\Lambda$CDM models at higher redshifts and we compare different mass function prescriptions. Finally, we compare the results of the HDE models with observations.

Dark matter and dark energy are usually assumed to couple only gravitationally. An extension to this picture is to model dark energy as a scalar field coupled directly to cold dark matter. This coupling leads to new physical effects, such as a fifth-force and a time-dependent dark matter particle mass. In this work we examine the impact that coupling has on weak lensing statistics by constructing realistic simulated weak-lensing maps using raytracing techniques through N-body cosmological simulations. We construct maps for different lensing quantities, covering a range of scales from a few arcminutes to several degrees. The concordance $\Lambda$CDM model is compared to different coupled dark energy models, described either by an exponential scalar field potential (standard coupled dark energy scenario) or by a SUGRA potential (bouncing model). We analyse several statistical quantities and our results, with sources at low redshifts are largely consistent with previous work on CMB lensing by Carbone et al., 2013. The most significant differences from the $\Lambda$CDM model are due to the enhanced growth of the perturbations and to the effective friction term in non-linear dynamics. For the most extreme models, we see differences in the power spectra up to 40% compared to the $\Lambda$CDM model. The different time evolution of the linear matter overdensity can account for most of the differences, but when controlling for this using a $\Lambda$CDM model having the same normalization, the overall signal is smaller due to the effect of the friction term appearing in the equation of motion for dark matter particles.

Under the commonly used assumption that clumped objects can be well described by a spherical top-hat matter density profile, we investigate the evolution of the cosmic growth index in clustering dark energy (CDE) scenarios on sub-horizon scales. We show that the evolution of the growth index $\gamma(z)$ strongly depends on the equation-of-state (EoS) parameter and on the clustering properties of the dark energy (DE) component. Performing a $\chi^2$ analysis, we show that CDE models have a better fit to observational growth rate data points with respect to the concordance $\Lambda$CDM model. We finally determine $\gamma(z)$ using an exponential parametrization and demonstrate that the growth index in CDE models presents large variations with cosmic redshift. In particular it is smaller (larger) than the theoretical value for the $\Lambda$CDM model, $\gamma_{\Lambda}\simeq0.55$, in the recent past (at the present time).

In the framework of the spherical collapse model we study the influence of shear and rotation terms for dark matter fluid in clustering dark energy models. We evaluate, for different equations of state, the effects of these terms on the linear overdensity threshold parameter, $\delta_{\rm c}$, and on the virial overdensity, $\Delta_{\rm V}$. The evaluation of their effects on $\delta_{\rm c}$ allows us to infer the modifications occurring on the mass function. Due to ambiguities in the definition of the halo mass in the case of clustering dark energy, we consider two different situations: the first is the classical one where the mass is of the dark matter halo only, while the second one is given by the sum of the mass of dark matter and dark energy. As previously found, the spherical collapse model becomes mass dependant and the two additional terms oppose to the collapse of the perturbations, especially on galactic scales, with respect to the spherical non-rotating model, while on clusters scales the effects of shear and rotation become negligible. The values for $\delta_{\rm c}$ and $\Delta_{\rm V}$ are higher than the standard spherical model. Regarding the effects of the additional non-linear terms on the mass function, we evaluate the number density of halos. As expected, major differences appear at high masses and redshifts. In particular, quintessence (phantom) models predict more (less) objects with respect to the $\Lambda$CDM model and the mass correction due to the contribution of the dark energy component has negligible effects on the overall number of structures.

We study the impact of Early Dark Energy fluctuations in the linear and non-linear regimes of structure formation. In these models the energy density of dark energy is non-negligible at high redshifts and the fluctuations in the dark energy component can have the same order of magnitude of dark matter fluctuations. Since two basic approximations usually taken in the standard scenario of quintessence models, that both dark energy density during the matter dominated period and dark energy fluctuations on small scales are negligible, are not valid in such models, we first study approximate analytical solutions for dark matter and dark energy perturbations in the linear regime. This study is helpful to find consistent initial conditions for the system of equations and to analytically understand the effects of Early Dark Energy and its fluctuations, which are also verified numerically. In the linear regime we compute the matter growth and variation of the gravitational potential associated with the Integrated Sachs-Wolf effect, showing that these observables present important modifications due to Early Dark Energy fluctuations, though making them more similar to $\Lambda$CDM model. We also make use of the Spherical Collapse model to study the influence of Early Dark Energy fluctuations in the nonlinear regime of structure formation, especially on $\delta_c$ parameter, and their contribution to the halo mass, which we show can be of the order of 10%. We finally compute how the number density of halos is modified in comparison to $\Lambda$CDM model and address the problem of how to correct the mass function in order to take into account the contribution of clustered dark energy. We conclude that the inhomogeneous Early Dark Energy models are more similar to $\Lambda$CDM model than its homogeneous counterparts.

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

The impacts of Compton scattering of hot cosmic gas with the cosmic microwave background radiation (Sunyaev-Zel'dovich effect, SZ) are consistently quantified in Gaussian and non-Gaussian scenarios, by means of 3D numerical, N-body, hydrodynamic simulations, including cooling, star formation, stellar evolution and metal pollution (He, C, O, Si, Fe, S, Mg, etc.) from different stellar phases, according to proper yields for individual metal species and mass-dependent stellar lifetimes. Light cones are built through the simulation outputs and samples of one hundred maps for the resulting temperature fluctuations are derived for both Gaussian and non-Gaussian primordial perturbations. From them, we estimate the possible changes due to early non-Gaussianities on: SZ maps, probability distribution functions, angular power spectra and corresponding bispectra. We find that the different growth of structures in the different cases induces significant spectral distortions only in models with large non-Gaussian parameters, $f_{\rm NL}$. In general, the overall trends are covered by the non-linear, baryonic evolution, whose feedback mechanisms tend to randomize the gas behaviour and homogenize its statistical features, quite independently from the background matter distribution. Deviations due to non-Gaussianity are almost undistinguishable for $f_{\rm NL}\lesssim 100$, remaining always at few-per-cent level, within the error bars of the Gaussian scenario. Rather extreme models with $f_{\rm NL}\sim1000$ present more substantial deviations from the Gaussian case, overcoming baryon contaminations and showing discrepancies up to a factor of a few in the spectral properties.

We study the effect of shear and rotation on results previously obtained dealing with the application of the spherical collapse model (SCM) to generalized Chaplygin gas (gCg) dominated universes. The system is composed of baryons and gCg and the collapse is studied for different values of the parameter $\alpha$ of the gCg. We show that the joint effect of shear and rotation is that of slowing down the collapse with respect to the simple SCM. This result is of utmost importance for the so-called unified dark matter models, since the described slow down in the growth of density perturbation can solve one of the main problems of the quoted models, namely the instability described in previous papers [e.g., H. B. Sandvik {\it et al.}, Phys. Rev. D {\bf 69}, 123524 (2004)] at the linear perturbation level.

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

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

We study, for the first time, how shear and angular momentum modify typical parameters of the spherical collapse model, in dark energy dominated universes. In particular, we study the linear density threshold for collapse $\delta_\mathrm{c}$ and the virial overdensity $\Delta_\mathrm{V}$, for several dark-energy models and its influence on the cumulative mass function. The equations of the spherical collapse are those obtained in Pace et al. (2010), who used the fully nonlinear differential equation for the evolution of the density contrast derived from Newtonian hydrodynamics, and assumed that dark energy is present only at the background level. With the introduction of the shear and rotation terms, the parameters of the spherical collapse model are now mass-dependant. The results of the paper show, as expected, that the new terms considered in the spherical collapse model oppose the collapse of perturbations on galactic scale giving rise to higher values of the linear overdensity parameter with respect to the non-rotating case. We find a similar effect also for the virial overdensity parameter. For what concerns the mass function, we find that its high mass tail is suppressed, while the low mass tail is slightly affected except in some cases, e.g. the Chaplygin gas case.

The influence of the shear stress and angular momentum on the nonlinear spherical collapse model is discussed in the framework of the Einstein-de Sitter (EdS) and $\Lambda$CDM models. By assuming that the vacuum component is not clustering within the homogeneous nonspherical overdensities, we show how the local rotation and shear affects the linear density threshold for collapse of the non-relativistic component ($\delta_\mathrm{c}$) and its virial overdensity ($\Delta_\mathrm{V}$). It is also found that the net effect of shear and rotation in galactic scale is responsible for higher values of the linear overdensity parameter as compared with the standard spherical collapse model (no shear and rotation).

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

(abridged) We explore the imprints of deviations from Gaussian primordial density fluctuations on the skeleton of the large-scale matter distribution as mapped through cosmological weak lensing. We computed the skeleton length of simulated effective convergence maps covering $\sim 35$ sq. deg each, extracted from a suite of cosmological $n-$body runs with different levels of local primordial non-Gaussianity. The latter is expected to alter the structure formation process with respect to the fiducial Gaussian scenario, and thus to leave a signature on the cosmic web. We found that alterations of the initial conditions consistently modify both the cumulative and the differential skeleton length, although the effect is generically smaller than the cosmic variance and depends on the smoothing of the map prior to the skeleton computation. Nevertheless, the qualitative shape of these deviations is rather similar to their primordial counterparts, implying that skeleton statistics retain good memory of the initial conditions. We performed a statistical analysis in order to find out at what Confidence Level primordial non-Gaussianity could be constrained by the skeleton test on cosmic shear maps of the size we adopted. At 68.3% Confidence Level we found an error on the measured level of primordial non-Gaussianity of $\Delta f_\mathrm{NL}\sim 300$, while at 90% Confidence Level it is of $\Delta f_\mathrm{NL}\sim 500$. While these values by themselves are not competitive with the current constraints, weak lensing maps larger than those used here would have a smaller field-to-field variance, and thus would likely lead to tighter constraints. A rough estimate indicates $\Delta f_\mathrm{NL} \sim$ a few tens at 68.3% Confidence Level for an all-sky weak lensing survey.

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

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

Strong lensing is one of the most direct probes of the mass distribution in the inner regions of galaxy clusters. It can be used to constrain the density profiles and to measure the mass of the lenses. Moreover, the abundance of strong lensing events can be used to constrain the structure formation and the cosmological parameters through the so-called "arc-statistics" approach. However, several issues related to the usage of strong lensing clusters in cosmological applications are still controversial, leading to the suspect that several biases may affect this very peculiar class of objects. With this study we aim at better understanding the properties of galaxy clusters which can potentially act as strong lenses. We do so by investigating the properties of a large sample of galaxy clusters extracted from the N-body/hydrodynamical simulation MareNostrum Universe. We explore the correlation between the cross section for lensing and many properties of clusters, like the mass, the three-dimensional and projected shapes, their concentrations, the X-ray luminosity and the dynamical activity. We find that the probability of strong alignments between the major axes of the lenses and the line of sight is a growing function of the lensing cross section. In projection, the strong lenses appear rounder within R200, but we find that their cores tend to be more elliptical as the lensing cross section increases. We also find that the cluster concentrations estimated from the projected density profiles tend to be biased high. The X-ray luminosity of strong lensing clusters is higher than that of normal lenses of similar mass and redshift. This is particular significant for the least massive lenses. Finally, we find that the strongest lenses generally exhibit an excess of kinetic energy within the virial radius, indicating that they are more dynamically active than usual clusters.

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

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

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