While the accelerated cosmological expansion suggests the presence of a dark energy component, the possibility of its interaction with other fluid components is not entirely ruled out. We investigate the cosmological dynamics of coupled quintessence (where the quintessence couples and interacts with a generic fluid) in a disformal framework. The novelty of this work is: (i) extend the idea of coupled quintessence mimicking the cosmological constant with an uncoupled fluid, by including a disformal coupling, and (ii) to bring these conformal and disformal couplings with a LCDM background, to a dynamical system analysis. Using this approach, the quintessence potential does not have to be specified. Three scenarios are studied. Scenario I consists of a pressureless fluid (e.g dark matter), which is coupled to quintessence, and both fluids are mimicking exactly a LCDM background. Scenario II consists of a relativistic fluid, which is coupled to quintessence and both fluids are also mimicking exactly the LCDM background. Scenario III consists of a relativistic fluid (e.g. neutrinos), which is coupled to quintessence and both fluids are mimicking a background with a cosmological constant and uncoupled radiation. The expansion history of the coupled quintessence is investigated for scenario I, which confirms that the mass scale of the quintessence influences the disformal characteristics of the dynamical system, which is portrayed by the evolution of an effective conformal coupling.

The key probes of the growth of large-scale structure are its rate $f$ and amplitude $\sigma_8$. Redshift space distortions in the galaxy power spectrum allow us to measure only the combination $f\sigma_8$, which can be used to constrain the standard cosmological model or alternatives. By using measurements of the galaxy-galaxy lensing cross-correlation spectrum or of the galaxy bispectrum, it is possible to break the $f\sigma_8$ degeneracy and obtain separate estimates of $f$ and $\sigma_8$ from the same galaxy sample. Currently there are only a handful of such separate measurements, but even this allows for improved constraints on cosmological models. We consider what can be achieved by a future nominal sample that delivers a $\sim 1\%$ constraint on $f$ and $\sigma_8$ separately, compared to the case with a similar precision on the combination $f\sigma_8$. For the six cosmological parameters of $\Lambda$CDM, we find improvements of $\sim\! 5$--$50\%$ on their constraints. For modified gravity models in the Horndeski class, the improvements on these standard parameters are $\sim\! 0$--$15\%$. However, the precision on the sum of neutrino masses improves by 65\% and there is a significant increase in the precision on the background and perturbation Horndeski parameters.

In the era of radio astronomy, the high sensitivity of the Square Kilometre Array (SKA) could play a decisive role in the detection of new radio sources. In this work, we study the SKA sensitivity to the synchrotron radio emission expected by the annihilation of TeV DM candidate in the Draco dwarf spheroidal galaxy. On the one hand, we consider model-independent DM candidates: we find out that with 1000 hours of data-taking, SKA1-MID will be able to exclude up to 10 TeV thermal DM candidates that annihilate in $W^+W^-$ and $b\bar b$ channels. We also study as these constraints improve by including a density enhancement due to a DM-spike associated with an intermediate-mass black hole in Draco. On the other hand, we consider extra-dimensional brane-world DM candidates, dubbed branons. In this specific scenario, SKA allows us to set constraints on the branon parameter space ($f$,$M$), where $f$ is related to the coupling of the branon to the Standard Model particles and $M$ is the mass of the branon itself. In particular, we consider two different branon DM candidates. We find out that SKA will be able to set more stringent constraints on the branon DM candidate required in order to fit the AMS-02 data, yet the sensitivity of the instrument should be improved in order to study the branon candidate for the Galactic Centre. Nonetheless, we show that SKA represents - among other detectors - the most promising instrument for multi-wavelength detection of synchrotron radio emission by annihilating multi-TeV DM.

We use growth of structure data to constrain the effective field theory of dark energy. Considering as case study Horndeski theories with the speed of gravitational waves equal to that of light, we show how constraints on the free parameters and the large-scale structure phenomenological functions can be improved by two ingredients: firstly by complementing the set of redshift-space distortions data with the three recent measurements of the growth rate $f$ and the amplitude of matter fluctuations $\sigma_8$ from the VIPERS and SDSS collaborations; secondly by applying a local Solar System bound on the variation of the Newton constant. This analysis allows us to conclude that: $i)$ despite firmly restricting the predictions of weaker gravity, the inclusion of the Solar System bound does not prevent suppressed growth relative to the standard model $\Lambda$CDM at low redshifts; $ii)$ the same bound in conjunction with the growth of structure data strongly restricts the redshift evolution of the gravitational slip parameter to be close to unity and the present value is constrained to one at the $10^{-3}$ level; $iii)$ the growth of structure data favours a fifth force contribution to the effective gravitational coupling at low redshifts and at more than one sigma at present time.

