Sub-horizon perturbations in scalar-tensor theories have been shown to grow generically faster than in uncoupled models due to a positive, additive Yukawa force. In such cases, the amount of clustering becomes larger than in the standard cosmological model, exacerbating the observed tension in the $\sigma_8$ parameter. Here we show instead that in some simple cases of conformally coupled dark energy one can obtain a transient regime of negative Yukawa force, without introducing ghosts or other instabilities.

We address the issue of the existence of inequivalent definitions of gravitational mass in $R^{2}$-gravity. We present several definitions of gravitational mass, and discuss the formal relations between them. We then consider the concrete case of a static and spherically symmetric neutron star, and solve numerically the equations of motion for several values of the free parameter of the model. We compare the features of the mass-radius relations obtained for each definition of gravitational mass, and we comment on their dependence on the free parameter. We then argue that $R^{2}$-gravity is a valuable proxy to discuss the existence of inequivalent definitions of gravitational mass in a generic modified gravity theory, and present some comments on the general case.

We review the features of Dark Matter as a particle, presenting some old and new instructive models, and looking for their physical implications in the early universe and in the process of structure formation. We also present a schematic of Dark Matter searches and introduce the most promising candidates to the role of Dark Matter particle.

We investigate the behavior of the asymptotic late-times effective equation of state for a class of nonlocal theories of gravity. These theories modify the Einstein-Hilbert Lagrangian introducing terms containing negative powers of the d'Alembert operator acting on the Ricci scalar. We find that imposing vanishing initial conditions for the nonlocal content during the radiation-dominated epoch implies the same asymptotic late-times behavior for most of these models. In terms of the effective equation of state of the universe, we find that asymptotically $\omega_{ \rm eff} \rightarrow -1$, approaching the value given by a cosmological constant. On the other hand, unlike in the case of $\Lambda$CDM, the Hubble factor is a monotonic growing function that diverges asymptotically. We argue that this behavior is not a coincidence and discuss under which conditions this is to be expected in these nonlocal models.

In this paper we establish \textit{formulae} for the inflationary slow-roll parameters $\epsilon$, $\eta$ and $\zeta$ as functions of the Ricci scalar $R$ for $f(R)$ theories of gravity. As examples, we present the analytic and numerical solutions of $\epsilon$, $\eta$ and $\zeta$ as functions of the number of e-folds $N$ in two important instances: for the Starobinsky model and for a $f(R)$ reconstruction of the $\alpha$-Attractors. The highlight of our proposal is to rewrite the slow-roll parameters in terms of $f(R)$, which allows to find directly $n_{\rm s}$, $r$, $\alpha_{\rm s}$ and $f^{\text{equil}}_{\text{NL}}$ as functions of $R$ itself. We obtain that both models indicate a small contribution to the non-Gaussianity parameters, which are in good agreement with current observational constraints.

We investigate the linear regime of f(R) = R + alpha R^2 gravity for static, spherically symmetric and asymptotically flat configurations of matter. We show that, in vacuum and deep inside the range of the extra scalar degree of freedom, the post-Newtonian parameter gamma is not equal to 1/2, as established in the literature, but it assumes larger values depending on the pressure of the star. We provide an explicit expression for gamma in terms of the mass, of the integrated pressure of the star and of the ratio between the star's radius and the range of the extra degree of freedom. We corroborate our results by providing numerical solutions for the case of a neutron star.

These lecture notes are based on the hand-written notes which I prepared for the cosmology course taught to graduate students of PPGFis and PPGCosmo at the Federal University of Esp\'irito Santo (UFES), starting from 2014. This course covers topics ranging from the evidence of the expanding universe to Cosmic Microwave Background anisotropies. They can be found also on my personal webpage http://ofp.cosmo-ufes.org/ and shall be published by Springer in 2018.

