In this work, we consider an extension of symmetric teleparallel gravity, namely, $f(Q)$ gravity, where the fundamental block to describe spacetime is the nonmetricity, $Q$. Within this formulation of gravitation, we perform an observational analysis of several modified $f(Q)$ models using the redshift approach, where the $f(Q)$ Lagrangian is reformulated as an explicit function of the redshift, $f(z)$. Various different polynomial parameterizations of $f(z)$ are proposed, including new terms which would allow for deviations from the $\Lambda$CDM model. Given a variety of observational probes, such as the expansion rate data from early-type galaxies, Type Ia Supernovae, Quasars, Gamma Ray Bursts, Baryon Acoustic Oscillations data and Cosmic Microwave Background distance priors, we have checked the validity of these models at the background level in order to verify if this new formalism provides us with plausible alternative models to explain the late time acceleration of the universe. Indeed, this novel approach provides a different perspective on the formulation of observationally reliable alternative models of gravity.

The clustering amplitude of 7143 clusters from the Sloan Digital Sky Survey (SDSS) is found to increase with cluster mass, closely agreeing with the Gaussian random field hypothesis for structure formation. The amplitude of the observed cluster correlation exceeds the predictions from pure cold dark matter (CDM) simulation by $\simeq 6\%$ for the standard Planck-based values of the cosmological parameters. We show that this excess can be naturally accounted for by free streaming of light neutrinos, which opposes gravitational growth, so clusters formed at fixed mass are fewer and hence more biased than for a pure CDM density field. An enhancement of the cluster bias by 7\% matches the observations, corresponding to a total neutrino mass, $m_{\nu} = 0.119 \pm 0.034$ eV at 67\% confidence level, for the standard relic neutrino density. If ongoing laboratory experiments favor a normal neutrino mass hierarchy then we may infer a somewhat larger total mass than the minimum oscillation based value, $\sum m_{\nu} \simeq 0.056eV$, with 90\% confidence. Much higher precision can be achieved by applying our method to a larger sample of more distant clusters with weak lensing derived masses.

In this work we study in detail a phenomenological generalization of the Chaplygin cosmological model, which we call as \textit{umami} Chaplygin model. We consider three different cosmological background scenarios in which our fluid can play three different roles: only as a dark energy component; as a dark matter and dark energy component; and as a dark plus baryonic matter and dark energy component. With such analysis we explore the possibility to unify the dark fluids under one single component within the context of General Relativity. We test this hypothesis against the main available data related to the cosmological background, namely: Type Ia supernovae; Baryon Acoustic Oscillations; Quasars; Gamma Ray Bursts; Hubble data from cosmic chronometers; and Cosmic Microwave Background. We eventually compare the statistically efficiency and reliability of our model to describe observational data with respect to the standard $\Lambda$CDM model by means of the Bayesian Evidence. Final results point to a positive (albeit not strong) evidence in favor of a possible unification of dark energy and dark matter with a the \textit{umami} fluid.

