We pioneer the black hole mass gap as a powerful new tool for constraining new particles. A new particle that couples to the Standard Model---such as an axion---acts as an additional source of loss in the cores of population-III stars, suppressing mass lost due to winds and quenching the pair-instability. This results in heavier astrophysical black holes. As an example, using stellar simulations we show that the solar axion explanation of the recent XENON1T excess implies astrophysical black holes of ~ 56 MS, squarely within the black hole mass gap predicted by the Standard Model.

New light particles produced in supernovae can lead to additional energy loss and a consequent deficit in neutrino production in conflict with the neutrinos observed from Supernova 1987A (SN1987A). Contrary to the majority of previous SN1987A studies, we examine the impact of $Z'$ bosons, axions, and axion-like particles (ALPs) interacting with the muons produced in SN1987A. For the first time, we find constraints on generic $Z'$ bosons coupled to muons, and apply our results to particle models including gauged $L_\mu-L_\tau$ number, $U(1)_{L_\mu-L_\tau}$, and gauged $B-L$ number, $U(1)_{B-L}$. We constrain $Z'$ bosons with masses up to about 250-500 MeV, and down to about $10^{-9}$ in $Z'$-muon coupling. We also extend previous work on axion-muon couplings by examining the importance of loop-level interactions, as well as performing calculations over a wider range of axion masses. We constrain muon-coupled axions from arbitrarily low masses up to about 200-500 MeV, with bounds extending down to axion-muon couplings of approximately $10^{-8}$ GeV$^{-1}$. We conclude that supernovae broadly provide a sensitive probe of new lightly-coupled particles interacting with muons.

We study the gravitational wave phenomenology in models of solitosynthesis. In such models, a first order phase transition is precipitated by a period in which non-topological solitons with a conserved global charge (Q-balls) accumulate charge. As such, the nucleation rate of critical bubbles differs significantly from thermal phase transitions. In general we find that the peak amplitude of the gravitational wave spectrum resulting from solitosynthesis is stronger than that of a thermal phase transition and the timescale of the onset of nonlinear plasma dynamics is comparable to Hubble. We demonstrate this explicitly in an asymmetric dark matter model, and discuss current and future constraints in this scenario.

Dark matter may be in the form of non-baryonic structures such as compact subhalos and boson stars. Structures weighing between asteroid and solar masses may be discovered via gravitational microlensing, an astronomical probe that has in the past helped constrain the population of primordial black holes and baryonic MACHOs. We investigate the non-trivial effect of the size of and density distribution within these structures on the microlensing signal, and constrain their populations using the EROS-2 and OGLE-IV surveys. Structures larger than a solar radius are generally constrained more weakly than point-like lenses, but stronger constraints may be obtained for structures with mass distributions that give rise to caustic crossings or produce larger magnifications.

We explore a simple model which naturally explains the observed baryon asymmetry of the Universe. In this model the strong coupling is promoted to a dynamical quantity, which evolves through the vacuum expectation value of a singlet scalar field that mixes with the Higgs field. In the resulting cosmic history, QCD confinement and electroweak symmetry breaking initially occur simultaneously close to the weak scale. The early confinement triggers the axion to roll toward its minimum, which creates a chemical potential between baryons and antibaryons through the interactions of the $\eta'$ meson, resulting in spontaneous baryogenesis. The electroweak sphalerons are sharply switched off after confinement and the baryon asymmetry is frozen in. Subsequently, evolution of the Higgs vacuum expectation value (which is modified in the confined phase) triggers a relaxation to a Standard Model-like vacuum. We identify viable regions of parameter space, and describe various experimental probes, including current and future collider constraints, and gravitational wave phenomenology.

In this white paper, we discuss the prospects for characterizing and identifying dark matter using gravitational waves, covering a wide range of dark matter candidate types and signals. We argue that present and upcoming gravitational wave probes offer unprecedented opportunities for unraveling the nature of dark matter and we identify the most urgent challenges and open problems with the aim of encouraging a strong community effort at the interface between these two exciting fields of research.