We show an analogy at high curvature between a $f(R) = R + aR^{n - 1} + bR^2$ theory and the $\alpha$-Attractors. We calculate the expressions of the parameters $a$, $b$ and $n$ as functions of $\alpha$ and the predictions of the model $f(R) = R + aR^{n - 1} + bR^2$ on the scalar spectral index $n_{\rm s}$ and the tensor-to-scalar ratio $r$. We find that the power law correction $R^{n - 1}$ allows for a production of gravitational waves enhanced with respect to the one in the Starobinsky model, while maintaining a viable prediction on $n_{\rm s}$. We numerically reconstruct the full $\alpha$-Attractors class of models testing the goodness of our high-energy approximation $f(R) = R + aR^{n - 1} + bR^2$. Moreover, we also investigate the case of a single power law $f(R) = \gamma R^{2 - \delta}$ theory, with $\gamma$ and $\delta$ free parameters. We calculate analytically the predictions of this model on the scalar spectral index $n_{\rm s}$ and the tensor-to-scalar ratio $r$ and the values of $\delta$ which are allowed from the current observational results. We find that $-0.015 < \delta < 0.016$, confirming once again the excellent agreement between the Starobinsky model and observation.

We investigate the effect of the redshift drift in strong gravitational lensing. The redshift drift produces a time variation of $i)$ the apparent position of a lensed source and $ii)$ the time delay among incoming signals from different images. We dub these effects as angular drift and time delay drift, respectively, and analyze their relevance in cosmology.

The k-essence theory with a power-law function of $(\partial\phi)^2$ and Rastall's non-conservative theory of gravity with a scalar field are shown to have the same solutions for the metric under the assumption that both the metric and the scalar fields depend on a single coordinate. This equivalence (called k-R duality) holds for static configurations with various symmetries (spherical, plane, cylindrical, etc.) and all homogeneous cosmologies. In the presence of matter, Rastall's theory requires additional assumptions on how the stress-energy tensor non-conservation is distributed between different contributions. Two versions of such non-conservation are considered in the case of isotropic spatially flat cosmological models with a perfect fluid: one (R1) in which there is no coupling between the scalar field and the fluid, and another (R2) in which the fluid separately obeys the usual conservation law. In version R1 it is shown that k-R duality holds not only for the cosmological models themselves but also for their adiabatic perturbations. In version R2, among other results, a particular model is singled out that reproduces the same cosmological expansion history as the standard $\Lambda$CDM model but predicts different behaviors of small fluctuations in the k-essence and Rastall frameworks.

We consider cosmological backreaction effects in Buchert's averaging formalism on the basis of an explicit solution of the Lema\^itre-Tolman-Bondi (LTB) dynamics which is linear in the LTB curvature parameter and has an inhomogeneous bang time. The volume Hubble rate is found in terms of the volume scale factor which represents a derivation of the simplest phenomenological solution of Buchert's equations in which the fractional densities corresponding to average curvature and kinematic backreaction are explicitly determined by the parameters of the underlying LTB solution at the boundary of the averaging volume. This configuration represents an exactly solvable toy model but it does not adequately describe our "real" Universe.

We analyse $f(R)$ theories of gravity from a dynamical system perspective, showing how the $R^2$ correction in Starobinsky's model plays a crucial role from the viewpoint of the inflationary paradigm. Then, we propose a modification of Starobinsky's model by adding an exponential term in the $f(R)$ Lagrangian. We show how this modification could allow to test the robustness of the model by means of the predictions on the scalar spectral index $n_s$.

There are some approaches, either based on General Relativity (GR) or modified gravity, that use galaxy rotation curves to derive the matter density of the corresponding galaxy, and this procedure would either indicate a partial or a complete elimination of dark matter in galaxies. Here we review these approaches, clarify the difficulties on this inverted procedure, present a method for evaluating them, and use it to test two specific approaches that are based on GR: the Cooperstock-Tieu (CT) and the Balasin-Grumiller (BG) approaches. Using this new method, we find that neither of the tested approaches can satisfactorily fit the observational data without dark matter. The CT approach results can be significantly improved if some dark matter is considered, while for the BG approach no usual dark matter halo can improve its results.

We investigate the effect of the cosmological expansion on the bending of light due to an isolated point-like mass. We adopt McVittie metric as the model for the geometry of the lens. Assuming a constant Hubble factor we find an analytic expression involving the bending angle, which turns out to be unaffected by the cosmological expansion at the leading order.

