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.

$SU(2)_L \times SU(2)_R$ gauge symmetry requires three right-handed neutrinos ($ N _i $), one of which, $N_1$, can be sufficiently stable to be dark matter. In the early universe, $ W _R $ exchange with the Standard Model thermal bath keeps the right-handed neutrinos in thermal equilibrium at high temperatures. $N_1$ can make up all of dark matter if they freeze-out while relativistic and are mildly diluted by subsequent decays of a long-lived and heavier right-handed neutrino, $N_2$. We systematically study this parameter space, constraining the symmetry breaking scale of $SU(2)_R$ and the mass of $N_1$ to a triangle in the $(v_R,M_1)$ plane, with $v_R = (10^6 - 3 \times 10^{12})$ GeV and $M_1 = (2\, {\rm keV} - 1 \, {\rm MeV)}$. Much of this triangle can be probed by signals of warm dark matter, especially if leptogenesis from $N_2$ decay yields the observed baryon asymmetry. The minimal value of $v_R$ is increased to $10^8 \, {\rm GeV}$ for doublet breaking of $SU(2)_R$, and further to $10^9 \, {\rm GeV}$ if leptogenesis occurs via $N_2$ decay, while the upper bound on $M_1$ is reduced to 100 keV. In addition, there is a component of hot $N_1$ dark matter resulting from the late decay of $N_2 \rightarrow N_1 \ell^+ \ell^-$ that can be probed by future cosmic microwave background observations. Interestingly, the range of $v_R$ allows both precision gauge coupling unification and the Higgs Parity understanding of the vanishing of the Standard Model Higgs quartic at scale $v_R$. Finally, we study freeze-in production of $N_1$ dark matter via the $W_R$ interaction, which allows a much wider range of $(v_R,M_1)$.

We present the possibility that the seesaw mechanism with thermal leptogenesis can be tested using the stochastic gravitational background. Achieving neutrino masses consistent with atmospheric and solar neutrino data, while avoiding non-perturbative couplings, requires right-neutrinos lighter than the typical scale of grand unification. This scale separation suggests a symmetry protecting the right handed neutrinos from getting a mass. Thermal leptogenesis would then require that such a symmetry be broken below the reheating temperature. We enumerate all such possible symmetries consistent with these minimal assumptions and their corresponding defects, finding that in many cases, gravitational waves from the network of cosmic strings should be detectable. Estimating the predicted gravitational wave background we find that future space-borne missions could probe the entire range relevant for thermal leptogenesis.

Absorption of fermionic dark matter leads to a range of distinct and novel signatures at dark matter direct detection and neutrino experiments. We study the possible signals from fermionic absorption by nuclear targets, which we divide into two classes of four Fermi operators: neutral and charged current. In the neutral current signal, dark matter is absorbed by a target nucleus and a neutrino is emitted. This results in a characteristically different nuclear recoil energy spectrum from that of elastic scattering. The charged current channel leads to induced $\beta$ decays in isotopes which are stable in vacuum as well as shifts of the kinematic endpoint of $ \beta$ spectra in unstable isotopes. To confirm the possibility of observing these signals in light of other constraints, we introduce UV completions of example higher dimensional operators that lead to fermionic absorption signals and study their phenomenology. Most prominently, dark matter which exhibits fermionic absorption signals is necessarily unstable leading to stringent bounds from indirect detection searches. Nevertheless, we find a large viable parameter space in which dark matter is sufficiently long lived and detectable in current and future experiments.

Pulsars act as accurate clocks, sensitive to gravitational redshift and acceleration induced by transiting clumps of matter. We study the sensitivity of pulsar timing arrays (PTAs) to single transiting compact objects, focusing on primordial black holes and compact subhalos in the mass range from $10^{-12} M _{\odot}$ to well above $100~M_\odot$. We find that the Square Kilometer Array can constrain such objects to be a subdominant component of the dark matter over this entire mass range, with sensitivity to a dark matter sub-component reaching the sub-percent level over significant parts of this range. We also find that PTAs offer an opportunity to probe substantially less dense objects than lensing because of the large effective radius over which such objects can be observed, and we quantify the subhalo concentration parameters which can be constrained.

Vector bosons heavier than $10^{-22}$ eV can be viable dark matter candidates with distinctive experimental signatures. Ultralight dark matter generally requires a non-thermal origin to achieve the observed density, while still behaving like a pressureless fluid at late times. We show that such a production mechanism naturally occurs for vectors whose mass originates from a dark Higgs. If the dark Higgs has a large field value after inflation, the energy in the Higgs field can be efficiently transferred to vectors through parametric resonance. Computing the resulting abundance and spectra requires careful treatment of the transverse and longitudinal components, whose dynamics are governed by distinct differential equations. We study these equations in detail and find that the mass of the vector may be as low as $10 ^{ - 18 }$ eV, while making up the dominant dark matter abundance. This opens up a wide mass range of vector dark matter as cosmologically viable, further motivating their experimental search.

If supersymmetry exists at any scale, regardless of whether it is restored around the weak scale, it may be a good symmetry of the dark sector, enforcing a degeneracy between its lowest lying fermions and bosons. We explore the implications of this scenario for the early universe and dark matter, as well as the corresponding signatures. In particular we show that the thermal history of the dark sector results in co-decaying dark matter in much of the parameter space. This implies new phenomenological signatures and presents a new way to discover high scale supersymmetry.

We study the cosmological consequences of co-decaying dark matter - a recently proposed mechanism for depleting the density of dark matter through the decay of nearly degenerate particles. A generic prediction of this framework is an early dark matter dominated phase in the history of the universe, that results in the enhanced growth of dark matter perturbations on small scales. We compute the duration of the early matter dominated phase and show that the perturbations are robust against washout from free-streaming. The enhanced small scale structure is expected to survive today in the form of compact micro-halos and can lead to significant boost factors for indirect detection experiments, such as FERMI, where dark matter would appear as point sources.

We propose a new mechanism for thermal dark matter freezeout, termed Co-Decaying Dark Matter. Multi-component dark sectors with degenerate particles and out-of-equilibrium decays can co-decay to obtain the observed relic density. The dark matter density is exponentially depleted through the decay of nearly degenerate particles, rather than from Boltzmann suppression. The relic abundance is set by the dark matter annihilation cross-section, which is predicted to be boosted, and the decay rate of the dark sector particles. The mechanism is viable in a broad range of dark matter parameter space, with a robust prediction of an enhanced indirect detection signal. Finally, we present a simple model that realizes co-decaying dark matter.