Astrophysical observations of spinning BHs, which span $ 5M_\odot\lesssim M_{\rm BH}\lesssim 5\times 10^8 M_\odot$, can be used to exclude the existence of certain massive bosons via the superradiance phenomenon. In this work, we explore for the first time how these measurements can be used to constrain properties of statistical distributions for the masses of multiple bosonic fields. Quite generally, our methodology excludes $N_{\rm ax}\gtrsim 30$ scalar fields with a range of mass distribution widths and central values spanning many orders of magnitude. We demonstrate this for the specific example of axions in string theory and M-theory, where the mass distributions in certain cases take universal forms. We place upper bounds on $N_{\rm ax}$ for certain scenarios of interest realised approximately as mass distributions in M-theory, including the QCD axion, grand unified theories, and fuzzy dark matter.

Axions arise in many theoretical extensions of the Standard Model of particle physics, in particular the "string axiverse". If the axion masses, $m_a$, and (effective) decay constants, $f_a$, lie in specific ranges, then axions contribute to the cosmological dark matter and dark energy densities. We compute the background cosmological (quasi-)observables for models with a large number of axion fields, $n_{\rm ax}\sim \mathcal{O}(10-100)$, with the masses and decay constants drawn from statistical distributions. This reduces the number of parameters from $2n_{\rm ax}$ to a small number of "hyperparameters". We consider a number of distributions, from those motivated purely by statistical considerations, to those where the structure is specified according to a class of M-theory models. Using Bayesian methods we are able to constrain the hyperparameters of the distributions. In some cases the hyperparameters can be related to string theory, e.g. constraining the number ratio of axions to moduli, or the typical decay constant scale needed to provide the correct relic densities. Our methodology incorporates the use of both random matrix theory and Bayesian networks.

We study the relative contribution of cusps and pseudocusps, on cosmic (super)strings, to the emitted bursts of gravitational waves. The gravitational wave emission in the vicinity of highly relativistic points on the string follows, for a high enough frequency, a logarithmic decrease. The slope has been analytically found to be $^{-4}/_3$ for points reaching exactly the speed of light in the limit $c=1$. We investigate the variations of this high frequency behaviour with respect to the velocity of the points considered, for strings formed through a numerical simulation, and we then compute numerically the gravitational waves emitted. We find that for string points moving with velocities as far as $10^{-3}$ from the theoretical (relativistic) limit $c=1$, gravitational wave emission follows a behaviour consistent with that of cusps, effectively increasing the number of cusps on a string. Indeed, depending on the velocity threshold chosen for such behaviour, we show the emitting part of the string worldsheet is enhanced by a factor ${\cal O}(10^3)$ with respect to the emission of cusps only.