We revisit the analysis made by Hwang and Noh [JCAP 1310 (2013)] aiming the construction of a Newtonian set of equations incorporating pressure effects typical of the General Relativity theory. We explicitly derive the Hwang-Noh equations, comparing them with similar computations found in the literature. Then, we investigate $i)$ the cosmological expansion, $ii)$ linear cosmological perturbations theory and $iii)$ stellar equilibrium by using the new set of equations and comparing the results with those coming from the usual Newtonian theory, from the Neo-Newtonian theory and from the General Relativity theory. We show that the predictions for the background evolution of the Universe are deeply changed with respect to the General Relativity theory: the acceleration of the Universe is achieved with positive pressure. On the other hand, the behaviour of small cosmological perturbations reproduces the one found in the relativistic context, even if only at small scales. We argue that this last result may open new possibilities for numerical simulations for structure formation in the Universe. Finally, the properties of neutron stars are qualitatively reproduced by Hwang-Noh equations, but the upper mass limit is at least one order of magnitude higher than the one obtained in General Relativity.

We present a new approximation to include fully general relativistic pressure and velocity in Newtonian hydrodynamics. The energy conservation, momentum conservation and two Poisson's equations are consistently derived from Einstein's gravity in the zero-shear gauge assuming weak gravity and action-at-a-distance limit. The equations show proper special relativity limit in the absence of gravity. Our approximation is complementary to the post-Newtonian approximation and the equations are valid in fully nonlinear situations.

We consider an effective viscous pressure as the result of a backreaction of inhomogeneities within Buchert's formalism. The use of an effective metric with a time-dependent curvature radius allows us to calculate the luminosity distance of the backreaction model. This quantity is different from its counterpart for a "conventional" spatially flat bulk viscous fluid universe. Both expressions are tested against the SNIa data of the Union2.1 sample with only marginally different results for the distance-redshift relation and in accordance with the $\Lambda$CDM model. Future observations are expected to be able to discriminate among these models on the basis of indirect measurements of the curvature evolution.

The Brans-Dicke action is one of the most natural extensions of the Einstein-Hilbert action. It is based on the introduction of a fundamental scalar field that effectively incorporates a dynamics to the gravitational coupling $G$. In spite of the diverse motivations and the rich phenomenology that comes from its solutions, Solar System tests impose strong constraints on the Brans-Dicke theory, rendering it indistinguishable from General Relativity. In the present text, new perspectives for the Brans-Dicke theory are presented, based on the possibility that the scalar field presented in the BD theory can be external, as well as on the applications to black hole physics and the primordial universe.

Effects of velocity dispersion of dark matter particles on the CMB TT power spectrum and on the matter linear power spectrum are investigated using a modified CAMB code. Cold dark matter originated from thermal equilibrium processes does not produce appreciable effects but this is not the case if particles have a non-thermal origin. A cut-off in the matter power spectrum at small scales, similar to that produced by warm dark matter or that produced in the late forming dark matter scenario, appears as a consequence of velocity dispersion effects, which act as a pressure perturbation.

We discuss the possibility to implement a viscous cosmological model, attributing to the dark matter component a behaviour described by bulk viscosity. Since bulk viscosity implies negative pressure, this rises the possibility to unify the dark sector. At the same time, the presence of dissipative effects may alleviate the so called small scale problems in the $\Lambda$CDM model. While the unified viscous description for the dark sector does not lead to consistent results, the non-linear behaviour indeed improves the situation with respect to the standard cosmological model.

We study the effects of sample variance in N--body simulations, as a function of the size of the simulation box, namely in connection with predictions on tomographic shear spectra. We make use of a set of 8 $\Lambda$CDM simulations in boxes of 128, 256, 512 $h^{-1}$Mpc aside, for a total of 24, differing just by the initial seeds. Among the simulations with 128 and 512 $h^{-1}$Mpc aside, we suitably select those closest and farthest from {\it average}. Numerical and linear spectra $P(k,z)$ are suitably connected at low $k$ so to evaluate the effects of sample variance on shear spectra $C_{ij}(\ell)$ for 5 or 10 tomographic bands. We find that shear spectra obtained by using 128 $h^{-1}$Mpc simulations can vary up to $\sim 25\, \%$, just because of the seed. Sample variance lowers to $\sim 3.3\, \%$, when using 512 $h^{-1}$Mpc. These very percentages could however slightly vary, if other sets of the same number of realizations were considered. Accordingly, in order to match the $\sim 1\, \%$ precision expected for data, if still using 8 boxes, we require a size $\sim 1300$ --$ 1700 \, h^{-1}$ Mpc for them.

