It is usually assumed that a healthy EFT should not allow superluminal propagation. In the presence of gravity, however, the notion of superluminality becomes subtle, since there is no invariant way to compare with an underlying Minkowski light cone. One can instead resort to an asymptotic criterion: whether the EFT can induce signal propagation faster than what allowed by the asymptotic structure of spacetime. In this work we study the asymptotic causal structure of gravitational EFTs by analysing signal propagation in black-hole backgrounds in the presence of higher-derivative operators. We show that in spacetime dimensions D>4 the effective light cones can lead to genuine asymptotic superluminality, which can be used to constrain the regime of validity of the EFT. By contrast, in D=4 the asymptotic causal structure is universally identical to that of Schwarzschild: prompt null curves remain insensitive to higher-derivative corrections and no asymptotic time advance is possible. We first study the representative operator R_μνρσF^μνF^ρσ, then show that this conclusion is true for any EFT, as it relies only on the asymptotic behaviour of the metric. Finally, we discuss two ways to define superluminality in D=4 spacetimes: introducing a covariant cut-off by putting the theory in an asymptotically-AdS background, or imposing a hard cut-off by working at finite distance.

Quadratic quasi-normal modes, generated at second order in black hole perturbation theory, are a promising target for testing gravity in the nonlinear regime with next-generation gravitational wave detectors. While their frequencies have long been known, their amplitudes remain poorly studied. We introduce regular variables and compute amplitudes for Schwarzschild black holes with the Leaver algorithm. We find a nonlinear ratio ℛ≃0.154e^-0.068i for the most excited ℓ=4 mode, matching results from Numerical Relativity. We also predict new low-frequency ℓ=2 quadratic modes.

General Relativity predicts the existence of quadratic quasi-normal modes at second order in perturbation theory. Building on our recent work, we compute the amplitudes and polarizations of these modes for non-rotating black holes, showing that they are completely determined by the amplitudes and polarizations of linear modes. We obtain the ratio of quadratic to linear amplitudes, which still depends on the initial conditions of the merger through the polarization of linear modes. However, we demonstrate that this dependence is captured by four fundamental numbers, independent of initial conditions, representing four different combinations of linear modes parities. Additionally, we prove two selection rules regarding the vanishing of classes of quadratic modes. Our results are available online as a package which provides the ratio of amplitudes across a broad spectrum of angular momenta.

Quadratic quasinormal mode, generated at second order in black hole perturbation theory, are a promising target for testing gravity in the nonlinear regime with next-generation gravitational wave detectors. While their frequencies have long been known, their amplitudes remain poorly studied. We introduce regular variables and compute amplitudes for Schwarzschild black holes with the Leaver algorithm. We find a nonlinear ratio R similar or equal to 0.154 e - 0.068 i for the most excited ( 1 / 4 4 mode, matching results from numerical relativity. We also predict new low-frequency ( 1 / 4 2 quadratic modes.

Studies of black hole superradiance often focus on the growth of a cloud in isolation, accompanied by the spin-down of the black hole. In this paper, we consider the additional effect of the accretion of matter and angular momentum from the environment. We show that, in many cases, the black hole evolves by drifting along the superradiance threshold, in which case the evolution of its parameters can be described analytically or semi-analytically. We quantify the conditions under which accretion can serve as a mechanism to increase the cloud-to-black hole mass ratio, beyond the standard maximum of about 10%. This occurs by a process we call over-superradiance, whereby accretion effectively feeds the superradiance cloud, by way of the black hole. We give two explicit examples: accretion from a vortex expected in wave dark matter and accretion from a baryonic disk. In the former case, we estimate the accretion rate by using an analytical fit to the asymptotic behavior of the confluent Heun function. Level transition, whereby one cloud level grows while the other shrinks, can be understood in a similar way.

Many binary systems of interest for gravitational-wave astronomy are orbited by a third distant body, which can considerably alter their relativistic dynamics. Precision computations are needed to understand the interplay between relativistic corrections and three-body interactions. We use an effective field theory approach to derive the effective action describing the long time-scale dynamics of hierarchical three-body systems up to 1PN quadrupole order. At this level of approximation, computations are complicated by the backreaction of small oscillations on orbital time-scales as well as deviations from the adiabatic approximation. We address these difficulties by eliminating the fast modes through the method of near-identity transformations. This allows us to compute for the first time the complete expression of the 1PN quadrupole cross-terms in generic configurations of three-body systems. We numerically integrate the resulting equations of motion and show that 1PN quadrupole terms can affect the long term dynamics of relativistic three-body systems.

Dark matter (DM) can consist of a scalar field so light that DM particles in the galactic halo are best described by classical waves. We investigate how these classical solutions are influenced by the presence of a non-rotating supermassive black hole at the center of the galaxy, using an analytical, albeit approximate, approach. Relying on this analytic control, we examine the consequences of imposing causal boundary conditions at the horizon, which are typically overlooked. First, we examine the scenario where the backreaction of dark matter can be neglected. The scalar field decays like a power law at large distances, thus endowing the black hole with "hair". We derive solutions for the field profile over a wide range of parameters, including cases with rotating dark matter. As a by-product, we extract the dynamical Love numbers for scalar perturbations. Next, we determine the spectrum of bound states and their behaviour. Finally, we incorporate the self-gravity of the scalar field, with a focus on the situation where dark matter forms a soliton (boson star) at the center of the galaxy. We derive an analytical expression for the soliton at every distance from the center. With a solution that remains applicable even at horizon scales, we can reliably compute the accretion rate of the black hole.

Abstract Quasinormal modes describe the ringdown of compact objects deformed by small perturbations. In generic theories of gravity that extend General Relativity, the linearized dynamics of these perturbations is described by a system of coupled linear differential equations of second order. We first show, under general assumptions, that such a system can be brought to a Schrödinger-like form. We then devise an analytic approximation scheme to compute the spectrum of quasinormal modes. We validate our approach using a toy model with a controllable mixing parameter ε and showing that the analytic approximation for the fundamental mode agrees with the numerical computation when the approximation is justified. The accuracy of the analytic approximation is at the (sub-) percent level for the real part and at the level of a few percent for the imaginary part, even when ε is of order one. Our approximation scheme can be seen as an extension of the approach of Schutz and Will [1] to the case of coupled systems of equations, although our approach is not phrased in terms of a WKB analysis, and offers a new viewpoint even in the case of a single equation.

Recent works have suggested that nonlinear (quadratic) effects in black hole perturbation theory may be important for describing a black hole ringdown. We show that the technique of uniform approximations can be used to accurately compute 1) nonlinear amplitudes at large distances in terms of the linear ones, 2) linear (and nonlinear) quasi-normal mode frequencies, 3) the wavefunction for both linear and nonlinear modes. Our method can be seen as a generalization of the WKB approximation, with the advantages of not losing accuracy at large overtone number and not requiring matching conditions. To illustrate the effectiveness of this method we consider a simplified source for the second-order Zerilli equation, which we use to numerically compute the amplitude of nonlinear modes for a range of values of the angular momentum number.

