We develop a full four-dimensional numerical code to study scalar gravitational radiation emitted from binary systems and probe the Vainshtein mechanism in situations that break the static and spherical symmetry, relevant for binary pulsars as well as black holes and neutron stars binaries. The present study focuses on the cubic Galileon which arises as the decoupling limit of massive theories of gravity. Limitations associated with the numerical methods prevent us from reaching a physically realistic hierarchy of scales; nevertheless, within this context we observe the same power law scaling of the radiated power as previous analytic estimates, and confirm a strong suppression of the power emitted in the monopole and dipole as compared with quadrupole radiation. Following the trend to more physically realistic parameters, we confirm the suppression of the power emitted in scalar gravitational radiation and the recovery of General Relativity with good accuracy. This paves the way for future numerical work, probing more generic, physically relevant situations and sets of interactions that may exhibit the Vainshtein mechanism.

Recently, aLIGO has announced the first direct detections of gravitational waves, a direct manifestation of the propagating degrees of freedom of gravity. The detected signals GW150914 and GW151226 have been used to examine the basic properties of these gravitational degrees of freedom, particularly setting an upper bound on their mass. It is timely to review what the mass of these gravitational degrees of freedom means from the theoretical point of view, particularly taking into account the recent developments in constructing consistent massive gravity theories. Apart from the GW150914 mass bound, a few other observational bounds have been established from the effects of the Yukawa potential, modified dispersion relation and fifth force that are all induced when the fundamental gravitational degrees of freedom are massive. We review these different mass bounds and examine how they stand in the wake of recent theoretical developments and how they compare to the bound from GW150914.

Here we consider the possibility of preheating the Universe via the parametric amplification of a massless, U(1) abelian gauge field. We assume that the gauge field is coupled to the inflaton via a conformal factor with one free parameter. We present the results of high-resolution three-dimensional simulations of this model and show this mechanism efficiently preheats the Universe to a radiation-dominated final state.

We present the simplest model for classical transitions in flux vacua. A complex field with a spontaneously broken U(1) symmetry is embedded in $M_2\times S_1$. We numerically construct different winding number vacua, the vortices interpolating between them, and simulate the collisions of these vortices. We show that classical transitions are generic at large boosts, independent of whether or not vortices miss each other in the compact $S_1$.