In this paper, we explore two different ways of implementing quantum effects in a classical structure. The first one is through an external field. The other one is modifying the classical conservation laws. In both cases, the consequences for the description of the evolution of the universe are discussed.

The cosmic microwave background (CMB) radiation provides a remarkable window onto the early universe, revealing its composition and structure. In these lectures we review and discuss the physics underlying the main features of the CMB.

We formulate a theory combining the principles of a scalar-tensor gravity and the Rastall proposal of a violation of the usual conservation laws. In the resulting Brans-Dicke-Rastall (BDR) theory the only exact, static, spherically symmetric solution is a Robinson-Bertotti type solution besides the trivial Schwarzschild one. The PPN constraints can be completely satisfied for some values of the free parameters.The cosmological solutions display, among others, a decelerate-accelerate transition in the matter dominated phase.

We formulate a theory combining the principles of a scalar-tensor gravity and Rastall's proposal of a violation of the usual conservation laws. We obtain a scalar-tensor theory with two parameters $\omega$ and $\lambda$, the latter quantifying the violation of the usual conservation laws. The only exact spherically symmetric solution is that of Robinson-Bertotti besides Schwarzschild solution. A PPN analysis reveals that General Relativity results are reproduced when $\lambda = 0$. The cosmological case displays a possibility of deceleration/acceleration or acceleration/deceleration transitions during the matter dominated phase depending on the values of the free parameters.

Rastall's theory is a modification of Einstein's theory of gravity where the covariant divergence of the stress-energy tensor is no more vanishing, but proportional to the gradient of the Ricci scalar. The motivation of this theory is to investigate a possible non-minimal coupling of the matter fields to geometry which, being proportional to the curvature scalar, may represent an effective description of quantum gravity effects. Non-conservation of the stress-energy tensor, via Bianchi identities, implies new field equations which have been recently used in a cosmological context, leading to some interesting results. In this paper we adopt Rastall's theory to reproduce some features of the effective Friedmann's equation emerging from loop quantum cosmology. We determine a class of bouncing cosmological solutions and comment about the possibility of employing these models as effective descriptions of the full quantum theory.

We consider a cosmological model of the late universe constituted by standard cold dark matter plus a dark energy component with constant equation of state w and constant effective speed of sound. Neglecting fluctuations in the dark energy component we obtain an equation describing the evolution of sub-horizon cold dark matter perturbations through the epoch of dark matter-dark energy equality. We explore its analytic solutions and calculate an exact w-dependent correction for the dark matter growth function, logarithmic growth function and growth index parameter through the epoch considered. We test our analytic approximation with the numerical solution and find that the discrepancy is less than 1% for k =0 in the epoch of interest.

The matter creation cosmology is revisited, including the evolution of baryons and dark matter particles. The creation process affects only dark matter and not baryons. The dynamics of the $\Lambda$CDM model can be reproduced only if two conditions are satisfied: 1) the entropy density production rate and the particle density variation rate are equal and 2) the (negative) pressure associated to the creation process is constant. However, the matter creation model predicts a present dark matter-to-baryon ratio much larger than that observed in massive X-ray clusters of galaxies, representing a potential difficulty for the model. In the linear regime, a fully relativistic treatment indicates that baryons are not affected by the creation process but this is not the case for dark matter. Both components evolve together at early phases but lately the dark matter density contrast decreases since the background tends to a constant value. This behaviour produces a negative growth factor, in disagreement with observations, being a further problem for this cosmology.