Abstract We derive causality constraints on the simplest scalar-tensor theories in which black holes differ from what General Relativity predicts, a scalar coupled to the Gauss-Bonnet or the Chern-Simons terms. Demanding that time advances are unobservable within the regime of validity of these effective field theories, we find their cutoff must be parametrically of the same size as the inverse Schwarzschild radius of the black holes for which the non-standard effects are of order one. For astrophysical black holes within the range of current gravitational wave detectors, this means a cutoff length of the order of kilometers. We further explore the leading additional higher-dimensional operators potentially associated with the scale of UV completion and discuss their phenomenological implications for gravitational wave science.

Scalar-tensor theories are frequently only consistent with fifth force constraints in the presence of a screening mechanism, namely in order to suppress an otherwise unacceptably large coupling between the scalar and ordinary matter. Here we investigate precisely which subsets of Horndeski theories do not give rise to and/or require such a screening mechanism. We investigate these subsets in detail, deriving their form and discussing how they are restricted upon imposing additional bounds from the speed of gravitational waves, solar system tests and cosmological observables. Finally, we also identify what subsets of scalar-tensor theories precisely recover the predictions of standard (linearised) ΛCDM cosmologies in the quasi-static limit.

The motion of three bodies can be solved perturbatively when a tightly bound inner binary is orbited by a distant perturber, giving rise for example to the well-known Kozai-Lidov oscillations. We propose to study the relativistic hierarchical three-body orbits by adapting the Effective Field Theory techniques used in the two-body problem. This allows us to conveniently treat the inner binary as an effective point-particle, thus reducing the complexity of the three-body problem to a simpler spinning two-body motion. We present in details the mapping between the inner binary osculating elements and the resulting spin of the effective point-particle. Our study builds towards a derivation of three-body analytic waveforms.

A bstract We develop the effective theory for perturbations around black holes with scalar hair, in two directions. First, we show that the scalar-Gauss-Bonnet theory, often used as an example exhibiting scalar black hole hair, can be deformed by galileon operators leading to order unity changes to its predictions. The effective theory for perturbations thus provides an efficient framework for describing and constraining broad classes of scalar-tensor theories, of which the addition of galileon operators is an example. Second, we extend the effective theory to perturbations around an axisymmetric, slowly rotating black hole, at linear order in the black hole spin. We also discuss the inclusion of parity-breaking operators in the effective theory.