The study of the dynamics of a two-body system in modified gravity constitutes a more complex problem than in Newtonian gravity. Numerical methods are typically needed to solve the equations of geodesics. Despite the complexity of the problem, the study of a two-body system in $f(R)$ gravity leads to a new exciting perspective hinting the right strategy to adopt in order to probe modified gravity. Our results point out some differences between the {\em semiclassical} (Newtonian) approach, and the {\em relativistic} (geodesic) one thus suggesting that the latter represents the best strategy for future tests of modified theories of gravity. { Finally, we have also highlighted the capability of forthcoming observations to serve as smoking gun of modified gravity revealing a departure from GR or further reducing the parameter space of $f(R)$ gravity}. \keywords{$f(R)$ gravity \and binary system \and geodesics \and precession

Alternative theories of gravity may serve to overcame several shortcomings of the standard cosmological model but, in their weak field limit, General Relativity must be recovered so as to match the tight constraints at the Solar System scale. Therefore, testing such alternative models at scales of stellar systems could give a unique opportunity to confirm or rule them out. One of the most straightforward modifications is represented by analytical $f(R)$-gravity models that introduce a Yukawa-like modification to the Newtonian potential thus modifying the dynamics of particles. Using the geodesics equations, we have illustrated the amplitude of these modifications. First, we have integrated numerically the equations of motion showing the orbital precession of a particle around a massive object. Second, we have computed an analytic expression for the periastron advance of systems having their semi-major axis much shorter than the Yukawa-scale length. Finally, we have extended our results to the case of a binary system composed of two massive objects. Our analysis provides a powerful tool to obtain constraints on the underlying theory of gravity using current and forthcoming datasets.

It is a very well established matter nowadays that many modified gravity models can offer a sound alternative to General Relativity for the description of the accelerated expansion of the universe. But it is also equally well known that no clear and sharp discrimination between any alternative theory and the classical one has been found so far. In this work, we attempt at formulating a different approach starting from the general class of $f(R)$ theories as test probes: we try to reformulate $f(R)$ Lagrangian terms as explicit functions of the redshift, i.e., as $f(z)$. In this context, the $f(R)$ setting to the consensus cosmological model, the $\Lambda$CDM model, can be written as a polynomial including just a constant and a third-order term. Starting from this result, we propose various different polynomial parameterizations $f(z)$, including new terms which would allow for deviations from $\Lambda$CDM, and we thoroughly compare them with observational data. While on the one hand we have found no statistically preference for our proposals (even if some of them are as good as $\Lambda$CDM by using Bayesian Evidence comparison), we think that our novel approach could provide a different perspective for the development of new and observationally reliable alternative models of gravity.

We have selected a sample of eleven massive clusters of galaxies observed by the Hubble Space Telescope in order to study the impact of the dynamical state on the IntraCluster Light (ICL) fraction, the ratio of total integrated ICL to the total galaxy member light. With the exception of the Bullet cluster, the sample is drawn from the Cluster Lensing and Supernova Survey and the Frontier Fields program, containing five relaxed and six merging clusters. The ICL fraction is calculated in three optical filters using the CHEFs IntraCluster Light Estimator, a robust and accurate algorithm free of a priori assumptions. We find that the ICL fraction in the three bands is, on average, higher for the merging clusters, ranging between $\sim7-23\%$, compared with the $\sim 2-11\%$ found for the relaxed systems. We observe a nearly constant value (within the error bars) in the ICL fraction of the regular clusters at the three wavelengths considered, which would indicate that the colors of the ICL and the cluster galaxies are, on average, coincident and, thus, their stellar populations. However, we find a higher ICL fraction in the F606W filter for the merging clusters, consistent with an excess of lower-metallicity/younger stars in the ICL, which could have migrated violently from the outskirts of the infalling galaxies during the merger event.

The cosmological redshift drift could lead to the next step in high-precision cosmic geometric observations, becoming a direct and irrefutable test for cosmic acceleration. In order to test the viability and possible properties of this effect, also called Sandage-Loeb (SL) test, we generate a model independent mock data set so as to compare its constraining power with that of the future mock data sets of Type Ia Supernovae (SNe) and Baryon Acoustic Oscillations (BAO). The performance of those data sets is analyzed by testing several cosmological models with the Markov chain Monte Carlo (MCMC) method, both independently and combining all data sets. Final results show that, in general, SL data sets allow for remarkable constraints on the matter density parameter today $\Omega_m$ on every tested model, showing also a great complementarity with SNe and BAO data regarding dark energy (DE) parameters.

Light Axionic Dark Matter, motivated by string theory, is increasingly favored for the "no-WIMP era". Galaxy formation is suppressed below a Jeans scale, of $\simeq 10^8 M_\odot$ by setting the axion mass to, $m_B \sim 10^{-22}$eV, and the large dark cores of dwarf galaxies are explained as solitons on the de-Broglie scale. This is persuasive, but detection of the inherent scalar field oscillation at the Compton frequency, $\omega_B= (2.5{\rm \, months})^{-1}(m_B/10^{-22}eV)$, would be definitive. By evolving the coupled Schr\"odinger-Poisson equation for a Bose-Einstein condensate, we predict the dark matter is fully modulated by de-Broglie interference, with a dense soliton core of size $\simeq 150pc$, at the Galactic center. The oscillating field pressure induces General Relativistic time dilation in proportion to the local dark matter density and pulsars within this dense core have detectably large timing residuals, of $\simeq 400nsec/(m_B/10^{-22}eV)$. This is encouraging as many new pulsars should be discovered near the Galactic center with planned radio surveys. More generally, over the whole Galaxy, differences in dark matter density between pairs of pulsars imprints a pairwise Galactocentric signature that can be distinguished from an isotropic gravitational wave background.

We present accurate photometric redshifts for galaxies observed by the Cluster Lensing and Supernova survey with Hubble (CLASH). CLASH observed 25 massive galaxy cluster cores with the Hubble Space Telescope in 16 filters spanning 0.2 - 1.7 $\mu$m. Photometry in such crowded fields is challenging. Compared to our previously released catalogs, we make several improvements to the photometry, including smaller apertures, ICL subtraction, PSF matching, and empirically measured uncertainties. We further improve the Bayesian Photometric Redshift (BPZ) estimates by adding a redder elliptical template and by inflating the photometric uncertainties of the brightest galaxies. The resulting photometric redshift accuracies are dz/(1+z) $\sim$ 0.8\%, 1.0\%, and 2.0\% for galaxies with I-band F814W AB magnitudes $<$ 18, 20, and 23, respectively. These results are consistent with our expectations. They improve on our previously reported accuracies by a factor of 4 at the bright end and a factor of 2 at the faint end. Our new catalog includes 1257 spectroscopic redshifts, including 382 confirmed cluster members. We also provide stellar mass estimates. Finally, we include lensing magnification estimates of background galaxies based on our public lens models. Our new catalog of all 25 CLASH clusters is available via MAST. The analysis techniques developed here will be useful in other surveys of crowded fields, including the Frontier Fields and surveys carried out with J-PAS and JWST.

We construct a large, redshift complete sample of distant galaxy clusters by correlating Sloan Digital Sky Survey (SDSS) Data Release 12 (DR12) redshifts with clusters identified with the red-sequence Matched-filter Probabilistic Percolation (redMaPPer) algorithm. Our spectroscopic completeness is 97% for ~ 7,000 clusters within the redMaPPer selection limit, $z \leqslant$ 0.325, so that our cluster correlation functions are much more precise than earlier work and not suppressed by photometric redshifts. We derive an accurate power-law mass-richness relation from the observed abundance with respect to the mass function from Millennium XXL (MXXL) simulation, adjusted to the Planck weighted cosmology. The number density of clusters is found to decline by 20% over the range 0.1 $< z <$ 0.3, in good agreement with the evolution predicted by MXXL. Our projected three-dimensional correlation function scales with richness, $\lambda$, rising from $r_0=$ 14 h$^{-1}$ Mpc at $\lambda\simeq$ 25, to $r_0=$ 22 h$^{-1}$ Mpc at $\lambda\simeq$ 60, with a gradient that matches MXXL when applying our mass-richness relation, whereas the observed amplitude of the correlation function at $\left<z\right>=$ 0.24 exceeds the MXXL prediction by 20% at the $\simeq$ 2.5$\sigma$ level. This tension cannot be blamed on spurious, randomly located clusters as this would reduce the correlation amplitude. Full consistency between the correlation function and the abundances is achievable for the pre-Planck values of $\sigma_8=$ 0.9, $\Omega_m=$ 0.25, and h = 0.73, matching the improved distance ladder estimate of the Hubble constant.

Using Planck satellite data, we construct SZ gas pressure profiles for a large, volume-complete sample of optically selected clusters. We have defined a sample of over 8,000 redMaPPer clusters from the Sloan Digital Sky Survey (SDSS), within the volume-complete redshift region 0.100 < z < 0.325, for which we construct Sunyaev-Zel'dovich (SZ) effect maps by stacking Planck data over the full range of richness. Dividing the sample into richness bins we simultaneously solve for the mean cluster mass in each bin together with the corresponding radial pressure profile parameters, employing an MCMC analysis. These profiles are well detected over a much wider range of cluster mass and radius than previous work, showing a clear trend towards larger break radius with increasing cluster mass. Our SZ-based masses fall ~24% below the mass-richness relations from weak lensing, in a similar fashion as the "hydrostatic bias" related with X-ray derived masses. We correct for this bias to derive an optimal mass-richness relation finding a slope 1.22 +/- 0.04 and a pivot mass log(M_500/M_0)= 14.432 +/- 0.041, evaluated at a richness lambda=60. Finally, we derive a tight Y_500-M_500 relation over a wide range of cluster mass, with a power law slope equal to 1.72 +/- 0.07, that agrees well with the independent slope obtained by the Planck team with an SZ-selected cluster sample, but extends to lower masses with higher precision.

In this work we consider a family of cosmological models featuring future singularities. This type of cosmological evolution is typical of dark energy models with an equation of state violating some of the standard energy conditions (e.g. the null energy condition). Such kind of behavior, widely studied in the literature, may arise in cosmologies with phantom fields, theories of modified gravity or models with interacting dark matter/dark energy. We briefly review the physical consequences of these cosmological evolution regarding geodesic completeness and the divergence of tidal forces in order to emphasize under which circumstances the singularities in some cosmological quantities correspond to actual singular spacetimes. We then introduce several phenomenological parameterizations of the Hubble expansion rate to model different singularities existing in the literature and use SN Ia, BAO and H(z) data to constrain how far in the future the singularity needs to be (under some reasonable assumptions on the behaviour of the Hubble factor). We show that, for our family of parameterizations, the lower bound for the singularity time can not be smaller than about 1.2 times the age of the universe, what roughly speaking means about 2.8 Gyrs from the present time.

In [1] a new method to measure the speed of light through Baryon Acoustic Oscillations (BAO) was introduced. Here, we describe in much more detail the theoretical basis of that method, its implementation, and give some newly updated results about its application to the forecast data. In particular, we will show that SKA will be able to detect a 1% variation (if any) in the speed of light at 3$\sigma$ level. Smaller signals will be hardly detected by already-planned future galaxy surveys, but we give indications about what sensitivity requirements should a survey fulfill in order to be successful.

We explore the observational adequacy of a class of Unified Dark Energy/Matter (UDE/M) models with a fast transition. Our constraints are set using a combination of geometric probes, some low redshift ones, and some high redshift ones (CMB related included). The transition is phenomenologically modeled by two different transition functions corresponding to a fast and to an ultra-fast transition respectively. We find that the key parameters governing the transition can be well constrained, and from the statistical point of view it follows that the models cannot be discarded when compared to $\Lambda$CDM. We find the intriguing result that standard/input parameters such as $\Omega_m$ and $\Omega_b$ are far better constrained than in $\Lambda$CDM, and this is the case for the derived/output parameter measuring the deceleration value at present, $q_0$.

We review the appropriateness of using SNIa observations to detect potential signatures of anisotropic expansion in the Universe. We focus on Union2 and SNLS3 SNIa datasets and use the hemispherical comparison method to detect possible anisotropic features. Unlike some previous works where non-diagonal elements of the covariance matrix were neglected, we use the full covariance matrix of the SNIa data, thus obtaining more realistic and not underestimated errors. As a matter of fact, the significance of previously claimed detections of a preferred direction in the Union2 dataset completely disappears once we include the effects of using the full covariance matrix. Moreover, we also find that such a a preferred direction is aligned with the orthogonal direction of the SDSS observational plane and this suggests a clear indication that the SDSS subsample of the Union2 dataset introduces a significant biased, making the detected preferred direction unphysical. We thus find that current SNIa surveys are inappropriate to test anisotropic features due to their highly non-homogeneous angular distribution in the sky. In addition, after removal of the highest inhomogeneous sub-samples, the number of SNIa is too low. Finally, we take advantage of the particular distribution of SNLS SNIa sub-sample in the SNLS3 data set, in which the observations were taken along 4 different directions. We fit each direction independently and find consistent results at the 1$\sigma$ level. Although the likelihoods peak at relatively different values of $\Omega_m$,the low number of data along each direction gives rise to large errors so that the likelihoods are sufficiently broad as to overlap within 1$\sigma$.

In this letter we describe a new method to use Baryon Acoustic Oscillations (BAO) to derive a constraint on the possible variation of the speed of light. The method relies on the fact that there is a simple relation between the angular diameter distance $(D_{A})$ maximum and the Hubble function $(H)$ evaluated at the same maximum-condition redshift, which includes speed of light $c$. We note the close analogy of the BAO probe with a laboratory experiment: here we have $D_{A}$ which plays the role of a standard (cosmological) ruler, and $H^{-1}$, with the dimension of time, as a (cosmological) clock. We evaluate if current or future missions such as Euclid can be sensitive enough to detect any variation of $c$.

In this letter we propose a test to detect the linearity of the dark energy equation of state, and apply it to two different Type Ia Supernova (SN Ia) data sets, Union2.1 and SNLS3. We find that: a. current SN Ia data are well described by a dark energy equation of state linear in the cosmic scale factor a, at least up to a redshift z = 1, independent of the pivot points chosen for the linear relation; b. there is no significant evidence of any deviation from linearity. This apparent linearity may reflect the limit of dark energy information extractable from current SN Ia data.

A clean measurement of the evolution of the galaxy cluster mass function can significantly improve our understanding of cosmology from the rapid growth of cluster masses below z < 0.5. Here we examine the consistency of cluster catalogues selected from the SDSS by applying two independent gravity-based methods using all available spectroscopic redshifts from the DR10 release. First, we detect a gravitational redshift related signal for 20,119 and 13,128 clusters with spectroscopic redshifts contained in the GMBCG and redMaPPer catalogues, respectively, at a level of $\sim - 10$ km s$^{-1}$. This we show is consistent with the magnitude expected using the richness-mass relations provided by the literature and after applying recently clarified relativistic and flux bias corrections. This signal is also consistent with the richest clusters in the larger catalogue of Wen et al. (2012), corresponding to $M_{200m} \gtrsim 2 \times 10^{14}\,\mathrm{M}_\odot\,h^{-1}$, however we find no significant detection of gravitational redshift signal for less riched clusters, which may be related to bulk motions from substructure and spurious cluster detections. Second, we find all three catalogues generate mass-dependent levels of lensing magnification bias, which enhances the mean redshift of flux-selected background galaxies from the BOSS survey. The magnitude of this lensing effect is generally consistent with the corresponding richness-mass relations advocated for the surveys. We conclude that all catalogues comprise a high proportion of reliable clusters, and that the GMBCG and redMaPPer cluster finder algorithms favor more relaxed clusters with a meaningful gravitational redshift signal, as anticipated by the red-sequence colour selection of the GMBCG and redMaPPer samples.

A fast transition between a standard matter-like era and a late $\Lambda$CDM-like epoch generated by a single Unified Dark Matter component can explain the observed acceleration of the Universe. UDM models with a fast transition should be clearly distinguishable from $\Lambda$CDM (and alternatives) through observations. Here we focus on a particularly simple model and analyse its viability by studying features of the background model and properties of the adiabatic UDM perturbations.

In this work we investigate if a small fraction of quarks and gluons, which escaped hadronization and survived as a uniformly spread perfect fluid, can play the role of both dark matter and dark energy. This fluid, as developed in \citep{Brilenkov}, is characterized by two main parameters: $\beta$, related to the amount of quarks and gluons which act as dark matter; and $\gamma$, acting as the cosmological constant. We explore the feasibility of this model at cosmological scales using data from type Ia Supernovae (SNeIa), Long Gamma-Ray Bursts (LGRB) and direct observational Hubble data. We find that: (i) in general, $\beta$ cannot be constrained by SNeIa data nor by LGRB or H(z) data; (ii) $\gamma$ can be constrained quite well by all three data sets, contributing with $\approx78\%$ to the energy-matter content; (iii) when a strong prior on (only) baryonic matter is assumed, the two parameters of the model are constrained successfully.

In this work we present a nonparametric approach, which works on minimal assumptions, to reconstruct the cosmic expansion of the Universe. We propose to combine a locally weighted scatterplot smoothing method and a simulation-extrapolation method. The first one (Loess) is a nonparametric approach that allows to obtain smoothed curves with no prior knowledge of the functional relationship between variables nor of the cosmological quantities. The second one (Simex) takes into account the effect of measurement errors on a variable via a simulation process. For the reconstructions we use as raw data the Union2.1 Type Ia Supernovae compilation, as well as recent Hubble parameter measurements. This work aims to illustrate the approach, which turns out to be a self-sufficient technique in the sense we do not have to choose anything by hand. We examine the details of the method, among them the amount of observational data needed to perform the locally weighted fit which will define the robustness of our reconstruction. In view of our results, we believe that our proposal offers a promising alternative for reconstructing global trends of cosmological data when there is little intuition on the relationship between the variables and we also think it even presents good prospects to generate reliable mock data points where the original sample is poor.

The Cluster Lensing And Supernovae survey with Hubble (CLASH) is an Hubble Space Telescope (HST) Multi-Cycle Treasury program observing 25 massive galaxy clusters. CLASH observations are carried out in 16 bands from UV to NIR to derive accurate and reliable estimates of photometric redshifts. We present the CLASH photometric redshifts and study the photometric redshift accuracy of the arcs in more detail for the case of MACS1206.2-0847. We use the publicly available Le Phare and BPZ photometric redshift codes on 17 CLASH galaxy clusters. Using Le Phare code for objects with StoN>=10, we reach a precision of 3%(1+z) for the strong lensing arcs, which is reduced to 2.4%(1+z) after removing outliers. For galaxies in the cluster field the corresponding values are 4%(1+z) and 3%(1+z). Using mock galaxy catalogues, we show that 3%(1+z) precision is what one would expect from the CLASH photometry when taking into account extinction from dust, emission lines and the finite range of SEDs included in the photo-z template library. We study photo-z results for different aperture photometry and find that the SExtractor isophotal photometry works best. Le Phare and BPZ give similar photo-z results for the strong lensing arcs as well as galaxies of the cluster field. Results are improved when optimizing the photometric aperture shape showing an optimal aperture size around 1" radius giving results which are equivalent to isophotal photometry. Tailored photometry of the arcs improve the photo-z results.

Cosmography provides a model-independent way to map the expansion history of the Universe. In this paper we simulate a Euclid-like survey and explore cosmographic constraints from future Baryonic Acoustic Oscillations (BAO) observations. We derive general expressions for the BAO transverse and radial modes and discuss the optimal order of the cosmographic expansion that provide reliable cosmological constraints. Through constraints on the deceleration and jerk parameters, we show that future BAO data have the potential to provide a model-independent check of the cosmic acceleration as well as a discrimination between the standard $\Lambda$CDM model and alternative mechanisms of cosmic acceleration.

Aims. We investigate the degree of improvement in dark energy constraints that can be achieved by extending Type Ia Supernova (SN Ia) samples to redshifts z > 1.5 with the Hubble Space Telescope (HST), particularly in the ongoing CANDELS and CLASH multi-cycle treasury programs. Methods. Using the popular CPL parametrization of the dark energy, w = w0 +wa(1-a), we generate mock SN Ia samples that can be projected out to higher redshifts. The synthetic datasets thus generated are fitted to the CPL model, and we evaluate the improvement that a high-z sample can add in terms of ameliorating the statistical and systematic uncertainties on cosmological parameters. Results. In an optimistic but still very achievable scenario, we find that extending the HST sample beyond CANDELS+CLASH to reach a total of 28 SN Ia at z > 1.0 could improve the uncertainty in the wa parameter by up to 21%. The corresponding improvement in the figure of merit (FoM) would be as high as 28%. Finally, we consider the use of high-redshift SN Ia samples to detect non-cosmological evolution in SN Ia luminosities with redshift, finding that such tests could be undertaken by future spacebased infrared surveys using the James Webb Space Telescope (JWST).

In the strong lensing regime non-parametric lens models struggle to achieve sufficient angular resolution for a meaningful derivation of the central cluster mass distribution. The problem lies mainly with cluster members which perturb lensed images and generate additional images, requiring high resolution modeling, even though we mainly wish to understand the relatively smooth cluster component. The required resolution is not achievable because the separation between lensed images is several times larger than the deflection angles by member galaxies, even for the deepest data. Here we bypass this limitation by incorporating a simple physical prior for member galaxies, using their observed positions and their luminosity scaled masses. This galaxy contribution is added to a relatively coarse Gaussian pixel grid for modeling the cluster mass distribution, extending our established WSLAP code (Diego et al. 2007). We test this new code with a simulation based on A1689, using the pixels belonging to multiply-lensed images and the observed member galaxies. Dealing with the cluster members this way leads to convergent solutions, without resorting to regularization, reproducing well the input cluster and substructures. We highlight the ability of this method to recover dark sub-components of the cluster, unrelated to member galaxies. Such anomalies can provide clues to the nature of invisible dark matter, but are hard to discover using parametrized models where substructures are defined by the visible data. With our increased resolution and stability we show, for the first time, that non-parametric models can be made sufficiently precise to locate multiply-lensed systems, thereby achieving fully self-consistent solutions without reliance on input systems from less objective means.

A fast transition between a standard matter-like era and a late $\Lambda$ cold dark matter ($\Lambda$CDM)-like epoch (or more in general, a CDM+DE era), generated by a single Unified Dark Matter component, can provide a new interesting paradigm in the context of general relativity, alternative to $\Lambda$CDM itself or other forms of DE or modified gravity theories invoked to explain the observed acceleration of the Universe. UDM models with a fast transition have interesting features, leading to measurable predictions, thus they should be clearly distinguishable from $\Lambda$CDM (and alternatives) through observations. Here we look at different ways of prescribing phenomenological UDM models with fast transition, then focusing on a particularly simple model. We analyse the viability of this model by studying features of the background model and properties of the adiabatic UDM perturbations, which depend on the effective speed of sound and the functional form of the Jeans scale. As a result, theoretical constraints on the parameters of the models are found that allow for a behaviour compatible with observations.

In this work we set observational constraints of the Superfluid Chaplygin gas model, which gives a unified description of the dark sector of the Universe as a Bose-Einstein condensate (BEC) that behaves as dark energy (DE) while it is in the ground state and as dark matter (DM) when it is in the excited state. We first show and perform the various steps leading to a form of the equations suitable for the observational tests to be carried out. Then, by using a Markov Chain Monte Carlo (MCMC) code, we constrain the model with a sample of cosmology-independent long gamma-ray bursts (LGRBs) calibrated using their Type I Fundamental Plane, as well as the Union2.1 set and observational Hubble parameter data. In this analysis, using our cosmological constraints, we sketch the effective equation of state parameter and deceleration parameter, and we also obtain the redshift of the transition from deceleration to acceleration: $z_t$.

The cosmological model consisting of an electromagnetic Born-Infeld (BI) field coupled to a Robertson-Walker geometry is tested with the standard probes of SNIa, GRBs and direct Hubble parameter. The analysis shows that the inclusion of the nonlinear electromagnetic component does not contribute in a significative way to the observed expansion. The BI electromagnetic matter is considered with an abundance of $\Omega_{BI}$, that our best fit leads to $\Omega_{BI}=0.037$ when tested with SNIa and the Hubble parameter data (0.1<z<1.75); while when tested with GRBs the result is of $\Omega_{BI}=0.304$, which may indicate that this electrodynamics was important at epochs close to the appearance of large structure (redshifts of approximately 7), although this late result has not as much reliability as that corresponding to the first two probes, since we know that the dispersion in GRBs data is still considerable. In view of these results we can rule out the electromagnetic Born-Infeld matter as the origin of the present accelerated expansion, this conclusion concerns exclusively the Born-Infeld theory.

We examine the inner mass distribution of the relaxed galaxy cluster Abell 383 in deep 16-band HST/ACS+WFC3 imaging taken as part of the CLASH multi-cycle treasury program. Our program is designed to study the dark matter distribution in 25 massive clusters, and balances depth with a wide wavelength coverage to better identify lensed systems and generate precise photometric redshifts. This information together with the predictive strength of our strong-lensing analysis method identifies 13 new multiply-lensed images and candidates, so that a total of 27 multiple-images of 9 systems are used to tightly constrain the inner mass profile, $d\log \Sigma/d\log r\simeq -0.6\pm 0.1$ (r<160 kpc). We find consistency with the standard distance-redshift relation for the full range spanned by the lensed images, 1.01<z<6.03, with the higher redshift sources deflected through larger angles as expected. The inner mass profile derived here is consistent with the results of our independent weak-lensing analysis of wide-field Subaru images, with good agreement in the region of overlap. The overall mass profile is well fitted by an NFW profile with M_{vir}=(5.37^{+0.70}_{-0.63}\pm 0.26) x 10^{14}M_{\odot}/h and a relatively high concentration, c_{vir}=8.77^{+0.44}_{-0.42}\pm 0.23, which lies above the standard c-M relation similar to other well-studied clusters. The critical radius of Abell 383 is modest by the standards of other lensing clusters, r_{E}\simeq16\pm2\arcsec (for z_s=2.55), so the relatively large number of lensed images uncovered here with precise photometric redshifts validates our imaging strategy for the CLASH survey. In total we aim to provide similarly high-quality lensing data for 25 clusters, 20 of which are X-ray selected relaxed clusters, enabling a precise determination of the representative mass profile free from lensing bias. (ABRIDGED)

We show that a large-area imaging survey using narrow-band filters could detect quasars in sufficiently high number densities, and with more than sufficient accuracy in their photometric redshifts, to turn them into suitable tracers of large-scale structure. If a narrow-band optical survey can detect objects as faint as i=23, it could reach volumetric number densities as high as 10^{-4} h^3 Mpc^{-3} (comoving) at z~1.5 . Such a catalog would lead to precision measurements of the power spectrum up to z~3-4. We also show that it is possible to employ quasars to measure baryon acoustic oscillations at high redshifts, where the uncertainties from redshift distortions and nonlinearities are much smaller than at z<1. As a concrete example we study the future impact of J-PAS, which is a narrow-band imaging survey in the optical over 1/5 of the unobscured sky with 42 filters of ~100 A full-width at half-maximum. We show that J-PAS will be able to take advantage of the broad emission lines of quasars to deliver excellent photometric redshifts, \sigma_{z}~0.002(1+z), for millions of objects.

In this work we introduce two new polynomial parametrizations of dark energy and explore their correlation properties. The parameters to fit are the equation of state values at z=0 and z=0.5, which have naturally low correlation and have already been shown to improve the popular Chevallier-Polarski-Linder (CPL) parametrization. We test our models with low redshift astronomical probes: type Ia supernovae and baryon acoustic oscillations (BAO), in the form of both current and synthetic data. Specifically, we present simulations of measurements of the radial and transversal BAO scales similar to those expected in a BAO high precision spectroscopic redshift survey similar to EUCLID. According to the Bayesian deviance information criterion (DIC), which penalizes large errors and correlations, we show that our models perform better than the CPL re-parametrization proposed by Wang (in terms of z=0 and z=0.5). This is due to the combination of a lower correlation and smaller relative errors. The same holds for a frequentist perspective: our Figure-of-Merit is larger for our parametrizations.

Model independent reconstructions of dark energy have received some attention. The approach that addresses the reconstruction of the dimensionless coordinate distance and its two first derivatives using a polynomial fit in different redshift windows is well developed \cite{DalyDjorgovski1,DalyDjorgovski2,DalyDjorgovski3}. In this work we offer new insights into the problem by focusing on two types of observational probes: SNeIa and GRBs. Our results allow to highlight some of the intrinsic weaknesses of the method. One of the directions we follow is to consider updated observational samples. Our results indicate than conclusions on the main dark energy features as drawn from this method are intimately related to the features of the samples themselves (which are not quite ideal). This is particularly true of GRBs, which manifest themselves as poor performers in this context. In contrast to original works, we conclude they cannot be used for cosmological purposes, and the state of the art does not allow to regard them on the same quality basis as SNeIa. The next direction we contribute to is the question of how the adjusting of some parameters (window width, overlap, selection criteria) affect the results. We find again there is a considerable sensitivity to these features. Then, we try to establish what is the current redshift range for which one can make solid predictions on dark energy evolution. Finally, we strengthen the former view that this model is modest in the sense it provides only a picture of the global trend. But, on the other hand, we believe it offers an interesting complement to other approaches given that it works on minimal assumptions.

We study the possibility to extract model independent information about the dynamics of the universe by using Cosmography. We intend to explore it systematically, to learn about its limitations and its real possibilities. Here we are sticking to the series expansion approach on which Cosmography is based. We apply it to different data sets: Supernovae Type Ia (SNeIa), Hubble parameter extracted from differential galaxy ages, Gamma Ray Bursts (GRBs) and the Baryon Acoustic Oscillations (BAO) data. We go beyond past results in the literature extending the series expansion up to the fourth order in the scale factor, which implies the analysis of the deceleration, q_{0}, the jerk, j_{0} and the snap, s_{0}. We use the Markov Chain Monte Carlo Method (MCMC) to analyze the data statistically. We also try to relate direct results from Cosmography to dark energy (DE) dynamical models parameterized by the Chevalier-Polarski-Linder (CPL) model, extracting clues about the matter content and the dark energy parameters. The main results are: a) even if relying on a mathematical approximate assumption such as the scale factor series expansion in terms of time, cosmography can be extremely useful in assessing dynamical properties of the Universe; b) the deceleration parameter clearly confirms the present acceleration phase; c) the MCMC method can help giving narrower constraints in parameter estimation, in particular for higher order cosmographic parameters (the jerk and the snap), with respect to the literature; d) both the estimation of the jerk and the DE parameters, reflect the possibility of a deviation from the LCDM cosmological model.

