A bouncing universe is a viable candidate to solve the initial singularity problem. Here we consider bouncing solutions in the context of $f(R,\mathcal{G})$ gravity by using an order reduction technique which allows one to find solutions that are perturbatively close to General Relativity. This procedure also acts as a model selection approach. Indeed, several covariant gravitational actions leading to a bounce are directly selected by demanding that the Friedmann equation derived within such gravity theories coincides with the one emerging from Loop Quantum Cosmology.

We study non-rotating and isotropic strange quark stars in Lorentz-violating theories of gravity, and in particular in Ho\v{r}ava gravity and Einstein-{\ae}ther theory. For quark matter we adopt both linear and non-linear equations of state, corresponding to the MIT bag model and colour flavour locked state, respectively. The new structure equations describing hydrostatic equilibrium generalize the usual Tolman-Oppenheimer-Volkoff (TOV) equations of Einstein's general relativity. A dimensionless parameter $\nu$ measures the deviation from the standard TOV equations, which are recovered in the limit $\nu \rightarrow 0$. We compute the mass, the radius as well as the compactness of the stars, and we show graphically the impact of the parameter $\nu$ on the mass-to-radius profiles for different equations of state describing quark matter. The energy conditions and stability criteria are also considered, and they are all found to be fulfilled.

Cyclic universes with bouncing solutions are candidates for solving the big bang initial singularity problem. Here we seek bouncing solutions in a modified Gauss-Bonnet gravity theory, of the type $R+f(G)$, where $R$ is the Ricci scalar, $G$ is the Gauss-Bonnet term, and $f$ some function of it. In finding such a bouncing solution we resort to a technique that reduces the order of the differential equations of the $R+f(G)$ theory to second order equations. As general relativity is a theory whose equations are of second order, this order reduction technique enables one to find solutions which are perturbatively close to general relativity. We also build the covariant action of the order reduced theory.

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 quadratic-order degenerate higher-order scalar-tensor (DHOST) theories compatible with gravitational-wave constraints, we derive the most general Lagrangian allowing for tracker solutions characterized by $\dot{\phi}/H^p={\rm constant}$, where $\dot{\phi}$ is the time derivative of a scalar field $\phi$, $H$ is the Hubble expansion rate, and $p$ is a constant. While the tracker is present up to the cubic-order Horndeski Lagrangian $L=c_2X-c_3X^{(p-1)/(2p)} \square \phi$, where $c_2, c_3$ are constants and $X$ is the kinetic energy of $\phi$, the DHOST interaction breaks this structure for $p \neq 1$. Even in the latter case, however, there exists an approximate tracker solution in the early cosmological epoch with the nearly constant field equation of state $w_{\phi}=-1-2p\dot{H}/(3H^2)$. The scaling solution, which corresponds to $p=1$, is the unique case in which all the terms in the field density $\rho_{\phi}$ and the pressure $P_{\phi}$ obey the scaling relation $\rho_{\phi} \propto P_{\phi} \propto H^2$. Extending the analysis to the coupled DHOST theories with the field-dependent coupling $Q(\phi)$ between the scalar field and matter, we show that the scaling solution exists for $Q(\phi)=1/(\mu_1 \phi+\mu_2)$, where $\mu_1$ and $\mu_2$ are constants. For the constant $Q$, i.e., $\mu_1=0$, we derive fixed points of the dynamical system by using the general Lagrangian with scaling solutions. This result can be applied to the model construction of late-time cosmic acceleration preceded by the scaling $\phi$-matter-dominated epoch.

We generalise the covariant Tolman-Oppenheimer-Volkoff equations proposed in arXiv:1709.02818 [gr-qc] to the case of static and spherically symmetric spacetimes with anisotropic sources. The extended equations allow a detailed analysis of the role of the anisotropic terms in the interior solution of relativistic stars and lead to the generalisation of some well known solutions of this type. We show that, like in the isotropic case, one can define generating theorems for the anisotropic Tolman-Oppenheimer-Volkoff equations. We also find that it is possible to define a reconstruction algorithm able to generate a double infinity of interior solutions. Among these, we derive a class of solutions that can represent "quasi-isotropic" stars.

