The ratio between the CMB lensing/galaxy counts and the galaxy shear/galaxy counts cross-correlations combines the information from different cosmological probes to infer cosmographic measurements that are less dependent on astrophysical uncertainties and can constrain the geometry of the Universe. We discuss the future perspectives for the measurement of this lensing ratio as previously introduced, i.e. with the use of the Limber and flat-sky approximations and neglecting all the effects on the galaxy survey from observing on the past lightcone. We then show how the cosmological information in this estimator is affected by the Limber approximation and by the inclusion of the redshift space distorsions (RSD) and magnification lensing contributions to the galaxy number counts. We find that the magnification lensing contribution induces a multipole dependence of the lensing ratio that we show to be detectable at a statistical significant level combining post-$Planck$ CMB surveys and a Euclid-like experiment. We propose an improved estimator which takes into account this angular scale dependence. Using this extended formalism, we present forecasts for upcoming and future cosmological surveys and we show at which extent the lensing ratio information can improve the CMB constraints on cosmological parameters. We get that for extended cosmological models where the neutrino mass, the spatial curvature and the dark energy equation of state are allowed to vary, the constraints from $Planck$ on these parameters and on $H_0$ can be reduced by $\sim 40\%$ with the inclusion of a single lensing ratio and by $\sim 60-70\%$ adding the joint measurement of 9 lensing ratios with a Euclid-like survey. We also find that neglegcting the contribution from magnification lensing can induce a bias on the derived cosmological parameters in a combined analysis.

The data from the Euclid mission will enable the measurement of the photometric redshifts, angular positions, and weak lensing shapes for over a billion galaxies. This large dataset will allow for cosmological analyses using the angular clustering of galaxies and cosmic shear. The cross-correlation (XC) between these probes can tighten constraints and it is therefore important to quantify their impact for Euclid. In this study we carefully quantify the impact of XC not only on the final parameter constraints for different cosmological models, but also on the nuisance parameters. In particular, we aim at understanding the amount of additional information that XC can provide for parameters encoding systematic effects, such as galaxy bias or intrinsic alignments (IA). We follow the formalism presented in Euclid Collaboration: Blanchard et al. (2019) and make use of the codes validated therein. We show that XC improves the dark energy Figure of Merit (FoM) by a factor $\sim 5$, whilst it also reduces the uncertainties on galaxy bias by $\sim 17\%$ and the uncertainties on IA by a factor $\sim 4$. We observe that the role of XC on the final parameter constraints is qualitatively the same irrespective of the galaxy bias model used. We also show that XC can help in distinguishing between different IA models, and that if IA terms are neglected then this can lead to significant biases on the cosmological parameters. We find that the XC terms are necessary to extract the full information content from the data in future analyses. They help in better constraining the cosmological model, and lead to a better understanding of the systematic effects that contaminate these probes. Furthermore, we find that XC helps in constraining the mean of the photometric-redshift distributions, but it requires a more precise knowledge of this mean in order not to degrade the final FoM. [Abridged]

Stage IV weak lensing experiments will offer more than an order of magnitude leap in precision. We must therefore ensure that our analyses remain accurate in this new era. Accordingly, previously ignored systematic effects must be addressed. In this work, we evaluate the impact of the reduced shear approximation and magnification bias, on the information obtained from the angular power spectrum. To first-order, the statistics of reduced shear, a combination of shear and convergence, are taken to be equal to those of shear. However, this approximation can induce a bias in the cosmological parameters that can no longer be neglected. A separate bias arises from the statistics of shear being altered by the preferential selection of galaxies and the dilution of their surface densities, in high-magnification regions. The corrections for these systematic effects take similar forms, allowing them to be treated together. We calculated the impact of neglecting these effects on the cosmological parameters that would be determined from Euclid, using cosmic shear tomography. To do so, we employed the Fisher matrix formalism, and included the impact of the super-sample covariance. We also demonstrate how the reduced shear correction can be calculated using a lognormal field forward modelling approach. These effects cause significant biases in Omega_m, sigma_8, n_s, Omega_DE, w_0, and w_a of -0.53 sigma, 0.43 sigma, -0.34 sigma, 1.36 sigma, -0.68 sigma, and 1.21 sigma, respectively. We then show that these lensing biases interact with another systematic: the intrinsic alignment of galaxies. Accordingly, we develop the formalism for an intrinsic alignment-enhanced lensing bias correction. Applying this to Euclid, we find that the additional terms introduced by this correction are sub-dominant.

Stage IV lensing surveys promise to make available an unprecedented amount of excellent data which will represent a huge leap in terms of both quantity and quality. This will open the way to the use of novel tools, which go beyond the standard second order statistics probing the high order properties of the convergence field. We discuss the use of Minkowski Functionals (MFs) as complementary probes to increase the lensing Figure of Merit (FoM), for a survey made out of a wide total area $A_{\rm{tot}}$ imaged at a limiting magnitude $\rm{mag_{W}}$ containing a subset of area $A_{\rm{deep}}$ where observations are pushed to a deeper limiting magnitude $\rm{mag_{D}}$. We present an updated procedure to match the theoretically predicted MFs to the measured ones, taking into account the impact of map reconstruction from noisy shear data. We validate this renewed method against simulated data sets with different source redshift distributions and total number density, setting these quantities in accordance with the depth of the survey. We can then rely on a Fisher matrix analysis to forecast the improvement in the FoM due to the joint use of shear tomography and MFs under different assumptions on $(A_{\rm{tot}},\,A_{\rm{deep}},\,\rm{mag_{D}})$, and the prior on the MFs nuisance parameters. It turns out that MFs can provide a valuable help in increasing the FoM of the lensing survey, provided the nuisance parameters are known with a non negligible precision. What is actually more interesting is the possibility to compensate for the loss of FoM due to a cut in the multipole range probed by shear tomography, which makes the results more robust against uncertainties in the modeling of nonlinearities. This makes MFs a promising tool to both increase the FoM and make the constraints on the cosmological parameters less affected by theoretical systematic effects.

Our aim is to quantify the impact of systematic effects on the inference of cosmological parameters from cosmic shear. We present an end-to-end approach that introduces sources of bias in a modelled weak lensing survey on a galaxy-by-galaxy level. Residual biases are propagated through a pipeline from galaxy properties (one end) through to cosmic shear power spectra and cosmological parameter estimates (the other end), to quantify how imperfect knowledge of the pipeline changes the maximum likelihood values of dark energy parameters. We quantify the impact of an imperfect correction for charge transfer inefficiency (CTI) and modelling uncertainties of the point spread function (PSF) for Euclid, and find that the biases introduced can be corrected to acceptable levels.

The Euclid space telescope will measure the shapes and redshifts of galaxies to reconstruct the expansion history of the Universe and the growth of cosmic structures. Estimation of the expected performance of the experiment, in terms of predicted constraints on cosmological parameters, has so far relied on different methodologies and numerical implementations, developed for different observational probes and for their combination. In this paper we present validated forecasts, that combine both theoretical and observational expertise for different cosmological probes. This is presented to provide the community with reliable numerical codes and methods for Euclid cosmological forecasts. We describe in detail the methodology adopted for Fisher matrix forecasts, applied to galaxy clustering, weak lensing and their combination. We estimate the required accuracy for Euclid forecasts and outline a methodology for their development. We then compare and improve different numerical implementations, reaching uncertainties on the errors of cosmological parameters that are less than the required precision in all cases. Furthermore, we provide details on the validated implementations that can be used by the reader to validate their own codes if required. We present new cosmological forecasts for Euclid. We find that results depend on the specific cosmological model and remaining freedom in each setup, i.e. flat or non-flat spatial cosmologies, or different cuts at nonlinear scales. The validated numerical implementations can now be reliably used for any setup. We present results for an optimistic and a pessimistic choice of such settings. We demonstrate that the impact of cross-correlations is particularly relevant for models beyond a cosmological constant and may allow us to increase the dark energy Figure of Merit by at least a factor of three.

Minkowski functionals (MFs) quantify the topological properties of a given field probing its departure from Gaussianity. We investigate their use on lensing convergence maps in order to see whether they can provide further insights on the underlying cosmology with respect to the standard second-order statistics, i.e., cosmic shear tomography. To this end, we first present a method to match theoretical predictions with measured MFs taking care of the shape noise, imperfections in the map reconstruction, and inaccurate description of the nonlinearities in the matter power spectrum and bispectrum. We validate this method against simulated maps reconstructed from shear fields generated by the MICE simulation. We then perform a Fisher matrix analysis to forecast the accuracy on cosmological parameters from a joint MFs and shear tomography analysis. It turns out that MFs are indeed helpful to break the $\Omega_{\rm m}$--$\sigma_8$ degeneracy thus generating a sort of chain reaction leading to an overall increase of the Figure of Merit.

Accurate shape measurements are essential to infer cosmological parameters from large area weak gravitational lensing studies. The compact diffraction-limited point-spread function (PSF) in space-based observations is greatly beneficial, but its chromaticity for a broad band observation can lead to new subtle effects that could hitherto be ignored: the PSF of a galaxy is no longer uniquely defined and spatial variations in the colours of galaxies result in biases in the inferred lensing signal. Taking Euclid as a reference, we show that this colourgradient bias (CG bias) can be quantified with high accuracy using available multi-colour Hubble Space Telescope (HST) data. In particular we study how noise in the HST observations might impact such measurements and find this to be negligible. We determine the CG bias using HST observations in the F606W and F814W filters and observe a correlation with the colour, in line with expectations, whereas the dependence with redshift is weak. The biases for individual galaxies are generally well below 1%, which may be reduced further using morphological information from the Euclid data. Our results demonstrate that CG bias should not be ignored, but it is possible to determine its amplitude with sufficient precision, so that it will not significantly bias the weak lensing measurements using Euclid data.

The unprecedented quality, the increased dataset, and the wide area of ongoing and near future weak lensing surveys allows to move beyond the standard two points statistics thus making worthwhile to investigate higher order probes. As an interesting step towards this direction, we expolore the use of higher order moments (HOM) of the convergence field as a way to increase the lensing Figure of Merit (FoM). To this end, we rely on simulated convergence to first show that HOM can be measured and calibrated so that it is indeed possible to predict them for a given cosmological model provided suitable nuisance parameters are introduced and then marginalized over. We then forecast the accuracy on cosmological parameters from the use of HOM alone and in combination with standard shear power spectra tomography. It turns out that HOM allow to break some common degeneracies thus significantly boosting the overall FoM. We also qualitatively discuss possible systematics and how they can be dealt with.