In this paper we explore the possibility of observable gravitational waves as a manifestation of the QCD axion dynamics. In particular, we focus on dynamical axion models which solve the strong CP problem, and include the confinement of a QCD-like gauge group at the TeV scale. We study the resulting chiral symmetry breaking phase transition for models with $N_F=3$ and $N_F=4$ light flavors using the linear sigma model. This model describes the scalar meson spectrum and its interactions, with the diagonal field $\varphi$ as the order parameter. We find that the amplitude of the gravitational wave spectrum depends on the mass of the dynamical axion $\eta'$ via the ratio $m_{\eta'}/m_\varphi$. The resulting spectra may be observed at future mid-range gravitational wave experiments such as AION/MAGIS, DECIGO, and BBO. Moreover, the TeV states can be searched for at colliders and their quantum numbers characterized, providing a unique connection between axion physics, gravitational waves and collider searches.

We study properties of boson stars consisting of ultra-light scalar dark matter with repulsive self-interactions. We investigate the origin of the maximum mass of spherically symmetric stable stars which emerges only when solving the full equations of motion in curved space-time, but not when solving the approximated Schr\"odinger-Newton equations. When the repulsion is weak, the backreaction of the curvature on the scalars acts as an additional source of attraction and can overcome the repulsion, resulting in a maximum star mass and compactness. We also point out that the potential in a UV completed particle physics model of light scalar dark matter is generally more complicated than the widely used $\phi^4$ interaction. Additional interactions beyond $\phi^4$ in the potential can dramatically change the properties of boson stars as well as modify the prospect of LIGO gravitational wave detection for binary mergers of boson stars.

We analyse the gravitational wave and low energy signatures of a Pati-Salam phase transition. For a Pati-Salam scale of $M_{PS} \sim 10^5$ GeV, we find a stochastic power spectrum within reach of the next generation of ground-based interferometer experiments such as the Einstein Telescope, in parts of the parameter space. We study the lifetime of the proton in this model, as well as complementarity with low energy constraints including electroweak precision data, neutrino mass measurements, lepton flavour violation, and collider constraints.

Models in which the dark matter is produced at extremely low rates from the annihilation of Standard Model particles in the early Universe allow us to explain the current dark matter relic density while easily evading the traditional experimental constraints. In scenarios where the dark matter interacts with the Standard Model via a new physics mediator, the early Universe dynamics of the dark sector can be particularly complex, as the dark matter and the mediator could be in thermal and chemical equilibrium with each other. This equilibration takes place via number-changing processes such as double Compton scattering and bremsstrahlung, whose amplitudes are cumbersome to calculate. In this paper, we show that in large regions of the parameter space, these equilibration mechanisms do not significantly affect the final dark matter relic density. In particular, for a model with a light dark photon mediator, the relic density can be reasonably estimated by considering that the dark matter is solely produced through the annihilation of Standard Model particles. This result considerably simplifies the treatment of a large class of dark matter theories, facilitating in particular the superimposition of the relic density constraints on the current and future experimental bounds.

We demonstrate that the relic gravitational wave background from a multi-step phase transition may deviate from the simple sum of the single spectra, for phase transitions with similar nucleation temperatures $T_N$. We demonstrate that the temperature range $\Delta T$ between the volume fractions $f(T)=0.1$ and $f(T)=0.9$ occupied by the vacuum bubbles can span $\sim 20$ GeV. This allows for a situation in which phase transitions overlap, such that the later bubbles may nucleate both in high temperature and intermediate temperature phases. Such scenarios may lead to more exotic gravitational wave spectra, which cannot be fitted that of a consecutive PTs. We demonstrate this explicitly in the singlet extension of the Standard Model. Finally, we comment on potential additional effects due to the more exotic dynamics of overlapping phase transitions.

In anticipation of upcoming gravitational wave experiments, we provide a comprehensive overview of the spectra predicted by phase transitions triggered by states from a large variety of dark sector models. Such spectra are functions of the quantum numbers and (self-) couplings of the scalar that triggers the dark phase transition. We classify dark sectors that give rise to a first order phase transition and perform a numerical scan over the thermal parameter space. We then characterize scenarios in which a measurement of a new source of gravitational waves could allow us to discriminate between models with differing particle content.

