Stepping up superradiance constraints on axions

Light feebly-coupled bosonic particles can efficiently extract the rotational energy of rapidly spinning black holes on subastrophysical timescales via a phenomenon known as black hole superradiance. In the case of light axions, the feeble self-interactions of these particles can lead to a nonlinear coupled evolution of many superradiant quasibound states, dramatically altering the rate at which the black hole is spun down. In this work, we extend the study of axion superradiance to higher order states, solving for the first time the coupled evolution of all states with <a:math xmlns:a="http://www.w3.org/1998/Math/MathML" display="inline"><a:mi>n</a:mi><a:mo>≤</a:mo><a:mn>5</a:mn></a:math> in the fully relativistic limit (with <c:math xmlns:c="http://www.w3.org/1998/Math/MathML" display="inline"><c:mi>n</c:mi></c:math> being the principal quantum number). Using a Bayesian framework, we rederive constraints on axions using the inferred spins of solar mass black holes, demonstrating that previously adopted limit-setting procedures have underestimated current sensitivity to the axion decay constant <e:math xmlns:e="http://www.w3.org/1998/Math/MathML" display="inline"><e:msub><e:mi>f</e:mi><e:mi>a</e:mi></e:msub></e:math> by around 1 order of magnitude and that the inclusion to higher order states allows one to reasonably capture the evolution of typical high-spin black holes across a much wider range of parameter space, thereby allowing constraints to be extended to more massive axions. We conclude with an extensive discussion on the systematics associated with spin inference from x-ray observations.