Numerical computation of second-order vacuum perturbations of Kerr black holes

Motivated by the desire to understand the leading-order nonlinear gravitational wave interactions around arbitrarily rapidly rotating Kerr black holes, we describe a numerical code designed to compute secondorder vacuum perturbations on such spacetimes. A general discussion of the formalism we use is presented in [N. Loutrel et al., Phys. Rev. D 103, 104017 (2021)]; here we show how we numerically implement that formalism with a particular choice of coordinates and tetrad conditions and give example results for black holes with dimensionless spin parameters a = 0.7 and a = 0.998. We first solve the Teukolsky equation for the linearly perturbedWeyl scalar Psi((1))(4), followed by direct reconstruction of the spacetime metric from Psi((1))(4), and then solve for the dynamics of the second-order perturbed Weyl scalar Psi((1))(4). This code is a first step toward a more general purpose second-order code, and we outline how our basic approach could be further developed to address current questions of interest, including extending the analysis of ringdown in black hole mergers to before the linear regime, exploring gravitational wave "turbulence" around near-extremal Kerr black holes, and studying the physics of extreme mass ratio inspiral.