In the following work, we compute the positron production from branon dark matter annihilations in order to constrain extra-dimensional theories. By having assumed that the positron fraction measured by AMS-02 is well explained just with astrophysical sources, exclusion diagrams for the branon mass and the tension of the brane, the two parameters characterising the branon phenomenology become possible. Our analysis has been performed for a minimal and a medium diffusion model in one extra dimension for both pseudo-isothermal and Navarro Frenck White dark matter haloes. Our constraints in the dark matter mass candidate range between 200 GeV and 100 TeV. Combined with previous cosmological analyses and experimental data in colliders, it allows us to set bounds on the parameter space of branons. In particular, we have discarded regions in the mass-tension diagram up to a branon mass of 28 TeV for the pseudo-Isothermal prole and minimal diffusion, and 63 TeV for the Navarro-Frenck-White profile and medium diffusion.

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.

In the context of f(R)=R + alpha R^2 gravity, we study the existence of neutron and quark stars with no intermediate approximations in the generalised system of Tolman-Oppenheimer-Volkov equations. Analysis shows that for positive alpha's the scalar curvature does not drop to zero at the star surface (as in General Relativity) but exponentially decreases with distance. Also the stellar mass bounded by star surface decreases when the value alpha increases. Nonetheless distant observers would observe a gravitational mass due to appearance of a so-called gravitational sphere around the star. The non-zero curvature contribution to the gravitational mass eventually is shown to compensate the stellar mass decrease for growing alpha's. We perform our analysis for several equations of state including purely hadronic configurations as well as hyperons and quark stars. In all cases, we assess that the relation between the parameter $\alpha$ and the gravitational mass weakly depend upon the chosen equation of state. Another interesting feature is the increase of the star radius in comparison to General Relativity for stars with masses close to maximal, whereas for intermediate masses around 1.4-1.6 solar masses, the radius of star depends upon alpha very weakly. Also the decrease in the mass bounded by star surface may cause the surface redshift to decrease in R^2-gravity when compared to Einsteinian predictions. This effect is shown to hardly depend upon the observed gravitational mass. Finally, for negative values of alpha our analysis shows that outside the star the scalar curvature has damped oscillations but the contribution of the gravitational sphere into the gravitational mass increases indefinitely with radial distance putting into question the very existence of such relativistic stars.

In the context of f(R) gravity theories, we show that the apparent mass of a neutron star as seen from an observer at infinity is numerically calculable but requires careful matching, first at the star's edge, between interior and exterior solutions, none of them being totally Schwarzschild-like but presenting instead small oscillations of the curvature scalar R; and second at large radii, where the Newtonian potential is used to identify the mass of the neutron star. We find that for the same equation of state, this mass definition is always larger than its general relativistic counterpart. We exemplify this with quadratic $R^2$ and Hu-Sawicki-like modifications of the standard General Relativity action. Therefore, the finding of two-solar mass neutron stars basically imposes no constraint on stable f(R) theories. However, star radii are in general smaller than in General Relativity, which can give an observational handle on such classes of models at the astrophysical level. Both larger masses and smaller matter radii are due to much of the apparent effective energy residing in the outer metric for scalar-tensor theories. Finally, because the f(R) neutron star masses can be much larger than General Relativity counterparts, the total energy available for radiating gravitational waves could be of order several solar masses, and thus a merger of these stars constitutes an interesting wave source.

Model-independent methods in cosmology has become an essential tool in order to deal with an increasing number of theoretical alternatives for explaining the late-time acceleration of the Universe. In principle, this provides a way of testing the Cosmological Concordance (or $\Lambda$CDM) model under different assumptions and ruling out whole classes of competing theories. One such model-independent method is the so-called cosmographic approach, which relies only on the homogeneity and isotropy of the Universe on large scales. We show that this method suffers from many shortcomings, providing biased results depending on the auxiliary variable used in the series expansion and is unable to rule out models or adequately reconstruct theories with higher-order derivatives in either the gravitational or matter sector. Consequently, in its present form, this method seems unable to provide reliable or useful results for cosmological applications.

