We present the results of a full cosmological simulation with the new code SCALAR, where dark matter is in form of fuzzy dark matter, described by a light scalar field with a mass of $m_{\rm B} = 2.5 \times 10^{-22}$ eV and evolving according to the Schr\"{o}dinger-Poisson system of equations. In comoving units, the simulation volume is $2.5 ~ h^{-1} {\rm Mpc}$ on a side, with a resolution of $20~h^{-1}{\rm pc}$ at the finest refinement level. We analyse the formation and the evolution of central solitonic cores, which are found to leave their imprints on dark matter density profiles, resulting in shallower central densities, and on rotation curves, producing an additional circular velocity peak at small radii from the center. We find that the suppression of structures due to the quantum nature of the scalar field results in an shallower halo mass function in the low-mass end compared to the case of a $\Lambda$CDM simulation, in which dark matter is expected to cluster at all mass scales even if evolved with the same initial conditions used for fuzzy dark matter. Furthermore, we verify the scaling relations characterising the solution to the Schr\"{o}dinger-Poisson system, for both isolated and merging halos, and we find that they are preserved by merging processes. We characterise each fuzzy dark matter halo in terms of the dimensionless quantity $\Xi \propto \left | E_{\rm halo} \right |/M_{\rm halo}^3$ and we show that the core mass is tightly linked to the halo mass by the core-halo mass relation $M_{\rm core}/M_{\rm halo} \propto \Xi^{1/3}$. We also show that the core surface density of the simulated fuzzy dark matter halos does not follow the scaling with the core radius as observed for dwarf galaxies, representing a big challenge for the fuzzy dark matter model as the sole explanation of core formation.

We present a new code, SCALAR, based on the high-resolution hydrodynamics and N-body code RAMSES, to solve the Schr\"odinger equation on adaptive refined meshes. The code is intended to be used to simulate axion or fuzzy dark matter models where the evolution of the dark matter component is determined by a coupled Schr\"odinger-Poisson equation, but it can also be used as a standalone solver for both linear and non-linear Schr\"odinger equations with any given external potential. This paper describes the numerical implementation of our solver and presents tests to demonstrate how accurately it operates.

It is always interesting to investigate how well can a future experiment perform with respect to others (present or future ones). Cosmology is really an exciting field where a lot of puzzles are still unknown. In this article we consider a generalized dark energy (DE) scenario where anisotropic stress is present. We constrain this generalized cosmic scenario with an aim to investigate how gravitational waves standard sirens (GWSS) may constrain the anisotropic stress, which according to the standard cosmological probes, remains unconstrained. In order to do this, we generate the luminosity distance measurements from $\mathcal{O} (10^3)$ mock GW events which match the expected sensitivity of the Einstein Telescope. Our analyses report that, first of all, GWSS can give better constraints on various cosmological parameters compared to the usual cosmological probes, but the viscous sound speed appearing due to the dark energy anisotropic stress, is totally unconstrained even after the inclusion of GWSS.

Searching for violations of the no-hair theorem (NHT) is a powerful way to test gravity, and more generally fundamental physics, particularly with regards to the existence of additional scalar fields. The first observation of a black hole (BH) shadow by the Event Horizon Telescope (EHT) has opened a new direct window onto tests of gravity in the strong-field regime, including probes of violations of the NHT. We consider two scenarios described by the Einstein-Maxwell equations of General Relativity and electromagnetism, to which we add a scalar field. In the first case we consider a minimally-coupled scalar field with a potential, whereas in the second case the field is conformally-coupled to curvature. In both scenarios we construct charged BH solutions, which are found to carry primary scalar hair. We then compute the shadows cast by these two BHs as a function of their electric charge and scalar hair parameter. Comparing these shadows to the shadow of M87* recently imaged by the EHT collaboration, we set constraints on the amount of scalar hair carried by these two BHs. The conformally-coupled case admits a regime for the hair parameter, compatible with EHT constraints, describing a so-called mutated Reissner-Nordstr\"{o}m BH: this solution was recently found to effectively mimic a wormhole. Our work provides novel constraints on fundamental physics, and in particular on violations of the no-hair theorem and the existence of additional scalar fields, from the shadow of M87*.

Motivated by the discrepancy in measurements of $H_0$ between local and global probes, we investigate whether teleparallel gravities could be a better model to describe the present days observations or at least to alleviate the $H_0$ tension. Specifically, in this work we study and place constraints on three popular $f(T)$ models in light of the Planck-2018 CMB data release. We find that the $f(T)$ power-law model can alleviate the $H_0$ tension from $4.4\sigma$ to $1.9\sigma$ level, while the $f(T)$ model of two exponential fail to resolve this inconsistency. Moreover, for the first time, we obtain constraints on the effective number of relativistic species $N_{eff}$ and on the sum of the neutrino masses $\Sigma m_\nu$ in $f(T)$ gravity. We find that the constraints obtained are looser than in $\Lambda$CDM. However, the introduction of massive neutrinos into the cosmological model alleviate the $H_0$ tension for the power-law model. Finally, we find that whether a viable $f(T)$ theory can mitigate the $H_0$ tension depends on the mathematical structure of the distortion factor $y(z,\,b)$. These results could provide a clue for theoreticians to write a more physical-motivated expression of $f(T)$ function.

We investigate a well known scenario of interaction in the dark sector where the vacuum energy is interacting with cold dark matter throughout the cosmic evolution in light of the cosmic microwave background (CMB) data from final Planck 2018 release. In addition to this minimal scenario, we generalize the model baseline by including the properties of neutrinos, such as the neutrino mass scale ($M_{\nu}$) and the effective number of neutrino species ($N_{\rm eff}$) as free parameters, in order to verify the possible effects that such parameters might generate on the coupling parameter, and vice versa. As already known, we again confirm that in light of the Planck 2018 data, such dark coupling can successfully solve the $H_0$ tension (with and without the presence of neutrinos). Concerning the properties of neutrinos, we find that $M_{\nu}$ may be wider than expected within the $\Lambda$CDM model and $N_{\rm eff}$ is fully compatible with three neutrino species (similar to $\Lambda$CDM prevision). The parameters characterizing the properties of neutrinos do not correlate with the coupling parameter of the interaction model. When considering the joint analysis of CMB from Planck 2018 and an estimate of $H_0$ from Hubble Space Telescope 2019 data, {\it we find an evidence for a non-null value of the coupling parameter at more than 3$\sigma$ confidence-level.} We also discuss the possible effects on the interacting scenario due to the inclusion of baryon acoustic oscillations data with Planck 2018. Our main results updating the dark sectors' interaction and neutrino properties in the model baseline, represent a new perspective in this direction. Clearly, a possible new physics in light of some dark interaction between dark energy and dark matter can serve as an alternative to $\Lambda$CDM scenario to explain the observable Universe, mainly in light of the current tension on $H_0$.

We consider the possibility that dark energy and baryons might scatter off each other. The type of interaction we consider leads to a pure momentum exchange, and does not affect the background evolution of the expansion history. We parametrize this interaction in an effective way at the level of Boltzmann equations. We compute the effect of dark energy-baryon scattering on cosmological observables, focusing on the Cosmic Microwave Background (CMB) temperature anisotropy power spectrum and the matter power spectrum. Surprisingly, we find that even huge dark energy-baryon cross-sections $\sigma_{xb} \sim {\cal O}({\rm b})$, which are generically excluded by non-cosmological probes such as collider searches or precision gravity tests, only leave an insignificant imprint on the observables considered. In the case of the CMB temperature power spectrum, the only imprint consists in a sub-percent enhancement or depletion of power (depending whether or not the dark energy equation of state lies above or below $-1$) at very low multipoles, which is thus swamped by cosmic variance. These effects are explained in terms of differences in how gravitational potentials decay in the presence of a dark energy-baryon scattering, which ultimately lead to an increase or decrease in the late-time integrated Sachs-Wolfe power. Even smaller related effects are imprinted on the matter power spectrum. The imprints on the CMB are not expected to be degenerate with the effects due to altering the dark energy sound speed. We conclude that, while strongly appealing, the prospects for a direct detection of dark energy through cosmology do not seem feasible when considering realistic dark energy-baryon cross-sections. As a caveat, our results hold to linear order in perturbation theory.

Non-linear electrodynamics (NLED) theories are well-motivated extensions of QED in the strong field regime, and have long been studied in the search for regular black hole (BH) solutions. We consider two well-studied and well-motivated NLED models coupled to General Relativity: the Euler-Heisenberg model and the Bronnikov model. After carefully accounting for the effective geometry induced by the NLED corrections, we determine the shadows of BHs within these two models. We then compare these to the shadow of the supermassive BH M87* recently imaged by the Event Horizon Telescope collaboration. In doing so, we are able to extract upper limits on the black hole magnetic charge, thus providing novel constraints on fundamental physics from this new extraordinary probe.

The Lambda-Cold Dark Matter (LCDM) model agrees with most of the cosmological observations, but has some hindrances from observed data at smaller scales such as galaxies. Recently, Berezhiani and Khoury (2015) proposed a new theory involving interacting superfluid dark matter with three model parameters, which explains galactic dynamics with great accuracy. In the present work, we study the cosmological behaviour of this model in the linear regime of cosmological perturbations. In particular, we compute both analytically and numerically the matter linear growth factor and obtain new bounds for the model parameters which are significantly stronger than previously found. These new constraints come from the fact that structures within the superfluid dark matter framework grow quicker than in LCDM, and quite rapidly when the DM-baryon interactions are strong.

We intend to understand cosmological structure formation within the framework of superfluid models of dark matter with finite temperatures. Of particular interest is the evolution of small-scale structures where the pressure and superfluid properties of the dark matter fluid are prominent. We compare the growth of structures in the superfluid dark matter with that of the standard cold dark matter paradigm and non-superfluid dark matter. The equations for superfluid hydrodynamics are computed numerically in an expanding $\Lambda$CDM background with spherical symmetry, and the effect of various superfluid fractions, temperatures, interactions, and masses on the collapse of structures is considered. The linear perturbation of the superfluid equations is derived, which gives further insight of the dynamics of the superfluid dark matter collapse. It is found that while a conventional dark matter fluid with self-interactions and finite temperatures experiences a suppression in the growth of structures at smaller scales, as expected due to the presence of pressure terms, a superfluid can collapse much more efficiently due to its ability to suppress the growth of entropy perturbations and thus gradients in the thermal pressure. It is also found that the cores of the dark matter halos initially become more superfluid during the collapse, but eventually reach a point where the superfluid fraction falls sharply. The formation of superfluid dark matter halos surrounded by a normal fluid dark matter background is therefore disfavored by the present work.

While the origin and composition of dark matter and dark energy remains unknown, it is possible that they might represent two manifestations of a single entity, as occurring in unified dark sector models. On the other hand, advances in our understanding of the dark sector of the Universe might arise from Cosmic Dawn, the epoch when the first stars formed. In particular, the first detection of the global 21-cm absorption signal at Cosmic Dawn from the EDGES experiment opens up a new arena wherein to test models of dark matter and dark energy. Here, we consider generalized and modified Chaplygin gas models as candidate unified dark sector models. We first constrain these models against Cosmic Microwave Background data from the \textit{Planck} satellite, before exploring how the inclusion of the global 21-cm signal measured by EDGES can improve limits on the model parameters, finding that the uncertainties on the parameters of the Chaplygin gas models can be reduced by a factor between $1.5$ and $10$. We also find that within the generalized Chaplygin gas model, the tension between the CMB and local determinations of the Hubble constant $H_0$ is reduced from $\approx 4\sigma$ to $\approx 1.3\sigma$. In conclusion, we find that the global 21-cm signal at Cosmic Dawn can provide an extraordinary window onto the physics of unified dark sectors.

Modified gravity (MG) theories aim to reproduce the observed acceleration of the Universe by reducing the dark sector while simultaneously recovering General Relativity (GR) within dense environments. Void studies appear to be a suitable scenario to search for imprints of alternative gravity models on cosmological scales. Voids cover an interesting range of density scales where screening mechanisms fade out, which reaches from a density contrast $\delta \approx -1$ close to their centers to $\delta \approx 0$ close to their boundaries. We present an analysis of the level of distinction between GR and two modified gravity theories, the Hu-Sawicki $f(R)$ and the symmetron theory. This study relies on the abundance, linear bias, and density profile of voids detected in n-body cosmological simulations. We define voids as connected regions made up of the union of spheres with a {\it \textup{mean}} density given by $\overline\rho_v=0.2\,\overline\rho_m$, but disconnected from any other voids. We find that the height of void walls is considerably affected by the gravitational theory, such that it increases for stronger gravity modifications. Finally, we show that at the level of dark matter n-body simulations, our constraints allow us to distinguish between GR and MG models with $|f_{R0}| > 10^{-6}$ and $z_{SSB} > 1$. Differences of best-fit values for MG parameters that are derived independently from multiple void probes may indicate an incorrect MG model. This serves as an important consistency check.

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

Since the very beginning of astronomy the location of objects on the sky has been a fundamental observational quantity that has been taken for granted. While precise two dimensional positional information is easy to obtain for observations in the electromagnetic spectrum, the positional accuracy of current and near future gravitational wave detectors is limited to between tens and hundreds of square degrees, which makes it extremely challenging to identify the host galaxies of gravitational wave events or to confidently detect any electromagnetic counterparts. Gravitational wave observations provide information on source properties and distances that is complementary to the information in any associated electromagnetic emission and that is very hard to obtain in any other way. Observing systems with multiple messengers thus has scientific potential much greater than the sum of its parts. A gravitational wave detector with higher angular resolution would significantly increase the prospects for finding the hosts of gravitational wave sources and triggering a multi-messenger follow-up campaign. An observatory with arcminute precision or better could be realised within the Voyage 2050 programme by creating a large baseline interferometer array in space and would have transformative scientific potential. Precise positional information of standard sirens would enable precision measurements of cosmological parameters and offer new insights on structure formation; a high angular resolution gravitational wave observatory would allow the detection of a stochastic background and resolution of the anisotropies within it; it would also allow the study of accretion processes around black holes; and it would have tremendous potential for tests of modified gravity and the discovery of physics beyond the Standard Model.

We consider a cosmological model with an interaction between dark matter and dark energy which leaves the background cosmology unaffected and only affects the evolution of the perturbations. This is achieved by introducing a coupling given in terms of the relative velocities of dark matter and dark energy. This interaction has the distinctive feature of appearing predominantly on small scales, where peculiar velocities can become important. We confront the predictions of the model to cosmological observations and find a potential alleviation of the known tension in the amplitude of density perturbations as measured by low redshift galaxy surveys and the Planck data. The model also predicts a shift in the turnover of the matter power spectrum which does not depend on the horizon at equality (fixed by the background cosmology and, thus, unaffected by the perturbations) and is entirely due to the interaction between dark matter and dark energy. A bias in the peculiar velocity between baryons and dark matter is also shown to be a unique feature of this type of interactions in the dark sector.

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

The detection of gravitational waves (GW) by the LIGO and Virgo collaborations offers a whole new range of possible tests and opens up a new window which may shed light on the nature of dark energy and dark matter. In the present work we investigate how future gravitational waves data could help to constrain different dynamical dark energy models. In particular, we perform cosmological forecastings of a class of well known and most used dynamical dark energy models using the third-generation gravitational wave detector, the Einstein Telescope. We have considered 1000 simulated GW events in order to constrain the parameter space of the dynamical dark energy models. Our analyses show that the inclusion of the GW data from the Einstein Telescope, significantly improves the parameter space of the dynamical dark energy models compared to their constraints extracted from the standard cosmological probes, namely, the cosmic microwave observations, baryon acoustic oscillations distance measurements, Supernove type Ia, and the Hubble parameter measurements.