Euclid and the Large Synoptic Survey Telescope (LSST) are poised to dramatically change the astronomy landscape early in the next decade. The combination of high cadence, deep, wide-field optical photometry from LSST with high resolution, wide-field optical photometry and near-infrared photometry and spectroscopy from Euclid will be powerful for addressing a wide range of astrophysical questions. We explore Euclid/LSST synergy, ignoring the political issues associated with data access to focus on the scientific, technical, and financial benefits of coordination. We focus primarily on dark energy cosmology, but also discuss galaxy evolution, transient objects, solar system science, and galaxy cluster studies. We concentrate on synergies that require coordination in cadence or survey overlap, or would benefit from pixel-level co-processing that is beyond the scope of what is currently planned, rather than scientific programs that could be accomplished only at the catalog level without coordination in data processing or survey strategies. We provide two quantitative examples of scientific synergies: the decrease in photo-z errors (benefitting many science cases) when high resolution Euclid data are used for LSST photo-z determination, and the resulting increase in weak lensing signal-to-noise ratio from smaller photo-z errors. We briefly discuss other areas of coordination, including high performance computing resources and calibration data. Finally, we address concerns about the loss of independence and potential cross-checks between the two missions and potential consequences of not collaborating.

We present a conjugate analysis of two different dark energy models investigating both their agreement with recent data and their thermodynamical properties. The successful match with the data allows to both constrain the model parameters and characterize their kinematical properties. As a novel step, we exploit the strong connection between gravity and thermodynamics to further check models viability by investigating their thermodynamical quantities. In particular, we study whether the cosmological scenario fulfills the generalized second law of thermodynamics and, moreover, we contrast the two model asking whether the evolution of the total entropy is in agreement with the expectation for a closed system. As a general result, we discuss whether thermodynamic constraints can be a valid complementary way to both constrain dark energy models and differentiate among rival scenarios.

Gravitational lensing has long been considered as a valuable tool to determine the total mass of galaxy clusters. The shear profile as inferred from the statistics of ellipticity of background galaxies allows to probe the cluster intermediate and outer regions thus determining the virial mass estimate. However, the mass sheet degeneracy and the need for a large number of background galaxies motivate the search for alternative tracers which can break the degeneracy among model parameters and hence improve the accuracy of the mass estimate. Lensing flexion, i.e. the third derivative of the lensing potential, has been suggested as a good answer to the above quest since it probes the details of the mass profile. We investigate here whether this is indeed the case considering jointly using weak lensing, magnification and flexion. We use a Fisher matrix analysis to forecast the relative improvement in the mass accuracy for different assumptions on the shear and flexion signal - to - noise (S/N) ratio also varying the cluster mass, redshift, and ellipticity. It turns out that the error on the cluster mass may be reduced up to a factor 2 for reasonable values of the flexion S/N ratio. As a general result, we get that the improvement in mass accuracy is larger for more flattened haloes, but extracting general trends is a difficult because of the many parameters at play. We nevertheless find that flexion is as efficient as magnification to increase the accuracy in both mass and concentration determination.

Large area lensing surveys are expected to make it possible to use cosmic shear tomography as a tool to severely constrain cosmological parameters. To this end, one typically relies on second order statistics such as the two - point correlation fucntion and its Fourier counterpart, the power spectrum. Moving a step forward, we wonder whether and to which extent higher order stastistics can improve the lensing Figure of Merit (FoM). In this first paper of a series, we investigate how second, third and fourth order lensing convergence moments can be measured and use as probe of the underlying cosmological model. We use simulated data and investigate the impact on moments estimate of the map reconstruction procedure, the cosmic variance, and the intrinsic ellipticity noise. We demonstrate that, under realistic assumptions, it is indeed possible to use higher order moments as a further lensing probe.

Strong lensing provides popular techniques to investigate the mass distribution of intermediate redshift galaxies, testing galaxy evolution and formation scenarios. It especially probes the background cosmic expansion, hence constraining cosmological parameters. The measurement of Einstein radii and central velocity dispersions indeed allows to trace the ratio D_s/D_ls between the distance D_s from the observer to the source and the distance D_ls from the lens to the source. We present an improved method to explicitly include the two - component structure in the galaxy lens modeling, in order to analyze the role played by the redshift and the model dependence on a nuisance parameter, F_E, which is usually marginalized in the cosmological applications. We show how to deal with these problems and carry on a Fisher matrix analysis to infer the accuracy on cosmological parameters achieved by this method.

Euclid is the next ESA mission devoted to cosmology. It aims at observing most of the extragalactic sky, studying both gravitational lensing and clustering over $\sim$15,000 square degrees. The mission is expected to be launched in year 2020 and to last six years. The sheer amount of data of different kinds, the variety of (un)known systematic effects and the complexity of measures require efforts both in sophisticated simulations and techniques of data analysis. We review the mission main characteristics, some aspects of the the survey and highlight some of the areas of interest to this meeting

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.

Alterations of the gravity Lagrangian introduced in modified torsion gravity theories---also referred to as $f(T)$ gravity---allows for an accelerated expansion in a matter dominated Universe. In this framework, the cosmic speed up is driven by an effective `torsion fluid'. Besides the background evolution of the Universe, structure formation is also modified because of a time dependent effective gravitational constant. Here, we investigate the imprints of $f(T)$ gravity on galaxy clustering and weak gravitational lensing to the aim of understanding whether future galaxy surveys could constrain torsion gravity and discriminate amongst it and standard general relativity. Specifically, we compute Fisher matrix forecasts for two viable $f(T)$ models to both infer the accuracy on the measurement of the model parameters and evaluate the power that a combined clustering and shear analysis will have as a tool for model selection. We find that with such a combination of probes it will indeed be possible to tightly constrain $f(T)$ model parameters. Moreover, the Occam's razor provided by the Bayes factor will allow us to confirm an $f(T)$ power-law extension of the concordance $\Lambda{\rm CDM}$ model, were a value larger than 0.02 of its power-law slope measured, whereas in $\Lambda{\rm CDM}$ it is exactly 0.

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.

Beyond the standard cosmological model the late-time accelerated expansion of the universe can be reproduced by the introduction of an additional dynamical scalar field. In this case, the field is expected to be naturally coupled to the rest of the theory's fields, unless a (still unknown) symmetry suppresses this coupling. Therefore, this would possibly lead to some observational consequences, such as space-time variations of nature's fundamental constants. In this paper we investigate the coupling between a dynamical Dark Energy model and the electromagnetic field, and the corresponding evolution of the fine structure constant ($\alpha$) with respect to the standard local value $\alpha_0$. In particular, we derive joint constraints on two dynamical Dark Energy model parametrizations (the Chevallier-Polarski-Linder and Early Dark Energy model) and on the coupling with electromagnetism $\zeta$, forecasting future low-medium redshift observations. We combine supernovae and weak lensing measurements from the Euclid experiment with high-resolution spectroscopy measurements of fundamental couplings and the redshift drift from the European Extremely Large Telescope, highlighting the contribution of each probe. Moreover, we also consider the case where the field driving the $\alpha$ evolution is not the one responsible for cosmic acceleration and investigate how future observations can constrain this scenario.

MOdified Newtonian dynamics (MOND) represents a phenomenological alternative to dark matter (DM) for the missing mass problem in galaxies and clusters of galaxies. We analyze the central regions of a local sample of $\sim 220$ early-type galaxies from the $\rm ATLAS^{3D}$ survey, to see if the data can be reproduced without recourse to DM. We estimate dynamical masses in the MOND context through Jeans analysis, and compare to $\rm ATLAS^{3D}$ stellar masses from stellar population synthesis. We find that the observed stellar mass--velocity dispersion relation is steeper than expected assuming MOND with a fixed stellar initial mass function (IMF) and a standard value for the acceleration parameter $a_{\rm 0}$. Turning from the space of observables to model space, a) fixing the IMF, a universal value for $a_{\rm 0}$ cannot be fitted, while, b) fixing $a_{\rm 0}$ and leaving the IMF free to vary, we find that it is "lighter" (Chabrier-like) for low-dispersion galaxies, and "heavier" (Salpeter-like) for high dispersions. This MOND-based trend matches inferences from Newtonian dynamics with DM, and from detailed analysis of spectral absorption lines, adding to the converging lines of evidence for a systematically-varying IMF.

Gamma -ray bursts (GRBs) observed up to redshifts $z>9.4$ can be used as possible probes to test cosmological models. Here we show how changes of the slope of the {\it luminosity $L^*_X$ -break time $T^*_a$} correlation in GRB afterglows, hereafter the LT correlation, affect the determination of the cosmological parameters. With a simulated data set of 101 GRBs with a central value of the correlation slope that differs on the intrinsic one by a $5\sigma$ factor, we find an overstimated value of the matter density parameter, $\Omega_M$, compared to the value obtained with SNe Ia, while the Hubble constant, $H_0$, best fit value is still compatible in 1$\sigma$ compared to other probes. We show that this compatibility of $H_0$ is due to the large intrinsic scatter associated with the simulated sample. Instead, if we consider a subsample of high luminous GRBs ($HighL$), we find that both the evaluation of $H_0$ and $\Omega_M$ are not more compatible in 1$\sigma$ and $\Omega_M$ is underestimated by the $13\%$. However, the $HighL$ sample choice reduces dramatically the intrinsic scatter of the correlation, thus possibly identifying this sample as the standard canonical `GRBs' confirming previous results presented in Dainotti et al. (2010,2011). Here, we consider the LT correlation as an example, but this reasoning can be extended also for all other GRB correlations. In literature so far GRB correlations are not corrected for redshift evolution and selection biases, therefore we are not aware of their intrinsic slopes and consequently how far the use of the observed correlations can influence the derived `best' cosmological settings. Therefore, we conclude that any approach that involves cosmology should take into consideration only intrinsic correlations not the observed ones.

Cosmic shear tomography has emerged as one of the most promising tools to both investigate the nature of dark energy and discriminate between General Relativity and modified gravity theories. In order to successfully achieve these goals, systematics in shear measurements have to be taken into account; their impact on the weak lensing power spectrum has to be carefully investigated in order to estimate the bias induced on the inferred cosmological parameters. To this end, we develop here an efficient tool to compute the power spectrum of systematics by propagating, in a realistic way, shear measurement, source properties and survey setup uncertainties. Starting from analytical results for unweighted moments and general assumptions on the relation between measured and actual shear, we derive analytical expressions for the multiplicative and additive bias, showing how these terms depend not only on the shape measurement errors, but also on the properties of the source galaxies (namely, size, magnitude and spectral energy distribution). We are then able to compute the amplitude of the systematics power spectrum and its scaling with redshift, while we propose a multigaussian expansion to model in a non-parametric way its angular scale dependence. Our method allows to self-consistently propagate the systematics uncertainties to the finally observed shear power spectrum, thus allowing us to quantify the departures from the actual spectrum. We show that even a modest level of systematics can induce non-negligible deviations, thus leading to a significant bias on the recovered cosmological parameters.

