We present a method to introduce relativistic corrections including linear dark energy perturbations in Horndeski theory into Newtonian simulations based on the N-body gauge approach. We assume that standard matter species (cold dark matter, baryons, photons and neutrinos) are only gravitationally-coupled with the scalar field and we then use the fact that one can include modified gravity effects as an effective dark energy fluid in the total energy-momentum tensor. In order to compute the scalar field perturbations, as well as the cosmological background and metric perturbations, we use the Einstein-Boltzmann code \hiclass. As an example, we study the impact of relativistic corrections on the matter power spectrum in k-essence, a subclass of Horndeski theory, including the effects of massless and massive neutrinos. For massive neutrinos with $\sum m_{\nu} = 0.1$ eV, the corrections due to relativistic species (photons, neutrinos and dark energy) can introduce a maximum deviation of approximately $8\%$ to the power spectrum at $k \sim 10^{-3} \ \textrm{Mpc}^{-1}$ at $z=0$, for a scalar field with sound speed $c_{s}^{2}\sim0.013$ during matter domination epoch. Our formalism makes it possible to test beyond \lcdm \ models probed by upcoming large-scale structure surveys on very large scales.

Perturbations from inflation evolve into large scale structure of the late universe, and encode abundant cosmic structure formation physics. We allow freedom in the primordial power spectrum, rather than assuming a power law scale dependence, to study its impact on cosmological parameter determination. Combining various generations of cosmic microwave background (CMB) data and galaxy redshift survey data, we investigate the constraints on reconstruction of the primordial curvature perturbation power spectrum and the late time cosmology, especially the sum of neutrino masses. We quantify how each successive generation, in CMB and galaxy surveys, provides significant improvements, often by factors of several. By using CMB polarization information over a broad range of angular scales, and galaxy redshift data in many bins of redshift, one can allow inflationary freedom and still constrain parameters comparably to assuming power law dependence. The primordial power spectrum can be reconstructed at the subpercent level in a dozen wavenumber bins, while simultaneously fitting the sum of neutrino masses to 14 meV.

Within the effective field theory approach to cosmic acceleration, the background expansion can be specified separately from the gravitational modifications. We explore the impact of modified gravity in a background different from a cosmological constant plus cold dark matter ($\Lambda$CDM) on the stability and cosmological observables, including covariance between gravity and expansion parameters. In No Slip Gravity the more general background allows more gravitational freedom, including both positive and negative Planck mass running. We examine the effects on cosmic structure growth, as well as showing that a viable positive integrated Sachs-Wolfe effect crosscorrelation easily arises from this modified gravity theory. Using current data we constrain parameters with a Monte Carlo analysis, finding a maximum running $|\alpha_M|\lesssim 0.03$. We provide the modified {\tt hi\_class} code publicly on GitHub, now enabling computation and inclusion of the redshift space distortion observable $f\sigma_8$ as well as the No Slip Gravity modifications.

Generally the Brans-Dicke theory reduces to General Relativity in the limit $\omega\rightarrow\infty$ if the scalar field goes as $\phi\propto1/\omega$. However, it is also known that there are examples with $\phi\propto1/\sqrt{\omega}$ that does not tend to GR. We discuss another case: a homogeneous and isotropic universe filled with stiff matter. The power of time dependence of these solutions do not depend on $\omega$, and there is no General Relativity limit even though we have $\phi\propto1/\omega$. A perturbative analysis of this exotic case is carried out.

We study the formation of classical singularities in Generalized Brans-Dicke theories that are natural extensions to Brans-Dicke where the kinetic term is modified by a new coupling function $\omega(\varphi)$. We discuss the asymptotic limit $\omega(\varphi)\rightarrow\infty$ and show that the system generically does not approach General Relativity. Given the arbitrariness of $\omega(\varphi)$, one can search for coupling functions chosen specifically to avoid classical singularities. However, we prove that this is not the case. Homogeneous and spherically symmetric collapsing objects form singularities for arbitrary coupling functions. On the other hand, expanding cosmological scenarios are completely free of Big Rip type singularities. In an expanding universe, the scalar field behaves at most as stiff matter, which makes these cosmological solutions asymptotically approach General Relativity.