The Rastall's theory is a modification of General Relativity touching one of the cornestone of gravity theory: the conservation laws. In Rastall's theory, the energy-momentum tensor is not conserved anymore, depending now on the gradient of the Ricci curvature. In this sense, this theory can be seen as a classical implementation of quantum effects in a curved background space-time. We exploit this structure in order to reproduce some results of an effective theory of quantum loop cosmology. Later, we propose a model for the dark sector of the universe. In this case, the corresponding $\Lambda$CDM model appears as the only model consistent with observational data.

We investigate the correspondence between a perfect fluid and a scalar field and show a possible way of expressing thermodynamic quantities such as entropy, particle number density, temperature and chemical potential in terms of the scalar field phi and its kinetic term X. We prove a theorem which relates isentropy with purely kinetic Lagrangian. As an application, we study the evolution of the gravitational potential in cosmological perturbation theory.

Newtonian fluid dynamics allows the construction of acoustic metrics from which black hole configurations can be studied. However, relativistic pressure effects are neglected within Newtonian theory. We study acoustic black holes in the framework of neo-Newtonian hydrodynamics, which is designed to take into account relativistic inertial effects of the pressure $p$. Within this new hydrodynamical context we show how $p$ can influence the formation of the acoustic horizons.

We discuss solutions of Vlasov-Einstein equation for collisionless dark matter particles in the context of a flat Friedmann universe. We show that, after decoupling from the primordial plasma, the dark matter phase-space density indicator Q remains constant during the expansion of the universe, prior to structure formation. This well known result is valid for non-relativistic particles and is not "observer dependent" as in solutions derived from the Vlasov-Poisson system. In the linear regime, the inclusion of velocity dispersion effects permits to define a physical Jeans length for collisionless matter as function of the primordial phase-space density indicator: \lambda_J = (5\pi/G)^(1/2)Q^(-1/3)\rho_dm^(-1/6). The comoving Jeans wavenumber at matter-radiation equality is smaller by a factor of 2-3 than the comoving wavenumber due to free-streaming, contributing to the cut-off of the density fluctuation power spectrum at the lowest scales. We discuss the physical differences between these two scales. For dark matter particles of mass equal to 200 GeV, the derived Jeans mass is 4.3 x 10^(-6) solar masses.

We investigate $f(R,T)$ gravity models ($R$ is the curvature scalar and $T$ is the trace of the stress-energy tensor of ordinary matter) that are able to reproduce the four known types of future finite-time singularities. We choose a suitable expression for the Hubble parameter in order to realise the cosmic acceleration and we introduce two parameters, $\alpha$ and $H_s$, which characterise each type of singularity. We address conformal anomaly and we observe that it cannot remove the sudden singularity or the type IV one, but, for some values of $\alpha$, the big rip and the type III singularity may be avoided. We also find that, even without taking into account conformal anomaly, the big rip and the type III singularity may be removed thanks to the presence of the $T$ contribution of the $f(R,T)$ theory.

We investigate the unification scenario provided by the generalised Chaplygin gas model (a perfect fluid characterized by an equation of state p = -A/\rho^{\alpha}). Our concerns lie with a possible tension existing between background kinematic tests and those related to the evolution of small perturbations. We analyse data from the observation of the differential age of the universe, type Ia supernovae, baryon acoustic oscillations and the position of the first peak of the angular spectrum of the cosmic background radiation. We show that these tests favour negative values of the parameter \alpha: we find \alpha = - 0.089^{+0.161}_{-0.128} at the 2\sigma\ level and that \alpha < 0 with 85% confidence. These would correspond to negative values of the square speed of sound which are unacceptable from the point of view of structure formation. We discuss a possible solution to this problem, when the generalised Chaplygin gas is framed in the modified theory of gravity proposed by Rastall. We show that a fluid description within this theory does not serve the purpose, but it is necessary to frame the generalised Chaplygin gas in a scalar field theory. Finally, we address the standard general relativistic unification picture provided by the generalised Chaplygin gas in the case \alpha = 0: this is usually considered to be undistinguishable from the standard \Lambda CDM model, but we show that the evolution of small perturbations, governed by the M\'esz\'aros equation, is indeed different and the formation of sub-horizon GCG matter halos may be importantly affected in comparison with the \Lambda CDM scenario.

