Perturbations of spinning black holes in dynamical Chern-Simons gravity: Slow rotation equations

The detection of gravitational waves resulting from the coalescence of binary black holes by the LIGO-Virgo-Kagra Collaboration has inaugurated a new era in gravitational physics. These gravitational waves provide a unique opportunity to test Einstein's general relativity and its modifications in the regime of extreme gravity. A significant aspect of such tests involves the study of the ringdown phase of gravitational waves from binary black hole coalescence, which can be decomposed into a superposition of various quasinormal modes. In general relativity, the spectra of quasinormal modes depend on the mass, spin, and charge of the final black hole, but they can also be influenced by additional properties of the black hole spacetime, as well as corrections to the general theory of relativity. In this work, we focus on a specific modified theory known as dynamical Chern-Simons gravity. We employ the modified Teukolsky formalism developed in a previous study and lay down the foundations to investigate perturbations of slowly rotating black holes admitted by the theory. Specifically, we derive the master equations for the ψ0 and ψ4 Weyl scalar perturbations that characterize the radiative part of gravitational perturbations, as well as the master equation for the scalar field perturbations. We employ metric reconstruction techniques to obtain explicit expressions for all relevant quantities. Finally, by leveraging the properties of spin-weighted spheroidal harmonics to eliminate the angular dependence from the evolution equations, we derive two, radial, second-order, ordinary differential equations for ψ0 and ψ4, respectively. These two equations are coupled to another radial, second-order, ordinary differential equation for the scalar field perturbations. This work is the first attempt to derive a master equation for black holes in dynamical Chern-Simons gravity using curvature perturbations. The master equations we obtain can then be numerically integrated to obtain the quasinormal mode spectrum of slowly rotating black holes in this theory, making progress in the study of ringdown in dynamical Chern-Simons gravity.