As soon as their extragalactic origins were established, the hope to make Gamma - Ray Bursts (GRBs) standardizeable candles to probe the very high - z universe has opened the search for scaling relations between redshift independent observable quantities and distance dependent ones. Although some remarkable success has been achieved, the empirical correlations thus found are still affected by a significant intrinsic scatter which downgrades the precision in the inferred GRBs Hubble diagram. We investigate here whether this scatter may come from fitting together objects belonging to intrinsically different classes. To this end, we rely on a cladistics analysis to partition GRBs in homogenous families according to their rest frame properties. Although the poor statistics prevent us from drawing a definitive answer, we find that both the intrinsic scatter and the coefficients of the $E_{peak}$\,-\,$E_{iso}$ and $E_{peak}$\,-\,$L$ correlations significantly change depending on which subsample is fitted. It turns out that the fit to the full sample leads to a scaling relation which approximately follows the diagonal of the region delimited by the fits to each homogenous class. We therefore argue that a preliminary identification of the class a GRB belongs to is necessary in order to select the right scaling relation to be used in order to not bias the distance determination and hence the Hubble diagram.

Weak gravitational lensing by galaxy clusters on faint higher redshift galaxies has been traditionally used to study the cluster mass distribution and as a tool to identify clusters as peaks in the shear maps. However, it becomes soon clear that peaks statistics can also be used as a way to constrain the underlying cosmological model due to its dependence on both the cosmic expansion rate and the growth rate of structures. This feature makes peak statistics particularly interesting from the point of view of discriminating between General Relativity and modified gravity. Here we consider a general class of $f(R)$ theories and compute the observable mass function based on the aperture mass statistics. We complement our theoretical analysis with a Fisher matrix forecast of the constraints that an Euclid\,-\,like survey can impose on the $f(R)$ model parameters. We show that peak statistics alone can in principle discriminate between General Relativity and $f(R)$ models and strongly constrain the $f(R)$ parameters that are sensitive to the non-linear growth of structure. However, further analysis is needed in order to include possible selection function in the peaks redshift determination.

In this paper, in the framework of the secondary infall model, the correlation between the central surface density and the halo core radius of galaxy, and cluster of galaxies, dark matter haloes was analyzed, this having recently been studied on a wide range of scales. We used Del Popolo (2009) secondary infall model taking into account ordered and random angular momentum, dynamical friction, and dark matter (DM) adiabatic contraction to calculate the density profile of haloes, and then these profiles are used to determine the surface density of DM haloes. The main result is that $r_\ast$ (the halo characteristic radius) is not an universal quantity as claimed by Donato et al. (2009) and Gentile et al. (2009). On the contrary, we find a correlation with the halo mass $M_{200}$ in agreement with Cardone & Tortora (2010), Boyarsky at al. (2009) and Napolitano et al. (2010), but with a significantly smaller scatter, namely $0.16 \pm 0.05$. We also consider the baryon column density finding this latter being indeed a constant for low mass systems such as dwarfs, but correlating with mass with a slope $\alpha= 0.18 \pm 0.05$. In the case of the surface density of dark matter for a system composed only of dark matter, as in dissipationless simulations, we get $\alpha=0.20 \pm 0.05$. These results leave little room for the recently claimed universality of (dark and stellar) column density.

The spatial variation of the colour of a galaxy may introduce a bias in the measurement of its shape if the PSF profile depends on wavelength. We study how this bias depends on the properties of the PSF and the galaxies themselves. The bias depends on the scales used to estimate the shape, which may be used to optimise methods to reduce the bias. Here we develop a general approach to quantify the bias. Although applicable to any weak lensing survey, we focus on the implications for the ESA Euclid mission. Based on our study of synthetic galaxies we find that the bias is a few times 10^-3 for a typical galaxy observed by Euclid. Consequently, it cannot be neglected and needs to be accounted for. We demonstrate how one can do so using spatially resolved observations of galaxies in two filters. We show that HST observations in the F606W and F814W filters allow us to model and reduce the bias by an order of magnitude, sufficient to meet Euclid's scientific requirements. The precision of the correction is ultimately determined by the number of galaxies for which spatially-resolved observations in at least two filters are available. We use results from the Millennium Simulation to demonstrate that archival HST data will be sufficient for the tomographic cosmic shear analysis with the Euclid dataset.

We revisit the issue of the constancy of the dark matter (DM) and baryonic Newtonian acceleration scales within the DM scale radius by considering a large sample of late - type galaxies. We rely on a Markov Chain Monte Carlo (MCMC) method to estimate the parameters of the halo model and the stellar mass - to - light ratio and then propagate the uncertainties from the rotation curve data to the estimate of the acceleration scales. This procedure allows us to compile a catalog of 58 objects with estimated values of the $B$ band absolute magnitude $M_B$, the virial mass $M_{vir}$, the DM and baryonic Newtonian accelerations (denoted as $g_{DM}(r_0)$ and $g_{bar}(r_0)$, respectively) within the scale radius $r_0$ which we use to investigate whether it is possible to define a universal acceleration scale. We find a weak but statistically meaningful correlation with $M_{vir}$ thus making us argue against the universality of the acceleration scales. However, the results somewhat depend on the sample adopted so that a careful analysis of selection effects should be carried out before any definitive conclusion can be drawn.

The assumptions that "light propagates along null geodesics of the spacetime metric" and "the number of photons is conserved along the light path" lead to the distance duality relation (DDR), $\eta = D_L(z) (1 + z)^{-2}/D_A(z) = 1$, with $D_L(z)$ and $D_A(z)$ the luminosity and angular diameter distances to a source at redshift $z$. In order to test the DDR, we follow the usual strategy comparing the angular diameter distances of a set of clusters, inferred from X - ray and radio data, with the luminosity distance at the same cluster redshift using the local regression technique to estimate $D_L(z)$ from Type Ia Supernovae (SNeIa) Hubble diagram. In order to both strengthen the constraints on the DDR and get rid of the systematics related to the unknown cluster geometry, we also investigate the possibility to use Baryon Acoustic Oscillations (BAO) to infer $D_A(z)$ from future BAO surveys. As a test case, we consider the proposed Euclid mission investigating the precision can be afforded on $\eta(z)$ from the expected SNeIa and BAO data. We find that the combination of BAO and the local regression coupled allows to reduce the errors on $\eta_a = d\eta/dz|_{z = 0}$ by a factor two if one $\eta_0 = \eta(z = 0) = 1$ is forced and future data are used. On the other hand, although the statistical error on $\eta_0$ is not significantly reduced, the constraints on this quantity will be nevertheless ameliorated thanks to the reduce impact of systematics.

A covariant formulation of a theory with a massive graviton and no negative energy state has been recently proposed as an alternative to the usual General Relativity framework. For a spatially flat homogenous and isotropic universe, the theory introduces modified Friedmann equations where the standard matter term is supplemented by four effective fluids mimicking dust, cosmological constant, quintessence and stiff matter, respectively. We test the viability of this massive gravity formulation by contrasting its theoretical prediction to the Hubble diagram as traced by Type Ia Supernovae (SNeIa) and Gamma Ray Bursts (GRBs), the $H(z)$ measurements from passively evolving galaxies, Baryon Acoustic Oscillations (BAOs) from galaxy surveys and the distance priors from the Cosmic Microwave Background Radiation (CMBR) anisotropy spectrum. It turns out that the model is indeed able to very well fit this large dataset thus offering a viable alternative to the usual dark energy framework. We finally set stringent constraints on its parameters also narrowing down the allowed range for the graviton mass.

Generalised Teleparallel gravity, also referred to as f(T) gravity, has been recently proposed as an extended theory of gravitation able to give rise to an accelerated expansion in a matter only universe. The cosmic speed up is driven by an effective torsion fluid whose equation of state depend on the f(T) function entering the modified gravity Lagrangian. We focus on two particular choices for f(T) which share the nice property to emulate a phantom divide crossing as suggested by some recent data. We check their viability contrasting the predicted background dynamics to the Hubble diagram as traced by both Type Ia Supernovae (SNeIa) and Gamma Ray Bursts (GRBs), the measurement of the rate expansion H(z), the Baryon Acoustic Oscillations (BAOs) at different redshifts, and the Cosmic Microwave Background Radiation (CMBR) distance priors. Both f(T) models turn out to be in very good agreement with this large dataset so that we also investigate whether it is possible to discriminate among them relying on the different growth factors.

According to the now strongly supported concordance $\Lambda$CDM model, galaxies may be grossly described as a luminous component embedded in a dark matter halo. The density profile of these mass dominating haloes may be determined by N - body simulations which mimic the evolution of the tiny initial density perturbations during the process leading to the structures we observe today. Unfortunately, when the effect of baryons is taken into account, the situation gets much more complicated due to the difficulties in simulating their physics. As a consequence, a definitive prediction of how dark matter haloes should presently look like is still missing. We revisit here this issue from an observational point of view devoting our attention to dwarf galaxies. Being likely dark matter dominated, these systems are ideal candidates to investigate the present day halo density profiles and check whether dark matter related quantities correlate with the stellar ones or the environment. By fitting a large sample of well measured rotation curves, we infer constraints on both halo structural parameters (such as the logarithmic slope of the density profile and its concentration) and derived quantities (e.g., the mass fraction and the Newtonian acceleration) which could then be used to constrain galaxy formation scenarios. Moreover, we investigate whether the halo properties correlates with the environment the galaxy lives in thus offering a new tool to deepen our understanding of galaxy formation.

We focus on a modified version of Horava - Lifschitz theory and, in particular, we consider the impact of its weak - field static spherically symmetric limit on the galaxy dynamics. In a previous paper, we used the modified gravitational potential obtained in this theory to evaluate the Milky Way rotation curve using a spheroidal truncated power - law bulge and a double exponential disc as the only sources of the gravitational field and showed that the modified rotation curved is not in agreement with the data. Making a step forward, we here include also the contribution from a dark matter halo in order to see whether this helps fitting the rotation curve data. As a test case, we consider a sample of spiral galaxies with smooth baryon matter distribution and well measured circular velocity profiles. It turns out that, although a marginal agreement with the data can be found, this can only be obtained if the dark matter halo has an unrealistically small virial mass and incredibly large concentration. Such results can be interpreted as a strong evidence against the reliability of the gravitational potential obtained in the modified version of Horava -Lifschitz theory that we consider.