The Cluster Lensing And Supernova survey with Hubble (CLASH) is a 524-orbit multi-cycle treasury program to use the gravitational lensing properties of 25 galaxy clusters to accurately constrain their mass distributions. The survey, described in detail in this paper, will definitively establish the degree of concentration of dark matter in the cluster cores, a key prediction of CDM. The CLASH cluster sample is larger and less biased than current samples of space-based imaging studies of clusters to similar depth, as we have minimized lensing-based selection that favors systems with overly dense cores. Specifically, twenty CLASH clusters are solely X-ray selected. The X-ray selected clusters are massive (kT > 5 keV; 5 - 30 x 10^14 M_solar) and, in most cases, dynamically relaxed. Five additional clusters are included for their lensing strength (Einstein radii > 35 arcsec at z_source = 2) to further quantify the lensing bias on concentration, to yield high resolution dark matter maps, and to optimize the likelihood of finding highly magnified high-redshift (z > 7) galaxies. The high magnification, in some cases, provides angular resolutions unobtainable with any current UVOIR facility and can yield z > 7 candidates bright enough for spectroscopic follow-up. A total of 16 broadband filters, spanning the near-UV to near-IR, are employed for each 20-orbit campaign on each cluster. These data are used to measure precise (sigma_phz < 0.02(1+z)) photometric redshifts for dozens of newly discovered multiply-lensed images per cluster. Observations of each cluster are spread over 8 epochs to enable a search, primarily in the parallel fields, for Type Ia supernovae at z > 1 to improve constraints on the time dependence of the dark energy equation of state and the evolution of such supernovae in an epoch when the universe is matter dominated.