Under the assumption that a dynamical scalar field is responsible for the current acceleration of the Universe, we explore the possibility of probing its physics in black hole merger processes with gravitational wave interferometers. Remaining agnostic about the microscopic physics, we use an effective field theory approach to describe the scalar dynamics. We investigate the case in which some of the higher-derivative operators, that are highly suppressed on cosmological scales, instead become important on typical distances for black holes. If a coupling to the Gauss-Bonnet operator is one of them, a nontrivial background profile for the scalar field can be sourced in the surroundings of the black hole, resulting in a potentially large amount of "hair." In turn, this can induce sizeable modifications to the spacetime geometry or a mixing between the scalar and the gravitational perturbations. Both effects will ultimately translate into a modification of the quasinormal mode spectrum in a way that is also sensitive to other operators besides the one sourcing the scalar background. The presence of deviations from the predictions of general relativity in the observed spectrum can therefore serve as a window onto dark energy physics.

We reconsider the issue of whether scalar-tensor theories can admit stable wormhole configurations supported by a non-trivial radial profile for the scalar field. Using a recently proposed effective theory for perturbations around static, spherically symmetric backgrounds, we show that scalar-tensor theories of “beyond Horndeski” type can have wormhole solutions that are free of ghost and gradient instabilities. Such solutions are instead forbidden within the more restrictive “Horndeski” class of theories.

The final ringdown phase in a coalescence process is a valuable laboratory to test General Relativity and potentially constrain additional degrees of freedom in the gravitational sector. We introduce here an effective description for perturbations around spherically symmetric spacetimes in the context of scalar-tensor theories, which we apply to study quasi-normal modes for black holes with scalar hair. We derive the equations of motion governing the dynamics of both the polar and the axial modes in terms of the coefficients of the effective theory. Assuming the deviation of the background from Schwarzschild is small, we use the WKB method to introduce the notion of "light ring expansion". This approximation is analogous to the slow-roll expansion used for inflation, and it allows us to express the quasinormal mode spectrum in terms of a small number of parameters. This work is a first step in describing, in a model independent way, how the scalar hair can affect the ringdown stage and leave signatures on the emitted gravitational wave signal. Potential signatures include the shifting of the quasi-normal spectrum, the breaking of isospectrality between polar and axial modes, and the existence of scalar radiation.

Higher derivative scalar interactions can give rise to interesting cosmological scenarios. We present a complete classification of such operators that can yield sizeable effects without introducing ghosts and, at the same time, define an effective field theory robust under the inclusion of quantum corrections. A set of rules to power count consistently the coefficients of the resulting Lagrangian is provided by the presence of an approximate global symmetry. The interactions that we derive in this way contain a subset of the so-called Horndeski and beyond Horndeski theories. Our construction therefore provides a structurally robust context to study their phenomenology. Applications to dark energy/modified gravity and geodesically complete cosmologies are briefly discussed.

We study the stability of spatially flat FRW solutions which are geodesically complete, i.e. for which one can follow null (graviton) geodesics both in the past and in the future without ever encountering singularities. This is the case of NEC-violating cosmologies such as smooth bounces or solutions which approach Minkowski in the past. We study the EFT of linear perturbations around a solution of this kind, including the possibility of multiple fields and fluids. One generally faces a gradient instability which can be avoided only if the operator (3) R δ N is present and its coefficient changes sign along the evolution. This operator (typical of beyond-Horndeski theories) does not lead to extra degrees of freedom, but cannot arise starting from any theory with second-order equations of motion. The change of sign of this operator prevents to set it to zero with a generalised disformal transformation.

Many alternatives to canonical slow-roll inflation have been proposed over the years, one of the main motivations being to have a model, capable of generating observable values of non-Gaussianity. In this work, we (re-)explore the physical implications of a great majority of such models within a single, effective field theory framework (including novel models with large non-Gaussianity discussed for the first time below). The constraints we apply—both theoretical and experimental—are found to be rather robust, determined to a great extent by just three parameters: the coefficients of the quadratic EFT operators (δN)2 and δNδE, and the slow-roll parameter ε. This allows to significantly limit the majority of single-field alternatives to canonical slow-roll inflation. While the existing data still leaves some room for most of the considered models, the situation would change dramatically if the current upper limit on the tensor-to-scalar ratio decreased down to r < 10−2. Apart from inflationary models driven by plateau-like potentials, the single-field model that would have a chance of surviving this bound is the recently proposed slow-roll inflation with weakly-broken galileon symmetry. In contrast to canonical slow-roll inflation, the latter model can support r < 10−2 even if driven by a convex potential, as well as generate observable values for the amplitude of non-Gaussianity.

Canonical models of single-field, slow-roll inflation do not lead to appreciable non-Gaussianity, unless derivative interactions of the inflaton become uncontrollably large. We propose a novel slow-roll scenario where scalar perturbations propagate at a subluminal speed, leading to sizeable equilateral non-Gaussianity, f NL equil  ∝ 1/c s 4 , largely insensitive to the ultraviolet physics. The model is based on a low-energy effective theory characterized by weakly broken invariance under internal galileon transformations, ϕ → ϕ + b μ x μ , which protects the properties of perturbations from large quantum corrections. This provides the unique alternative to models such as DBI inflation in generating strongly subluminal/non-Gaussian scalar perturbations.