Galileon radiation in the collapse of a thin spherical shell of matter is analyzed. In the framework of a cubic Galileon theory, we compute the field profile produced at large distances by a short collapse, finding that the radiated field has two peaks traveling ahead of light fronts. The total energy radiated during the collapse follows a power law scaling with the shell's physical width and results from two competing effects: a Vainshtein suppression of the emission and an enhancement due to the thinness of the shell.

In Ref. [1, 2] a formalism to deal with high-order perturbations of a general spherical background was developed. In this article, we apply it to the particular case of a perfect fluid background. We have expressed the perturbations of the energy-momentum tensor at any order in terms of the perturbed fluid's pressure, density and velocity. In general, these expressions are not linear and have sources depending on lower order perturbations. For the second-order case we make the explicit decomposition of these sources in tensor spherical harmonics. Then, a general procedure is given to evolve the perturbative equations of motions of the perfect fluid for any value of the harmonic label. Finally, with the problem of a spherical collapsing star in mind, we discuss the high-order perturbative matching conditions across a timelike surface, in particular the surface separating the perfect fluid interior from the exterior vacuum.

Within linearized perturbation theory, black holes decay to their final stationary state through the well-known spectrum of quasinormal modes. Here we numerically study whether nonlinearities change this picture. For that purpose we study the ringdown frequencies of gauge-invariant second-order gravitational perturbations induced by self-coupling of linearized perturbations of Schwarzschild black holes. We do so through high-accuracy simulations in the time domain of first and second-order Regge-Wheeler-Zerilli type equations, for a variety of initial data sets. We consider first-order even-parity $(\ell=2,m=\pm 2)$ perturbations and odd-parity $(\ell=2,m=0)$ ones, and all the multipoles that they generate through self-coupling. For all of them and all the initial data sets considered we find that ---in contrast to previous predictions in the literature--- the numerical decay frequencies of second-order perturbations are the same ones of linearized theory, and we explain the observed behavior. This would indicate, in particular, that when modeling or searching for ringdown gravitational waves, appropriately including the standard quasinormal modes already takes into account nonlinear effects.

We show how the interplay of non-linear dynamics, self-gravity, and fluctuations leads to self-affine behavior of matter density correlations quite generically, i.e., with a power-law exponent whose value does not depend in a very direct way on the specific details of the dynamics. This we do by means of the Renormalization Group, a powerful analytic tool for extracting asymptotic behavior of many-body systems