The gravitational-wave astronomical revolution began in 2015 with LIGO's observation of the coalescence of two stellar-mass black holes. Over the coming decades, ground-based detectors like LIGO will extend their reach, discovering thousands of stellar-mass binaries. In the 2030s, the space-based LISA will enable gravitational-wave observations of the massive black holes in galactic centres. Between LISA and ground-based observatories lies the unexplored decihertz gravitational-wave frequency band. Here, we propose a Decihertz Observatory to cover this band, and complement observations made by other gravitational-wave observatories. The decihertz band is uniquely suited to observation of intermediate-mass ($\sim 10^2$-$10^4 M_\odot$) black holes, which may form the missing link between stellar-mass and massive black holes, offering a unique opportunity to measure their properties. Decihertz observations will be able to detect stellar-mass binaries days to years before they merge and are observed by ground-based detectors, providing early warning of nearby binary neutron star mergers, and enabling measurements of the eccentricity of binary black holes, providing revealing insights into their formation. Observing decihertz gravitational-waves also opens the possibility of testing fundamental physics in a new laboratory, permitting unique tests of general relativity and the Standard Model of particle physics. Overall, a Decihertz Observatory will answer key questions about how black holes form and evolve across cosmic time, open new avenues for multimessenger astronomy, and advance our understanding of gravitation, particle physics and cosmology.

We investigate the ability of the Laser Interferometer Space Antenna (LISA) to measure the center of mass acceleration of stellar-origin black hole binaries emitting gravitational waves. Our analysis is based on the idea that the acceleration of the center of mass induces a time variation in the redshift of the gravitational wave, which in turn modifies its waveform. We confirm that while the cosmological acceleration is too small to leave a detectable imprint on the gravitational waveforms observable by LISA, larger peculiar accelerations may be measurable for sufficiently long lived sources. We focus on stellar mass black hole binaries, which will be detectable at low frequencies by LISA and near coalescence by ground based detectors. These sources may have large peculiar accelerations, for instance, if they form in nuclear star clusters or in AGN accretion disks. If that is the case, we find that in an astrophysical population calibrated to the LIGO-Virgo observed merger rate, LISA will be able to measure the peculiar acceleration of a small but significant fraction of the events if the mission lifetime is extended beyond the nominal duration of 4 years. In this scenario LISA will be able to assess whether black hole binaries form close to galactic centers, particularly in AGN disks, and will thus help discriminate between different formation mechanisms. Although for a nominal 4 years LISA mission the peculiar acceleration effect cannot be measured, a consistent fraction of events may be biased by strong peculiar accelerations which, if present, may imprint large systematic errors on some waveform parameters. In particular, estimates of the luminosity distance could be strongly biased and consequently induce large systematic errors on LISA measurements of the Hubble constant with stellar mass black hole binaries.

We explore the prospects for the detection of giant circumbinary exoplanets and brown dwarfs (BDs) orbiting Galactic double white dwarfs binaries (DWDs) with the Laser Interferometer Space Antenna (LISA). By assuming an occurrence rate of 50%, motivated by white dwarf pollution observations, we built a Galactic synthetic population of P-type giant exoplanets and BDs orbiting DWDs. We carried this out by injecting different sub-stellar populations, with various mass and orbital separation characteristics, into the DWD population used in the LISA mission proposal. We then performed a Fisher matrix analysis to measure how many of these three-body systems show a periodic Doppler-shifted gravitational wave perturbation detectable by LISA. We report the number of circumbinary planets (CBPs) and (BDs) that can be detected by LISA for various combinations of mass and semi-major axis distributions. We identify pessimistic and optimistic scenarios corresponding, respectively, to 3 and 83 (14 and 2218) detections of CBPs (BDs), observed during the length of the nominal LISA mission. These detections are distributed all over the Galaxy following the underlying DWD distribution, and they are biased towards DWDs with higher LISA signal-to-noise ratio and shorter orbital period. Finally, we show that if LISA were to be extended for four more years, the number of systems detected will be more than doubled in both the optimistic and pessimistic scenarios. Our results present promising prospects for the detection of post-main sequence exoplanets and BDs, showing that gravitational waves can prove the existence of these populations over the totality of the Milky Way. Detections by LISA will deepen our knowledge on the life of exoplanets subsequent to the most extreme evolution phases of their hosts, clarifying whether new phases of planetary formation take place later in the life of the stars.