We consider two stable Interacting Dark Matter -- Dark Energy models and confront them against current Cosmic Microwave Background data from the \textit{Planck} satellite. We then generate luminosity distance measurements from ${\cal O}(10^3)$ mock Gravitational Wave events matching the expected sensitivity of the proposed Einstein Telescope. We use these to forecast how the addition of Gravitational Wave standard sirens data can improve current limits on the Dark Matter -- Dark Energy coupling strength ($\xi$). We find that the addition of Gravitational Waves data can reduce the current uncertainty by a factor of $5$. Moreover, if the underlying cosmological model truly features Dark Matter -- Dark Energy interactions with a value of $\xi$ within the currently allowed $1\sigma$ upper limit, the addition of Gravitational Wave data would help disentangle such an interaction from the standard case of no interaction at a significance of more than $3\sigma$.

There are currently no reliable methods to measure transverse velocities of galaxies. This is an important piece of information that could allow us to probe the physics of structure formation as well as testing the underlying theory of gravity. The slingshot effect, a special case of the Integrated Sachs-Wolfe effect, is expected to create dipole signals in the temperature fluctuations of the Cosmic Microwave Background Radiation (CMB). This effect creates a hot spot behind and a cold spot in front of moving massive objects. The dipole signal created by the slingshot effect can be used to measure transverse velocities, but because the signal is expected to be weak, the effect has not been measured yet. The aim is to show that the slingshot effect can be measured by stacking the signals of galaxies falling into a collapsing cluster. We evaluate if the effect can probe modified gravity. We use data from a simulated galaxy catalogue (MDPL2) to mimic observations. We identify a massive galaxy cluster, and make maps of the slingshot effect around infalling galaxies. We add uncorrelated Gaussian noise to each map. The maps are rotated according to the direction to the cluster centre, such that the dipole signal will add up constructively when stacking. We compare each stack to a dipole stencil and we find the probability for a false positive in the absence of the slingshot signal. Each galaxy gives a signal of around $\Delta T/T \approx 10^{-9}$, while the precision of CMB experiments of today are $\Delta T/T \approx 4 \times 10^{-6}$. By stacking around 10 000 galaxies, the slingshot signal can be over the detectable threshold with experiments of today. However, future CMB experiments must be used to be certain of the strength of the observed signal.

The Quintessence model is one of the simplest and better known alternatives to Einstein's theory for gravity. The properties of the solutions have been studied in great detail in the background, linear and non-linear contexts in cosmology. Here we discuss new phenomenology that is induced by adding disformal terms to the interactions. Among other results, we show analytically and using cosmological simulations ran with the code \texttt{Isis} that the model posses a mechanism through which is it possible to obtain repulsive fifth forces, which are opposite to gravity. Although the equations are very complex, we also find that most of the new phenomenology can be explained by studying background quantities. We used our simulation data to test approximate relations that exist between the metric and scalar field perturbations as well as between the fifth force and gravity. Excellent agreement was found between exact and approximated solutions, which opens the way for running disformal gravity cosmological simulations using simply a Newtonian solver. These results could not only help us to find new ways of testing gravity, but also provide new motivations for building alternative models.

Searches for modified gravity in the large-scale structure try to detect the enhanced amplitude of density fluctuations caused by the fifth force present in many of these theories. Neutrinos, on the other hand, suppress structure growth below their free-streaming length. Both effects take place on comparable scales, and uncertainty in the neutrino mass leads to a degeneracy with modified gravity parameters for probes that are measuring the amplitude of the matter power spectrum. We explore the possibility to break the degeneracy between modified gravity and neutrino effects in the growth of structures by considering kinematic information related to either the growth rate on large scales or the virial velocities inside of collapsed structures. In order to study the degeneracy up to fully non-linear scales, we employ a suite of $N$-body simulations including both $f(R)$ modified gravity and massive neutrinos. Our results indicate that velocity information provides an excellent tool to distinguish massive neutrinos from modified gravity. Models with different values of neutrino masses and modified gravity parameters possessing a comparable matter power spectrum at a given time have different growth rates. This leaves imprints in the velocity divergence, which is therefore better suited than the amplitude of density fluctuations to tell the models apart. In such models with a power spectrum comparable to $\Lambda$CDM today, the growth rate is strictly enhanced. We also find the velocity dispersion of virialised clusters to be well suited to constrain deviations from general relativity without being affected by the uncertainty in the sum of neutrino masses.

In this paper, we study the chameleon profile in inhomogeneous density distributions and find that the fifth force in thin shell near the surface is weaker from what expected in homogeneous density distributions. Also, we check the validity of quasi-static approximation for the chameleon scalar field in the astrophysical time scales. We have investigated the rolling down behaviour of the scalar field on its effective potential inside a one solar mass red giant star by using MESA code. We have found that the scalar field is fast enough to follow the minimum of the potential. This adiabatic behaviour reduces the fifth force and extends the screened regions to lower densities where the field has smaller mass and was expected to be unscreened. As a consequence, the star evolution is similar to what expected from standard general relativity. In addition, considering the stability of star, an approximate constraint on the coupling constant $\beta$ is found.

The grand challenges of contemporary fundamental physics---dark matter, dark energy, vacuum energy, inflation and early universe cosmology, singularities and the hierarchy problem---all involve gravity as a key component. And of all gravitational phenomena, black holes stand out in their elegant simplicity, while harbouring some of the most remarkable predictions of General Relativity: event horizons, singularities and ergoregions. The hitherto invisible landscape of the gravitational Universe is being unveiled before our eyes: the historical direct detection of gravitational waves by the LIGO-Virgo collaboration marks the dawn of a new era of scientific exploration. Gravitational-wave astronomy will allow us to test models of black hole formation, growth and evolution, as well as models of gravitational-wave generation and propagation. It will provide evidence for event horizons and ergoregions, test the theory of General Relativity itself, and may reveal the existence of new fundamental fields. The synthesis of these results has the potential to radically reshape our understanding of the cosmos and of the laws of Nature. The purpose of this work is to present a concise, yet comprehensive overview of the state of the art in the relevant fields of research, summarize important open problems, and lay out a roadmap for future progress.

We reconsider the dynamics of the Universe in the presence of interactions in the cosmological dark sector. A class of interacting models is introduced via a real function $f\left(r\right)$ of the ratio $r$ between the energy densities of the (pressureless) cold dark matter (CDM) and dark energy (DE). The subclass of models for which the ratio $r$ depends only on the scale factor is shown to be equivalent to unified models of the dark sector, i.e. models for which the CDM and DE components can be combined in order to form a unified dark fluid. For specific choices of the function $f\left(r\right)$ we recover several models already studied in the literature. We analyse various special cases of this type of interacting models using a suitably modified version of the CLASS code combined with MontePython in order to constrain the parameter space with the data from supernova of type SNe Ia (JLA), the Hubble constant $H_{0}$, cosmic chronometers (CC), baryon acoustic oscilations (BAO) and data from the Planck satellite (Planck TT). Our analysis shows that even if data from the late Universe ($H_{0}$, SNe Ia and CC) indicate an interaction in the dark sector, the data related to the early Universe (BAO and Planck TT) constrain this interaction substantially, in particular for cases in which the background dynamics is strongly affected.

We study a novel interacting dark energy $-$ dark matter scenario where the anisotropic stress of the large scale inhomogeneities is considered. The dark energy has a constant equation of state and the interaction model produces stable perturbations. The resulting picture is constrained using different astronomical data aiming to measure the impact of the anisotropic stress on the cosmological parameters. Our analyses show that a non-zero interaction in the dark sector is allowed while a non-interaction scenario is recovered within 68\% CL. The anisotropic stress is also constrained to be small, and its zero value is permitted within 68\% CL. The dark energy equation of state, $w_x$, is also found to be close to `$-1$' boundary. However, from the ratio of the CMB TT spectra, we see that the model has a mild deviation from the $\Lambda$CDM cosmology while such deviation is almost forbidden from the CMB TT spectra alone. Although the deviation is not much significant, but from the present data, we cannot exclude such deviation. Overall, at the background level, the model is close to the $\Lambda$CDM cosmology while at the level of perturbations, a non-zero but a very small interaction in the dark sector is permitted. Perhaps, a more accurate conclusion can be made with the next generation of surveys. We also found that the region $w_x < -1$, is found to be effective to release the tension on $H_0$. Finally, from the Bayesian analysis, we find that $\Lambda$CDM remains in still preferred over the interacting scenarios.

We present a formalism to study screening mechanisms in modified theories of gravity via perturbative methods in different cosmological scenarios. We consider Einstein frame posed theories that are recast as Jordan frame theories, where a known formalism is employed, though the resulting non-linearities of the Klein-Gordon equation acquire an explicit coupling between matter and the scalar field, which is not present in Jordan frame theories. The obtained growth functions are then separated in screening and non-screened contributions to facilitate its analysis. This allows us to compare several theoretical models and to recognize patterns which can be used to differentiate models and their screening mechanisms. In particular, we find anti-screening features in the Symmetron model. In opposition, chameleon type theories, both in the Jordan and in the Einstein frame, always present a screening behaviour. Up to third order in perturbation, we find no anti-screening behaviour in theories with a Vainshtein mechanism, such as the DGP and the cubic Galileon.

Multiple modifications of general relativity (GR) have been proposed in the literature in order to understand the nature of the accelerated expansion of the Universe. However, thus far all the predictions of GR have been confirmed with constantly increasing accuracy. In this work, we study the imprints of a particular class of models -- "screened" modified gravity theories -- on cosmic filaments. We have utilized the $N$-body code ISIS/RAMSES to simulate the symmetron model and the Hu-Sawicky $f(R)$ model, and we post-process the output with DisPerSE to identify the filaments of the cosmic web. We investigated how the global properties of the filaments -- such as their lengths, masses, and thicknesses -- as well as their radial density and speed profiles change under different gravity theories. We find that filaments are, on average, shorter and denser in modified gravity models compared to in $\Lambda$CDM. We also find that the speed profiles of the filaments are enhanced, consistent with theoretical expectations. Overall, our results suggest that cosmic filaments can be an effective complementary probe of screened modified gravity theories on Mpc scales.

We investigate the observational consequences of a novel class of stable interacting dark energy (IDE) models, featuring interactions between dark matter (DM) and dark energy (DE). In the first part of our work, we start by considering two IDE models which are known to present early-time linear perturbation instabilities. Applying a transformation depending on the dark energy equation of state (EoS) to the DM-DE coupling, we then obtain two novel stable IDE models. Subsequently, we derive robust and accurate constraints on the parameters of these models, assuming a constant EoS $w_x$ for the DE fluid, in light of some of the most recent publicly available cosmological data. These include Cosmic Microwave Background (CMB) temperature and polarization anisotropy measurements from the \textit{Planck} satellite, a selection of Baryon Acoustic Oscillation measurements, Supernovae Type-Ia luminosity distance measurements from the JLA sample, and measurements of the Hubble parameter up to redshift $2$ from cosmic chronometers. Our analysis displays a mild preference for the DE fluid residing in the phantom region ($w_x<-1$), with significance up to 95\% confidence level, while we obtain new upper limits on the coupling parameter between the dark components. The preference for a phantom DE suggests a coupling function $Q<0$, thus a scenario where energy flows from the DE to the DM. We also examine the possibility of addressing the $H_0$ and $\sigma_8$ tensions, finding that only the former can be partially alleviated. Finally, we perform a Bayesian model comparison analysis to quantify the possible preference for the two IDE models against the standard concordance $\Lambda$CDM model, finding that the latter is always preferred with the strength of the evidence ranging from positive to very strong.

We investigate the effects of a disformal coupling between dark energy and dark matter in the predictions of the spherical collapse and its signatures in galaxy cluster number counts. We find that the disformal coupling has no significant effects on spherical collapse at high redshifts, and in particular during matter domination epoch. However, at lower redshifts, the extrapolated linear density contrast at collapse close to redshift $z \lesssim 1$ and overdensity at virialization can be strongly suppressed by a disformal coupling between dark energy and dark matter. We also find that disformal coupling can have different imprints on cluster number counts compared with conformal coupling, such that the disformal coupling can strongly suppress the predicted number of clusters per redshift interval at $z > 0.1$ while enhance the number of cluster at $z < 0.05$. Using the specifications of eROSITA survey, we find that the disformal coupling between dark energy and dark matter can be tightly constrained by cluster number counts.

In order to determine the observable signatures of modified gravity theories, it is important to consider the effect of baryonic physics. We use a modified version of the ISIS code to run cosmological hydrodynamic simulations to study degeneracies between modified gravity and radiative hydrodynamical processes. Of these, one was the standard $\Lambda$CDM model and four were variations of the Symmetron model. For each model we ran three variations of baryonic processes: non-radiative hydrodynamics; cooling and star formation; and cooling, star formation, and supernova feedback. We construct stacked gas density, temperature, and dark matter density profiles of the halos in the simulations, and study the differences between them. We find that both radiative variations of the models show degeneracies between their processes and at least two of the three parameters defining the Symmetron model.

We investigate the efficiency of screening mechanisms in the hybrid metric-Palatini gravity. The value of the field is computed around spherical bodies embedded in a background of constant density. We find a thin shell condition for the field depending on the background field value. In order to quantify how the thin shell effect is relevant, we analyze how it behaves in the neighborhood of different astrophysical objects (planets, moons or stars). We find that the condition is very well satisfied except only for some peculiar objects. Furthermore we establish bounds on the model using data from solar system experiments such as the spectral deviation measured by the Cassini mission and the stability of the Earth-Moon system, which gives the best constraint to date on $f(R)$ theories. These bounds contribute to fix the range of viable hybrid gravity models.

Consistency relations of large-scale structures provide exact nonperturbative results for cross-correlations of cosmic fields in the squeezed limit. They only depend on the equivalence principle and the assumption of Gaussian initial conditions, and remain nonzero at equal times for cross-correlations of density fields with velocity or momentum fields, or with the time derivative of density fields. We show how to apply these relations to observational probes that involve the integrated Sachs-Wolfe effect or the kinematic Sunyaev-Zeldovich effect. In the squeezed limit, this allows us to express the three-point cross-correlations, or bispectra, of two galaxy or matter density fields, or weak lensing convergence fields, with the secondary Cosmic Microwave Background (CMB) distortion in terms of products of a linear and a nonlinear power spectrum. In particular, we find that cross-correlations with the integrated Sachs-Wolfe effect show a specific angular dependence. These results could be used to test the equivalence principle and the primordial Gaussianity, or to check the modeling of large-scale structures.

Stability analysis of interacting dark energy models generally divides its parameters space into two regions: (i) $w_x \geq -1$ and $\xi \geq 0$ and (ii) $w_x \leq -1$ and $\xi \leq 0$, where $w_x$ is the dark energy equation of state and $\xi$ is the coupling strength of the interaction. Due to this separation, crucial information about the cosmology and phenomenology of these models may be lost. In a recent study it has been shown that one can unify the two regions with a coupling function which depends on the dark energy equation of state. In this work we introduce a new coupling function which also unifies the two regions of the parameter space and generalises the previous proposal. We analyse this scenario considering the equation of state of DE to be either constant or dynamical. We study the cosmology of such models and constrain both scenarios with the use of latest astronomical data from both background evolution as well as large scale structures. Our analysis shows that a non-zero value of the coupling parameter $\xi$ as well as the dark energy equation of state other than `$-1$' are allowed. However, within $1\sigma$ confidence level, $\xi = 0$, and the dark energy equation of state equal to `$-1$' are compatible with the current data. In other words, the observational data allow a very small but nonzero deviation from the $\Lambda$-cosmology, however, within $1\sigma$ confidence-region the interacting models can mimick the $\Lambda$-cosmology. In fact we observe that the models both at background and perturbative levels are very hard to distinguish form each other and from $\Lambda$-cosmology as well. Finally, we offer a rigorous analysis on the current tension on $H_0$ allowing different regions of the dark energy equation of state which shows that interacting dark energy models reasonably solve the current tension on $H_0$.