Rastall's theory is a modification of General Relativity, based on the non-conservation of the stress-energy tensor. The latter is encoded in a parameter $\gamma$ such that $\gamma = 1$ restores the usual $\nabla_\nu T^{\mu\nu} = 0$ law. We test Rastall's theory in cosmology, on a flat Robertson-Walker metric, investigating a two-fluid model and using the type Ia supernovae Constitution dataset. One of the fluids is pressureless and obeys the usual conservation law, whereas the other is described by an equation of state $p_x = w_x\rho_x$, with $w_x$ constant. The Bayesian analysis of the Constitution set does not strictly constrain the parameter $\gamma$ and prefers values of $w_x$ close to -1. We then address the evolution of small perturbations and show that they are dramatically unstable if $w_x \neq -1$ and $\gamma \neq 1$, i.e. General Relativity is the favored configuration. The only alternative is $w_x = -1$, for which the dynamics becomes independent from $\gamma$.

We explore the phenomenology of nontrivial quantum effects on low-energy gravity. These effects come from the running of the gravitational coupling parameter G and the cosmological constant L in the Einstein-Hilbert action, as induced by the Renormalization Group (RG). The Renormalization Group corrected General Relativity (RGGR model) is used to parametrize these quantum effects, and it is assumed that the dominant dark matter-like effects inside galaxies is due to these nontrivial RG effects. Here we present additional details on the RGGR model application, in particular on the Poisson equation extension that defines the effective potential, also we re-analyse the ordinary elliptical galaxy NGC 4494 using a slightly different model for its baryonic contribution, and explicit solutions are presented for the running of G and L. The values of the NGC 4494 parameters as shown here have a better agreement with the general RGGR picture for galaxies, and suggest a larger radial anisotropy than the previously published result.

We review the difficulties of the generalized Chaplygin gas model to fit observational data, due to the tension between background and perturbative tests. We argue that such issues may be circumvented by means of a self-interacting scalar field representation of the model. However, this proposal seems to be successful only if the self-interacting scalar field has a non-canonical form. The latter can be implemented in Rastall's theory of gravity, which is based on a modification of the usual matter conservation law. We show that, besides its application to the generalized Chaplygin gas model, other cosmological models based on Rastall's theory have many interesting and unexpected new features.

We investigate the evolution of the gravitational potential in Rastall scalar field theories. In a single component model a consistent perturbation theory, formulated in the newtonian gauge, is possible only for $\gamma = 1$, which is the General Relativity limit. On the other hand, the addition of another canonical fluid component allows also to consider the case $\gamma \neq 1$.

We review the difficulties of the generalized Chaplygin gas model to fit observational data, due to the tension between background and perturbative tests. We argue that such issues may be circumvented by means of a self-interacting scalar field representation of the model. However, this proposal seems to be successful only if the self-interacting scalar field has a non-canonical form. The latter can be implemented in Rastall's theory of gravity.

Rastall's theory is based on the non-conservation of the energy-momentum tensor. We show that, in this theory, if we introduce a two-fluid model, one component representing vacuum energy whereas the other pressureless matter (e.g. baryons plus cold dark matter), the cosmological scenario is the same as for the \Lambda CDM model, both at background and linear perturbative levels, except for one aspect: now dark energy may cluster. We speculate that this can lead to a possibility of distinguishing the models at the non-linear perturbative level.

We consider cosmological scenarios based on $f(R,T)$ theories of gravity ($R$ is the Ricci scalar and $T$ is the trace of the energy-momentum tensor) and numerically reconstruct the function $f(R,T)$ which is able to reproduce the same expansion history generated, in the standard General Relativity theory, by dark matter and holographic dark energy. We consider two special $f(R,T)$ models: in the first instance, we investigate the modification $R + 2f(T)$, i.e. the usual Einstein-Hilbert term plus a $f(T)$ correction. In the second instance, we consider a $f(R)+\lambda T$ theory, i.e. a $T$ correction to the renown $f(R)$ theory of gravity.