Modified gravity has attracted much attention over the last few years and remains a potential candidate for dark energy. In particular, the so-called viable f(R) gravity theories, which are able to both recover General Relativity (GR) and produce late-time cosmic acceleration, have been widely studied in recent literature. Nevertheless, extended theories of gravity suffer from several shortcomings which compromise their ability to provide realistic alternatives to the standard cosmological Lambda CDM Concordance model. We address the existence of cosmological singularities and the conditions that guarantee late-time acceleration,assuming reasonable energy conditions for standard matter in the so-called Hu-Sawicki f(R) model, currently among the most widely studied modifications to General Relativity. Then using the Supernovae Ia Union 2.1 catalogue, we further constrain the free parameters of this model. The combined analysis of both theoretical and observational constraints sheds some light on the viable parameter space of these models and the form of the underlying effective theory of gravity.

We present a complete analysis of the imprint of tensor anisotropies on the Cosmic Microwave Background for a class of f(R) gravity theories within the PPF-CAMB framework. We derive the equations, both for the cosmological background and gravitational wave perturbations, required to obtain the standard temperature and polarization power spectra, taking care to include all effects which arise from f(R) modifications of both the background and the perturbation equations. For R^n gravity, we show that for n different from 2, the initial conditions in the radiation dominated era are the same as those found in General Relativity. We also find that by doing simulations which involve either modifying the background evolution while keeping the perturbation equations fixed or fixing the background to be the Lambda-CDM model and modifying the perturbation equations, the dominant contribution to deviations from General Relativity in the temperature and polarization spectra can be attributed to modifications in the background. This demonstrates the importance of using the correct background in perturbative studies of f(R) gravity. Finally an enhancement in the B-modes power spectra is observed which may allow for lower inflationary energy scales.

One possible explanation for the present observed acceleration of the Universe is the breakdown of homogeneity and isotropy due to the formation of non-linear structures. How inhomogeneities affect the averaged cosmological expansion rate and lead to late-time acceleration is generally considered to be due to some backreaction mechanism. General Relativity together with pressure-free matter have until recently been considered as the sole ingredients for averaged calculations. In this communication we focus our attention on more general scenarios, including imperfect fluids as well as alternative theories of gravity, and apply an averaging procedure to them in order to determine possible backreaction effects. For illustrative purposes, we present our results for dark energy models, quintessence and Brans-Dicke theories. We also provide a discussion about the limitations of frame choices in the averaging procedure.

We present a study of the Galactic Center region as a possible source of both secondary gamma-ray and neutrino fluxes from annihilating dark matter. We have studied the gamma-ray flux observed by the High Energy Stereoscopic System (HESS) from the J1745-290 Galactic Center source. The data are well fitted as annihilating dark matter in combination with an astrophysical background. The analysis was performed by means of simulated gamma spectra produced by Monte Carlo event generators packages. We analyze the differences in the spectra obtained by the various Monte Carlo codes developed so far in particle physics. We show that, within some uncertainty, the HESS data can be fitted as a signal from a heavy dark matter density distribution peaked at the Galactic Center, with a power-law for the background with a spectral index which is compatible with the Fermi-Large Area Telescope (LAT) data from the same region. If this kind of dark matter distribution generates the gamma-ray flux observed by HESS, we also expect to observe a neutrino flux. We show prospective results for the observation of secondary neutrinos with the Astronomy with a Neutrino Telescope and Abyss environmental RESearch project (ANTARES), Ice Cube Neutrino Observatory (Ice Cube) and the Cubic Kilometer Neutrino Telescope (KM3NeT). Prospects solely depend on the device resolution angle when its effective area and the minimum energy threshold are fixed.

We study the differences in the gamma ray spectra simulated by four Monte Carlo event generator packages developed in particle physics. Two different versions of PYTHIA and two of HERWIG are analyzed, namely PYTHIA 6.418 and HERWIG 6.5.10 in Fortran and PYTHIA 8.165 and HERWIG 2.6.1 in C++. For all the studied channels, the intrinsic differences between them are shown to be significative and may play an important role in misunderstanding dark matter signals.