We present a strong-lensing analysis of the galaxy cluster MACS J1206.2-0847 ($z$=0.44) using UV, Optical, and IR, HST/ACS/WFC3 data taken as part of the CLASH multi-cycle treasury program, with VLT/VIMOS spectroscopy for some of the multiply-lensed arcs. The CLASH observations, combined with our mass-model, allow us to identify 47 new multiply-lensed images of 12 distant sources. These images, along with the previously known arc, span the redshift range $1\la z\la5.5$, and thus enable us to derive a detailed mass distribution and to accurately constrain, for the first time, the inner mass-profile of this cluster. We find an inner profile slope of $d\log \Sigma/d\log \theta\simeq -0.55\pm 0.1$ (in the range [1\arcsec, 53\arcsec], or $5\la r \la300$ kpc), as commonly found for relaxed and well-concentrated clusters. Using the many systems uncovered here we derive credible critical curves and Einstein radii for different source redshifts. For a source at $z_{s}\simeq2.5$, the critical curve encloses a large area with an effective Einstein radius of $\theta_{E}=28\pm3\arcsec$, and a projected mass of $1.34\pm0.15\times10^{14} M_{\odot}$. From the current understanding of structure formation in concordance cosmology, these values are relatively high for clusters at $z\sim0.5$, so that detailed studies of the inner mass distribution of clusters such as MACS J1206.2-0847 can provide stringent tests of the $\Lambda$CDM paradigm.