The generalized Chaplygin gas model represents an attempt to unify dark matter and dark energy. It is characterized by a fluid with an equation of state $p = - A/\rho^\alpha$. It can be obtained from a generalization of the DBI action for a scalar, tachyonic field. At background level, this model gives very good results, but it suffers from many drawbacks at perturbative level. We show that, while for background analysis it is possible to consider any value for $\alpha$, the perturbative analysis must be restricted to positive values of $\alpha$. This restriction can be circumvented if the origin of the generalized Chaplygin gas is traced back to a self-interacting scalar field, instead of the DBI action. But, in doing so, the predictions coming from formation of large scale structures reduce the generalized Chaplygin gas model to a kind of quintessence model, and the unification scenario is lost, if the scalar field is the canonical one. However, if the unification condition is imposed from the beginning as a prior, the model may remain competitive. More interesting results, concerning the unification program, are obtained if a non-canonical self-interacting scalar field, inspired by Rastall's theory of gravity, is imposed. In this case, an agreement with the background tests is possible.

We present an analysis of the cross-correlation between the CMB and the large-scale structure (LSS) of the Universe in Unified Dark Matter (UDM) scalar field cosmologies. We work out the predicted cross-correlation function in UDM models, which depends on the speed of sound of the unified component, and compare it with observations from six galaxy catalogues (NVSS, HEAO, 2MASS, and SDSS main galaxies, luminous red galaxies, and quasars). We sample the value of the speed of sound and perform a likelihood analysis, finding that the UDM model is as likely as the LambdaCDM, and is compatible with observations for a range of values of c_\infinity (the value of the sound speed at late times) on which structure formation depends. In particular, we obtain an upper bound of c_\infinity^2 \leq 0.009 at 95% confidence level, meaning that the LambdaCDM model, for which c_\infinity^2 = 0, is a good fit to the data, while the posterior probability distribution peaks at the value c_\infinity^2=10^(-4) . Finally, we study the time dependence of the deviation from LambdaCDM via a tomographic analysis using a mock redshift distribution and we find that the largest deviation is for low-redshift sources, suggesting that future low-z surveys will be best suited to constrain UDM models.