Dimensional flow, the scale dependence of the dimensionality of spacetime, is a feature shared by many theories of quantum gravity (QG). We present the first study of the consequences of QG dimensional flow for the luminosity distance scaling of gravitational waves in the frequency ranges of LIGO and LISA. We find generic modifications with respect to the standard general-relativistic scaling, largely independent of specific QG proposals. We constrain these effects using two examples of multimessenger standard sirens, the binary neutron-star merger GW170817 and a simulated supermassive black-hole merger event detectable with LISA. We apply these constraints to various QG candidates, finding that the quantum geometries of group field theory, spin foams and loop quantum gravity can give rise to observable signals in the gravitational-wave spin-2 sector. Our results complement and improve GW propagation-speed bounds on modified dispersion relations. Under more model-dependent assumptions, we also show that bounds on quantum geometry can be strengthened by solar-system tests.

We investigate possible signatures of quantum gravity which could be tested with current and future gravitational-wave (GW) observations. In particular, we analyze how quantum gravity can influence the GW luminosity distance, the time dependence of the effective Planck mass and the instrumental strain noise of interferometers. Using both model-dependent and model-independent formulae, we show that these quantities can encode a non-perturbative effect typical of all quantum-gravity theories, namely the way the dimension of spacetime changes with the probed scale. Effects associated with such dimensional flow might be tested with GW observations and constrained significantly in theories with a microscopically discrete spacetime geometry, more strongly than from propagation-speed constraints. Making use of public LIGO data as well as of a simulated higher-redshift LISA source, we impose the first, respectively, actual and mock constraints on quantum-gravity parameters affecting the GW luminosity distance and discuss specific theoretical examples. If also the Newtonian potential is modified but light geodesics are not, then solar-system bounds may be stronger than GW ones. Yet, for some theories GW astronomy can give unique information not available from solar-system tests.

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 propose a space-based interferometer surveying the gravitational wave (GW) sky in the milli-Hz to $\mu$-Hz frequency range. By the 2040s', the $\mu$-Hz frequency band, bracketed in between the Laser Interferometer Space Antenna (LISA) and pulsar timing arrays, will constitute the largest gap in the coverage of the astrophysically relevant GW spectrum. Yet many outstanding questions related to astrophysics and cosmology are best answered by GW observations in this band. We show that a $\mu$-Hz GW detector will be a truly overarching observatory for the scientific community at large, greatly extending the potential of LISA. Conceived to detect massive black hole binaries from their early inspiral with high signal-to-noise ratio, and low-frequency stellar binaries in the Galaxy, this instrument will be a cornerstone for multimessenger astronomy from the solar neighbourhood to the high-redshift Universe.

Modifications of General Relativity leave their imprint both on the cosmic expansion history through a non-trivial dark energy equation of state, and on the evolution of cosmological perturbations in the scalar and in the tensor sectors. In particular, the modification in the tensor sector gives rise to a notion of gravitational-wave (GW) luminosity distance, different from the standard electromagnetic luminosity distance, that can be studied with standard sirens at GW detectors such as LISA or third-generation ground based experiments. We discuss the predictions for modified GW propagation from some of the best studied theories of modified gravity, such as Horndeski or the more general degenerate higher order scalar-tensor (DHOST) theories, non-local infrared modifications of gravity, bigravity theories and the corresponding phenomenon of GW oscillation, as well as theories with extra or varying dimensions. We show that modified GW propagation is a completely generic phenomenon in modified gravity. We then use a simple parametrization of the effect in terms of two parameters $(\Xi_0,n)$, that is shown to fit well the results from a large class of models, to study the prospects of observing modified GW propagation using supermassive black hole binaries as standard sirens with LISA. We construct mock source catalogs and perform detailed Markov Chain Monte Carlo studies of the likelihood obtained from LISA standard sirens alone, as well as by combining them with CMB, BAO and SNe data to reduce the degeneracies between cosmological parameters. We find that the combination of LISA with the other cosmological datasets allows one to measure the parameter $\Xi_0$ that characterizes modified GW propagation to the percent level accuracy, sufficient to test several modified gravity theories. [Abridged]