We analyse dark energy models where self-interacting three-forms or phantom fields drive the accelerated expansion of the Universe. The dynamics of such models is often studied by rewriting the cosmological field equations in the form of a system of autonomous differential equations, or simply a dynamical system. Properties of these systems are usually studied via linear stability theory. In situations where this method fails, for instance due to the presence of zero eigenvalues in the Jacobian, centre manifold theory can be applied. We present a concise introduction and show explicitly how to use this theory in two concrete examples.

Using real and synthetic Type Ia SNe (SNeIa) and baryon acoustic oscillations (BAO) data representing current observations forecasts, this paper investigates the tension between those probes in the dark energy equation of state (EoS) reconstruction considering the well known CPL model and Wang's low correlation reformulation. In particular, here we present simulations of BAO data from both the the radial and transverse directions. We also explore the influence of priors on Omega_m and Omega_b on the tension issue, by considering 1-sigma deviations in either one or both of them. Our results indicate that for some priors there is no tension between a single dataset (either SNeIa or BAO) and their combination (SNeIa+BAO). Our criterion to discern the existence of tension (sigma-distance) is also useful to establish which is the dataset with most constraining power; in this respect SNeIa and BAO data switch roles when current and future data are considered, as forecasts predict and spectacular quality improvement on BAO data. We also find that the results on the tension are blind to the way the CPL model is addressed: there is a perfect match between the original formulation and that by the low correlation optimized, but the errors on the parameters are much narrower in all cases of our exhaustive exploration, thus serving the purpose of stressing the convenience of this reparametrization.

