Accurate model for recoil velocity in binary black hole mergers up to mass ratio 200:1

Using a combination of numerical relativity and black hole perturbation theory simulations for binary black hole mergers spanning mass ratios from 1:1 up to 200:1, we present an accurate model for the recoil velocity of the remnant black hole. Our approach integrates post-Newtonian–based analytical prescriptions with data-driven Gaussian process regression to enforce the correct theoretical limiting behaviors near the equal-mass and extreme–mass-ratio regimes while maintaining high accuracy across the parameter space. We make our model, gwModel_kick_q200, publicly available through the gwModels package.