This Astro-2020 White Paper deals with what we might learn from future gravitational wave observations about the early universe phase transitions and their energy scales, primordial black holes, Hubble parameter, dark matter and dark energy, modified theories of gravity and extra dimensions.

There are two big questions cosmologists would like to answer -- How does the Universe work, and what are its origin and destiny? A long wavelength gravitational wave detector -- with million km interferometer arms, achievable only from space -- gives a unique opportunity to address both of these questions. A sensitive, mHz frequency observatory could use the inspiral and merger of massive black hole binaries as standard sirens, extending our ability to characterize the expansion history of the Universe from the onset of dark energy-domination out to a redshift z ~ 10. A low-frequency detector, furthermore, offers the best chance for discovery of exotic gravitational wave sources, including a primordial stochastic background, that could reveal clues to the origin of our Universe.

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.

To date more than 3500 exoplanets have been discovered orbiting a large variety of stars. Due to the sensitivity limits of the currently used detection techniques, these planets populate zones restricted either to the solar neighbourhood or towards the Galactic bulge. This selection problem prevents us from unveiling the true Galactic planetary population and is not set to change for the next two decades. Here we present a new detection method that overcomes this issue and that will allow us to detect gas giant exoplanets using gravitational wave astronomy. We show that the Laser Interferometer Space Antenna (LISA) mission can characterise hundreds of new circumbinary exoplanets orbiting white dwarf binaries everywhere in our Galaxy - a population of exoplanets so far completely unprobed - as well as detecting extragalactic bound exoplanets in the Magellanic Clouds. Such a method is not limited by stellar activity and, in extremely favourable cases, will allow LISA to detect super-Earths down to 10 Earth masses.

The Nobel Prize winning confirmation in 1998 of the accelerated expansion of our Universe put into sharp focus the need of a consistent theoretical model to explain the origin of this acceleration. As a result over the past two decades there has been a huge theoretical and observational effort into improving our understanding of the Universe. The cosmological equations describing the dynamics of a homogeneous and isotropic Universe are systems of ordinary differential equations, and one of the most elegant ways these can be investigated is by casting them into the form of dynamical systems. This allows the use of powerful analytical and numerical methods to gain a quantitative understanding of the cosmological dynamics derived by the models under study. In this review we apply these techniques to cosmology. We begin with a brief introduction to dynamical systems, fixed points, linear stability theory, Lyapunov stability, centre manifold theory and more advanced topics relating to the global structure of the solutions. Using this machinery we then analyse a large number of cosmological models and show how the stability conditions allow them to be tightly constrained and even ruled out on purely theoretical grounds. We are also able to identify those models which deserve further in depth investigation through comparison with observational data. This review is a comprehensive and detailed study of dynamical systems applications to cosmological models focusing on the late-time behaviour of our Universe, and in particular on its accelerated expansion. In self contained sections we present a large number of models ranging from canonical and non-canonical scalar fields, interacting models and non-scalar field models through to modified gravity scenarios. Selected models are discussed in detail and interpreted in the context of late-time cosmology.

A significant fraction of stars are members of gravitationally bound hierarchies containing three or more components. Almost all low mass stars in binaries with periods shorter three days are part of a hierarchical system. We therefore anticipate that a large fraction of compact galactic binaries detected by the Laser Interferometer Space Antenna (LISA) will be members of hierarchical triple or quadruple system. The acceleration imparted by the hierarchical companions can be detected in the gravitational wave signal for outer periods as large as 100 years. For systems with periods that are shorter than, or comparable to, the mission lifetime, it will be possible to measure the period and eccentricity of the outer orbit. LISA observations of hierarchical stellar systems will provide insight into stellar evolution, including the role that Kozai-Lidov oscillations play in driving systems towards merger.

