The superradiant instability of black hole space-times has been used to place limits on ultra-light bosonic particles. We show that these limits are model dependent. While the initial growth of the mode is gravitational and thus model independent, the ability to place a limit on new particles requires the mode to grow unhindered to a large number density. Non-linear interactions between the particle and other light degrees of freedom that are mediated through higher dimension operators can damp this growth, eliminating the limit. However, these non-linearities may also destroy a cosmic abundance of these light particles, an attractive avenue for their discovery in several experiments. We study the specific example of the QCD axion and show that it is easy to construct models where these non-linearities eliminate limits from superradiance while preserving their cosmic abundance.

Gravitational waves (GW) produced in the early Universe contribute to the number of relativistic degrees of freedom, $N_{\rm eff}$, during Big Bang Nucleosynthesis (BBN). By using the constraints on $N_{\rm eff}$, we present a new bound on how much the observable Universe could have expanded during cosmic inflation. The new bound is roughly four orders of magnitude more stringent than previous bounds. We also discuss the sensitivities of the current and planned GW observatories such as LIGO and LISA, and show that the constraints they could impose are always less stringent than the BBN bound.

From a model-building perspective, the inflationary sector might very well have no direct couplings to other species, apart from inevitable gravitational interactions. Within the context of General Relativity, a thermal universe can still emerge after inflation if: $i)$ some radiation sector is excited towards the end of inflation, and $ii)$ the post-inflationary equation of state becomes sufficiently stiff $w \geq w_{\rm RD}\gtrsim 0.57$, with $w_{\rm RD}$ a threshold depending on the inflationary scale $H_*$ and the initial radiation-to-inflaton energy ratio $\Delta_*$. Furthermore, a stiff period in the expansion history enhances significantly the inflationary gravitational wave (GW) background, making this signal (potentially) observable by aLIGO, LISA and other experiments. The very same enhancement leads however to an inconsistency of the scenario: the energy of the GWs becomes too large compared to the rest of the radiation sector, violating standard BBN and CMB bounds on GW backgrounds. Except for very special scenarios where the initial radiation sector comprises hundreds of fields with couplings tuned to specific values, our result applies independently of $w$, $H_*$ and $\Delta_*$. This suggests that in order to reheat the universe, the inflationary sector should be coupled directly to other particle species. Alternatively the inflationary sector could be implemented in modified gravity theories.

The expansion history of the Universe between the end of inflation and the onset of radiation-domination (RD) is currently unknown. If the equation of state during this period is stiffer than that of radiation, $w > 1/3$, the gravitational wave (GW) background from inflation acquires a blue-tilt ${d\log\rho_{\rm GW}\over d\log f} = {2(w-1/3)\over (w+1/3)} > 0$ at frequencies $f \gg f_{\rm RD}$ corresponding to modes re-entering the horizon during the stiff-domination (SD), where $f_{\rm RD}$ is the frequency today of the horizon scale at the SD-to-RD transition. We characterized in detail the transfer function of the GW energy density spectrum, considering both 'instant' and smooth modelings of the SD-to-RD transition. The shape of the spectrum is controlled by $w$, $f_{\rm RD}$, and $H_{\rm inf}$ (the Hubble scale of inflation). We determined the parameter space compatible with a detection of this signal by LIGO and LISA, including possible changes in the number of relativistic degrees of freedom, and the presence of a tensor tilt. Consistency with upper bounds on stochastic GW backgrounds, however, rules out a significant fraction of the observable parameter space. We find that this renders the signal unobservable by Advanced LIGO, in all cases. The GW background remains detectable by LISA, though only in a small island of parameter space, corresponding to scenarios with an equation of state in the range $0.46 \lesssim w \lesssim 0.56$ and a high inflationary scale $H_{\rm inf} \gtrsim 10^{13}~{\rm GeV}$, but low reheating temperature $1~{\rm MeV} \lesssim T_{\rm RD} \lesssim 150~{\rm MeV}$ (equivalently, $10^{-11}~{\rm Hz} \lesssim f_{\rm RD} \lesssim 3.6\cdot10^{-9}~{\rm Hz}$). Implications for early Universe scenarios resting upon an SD epoch are briefly discussed.

The damping process of a homogeneous oscillating scalar field that indirectly interacts with a thermal bath through a mediator field is investigated over a wide range of model parameters. We consider two types of mediator fields, those that can decay to the thermal bath and those that are individually stable but pair annihilate. The former case has been extensively studied in the literature by treating the damping as a local effect after integrating out the assumed close-to-equilibrium mediator field. The same approach does not apply if the mediator field is stable and freezes out of equilibrium. To account for the latter case, we adopt a non-local description of damping that is only meaningful when we consider full half-oscillations of the field being damped. The damping rates of the oscillating scalar field and the corresponding heating rate of the thermal bath in all bulk parameter regions are calculated in both cases, corroborating previous results in the direct decay case. Using the obtained results, the time it takes for the amplitude of the scalar field to be substantially damped is estimated.