There are currently no reliable methods to measure transverse velocities of galaxies. This is an important piece of information that could allow us to probe the physics of structure formation as well as testing the underlying theory of gravity. The slingshot effect, a special case of the Integrated Sachs-Wolfe effect, is expected to create dipole signals in the temperature fluctuations of the Cosmic Microwave Background Radiation (CMB). This effect creates a hot spot behind and a cold spot in front of moving massive objects. The dipole signal created by the slingshot effect can be used to measure transverse velocities, but because the signal is expected to be weak, the effect has not been measured yet. The aim is to show that the slingshot effect can be measured by stacking the signals of galaxies falling into a collapsing cluster. We evaluate if the effect can probe modified gravity. We use data from a simulated galaxy catalogue (MDPL2) to mimic observations. We identify a massive galaxy cluster, and make maps of the slingshot effect around infalling galaxies. We add uncorrelated Gaussian noise to each map. The maps are rotated according to the direction to the cluster centre, such that the dipole signal will add up constructively when stacking. We compare each stack to a dipole stencil and we find the probability for a false positive in the absence of the slingshot signal. Each galaxy gives a signal of around $\Delta T/T \approx 10^{-9}$, while the precision of CMB experiments of today are $\Delta T/T \approx 4 \times 10^{-6}$. By stacking around 10 000 galaxies, the slingshot signal can be over the detectable threshold with experiments of today. However, future CMB experiments must be used to be certain of the strength of the observed signal.

The Quintessence model is one of the simplest and better known alternatives to Einstein's theory for gravity. The properties of the solutions have been studied in great detail in the background, linear and non-linear contexts in cosmology. Here we discuss new phenomenology that is induced by adding disformal terms to the interactions. Among other results, we show analytically and using cosmological simulations ran with the code \texttt{Isis} that the model posses a mechanism through which is it possible to obtain repulsive fifth forces, which are opposite to gravity. Although the equations are very complex, we also find that most of the new phenomenology can be explained by studying background quantities. We used our simulation data to test approximate relations that exist between the metric and scalar field perturbations as well as between the fifth force and gravity. Excellent agreement was found between exact and approximated solutions, which opens the way for running disformal gravity cosmological simulations using simply a Newtonian solver. These results could not only help us to find new ways of testing gravity, but also provide new motivations for building alternative models.

Extending general relativity by adding extra degrees of freedom is a popular approach for explaining the accelerated expansion of the Universe and to build high energy completions of the theory of gravity. The presence of such new degrees of freedom is, however, tightly constrained from several observations and experiments that aim to test general relativity in a wide range of scales. The viability of a given modified theory of gravity, therefore, strongly depends on the existence of a screening mechanism that suppresses the extra degrees of freedom. We perform simulations, and find that waves propagating in the new degrees of freedom can significantly impact the efficiency of some screening mechanisms, thereby threatening the viability of these modified gravity theories. Specifically, we show that the waves produced in the symmetron model can increase the amplitude of the fifth force and the parametrized post-Newtonian parameters by several orders of magnitude.

We investigate statistical properties of the distribution of matter at redshift zero in disformal gravity by using N-body simulations. The disformal model studied here consists of a conformally coupled symmetron field with an additional exponential disformal term. We conduct cosmological simulations to discover the impact of the new disformal terms in the matter power spectrum, halo mass function, and radial profile of the scalar field. We calculated the disformal geodesic equation and the equation of motion for the scalar field. We then implemented these equations into the N-body code ISIS, which is a modified gravity version of the code RAMSES. The presence of a conformal symmetron field increases both the power spectrum and mass function compared to standard gravity on small scales. Our main finding is that the newly added disformal terms tend to counteract these effects and can make the evolution slightly closer to standard gravity. We finally show that the disformal terms give rise to oscillations of the scalar field in the centre of the dark matter haloes.