Multi-frequency gravitational wave (GW) observations are useful probes of the formation processes of coalescing stellar-mass binary black holes (BBHs). We discuss the phase drift in the GW inspiral waveform of the merging BBH caused by its center-of-mass acceleration. The acceleration strongly depends on the location where a BBH forms within a galaxy, allowing observations of the early inspiral phase of LIGO-like BBH mergers by the Laser Interferometer Space Antenna (LISA) to test the formation mechanism. In particular, BBHs formed in dense nuclear star clusters or via compact accretion disks around a nuclear supermassive black hole in active galactic nuclei would suffer strong acceleration, and produce large phase drifts measurable by LISA. The host galaxies of the coalescing BBHs in these scenarios can also be uniquely identified in the LISA error volume, without electromagnetic counterparts. A non-detection of phase drifts would rule out or constrain the contribution of the nuclear formation channels to the stellar-mass BBH population.

We consider scalar field models of dark energy interacting with dark matter through a coupling proportional to the contraction of the four-derivative of the scalar field with the four-velocity of the dark matter fluid. The coupling is realized at the Lagrangian level employing the formalism of Scalar-Fluid theories, which use a consistent Lagrangian approach for relativistic fluid to describe dark matter. This framework produces fully covariant field equations, from which we can derive unequivocal cosmological equations at both background and linear perturbations levels. The background evolution is analyzed in detail applying dynamical systems techniques, which allow us to find the complete asymptotic behavior of the universe given any set of model parameters and initial conditions. Furthermore we study linear cosmological perturbations investigating the growth of cosmic structures within the quasi-static approximation. We find that these interacting dark energy models give rise to interesting phenomenological dynamics, including late-time transitions from dark matter to dark energy domination, matter and accelerated scaling solutions and dynamical crossing of the phantom barrier. Moreover we obtain possible deviations from standard $\Lambda$CDM behavior at the linear perturbations level, which have an impact on the dynamics of structure formation and might provide characteristic observational signatures.

We perform a forecast analysis of the ability of the LISA space-based interferometer to reconstruct the dark sector interaction using gravitational wave (GW) standard sirens at high redshift. We employ Gaussian process methods to reconstruct the distance-redshift relation in a model independent way. We adopt simulated catalogues of standard sirens given by merging massive black hole binaries (MBHBs) visible by LISA, with an electromagnetic (EM) counterpart detectable by future telescopes. The catalogues are constructed considering three different astrophysical scenarios for the evolution of MBHB mergers based on the semi-analytic model of E. Barausse (2012). We first use these standard siren datasets to assess the potential of LISA in reconstructing a possible interaction between vacuum dark energy and dark matter. Then we combine the LISA cosmological data with supernovae data simulated for the Dark Energy Survey (DES). We consider two scenarios distinguished by the time duration of the LISA mission: 5 and 10 years. Using only LISA standard siren data, the dark sector interaction can be well reconstructed from redshift $z\sim1$ to $z\sim3$ (5 yr) and $z\sim1$ up to $z\sim5$ (10 yr), though the reconstruction is inefficient at lower redshift. When combined with the DES datasets, the interaction is well reconstructed in the whole redshift region from $z\sim0$ to $z\sim3$ (5 yr) and $z\sim0$ to $z\sim5$ (10 yr). MBHB standard sirens can thus be used to constrain the dark sector interaction at redshift ranges not reachable by usual supernovae datasets which probe only the $z\lesssim 1.5$ range. GW standard sirens will not only constitute a complementary and alternative way, with respect to familiar EM observations, to probe the cosmic expansion, but will also provide new tests to constrain possible deviations from the standard $\Lambda$CDM dynamics, especially at high redshift.

We expand the dynamical systems investigation of cosmological scalar fields characterised by kinetic corrections presented in [N. Tamanini, Phys. Rev. D 89 (2014) 083521]. In particular we do not restrict the analysis to exponential potentials only, but we consider arbitrary scalar field potentials and derive general results regarding the corresponding cosmological dynamics. Two specific potentials are then used as examples to show how these models can be employed not only to describe dark energy, but also to achieve dynamical crossing of the phantom barrier at late times. Stability and viability issues at the classical level are also discussed.

