Growing black-hole hair in nonminimally coupled biscalar gravity

Black holes offer a unique laboratory for fundamental physics and are crucial for understanding theories beyond Einstein's theory of general relativity. In this paper, we focus on 4D effective field theories and string-theory inspired models that include scalar fields. We focus on one such model, axi-dilaton gravity, a quadratic gravity theory with two kinetically coupled scalar fields, an axion and a dilaton. To study the evolution and structure of these fields around black holes, we introduce canuda-axidil, the first open-source, parametrized numerical relativity code for quadratic and biscalar gravity. Using this code, we perform single black hole simulations to show the dynamical formation of axion and dilaton hairs and quantify the effect of higher-order terms in coupling and spin. Through these simulations, we measure the impact of black hole spin and curvature coupling strength on the profiles of the axion and dilaton and show that including kinetic coupling between the fields increases the observed deviations from general relativity. Furthermore, we simulate the axion and dilaton fields around a binary black hole coalescence demonstrating the growth of axion hair during the inspiral and the production of radiative modes for both fields.