We study a theory of gravity in which the action is a result from the general purely disformal transformation on the Einstein-Hilbert action. This theory is a sub-class of GLPV theory which is the the generalization of covariant Galileon. Nevertheless, we find that the self accelerating solution for the background universe disappears in this theory. We also find that, for this theory, the Vainshtein mechanism is absent. However, the Vainshtein mechanism is not necessary for this theory, because this theory can nearly mimic the Einstein theory of gravity at all scales inside the Huble radius without this mechanism.

We put constraints on dark energy properties using the PADE parameterisation, and compare it to the same constraints using Chevalier-Polarski-Linder (CPL) and $\Lambda$CDM, at both the background and the perturbation levels. The dark energy equation of state parameter of the models is derived following the mathematical treatment of PADE expansion. Unlike CPL parameterisation, the PADE approximation provides different forms of the equation of state parameter which avoid the divergence in the far future. Initially, we perform a likelihood analysis in order to put constraints on the model parameters using solely background expansion data and we find that all parameterisations are consistent with each other. Then, combining the expansion and the growth rate data we test the viability of PADE parameterisations and compare them with CPL and $\Lambda$CDM models respectively. Specifically, we find that the growth rate of the current PADE parameterisations is lower than $\Lambda$CDM model at low redshifts, while the differences among the models are negligible at high redshifts. In this context, we provide for the first time growth index of linear matter perturbations in PADE cosmologies. Considering that dark energy is homogeneous we recover the well known asymptotic value of the growth index, namely $\gamma_{\infty}=\frac{3(w_{\infty}-1)}{6w_{\infty}-5}$, while in the case of clustered dark energy we obtain $\gamma_{\infty}\simeq \frac{3w_{\infty}(3w_{\infty}-5)}{(6w_{\infty}-5)(3w_{\infty}-1)}$. Finally, we generalize the growth index analysis in the case where $\gamma$ is allowed to vary with redshift and we find that the form of $\gamma(z)$ in PADE parameterisation extends that of the CPL and $\Lambda$CDM cosmologies respectively.

We investigate how three different possibilities of neutrino mass hierarchies, namely normal, inverted, and degenerate, can affect the observational constraints on three well known dynamical dark energy models, namely the Chevallier-Polarski-Linder, logarithmic, and the Jassal-Bagla-Padmanabhan parametrizations. In order to impose the observational constraints on the models, we performed a robust analysis using Planck 2015 temperature and polarization data, Supernovae type Ia from Joint Light curve analysis, baryon acoustic oscillations distance measurements, redshift space distortion characterized by $f(z)\sigma_8(z)$ data, weak gravitational lensing data from Canada-France-Hawaii Telescope Lensing Survey, and cosmic chronometers data plus the local value of the Hubble parameter. We find that different neutrino mass hierarchies return similar fit on almost all model parameters and mildly change the dynamical dark energy properties.

We propose a new method to test modified gravity theories, taking advantage of the available data on extrasolar planets. We computed the deviations from the Kepler third law and use that to constrain gravity theories beyond General Relativity. We investigate gravity models which incorporate three screening mechanisms: the Chameleon, the Symmetron and the Vainshtein. We find that data from exoplanets orbits are very sensitive to the screening mechanisms putting strong constraints in the parameter space for the Chameleon models and the Symmetron, complementary and competitive to other methods, like interferometers and solar system. With the constraints on Vainshtein we are able to work beyond the hypothesis that the crossover scale is of the same order of magnitude than the Hubble radius $r_c \sim H_0^{-1}$, which makes the screening work automatically, testing how strong this hypothesis is and the viability of other scales.

We study the properties of dark matter haloes in a wide range of modified gravity models, namely, $f(R)$, DGP, and interacting dark energy models. We study the effects of modified gravity and dark energy on the internal properties of haloes, such as the spin and the structural parameters. We find that $f(R)$ gravity enhance the median value of the Bullock spin parameter, but could not detect such effects for DGP and coupled dark energy. $f(R)$ also yields a lower median sphericity and oblateness, while coupled dark energy has the opposite effect. However, these effects are very small. We then study the interaction rate of haloes in different gravity, and find that only strongly coupled dark energy models enhance the interaction rate. We then quantify the enhancement of the alignment of the spins of interacting halo pairs by modified gravity. Finally, we study the alignment of the major axes of haloes with the large-scale structures. The alignment of the spins of interacting pairs of haloes in DGP and coupled dark energy models show no discrepancy with GR, while $f(R)$ shows a weaker alignment. Strongly coupled dark energy shows a stronger alignment of the halo shape with the large-scale structures.

Extending general relativity by adding extra degrees of freedom is a popular approach for explaining the accelerated expansion of the Universe and to build high energy completions of the theory of gravity. The presence of such new degrees of freedom is, however, tightly constrained from several observations and experiments that aim to test general relativity in a wide range of scales. The viability of a given modified theory of gravity, therefore, strongly depends on the existence of a screening mechanism that suppresses the extra degrees of freedom. We perform simulations, and find that waves propagating in the new degrees of freedom can significantly impact the efficiency of some screening mechanisms, thereby threatening the viability of these modified gravity theories. Specifically, we show that the waves produced in the symmetron model can increase the amplitude of the fifth force and the parametrized post-Newtonian parameters by several orders of magnitude.

We present exact kinematic consistency relations for cosmological structures that do not vanish at equal times and can thus be measured in surveys. These rely on cross-correlations between the density and velocity, or momentum, fields. Indeed, the uniform transport of small-scale structures by long wavelength modes, which cannot be detected at equal times by looking at density correlations only, gives rise to a shift in the amplitude of the velocity field that could be measured. These consistency relations only rely on the weak equivalence principle and Gaussian initial conditions. They remain valid in the non-linear regime and for biased galaxy fields. They can be used to constrain non-standard cosmological scenarios or the large-scale galaxy bias.

Yes, but only for a parameter value that makes it almost coincide with the standard model. We reconsider the cosmological dynamics of a generalized Chaplygin gas (gCg) which is split into a cold dark matter (CDM) part and a dark energy (DE) component with constant equation of state. This model, which implies a specific interaction between CDM and DE, has a $\Lambda$CDM limit and provides the basis for studying deviations from the latter. Including matter and radiation, we use the (modified) CLASS code \cite{class} to construct the CMB and matter power spectra in order to search for a gCg-based concordance model that is in agreement with the SNIa data from the JLA sample and with recent Planck data. The results reveal that the gCg parameter $\alpha$ is restricted to $|\alpha|\lesssim 0.05$, i.e., to values very close to the $\Lambda$CDM limit $\alpha =0$. This excludes, in particular, models in which DE decays linearly with the Hubble rate.

In this work we will test an alternative model of gravity belonging to the large family of galileon models. It is characterized by an intrinsic breaking of the Vainshtein mechanism inside large astrophysical objects, thus having possibly detectable observational signatures. We will compare theoretical predictions from this model with the observed total mass profile for a sample of clusters of galaxies. The profiles are derived using two complementary tools: X-ray hot intra-cluster gas dynamics, and strong and weak gravitational lensing. We find that a dependence with the dynamical internal status of each cluster is possible; for those clusters which are very close to be relaxed, and thus less perturbed by possible astrophysical local processes, the galileon model gives a quite good fit to both X-ray and lensing observations. Both masses and concentrations for the dark matter halos are consistent with earlier results found in numerical simulations and in the literature, and no compelling statistical evidence for a deviation from general relativity is detectable from the present observational state. Actually, the characteristic galileon parameter $\Upsilon$ is always consistent with zero, and only an upper limit ($\lesssim0.086$ at $1\sigma$, $\lesssim0.16$ at $2\sigma$, and $\lesssim0.23$ at $3\sigma$) can be established. Some interesting distinctive deviations might be operative, but the statistical validity of the results is far from strong, and better data would be needed in order to either confirm or reject a potential tension with general relativity.

In this work we study the dynamics of the Local Group (LG) within the context of cosmological models beyond General Relativity (GR). Using observable kinematic quantities to identify candidate pairs we build up samples of simulated LG-like objects drawing from $f(R)$, symmetron, DGP and quintessence N-body simulations together with their $\Lambda$CDM counterparts featuring the same initial random phase realisations. The variables and intervals used to define LG-like objects are referred to as Local Group model; different models are used throughout this work and adapted to study their dynamical and kinematic properties. The aim is to determine how well the observed LG-dynamics can be reproduced within cosmological theories beyond GR, We compute kinematic properties of samples drawn from alternative theories and $\Lambda$CDM and compare them to actual observations of the LG mass, velocity and position. As a consequence of the additional pull, pairwise tangential and radial velocities are enhanced in modified gravity and coupled dark energy with respect to $\Lambda$CDM, inducing significant changes to the total angular momentum and energy of the LG. For example, in models such as $f(R)$ and the symmetron this increase can be as large as $60\%$, peaking well outside of the $95\%$ confidence region allowed by the data. This shows how simple considerations about the LG dynamics can lead to clear small-scale observational signatures for alternative scenarios, without the need of expensive high-resolution simulations.

We use cosmological hydrodynamical simulations to study the effect of screened modified gravity models on the mass estimates of galaxy clusters. In particular, we focus on two novel aspects: (i) we study modified gravity models in which baryons and dark matter are coupled with different strengths to the scalar field, and, (ii) we put the simulation results into the greater context of a general screened-modified gravity parametrization. We have compared the mass of clusters inferred via lensing versus the mass inferred via kinematical measurements as a probe of violations of the equivalence principle at Mpc scales. We find that estimates of cluster masses via X-ray observations is mainly sensitive to the coupling between the scalar degree of freedom and baryons - while the kinematical mass is mainly sensitive to the coupling to dark matter. Therefore, the relation between the two mass estimates is a probe of a possible non-universal coupling between the scalar field, the standard model fields, and dark matter. Finally, we used observational data of kinetic, thermal and lensing masses to place constraints on deviations from general relativity on cluster scales for a general parametrization of screened modified gravity theories which contains $f(R)$ and Symmetron models. We find that while the kinematic mass can be used to place competitive constraints, using thermal measurements is challenging as a potential non-thermal contribution is degenerate with the imprint of modified gravity.

Modified gravity theories with a screening mechanism have acquired much interest recently in the quest for a viable alternative to General Relativity on cosmological scales, given their intrinsic property of being able to pass Solar System scale tests and, at the same time, to possibly drive universe acceleration on much larger scales. Here, we explore the possibility that the same screening mechanism, or its breaking at a certain astrophysical scale, might be responsible of those gravitational effects which, in the context of general relativity, are generally attributed to Dark Matter. We consider a recently proposed extension of covariant Galileon models in the so-called "beyond Horndeski" scenario, where a breaking of the Vainshtein mechanism is possible and, thus, some peculiar observational signatures should be detectable and make it distinguishable from general relativity. We apply this model to a sample of clusters of galaxies observed under the \textit{CLASH} survey, using both new data from gravitational lensing events and archival data from X-ray intra-cluster hot gas observations. In particular, we use the latter to model the gas density, and then use it as the only ingredient in the matter clusters' budget to calculate the expected lensing convergence map. Results show that, in the context of this extended Galileon, the assumption of having only gas and no Dark Matter at all in the clusters is able to match observations. We also obtain narrow and very interesting bounds on the parameters which characterize this model. In particular, we find that, at least for one of them, the general relativity limit is excluded at $2\sigma$ confidence level, thus making this model clearly statistically different and competitive with respect to general relativity.

Alternative theories of gravity typically invoke an environment-dependent screening mechanism to allow phenomenologically interesting deviations from general relativity (GR) to manifest on larger scales, while reducing to GR on small scales. The observation of the transition from screened to unscreened behavior would be compelling evidence for beyond-GR physics. We show that pairwise peculiar velocity statistics, in particular the relative radial velocity dispersion, $\sigma_\parallel$ , can be used to observe this transition when they are binned by some measure of halo environment. We established this by measuring the radial velocity dispersion between pairs of halos in N-body simulations for three $f(R)$ gravity and four symmetron models. We developed an estimator involving only line-of-sight velocities to show that this quantity is observable, and binned the results in halo mass, ambient density, and the isolatedness of halos. Ambient density is found to be the most relevant measure of environment; it is distinct from isolatedness, and correlates well with theoretical expectations for the symmetron model. By binning $\sigma_\parallel$ in ambient density, we find a strong environment-dependent signature for the symmetron models, with the velocities showing a clear transition from GR to non-GR behavior. No such transition is observed for $f(R)$, as the relevant scales are deep in the unscreened regime. Observations of the relative radial velocity dispersion in forthcoming peculiar velocity surveys, if binned appropriately by environment, therefore offer a valuable way of detecting the screening signature of modified gravity.

We use a spherical model and an extended excursion set formalism with drifting diffusive barriers to predict the abundance of cosmic voids in the context of general relativity as well as f(R) and symmetron models of modified gravity. We detect spherical voids from a suite of N-body simulations of these gravity theories and compare the measured void abundance to theory predictions. We find that our model correctly describes the abundance of both dark matter and galaxy voids, providing a better fit than previous proposals in the literature based on static barriers. We use the simulation abundance results to fit for the abundance model free parameters as a function of modified gravity parameters, and show that counts of dark matter voids can provide interesting constraints on modified gravity. For galaxy voids, more closely related to optical observations, we find that constraining modified gravity from void abundance alone may be significantly more challenging. In the context of current and upcoming galaxy surveys, the combination of void and halo statistics including their abundances, profiles and correlations should be effective in distinguishing modified gravity models that display different screening mechanisms.

We present measurements of the spatial clustering statistics in redshift space of various scalar field modified gravity simulations. We utilise the two-point and the three-point correlation functions to quantify the spatial distribution of dark matter halos within these simulations and thus discern between the models. We compare $\Lambda$CDM simulations to various modified gravity scenarios and find consistency with previous work in terms of 2-point statistics in real and redshift-space. However using higher order statistics such as the three-point correlation function in redshift space we find significant deviations from $\Lambda$CDM hinting that higher order statistics may prove to be a useful tool in the hunt for deviations from General Relativity.