We construct a covariant version of the Tolman-Oppenheimer-Volkoff equations in the case of isotropic sources. The new equations make evident the mathematical problems in the determination of interior solutions of relativistic stellar objects. Using a reconstruction algorithm we find two physically interesting generalisations of previously known stellar interior solutions. The variables that we use also allow an easier formulation of known generating theorems for solutions associated to relativistic stellar objects.

We consider Horava gravity within the framework of the effective field theory (EFT) of dark energy and modified gravity. We work out a complete mapping of the theory into the EFT language for an action including all the operators which are relevant for linear perturbations with up to sixth order spatial derivatives. We then employ an updated version of the EFTCAMB/EFTCosmoMC package to study the cosmology of the low-energy limit of Horava gravity and place constraints on its parameters using several cosmological data sets. In particular we use cosmic microwave background (CMB) temperature-temperature and lensing power spectra by Planck 2013, WMAP low-l polarization spectra, WiggleZ galaxy power spectrum, local Hubble measurements, Supernovae data from SNLS, SDSS and HST and the baryon acoustic oscillations measurements from BOSS, SDSS and 6dFGS. We get improved upper bounds, with respect to those from Big Bang Nucleosynthesis, on the deviation of the cosmological gravitational constant from the local Newtonian one. At the level of the background phenomenology, we find a relevant rescaling of the Hubble rate at all epoch, which has a strong impact on the cosmological observables; at the level of perturbations, we discuss in details all the relevant effects on the observables and find that in general the quasi-static approximation is not safe to describe the evolution of perturbations. Overall we find that the effects of the modifications induced by the low-energy Horava gravity action are quite dramatic and current data place tight bounds on the theory parameters.

Modifications of General Relativity usually include extra dynamical degrees of freedom, which to date remain undetected. Here we explore the possibility of modifying Einstein's theory by adding solely nondynamical fields. With the minimal requirement that the theory satisfies the weak equivalence principle and admits a covariant Lagrangian formulation, we show that the field equations generically have to include higher-order derivatives of the matter fields. This has profound consequences for the viability of these theories. We develop a parametrization based on a derivative expansion and show that - to next to leading order - all theories are described by just two parameters. Our approach can be used to put stringent, theory-independent constraints on such theories, as we demonstrates using the Newtonian limit as an example.

Dwarf galaxies are good candidates to investigate the nature of Dark Matter, because their kinematics are dominated by this component down to small galactocentric radii. We present here the results of detailed kinematic analysis and mass modelling of the Orion dwarf galaxy, for which we derive a high quality and high resolution rotation curve that contains negligible non-circular motions and we correct it for the asymmetric drift. Moreover, we leverage the proximity (D = 5.4 kpc) and convenient inclination (47{\deg}) to produce reliable mass models of this system. We find that the Universal Rotation Curve mass model (Freeman disk + Burkert halo + gas disk) fits the observational data accurately. In contrast, the NFW halo + Freeman disk + gas disk mass model is unable to reproduce the observed Rotation Curve, a common outcome in dwarf galaxies. Finally, we attempt to fit the data with a MOdified Newtonian Dynamics (MOND) prescription. With the present data and with the present assumptions on distance, stellar mass, constant inclination and reliability of the gaseous mass, the MOND "amplification" of the baryonic component appears to be too small to mimic the required "dark component". The Orion dwarf reveals a cored DM density distribution and a possible tension between observations and the canonical MOND formalism.