We perform a forecast analysis of the capability of the eLISA space-based interferometer to constrain models of early and interacting dark energy using gravitational wave standard sirens. We employ simulated catalogues of standard sirens given by merging massive black hole binaries visible by eLISA, with an electromagnetic counterpart detectable by future telescopes. We consider three-arms mission designs with arm length of 1, 2 and 5 million km, 5 years of mission duration and the best-level low frequency noise as recently tested by the LISA Pathfinder. Standard sirens with eLISA give access to an intermediate range of redshift $1\lesssim z \lesssim 8$, and can therefore provide competitive constraints on models where the onset of the deviation from $\Lambda$CDM (i.e. the epoch when early dark energy starts to be non-negligible, or when the interaction with dark matter begins) occurs relatively late, at $z\lesssim 6$. If instead early or interacting dark energy is relevant already in the pre-recombination era, current cosmological probes (especially the cosmic microwave background) are more efficient than eLISA in constraining these models, except possibly in the interacting dark energy model if the energy exchange is proportional to the energy density of dark energy.

We investigate the cosmological dynamics of the recently proposed extended chameleon models at both background and linear perturbation levels. Dynamical systems techniques are employed to fully characterize the evolution of the universe at the largest distances, while structure formation is analysed at sub-horizon scales within the quasi-static approximation. The late time dynamical transition from dark matter to dark energy domination can be well described by almost all extended chameleon models considered, with no deviations from $\Lambda$CDM results at both background and perturbation levels. The results obtained in this work confirm the cosmological viability of extended chameleons as alternative dark energy models.

We extend the chameleon models by considering Scalar-Fluid theories where the coupling between matter and the scalar field can be represented by a quadratic effective potential with density-dependent minimum and mass. In this context, we study the effects of the scalar field on Solar System tests of gravity and show that models passing these stringent constraints can still induce large modifications of Newton's law on galactic scales. On these scales we analyse models which could lead to a percent deviation of Newton's law outside the virial radius. We then model the dark matter halo as a Navarro-Frenk-White profile and explicitly find that the fifth force can give large contributions around the galactic core in a particular model where the scalar field mass is constant and the minimum of its potential varies linearly with the matter density. At cosmological distances, we find that this model does not alter the growth of large scale structures and therefore would be best tested on galactic scales, where interesting signatures might arise in the galaxy rotation curves.

We investigate the capability of various configurations of the space interferometer eLISA to probe the late-time background expansion of the universe using gravitational wave standard sirens. We simulate catalogues of standard sirens composed by massive black hole binaries whose gravitational radiation is detectable by eLISA, and which are likely to produce an electromagnetic counterpart observable by future surveys. The main issue for the identification of a counterpart resides in the capability of obtaining an accurate enough sky localisation with eLISA. This seriously challenges the capability of four-link (2 arm) configurations to successfully constrain the cosmological parameters. Conversely, six-link (3 arm) configurations have the potential to provide a test of the expansion of the universe up to $z\sim 8$ which is complementary to other cosmological probes based on electromagnetic observations only. In particular, in the most favourable scenarios, they can provide a significant constraint on $H_0$ at the level of 0.5%. Furthermore, $(\Omega_M, \Omega_\Lambda)$ can be constrained to a level competitive with present SNIa results. On the other hand, the lack of massive black hole binary standard sirens at low redshift allows to constrain dark energy only at the level of few percent.

Following the selection of The Gravitational Universe by ESA, and the successful flight of LISA Pathfinder, the LISA Consortium now proposes a 4 year mission in response to ESA's call for missions for L3. The observatory will be based on three arms with six active laser links, between three identical spacecraft in a triangular formation separated by 2.5 million km. LISA is an all-sky monitor and will offer a wide view of a dynamic cosmos using Gravitational Waves as new and unique messengers to unveil The Gravitational Universe. It provides the closest ever view of the infant Universe at TeV energy scales, has known sources in the form of verification binaries in the Milky Way, and can probe the entire Universe, from its smallest scales near the horizons of black holes, all the way to cosmological scales. The LISA mission will scan the entire sky as it follows behind the Earth in its orbit, obtaining both polarisations of the Gravitational Waves simultaneously, and will measure source parameters with astrophysically relevant sensitivity in a band from below $10^{-4}\,$Hz to above $10^{-1}\,$Hz.