Euclid is a European Space Agency medium class mission selected for launch in 2020 within the Cosmic Vision 2015 2025 program. The main goal of Euclid is to understand the origin of the accelerated expansion of the universe. Euclid will explore the expansion history of the universe and the evolution of cosmic structures by measuring shapes and redshifts of galaxies as well as the distribution of clusters of galaxies over a large fraction of the sky. Although the main driver for Euclid is the nature of dark energy, Euclid science covers a vast range of topics, from cosmology to galaxy evolution to planetary research. In this review we focus on cosmology and fundamental physics, with a strong emphasis on science beyond the current standard models. We discuss five broad topics: dark energy and modified gravity, dark matter, initial conditions, basic assumptions and questions of methodology in the data analysis. This review has been planned and carried out within Euclid's Theory Working Group and is meant to provide a guide to the scientific themes that will underlie the activity of the group during the preparation of the Euclid mission.

We propose that the mass-temperature relation of galaxy clusters is a prime candidate for testing gravity theories beyond Einstein's general relativity, for modified gravity models with universal coupling between matter and the scalar field. For non-universally coupled models we discover that the impact of modified gravity can remain hidden from the mass-temperature relation. Using cosmological simulations, we find that in modified gravity the mass-temperature relation varies significantly from the standard gravity prediction of $M \propto T^{1.73}$. To be specific, for symmetron models with a coupling factor of $\beta=1$ we find a lower limit to the power law as $M\propto T^{1.6}$; and for f(R) gravity we compute predictions based on the model parameters. We show that the mass-temperature relation, for screened modified gravities, is significantly different from that of standard gravity for the less massive and colder galaxy clusters, while being indistinguishable from Einstein's gravity at massive, hot galaxy clusters. We further investigate the mass-temperature relation for other mass estimates than the thermal mass estimate, and discover that the gas mass-temperature results show an even more significant deviations from Einstein's gravity than the thermal mass-temperature.

We have investigated structure formation in the $\gamma$ gravity $f(R)$ model with {\it N}-body simulations. The $\gamma$ gravity model is a proposal which, unlike other viable $f(R)$ models, not only changes the gravitational dynamics, but can in principle also have signatures at the background level that are different from those obtained in $\Lambda$CDM (Cosmological constant, Cold Dark Matter). The aim of this paper is to study the nonlinear regime of the model in the case where, at late times, the background differs from $\Lambda$CDM. We quantify the signatures produced on the power spectrum, the halo mass function, and the density and velocity profiles. To appreciate the features of the model, we have compared it to $\Lambda$CDM and the Hu-Sawicki $f(R)$ models. For the considered set of parameters we find that the screening mechanism is ineffective, which gives rise to deviations in the halo mass function that disagree with observations. This does not rule out the model per se, but requires choices of parameters such that $|f_{R0}|$ is much smaller, which would imply that its cosmic expansion history cannot be distinguished from $\Lambda$CDM at the background level.

We investigate the propagation of scalar waves induced by matter sources in the context of scalar-tensor theories of gravity which include screening mechanisms for the scalar degree of freedom. The usual approach when studying these theories in the non-linear regime of cosmological perturbations is based on the assumption that scalar waves travel at the speed of light. Within General Relativity such approximation is good and leads to no loss of accuracy in the estimation of observables. We find, however, that mass terms and non-linearities in the equations of motion lead to propagation and dispersion velocities significantly different from the speed of light. As the group velocity is the one associated to the propagation of signals, a reduction of its value has direct impact on the behavior and dynamics of nonlinear structures within modified gravity theories with screening. For instance, the internal dynamics of galaxies and satellites submerged in large dark matter halos could be affected by the fact that the group velocity is smaller than the speed of light. It is therefore important, within such framework, to take into account the fact that different part of a galaxy will see changes in the environment at different times. A full non-static analysis may be necessary under those conditions.

We investigate statistical properties of the distribution of matter at redshift zero in disformal gravity by using N-body simulations. The disformal model studied here consists of a conformally coupled symmetron field with an additional exponential disformal term. We conduct cosmological simulations to discover the impact of the new disformal terms in the matter power spectrum, halo mass function, and radial profile of the scalar field. We calculated the disformal geodesic equation and the equation of motion for the scalar field. We then implemented these equations into the N-body code ISIS, which is a modified gravity version of the code RAMSES. The presence of a conformal symmetron field increases both the power spectrum and mass function compared to standard gravity on small scales. Our main finding is that the newly added disformal terms tend to counteract these effects and can make the evolution slightly closer to standard gravity. We finally show that the disformal terms give rise to oscillations of the scalar field in the centre of the dark matter haloes.

In this paper we study the effects of letting the dark matter and the gas in the Universe couple to the scalar field of the symmetron model, a modified gravity theory, with varying coupling strength. We also search for a way to distinguish between universal and non-universal couplings in observations. The research is performed utilising a series of hydrodynamic, cosmological N-Body simulations, studying the resulting power spectra and galaxy halo properties, such as the density and temperature profiles. Results show that in the cases of universal couplings, the deviations in the baryon fraction from $\Lambda$CDM are smaller than in the cases of non-universal couplings throughout the halos. The same is apparent in the power spectrum baryon bias, defined as the ratio of gas to dark matter power spectrum. Deviations of the density profiles and power spectra from the $\Lambda$CDM reference values can differ significantly between dark matter and gas because the dark matter deviations are mostly larger than the deviations in the gas.

Euclid is a European Space Agency medium class mission selected for launch in 2019 within the Cosmic Vision 2015-2025 programme. The main goal of Euclid is to understand the origin of the accelerated expansion of the Universe. Euclid will explore the expansion history of the Universe and the evolution of cosmic structures by measuring shapes and redshifts of galaxies as well as the distribution of clusters of galaxies over a large fraction of the sky. Although the main driver for Euclid is the nature of dark energy, Euclid science covers a vast range of topics, from cosmology to galaxy evolution to planetary research. In this review we focus on cosmology and fundamental physics, with a strong emphasis on science beyond the current standard models. We discuss five broad topics: dark energy and modified gravity, dark matter, initial conditions, basic assumptions and questions of methodology in the data analysis. This review has been planned and carried out within Euclid's Theory Working Group and is meant to provide a guide to the scientific themes that will underlie the activity of the group during the preparation of the Euclid mission.

Modified gravity models require a screening mechanism to be able to evade the stringent constraints from local gravity experiments and, at the same time, give rise to observable astrophysical and cosmological signatures. Such screened modified gravity models necessarily have dynamics determined by complex nonlinear equations that usually need to be solved on a model-by-model basis to produce predictions. This makes testing them a cumbersome process. In this paper, we investigate whether there is a common signature for all the different models that is suitable to testing them on cluster scales. To do this we propose an observable related to the fifth force, which can be observationally related to the ratio of dynamical-to-lensing mass of a halo, and then show that the predictions for this observable can be rescaled to a near universal form for a large class of modified gravity models. We demonstrate this using the Hu-Sawicki $f(R)$, the Symmetron, the nDGP, and the Dilaton models, as well as unifying parametrizations. The universal form is determined by only three quantities: a strength, a mass, and a width parameter. We also show how these parameters can be derived from a specific theory. This self-similarity in the predictions can hopefully be used to search for signatures of modified gravity on cluster scales in a model-independent way.

Self-consistent ${\it N}$-body simulations of modified gravity models are a key ingredient to obtain rigorous constraints on deviations from General Relativity using large-scale structure observations. This paper provides the first detailed comparison of the results of different ${\it N}$-body codes for the $f(R)$, DGP, and Symmetron models, starting from the same initial conditions. We find that the fractional deviation of the matter power spectrum from $\Lambda$CDM agrees to better than $1\%$ up to $k \sim 5-10~h/{\rm Mpc}$ between the different codes. These codes are thus able to meet the stringent accuracy requirements of upcoming observational surveys. All codes are also in good agreement in their results for the velocity divergence power spectrum, halo abundances and halo profiles. We also test the quasi-static limit, which is employed in most modified gravity ${\it N}$-body codes, for the Symmetron model for which the most significant non-static effects among the models considered are expected. We conclude that this limit is a very good approximation for all of the observables considered here.

In this paper we present the first results from implementing two scalar-tensor modified gravity theories, the symmetron and the Hu-Sawicki $f(R)$-gravity model, into a hydrodynamic N-body code with dark matter particles and a baryonic ideal gas. The study is a continuation of previous work where the symmetron and $f(R)$ have been successfully implemented in the RAMSES code, but for dark matter only. By running simulations, we show that the deviation from $\Lambda$CDM in these models for the gas density profiles are significantly lower than the dark matter equivalents. When it comes to the matter power-spectrum we find that hydrodynamic simulations agree very well with dark matter only simulations as long as we consider scales larger than $k\sim 0.5$ h/Mpc. In general the effects of modified gravity on the baryonic gas is found to not always mirror the effects it has on the dark matter. The largest signature is found when considering temperature profiles. We find that the gas temperatures in the modified gravity model studied here show deviations, when compared to $\Lambda$CDM, that can be a factor of a few larger than the deviations found in density profiles and power spectra.

Screened modified gravity predicts potentially large signatures in the peculiar velocity field that makes it an interesting probe to test gravity on cosmological scales. We investigate the signatures induced by the Symmetron and a Chameleon $f(R)$ model in the peculiar velocity field using $N$-body simulations. By studying fifth force and halo velocity profiles we identify three general categories of effects found in screened modified gravity models: a fully screened regime where we recover $\Lambda$CDM to high precision, an unscreened regime where the fifth force is in full operation, and, a partially screened regime where screening occurs in the inner part of a halo, but the fifth force is active at larger radii. These three regimes can be pointed out very clearly by analyzing the deviation in the maximum cluster velocity. Observationally, the partially screened regime is of particular interest since an uniform increase of the gravitational force - as present in the unscreened regime - is degenerate with the (dynamical) halo mass estimate, and, thus, hard to detect.

Scalar fields are crucial components in high energy physics and extensions of General Relativity. The fact they are not observed in the solar system may be due to a mechanism which screens their presence in high dense regions. We show how observations of the ellipticity of galaxy clusters can discriminate between models with and without scalar fields and even between different screening mechanisms. Using nowadays X-ray observations we put novel constraints on the different models.

Several extensions of General Relativity and high energy physics include scalar fields as extra degrees of freedom. In the search for predictions in the non-linear regime of cosmological evolution, the community makes use of numerical simulations in which the quasi-static limit is assumed when solving the equation of motion of the scalar field. In this Letter, we propose a method to solve the full equations of motion for scalar degrees of freedom coupled to matter. We run cosmological simulations which track the full time and space evolution of the scalar field, and find striking differences with respect to the commonly used quasi-static approximation. This novel procedure reveals new physical properties of the scalar field and uncovers concealed astrophysical phenomena which were hidden in the old approach.

Several extensions of the standard cosmological model include scalar fields as new degrees of freedom in the underlying gravitational theory. A particular class of these scalar field theories include screening mechanisms intended to hide the scalar field below observational limits in the solar system, but not on galactic scales, where data still gives freedom to find possible signatures of their presence. In order to make predictions to compare with observations coming from galactic and clusters scales (i.e. in the non-linear regime of cosmological evolution), cosmological N-body simulations are needed, for which codes that can solve for the scalar field must be developed. We present a new implementation of scalar-tensor theories of gravity which include screening mechanisms. The code is based in the already existing code RAMSES, to which we have added a non-linear multigrid solver that can treat a large class of scalar tensor theories of modified gravity. We present details of the implementation and the tests that we made to it. As application of the new code, we have studied the influence that two particular modified gravity theories, the symmetron and $f(R)$ gravity, have on the shape of cluster sized dark matter halos and found consistent results with previous estimations made with a static analysis.

In the context of scalar tensor theories for gravity, there is a universally adopted hypothesis when running N-body simulations that time derivatives in the equation of motion for the scalar field are negligible. In this work we propose to test this assumption for one specific scalar-tensor model with a gravity screening mechanism: the symmetron. To this end, we implemented the necessary modifications to include the non-static terms in the N-body code Ramses. We present test cases and results from cosmological simulations. Our main finding when comparing static vs. non-static simulations is that the global power spectrum is only slightly modified when taking into account the inclusion of non-static terms. On the contrary, we find that the calculation of the local power spectrum gives different measurements. Such results imply one must be careful when assuming the quasi-static approximation when investigating the environmental effects of modified gravity and screening mechanisms in structure formation of halos and voids distributions.

Primordial statistical anisotropy is a key indicator to investigate early Universe models and has been probed by the cosmic microwave background (CMB) anisotropies. In this paper, we examine tensor-mode CMB fluctuations generated from anisotropic gravitational waves, parametrised by $P_h({\bf k}) = P_h^{(0)}(k) [ 1 + \sum_{LM} f_L(k) g_{LM} Y_{LM} (\hat{\bf k}) ]$, where $P_h^{(0)}(k)$ is the usual scale-invariant power spectrum. Such anisotropic tensor fluctuations may arise from an inflationary model with noncommutativity of fields. It is verified that in this model, an isotropic component and a quadrupole asymmetry with $f_0(k) = f_2(k) \propto k^{-2}$ are created and hence highly red-tilted off-diagonal components arise in the CMB power spectra, namely $\ell_2 = \ell_1 \pm 2$ in $TT$, $TE$, $EE$ and $BB$, and $\ell_2 = \ell_1 \pm 1$ in $TB$ and $EB$. We find that B-mode polarisation is more sensitive to such signals than temperature and E-mode polarisation due to the smallness of large-scale cosmic variance and we can potentially measure $g_{00} = 30$ and $g_{2M} = 58$ at 68% CL in a cosmic-variance-limited experiment. Such a level of signal may be measured in a PRISM like experiment, while the instrumental noise contaminates it in the $Planck$ experiment. These results imply that it is impossible to measure the noncommutative parameter if it is small enough for the perturbative treatment to be valid. Our formalism and methodology for dealing with the CMB tensor statistical anisotropy are general and straightforwardly applicable to other early Universe models.

We study the statistics of the primordial power spectrum in models where massless gauge vectors are coupled to the inflaton, paying special attention to observational implications of having fundamental or effective horizons embedded in a bath of infrared fluctuations. As quantum infrared modes cross the horizon, they classicalize and build a background vector field. We find that the vector experiences a statistical precession phenomenon. Implications for primordial correlators and the interpretation thereof are considered. Firstly, we show how in general two, not only one, additional observables, a quadrupole amplitude and an intrinsic shape parameter, are necessary to fully describe the correction to the curvature power spectrum, and develop a unique parametrization for them. Secondly, we show that the observed anisotropic amplitude and the associated preferred direction depend on the volume of the patch being probed. We calculate non-zero priors for the expected deviations between detections based on microwave background data (which probes the entire Hubble patch) and large scale structure (which only probes a fraction of it).

Cosmological local observables are at best statistically determined by the fundamental theory describing inflation. When the scalar inflaton is coupled uniformly to a collection of subdominant massless gauge vectors, rotational invariance is obeyed locally. However, the statistical isotropy of fluctuations is spontaneously broken by gauge modes whose wavelength exceed our causal horizon. This leads to a landscape picture where primordial correlators depend on the position of the observer. We compute the stochastic corrections to the curvature power spectrum, show the existence of a new local observable (the shape of the quadrupole), and constrain the theory using Planck limits.

