In this work we study some features of head-on mergers of equilibrium solutions of the Gross-Pitaevskii-Poisson system that rules the dynamics of the ultralight bosonic dark matter model. The importance of equilibrium solutions is that they play the role of halo cores in structure formation simulations. We consider a given range of initial conditions in order to sample the parameter space in terms of mass ratio and head-on momentum. In each case we analyze the relaxation process induced by gravitational cooling in the high and low momentum regimes and estimate the relaxation time scales in each case. We detect a low frequency mode in the whole parameters space and it was found that the resulting configuration oscillates under this mode with amplitude that depends on the mass ratio and head-on momentum. In some cases the resulting configuration oscillates with changes in density of two orders of magnitude and with a matter distribution that is far from isotropic. These results could contribute to the collection of possible mass distributions considered in the reconstruction of mass profiles obtained in structure formation simulations.

We compare Einstein-Boltzmann solvers that include modifications to General Relativity and find that, for a wide range of models and parameters, they agree to a high level of precision. We look at three general purpose codes that primarily model general scalar-tensor theories, three codes that model Jordan-Brans-Dicke (JBD) gravity, a code that models f(R) gravity, a code that models covariant Galileons, a code that models Ho\v{r}ava-Lifschitz gravity and two codes that model non-local models of gravity. Comparing predictions of the angular power spectrum of the cosmic microwave background and the power spectrum of dark matter for a suite of different models, we find agreement at the sub-percent level. This means that this suite of Einstein-Boltzmann solvers is now sufficiently accurate for precision constraints on cosmological and gravitational parameters.

Light degrees of freedom that modify gravity on cosmological scales must be "screened" on solar system scales in order to be compatible with data. The Vainshtein mechanism achieves this through a breakdown of classical perturbation theory, as large interactions involving new degrees of freedom become important below the so-called Vainshtein radius. We begin to develop an extension of the Parameterized Post-Newtonian (PPN) formalism that is able to handle Vainshteinian corrections. We argue that theories with a unique Vainshtein scale must be expanded using two small parameters. In this Parameterized Post-Newtonian-Vainshteinian (PPNV) expansion, the primary expansion parameter that controls the PPN order is, as usual, the velocity $v$. The secondary expansion parameter, $\alpha$, controls the strength of the Vainshteinian correction and is a theory-specific combination of the Schwarzschild radius and the Vainshtein radius of the source that is independent of its mass. We present the general framework and apply it to Cubic Galileon theory both inside and outside the Vainshtein radius. The PPNV framework can be used to determine the compatibility of such theories with solar system and other strong-field data.