Third generation ground-based interferometers as well as the planned space-based interferometer LISA are expected to detect a plethora of gravitational wave signals from coalescing binaries at cosmological distance. The emitted gravitational waves propagate in the expanding universe through the inhomogeneous distribution of matter. Here we show that the acceleration of the universe and the peculiar acceleration of the binary with respect to the observer distort the gravitational chirp signal from the simplest General Relativity prediction beyond a mere time independent rescaling of the chirp mass, affecting intrinsic parameter estimations for the binaries visible by LISA. We find that the effect due to the peculiar acceleration can be much larger than the one due to the universe acceleration. Moreover, peculiar accelerations can introduce a bias in the estimation of parameters such as the time of coalescence and the individual masses of the binary. An error in the estimation of the arrival time will have an impact in the case of sources visible first by LISA and later by ground based interferometers.

We extend the dynamical systems analysis of Scalar-Fluid interacting dark energy models performed in C. G. Boehmer et al, Phys. Rev. D 91, 123002 (2015), by considering scalar field potentials beyond the exponential type. The properties and stability of critical points are examined using a combination of linear analysis, computational methods and advanced mathematical techniques, such as centre manifold theory. We show that the interesting results obtained with an exponential potential can generally be recovered also for more complicated scalar field potentials. In particular, employing power-law and hyperbolic potentials as examples, we find late time accelerated attractors, transitions from dark matter to dark energy domination with specific distinguishing features, and accelerated scaling solutions capable of solving the cosmic coincidence problem.

This paper summarises the potential of the LISA mission to constrain the expansion history of the universe using massive black hole binary mergers as gravitational wave standard sirens. After briefly reviewing the concept of standard siren, the analysis and methodologies of Ref. [1] are briefly outlined to show how LISA can be used as a cosmological probe, while a selection of results taken from Refs. [1,2] is presented in order to estimate the power of LISA in constraining cosmological parameters.

A new Lagrangian framework has recently been proposed to describe interactions between relativistic perfect fluids and scalar fields. In this paper we investigate the Einstein static universe in this new class of theories, which have been named Scalar-Fluid theories. The stability of the static solutions to both homogeneous and inhomogeneous perturbations is analysed deriving the relevant cosmological perturbation equations at the linear order. We can find several configurations corresponding to an Einstein static universes which are stable against inhomogeneous perturbations, but unstable against homogeneous perturbations. This shows the possible applications of Scalar-Fluid theories to the inflationary emergent universe scenario.

An interaction between dark matter and dark energy is usually introduced by a phenomenological modification of the matter conservation equations, while the Einstein equations are left unchanged. Starting from some general and fundamental considerations, in this work it is shown that a coupling in the dark sector is likely to introduce new terms also in the gravitational dynamics. Specifically in the cosmological background equations a bulk dissipative pressure, characterizing viscous effects and able to suppress structure formation at small scales, should appear from the dark coupling. At the level of the perturbations the analysis presented in this work reveals instead the difficulties in properly defining the dark sector interaction from a phenomenological perspective.

We present a new approach to build models of quintessence interacting with dark or baryonic matter. We use a variational approach for relativistic fluids to realize an effective description of matter fields at the Lagrangian level. The coupling is introduced directly in the action by considering a single function mixing the dynamical degrees of freedom of the theory. The resulting gravitational field equations are derived by variations with respect to the independent variables. New interesting phenomenology can be obtained at both small scales, where new screening mechanisms for scalar fields can be realized, and large scales, where one finds an original and rich class of interacting quintessence models. The background cosmology of two of these models is studied in detail using dynamical system techniques. We find a variety of interesting results: for instance, these models contain dark energy dominated late time attractors and scaling solutions, both with early time matter dominated epochs and a possible inflationary origin. In general this new approach provides the starting point for future in depth studies on new interacting quintessence models.

We consider an original variational approach for building new models of quintessence interacting with dark or baryonic matter. The coupling is introduced at the Lagrangian level using a variational formulation for relativistic fluids, where the interacting term generally depends on both the dynamical degrees of freedom of the theory and their spacetime derivatives. After deriving the field equations from the action, we consider applications in the context of cosmology. Two simple models are studied using dynamical system techniques showing the interesting phenomenology arising in this framework. We find that these models contain dark energy dominated late time attractors with early time matter dominated epochs and also obtain a possible dynamical crossing of the phantom barrier. The formulation and results presented here complete and expand the analysis exposed in the first part of this work, where only algebraic couplings, without spacetime derivatives, were considered.