In spiral galaxies, we explain their non-Keplerian rotation curves (RCs) by means of a non-luminous component embedding their stellar-gaseous disks. Understanding the detailed properties of this component (labelled Dark Matter, DM) is one of the most pressing issues of Cosmology. We investigate the recent relationship (claimed by Walker et al. 2010, hereafter W+10) between $r$, the galaxy radial coordinate, and $V_h(r)$, the dark halo contribution to the circular velocity at $r$, {\it a}) in the framework of the Universal Rotation Curve (URC) paradigm and directly {\it b}) by means of the kinematics of a large sample of DM dominated spirals. We find a general agreement between the W+10 claim, the distribution of DM emerging from the URC and that inferred in the (low luminosity) objects of our sample. We show that such a phenomenology, linking the spiral's luminosity, radii and circular velocities, implies an evident inconsistency with (naive) predictions in the $\Lambda$ Cold Dark Matter ($\Lambda$CDM) scenario.

Spherical symmetry for f(R)-gravity is discussed by searching for Noether symmetries. The method consists in selecting conserved quantities in form of currents that reduce dynamics of f(R)-models compatible with symmetries. In this way we get a general method to obtain constants of motion without setting a priori the form of f(R). In this sense, the Noether symmetry results a physical criterium. Relevant cases are discussed.

We attempt a critical reconsideration of "detailed balance" as a principle that can be used to restrict the proliferation of couplings in Horava-Lifshitz gravity. We re-examine the shortcomings that have been usually associated with it in the literature and we argue that easy remedies can be found for all of them within the framework of detailed balance, and that the most persistent of them are actually related to projectability. We show that, once projectability is abandoned, detailed balance reduces the number of independent couplings by roughly an order of magnitude and imposes only one restriction that constitutes a phenomenological concern: the size of the (bare) cosmological constant is unacceptably large. Remarkably, this restriction (which is present in the projectable version as well) has been so far under-appreciated in the literature. Optimists might prefer to interpret it as a potential blessing in disguise, as it allows one to entertain the idea of a miraculous cancelation between the bare cosmological constant and the (still poorly understood) vacuum energy contribution.

We discuss the possibility to obtain an electromagnetic emission accompanying the gravitational waves emitted in the coalescence of a compact binary system. Motivated by the existence of black hole configurations with open magnetic field lines along the rotation axis, we consider a magnetic dipole in the system, the evolution of which leads to (i) electromagnetic radiation, and (ii) a contribution to the gravitational radiation, the luminosity of both being evaluated. Starting from the observations on magnetars, we impose upper limits for both the electromagnetic emission and the contribution of the magnetic dipole to the gravitational wave emission. Adopting this model for the evolution of neutron star binaries leading to short gamma ray bursts, we compare the correction originated by the electromagnetic field to the gravitational waves emission, finding that they are comparable for particular values of the magnetic field and of the orbital radius of the binary system. Finally we calculate the electromagnetic and gravitational wave energy outputs which result comparable for some values of magnetic field and radius.

We discuss the conformal symmetry between Jordan and Einstein frames considering their relations with the metric and Palatini formalisms for modified gravity. Appropriate conformal transformations are taken into account leading to the evident connection between the gravitational actions in the two mentioned frames and the Hilbert-Einstein action with a cosmological constant. We show that the apparent differences between Palatini and metric formalisms strictly depend on the representation while the number of degrees of freedom is preserved. This means that the dynamical content of both formalism is identical.

The equivalence between metric and Palatini formalisms in f(R)-gravity can be achieved in the general context of theories with divergence free current. This equivalence is a necessary result of a symmetry which is included in a particular conservation equation of the current. In fact the conservation equation, by an appropriate redefinition of the introduced auxiliary field, may be encoded in a massless scalar field equation.

The gravitational phase shift of neutrino oscillation can be discussed in the framework of f(R)-gravity. We show that the shift of quantum mechanical phase can depend on the given f(R)-theory that we choose. This fact is general and could constitute a fundamental test to discriminate among the various alternative relativistic theories of gravity. Estimations of ratio between the gravitational phase shift and the standard phase are carried out for the electronic Solar neutrinos.