We calculate the gravitational radiation background generated from boson star binaries formed in locally dense clusters with formation rate tracked by the regular star formation rate. We compute how the the frequency window in gravitational waves is affected by the boson field mass and repulsive self-coupling, anticipating constraints from EPTA and LISA. We also comment on the possible detectability of these binaries.

The flatness of the inflaton potential and lightness of the Higgs could have the common origin of the breaking of a global symmetry. This scenario provides a unified framework of Goldstone Inflation and Composite Higgs, where the inflaton and the Higgs both have a pseudo--Goldstone boson nature. The inflaton reheats the Universe via decays to the Higgs and subsequent secondary production of other SM particles via the top and massive vector bosons. We find that inflationary predictions and perturbative reheating conditions are consistent with CMB data for sub--Planckian values of the fields, as well as opening up the possibility of inflation at the TeV scale. We explore this exciting possibility, leading to an interplay between collider data and cosmological constraints.

Identifying the inflaton with a pseudo-Goldstone boson explains the flatness of its potential. Successful Goldstone Inflation should also be robust against UV corrections, such as from quantum gravity: in the language of the effective field theory this implies that all scales are sub-Planckian. In this paper we present scenarios which realise both requirements by examining the structure of Goldstone potentials arising from Coleman-Weinberg contributions. We focus on single-field models, for which we notice that both bosonic and fermionic contributions are required and that spinorial fermion representations can generate the right potential shape. We then evaluate the constraints on non-Gaussianity from higher-derivative interactions, finding that axiomatic constraints on Goldstone boson scattering prevail over the current CMB measurements. The fit to CMB data can be connected to the UV completions for Goldstone Inflation, finding relations in the spectrum of new resonances. Finally, we show how hybrid inflation can be realised in the same context, where both the inflaton and the waterfall fields share a common origin as Goldstones.

We show that General Relativity coupled to a quantum field theory generically leads to non-local effects in the matter sector. These non-local effects can be described by non-local higher dimensional operators which remarkably have an approximate shift symmetry. When applied to inflationary models, our results imply that small non-Gaussianities are a generic feature of models based on General Relativity coupled to matter fields. However, these effects are too small to be observable in the Cosmic Microwave Background.

Slow-roll inflation requires the inflaton field to have an exceptionally flat potential, which combined with measurements of the scale of inflation demands some degree of fine-tuning. Alternatively, the flatness of the potential could be due to the inflaton's origin as a pseudo-Goldstone boson, as in Natural Inflation. Alas, consistency with Planck data places the original proposal of Natural Inflation in a tight spot, as it requires a trans-Planckian excursion of the inflaton. Although one can still tune the renormalizable potential to sub-Planckian values, higher order corrections from quantum gravity or sources of breaking of the Goldstone symmetry would ruin the predictivity of the model. In this paper we show how in more realistic models of Natural Inflation one could achieve inflation without a trans-Planckian excursion of the field. We show how a variant of Extra-natural inflation with bulk fermions can achieve the desired goal and discuss its four-dimensional duals. We also present a new type of four dimensional models inspired in Little Higgs and Composite Higgs models which can lead to sub-Planckian values of the inflaton field.

We discuss cosmological inflation in the minimal Wess-Zumino model with a single massive chiral supermultiplet. With suitable parameters and assuming a plausible initial condition at the start of the inflationary epoch, the model can yield scalar perturbations in the Cosmic Microwave Background (CMB) of the correct strength with a spectral index n_s ~ 0.96 and a tensor-to-scalar perturbation ratio r < 0.1, consistent with the Planck CMB data. We also discuss the possibility of topological inflation within the Wess-Zumino model, and the possibility of combining it with a seesaw model for neutrino masses. This would violate R-parity, but at such a low rate that the lightest supersymmetric particle would have a lifetime long enough to constitute the astrophysical cold dark matter.