We consider the Einstein equations within the DBI scenario for a spatially flat Friedmann-Robertson-Walker (FRW) spacetime without a cosmological constant. We derive the inflationary scenario by applying the symmetry transformations which preserve the form of the Friedmann and conservation equations. These form invariance transformations generate a symmetry group parametrized by the Lorentz factor $\ga$. We explicitly obtain an inflationary scenario by the cooperative effect of adding energy density into the Friedmann equation. For the case of a constant Lorentz factor, and under the slow roll assumption, we find the transformation rules for the scalar and tensor power spectra of perturbations as well as their ratio under the action of the form invariance symmetry group. Within this case and due to its relevance for the inflationary paradigm, we find the general solution of the dynamical equations for a DBI field driven by an exponential potential and show a broad set of inflationary solutions. The general solution can be split into three subsets and all these behave asymptotically as a power law solution at early and at late times.

We study the possibility of detecting oscillating patterns in the equation of state (EoS) of the dark energy using different cosmological datasets. We follow a phenomenological approach and study three different oscillating models for the EoS, one of them periodic and the other two damped (proposed here for the first time). All the models are characterised by the amplitude value, the centre and the frequency of oscillations. In contrast to previous works in the literature, we do not fix the value of the frequency to a fiducial value related to the time extension of chosen datasets, but consider a discrete set of values, so to avoid arbitrariness and try and detect any possible time period in the EoS. We test the models using a recent collection of SNeIa, direct Hubble data and Gamma Ray Bursts data. Main results are: I. even if constraints on the amplitude are not too strong, we detect a trend of it versus the frequency, i.e. decreasing (and even negatives) amplitudes for higher frequencies; II. the centre of oscillation (which corresponds to the present value of the EoS parameter) is very well constrained, phantom behaviour is excluded at $1\sigma$ level and trend which is in agreement with the one for the amplitude appears; III. the frequency is hard to constrain, showing similar statistical validity for all the values of the discrete set chosen, but the best fit of all the scenarios considered is associated with a period which is in the redshift range depicted by our cosmological data. The "best" oscillating models are compared with $\Lambda$CDM using dimensionally consistent a Bayesian approach based information criterion and the conclusion reached is the non existence of significant evidence against dark energy oscillations.

We introduce a new form of coupling between dark energy and dark matter that is quadratic in their energy densities. Then we investigate the background dynamics when dark energy is in the form of exponential quintessence. The three types of quadratic coupling all admit late-time accelerating critical points, but these are not scaling solutions. We also show that two types of coupling allow for a suitable matter era at early times and acceleration at late times, while the third type of coupling does not admit a suitable matter era.

It has been recently shown that the presence of a vector field over cosmological scales could explain the observed accelerated expansion of the universe without introducing neither new scales nor unnatural initial conditions in the early universe, thus avoiding the coincidence problem. Here, we present a detailed analysis of the constraints imposed by SNIa, CMB and BAO data on the vector dark energy model with general spatial curvature. We find that contrary to standard cosmology, CMB data excludes a flat universe for this model and, in fact, predicts a closed geometry for the spatial sections. We see that CMB and SNIa Gold data are perfectly compatible at the 1-sigma level, however SNIa Union dataset exhibits a 3-sigma tension with CMB. The same level of tension is also found between SNIa and BAO measurements.

We propose a model based on a DBI action for the unification of dark matter and dark energy. This is supported by the results of the study of its background behavior at early and late times, and reinforced by the analysis of the evolution of perturbations. We also perform a Bayesian analysis to set observational constraints on the parameters of the model using type Ia SN, CMB shift and BAO data. Finally, to complete the study we investigate its kinematics aspects, such as the effective equation of state parameter, acceleration parameter and transition redshift. Particularizing those parameters for the best fit one appreciates that an effective phantom is preferred.

We introduce the Lagrangian for a multi-scalar field configuration in a $N$-dimensional internal space endowed with a constant metric $Q_{ik}$ and generalize the quintom cosmological scenario. We find the energy momentum tensor of the model and show that the set of dual transformations, that preserve the form of the Einstein equations in the Friedmann-Robertson-Walker (FRW) cosmology, is enlarged. We show that the stability of the power law solutions leads to an exponential potential which is invariant under linear transformations in the internal space. Moreover, we obtain the general exact solution of the Einstein-Klein-Gordon equations for that potential. There exist solutions that cross the phantom divide and solutions that blow up at a finite time, exhibiting a superaccelerated behavior and ending in a big rip. We show that the quintom model with a separable potential can be identified with a mixture of several fluids. This framework includes the $\Lambda$CDM model, a bouncing model, and a setting sourced by a cosmic string network plus a cosmological constant. The we concentrate on the case where the dimension of the internal quintessence sector $N_{q}$ exceeds the dimension of the internal phantom sector $N_{ph}$. For $(N_q,N_{ph})=(2,1)$ the dark energy density derived from the 3-scalar field crosses the phantom divide and its negative component plays the role of the negative part of a classical Dirac Field.