In the context of f(R,T) theories of gravity, we study the evolution of scalar cosmological perturbations in the metric formalism. According to restrictions on the background evolution, a specific model within these theories is assumed in order to guarantee the standard continuity equation. Using a completely general procedure, we find the complete set of differential equations for the matter density perturbations. In the case of sub-Hubble modes, the density contrast evolution reduces to a second-order equation. We show that for well-motivated f(R,T) Lagrangians the quasistatic approximation yields to very different results from the ones derived in the frame of the concordance {\Lambda}CDM model constraining severely the viability of such theories.

An arguable aspect of the modified gravity theories is that many of them present the so-called degeneracy problem. For instance, the cosmological evolution, gravitational collapse and the main features of standard black-hole configurations, can be mimicked by many of those theories. In this communication we revise briefly the appropriate observable quantities to be measured in order to discard alternative theories to LCDM, such as the observed growth of scalar perturbations with Sloan data and the CMB tensor perturbations evolution.

Along this review, we focus on the study of several properties of modified gravity theories, in particular on black-hole solutions and its comparison with those solutions in General Relativity, and on Friedmann-Lemaitre-Robertson-Walker metrics. The thermodynamical properties of fourth order gravity theories are also a subject of this investigation with special attention on local and global stability of paradigmatic f(R) models. In addition, we revise some attempts to extend the Cardy-Verlinde formula, including modified gravity, where a relation between entropy bounds is obtained. Moreover, a deep study on cosmological singularities, which appear as a real possibility for some kind of modified gravity theories, is performed, and the validity of the entropy bounds is studied.

If the present dark matter in the Universe annihilates into Standard Model particles, it must contribute to the gamma ray fluxes detected on the Earth. Here we briefly review the present constraints for the detection of gamma ray photons produced in the annihilation of branon dark matter. We show that observations of dwarf spheroidal galaxies and the galactic center by EGRET, Fermi-LAT or MAGIC are below the sensitivity limits for branon detection. However,future experiments such as CTA could be able to detect gamma-ray photons from annihilating branons of masses above 150 GeV.

Branons are new degrees of freedom that appear in flexible brane-world models corresponding to brane fluctuations. These new fields can behave as standard weakly interacting massive particles (WIMPs) with a significant associated thermal relic density. We analyze the present constraints from their spontaneous annihilations into photons for EGRET, Fermi-LAT and MAGIC, and the prospects for detection in future Cherenkov telescopes. In particular, we focus on possible signals coming from the Galactic Center and different dwarf spheroidals, such as Draco, Sagittarius, Canis Major and SEGUE 1. We conclude that for those targets, present observations are below the sensitivity limits for branon detection by assuming standard dark matter distributions and no additional boost factors. However, future experiments such as CTA could be able to detect gamma-ray photons coming from the annihilation of branons with masses higher than 150 GeV.

Annihilation of different dark matter (DM) candidates into Standard Model (SM) particles could be detected through their contribution to the gamma ray fluxes that are measured on the Earth. The magnitude of such contributions depends on the particular DM candidate, but certain imprints of produced photon spectra may be analyzed in a model-independent fashion. In this work we provide the fitting formulae for the photon spectra generated by WIMP annihilation into quarks, leptons and gauge bosons channels in a wide range of WIMP masses.

In this work we provide the fitting formula valid for the simulated photon spectra from WIMP annihilation into light quark-anti quark (qq-) channels in a wide range of WIMP masses. We illustrate our results for the cc- channel.

If the present dark matter in the Universe annihilates into Standard Model particles, it must contribute to the fluxes of cosmic rays that are detected on the Earth, and in particular, to the observed gamma ray fluxes. The magnitude of such contribution depends on the particular dark matter candidate, but certain features of the produced photon spectra may be analyzed in a rather model-independent fashion. In this work we provide the complete photon spectra coming from WIMP annihilation into Standard Model particle-antiparticle pairs obtained by extensive Monte Carlo simulations. We present results for each individual annihilation channel and provide analytical fitting formulae for the different spectra for a wide range of WIMP masses.