The observed accelerated cosmic expansion can be a signature of fourth\,-\,order gravity theories, where the acceleration of the Universe is a consequence of departures from Einstein General Relativity, rather than the sign of the existence of a fluid with negative pressure. In the fourth\,-\,order gravity theories, the gravity Lagrangian is described by an analytic function $f(R)$ of the scalar curvature $R$ subject to the demanding conditions that no detectable deviations from standard GR is observed on the Solar System scale. Here we consider two classes of $f(R)$ theories able to pass Solar System tests and investigate their viability on cosmological scales. To this end, we fit the theories to a large dataset including the combined Hubble diagram of Type Ia Supernovae and Gamma Ray Bursts, the Hubble parameter $H(z)$ data from passively evolving red galaxies, Baryon Acoustic Oscillations extracted from the seventh data release of the Sloan Digital Sky Survey (SDSS) and the distance priors from the Wilkinson Microwave Anisotropy Probe seven years (WMAP7) data. We find that both classes of $f(R)$ fit very well this large dataset with the present\,-\,day values of the matter density, Hubble constant and deceleration parameter in agreement with previous estimates; however, the strong degeneracy among the $f(R)$ parameters prevents us from strongly constraining their values. We also derive the growth factor $g = d\ln{\delta}/d\ln{a}$, with $\delta = \delta \rho_M/\rho_M$ the matter density perturbation, and show that it can still be well approximated by $g(z) \propto \Omega_M(z)^{\gamma}$. We finally constrain $\gamma$ (on some representative scales) and investigate its redshift dependence to see whether future data can discriminate between these classes of $f(R)$ theories and standard dark energy models.

Euclid is a space-based survey mission from the European Space Agency designed to understand the origin of the Universe's accelerating expansion. It will use cosmological probes to investigate the nature of dark energy, dark matter and gravity by tracking their observational signatures on the geometry of the universe and on the cosmic history of structure formation. The mission is optimised for two independent primary cosmological probes: Weak gravitational Lensing (WL) and Baryonic Acoustic Oscillations (BAO). The Euclid payload consists of a 1.2 m Korsch telescope designed to provide a large field of view. It carries two instruments with a common field-of-view of ~0.54 deg2: the visual imager (VIS) and the near infrared instrument (NISP) which contains a slitless spectrometer and a three bands photometer. The Euclid wide survey will cover 15,000 deg2 of the extragalactic sky and is complemented by two 20 deg2 deep fields. For WL, Euclid measures the shapes of 30-40 resolved galaxies per arcmin2 in one broad visible R+I+Z band (550-920 nm). The photometric redshifts for these galaxies reach a precision of dz/(1+z) < 0.05. They are derived from three additional Euclid NIR bands (Y, J, H in the range 0.92-2.0 micron), complemented by ground based photometry in visible bands derived from public data or through engaged collaborations. The BAO are determined from a spectroscopic survey with a redshift accuracy dz/(1+z) =0.001. The slitless spectrometer, with spectral resolution ~250, predominantly detects Ha emission line galaxies. Euclid is a Medium Class mission of the ESA Cosmic Vision 2015-2025 programme, with a foreseen launch date in 2019. This report (also known as the Euclid Red Book) describes the outcome of the Phase A study.

LambdaCDM, for the currently preferred cosmological density Omega_0 and cosmological constant Omega_Lambda, predicts that the Universe expansion decelerates from early times to redshift z~0.9 and accelerates at later times. On the contrary, the cosmological model based on conformal gravity predicts that the cosmic expansion has always been accelerating. To distinguish between these two very different cosmologies, we resort to gamma-ray bursts (GRBs), which have been suggested to probe the Universe expansion history at z>1, where identified type Ia supernovae (SNe) are rare. We use the full Bayesian approach to infer the cosmological parameters and the additional parameters required to describe the GRB data available in the literature. For the first time, we use GRBs as cosmological probes without any prior information from other data. In addition, when we combine the GRB samples with SNe, our approach neatly avoids all the inconsistencies of most numerous previous methods that are plagued by the so-called circularity problem. In fact, when analyzed properly, current data are consistent with distance moduli of GRBs and SNe that can respectively be, in a variant of conformal gravity, ~15 and ~3 magnitudes fainter than in LambdaCDM. Our results indicate that the currently available SN and GRB samples are accommodated equally well by both LambdaCDM and conformal gravity and do not exclude a continuous accelerated expansion. We conclude that GRBs are currently far from being effective cosmological probes, as they are unable to distinguish between these two very different expansion histories.

Being based on the only assumption that the universe is homogenous and isotropic on large scales, cosmography is an ideal tool to investigate the cosmic expansion history in a almost model-independent way. Fitting the data on the luminosity distance and Baryon Acoustic Oscillations allows to determine the confidence ranges for the cosmographic parameters hence giving some quantitative constraints that a whatever theory has to fulfill. As an application, we consider here the case of teleparallel gravity (TEGR) also referred to as f(T)-gravity. To this end, we first work out analytical expressions to express the present day values of f(T)-derivatives as a function of the cosmographic parameters which hold under quite general and physically motivated conditions. We then use the constraints coming from cosmography to find out the confidence ranges for f(T)-derivatives up to the fifth order and show how these can be used to check the viability of given TEGR models without the need to explicitly solve the second order dynamic equations.

We present the cosmic shear signal predicted by two viable cosmological models in the framework of modified-action f(R) theories. We use f(R) models where the current accelerated expansion of the Universe is a direct consequence of the modified gravitational Lagrangian rather than Dark Energy (DE), either in the form of vacuum energy/cosmological constant or of a dynamical scalar field (e.g. quintessence). We choose Starobinsky's (St) and Hu & Sawicki's (HS) f(R) models, which are carefully designed to pass the Solar System gravity tests. In order to further support - or rule out - f(R) theories as alternative candidates to the DE hypothesis, we exploit the power of weak gravitational lensing, specifically of cosmic shear. We calculate the tomographic shear matrix as it would be measured by the upcoming ESA Cosmic Vision Euclid satellite. We find that in the St model the cosmic shear signal is almost completely degenerate with LCDM, but it is easily distinguishable in the HS model. Moreover, we compute the corresponding Fisher matrix for both the St and HS models, thus obtaining forecasts for their cosmological parameters. Finally, we show that the Bayes factor for cosmic shear will definitely favour the HS model over LCDM if Euclid measures a value larger than ~0.02 for the extra HS parameter n_HS.

We analyze the stellar mass-to-light ratio (M/L) gradients in a large sample of local galaxies taken from the Sloan Digital Sky Survey, spanning a wide range of stellar masses and morphological types. As suggested by the well known relationship between M/L ratios and colors, we show that M/L gradients are strongly correlated with colour gradients, which we trace to the effects of age variations. Stellar M/L gradients generally follow patterns of variation with stellar mass and galaxy type that were previous found for colour and metallicty gradients. In late-type galaxies M/L gradients are negative, steepening with increasing mass. In early-type galaxies M/L gradients are shallower while presenting a two-fold trend: they decrease with mass up to a characteristic mass of \M* \sim 10^10.3 M_sun and increase at larger masses. We compare our findings with other analyses and discuss some implications for galaxy formation and for dark matter estimates.

Thanks to their enormous energy release, Gamma Rays Bursts (GRBs) have recently attracted a lot of interest to probe the Hubble diagram (HD) deep into the matter dominated era and hence complement Type Ia Supernovae (SNeIa). We consider here three different calibration methods based on the use of a fiducial LCDM model, on cosmographic parameters and on the local regression on SNeIa to calibrate the scaling relations proposed as an equivalent to the Phillips law to standardize GRBs finding any significant dependence. We then investigate the evolution of these parameters with the redshift to obtain any statistical improvement. Under this assumption, we then consider possible systematics effects on the HDs introduced by the calibration method, the averaging procedure and the homogeneity of the sample arguing against any significant bias.

Scaling relations among dark matter (DM) and stellar quantities are a valuable tool to constrain formation scenarios and the evolution of galactic structures. However, most of the DM properties are actually not directly measured, but derived through model dependent mass mapping procedures. It is therefore crucial to adopt theoretically and observationally well founded models. We use here an updated version of the secondary infall model (SIM) to predict the halo density profile, taking into account the effects of angular momentum, dissipative friction and baryons collapse. The resulting family of halo profiles depends on one parameter only, the virial mass, and nicely fits the projected mass and aperture velocity dispersion of a sample of intermediate redshift lens galaxies. We derive DM related quantities (namely the column density and the Newtonian acceleration) and investigate their correlations with stellar mass, luminosity, effective radius and virial mass.

Thanks to their enormous energy release which allows to detect them up to very high redshift, Gamma Rays Bursts (GRBs) have recently attracted a lot of interest to probe the Hubble diagram (HD) deep into the matter dominated era and hence complement Type Ia Supernoave (SNeIa). However, lacking a local GRBs sample, calibrating the scaling relations proposed as an equivalent to the Phillips law to standardize GRBs is not an easy task because of the need to estimate the GRBs luminosity distance in a model independent way. We consider here three different calibration methods based on the use of a fiducial $\Lambda$CDM model, on cosmographic parameters and on the local regression on SNeIa. We find that the calibration coefficients and the intrinsic scatter do not significantly depend on the adopted calibration procedure. We then investigate the evolution of these parameters with the redshift finding no statistically motivated improvement in the likelihood so that the no evolution assumption is actually a well founded working hypothesis. Under this assumption, we then consider possible systematics effects on the HDs introduced by the calibration method, the averaging procedure and the homogeneity of the sample arguing against any significant bias. We nevertheless stress that a larger GRBs sample with smaller uncertainties is needed to definitely conclude that the different systematics considered here have indeed a negligible impact on the HDs thus strengthening the use of GRBs as cosmological tools.

Considering galaxies as self - gravitating systems of many collisionless particles allows to use methods of statistical mechanics inferring the distribution function of these stellar systems. Actually, the long range nature of the gravitational force contrasts with the underlying assumptions of Boltzmann statistics where the interactions among particles are assumed to be short ranged. A particular generalization of the classical Boltzmann formalism is available within the nonextensive context of Tsallis q -statistics, subject to non -additivity of the entropies of sub - systems. Assuming stationarity and isotropy in the velocity space, it is possible solving the generalized collsionless Boltzmann equation to derive the galaxy distribution function and density profile. We present a particular set of nonextensive models and investigate their dynamical and observable properties. As a test of the viability of this generalized context, we fit the rotation curve of M33 showing that the proposed approach leads to dark matter haloes in excellent agreement with the observed data.

