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 non-trivial 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 quasi-normal 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.

Scalar hair of black holes in theories with a shift symmetry are constrained by the no-hair theorem of Hui and Nicolis, assuming spherical symmetry, time-independence of the scalar field and asymptotic flatness. The most studied counterexample is a linear coupling of the scalar with the Gauss-Bonnet invariant. However, in this case the norm of the shift-symmetry current $J^2$ diverges at the horizon casting doubts on whether the solution is physically sound. We show that this is not an issue since $J^2$ is not a scalar quantity, since $J^\mu$ is not a diff-invariant current in the presence of Gauss-Bonnet. The same theory can be written in Horndeski form with a non-analytic function $G_5 \sim \log X$. In this case the shift-symmetry current is diff-invariant, but contains powers of $X$ in the denominator, so that its divergence at the horizon is again immaterial. We confirm that other hairy solutions in the presence of non-analytic Horndeski functions are pathological, featuring divergences of physical quantities as soon as one departs from time-independence and spherical symmetry. We generalise the no-hair theorem to Beyond Horndeski and DHOST theories, showing that the coupling with Gauss-Bonnet is necessary to have hair.

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.