This thesis will try to contribute to the understanding of open issues in cosmology by considering f(R) and brane-world theories. In Chapter 1, we shall summarize the main features of f(R) gravities in the metric formalism and we shall introduce both the notion of brane excitations, the branons, and the brane-skyrmions. We shall finish the chapter by providing some insight about the possibility of mini black holes detection in the LHC as a signature for the validity of these modified gravity theories. The Chapter 2 will deal with f(R) theories able to mimic Einstein-Hilbert plus cosmological constant solutions and f(R) theories will be shown to be able to mimic the cosmological evolution generated by any perfect fluid with constant equation of state. The Chapter 3 will be devoted to the computation of cosmological perturbations for f(R) theories. Special attention will be paid here to obtain a completely general differential equation for the evolution of perturbations and its particularization for the so-called sub-Hubble scales will be explicitly shown. In the Chapter 4 we shall focus on the study of black holes in f(R) gravities in an arbitrary number of dimensions. With this purpose we shall study constant curvature solutions for f(R) theories as well as perturbative solutions around the standard SAdS geometry. An important part of this chapter will be then devoted to the thermodynamics of SAdS black holes in f(R) theories. In the Chapter 5 we will thoroughly study brane-skyrmions. In this context, the recent claim of detection of an unexpected feature in the CMB, referred to as the cold spot, will be explained as a topological defect on the brane in complete agreement with those calculations in the literature that tried to explain that cold spot as a texture of a NLSM. Main conclusions are summarized all together in Chapter 6 .

A modified f(R) gravity model has been recently proposed in Phys. Rev. Lett. 102, 221101 (2009) [arXiv:0905.1941] whose cosmological behaviour is clearly distinguishable from LCDM. Contrary to previous opinions which consider that self-consistent f(R) gravity models distinct from LCDM are almost ruled out, the authors claim that the proposed model is cosmologically viable. Here we show that although the model satisfies some consistency conditions, precisely because of its departure from LCDM behaviour, it does not satisfy local gravity constraints and, in addition, the predicted matter power spectrum conflicts with SDSS data.

In the context of f(R) theories of gravity, we study the cosmological evolution of scalar perturbations by using a completely general procedure. We find that the exact fourth-order differential equation for the matter density perturbations in the longitudinal gauge, reduces to a second-order equation for sub-Hubble modes. This simplification is compared with the standard (quasi-static) equation used in the literature. We show that for general f(R) functions the quasi-static approximation is not justified. However for those f(R) adequately describing the present phase of accelerated expansion and satisfying local gravity tests, it does give a correct description for the evolution of perturbations.

One of the most striking features found in the cosmic microwave background data is the presence of an anomalous Cold Spot (CS) in the temperature maps made by the Wilkinson Microwave Anisotropy Probe (WMAP). This CS has been interpreted as the result of the presence of a collapsing texture, perhaps coming from some early universe Grand Unified Theory (GUT) phase transition. In this work we propose an alternative explanation based on a completely different kind of texture which appears in a natural way in a broad class of brane-world models. This type of textures known as brane-skyrmions can be understood as holes in the brane which make possible to pass through them along the extra-dimensional space. The typical scales needed for the proposed brane-skyrmions to correctly describe the observed CS can be as low as the electroweak scale.

In the context of $f(R)$ theories of gravity, we study the evolution of scalar cosmological perturbations in the metric formalism. Using a completely general procedure, we find the exact fourth-order differential equation for the matter density perturbations in the longitudinal gauge. In the case of sub-Hubble modes, the expression reduces to a second-order equation which is compared with the standard (quasi-static) equation used in the literature. We show that for general $f(R)$ functions the quasi-static approximation is not justified. However, for those functions adequately describing the present phase of accelerated expansion and satisfying local gravity tests, it provides a correct description for the evolution of perturbations.

It has been suggested that our universe could be a 3-dimensional brane where the SM fields live embedded in a D-dimensional space-time. In flexible raneworlds, in addition to the SM fields, new degrees of freedom appear on the brane associated to brane fluctuations. These new fields, known as branons, are standard WIMPs (Weakly Interacting Massive Particles) and therefore natural dark matter candidates, whose spontaneous annihilations can provide first evidences for this scenario.