A viable alternative to the dark energy as a solution of the cosmic speed up problem is represented by Extended Theories of Gravity. Should this be indeed the case, there will be an impact not only on cosmological scales, but also at any scale, from the Solar System to extragalactic ones. In particular, the gravitational potential can be different from the Newtonian one commonly adopted when computing the circular velocity fitted to spiral galaxies rotation curves. Phenomenologically modelling the modified point mass potential as the sum of a Newtonian and a Yukawa like correction, we simulate observed rotation curves for a spiral galaxy described as the sum of an exponential disc and a NFW dark matter halo. We then fit these curves assuming parameterized halo models (either with an inner cusp or a core) and using the Newtonian potential to estimate the theoretical rotation curve. Such a study allows us to investigate the bias on the disc and halo model parameters induced by the systematic error induced by forcing the gravity theory to be Newtonian when it is not. As a general result, we find that both the halo scale length and virial mass are significantly overestimated, while the dark matter mass fraction within the disc optical radius is typically underestimated. Moreover, should the Yukawa scale length be smaller than the disc half mass radius, then the logarithmic slope of the halo density profile would turn out to be shallower than the NFW one. Finally, cored models are able to fit quite well the simulated rotation curves, provided the disc mass is biased high in agreement with the results in literature, favoring cored haloes and maximal discs. Such results make us argue that the cusp/core controversy could actually be the outcome of an incorrect assumption about which theory of gravity must actually be used in computing the theoretical circular velocity.

A phenomenological extension of the well-known brane-world cosmology of Dvali, Gabadadze and Porrati (eDGP) has recently been proposed. In this model, a cosmological-constant-like term is explicitly present as a non-vanishing tension sigma on the brane, and an extra parameter alpha tunes the cross-over scale r_c, the scale at which higher dimensional gravity effects become non negligible. Since the Hubble parameter in this cosmology reproduces the same LCDM expansion history, we study how upcoming weak lensing surveys, such as Euclid and DES (Dark Energy Survey), can confirm or rule out this class of models. We perform Markov Chain Monte Carlo simulations to determine the parameters of the model, using Type Ia Supernov\ae, H(z) data, Gamma Ray Bursts and Baryon Acoustic Oscillations. We also fit the power spectrum of the temperature anisotropies of the Cosmic Microwave Background to obtain the correct normalisation for the density perturbation power spectrum. Then, we compute the matter and the cosmic shear power spectra, both in the linear and non-linear regimes. The latter is calculated with the two different approaches of Hu and Sawicki (2007) (HS) and Khoury and Wyman (2009) (KW). With the eDGP parameters coming from the Markov Chains, KW reproduces the LCDM matter power spectrum at both linear and non-linear scales and the LCDM and eDGP shear signals are degenerate. This result does not hold with the HS prescription: Euclid can distinguish the eDGP model from LCDM because their expected power spectra roughly differ by the 3sigma uncertainty in the angular scale range 700<l<3000; on the contrary, the two models differ at most by the 1sigma uncertainty over the range 500<l<3000 in the DES experiment and they are virtually indistinguishable.

Gamma Ray Bursts (GRBs) are the most energetic sources in the universe and among the farthest known astrophysical sources. These features make them appealing candidates as standard candles for cosmological applications so that studying the physical mechanisms for the origin of the emission and correlations among their observable properties is an interesting task. We consider here the luminosity L*X - break time Ta* (hereafter LT) correlation and investigate whether there are systematics induced by selection effects or redshift dependent calibra- tion. We perform this analysis both for the full sample of 77 GRBs with known redshift and for the subsample of GRBs having canonical X-ray light curves, hereafter called U0095 sample. We do not find any systematic bias thus con- firming the existence of physical GRB subclasses revealed by tight correlations of their afterglow properties. Furthermore, we study the possibility of applying the LT correlation as a redshift estimator both for the full distribution and for the canonical lightcurves. The large uncertainties and the non negligible intrin- sic scatter make the results not so encouraging, but there are nevertheless some hints motivating a further analysis with an increased U0095 sample.

Modified Newtonian Dynamics (MOND) has been shown to be able to fit spiral galaxy rotation curves as well as giving a theoretical foundation for empirically determined scaling relations, such as the Tully - Fisher law, without the need for a dark matter halo. As a complementary analysis, one should investigate whether MOND can also reproduce the dynamics of early - type galaxies (ETGs) without dark matter. As a first step, we here show that MOND can indeed fit the observed central velocity dispersion $\sigma_0$ of a large sample of ETGs assuming a simple MOND interpolating functions and constant anisotropy. We also show that, under some assumptions on the luminosity dependence of the Sersic n parameter and the stellar M/L ratio, MOND predicts a fundamental plane for ETGs : a log - linear relation among the effective radius $R_{eff}$, $\sigma_0$ and the mean effective intensity $\langle I_e \rangle$. However, we predict a tilt between the observed and the MOND fundamental planes.

In an attempt to look for a viable mechanism leading to a present-day accelerated expansion, we investigate the possibility that the observed cosmic speed up may be recovered in the framework of the Rastall's theory, relying on the non-conservativity of the stress-energy tensor, i.e. $T^{\mu}_{\nu ; \mu} \neq 0$. We derive the modified Friedmann equations and show that they correspond to Cardassian-like equations. We also show that, under suitable assumptions on the equation of state of the matter term sourcing the gravitational field, it is indeed possible to get an accelerated expansion, in agreement with the Hubble diagram of both Type Ia Supernovae (SNeIa) and Gamma Ray Bursts (GRBs). Unfortunately, to achieve such a result one has to postulate a matter density parameter much larger than the typical $\Omega_M \simeq 0.3$ value inferred from cluster gas mass fraction data.

Cosmography is a useful tool to constrain cosmological models, in particular dark energy models. In the case of modified theories of gravity, where the equations of motion are generally quite complicated, cosmography can contribute to select realistic models without imposing arbitrary choices a priori. Indeed, its reliability is based on the assumptions that the universe is homogeneous and isotropic on large scale and luminosity distance can be "tracked" by the derivative series of the scale factor a(t). We apply this approach to induced gravity brane-world models where an f(R)-term is present in the brane effective action. The virtue of the model is to self-accelerate the normal and healthy DGP branch once the f(R)-term deviates from the Hilbert-Einstein action. We show that the model, coming from a fundamental theory, is consistent with the LCDM scenario at low redshift. We finally estimate the cosmographic parameters fitting the Union2 Type Ia Supernovae (SNeIa) dataset and the distance priors from Baryon Acoustic Oscillations (BAO) and then provide constraints on the present day values of f(R) and its second and third derivatives.

We discuss the colour, age and metallicity gradients in a wide sample of local SDSS early- and late-type galaxies. From the fitting of stellar population models we find that metallicity is the main driver of colour gradients and the age in the central regions is a dominant parameter which rules the scatter in both metallicity and age gradients. We find a consistency with independent observations and a set of simulations. From the comparison with simulations and theoretical considerations we are able to depict a general picture of a formation scenario.