Aims. We investigate the gravitational redshift in clusters of galaxies in the symmetron and Hu-Sawicky $f(R)$ models. The characteristic feature of both models is the screening mechanism that hides the fifth force in dense environments recovering general relativity. Methods. We use $N$-body simulations that were run using the code \texttt{Isis}, which includes scalar fields, to analyse the deviation of observables in modified gravity models with respect to $\Lambda$CDM. Results. We find that due to the presence of the screening, the deviation is highly dependent on the halo mass. For instance, the $f(R)$ parameters $|f_{R0}|=10^{-5}$, $n=1$ cause an enhancement of the gravitational signal by up to 50% for halos with masses between $10^{13}\,M_\odot h^{-1}$ and $10^{14}\,M_\odot h^{-1}$. The characteristic mass range where the fifth force is most active varies with the model parameters. The usual assumption is that the presence of a fifth force leads to a deeper potential well and thus, a stronger gravitational redshift. However, we find that in cases in which only the central regions of the halos are screened, there could also be a weaker gravitational redshift.

We aim to constrain mass-varying neutrino models using large scale structure observations and produce forecast for the Euclid survey. We investigate two models with different scalar field potential and both positive and negative coupling parameters \beta. These parameters correspond to growing or decreasing neutrino mass, respectively. We explore couplings up to |\beta|<5. In the case of the exponential potential, we find an upper limit on $\Omega_\nu h^2$<0.004 at 2-$\sigma$ level. In the case of the inverse power law potential the null coupling can be excluded with more than 2-\sigma significance; the limits on the coupling are \beta>3 for the growing neutrino mass and \beta<-1.5 for the decreasing mass case. This is a clear sign for a preference of higher couplings. When including a prior on the present neutrino mass the upper limit on the coupling becomes |\beta|<3 at 2-$\sigma$ level for the exponential potential. Finally, we present a Fisher forecast using the tomographic weak lensing from an Euclid-like experiment and we also consider the combination with the cosmic microwave background (CMB) temperature and polarisation spectra from a Planck-like mission. If considered alone, lensing data is more efficient in constraining $\Omega_\nu$ with respect to CMB data alone. There is, however, a strong degeneracy in the \beta-$\Omega_\nu h^2$ plane. When the two data sets are combined, the latter degeneracy remains, but the errors are reduced by a factor ~2 for both parameters.

Aims. We investigated the gravitational effects of a scalar field within scalar-tensor gravity as an alternative to dark matter. Motivated by chameleon, symmetron and f(R)-gravity models, we studied a phenomenological scenario where the scalar field has both a mass (i.e. interaction length) and a coupling constant to the ordinary matter which scale with the local properties of the considered astrophysical system. Methods. We analysed the feasibility of this scenario using the modified gravitational potential obtained in its context and applied it to the galactic and hot gas/stellar dynamics in galaxy clusters and elliptical/spiral galaxies respectively. This is intended to be a first step in assessing the viability of this new approach in the context of "alternative gravity" models. Results. The main results are: 1. the velocity dispersion of elliptical galaxies can be fitted remarkably well by the suggested scalar field, with model significance similar to a classical Navarro-Frenk-White dark halo profile; 2. the analysis of the stellar dynamics and the gas equilibrium in elliptical galaxies has shown that the scalar field can couple with ordinary matter with different strengths (different coupling constants) producing and/or depending on the different clustering state of matter components; 3. elliptical and spiral galaxies, combined with clusters of galaxies, show evident correlations among theory parameters which suggest the general validity of our results at all scales and a way toward a possible unification of the theory for all types of gravitational systems we considered. All these results demonstrate that the proposed scalar field scenario can work fairly well as an alternative to dark matter.

A ghost-free theory of gravity with two dynamical metrics both coupled to matter is shown to be consistent and viable. Its cosmological implications are studied, and the models, in particular in the context of partially massless gravity, are found to explain the cosmic acceleration without resorting to dark energy.

We investigate the variation of the fine-structure constant, {\alpha}, in symmetron models using N-body simulations in which the full spatial distribution of {\alpha} at different redshifts has been calculated. In particular, we obtain simulated sky maps for this variation, and determine its power spectrum. We find that in high-density regions of space (such as deep inside dark matter halos) the value of {\alpha} approaches the value measured on Earth. In the low-density outskirts of halos the scalar field value can approach the symmetry breaking value and leads to significantly different values of {\alpha}. If the scalar-photon coupling strength {\beta}{\gamma} is of order unity we find that the variation of {\alpha} inside dark matter halos can be of the same magnitude as the recent claims by Webb et al. of a dipole variation. Importantly, our results also show that with low-redshift symmetry breaking these models exhibit some dependence of {\alpha} on lookback time (as opposed to a pure spatial dipole) which could in principle be detected by sufficiently accurate spectroscopic measurements, such as those of ALMA and the ELT-HIRES.

A screening mechanism for conformal vector-tensor modifications of general relativity is proposed. The conformal factor depends on the norm of the vector field and makes the field to vanish in high dense regions, whereas drives it to a non-null value in low density environments. Such process occurs due to a spontaneous symmetry breaking mechanism and gives rise to both the screening of fifth forces as well as Lorentz violations. The cosmology and local constraints are also computed.

Apparent similarities between non-local theories of gravity and the so-called C-theories are pointed out. It is shown that some simple C-theories can be mapped exactly into a previously considered type of ghost-free nonlocal gravity. This may introduce a useful tool to tackle some infinite-order derivative theories and raises the possibility of describing renormalisable gravity in a new context of D-theories.

The field equations in FRW background for the so called C-theories are presented and investigated. In these theories the usual Ricci scalar is substituted with $f(\mathcal{R})$ where $\mathcal{R}$ is a Ricci scalar related to a conformally scaled metric $\hat{g}_{\mu\nu} = \mathcal{C}(\mathcal{R})g_{\mu\nu}$, where the conformal factor itself depends on $\mathcal{R}$. It is shown that homogeneous perturbations of this Ricci scalar around general relativity FRW background of a large class of these theories are either inconsistent or unstable.

We study the cosmological evolution for a universe in the presence of a continuous tower of massive scalar fields which can drive the current phase of accelerated expansion of the universe and, in addition, can contribute as a dark matter component. The tower consists of a continuous set of massive scalar fields with a gaussian mass distribution. We show that, in a certain region of the parameter space, the {\it heavy} modes of the tower (those with masses much larger than the Hubble expansion rate) dominate at early times and make the tower behave like the usual single scalar field whose coherent oscillations around the minimum of the potential give a matter-like contribution. On the other hand, at late times, the {\it light} modes (those with masses much smaller than the Hubble expansion rate) overcome the energy density of the tower and they behave like a perfect fluid with equation of state ranging from 0 to -1, depending on the spectral index of the initial spectrum. This is a distinctive feature of the tower with respect to the case of quintessence fields, since a massive scalar field can only give acceleration with equation of state close to -1. Such unique property is the result of a synergy effect between the different mass modes. Interestingly, we find that, for some choices of the spectral index, the tower tracks the matter component at high redshifts (or it can even play the role of the dark matter) and eventually becomes the dominant component of the universe and give rise to an accelerated expansion.

It is shown that a disformally coupled theory in which the gravitational sector has the Einstein-Hilbert form is equivalent to a quartic DBI Galileon Lagrangian, possessing non-linear higher derivative interactions, and hence allowing for the Vainshtein effect. This Einstein Frame description considerably simplifies the dynamical equations and highlights the role of the different terms. The study of highly dense, non-relativistic environments within this description unravels the existence of a disformal screening mechanism, while the study of static vacuum configurations reveals the existence of a Vainshtein radius, at which the asymptotic solution breaks down. Disformal couplings to matter also allow the construction of Dark Energy models, which behave differently than conformally coupled ones and introduce new effects on the growth of Large Scale Structure over cosmological scales, on which the scalar force is not screened. We consider a simple Disformally Coupled Dark Matter model in detail, in which standard model particles follow geodesics of the gravitational metric and only Dark Matter is affected by the disformal scalar field. This particular model is not compatible with observations in the linearly perturbed regime. Nonetheless, disformally coupled theories offer enough freedom to construct realistic cosmological scenarios, which can be distinguished from the standard model through characteristic signatures.

In order to satisfy limits on the allowed variation of particle masses from big bang nucleosynthesis (BBN) until today, the chameleon scalar field is required to reach its attractor solution early on in its cosmological evolution. Brax et al. (2004) have shown this to be possible for certain specific initial conditions on the chameleon field at the end of inflation. However the extreme fine-tuning necessary to achieve this, poses a problem if the chameleon is to be viewed a natural candidate for dark energy. In this article we revisit the behaviour of the chameleon in the early Universe, including the additional coupling to electromagnetism proposed by Brax et al. (2011). Solving the chameleon evolution equations in the presence of a primordial magnetic field, we find that the strict initial conditions on the chameleon field at the end of inflation can be relaxed, and we determine the associated lower bound on the strength of the primordial magnetic field.

It is shown that extensions to General Relativity, which introduce a strongly coupled scalar field, can be viable if the interaction has a non-conformal form. Such disformal coupling depends upon the gradients of the scalar field. Thus, if the field is locally static and smooth, the coupling becomes invisible in the solar system: this is the disformal screening mechanism. A cosmological model is considered where the disformal coupling triggers the onset of accelerated expansion after a scaling matter era, giving a good fit to a wide range of observational data. Moreover, the interaction leaves signatures in the formation of large-scale structure that can be used to probe such couplings.

Motivated by the couplings of the dilaton in four-dimensional effective actions, we investigate the cosmological consequences of a scalar field coupled both to matter and a Maxwell-type vector field. The vector field has a background isotropy-violating component. New anisotropic scaling solutions which can be responsible for the matter and dark energy dominated epochs are identified and explored. For a large parameter region the universe expands almost isotropically. Using that the CMB quadrupole is extremely sensitive to shear, we constrain the ratio of the matter coupling to the vector coupling to be less than 10^(-5). Moreover, we identify a large parameter region, corresponding to a strong vector coupling regime, yielding exciting and viable cosmologies close to the LCDM limit.

This paper has been withdrawn by the author due to a mistake.

The scalar-tensor Dirac-Born-Infeld (DBI) inflation scenario provides a simple mechanism to reduce the large values of the boost factor associated with single field models with DBI action, whilst still being able to drive 60 efolds of inflation. Using a slow-roll approach, we obtain an analytical expression for the spectral index of the perturbations and, moreover, determine numerically the regions of the parameter space of the model capable of giving rise to a power spectrum with amplitude and spectral index within the observed bounds. We find that regions that exhibit significant DBI effects throughout the inflationary period can be discarded by virtue of a blue-tilted spectral index, however, there are a number of viable cases --- associated with a more red-tilted spectral index --- for which the boost factor is initially suppressed by the effect of the coupling between the fields, but increases later to moderate values.

Recently an inflationary model with a vector field coupled to the inflaton was proposed and the phenomenology studied for the Bianchi type I spacetime. It was found that the model demonstrates a counter-example to the cosmic no-hair theorem since there exists a stable anisotropically inflationary fix-point. One of the great triumphs of inflation, however, is that it explains the observed flatness and isotropy of the universe today without requiring special initial conditions. Any acceptable model for inflation should thus explain these observations in a satisfactory way. To check whether the model meets this requirement, we introduce curvature to the background geometry and consider axisymmetric spacetimes of Bianchi type II,III and the Kantowski-Sachs metric. We show that the anisotropic Bianchi type I fix-point is an attractor for the entire family of such spacetimes. The model is predictive in the sense that the universe gets close to this fix-point after a few e-folds for a wide range of initial conditions. If inflation lasts for N e-folds, the curvature at the end of inflation is typically of order exp(-2N). The anisotropy in the expansion rate at the end of inflation, on the other hand, while being small on the one-percent level, is highly significant. We show that after the end of inflation there will be a period of isotropization lasting for about 2N/3 e-folds. After that the shear scales as the curvature and becomes dominant around N e-folds after the end of inflation. For plausible bounds on the reheat temperature the minimum number of e-folds during inflation, required for consistency with the isotropy of the supernova Ia data, lays in the interval (21,48). Thus the results obtained for our restricted class of spacetimes indicates that inflation with anisotropic hair is cosmologically viable.

We investigate the environment dependence of dark matter halos in the symmetron modified gravity scenario. The symmetron is one of three known mechanisms for screening a fifth-force and thereby recovering General Relativity in dense environments. The effectiveness of the screening depends on both the mass of the object and the environment it lies in. Using high-resolution N-body simulations we find a significant difference, which depends on the halos mass and environment, between the lensing and dynamical masses of dark matter halos similar to the f(R) modified gravity. The symmetron can however yield stronger signatures due to a freedom in the strength of the coupling to matter.

We provide an exhaustive analysis of the Integrated Sach-Wolfe effect (ISW) in the context of coupled Dark Energy cosmologies where a component of massive neutrinos is also present. We focus on the effects of both the coupling between Dark Matter and Dark Energy and of the neutrino mass on the cross-correlation between galaxy/quasar distributions and ISW effect. We provide a simple expression to appropriately rescale the galaxy bias when comparing different cosmologies. Theoretical predictions of the cross-correlation function are then compared with observational data. We find that, while it is not possible to distinguish among the models at low redshifts, discrepancies between coupled models and $\Lambda$CDM increase with $z$. In spite of this, current data alone seems not able to distinguish between coupled models and $\Lambda$CDM. However, we show that upcoming galaxy surveys will permit tomographic analysis which allow to better discriminate among the models. We discuss the effects on cross-correlation measurements of ignoring galaxy bias evolution, b(z), and magnification bias correction and provide fitting formulae for b(z) for the cosmologies considered. We compare three different tomographic schemes and investigate how the expected signal to noise ratio, snr, of the ISW-LSS cross-correlation changes when increasing the number of tomographic bins. The dependence of snr on the area of the survey and the survey shot noise is also discussed.

We test anisotropic dark energy models with the 7-year WMAP temperature observations data. In the presence of imperfect sources, due to large-scale gradients or anisotropies in the dark energy field, the CMB sky will be distorted anisotropically on its way to us by the ISW effect. The signal covariance matrix then becomes nondiagonal for small multipoles, but at $\ell \gtrsim 20$ the anisotropy is negligible. We parametrize possible violations of rotational invariance in the late universe by the magnitude of a post-Friedmannian deviation from isotropy and its scale dependence. This allows to obtain hints on possible imperfect nature of dark energy and the large-angle anomalous features in the CMB. A robust statistical analysis, subjected to various tests and consistency checks, is performed to compare the predicted correlations with those obtained from the satellite-measured CMB full sky maps. The preferred axis point towards $(l,b) = (168^\circ, -31^\circ)$ and the amplitude of the anisotropy is $\varpi_0 = (0.51\pm 0.94)$ (1$\sigma$ deviation quoted). The best-fit model has a steep blue anisotropic spectrum ($n_{\mathrm{de}} = 3.1\pm1.5$).

Scalar fields, strongly coupled to matter, can be present in nature and still be invisible to local experiments if they are subject to a screening mechanism. The symmetron is one such mechanism which relies on restoration of a spontaneously broken symmetry in regions of high density to shield the scalar fifth force. We have investigated structure formation in the symmetron model by using N-body simulations and find strong observable signatures in both the linear and nonlinear matter power spectrum and on the halo mass function. The mechanism for suppressing the scalar fifth force in high density regions is also found to work very well.

We investigate whether there are any cosmological evidences for a scalar field with a mass and coupling to matter which change accordingly to the properties of the astrophysical system it "lives in", without directly focusing on the underlying mechanism that drives the scalar field scale-dependent properties. We assume a Yukawa type of coupling between the field and matter and also that the scalar field mass grows with density, in order to overcome all gravity constraints within the solar system. We analyse three different gravitational systems assumed as "cosmological indicators": supernovae type Ia, low surface brightness spiral galaxies and clusters of galaxies. Results show that: a) a quite good fit to the rotation curves of low surface brightness galaxies only using visible stellar and gas mass components is obtained; b) a scalar field can fairly well reproduce the matter profile in clusters of galaxies, estimated by X-ray observations and without the need of any additional dark matter; c) there is an intrinsic difficulty in extracting information about the possibility of a scale-dependent massive scalar field (or more generally about a varying gravitational constant) from supernovae type Ia.