We consider cosmological scenarios originating from a single imperfect fluid with bulk viscosity and apply Eckart's and both the full and the truncated M\"uller-Israel-Stewart's theories as descriptions of the non-equilibrium processes. Our principal objective is to investigate if the dynamical properties of Dark Matter and Dark Energy can be described by a single viscous fluid and how such description changes when a causal theory (M\"uller-Israel-Stewart's, both in its full and truncated forms) is taken into account instead of Eckart's non-causal theory. To this purpose, we find numerical solutions for the gravitational potential and compare its behaviour with the corresponding LambdaCDM case. Eckart's and the full causal theory seem to be disfavoured, whereas the truncated theory leads to results similar to those of the LambdaCDM model for a bulk viscous speed in the interval 10^{-11} << c_b^2 < 10^{-8}. Tentatively relating such value to a square propagation velocity of the order of T/m of perturbations in a non-relativistic gas of particles with mass m at the epoch of matter-radiation equality, this may be compatible with a mass range 0.1 GeV < m << 100 GeV.

We investigate the time evolution of the gravitational potential Phi for a special class of non-adiabatic Unified Dark matter Models described by scalar field lagrangians. These models predict the same background evolution as in the LambdaCDM and possess a non-vanishing speed of sound. We provide a very accurate approximation of Phi, valid after the recombination epoch, in the form of a Bessel function of the first kind. This approximation may be useful for a future deeper analysis of Unified Dark Matter scalar field models.

We investigate the general properties of Unified Dark Matter (UDM) scalar field models with Lagrangians with a non-canonical kinetic term, looking specifically for models that can produce a fast transition between an early Einstein-de Sitter CDM-like era and a later Dark Energy like phase, similarly to the barotropic fluid UDM models in JCAP1001(2010)014. However, while the background evolution can be very similar in the two cases, the perturbations are naturally adiabatic in fluid models, while in the scalar field case they are necessarily non-adiabatic. The new approach to building UDM Lagrangians proposed here allows to escape the common problem of the fine-tuning of the parameters which plague many UDM models. We analyse the properties of perturbations in our model, focusing on the the evolution of the effective speed of sound and that of the Jeans length. With this insight, we can set theoretical constraints on the parameters of the model, predicting sufficient conditions for the model to be viable. An interesting feature of our models is that what can be interpreted as w_{DE} can be <-1 without violating the null energy conditions.

Unified Dark Matter models describe Dark Matter and Dark Energy as a single entity which is, in the simplest case, embodied in a perfect barotropic fluid. It is a well-established fact that small adiabatic perturbations of Unified Dark Matter have an evolution characterised by oscillations and decay which provide predictions on the Cosmic Background Radiation anisotropies which are in poor agreement with observation. In this paper we investigate the generalised Chaplygin gas and we find that the Integrated Sachs-Wolfe effect excludes the model for $10^{-3} < \alpha < 350$. We discuss the implications on the background evolution of the Universe if large values of $\alpha$ are considered. In this case, the Universe expansion mimics a matter-dominated phase abruptly followed by a de Sitter one at the transition redshift $z_tr$. Thanks to an analysis of the type Ia supernovae Constitution set we are able to place $z_tr = 0.22$.

We investigate the general properties of Unified Dark Matter (UDM) fluid models where the pressure and the energy density are linked by a barotropic equation of state (EoS) $p = p(\rho)$ and the perturbations are adiabatic. The EoS is assumed to admit a future attractor that acts as an effective cosmological constant, while asymptotically in the past the pressure is negligible. UDM models of the dark sector are appealing because they evade the so-called "coincidence problem" and "predict" what can be interpreted as $w_{\rm DE} \approx -1$, but in general suffer the effects of a non-negligible Jeans scale that wreak havoc in the evolution of perturbations, causing a large Integrated Sachs-Wolfe effect and/or changing structure formation at small scales. Typically, observational constraints are violated, unless the parameters of the UDM model are tuned to make it indistinguishable from $\Lambda$CDM. Here we show how this problem can be avoided, studying in detail the functional form of the Jeans scale in adiabatic UDM perturbations and introducing a class of models with a fast transition between an early Einstein-de Sitter CDM-like era and a later $\Lambda$CDM-like phase. If the transition is fast enough, these models may exhibit satisfactory structure formation and CMB fluctuations. To consider a concrete case, we introduce a toy UDM model and show that it can predict CMB and matter power spectra that are in agreement with observations for a wide range of parameter values.

We investigate the Tolman-Oppenheimer-Volkoff equations for the generalized Chaplygin gas with the aim of extending the findings of V. Gorini, U. Moschella, A. Y. Kamenshchik, V. Pasquier, and A. A. Starobinsky [Phys. Rev. D {\bf 78}, 064064 (2008)]. We study both the standard case, where we reproduce some previous results, and the phantom case. In the phantom case we show that even a superluminal group velocity arising for $\alpha > 1$ cannot prevent the divergence of the pressure at a finite radial distance. Finally, we investigate how a modification of the generalized Chaplygin gas equation of state, required by causality arguments at densities very close to $\Lambda$, affects the results found so far.

We exploit the gauge-invariant formalism to analyse the perturbative behaviour of two cosmological models based on the generalized Chaplygin gas describing both dark matter and dark energy in the present Universe. In the first model we consider the generalized Chaplygin gas alone, while in the second one we add a baryon component to it. We extend our analysis also into the parameter range $\alpha > 1$, where the generalized Chaplygin gas sound velocity can be larger than that of light. In the first model we find that the matter power spectrum is compatible with the observed one only for $\alpha < 10^{-5}$, which makes the generalized Chaplygin gas practically indistinguishable from $\Lambda$CDM. In the second model we study the evolution of inhomogeneities of the baryon component. The theoretical power spectrum is in good agreement with the observed one for almost all values of $\alpha$. However, the growth of inhomogeneities seems to be particularly favoured either for sufficiently small values of $\alpha$ or for $\alpha \gtrsim 3$. Thus, it appears that the viability of the generalized Chaplygin gas as a cosmological model is stronger when its sound velocity is superluminal. We show that in this case the generalized Chaplygin gas equation of state can be changed in an unobservable region in such a way that its equivalent $k$-essence microscopical model has no problems with causality.