Gamma Ray Bursts (GRB) observed up to redshifts $z>8$ are fascinating objects to study due to their still unexplained relativistic outburst mechanisms and a possible use to test cosmological models. Our analysis of 77 GRB afterglows with known redshifts revealed a physical subsample of long GRBs with canonical {\it plateau breaking to power-law} light curves with a significant {\it luminosity $L^*_X$ - break time $T^*_a$} correlation in the GRB rest frame. This subsample forms approximately the {\it upper envelope} of the studied distribution. We have also found a similar relation for a small sample of GRB afterglows that belong to the intermediate class (IC) between the short and the long ones. It proves that within the full sample of afterglows there exist physical subclasses revealed here by tight correlations of their afterglow properties. The afterglows with regular (`canonical') light curves obey not only a mentioned tight physical scaling, but -- for a given $T^*_a$ -- the more regular progenitor explosions lead to preferentially brighter afterglows.

We investigate scaling relations between the dark matter (DM) halo model parameters for a sample of intermediate redshift early - type galaxies (ETGs) resorting to a combined analysis of Einstein radii and aperture velocity dispersions. Modeling the dark halo with a Navarro - Frenk - White profile and assuming a Salpeter initial mass function (IMF) to estimate stellar masses, we find that the column density ${\cal{S}}$ and the Newtonian acceleration within the halo characteristic radius $r_s$ and effective radius $R_{eff}$ are not universal quantities, but correlate with the luminosity $L_V$, the stellar mass $M_{\star}$ and the halo mass $M_{200}$, contrary to recent claims in the literature. We finally discuss a tight correlation among the DM mass $M_{DM}(R_{eff})$ within the effective radius $R_{eff}$, the stellar mass $M_{\star}(R_{eff})$ and $R_{eff}$ itself. The slopes of the scaling relations discussed here strongly depend, however, on the DM halo model and the IMF adopted so that these ingredients have to be better constrained in order to draw definitive conclusions on the DM scaling relations for ETGs.

We present the Hubble diagram (HD) of 66 Gamma Ray Bursts (GRBs) derived using only data from their X - ray afterglow lightcurve. To this end, we use the recently updated L_X - T_a correlation between the break time T_a and the X - ray luminosity L_X measured at T_a calibrated from a sample of Swift GRBs with lightcurves well fitted by the Willingale et al. (2007) model. We then investigate the use of this HD to constrain cosmological parameters when used alone or in combination with other data showing that the use of GRBs leads to constraints in agreement with previous results in literature. We finally argue that a larger sample of high luminosity GRBs can provide a valuable information in the search for the correct cosmological model.

We analyze the color gradients (CGs) of ~50000 nearby SDSS galaxies. From synthetic spectral models based on a simplified star formation recipe, we derive the mean spectral properties, and explain the observed radial trends of the color as gradients of the stellar population age and metallicity (Z). The most massive ETGs (M_* > 10^{11} Msun) have shallow CGs in correspondence of shallow (negative) Z gradients. In the stellar mass range 10^(10.3-10.5) < M_* < 10^(11) Msun, the Z gradients reach their minimum of ~ -0.5 dex^{-1}. At M_* ~ 10^{10.3-10.5} Msun, color and Z gradient slopes suddenly change. They turn out to anti-correlate with the mass, becoming highly positive at the very low masses. We have also found that age gradients anti-correlate with Z gradients, as predicted by hierarchical cosmological simulations for ETGs. On the other side, LTGs have gradients which systematically decrease with mass (and are always more negative than in ETGs), consistently with the expectation from gas infall and SN feedback scenarios. Z is found to be the main driver of the trend of color gradients, especially for LTGs, but age gradients are not negligible and seem to play a significant role too. We have been able to highlight that older galaxies have systematically shallower age and Z gradients than younger ones. Our results for high-mass galaxies are in perfect agreement with predictions based on the merging scenario, while the evolution of LTGs and younger and less massive ETGs seems to be mainly driven by infall and SN feedback. (Abridged)

The Horava - Lifshitz (HL) theory has recently attracted a lot of interest as a viable solution to some quantum gravity related problems and the presence of an effective cosmological constant able to drive the cosmic speed up. We show here that, in the weak field limit, the HL proposal leads to a modification of the gravitational potential because of two additive terms (scaling respectively as $r^2$ and $r^{-4}$) to the Newtonian $1/r$ potential. We then derive a general expression to compute the rotation curve of an extended system under the assumption that the mass density only depends on the cylindrical coordinates $(R, z)$ showing that the HL modification induces a dependence of the circular velocity on the mass function which is a new feature of the theory. As a first exploratory analysis, we then try fitting the Milky Way rotation curve using its visible components only in order to see whether the HL modified potential can be an alternative to the dark matter framework. This turns out not to be the case so that we argue that dark matter is still needed, but the amount of dark matter and the dark halo density profile have to be revised according to the new HL potential.

Generalized Einstein - Aether vector field models have been shown to provide, in the weak field regime, modifications to gravity which can be reconciled with the successfull MOND proposal. Very little is known, however, on the function F(K) defining the vector field Lagrangian so that an analysis of the viability of such theories at the cosmological scales has never been performed. As a first step along this route, we rely on the relation between F(K) and the MOND interpolating function $\mu(a/a_0)$ to assign the vector field Lagrangian thus obtaining what we refer to as "MONDian vector models". Since they are able by construction to recover the MOND successes on galaxy scales, we investigate whether they can also drive the observed accelerated expansion by fitting the models to the Type Ia Supernovae data. Should be this the case, we have a unified framework where both dark energy and dark matter can be seen as different manifestations of a single vector field. It turns out that both MONDian vector models are able to well fit the low redshift data on Type Ia Supernovae, while some tension could be present in the high z regime.

Fourth - order gravity theories have received much interest in recent years thanks to their ability to provide an accelerated cosmic expansion in a matter only universe. In these theories, the Lagrangian density of the gravitational field has the form R + f(R), and the explicit choice of the arbitrary function f(R) must meet the local tests of gravity and the constraints from the primordial abundance of the light elements. Two popular classes of f(R) models, which are expected to fulfill all the above requirements, have recently been proposed. However, neither of these models has ever been quantitatively tested against the available astrophysical data. Here, by combining Type Ia Supernovae and Gamma Ray Bursts, we investigate the ability of these models to reproduce the observed Hubble diagram over the redshift range (0, 7). We find that both models fit very well this dataset with the present day values of the matter density and deceleration parameters which agree with previous estimates. However, the strong degeneracy among the f(R) parameters prevents us from putting strong constraints on the values of these parameters; nevertheless, we can identify the regions of the parameter space that should, in principle, be carefully explored with future data and dynamical probes in order to discriminate among f(R) theories and standard dark energy models.

Gamma ray bursts (GRBs) have recently attracted much attention as a possible way to extend the Hubble diagram to very high redshift. To this aim, the luminosity (or isotropic emitted energy) of a GRB at redshift z must be evaluated from a correlation with a distance independent quantity so that one can then solve for the luminosity distance D_L(z) and hence the distance modulus mu(z). Averaging over five different two parameters correlations and using a fiducial cosmological model to calibrate them, Schaefer (2007) has compiled a sample of 69 GRBs with measured mu(z) which has since then been widely used to constrain cosmological parameters. We update here that sample by many aspects. First, we add a recently found correlation for the X - ray afterglow and use a Bayesian inspired fitting method to calibrate the different GRBs correlations known insofar assuming a fiducial LCDM model in agreement with the recent WMAP5 data. Averaging over six correlations, we end with a new GRBs Hubble diagram comprising 83 objects. We also extensively explore the impact of varying the fiducial cosmological model considering how the estimated mu(z) change as a function of the $(\Omega_M, w_0, w_a)$ parameters of the Chevallier - Polarski - Linder phenomenological dark energy equation of state. In order to avoid the need of assuming an {\it a priori} cosmological model, we present a new calibration procedure based on a model independent local regression estimate of mu(z) using the Union SNeIa sample to calibrate the GRBs correlations. This finally gives us a GRBs Hubble diagram made out of 69 GRBs whose estimated distance modulus mu(z) is almost independent on the underlying cosmological model.

We explore the predictions of a class of dark energy models, quinstant dark energy, concerning the structure formation in the Universe, both in the linear and non-linear regimes. Quinstant dark energy is considered to be formed by quintessence and a negative cosmological constant. We conclude that these models give good predictions for structure formation in the linear regime, but fail to do so in the non-linear one, for redshifts larger than one.

The dark matter content of early,- type galaxies (ETGs) is a hotly debated topic with contrasting results arguing in favour or against the presence of significant dark mass within the effective radius and the change with luminosity and mass. In order to address this question, we investigate here the global mass - to - light ratio $\Upsilon(r) = M(r)/L(r)$ of a sample of 21 lenses observed within the Sloan Lens ACS (SLACS) survey. We follow the usual approach of modeling the galaxy as a two component systems, but we use a phenomenological ansatz for $\Upsilon(r)$, proposed by some of us in Tortora et al. (2007), able to smoothly interpolate between constant $M/L$ models and a wide class of dark matter haloes. The resulting galaxy model is then fitted to the data on the Einstein radius and velocity dispersion. Our phenomenological model turns out to be in well agreement with the data suggesting the presence of massive dark matter haloes in order to explain the lensing and dynamics properties of the SLACS lenses. According to the values of the dark matter mass fraction, we argue that the halo may play a significant role in the inner regions probed by the data, but such a conclusion strongly depends on the adopted initial mass function of the stellar population. Finally, we find that the dark matter mass fraction within $R_{eff}$ scales with both the total luminosity and stellar mass in such a way that more luminous (and hence more massive) galaxies have a larger dark matter content.

Early - type galaxies (ETGs) define a narrow strip in the size - mass plane because of the observed correlation between the effective radius $R_{eff}$ and the total stellar mass $M_{\star}$. When expressed in logarithmic units, a linear relation, $\log{R_{eff}} \propto \gamma \log{M_{\star}}$, is indeed observationally found, but the slope $\gamma$ deviates from the canonical $\gamma = 1/2$ value which can be naively predicted for a spherically symmetric isotropic system. We propose here that a transfer of angular momentum to the stellar component induces an anisotropy in the velocity space thus leading to a modified distribution function (DF). Assuming an Osipkov - Merritt like anisotropic DF, we derive an analytic relation between the slope $\gamma$ of the size - mass relation and the slope $\alpha$ of the angular momentum term in the DF. With this simple model, we are then able to recover the observed $\gamma$ value provided $\alpha$ is suitably set. It turns out that an anisotropy profile which is tangential inside $\sim 0.6 r_a$ and radial outside, with $r_a$ the anisotropy radius, is able to reproduce the observed size - mass relation observed for massive ($M_{\star} \ge 2 \times 10^{10} \ h^{-1} {\rm M}_{\odot}$) elliptical galaxies.

Previous work in the literature had built a formalism for spatially averaged equations for the scale factor, giving rise to an averaged Raychaudhuri equation and averaged Hamiltonian constraint, which involve a backreaction source term. The present paper extends these equations to include models with variable Newton parameter and variable cosmological term, motivated by the nonperturbative renormalization program for quantum gravity based upon the Einstein-Hilbert action. We focus on the Brans-Dicke form of the renormalization-group improved action functional. The coupling between backreaction and spatially averaged three-dimensional scalar curvature is found to survive, and a variable-G cosmic quintet is found to emerge. Interestingly, under suitable assumptions, an approximate solution can be found where the early universe tends to a FLRW model, while keeping track of the original inhomogeneities through three effective fluids. The resulting qualitative picture is that of a universe consisting of baryons only, while inhomogeneities average out to give rise to the full dark-side phenomenology.

Gamma ray bursts (GRBs) have recently attracted much attention as a possible way to extend the Hubble diagram to very high redshift. However, the large scatter in their intrinsic properties prevents directly using them as distance indicator so that the hunt is open for a relation involving an observable property to standardize GRBs in the same way as the Phillips law makes it possible to use Type Ia Supernovae (SNeIa) as standardizable candles. We use here the data on the X - ray decay curve and spectral index of a sample of GRBs observed with the Swift satellite. These data are used as input to a Bayesian statistical analysis looking for a correlation between the X - ray luminosity L_X(T_a) and the time constant T_a of the afterglow curve. We find a linear relation between \log{[L_X(T_a)]} and \log{[T_a/(1+z)]} with an intrinsic scatter sigma_{int} = 0.33 comparable to previously reported relations. Remarkably, both the slope and the intrinsic scatter are almost independent on the matter density Omega_M and the constant equation of state w of the dark energy component thus suggesting that the circularity problem is alleviated for the $L_X - T_a$ relation.

It is nowadays accepted that the universe is undergoing a phase of accelerated expansion as tested by the Hubble diagram of Type Ia Supernovae (SNeIa) and several LSS observations. Future SNeIa surveys and other probes will make it possible to better characterize the dynamical state of the universe renewing the interest in cosmography which allows a model independent analysis of the distance - redshift relation. On the other hand, fourth order theories of gravity, also referred to as $f(R)$ gravity, have attracted a lot of interest since they could be able to explain the accelerated expansion without any dark energy. We show here how it is possible to relate the cosmographic parameters (namely the deceleration $q_0$, the jerk $j_0$, the snap $s_0$ and the lerk $l_0$ parameters) to the present day values of $f(R)$ and its derivatives $f^{(n)}(R) = d^nf/dR^n$ (with $n = 1, 2, 3$) thus offering a new tool to constrain such higher order models. Our analysis thus offers the possibility to relate the model independent results coming from cosmography to the theoretically motivated assumptions of $f(R)$ cosmology.

In order to solve the problem of eternal acceleration, a model has been recently proposed including both a negative cosmological constant $\Lambda$ and a scalar field evolving under the action of an exponential potential. We further explore this model by contrasting it against the Hubble diagram of Type Ia supernovae, the gas mass fraction in galaxy clusters and the acoustic peak and shift parameters. It turns out that the model is able to fit quite well this large dataset so that we conclude that a negative $\Lambda$ is indeed allowed and could represent a viable mechanism to halt eternal acceleration. In order to avoid problems with theoretical motivations for both a negative $\Lambda$ term and the scalar field, we reconstruct the gravity Lagrangian $f(R)$ of a fourth order theory of gravity predicting the same dynamics (scale factor and Hubble parameter) as the starting model. We thus end up with a $f(R)$ theory able to both fit the data and solve the problem of eternal acceleration without the need of unusual negative $\Lambda$ and ad hoc scalar fields.

Early-type galaxies (ETGs) are known to possess a number of quite useful scaling relations among their photometric and/or kinematical quantities. We propose a unified picture reducing both the fundamental plane and the photometric plane to suitable projections of a single relation. Modelling the ETG as a two component system, made out of a luminous Sersic profile and a NFW dark halo, and applying the virial theorem, we are able to express the velocity dispersion \sigma_0 as a function of the effective intensity I_e, the effective radius R_e and the Sersic index n. In a log-log plot, this relation reduces to a hyperplane (i.e., a plane in a four dimensional configuration space) which we dubbed the {\it Sersic Virial Hyperplane} (SVH). The tilt of the SVH can be fully explained in terms of a power-law scaling of the stellar (rather than the global) mass-to-light ratio Upsilon with the total luminosity L_T, while the scatter is determined by those on the c-M_{vir} relation between the concentration c and the virial mass M_{vir} of the dark halo component. In order to test whether such the observed SVH is consistent with the theoretical assumptions, we perform a detailed simulation of a large ETGs sample reproducing the same photometric properties of a SDSS low redshift ETGs catalog. It turns out that the simulated SVH is fully consistent with the observed one thus validating our theoretical approach.

We investigate the possibility that the observed flatness of the rotation curves of spiral galaxies is not an evidence for the existence of dark matter haloes, but rather a signal of the breakdown of General Relativity. To this aim, we consider power - law fourth order theories of gravity obtained by replacing the scalar curvature $R$ with $f(R) = f_0 R^n$ in the gravity Lagrangian. We show that, in the low energy limit, the gravitational potential generated by a pointlike source may be written as $\Phi(r) \propto r^{-1} \left [ 1 + (r/r_c)^{\beta} \right ]$ with $\beta$ a function of the slope $n$ of the gravity Lagrangian and $r_c$ a scalelength depending on the gravitating system properties. In order to apply the model to realistic systems, we compute the modified potential and the rotation curve for spherically symmetric and for thin disk mass distributions. It turns out that the potential is still asymptotically decreasing, but the corrected rotation curve, although not flat, is higher than the Newtonian one thus offering the possibility to fit rotation curves without dark matter. To test the viability of the model, we consider a sample of 15 low surface brightness (LSB) galaxies with combined HI and H$\alpha$ measurements of the rotation curve extending in the putative dark matter dominated region. We find a very good agreement between the theoretical rotation curve and the data using only stellar disk and interstellar gas.

Recent observations of the line of sight velocity profile of elliptical galaxies have furnished controversial results with some works favouring the presence of a large amount of dark matter in the outer regions and others arguing in favour of no dark matter at all. In order to shed new light on this controversy, we propose here a new phenomenological description of the total mass profile of galaxies. Under the hypothesis of spherical symmetry, we assume a double power-law expression for the global M/L ratio Upsilon(r)= Upsilon_0(r/r_0) ^{alpha}(1+r/r_0)^{beta}. In particular, Upsilon propto r^{alpha} for r/r_0<<1 so that alpha=0 mimics a constant M/L ratio in the inner regions, while, for (r/r_0>>1), Upsilon propto r^{alpha+beta} thus showing that models with alpha+beta=0 have an asymptotically constant M/L ratio. A wide range of possibilities is obtained by varying the slope parameters in the range we determine on the basis of physical considerations. Choosing a general expression for the luminosity density profile j(r), we work out an effective galaxy model that accounts for all the phenomenology observed in real elliptical galaxies. We derive the main dynamics and lensing properties of such an effective model. We analyze a general class of models, able to take into account different dynamical trends. We are able to obtain analytical expressions for the main dynamical and lensing quantities. We show that constraining the values of alpha+beta makes it possible to analyze the problem of the dark matter in elliptical galaxies. Indeed, positive values of alpha+beta would be a strong evidence for dark matter. Finally we indicate possible future approaches in order to face the observational data, in particular using velocity dispersion profiles and lensed quasar events.

The inflessence model has recently been proposed in an attempt to explain both early inflation and present day accelerated expansion within a single mechanism. The model has been successfully tested against the Hubble diagram of Type Ia Supernovae, the shift parameter, and the acoustic peak parameter. As a further mandatory test, we investigate here structure formation in the inflessence model determining the evolution of matter density contrast $\delta \equiv \delta \rho_M/\rho_M$ in the linear regime. We compare the growth factor $D(a) \equiv \delta/a$ and the growth index $f(z) \equiv d\ln{\delta}/d\ln{a}$ to these same quantities for the successful concordance $\Lambda$CDM model with a particular emphasis on the role of the inflessence parameters $(\gamma, z_Q)$. We also evaluate the anisotropy spectrum of the cosmic microwave background radiation (CMBR) to check whether the inflessence model may be in agreement with the observations. We find that, for large values of $(\gamma, z_Q)$, structure formation proceeds in a similar way to that in the $\Lambda$CDM scenario, and it is also possible to nicely fit the CMBR spectrum.

Gravitational lensing is investigated in the weak field limit of fourth order gravity in which the Lagrangian of the gravitational field is modified by replacing the Ricci scalar curvature R with an analytical expression $f(R)$. Considering the case of a pointlike lens, we study the behaviour of the deflection angle in the case of power law Lagrangians, i.e. with f(R) = f_0 R^n. In order to investigate possible detectable signatures, the position of the Einstein ring and the solutions of the lens equation are evaluated considering the change with respect to the standard case. Effects on the amplification of the images and the Paczynski curve in microlensing experiments are also estimated.

Astrophysical observations are pointing out huge amounts of dark matter and dark energy needed to explain the observed large scale structures and cosmic accelerating expansion. Up to now, no experimental evidence has been found, at fundamental level, to explain such mysterious components. The problem could be completely reversed considering dark matter and dark energy as shortcomings of General Relativity and claiming for the correct theory of gravity as that derived by matching the largest number of observational data. As a result, accelerating behavior of cosmic fluid and rotation curves of spiral galaxies are reproduced by means of curvature effects.

We investigate the effects of viscosity terms depending on the Hubble parameter and its derivatives in the dark energy equation of state. Such terms are possible if dark energy is a fictitious fluid originating from corrections to the Einstein general relativity as is the case for some braneworld inspired models or fourth order gravity. We consider two classes of models whose singularities in the early and late time universe have been studied by testing the models against the dimensionless coordinate distance to Type Ia Supernovae and radio-galaxies also including priors on the shift and the acoustic peak parameters. It turns out that both models are able to explain the observed cosmic speed up without the need of phantom (w < -1) dark energy.

We investigate the possibility that part of the dark matter is not made out of the usual cold dark matter (CDM) dustlike particles, but is under the form of a fluid of strings with barotropic factor $w_s = -1/3$ of cosmic origin. To this aim, we split the dark matter density parameter in two terms and investigate the dynamics of a spatially flat universe filled with baryons, CDM, fluid of strings and dark energy, modeling this latter as a cosmological constant or a negative pressure fluid with a constant equation of state $w < 0$. To test the viability of the models and to constrain their parameters, we use the Type Ia Supernovae Hubble diagram and the data on the gas mass fraction in galaxy clusters. We also discuss the weak field limit of a model comprising a significant fraction of dark matter in the form of a fluid of strings and show that this mechanism makes it possible to reduce the need for the elusive and up to now undetected CDM. We finally find that a model comprising both a cosmological constant and a fluid of strings fits very well the data and eliminates the need of phantom dark energy thus representing a viable candidate to alleviate some of the problems plaguing the dark side of the universe.

Abandoning the perfect fluid hypothesis, we investigate here the possibility that the dark energy equation of state (EoS) $w$ is a nonlinear function of the energy density $\rho$. To this aim, we consider four different EoS describing classical fluids near thermodynamical critical points and discuss the main features of cosmological models made out of dust matter and a dark energy term with the given EoS. Each model is tested against the data on the dimensionless coordinate distance to Type Ia Supernovae and radio galaxies, the shift and the acoustic peak parameters and the positions of the first three peaks in the anisotropy spectrum of the comic microwave background radation. We propose a possible interpretation of each model in the framework of scalar field quintessence determining the shape of the self interaction potential $V(\phi)$ that gives rise to each one of the considered thermodynamical EoS. As a general result, we demonstrate that replacing the perfect fluid EoS with more generar expressions gives both the possibility of successfully solving the problem of cosmic acceleration escaping the resort to phantom models.

It has been recently shown that, under some general conditions, it is always possible to find a fourth order gravity theory capable of reproducing the same dynamics of a given dark energy model. Here, we discuss this approach for a dark energy model with a scalar field evolving under the action of an exponential potential. In absence of matter, such a potential can be recovered from a fourth order theory via a conformal transformation. Including the matter term, the function f(R) entering the generalized gravity Lagrangian can be reconstructed according to the dark energy model.

We present a phenomenologically motivated model which is able to give rise both to an inflationary epoch and to the present day cosmic acceleration. We introduce an approach where the energy density depends on the scale factor $a$ in such a way that a smooth transition from the inflation to the matter dominated era and then to the nowadays accelerated expansion is achieved. We use the dimensionless coordinate distance $y(z)$ to the Gold Type Ia Supernovae dataset and to a sample comprising 20 radio galaxies, the shift parameter ${\cal{R}}$ and the acoustic peak parameter ${\cal{A}}$ to test whether the model is in agreement with observations and to constrain its main parameters. As an independent cross check, we also compare model predictions with the lookback time to galaxy clusters and the age of the universe. Although phenomenologically inspired, the model may be theoretically motivated either resorting to a scenario with a scalar and phantom fields (eventually interacting) or as an effective description arising from higher order gravity theories.

We discuss an alternative approach to quintessence modifying the usual equation of state of the cosmological fluid in order to see if going further than the approximation of perfect fluid allows to better reproduce the available data. We consider a cosmological model comprising only two fluids, namely baryons (modelled as dust) and dark matter with a Van der Waals equation of state. First, the general features of the model are presented and then the evolution of the energy density, the Hubble parameter and the scale factor are determined showing that it is possible to obtain accelerated expansion choosing suitably the model parameters. We use the estimated age of the universe and the data on the dimensionless coordinate distances to Type Ia supernovae and distant radio galaxies to see whether Van der Waals quintessence is viable to explain dark energy and to constrain its parameters.

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

Higher order theories of gravity have recently attracted a lot of interest as alternative candidates to explain the observed cosmic acceleration without the need of introducing any scalar field. A critical ingredient is the choice of the function f(R) of the Ricci scalar curvature entering the gravity Lagrangian and determining the dynamics of the universe. We describe an efficient procedure to reconstruct f(R) from the Hubble parameter $H$ depending on the redshift z. Using the metric formulation of f(R) theories, we derive a third order linear differential equation for f(R(z)) which can be numerically solved after setting the boundary conditions on the basis of physical considerations. Since H(z) can be reconstructed from the astrophysical data, the method we present makes it possible to determine, in principle, what is the f(R) theory which best reproduces the observed cosmological dynamics. Moreover, the method allows to reconcile dark energy models with f(R) theories finding out what is the expression of f(R) which leads to the same H(z) of the given quintessence model. As interesting examples, we consider "quiessence" (dark energy with constant equation of state) and the Chaplygin gas.

We present a new family of spherically symmetric models for the luminous components of elliptical and spiral galaxies and their dark matter haloes. Our starting point is a general expression for the logarithmic slope $\alpha(r) = d\log{\rho}/d\log{r}$ from which most of the cuspy models yet available in literature may be derived. We then dedicate our attention to a particular set of models whose logarithmic slope is a power law function of the radius $r$ investigating in detail their dynamics assuming isotropy in the velocity space. While the basic properties (such as the density profile and the gravitational potential) may be expressed analytically, both the distribution function and the observable quantities (surface brightness and line of sight velocity dispersion) have to be evaluated numerically. We also consider the extension to anisotropic models trying two different parameterization. Since the model recently proposed by Navarro et al. (2004) as the best fit to their sample of numerically simulated haloes belongs to the family we present here, analytical approximations are given for the most useful quantities.

Higher order curvature gravity has recently received a lot of attention due to the fact that it gives rise to cosmological models which seem capable of solving dark energy and quintessence issues without using "ad hoc" scalar fields. Such an approach is naturally related to fundamental theories of quantum gravity which predict higher order terms for loop expansions of quantum fields in curved spacetimes. In this framework, we obtain a class of cosmological solutions which are fitted against cosmological data. We reproduce reliable models able to fit high redshift supernovae and WMAP observations. The age of the universe and other cosmological parameters are recovered in this context. Furthermore, in the weak field limit, we obtain gravitational potentials which differ from the Newtonian one because of repulsive corrections increasing with distance. We evaluate the rotation curve of our Galaxy and compare it with the observed data in order to test the viability of these theories and to estimate the scale-length of the correction. It is remarkable that the Milky Way rotation curve is well fitted without the need of any dark matter halo and similar results hold also for other galaxies.

An impressive amount of different astrophysical data converges towards the picture of a spatially flat universe undergoing a today phase of accelerated expansion. The nature of the dark energy dominating the energy content of the universe is still unknown and a lot of different scenarios are viable candidates to explain cosmic acceleration. Most of the methods employed to test these cosmological models are essentially based on distance measurements to a particular class of objects. A different method, based on the lookback time to galaxy clusters and the age of the universe, is used here. In particular, we constrain the characterizing parameters of three classes of dark energy cosmological models to see whether they are in agreement with this kind of data, based on time measurements rather than distance observations.

We investigate the viability of f(R) theories in the framework of the Palatini approach as solutions to the problem of the observed accelerated expansion of the universe. Two physically motivated popular choices for f(R) are considered: power law, f(R) = \beta R^n, and logarithmic, f(R) = \alpha \ln{R}. Under the Palatini approach, both Lagrangians give rise to cosmological models comprising only standard matter and undergoing a present phase of accelerated expansion. We use the Hubble diagram of type Ia Supernovae and the data on the gas mass fraction in relaxed galaxy clusters to see whether these models are able to reproduce what is observed and to constrain their parameters. It turns out that they are indeed able to fit the data with values of the Hubble constant and of the matter density parameter in agreement with some model independent estimates, but the today deceleration parameter is higher than what is measured in the concordance LambdaCDM model.

The structure of the dark matter and the thermodynamical status of the hot gas in galaxy clusters is an interesting and widely discussed topic in modern astrophysics. Recently, Rasia et al. (2004) have proposed a new dynamical model for the mass density profile of clusters of galaxies as a result of a set of high resolution hydrodynamical simulations of structure formation. We investigate the lensing properties of this model evaluating the deflection angle, the lensing potential and the amplification of the images. We reserve particular attention to the structure and position of the critical curves in order to see whether this model is able to produce radial and tangential arcs. To this aim, we also investigate the effect of taking into account the brightest cluster galaxy in the lensing potential and the deviations from spherical symmetry mimicked by an external shear. We also analyze the implication of the gas density and temperature profiles of the Rasia et al.(2004) model on the properties of the X - ray emission and the comptonization parameter that determines the CMBR temperature decrement due to the Sunyaev - Zel'dovich effect.

We present a numerical method to estimate the lensing parameters and the Hubble constant H_0 from quadruply imaged gravitational lens systems. The lens galaxy is modeled using both separable deflection potentials and constant mass-to-light ratio profiles, while possible external perturbations have been taken into account introducing an external shear. The model parameters are recovered inverting the lens and the time delay ratio equations and imposing a set of physically motivated selection criteria. We investigate correlations among the model parameters and the Hubble constant. Finally, we apply the codes to the real lensed quasars PG 1115+080 and RX J0911+0551, and combine the results from these two systems to get H_0 = 56 +/- 23 km/(s Mpc). In addition, we are able to fit to the single systems a general elliptical potential with a non fixed angular part, and then we model the two lens systems with the same potential and a shared H_0: in this last case we estimate H_0=49_(-11)^(+6) km/(s Mpc).

A phenomenological approach is proposed to the problem of universe accelerated expansion and of the dark energy nature. A general class of models is introduced whose energy density depends on the redshift $z$ in such a way that a smooth transition among the three main phases of the universe evolution (radiation era, matter domination, asymptotical de Sitter state) is naturally achieved. We use the estimated age of the universe, the Hubble diagram of Type Ia Supernovae and the angular size - redshift relation for compact and ultracompact radio structures to test whether the model is in agreement with astrophysical observation and to constrain its main parameters. Although phenomenologically motivated, the model may be straightforwardly interpreted as a two fluids scenario in which the quintessence is generated by a suitably chosen scalar field potential. On the other hand, the same model may also be read in the context of unified dark energy models or in the framework of modified Friedmann equation theories.

It is well known that the surface brightness of elliptical galaxies and of bulges of spiral galaxies is best fitted by the Sersic r^{1/n} profile. It is thus interesting to explore the lensing properties of the Sersic model because of its wide range of applicability. To this aim, we evaluate the lensing potential, the deflection angle, the time delay between the images and the a mplification for a circularly symmetric lens whose surface density is described by the Sersic profile. We estimate the same quantities also adding a shear term to the deflection potential in order to study (at the lowest order) the impact of deviations from radial symmetry or the contribution from other nearby lenses. Moreover, we investigate the systematic errors due to the use of a de Vaucouleurs profile instead of the correct Sersic model also taking into account the presence of a dark halo.

Quintessence issues can be achieved by taking into account higher order curvature invariants into the effective action of gravitational field. Such an approach is naturally related to fundamental theories of quantum gravity which predict higher order terms in loop expansion of quantum fields in curved space-times. In this framework, we obtain a class of cosmological solutions which are fitted against cosmological data. We reproduce encouraging results able to fit high redshift supernovae and WMAP observations. The age of the universe and other cosmological parameters are discussed in this context.

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.

It has been suggested that a monitoring program employing the pixel lensing method to search for microlensing events towards galactic globular clusters may increase the statistics and discriminate among different halo models. Stimulated by this proposal, we evaluate an upper limit to the pixel lensing event rate for such a survey. Four different dark halo models have been considered changing both the flattening and the slope of the mass density profile. The lenses mass function has been modelled as a homogenous power - law for $\mu \in (\mu_l, \mu_u)$ and both the mass limits and the slope of the mass function have been varied to investigate their effect on the rate. The target globular clusters have been selected in order to minimize the disk contribution to the event rate. We find that a pixel lensing survey towards globular clusters is unable to discriminate among different halo models since the number of detectable events is too small to allow any reliable statistical analysis.

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

We present a new method to estimate the Hubble constant H_0 from the measured time delays in quadruply imaged gravitational lens systems. We show how it is possible to get an estimate of H_0 without the need to completely reconstruct the lensing potential thus avoiding any a priori hypotheses on the expression of the galaxy lens model. Our method only needs to assume that the lens potential may be expressed as r^{\alpha} F(\theta), whatever the shape function F(\theta) is, and it is thus able to fully explore the degeneracy in the mass models taking also into account the presence of an external shear. We test the method on simulated cases and show that it does work well in recovering the correct value of the slope \alpha of the radial profile and of the Hubble constant H_0. Then, we apply the same method to the real quadruple lenses PG1115+080 and B1422+231 obtaining H_0 = 58_{-15}^{+17} km/s/Mpc (68% CL).

We develop a semi - analytical method to reconstruct the lensing potential in quadruply imaged gravitational lens systems. Assuming that the potential belongs to a broad class of boxy non - elliptical models, we show how it is possible to write down a system of equations which can be numerically solved to recover the potential parameters directly from image positions and using physical constraints. We also describe a code developed to search for solutions of the system previously found and test it on simulated cases. Finally, we apply the method to the quadruple lens PG1115+080 which allows us to get also an estimate of the Hubble constant H_0 from the measured time delay as H_0 = 56_{-11}^{+17} km/s/Mpc.

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 dark halo lens mass function (MF) for a wide class of spheroidal non singular isothermal models comparing observed and observable microlensing quantities for MACHO observations towards LMC and taking into account the detection efficiency. We evaluate the microlensing observable quantities, i.e. observable optical depth, number of events and mean duration, for models with homogenous power - law MF changing the upper and lower mass limits and the flattening of the dark halo. By applying the simple technique of the inverse problem method we are then able to get some interesting constraints on the slope $\alpha$ of the MF and on the dark halo mass fraction f made out by MACHOs consistently with previous results.

SLOTT-AGAPE (Systematic Lensing Observation at Toppo Telescope - Andromeda Gravitational Amplification Pixel Lensing Experiment) is a new collaboration project among international partners from England, France, Germany, Italy and Switzerland that intends to perform microlensing observation by using M31 as target. The MACHOs search is made thanks to the pixel lensing technique.

We present in this paper a simple method to obtain informations on galaxy models parameters using the measured value of the microlensing optical depth. Assuming a 100 percent MACHO's dark halo, we ask the predicted optical depth for a given model to be the same as the observed one, in a given direction. Writing the optical depth in terms of the given halo model parameters and inverting this relation with respect to one of them, it is possible to get information on it, fixing under reasonable hypothesis the other parameters. This is what we call the "inverse problem in microlensing. We apply this technique to the class of power-law models with flat rotation curves, determining the range for the core radius R_c compatible with the constraints on the halo flattening q and the measures of the optical depth towards LMC. Next, we apply the same method to a simple triaxial model, evaluating the axial ratios.