The Dirac-Born-Infeld (DBI) action has been widely studied as an interesting example of a model of k-inflation in which the sound speed of the cosmological perturbations differs from unity. In this article we consider a scalar-tensor theory in which the matter component is a field with a DBI action. Transforming to the Einstein frame, we explore the effect of the resulting coupling on the background dynamics of the fields and the first-order perturbations. We find that the coupling forces the scalar field into the minimum of its effective potential. While the additional scalar field contributes significantly to the energy density during inflation, the dynamics are determined by the DBI field, which has the interesting effect of increasing the number of efolds of inflation and decreasing the boost factor of the DBI field. Focusing on this case, we show, with the benefit of numerical examples, that the power spectrum of the primordial perturbations is determined by the behaviour of the perturbations of the modified DBI field.

In chameleon field theories a scalar field can couple to matter with gravitational strength and still evade local gravity constraints due to a combination of self-interactions and the couplings to matter. Originally, these theories were proposed with a constant coupling to matter, however, the chameleon mechanism also extends to the case where the coupling becomes field-dependent. We study the cosmology of chameleon models with power-law couplings and power-law potentials. It is found that these generalized chameleons, when viable, have a background expansion very close to LCDM, but can in some special cases enhance the growth of the linear perturbations at low redshifts. For the models we consider it is found that this region of the parameter space is ruled out by local gravity constraints. Imposing a coupling to dark matter only, the local constraints are avoided, and it is possible to have observable signatures on the linear matter perturbations.

It has been recently proposed that the interpretation of gravity as an emergent, entropic force might have nontrivial implications to cosmology. Here two approaches are investigated: in one, the Friedman equation receives entropic contributions from the usually neglected surface terms, and in another, the extra terms are derived from quantum corrections to the entropy formula. UV terms may drive inflation, avoiding a recently derived no-go theorem, though in some cases leading to a graceful exit problem. IR terms can generate dark energy, alleviating the cosmological constant problem. The quantum corrections are bounded by their implications to the BBN, and the surface terms are constrained in addition by their effect upon the behavior of matter. Likelihood analyses are performed to constrain the modifications by the SNeIa, BAO and CMB data. It is found that a monomial correction to the area-entropy formula results in late acceleration in very good agreement with observations, which then turn out to be compatible with positive curvature. The evolution of perturbations is deduced by assuming the Jebsen-Birkhoff theorem. Distinct signatures can then be identified in the large scale structure formation. Furthermore, it is shown that the visible universe satisfies the Bekenstein bound.

Non-commutative geometry at inflation can give arise to parity violating modulations of the primordial power spectrum. We develop the statistical tools needed for investigating whether these modulations are evident in the Cosmic Microwave Background (CMB). The free parameters of the models are two directional parameters (theta,phi), the signal amplitude A*, and a tilt parameter n* that modulates correlation power on different scales. The signature of the model corresponds to a kind of hemispherical power asymmetry. When analyzing the 7-year WMAP data we find a weak signature for a preferred direction in the Q-, V-, and W bands with direction (l,b) = (-225 deg,-25 deg) +- (20 deg, 20 deg), which is close to another previously discovered hemispherical power asymmetry. Although these results are intriguing, the significance of the detection in the W-, V- and Q-bands are nonzero at about 2 sigma, suggesting that the simplest parameterization of the leading correction represents only partially the effects of the space-time non-commutativity possibly responsible for the hemispherical asymmetry. Our constraints on the presence of a dipole are independent of its physical origin and prefer a blue-tilted spectral index n* ~ 0 with the amplitude A* ~ 0.18.

In addition to shear and vorticity a homogeneous background may also exhibit anisotropic curvature. Here a class of spacetimes is shown to exist where the anisotropy is solely of the latter type, and the shear-free condition is supported by a canonical, massless 2-form field. Such spacetimes possess a preferred direction in the sky and at the same time a CMB which is isotropic at the background level. A distortion of the luminosity distances is derived and used to test the model against the CMB and supernovae (using the Union catalog), and it is concluded that the latter exhibit a higher-than-expected dependence on angular position. It is shown that future surveys could detect a possible preferred direction by observing ~ 20 / (\Omega_{k0}^2) supernovae over the whole sky.

Noncommutative geometry can provide effective description of physics at very short distances taking into account generic effects of quantum gravity. Inflation amplifies tiny quantum fluctuations in the early universe to macroscopic scales and may thus imprint high energy physics signatures in the cosmological perturbations that could be detected in the CMB. It is shown here that this can give rise to parity-violating modulations of the primordial spectrum and odd non-Gaussian signatures. The breaking of rotational invariance of the CMB provides constraints on the scale of noncommutativity that are competitive with the existing noncosmological bounds, and could explain the curious hemispherical asymmetry that has been claimed to be observed in the sky. This introduces also non-Gaussianity with peculiar shape- and scale-dependence, which in principle allows an independent cross-check of the presence of noncommutativity at inflation.

In chameleon dark energy models, local gravity constraints tend to rule out parameters in which observable cosmological signatures can be found. We study viable chameleon potentials consistent with a number of recent observational and experimental bounds. A novel chameleon field potential, motivated by f(R) gravity, is constructed where observable cosmological signatures are present both at the background evolution and in the growth-rate of the perturbations. We study the evolution of matter density perturbations on low redshifts for this potential and show that the growth index today gamma_0 can have significant dispersion on scales relevant for large scale structures. The values of gamma_0 can be even smaller than 0.2 with large variations of gamma on very low redshifts for the model parameters constrained by local gravity tests. This gives a possibility to clearly distinguish these chameleon models from the Lambda-Cold-Dark-Matter model in future high-precision observations.

In this paper we start from the original formulation of the galileon model with the original choice for couplings to gravity. Within this framework we find that there is still a subset of possible Lagrangians that give selfaccelerating solutions with stable spherically symmetric solutions. This is a certain constrained subset of the third order galileon which has not been explored before. We develop and explore the background cosmological evolution of this model drawing intuition from other even more restricted galileon models. The numerical results confirm the presence of selfacceleration, but also reveals a possible instability with respect to galileon perturbations.

We have studied the Bekenstein-Sandvik-Barrow-Magueijo (BSBM) model for the spatial and temporal variations of the fine structure constant, alpha, with the aid of full N-body simulations which explicitly and self-consistently solve for the scalar field driving the alpha-evolution. We focus on the scalar field (or equivalently alpha) inside the dark matter halos and find that the profile of the scalar field is essentially independent of the BSBM model parameter. This means that given the density profile of an isolated halo and the background value of the scalar field, we can accurately determine the scalar field perturbation in that halo. We also derive an analytic expression for the scalar-field perturbation using the Navarro-Frenk-White halo profile, and show that it agrees well with numerical results, at least for isolated halos; for non-isolated halos this prediction differs from numerical result by a (nearly) constant offset which depends on the environment of the halo.

We introduce the N-body simulation technique to follow structure formation in linear and nonlinear regimes for the extended quintessence models (scalar-tensor theories in which the scalar field has a self-interaction potential and behaves as dark energy), and apply it to a class of models specified by an inverse power-law potential and a non-minimal coupling. Our full solution of the scalar field perturbation confirms that, when the potential is not too nonlinear, the effects of the scalar field could be accurately approximated as a modification of background expansion rate plus a rescaling of the effective gravitational constant relevant for structure growth. For the models we consider, these have opposite effects, leading to a weak net effect in the linear perturbation regime. However, on the nonlinear scales the modified expansion rate dominates and could produce interesting signatures in the matter power spectrum and mass function, which might be used to improve the constraints on the models from cosmological data. We show that the density profiles of the dark matter halos are well described by the Navarro-Frenk-White formula, although the scalar field could change the concentration. We also derive an analytic formula for the scalar field perturbation inside halos assuming NFW density profile and sphericity, which agrees well with numerical results if the parameter is appropriately tuned. The results suggest that for the models considered, the spatial variation of the scalar field (and thus the locally measured gravitational constant) is very weak, and so local experiments could see the background variation of gravitational constant.

We introduce and carefully define an entire class of field theories based on non-standard spinors. Their dominant interaction is via the gravitational field which makes them naturally dark; we refer to them as Dark Spinors. We provide a critical analysis of previous proposals for dark spinors noting that they violate Lorentz invariance. As a working assumption we restrict our analysis to non-standard spinors which preserve Lorentz invariance, whilst being non-local and explicitly construct such a theory. We construct the complete energy-momentum tensor and derive its components explicitly by assuming a specific projection operator. It is natural to next consider dark spinors in a cosmological setting. We find various interesting solutions where the spinor field leads to slow roll and fast roll de Sitter solutions. We also analyse models where the spinor is coupled conformally to gravity, and consider the perturbations and stability of the spinor.

Certain scalar-tensor theories exhibit the so-called chameleon mechanism, whereby observational signatures of scalar fields are hidden by a combination of self-interactions and interactions with ambient matter. Not all scalar-tensor theories exhibit such a chameleon mechanism, which has been originally found in models with inverse power run-away potentials and field independent couplings to matter. In this paper we investigate field-theories with field-dependent couplings and a power-law potential for the scalar field. We show that the theory indeed is a chameleon field theory. We find the thin-shell solution for a spherical body and investigate the consequences for E\"ot-Wash experiments, fifth-force searches and Casimir force experiments. Requiring that the scalar-field evades gravitational tests, we find that the coupling is sensitive to a mass-scale which is of order of the Hubble scale today.

Disformal transformations have proven to be very useful to devise models of the dark sector. In the present paper we apply such transformation to a single scalar field theory as a way to drive the field into a slow roll phase. The canonical scalar field Lagrangian, when coupled to a disformal metric, turns out to have relations to bimetric dark matter theories and to describe many specific dark energy models at various limits, thus providing a surprisingly simple parametrisation of a wide variety of models including tachyon, Chaplygin gas, K-essence and dilatonic ghost condensate. We investigate the evolution of the background and linear perturbations in disformal quintessence in order to perform a full comparison of the predictions with the cosmological data. The dynamics of the expansion, in particular the mechanism of the transition to accelerating phase, is described in detail. We then study the effects of disformal quintessence on cosmic microwave background (CMB) anisotropies and large scale structures (LSS). A likelihood analysis using the latest data on wide-ranging SNIa, CMB and LSS observations is performed allowing variations in six cosmological parameters and the two parameters specifying the model. We find that while a large region of parameter space remains compatible with observations, models featuring either too much early dark energy or too slow transition to acceleration are ruled out.

The apparent alignment of the cosmic microwave background multipoles on large scales challenges the standard cosmological model. Scalar field inflation is isotropic and cannot account for the observed alignment. We explore the imprints, a non-standard spinor driven inflation would leave on the cosmic microwave background anisotropies. We show it is natural to expect an anisotropic inflationary expansion of the Universe which has the effect of suppressing the low multipole amplitude of the primordial power spectrum, while at the same time to provide the usual inflationary features.

A quintessence scalar field or cosmon interacting with neutrinos can have important effects on cosmological structure formation. Within growing neutrino models the coupling becomes effective only in recent times, when neutrinos become non-relativistic, stopping the evolution of the cosmon. This can explain why dark energy dominates the universe only in a rather recent epoch by relating the present dark energy density to the small mass of neutrinos. Such models predict the presence of stable neutrino lumps at supercluster scales (~200 Mpc and bigger), caused by an attractive force between neutrinos which is stronger than gravity and mediated by the cosmon. We present a method to follow the initial non-linear formation of neutrino lumps in physical space, by integrating numerically on a 3D grid non-linear evolution equations, until virialization naturally occurs. As a first application, we show results for cosmologies with final large neutrino average mass ~2 eV: in this case, neutrino lumps indeed form and mimic very large cold dark matter structures, with a typical gravitational potential 10^{-5} for a lump size ~10 Mpc, and reaching larger values for lumps of about 200 Mpc. A rough estimate of the cosmological gravitational potential at small k in the non-linear regime, Phi_nu = 10^{-6} (k/k_0)^{-2}, 1.2x10^{-2} h/Mpc < k_0 < 7.8x10^{-2} h/Mpc, turns out to be many orders of magnitude smaller than an extrapolation of the linear evolution of density fluctuations. The size of the neutrino-induced gravitational potential could modify the spectrum of CMB anisotropies for small angular momenta.

We study the evolution of cosmological perturbations in f(G) gravity, where the Lagrangian is the sum of a Ricci scalar R and an arbitrary function f in terms of a Gauss-Bonnet term G. We derive the equations for perturbations assuming matter to be described by a perfect fluid with a constant equation of state w. We show that density perturbations in perfect fluids exhibit negative instabilities during both the radiation and the matter domination, irrespective of the form of f(G). This growth of perturbations gets stronger on smaller scales, which is difficult to be compatible with the observed galaxy spectrum unless the deviation from General Relativity is very small. Thus f(G) cosmological models are effectively ruled out from this Ultra-Violet instability, even though they can be compatible with the late-time cosmic acceleration and local gravity constraints.

Anisotropic and isotropic slow-roll inflation supported by n-forms are sought. Canonical field strengths and their duals are taken into account, and they are allowed to have a potential and also, when necessary for slow-roll, a nonminimal curvature coupling. New isotropic solutions are found for three-forms. It is also shown that some n-form actions are equivalent to f(R) gravity and scalar field models with possible nonminimal couplings. Anisotropic solutions are found for two-forms, generalising vector inflation, however, also the two-form is unstable during inflation due to the nonminimal coupling to curvature. We also address the stability of the isotropic background solutions supported by a triad of vectors .

We investigate dark matter halo properties as a function of a time--varying dark energy equation of state. The dynamics of the collapse of the halo is governed by the form of the quintessence potential, the time evolution of its equation of state, the initial conditions of the field and its homogeneity nature in the highly non--linear regime. These have a direct impact on the turnaround, virialisation and collapse times, altering in consequence the non--linear density contrast and virial radius. We compute halo concentrations using the Eke, Navarro & Steinmetz algorithm, examining two extreme scenarios: first, we assume that the quintessence field does not exhibit fluctuations on cluster scales and below - homogeneous fluid; second, we assume that the field inside the overdensity collapses along with the dark matter - inhomogeneous fluid. The Eke, Navarro & Steinmetz prescription reveals, in general, higher halo concentrations in inhomogeneous dark energy models than in their homogeneous equivalents. Halo concentrations appear to be controlled by both changes in formation epochs of the halo cores as well as by differing virialisation overdensities. We derive physical halo properties in all models and discuss their observational implications. We examine two possible methods for comparing observations with theoretical predictions. The first method works on galaxy cluster scales and consists of fitting the observed X-ray cluster gas density distributions to those predicted for an NFW profile. The second method works on galaxy scales and involves the observational measurement of the so--called central density parameter.