Recently a new Lagrangian framework was introduced to describe interactions between scalar fields and relativistic perfect fluids. This allows two consistent generalizations of coupled quintessence models: non-vanishing pressures and a new type of derivative interaction. Here the implications of these to the formation of cosmological large-scale structure are uncovered at the linear order. The full perturbation equations in the two cases are derived in a unified formalism and their Newtonian, quasi-static limit is studied analytically. Requiring the absence of an effective sound speed for the coupled dark matter fluid restricts the Lagrangian to be a linear function of the matter number density. This still leaves new potentially viable classes of both algebraically and derivatively interacting models wherein the coupling may impact the background expansion dynamics and imprint signatures into the large-scale structure.

We investigate the cosmology of interacting spin-2 particles, formulating the multi-gravitational theory in terms of vierbeins and without imposing any Deser-van Nieuwenhuizen-like constraint. The resulting multi-vierbein theory represents a wider class of gravitational theories if compared to the corresponding multi-metric models. Moreover, as opposed to its metric counterpart which in general seems to contain ghosts, it has already been proved to be ghost-free. We outline a discussion about the possible matter couplings and we focus on the study of cosmological scenarios in the case of three and four interacting vierbeins. We find rich behavior, including de Sitter solutions with an effective cosmological constant arising from the multi-vierbein interaction, dark-energy solutions and nonsingular bouncing solutions.

The background dynamical evolution of a universe filled with matter and a cosmological scalar field is analyzed employing dynamical system techniques. After the phenomenology of a canonical scalar field with exponential potential is revised, square and square root kinetic corrections to the scalar field canonical Lagrangian are considered and the resulting dynamics at cosmological distances is obtained and studied. These noncanonical cosmological models imply new interesting phenomenology including early time matter dominated solutions, cosmological scaling solutions and late time phantom dominated solutions with dynamical crossing of the phantom barrier. Stability and viability issues for these scalar fields are presented and discussed.

We propose a new point of view for interpreting Newton's and Einstein's theories of gravity. By taking inspiration from Continuum Mechanics and its treatment of anisotropies, we formulate new gravitational actions for modified theories of gravity. These models are simple and natural generalisations with many interesting properties. Above all, their precise form can, in principle, be determined experimentally.

Theories with a non-minimal coupling between the space-time curvature and matter fields introduce an extra force due to the non-conservation of the matter energy momentum. In the present work the theoretical consistency of such couplings is studied using a scalar field Lagrangian to model the matter content. The conditions that the coupling does not introduce ghosts, classical instabilities or superluminal propagation of perturbations are derived. These consistency conditions are then employed to rule out or severely restrict the forms of the non-minimal coupling functions considered in the previous literature. For example, a power-law coupling is viable only for sublinear positive power of the curvature scalar.

We examine the existence of relativistic stars in f(T) modified gravity and explicitly construct several classes of static perfect fluid solutions. We derive the conservation equation from the complete f(T) gravity field equations and present the differences with its teleparallel counterpart. Firstly, we choose the tetrad field in the diagonal gauge and study the resulting field equations. Some exact solutions are explicitly constructed and it is noted that these solutions have to give a constant torsion scalar. Next, we choose a non diagonal tetrad field which results in field equations similar to those of general relativity. For specific models we are able to construct exact solutions of these field equations. Among those new classes of solutions, we find negative pressure solutions, and an interesting class of polynomial solutions.

We examine the generation of primordial perturbations during an inflationary epoch in generalised theories of gravity when the equations of motion are derived using the Palatini variational principle. Both f(R) and Scalar-Tensor theories are considered and we compare our results with those obtained under the conventional metric formalism. Non-linear generalisations of the action lead to different theories under the two variational choices and we obtain distinct results for scalar and tensor spectral indices and their ratio. We find the following general result; inflation driven solely by f(R) modifications alone do not result in suitable curvature perturbations whilst Scalar-Tensor theories generate nearly scalar invariant curvature perturbations but no tensor modes.