We show that gravitational wave emission from neutron star binaries can be used to discover any generic long-ranged muonic force due to the large inevitable abundance of muons inside neutron stars. As a minimal consistent example, we focus on a gauged U(1)$_{L_\mu - L_\tau}$ symmetry. In pulsar binaries, such U(1)$_{L_\mu - L_\tau}$ vectors induce an anomalously fast decay of the orbital period through the emission of dipole radiation. We study a range of different pulsar binaries, finding the most powerful constraints for vector masses below ${\cal O}(10^{-18} {\rm eV})$. For merging binaries the presence of muons in neutron stars can result in dipole radiation as well as a modification of the chirp mass during the inspiral phase. We make projections for a prospective search using both the GW170817 and S190814bv events and find that current data can discover light vectors with masses below ${\cal O}(10^{-10} {\rm eV})$. In both cases, the limits attainable with neutron stars reach gauge coupling $g^\prime\lesssim 10^{-20}$, which are many orders of magnitude stronger than previous constraints. We also show projections for next generation experiments, such as Einstein Telescope and Cosmic Explorer.

The mass hierarchy among the three generations of quarks and charged leptons is one of the greatest mysteries in particle physics. In various flavor models, the origin of this phenomenon is attributed to a series of hierarchical spontaneous symmetry breakings, most of which are beyond the reach of particle colliders. We point out that the observation of a multi-peaked stochastic gravitational wave signal from a series of cosmological phase transitions could well be a unique probe of the mechanism behind flavor hierarchies. To illustrate this point, we show how near future ground- and space-based gravitational wave observatories could detect up to three peaks in the recently proposed $PS^3$ model.

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.

We explore gravitational wave signals arising from first-order phase transitions occurring in a secluded hidden sector, allowing for the possibility that the hidden sector may have a different temperature than the Standard Model sector. We present the sensitivity to such scenarios for both current and future gravitational wave detectors in a model-independent fashion. Since secluded hidden sectors are of particular interest for dark matter models at the MeV scale or below, we pay special attention to the reach of pulsar timing arrays. Cosmological constraints on light degrees of freedom restrict the number of sub-MeV particles in a hidden sector, as well as the hidden sector temperature. Nevertheless, we find that observable first-order phase transitions can occur. To illustrate our results, we consider two minimal benchmark models: a model with two gauge singlet scalars and a model with a spontaneously broken $U(1)$ gauge symmetry in the hidden sector.

We present a model for viable gravitational reheating involving a scalar field directly coupled to the Ricci curvature scalar. Crucial to the model is a period of kination after inflation, which causes the Ricci scalar to change sign thus inducing a tachyonic effective mass $m^{2} \propto -H^2$ for the scalar field. The resulting tachyonic growth of the scalar field provides the energy for reheating, allowing for temperatures high enough for thermal leptogenesis. Additionally, the required period of kination necessarily leads to a blue-tilted primordial gravitational wave spectrum with the potential to be detected by future experiments. We find that for reheating temperatures $T_{\rm RH} \lesssim 1$ GeV, the possibility exists for the Higgs field to play the role of the scalar field.

We explore in detail the possibility that gravitational wave signals from binary inspirals are affected by a new force that couples only to dark matter particles. We discuss the impact of both the new force acting between the binary partners as well as radiation of the force carrier. We identify numerous constraints on any such scenario, ultimately concluding that observable effects on the dynamics of binary inspirals due to such a force are not possible if the dark matter is accrued during ordinary stellar evolution. Constraints arise from the requirement that the astronomical body be able to collect and bind at small enough radius an adequate number of dark matter particles, from the requirement that the particles thus collected remain bound to neutron stars in the presence of another neutron star, and from the requirement that the theory allows old neutron stars to exist and retain their charge. Thus, we show that any deviation from the predictions of general relativity observed in binary inspirals must be due either to the material properties of the inspiraling objects themselves, such as a tidal deformability, to a true fifth force coupled to baryons, or to a non-standard production mechanism for the dark matter cores of neutron stars. Viable scenarios of the latter type include production of dark matter in exotic neutron decays, or the formation of compact dark matter objects in the early Universe that later seed star formation or are captured by stars.

We propose an explanation of the galactic center gamma ray excess by supersymmetric WIMPs as heavy as 500 GeV. The lightest neutralino annihilates into vector-like leptons or quarks which cascade decay through intermediate Higgs bosons. Due to the long decay chains, the gamma ray spectrum is much softer than naively expected and peaks at GeV energies. The model predicts correlated diboson and dijet signatures to be tested at the LHC.