Dark energy and dark matter are the dominant sources in the evolution of the late universe. They are currently only indirectly detected via their gravitational effects, and there could be a coupling between them without violating observational constraints. We investigate the background dynamics when dark energy is modelled as exponential quintessence, and is coupled to dark matter via simple models of energy exchange. We introduce a new form of dark sector coupling, which leads to a more complicated dynamical phase space and has a better physical motivation than previous mathematically similar couplings.

We compare the dark energy model constraints obtained by using recent standard ruler data (Baryon Acoustic Oscillations (BAO) at z=0.2 and z=0.35 and Cosmic Microwave Background (CMB) shift parameters R and l_a) with the corresponding constraints obtained by using recent Type Ia Supernovae (SnIa) standard candle data (ESSENCE+SNLS+HST from Davis et. al.). We find that, even though both classes of data are consistent with LCDM at the 2\sigma level, there is a systematic difference between the two classes of data. In particular, we find that for practically all values of the parameters (\Omega_0m,\Omega_b) in the 2\sigma range of the the 3-year WMAP data (WMAP3) best fit, LCDM is significantly more consistent with the SnIa data than with the CMB+BAO data. For example for (\Omega_0m,\Omega_b)=(0.24,0.042) corresponding to the best fit values of WMAP3, the dark energy equation of state parametrization w(z)=w_0 + w_1 (z/(1+z)) best fit is at a 0.5\sigma distance from LCDM (w_0=-1,w_1=0) using the SnIa data and 1.7\sigma away from LCDM using the CMB+BAO data. There is a similar trend in the earlier data (SNLS vs CMB+BAO at z=0.35). This trend is such that the standard ruler CMB+BAO data show a mild preference for crossing of the phantom divide line w=-1, while the recent SnIa data favor LCDM. Despite of this mild difference in trends, we find no statistically significant evidence for violation of the cosmic distance duality relation \eta \equiv d_L(z)/(d_A(z) (1+z)^2)=1. For example, using a prior of \Omega_0m=0.24, we find \eta=0.95 \pm 0.025 in the redshift range 0<z<2, which is consistent with distance duality at the 2\sigma level.

We investigate the link between the warp function and the potential in DBI cosmologies in connection with the possibility they represent power-law solutions. A prescription is given to take advantage of the known result that given a warp factor there is always a choice of potential resulting in a constant ratio between pressure and energy density. The method is illustrated with examples with interesting models for either the warp factor or the potential. We complete this investigation by giving a recipe to exploit symmetries in order to generate new solutions from existing ones; this method can be applied, for instance, to the power-law cosmologies obtained using our prescription.

This contribution intends to give a pedagogical introduction to the topic of dark energy (the mysterious agent supposed to drive the observed late time acceleration of the Universe) and to various observational tests which require only assumptions on the geometry of the Universe. Those tests are the supernovae luminosity, the CMB shift, the direct Hubble data, and the baryon acoustic oscillations test. An historical overview of Cosmology is followed by some generalities on FRW spacetimes (the best large-scale description of the Universe), and then the test themselves are discussed. A convenient section on statistical inference is included as well.

Recently H(z) data obtained from differential ages of galaxies have been proposed as a new geometrical probe of dark energy. In this paper we use those data, combined with other background tests (CMB shift and SNIa data), to constrain a set of general relativistic dark energy models together with some other models motivated by extra dimensions. Our analysis rests mostly on Bayesian statistics, and we conclude that LCDM is at least substantially favoured, and that braneworld models are less favoured than general relativistic ones.

A new class of braneworld models displaying late-time phantom acceleration without resorting to a phantom fluid is presented. In this scenario expansion is fuelled by dark matter together with some effective dark energy capable of crossing the phantom divide. Unlike a previous proposal of this nature, in these models the effective phantom behaviour remains valid at all redshifts for some choices of the free parameters of the models. The construction is based on the generalised Chaplygin gas, and the cosmological history interpolates between a standard CDM-like behaviour at early times and a de Sitter-like behaviour at late times, so no future singularity is reached.

In this paper we investigate quintom cosmologies with arbitrary potentials from the dynamical systems perspective. The dynamical systems analysis is complete in the sense that it includes the asymptotic regime where both scalar fields diverge, which proves to be particularly relevant in connection with the existence of tracking phases. The results of the present study indicate that the existence of phantom attractors is not generic: for quintom models there may exist either de Sitter attractors associated with the saddle points of the potential, or tracking attractors in the asymptotic regime where the scalar fields diverge.

We investigate a simple braneworld model in which the universe contains only cold dark matter and a cosmological constant, but the effective dark energy is phantom-like because of extra-dimensional gravity effects. Modified gravity screens the cosmological constant Lambda, allowing for a larger Lambda. In practice, observations do not favour any significant screening. We use supernova data, the cosmic microwave background shift parameter, and the baryon oscillation peak in the galaxy distribution to constrain the model. We find the mean value of Omega_m with 68% confidence limits, and an upper limit on Omega_Lambda at the 68% confidence level. The best-fit model is very close to a standard LCDM model, but the LCDM model provides a better fit since it has one less parameter.

In this paper causal geodesic completeness of FLRW cosmological models is analysed in terms of generalised power expansions of the scale factor in coordinate time. The strength of the found singularities is discussed following the usual definitions due to Tipler and Krolak. It is shown that while classical cosmological models are both timelike and lightlike geodesically incomplete, certain observationally alllowed models which have been proposed recently are lightlike geodesically complete.

A class of braneworld models can lead to phantom-like acceleration of the late universe, but without the need for any phantom matter. In the simplest models, the universe contains only cold dark matter and a cosmological constant. We generalize these models by introducing a quintessence field. The new feature in our models is that quintessence leads to a crossing of the phantom divide, $w=-1$. This is a purely gravitational effect, and there is no phantom instability. Furthermore, the Hubble parameter is always decreasing, and there is no big rip singularity in the future.

We consider multiple k-essence sources and obtain the conditions their kinetic functions must satisfy so that purely kinetic k-essences lead to models with phantom barrier crossing. After that, we show that polynomial kinetic functions allow the integration of the dynamical equations determining the geometry and the k-fields. The models thus obtained have accelerated expansion at late cosmic time, representing universes ending at a finite time with a big rip singularity. In addition, these models begin to evolve from an initial singularity, so they have a finite time span.

In this paper we investigate the evolution of a class of cosmologies fuelled by quintom dark energy and dark matter. Quintom dark energy is a hybrid of quintessence and phantom which involves the participation of two reals scalar fields playing the roles of those two types of dark energy. In that framework we examine from a dynamical systems perspective the possibility that those fields are coupled among them by considering an exponential potential with an interesting functional dependence similar but not identical to others studied before. The model we consider represents a counterexample to the typical behavior of quintom models with exponential potentials because it admits either tracking attractors ($w=0$), or phantom attractors ($w<-1$).

Quintessence and k-essence have been proposed as candidates for the dark energy component of the universe that would be responsible of the currently observed accelerated expansion. In this paper we investigate the degree of resemblance between those two theoretical setups, and find that every quintessence model can be viewed as a k-essence model generated by a kinetic linear function. In addition, we show the true effects of k-essence begin at second order in the expansion of the kinetic function in powers of the kinetic energy.