We show that, within some modified gravity theories, such as the Palatini models, the non-linear nature of the field equations implies that the usual naive averaging procedure (replacing the microscopic energy-momentum by its cosmological average) could be invalid. As a consequence, the relative motion of particles in Palatini theories is actually indistinguishable from that predicted by General Relativity. Moreover, there is no WEP violation. Our new and most important result is that the cosmology and astrophysics, or put more generally, the behaviours on macroscopic scales, predicted by these two theories are the same, and as a result the naturalness problems associated with the cosmological constant are not alleviated. Palatini gravity does however predict alterations to the internal structure of particles and the particle physics laws, e.g., corrections to the hydrogen energy levels. Measurements of which place strong constraints on the properties of viable Palatini gravities.

We show that growing neutrino models, in which the growing neutrino mass stops the dynamical evolution of a dark energy scalar field, lead to a substantial neutrino clustering on the scales of superclusters. Nonlinear neutrino lumps form at redshift z ~ 1 and could partially drag the clustering of dark matter. If observed, large scale non-linear structures could be an indication for a new attractive force stronger than gravity and mediated by the 'cosmon' dark energy scalar field.

We consider several new classes of viable vector field alternatives to the inflaton and quintessence scalar fields. Spatial vector fields are shown to be compatible with the cosmological anisotropy bounds if only slightly displaced from the potential minimum while dominant, or if driving an anisotropic expansion with nearly vanishing quadropole today. The Bianchi I model with a spatial field and an isotropic fluid is studied as a dynamical system, and several types of scaling solutions are found. On the other hand, time-like fields are automatically compatible with large-scale isotropy. We show that they can be dynamically important if non-minimal gravity couplings are taken into account. As an example, we reconstruct a vector-Gauss-Bonnet model which generates the concordance model acceleration at late times and supports an inflationary epoch at high curvatures. The evolution of vortical perturbations is considered.

We show that, within modified gravity, the non-linear nature of the field equations implies that the usual naive averaging procedure (replacing the microscopic energy-momentum by its cosmological average) is invalid. We discuss then how the averaging should be performed correctly and show that, as a consequence, at classical level the physical masses and geodesics of particles, cosmology and astrophysics in Palatini modified gravity theories are all indistinguishable from the results of general relativity plus a cosmological constant. Palatini gravity is however a different theory from general relativity and predicts different internal structures of particles from the latter. On the other hand, and in contrast to classical particles, the electromagnetic field permeates in the space, hence a different averaging procedure should be applied here. We show that in general Palatini gravity theories would then affect the propagation of photons, thus changing the behaviour of a Universe dominated by radiation. Finally, Palatini theories also predict alterations to particle physics laws. For example, it can lead to sensitive corrections to the hydrogen energy levels, the measurements of which could be used to place very strong constraints on the properties of viable Palatini gravity theories

We propose a method to study and constrain modified gravity theories for dark matter using CMB temperature anisotropies and polarization. We assume that the theories considered here have already passed the matter power-spectrum test of large-scale structure. With this requirement met, we show that a modified gravity theory can be specified by parametrizing the time evolution of its dark-matter density contrast, which is completely controlled by the dark matter stress history. We calculate how the stress history with a given parametrization affects the CMB observables, and a qualitative discussion of the physical effects involved is supplemented with numerical examples. It is found that, in general, alternative gravity theories can be efficiently constrained by the CMB temperature and polarization spectra. There exist, however, special cases where modified gravity cannot be distinguished from the CDM model even by using both CMB and matter power spectrum observations, nor can they be efficiently restricted by other observables in perturbed cosmologies. Our results show how the stress properties of dark matter, which determine the evolutions of both density perturbations and the gravitational potential, can be effectively investigated using just the general conservation equations and without assuming any specific theoretical gravitational theory within a wide class.

We investigate cosmologies where the accelerated expansion of the Universe is driven by a field with an anisotropic equation of state. We model such scenarios within the Bianchi I framework, introducing two skewness parameters to quantify the deviation of pressure from isotropy. We study the dynamics of the background expansion in these models. A special case of anisotropic cosmological constant is analyzed in detail. The anisotropic expansion is then confronted with the redshift and angular distribution of the supernovae type Ia. In addition, we investigate the effects on the cosmic microwave background (CMB) anisotropies for which the main signature appears to be a quadrupole contribution. We find that the two skewness parameters can be very well constrained. Tightest bounds are imposed by the CMB quadrupole, but there are anisotropic models which avoid this bound completely. Within these bounds, the anisotropy can be beneficial as a potential explanation of various anomalous cosmological observations, especially in the CMB at the largest angles. We also consider the dynamics of linear perturbations in these models. The covariant approach is used to derive the general evolution equations for cosmological perturbations taking into account imperfect sources in an anisotropic background. The implications for the galaxy formation are then studied. These results might help to make contact between the observed anomalies in CMB and large scale structure and fundamental theories exhibiting Lorentz violation.

As a result of non-linear self-interactions, in chameleon theories where the field couples to matter much more strongly than gravity does, the fifth force between two bodies with thin-shell is independent of their coupling to the field. As a consequence the bounds on the coupling coming from terrestrial tests of gravity, measurements of the Casimir force and those constraints imposed by the physics of compact objects, big-bang nucleosynthesis and measurements of the cosmic microwave background anisotropies can be exponentially relaxed.

We study the cosmological implications of including angular motion in the DBI brane inflation scenario. The non-canonical kinetic terms of the Dirac-Born-Infeld action give an interesting alternative to slow roll inflation, and cycling branes can drive periods of accelerated expansion in the Universe. We present explicit numerical solutions demonstrating brane inflation in the Klebanov-Strassler throat. We find that demanding sufficient inflation takes place in the throat is in conflict with keeping the brane's total energy low enough so that local gravitational backreaction on the Calabi-Yau manifold can be safely ignored. We deduce that spinflation (brane inflation with angular momentum) can ease this tension by providing extra e-foldings at the start of inflation. Cosmological expansion rapidly damps the angular momentum causing an exit to a more conventional brane inflation scenario. Finally, we set up a general framework for cosmological perturbation theory in this scenario, where we have multi-field non-standard kinetic term inflation.

A growing neutrino mass can stop the dynamical evolution of a dark energy scalar field, thus explaining the 'why now' problem. We show that such models lead to a substantial neutrino clustering on the scales of superclusters. Nonlinear neutrino lumps form at redshift z \sim 1 and could partially drag the clustering of dark matter. If observed, large scale non-linear structures could be an indication for a new attractive 'cosmon force' stronger than gravity.

A coupling between a light scalar field and neutrinos has been widely discussed as a mechanism for linking (time varying) neutrino masses and the present energy density and equation of state of dark energy. However, it has been pointed out that the viability of this scenario in the non-relativistic neutrino regime is threatened by the strong growth of hydrodynamic perturbations associated with a negative adiabatic sound speed squared. In this paper we revisit the stability issue in the framework of linear perturbation theory in a model independent way. The criterion for the stability of a model is translated into a constraint on the scalar-neutrino coupling, which depends on the ratio of the energy densities in neutrinos and cold dark matter. We illustrate our results by providing meaningful examples both for stable and unstable models.

If the dark energy equation of state is anisotropic, the expansion rate of the universe becomes direction-dependent at late times. We show that such models are not only cosmologically viable but that they could explain some of the observed anomalies in the CMB, and shed some light into the coincidence problem. The possible anisotropy can then be constrained by studying its effects on the luminosity distance-redshift relation inferred from several observations. A vector field action for dark energy is also presented as an example of such possibility.

We investigate the clustering of dark energy within matter overdensities and voids. In particular, we derive an analytical expression for the dark energy density perturbations, which is valid both in the linear, quasi-linear and fully non-linear regime of structure formation. We also investigate the possibility of detecting such dark energy clustering through the ISW effect. In the case of uncoupled quintessence models, if the mass of the field is of order the Hubble scale today or smaller, dark energy fluctuations are always small compared to the matter density contrast. Even when the matter perturbations enter the non-linear regime, the dark energy perturbations remain linear. We find that virialised clusters and voids correspond to local overdensities in dark energy, with $\delta_{\phi}/(1+w) \sim \Oo(10^{-5})$ for voids, $\delta_{\phi}/(1+w) \sim \Oo(10^{-4})$ for super-voids and $\delta_{\phi}/(1+w) \sim \Oo(10^{-5})$ for a typical virialised cluster. If voids with radii of $100-300 {\rm Mpc}$ exist within the visible Universe then $\delta_{\phi}$ may be as large as $10^{-3}(1+w)$. Linear overdensities of matter and super-clusters generally correspond to local voids in dark energy; for a typical super-cluster: $\delta_{\phi}/(1+w) \sim \Oo(-10^{-5})$. The approach taken in this work could be straightforwardly extended to study the clustering of more general dark energy models.

We derive a semi-analytical extension of the spherical collapse model of structure formation that takes account of the effects of deviations from spherical symmetry and shell crossing which are important in the non-linear regime. Our model is designed so that it predicts a relation between the peculiar velocity and density contrast that agrees with the results of N-body simulations in the region where such a comparison can sensibly be made. Prior to turnaround, when the unmodified spherical collapse model is expect to be a good approximation, the predictions of the two models coincide almost exactly. The effects of a late time dominating dark energy component are also taken into account. The improved spherical collapse model is a useful tool when one requires a good approximation not just to the evolution of the density contrast but also its trajectory. Moreover, the analytical fitting formulae presented is simple enough to be used anywhere where the standard spherical collapse might be used but with the advantage that it includes a realistic model of the effects of virialisation.

We consider cosmology in the Einstein-aether theory (the generally covariant theory of gravitation coupled to a dynamical timelike Lorentz-violating vector field) with a linear aether-Lagrangian. The 3+1 spacetime splitting approach is used to derive covariant and gauge invariant perturbation equations which are valid for a general class of Lagrangians. Restricting attention to the parameter space of these theories which is consistent with local gravity experiments, we show that there are tracking behaviors for the aether field, both in the background cosmology and at linear perturbation level. The primordial power-spectrum of scalar perturbations in this model is shown to be the same that predicted by standard general relativity. However, the power-spectrum of tensor perturbation is different from that in general relativity, but has a smaller amplitude and so cannot be detected at present. We also study the implications for late-time cosmology and find that the evolution of photon and neutrino anisotropic stresses can source the aether field perturbation during the radiation and matter dominated epochs, and as a result the CMB and matter power spectra are modified. However these effects are degenerate with respect to other cosmological parameters, such as neutrino masses and the bias parameter in the observed galaxy spectrum.

We explore the effects of a positive cosmological constant on astrophysical and cosmological configurations described by a polytropic equation of state. We derive the conditions for equilibrium and stability of such configurations and consider some astrophysical examples where our analysis may be relevant. We show that in the presence of the cosmological constant the isothermal sphere is not a viable astrophysical model since the density in this model does not go asymptotically to zero. The cosmological constant implies that, for polytropic index smaller than five, the central density has to exceed a certain minimal value in terms of the vacuum density in order to guarantee the existence of a finite size object. We examine such configurations together with effects of $\Lambda$ in other exotic possibilities, such as neutrino and boson stars, and we compare our results to N-body simulations. The astrophysical properties and configurations found in this article are specific features resulting from the existence of a dark energy component. Hence, if found in nature would be an independent probe of a cosmological constant, complementary to other observations.

We consider the cosmology of the Ricci-tensor-squared gravity in the Palatini variational approach. The gravitational action of standard general relativity is modified by adding a function f(R^abR_ab) to the Einstein-Hilbert action, and the Palatini variation is used to derive the field equations. A general method of obtaining the background and first-order covariant and gauge-invariant perturbation equations is outlined. As an example, we consider the cosmological constraints on such theories arising from the supernova type Ia and cosmic microwave background observations. We find that the best fit to the data is a non-null leading-order correction to Einstein gravity, but the current data exhibit no significant preference over the concordance model. The growth of non-relativistic matter density perturbations at late times is also analyzed, and we find that a scale-dependent (positive or negative) sound-speed-squared term generally appears in the growth equation for small-scale density perturbations. We also estimate the observational bound imposed by the matter power spectrum for the model with f(R^abR_ab) = alpha(R^abR_ab)^beta to be roughly |\beta| \lesssim O(10^{-5}) so long as the dark matter does not possess compensating anisotropic stresses.

We investigate the possibility of using cosmological observations to probe and constrain an imperfect dark energy fluid. We consider a general parameterization of the dark energy component accounting for an equation of state, speed of sound and viscosity. We use present and future data from the cosmic microwave background radiation (CMB), large scale structures and supernovae type Ia. We find that both the speed of sound and viscosity parameters are difficult to nail down with the present cosmological data. Also, we argue that it will be hard to improve the constraints significantly with future CMB data sets. The implication is that a perfect fluid description might ultimately turn out to be a phenomenologically sufficient description of all the observational consequences of dark energy. The fundamental lesson is however that even then one cannot exclude, by appealing to observational evidence alone, the possibility of imperfectness in dark energy.

The effects of mass-varying neutrinos on cosmic microwave background (CMB) anisotropies and large scale structures (LSS) are studied. In these models, dark energy and neutrinos are coupled such that the neutrino masses are functions of the scalar field playing the role of dark energy. We begin by describing the cosmological background evolution of such a system. It is pointed out that, similar to models with a dark matter/dark energy interaction, the apparent equation of state measured with SNIa can be smaller than -1. We then discuss the effect of mass-varying neutrinos on the CMB anisotropies and the matter power spectrum. A suppression of power in the CMB power spectrum at large angular scales is usually observed. We give an explanation for this behaviour and discuss different couplings and quintessence potentials to show the generality of the results obtained. We perform a likelihood analysis using wide-ranging SNIa, CMB and LSS observations to assess whether such theories are viable. Treating the neutrino mass as a free parameter we find that the constraints on the coupling are weak, since CMB and LSS surveys give only upper bounds on the neutrino mass. However, fixing a priori the neutrino masses, we find that there is some evidence that the existence of such a coupling is actually preferred by current cosmological data over the standard LambdaCDM cosmology.

It was recently argued that the observed PVLAS anomaly can be explained by chameleon field theories in which large deviations from Newton's law can be avoided. Here we present the predictions for the dichroism and the birefringence induced in the vacuum by a magnetic field in these models. We show that chameleon particles behave very differently from standard axion-like particles (ALPs). We find that, unlike ALPs, the chameleon particles are confined within the experimental set-up. As a consequence, the birefringence is always bigger than the dichroism in PVLAS-type experiments.

We consider the cosmology where some function f(G) of the Gauss-Bonnet term G is added to the gravitational action to account for the late-time accelerating expansion of the universe. The covariant and gauge invariant perturbation equations are derived with a method which could also be applied to general f(R,R^abR_ab,R^abcdR_abcd) gravitational theories. It is pointed out that, despite their fourth-order character, such f(G) gravity models generally cannot reproduce arbitrary background cosmic evolutions; for example, the standard LCDM paradigm with Omega_DE = 0.76 cannot be realized in f(G) gravity theories unless f is a true cosmological constant because it imposes exclusionary constraints on the form of f(G). We analyze the perturbation equations and find that, as in f(R) model, the stability of early-time perturbation growth puts some constraints on the functional form of f(G), in this case d^2 f/d G^2 < 0. Furthermore, the stability of small-scale perturbations also requires that f not deviate significantly from a constant. These analyses are illustrated by numerically propagating the perturbation equations with a specific model reproducing a representative LCDM cosmic history. Our results show how the f(G) models are highly constrained by cosmological data.

We study the background cosmology governed by the Tensor-Vector-Scalar theory of gravity proposed by Bekenstein. We consider a broad family of potentials that lead to modified gravity and calculate the evolution of the field variables both numerically and analytically. We find a range of possible behaviors, from scaling to the late time domination of either the additional gravitational degrees of freedom or the background fluid.

We show that, as a result of non-linear self-interactions, it is feasible, at least in light of the bounds coming from terrestrial tests of gravity, measurements of the Casimir force and those constraints imposed by the physics of compact objects, big-bang nucleosynthesis and measurements of the cosmic microwave background, for there to exist, in our Universe, one or more scalar fields that couple to matter much more strongly than gravity does. These scalar fields behave like chameleons: changing their properties to fit their surroundings. As a result these scalar fields can be not only very strongly coupled to matter, but also remain relatively light over solar system scales. These fields could also be detected by a number of future experiments provided they are properly designed to do so. These results open up an altogether new window, which might lead to a completely different view of the role played by light scalar fields in particle physics and cosmology.

Cosmological consequences of a string-motivated dark energy scenario featuring a scalar field coupled to the Gauss-Bonnet invariant are investigated. We study the evolution of the universe in such a model, identifying its key properties. The evolution of the homogeneous background and cosmological perturbations, both at large and small scales, are calculated. The impact of the coupling on galaxy distributions and the cosmic microwave background is examined. We find the coupling provides a mechanism to viably onset the late acceleration, to alleviate the coincidence problem, and furthermore to effectively cross the phantom divide at the present while avoiding a Big Rip in the future. We show the model could explain the present cosmological observations, and discuss how various astrophysical and cosmological data, from the Solar system, supernovae Ia, cosmic microwave background radiation and large scale structure constrain it.

We investigate a string-inspired dark energy scenario featuring a scalar field with a coupling to the Gauss-Bonnet invariant. Such coupling can trigger the onset of late dark energy domination after a scaling matter era. The universe may then cross the phantom divide and perhaps also exit from the acceleration. We discuss extensively the cosmological and astrophysical implications of the coupled scalar field. Data from the Solar system, supernovae Ia, cosmic microwave background radiation, large scale structure and big bang nucleosynthesis is used to constrain the parameters of the model. A good Newtonian limit may require to fix the coupling. With all the data combined, there appears to be some tension with the nucleosynthesis bound, and the baryon oscillation scale seems to strongly disfavor the model. These possible problems might be overcome in more elaborate models. In addition, the validity of these constraints in the present context is not strictly established. Evolution of fluctuations in the scalar field and their impact to clustering of matter is studied in detail and more model-independently. Small scale limit is derived for the perturbations and their stability is addressed. A divergence is found and discussed. The general equations for scalar perturbations are also presented and solved numerically, confirming that the Gauss-Bonnet coupling can be compatible with the observed spectrum of cosmic microwave background radiation as well as the matter power spectrum inferred from large scale surveys.

We show that as a result of non-linear self-interactions, scalar field theories that couple to matter much more strongly than gravity are not only viable but could well be detected by a number of future experiments, provided these are properly designed to do so.

Cluster number counts can be used to test dark energy models. We investigate dark energy candidates which are coupled to dark matter. We analyze the cluster number counts dependence on the amount of dark matter coupled to dark energy. Further more, we study how dark energy inhomogeneities affect cluster abundances. It is shown that increasing the coupling reduces significantly the cluster number counts, and that dark energy inhomogeneities increases cluster abundances. Wiggles in cluster number counts are shown to be a specific signature of coupled dark energy models. Future observations will possibly detect such oscillations and discriminate among the different dark energy models.

This thesis describes a detailed investigation of the effects of matter inhomogeneities on the cosmological evolution of the fine structure constant using the Bekenstein-Sandvik-Barrow-Magueijo (BSBM) theory. We briefly review the observational and theoretical motivations to this work, together with the standard cosmological model. We start by analysing the phase space of the system of equations that describes a time-varying fine structure constant, in a homogeneous and isotropic background universe. We classify all the possible behaviours of the fine structure constant in ever-expanding universes and find exact solutions for it. Using a gauge-invariant formalism, we derive and solve the linearly perturbed Einstein cosmological equations for the BSBM theory. We calculate the time evolution of inhomogeneous perturbations of the fine structure constant on small and large scales with respect to the Hubble radius. We also investigate, in the non-linear regime of the large scale structure formation, the space-time evolution of the fine structure constant, inside evolving spherical overdensities. The dependence on the dark-energy equation of state is also analysed. Finally, we analyse the effects of the coupling of the field (that drives the variations in the fine structure constant) to the matter fields, on the space and time evolution of the fine structure constant.

We investigate f(R) theories of gravity within the Palatini approach and show how one can determine the expansion history, H(a), for an arbitrary choice of f(R). As an example, we consider cosmological constraints on such theories arising from the supernova type Ia, large scale structure formation and cosmic microwave background observations. We find that best fit to the data is a non-null leading order correction to the Einstein gravity, but the current data exhibits no significant preference over the concordance LCDM model. Our results show that the often considered 1/R models are not compatible with the data. The results demonstrate that the background expansion alone can act as a good discriminator between modified gravity models when multiple data sets are used.

We investigate the consequences of an imperfect dark energy component on the large scale structure. A phenomenological three parameter fluid description is used to study the effect of dark energy on the cosmic microwave background radiation (CMBR) and matter power spectrum. In addition to the equation of state and the sound speed, we allow a nonzero viscosity parameter for the fluid. Then anisotropic stress perturbations are generated in dark energy. In general, we find that this possibility is not excluded by the present day cosmological observations. In the simplest case when all of the three parameters are constant, we find that the observable effects of the anisotropic stress can be closely mimicked by varying the sound speed of perfect dark energy. However, now also negative values for the sound speed, as expected for adiabatic fluid model, are tolerable and in fact could explain the observed low quadrupole in the CMBR spectrum. We investigate also structure formation of imperfect fluid dark energy characterized by an evolving equation of state. In particular, we study unified models of dark energy with dark matter, such as the Chaplygin gas or the Cardassian expansion, with a shear perturbation included. This can stabilize the growth of inhomogeneities in these models, thus somewhat improving their compatibility with large scale structure observations.

Cosmological consequences of a coupling between massive neutrinos and dark energy are investigated. In such models, the neutrino mass is a function of a scalar field, which plays the role of dark energy. The background evolution, as well as the evolution of cosmological perturbations is discussed. We find that mass-varying neutrinos can leave a significant imprint on the anisotropies in the cosmic microwave background radiation and even lead to a reduction of power on large angular scales.

We investigate how inhomogeneous quintessence models may have a specific signature even in the linear regime of large scale structure formation. The dynamics of the collapse of a dark matter halo is governed by the value or the dynamical evolution of the dark energy equation of state, the energy density's initial conditions and its homogeneity nature in the highly non-linear regime. These have a direct impact on the redshift of collapse, altering in consequence the linearly extrapolated density threshold above which structures will end up collapsing. We compute this quantity for minimally coupled and coupled quintessence models, examining two extreme scenarios: first, when the quintessence field does not exhibit fluctuations on cluster scales and below - homogeneous dark energy; and second, when the field inside the overdensity collapses along with the dark matter - inhomogeneous dark energy. One shows that inhomogeneous dark energy models present distinct features which may be used to confront them with observational data, for instance galaxy number counting. Fitting formulae for the linearly extrapolated density threshold above which structures will end up collapsing are provided for models of dark energy with constant equation of state.

A relativistic theory of modified gravity has been recently proposed by Bekenstein. The tensor field in Einstein's theory of gravity is replaced by a scalar, a vector, and a tensor field which interact in such a way to give Modified Newtonian Dynamics (MOND) in the weak-field non-relativistic limit. We study the evolution of the universe in such a theory, identifying its key properties and comparing it with the standard cosmology obtained in Einstein gravity. The evolution of the scalar field is akin to that of tracker quintessence fields. We expand the theory to linear order to find the evolution of perturbations on large scales. The impact on galaxy distributions and the cosmic microwave background is calculated in detail. We show that it may be possible to reproduce observations of the cosmic microwave background and galaxy distributions with Bekenstein's theory of MOND.

Motivated by the string landscape we examine scenarios for which inflation is a two-step process, with a comparatively short inflationary epoch near the string scale and a longer period at a much lower energy (like the TeV scale). We quantify the number of $e$-foldings of inflation which are required to yield successful inflation within this picture. The constraints are very sensitive to the equation of state during the epoch between the two inflationary periods, as the extra-horizon modes can come back inside the horizon and become reprocessed. We find that the number of $e$-foldings during the first inflationary epoch can be as small as 12, but only if the inter-inflationary period is dominated by a network of cosmic strings (such as might be produced if the initial inflationary period is due to the brane-antibrane mechanism). In this case a further 20 $e$-foldings of inflation would be required at lower energies to solve the late universe's flatness and horizon problems.

Spatially flat loitering universe models have recently been shown to arise in the context of brane world scenarios. Such models allow more time for structure formation to take place at high redshifts, easing, e.g., the tension between the observed and predicted evolution of the quasar population with redshift. While having the desirable effect of boosting the growth of structures, we show that in such models the position of the first peak in the power spectrum of the cosmic microwave background anisotropies severely constrains the amount of loitering at high redshifts.

We study the inhomogeneous cosmological evolution of the Newtonian gravitational 'constant' G in the framework of scalar-tensor theories. We investigate the differences that arise between the evolution of G in the background universes and in local inhomogeneities that have separated out from the global expansion. Exact inhomogeneous solutions are found which describe the effects of masses embedded in an expanding FRW Brans-Dicke universe. These are used to discuss possible spatial variations of G in different regions. We develop the technique of matching different scalar-tensor cosmologies of different spatial curvature at a boundary. This provides a model for the linear and non-linear evolution of spherical overdensities and inhomogeneities in G. This allows us to compare the evolution of G and \dot{G} that occurs inside a collapsing overdense cluster with that in the background universe. We develop a simple virialisation criterion and apply the method to a realistic lambda-CDM cosmology containing spherical overdensities. Typically, far slower evolution of \dot{G} will be found in the bound virialised cluster than in the cosmological background. We consider the behaviour that occurs in Brans-Dicke theory and in some other representative scalar-tensor theories.

Dark energy might have an influence on the formation of non--linear structures during the cosmic history. For example, in models in which dark energy couples to dark matter, it will be non--homogeneous and will influence the collapse of a dark matter overdensity. We use the spherical collapse model to estimate how much influence dark energy might have.

We study the spherical collapse model in dark energy cosmologies, in which dark energy is modelled as a minimally coupled scalar field. We first follow the standard assumption that dark energy does not cluster on the scales of interest. Investigating four different popular potentials in detail, we show that the predictions of the spherical collapse model depend on the potential used. We also investigate the dependence on the initial conditions. Secondly, we investigate in how far perturbations in the quintessence field affect the predictions of the spherical collapse model. In doing so, we assume that the field collapses along with the dark matter. Although the field is still subdominant at the time of virialisation, the predictions are different from the case of a homogeneous dark energy component. This will in particular be true if the field is non--minimally coupled. We conclude that a better understanding of the evolution of dark energy in the highly non--linear regime is needed in order to make predictions using the spherical collapse model in models with dark energy.

We study the space-time evolution of the fine structure constant, $\alpha$, inside evolving spherical overdensities in a lambda-CDM Friedmann universe using the spherical infall model. We show that its value inside virialised regions will be significantly larger than in the low-density background universe. The consideration of the inhomogeneous evolution of the universe is therefore essential for a correct comparison of extragalactic and solar system limits on, and observations of, possible time variation in $\alpha$ and other constants. Time variation in $\alpha$ in the cosmological background can give rise to no locally observable variations inside virialised overdensities like the one in which we live, explaining the discrepancy between astrophysical and geochemical observations.

Using the BSBM varying-alpha theory, with dark matter dominated by magnetic energy, and the spherical collapse model for cosmological structure formation, we have studied the effects of the dark-energy equation of state and the coupling of alpha to the matter fields on the space and time evolution of alpha. We have compared its evolution inside virialised overdensities with that in the cosmological background, using the standard ($% \Lambda =0$) $CDM$ model of structure formation and the dark-energy modification, $wCDM$. We find that, independently of the model of structure formation one considers, there is always a difference between the value of alpha in an overdensity and in the background. In a $SCDM$ model, this difference is the same, independent of the virialisation redshift of the overdense region. In the case of a $wCDM$ model, especially at low redshifts, the difference depends on the time when virialisation occurs and the equation of state of the dark energy. At high redshifts, when the $wCDM$ model becomes asymptotically equivalent to the $SCDM$ one, the difference is constant. At low redshifts, when dark energy starts to dominate the cosmological expansion, the difference between alpha in a cluster and in the background grows. The inclusion of the effects of inhomogeneity leads naturally to no observable local time variations of alpha on Earth and in our Galaxy even though time variations can be significant on quasar scales. The inclusion of the effects of inhomogeneous cosmological evolution are necessary if terrestrial and solar-system bounds on the time variation of the fine structure 'constant' are to be correctly compared with extragalactic data.

Using a gauge-invariant formalism we derive and solve the perturbed cosmological equations for the BSBM theory of varying fine structure 'constant'. We calculate the time evolution of inhomogeneous perturbations of the fine structure constant, $\frac{\delta \alpha}{\alpha}$ on small and large scales with respect to the Hubble radius. In a radiation-dominated universe small inhomogeneities in $\alpha $ will decrease on large scales but on scales smaller than the Hubble radius they will undergo stable oscillations. In a dust-dominated universe small inhomogeneous perturbations in $\alpha $ will become constant on large scales and on small scales they will increase as $t^{2/3}$, and $\frac{\delta \alpha}{\alpha}$ will track $% \frac{\delta \rho _{m}}{\rho_{m}}$ . If the expansion accelerates, as in the case of a $\Lambda $ or quintessence-dominated phase, inhomogeneities in $\alpha $ will decrease on both large and small scales. The amplitude of perturbations in $\alpha $ will be much smaller than that of matter or radiation perturbations. We also present a numerical study of the non-linear evolution of spherical inhomogeneities in radiation and dust universes by means of a comparison between the evolution of flat and closed Friedmann models with time-varying $\alpha .$Various limitations of these simple models are also discussed.

Assuming a Friedmann universe which evolves with a power-law scale factor, $a=t^{n}$, we analyse the phase space of the system of equations that describes a time-varying fine structure 'constant', $\alpha$, in the Bekenstein-Sandvik-Barrow-Magueijo generalisation of general relativity. We have classified all the possible behaviours of $\alpha (t)$ in ever-expanding universes with different $n$ and find new exact solutions for $\alpha (t)$. We find the attractors points in the phase space for all $n$. In general, $\alpha $ will be a non-decreasing function of time that increases logarithmically in time during a period when the expansion is dust dominated ($n=2/3$), but becomes constant when $n>2/3$. This includes the case of negative-curvature domination ($n=1$). $\alpha $ also tends rapidly to a constant when the expansion scale factor increases exponentially. A general set of conditions is established for $\alpha $ to become asymptotically constant at late times in an expanding universe.

Spontaneous breaking of Lorentz invariance may take place in string theories, possibly endowing the photon with a mass. This leads to the breaking of the conformal symmetry of the electromagnetic action allowing for the generation within inflationary scenarios of magnetic fields over $Mpc$ scales. We show that the generated fields are consistent with amplification by the galactic dynamo processes and can be as large as to explain the observed galactic magnetic fields through the collapse of protogalactic clouds.