We use the SnIa Gold dataset to compare LCDM with 10 representative parametrizations of the recent Hubble expansion history $H(z)$. For the comparison we use two statistical tests; the usual $\chi_{min}^2$ which is insensitive to the parametrization number of parameters, and a statistic we call the p-test which depends on both the value of $\chi_{min}^2$ and the number $n$ of the parametrization parameters. The p-test measures the confidence level to which the parameter values corresponding to LCDM are excluded from the viewpoint of the parametrization tested. For example, for a linear equation of state parametrization $w(z)=w_0 + w_1 z$ the LCDM parameter values ($w_0=-1$, $w_1=0$) are excluded at 75% confidence level. We use a flat prior and $\Omega_{0m}=0.3$. All parametrizations tested are consistent with the Gold dataset at their best fit. According to both statistical tests, the worst fits among the 10 parametrizations, correspond to the Chaplygin gas, the brane world and the Cardassian parametrizations. The best fit is achieved by oscillating parametrizations which can exclude the parameter values corresponding to LCDM at 85% confidence level. Even though this level of significance does not provide a statistically significant exclusion of LCDM (it is less than $2\sigma$) and does not by itself constitute conclusive evidence for oscillations in the cosmological expansion, when combined with similar independent recent evidence for oscillations coming from the CMB and matter power spectra it becomes an issue worth of further investigation.

We extend the homogeneous modified Chaplygin cosmologies to large-scale perturbations by formulating a Zeldovich-like approximation. We show that the model interpolates between an epoch with a soft equation of state and a de Sitter phase, and that in the intermediate regime its matter content is simply the sum of dust and a cosmological constant. We then study how the large-scale inhomogeneities evolve and compare the results with cold dark matter (CDM), $\Lambda$CDM and generalized Chaplygin scenarios. Interestingly, we find that unlike the latter, our models do always resemble $\Lambda$CDM.

We consider modifications to the Friedmann equation motivated by recent proposals along these lines pursuing an explanation to the observed late time acceleration. Here we show those modifications can be framed within a theory with self-interacting gravity, where the term self-interaction refers here to the presence of functions of $\rho$ and $p$ in the right hand side of the Einstein equations. We then discuss the construction of the duals of the cosmologies generated within that framework. After that we investigate the modifications required to generate generalized and modified Chaplygin cosmologies and show that their duals belong to a larger family of cosmologies we call extended Chaplygin cosmologies. Finally, by letting the parameters of those models take values not earlier considered in the literature we show some representatives of that family of cosmologies display sudden future singularities, which indicates their behavior is rather different from generalized or modified Chaplygin gas cosmologies. This reinforces the idea that modifications of gravity can be responsible for unexpected evolutionary features in the universe.

In this paper we analyze the asymptotic behavior of Cardassian cosmological models filled with a perfect fluid and a scalar field with an exponential potential. Cardassian cosmologies arise from modifications of the Friedmann equation, and among the different proposals within that framework we will choose those of the form $3H^2-\rho\propto \rho^n$ with $n<1$. We construct two three dimensional dynamical systems arising from the evolution equations, respectively adapted for studying the high and low energy limits. Using standard dynamical systems techniques we find the fixed points and characterize the solutions they represent. We pay especial attention to the properties inherent to the modifications and compare with the (standard) unmodified scenario. Among other interesting results, we find there are no late-time tracking attractors.

We analyze the implications of having a divergent speed of sound in k-essence cosmological models. We first study a known theory of that kind, for which the Lagrangian density depends linearly on the time derivative of the k-field. We show that when k-essence is the only source consistency requires that the potential of the k-field be of the inverse square form. Then, we review the known result that the corresponding power-law solutions can be mapped to power-law solutions of theories with no divergence in the speed of sound. After that, we argue that the requirement of a divergent sound speed at some point fixes uniquely the form of the Lagrangian to be exactly the one considered earlier and prove the asymptotic stability of the most interesting solutions belonging to the divergent theory. Then, we discuss the implications having not just k-essence but also matter. This is interesting because introducing another component breaks the rigidity of the theory, and the form of the potential ceases to be unique as happened in the pure k-essence case. Finally, we show the finiteness of the effective sound speed under an appropiate definition.

We implement the transition from dark matter to dark energy in k-essence cosmologies for a very large set of kinetic functions $F$, in a way alternative to recent proposals which use generalized Chaplygin gas and transient models. Here we require that the pressure admits a power-law expansion around some value of the kinetic energy where the pressure vanishes. In addition, for suitable values of the parameters of the model, the speed of sound of the dark matter will be low. We first present the discussion in fairly general terms, and later consider for illustration two examples.

We present a general algorithm based on the concept of form-invariance which can be used for generating phantom cosmologies. It involves linear transformations between the kinetic energy and the potential of the scalar field, and transforms solutions of the Einstein-Klein-Gordon equations which preserve the weak energy condition into others which violate it, while keeping the energy density of the field positive. All known solutions representing phantom cosmologies are unified by this procedure. Using the general algorithm we obtain those solutions and show the relations between them. In addition, the scale factors of the product and seed solutions are related by a generalization of the well-known $a\to a^{-1}$ duality.

We devise a method to obtain a phantom version of FRW k-essence cosmologies with homogeneous k-fields by applying form-invariance transformations. It can be seen that the transformation performs the maps $H\toH$ and $\rho+p\to-(\rho+p)$, which in turn give $\gamma\to\gamma$ and $a\to a^{-1}$. The discussion is presented in a general setup, valid for FRW k-essence cosmologies, and then we discuss power-law solutions for illustration purposes. First, we deal with models such that the gradient of the k-field is not constant, which include standard and generalized tachyon cosmologies. We concentrate on the usual tachyon and show the phantom symmetry involves a change in the potential and that it generates an extended super-accelerated tachyon field. Then, we turn our attention to models for which the time derivative of the k-field is not constant, and we show the transformation can be implemented without changing the potential at all.

We perform a thorough phase-plane analysis of the flow defined by the equations of motion of a FRW universe filled with a tachyonic fluid plus a barotropic one. The tachyon potential is assumed to be of inverse square form, thus allowing for a two-dimensional autonomous system of equations. The Friedmann constraint, combined with a convenient choice of coordinates, renders the physical state compact. We find the fixed-point solutions, and discuss whether they represent attractors or not. The way the two fluids contribute at late-times to the fractional energy density depends on how fast the barotropic fluid redshifts. If it does it fast enough, the tachyonic fluid takes over at late times, but if the opposite happens, the situation will not be completely dominated by the barotropic fluid; instead there will be a residual non-negligible contribution from the tachyon subject to restrictions coming from nucleosynthesis.

We present here the general transformation that leaves unchanged the form of the field equations for perfect fluid Friedmann--Robertson--Walker and Bianchi V cosmologies. The symmetries found can be used as algorithms for generating new cosmological models from existing ones. A particular case of the general transformation is used to illustrate the crucial role played by the number of scalar fields in the occurrence of inflation. Related to this, we also study the existence and stability of Bianchi V power law solutions.

We construct exact solutions to Einstein equations which represent relativistic disks immersed into an expanding FRW Universe. It is shown that the expansion influences dynamical characteristics of the disks such as rotational curves, surface mass density, etc. The effects of the expansion is exemplified with non-static generalizations of Kuzmin-Curzon and generalized Schwarzschild disks.