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.

I present a short overview of the latest developments in indirect searches for dark matter using gamma rays, X-rays, charged cosmic rays, micro waves, radio waves, and neutrinos. I briefly outline key past, present, and future experiments and their search strategies. In several searches there are exciting anomalies which could potentially be emerging dark matter signals. I discuss these anomalous signals, and some future prospects to determine their origins.

Neutron stars serve as excellent next-generation thermal detectors of dark matter, heated by the scattering and annihilation of dark matter accelerated to relativistic speeds in their deep gravitational wells. However, the dynamics of neutron star cores are uncertain, making it difficult at present to unequivocally compute dark matter scattering in this region. On the other hand, the physics of an outer layer of the neutron star, the crust, is more robustly understood. We show that dark matter scattering solely with the low-density crust still kinetically heats neutron stars to infrared temperatures detectable by forthcoming telescopes. We find that for both spin-independent and spin-dependent scattering on nucleons, the crust-only cross section sensitivity is $10^{-43} - 10^{-41}$~cm$^2$ for dark matter masses of 100 MeV $-$ 1 PeV, with the best sensitivity arising from dark matter scattering with a crust constituent called nuclear pasta (including gnocchi, spaghetti, and lasagna phases). For dark matter masses from 10 eV to 1 MeV, the sensitivity is $10^{-39} - 10^{-34}$~cm$^2$, arising from exciting collective phonon modes in a neutron superfluid in the inner crust. Furthermore, for any $s$-wave or $p$-wave annihilating dark matter, we show that dark matter will efficiently annihilate by thermalizing just with the neutron star crust, regardless of whether the dark matter ever scatters with the neutron star core. This implies efficient annihilation in neutron stars for any electroweakly interacting dark matter with inelastic mass splittings of up to 200 MeV, including Higgsinos. We conclude that neutron star crusts play a key role in dark matter scattering and annihilation in neutron stars.

The Galactic Center GeV excess (GCE) has garnered great interest as a possible signal of either dark matter annihilation or some novel astrophysical phenomenon, such as a new population of gamma-ray emitting pulsars. In a companion paper, we showed that in a $10^\circ$ radius region of interest (ROI) surrounding the Galactic Center, apparent evidence for GCE point sources (PSs) from non-Poissonian template fitting (NPTF) is actually an artifact of unmodeled north-south asymmetry of the GCE. In this work, we develop a simplified analytic description of how signal mismodeling can drive an apparent preference for a PS population, and demonstrate how the behavior pointed out in the companion paper also appears in simpler simulated datasets that contain no PS signals at all. We explore the generality of this behavior in the real gamma-ray data, and discuss the implications for past and future studies using NPTF techniques. While the drop in PS preference once north-south asymmetry is included is not ubiquitous in larger ROIs, we show that any overly-rigid signal model is expected to yield a spurious PS signal that can appear very convincing: as well as apparent significance comparable to what one would expect from a true PS population, the signal can exhibit stability against a range of variations in the analysis, and a source count function that is very consistent with previous apparent NPTF-based detections of a GCE PS population. This contrasts with previously-studied forms of systematic mismodeling which are unlikely to mimic a PS population in the same way. In the light of this observation, and its explicit realization in the region where the GCE is brightest, we argue that a dominantly smooth origin for the GCE is not in tension with existing NPTF analyses.

We re-examine evidence that the Galactic Center Excess (GCE) originates primarily from point sources (PSs). We show that in our region of interest, non-Poissonian template fitting (NPTF) evidence for GCE PSs is an artifact of unmodeled north-south asymmetry of the GCE. This asymmetry is strongly favored by the fit (although it is unclear if this is physical), and when it is allowed, the preference for PSs becomes insignificant. We reproduce this behavior in simulations, including detailed properties of the spurious PS population. We conclude that NTPF evidence for GCE PSs is highly susceptible to certain systematic errors, and should not at present be taken to robustly disfavor a dominantly smooth GCE.

Statistical evidence has previously suggested that the Galactic Center GeV Excess (GCE) originates largely from point sources, and not from annihilating dark matter. We examine the impact of unmodeled source populations on identifying the true origin of the GCE using non-Poissonian template fitting (NPTF) methods. In a proof-of-principle example with simulated data, we discover that unmodeled sources in the Fermi Bubbles can lead to a dark matter signal being misattributed to point sources by the NPTF. We discover striking behavior consistent with a mismodeling effect in the real Fermi data, finding that large artificial injected dark matter signals are completely misattributed to point sources. Consequently, we conclude that dark matter may provide a dominant contribution to the GCE after all.

The Sun is an excellent laboratory for astroparticle physics but remains poorly understood at GeV--TeV energies. Despite the immense relevance for both cosmic-ray propagation and dark matter searches, only in recent years has the Sun become a target for precision gamma-ray astronomy with the Fermi-LAT instrument. Among the most surprising results from the observations is a hard excess of GeV gamma-ray flux that strongly anti-correlates with solar activity, especially at the highest energies accessible to Fermi-LAT. Most of the observed properties of the gamma-ray emission cannot be explained by existing models of cosmic-ray interactions with the solar atmosphere. GeV--TeV gamma-ray observations of the Sun spanning an entire solar cycle would provide key insights into the origin of these gamma rays, and consequently improve our understanding of the Sun's environment as well as the foregrounds for new physics searches, such as dark matter. These can be complemented with new observations with neutrinos and cosmic rays. Together these observations make the Sun a new testing ground for particle physics in dynamic environments.

We analyze the Sun as a source for the indirect detection of dark matter through a search for gamma rays from the solar disk. Capture of dark matter by elastic interactions with the solar nuclei followed by annihilation to long-lived mediators can produce a detectable gamma-ray flux. We search three years of data from the High Altitude Water Cherenkov Observatory and find no statistically significant detection of TeV gamma-ray emission from the Sun. Using this, we constrain the spin-dependent elastic scattering cross section of dark matter with protons for dark matter masses above 1 TeV, assuming an unstable mediator with a favorable lifetime. The results complement constraints obtained from Fermi-LAT observations of the Sun and together cover WIMP masses between 4 GeV and $10^6$ GeV. The cross section constraints for mediator decays to gamma rays can be as strong as $\sim10^{-45}$ cm$^{-2}$, which is more than four orders of magnitude stronger than current direct-detection experiments for 1 TeV dark matter mass. The cross-section constraints at higher masses are even better, nearly 7 orders of magnitude better than the current direct-detection constraints for 100 TeV dark matter mass. This demonstration of sensitivity encourages detailed development of theoretical models in light of these powerful new constraints.

Steady gamma-ray emission up to at least 200 GeV has been detected from the solar disk in the Fermi-LAT data, with the brightest, hardest emission occurring during solar minimum. The likely cause is hadronic cosmic rays undergoing collisions in the Sun's atmosphere after being redirected from ingoing to outgoing in magnetic fields, though the exact mechanism is not understood. An important new test of the gamma-ray production mechanism will follow from observations at higher energies. Only the High Altitude Water Cherenkov (HAWC) Observatory has the required sensitivity to effectively probe the Sun in the TeV range. Using three years of HAWC data from November 2014 to December 2017, just prior to the solar minimum, we search for 1--100 TeV gamma rays from the solar disk. No evidence of a signal is observed, and we set strong upper limits on the flux at a few $10^{-12}$ TeV$^{-1}$ cm$^{-2}$ s$^{-1}$ at 1 TeV. Our limit, which is the most constraining result on TeV gamma rays from the Sun, is $\sim10\%$ of the theoretical maximum flux (based on a model where all incoming cosmic rays produce outgoing photons), which in turn is comparable to the Fermi-LAT data near 100 GeV. The prospects for a first TeV detection of the Sun by HAWC are especially high during solar minimum, which began in early 2018.

Weakly Interacting Massive Particles (WIMPs) have long reigned as one of the leading classes of dark matter candidates. The observed dark matter abundance can be naturally obtained by freezeout of weak-scale dark matter annihilations in the early universe. This "thermal WIMP" scenario makes direct predictions for the total annihilation cross section that can be tested in present-day experiments. While the dark matter mass constraint can be as high as $m_\chi\gtrsim100$ GeV for particular annihilation channels, the constraint on the total cross section has not been determined. We construct the first model-independent limit on the WIMP total annihilation cross section, showing that allowed combinations of the annihilation-channel branching ratios considerably weaken the sensitivity. For thermal WIMPs with s-wave $2\rightarrow2$ annihilation to visible final states, we find the dark matter mass is only known to be $m_\chi\gtrsim20$ GeV. This is the strongest largely model-independent lower limit on the mass of thermal-relic WIMPs, together with the upper limit on the mass from the unitarity bound ($m_\chi\lesssim 100$ TeV), it defines what we call the "WIMP window". To probe the remaining mass range, we outline ways forward.

Many dark matter interaction types lead to annihilation processes which suffer from p-wave suppression or helicity suppression, rendering them subdominant to unsuppressed s-wave processes. We demonstrate that the natural inclusion of dark initial state radiation can open an unsuppressed s-wave annihilation channel, and thus provide the dominant dark matter annihilation process for particular interaction types. We illustrate this effect with the bremsstrahlung of a dark spin-0 or dark spin-1 particle from fermionic dark matter, $\overline{\chi}\chi\rightarrow \overline{f}f\phi$ or $\overline{f}fZ'$. The dark initial state radiation process, despite having a 3-body final state, proceeds at the same order in the new physics scale $\Lambda$ as the annihilation to the 2-body final state $\overline{\chi}\chi\rightarrow \overline{f}f$. This is lower order in $\Lambda$ than the well-studied lifting of helicity suppression via Standard Model final state radiation, or virtual internal bremsstrahlung. This dark bremsstrahlung process should influence LHC and indirect detection searches for dark matter.

Dark matter capture and annihilation in the Sun can produce detectable high-energy neutrinos, providing a probe of the dark matter-proton scattering cross section. We consider the case when annihilation proceeds via long-lived dark mediators, which allows gamma rays to escape the Sun and reduces the attenuation of neutrinos. For gamma rays, there are exciting new opportunities, due to detailed measurements of GeV solar gamma rays with Fermi, and unprecedented sensitivities in the TeV range with HAWC and LHAASO. For neutrinos, the enhanced flux, particularly at higher energies ($\sim$TeV), allows a more sensitive dark matter search with IceCube and KM3NeT. We show that these search channels can be extremely powerful, potentially improving sensitivity to the dark matter spin-dependent scattering cross section by several orders of magnitude relative to present searches for high-energy solar neutrinos, as well as direct detection experiments.

In the simplified dark matter models commonly studied, the mass generation mechanism for the dark fields is not typically specified. We demonstrate that the dark matter interaction types, and hence the annihilation processes relevant for relic density and indirect detection, are strongly dictated by the mass generation mechanism chosen for the dark sector particles, and the requirement of gauge invariance. We focus on the class of models in which fermionic dark matter couples to a spin-1 vector or axial-vector mediator. However, in order to generate dark sector mass terms, it is necessary in most cases to introduce a dark Higgs field and thus a spin-0 scalar mediator will also be present. In the case that all the dark sector fields gain masses via coupling to a single dark sector Higgs field, it is mandatory that the axial-vector coupling of the spin-1 mediator to the dark matter is non-zero; the vector coupling may also be present depending on the charge assignments. For all other mass generation options, only pure vector couplings between the spin-1 mediator and the dark matter are allowed. If these coupling restrictions are not obeyed, unphysical results may be obtained such as a violation of unitarity at high energies. These two-mediator scenarios lead to important phenomenology that does not arise in single mediator models. We survey two-mediator dark matter models which contain both vector and scalar mediators, and explore their relic density and indirect detection phenomenology.

We consider the indirect detection signals for a self-consistent hidden $U(1)$ model containing a Majorana dark matter candidate, $\chi$, a dark gauge boson, $Z'$, and a dark Higgs, $s$. Compared with a model containing only a dark matter candidate and $Z'$ mediator, the addition of the scalar provides a mass generation mechanism for the dark sector particles and is required in order to avoid unitarity violation at high energies. We find that the inclusion of the two mediators opens up a new two-body s-wave annihilation channel, $\chi\overline\chi\rightarrow sZ'$. This new process, which is missed in the usual single-mediator simplified model approach, can be the dominant annihilation channel. This provides rich phenomenology for indirect detection searches, allows indirect searches to explore regions of parameter space not accessible with other commonly considered s-wave annihilation processes, and enables both the $Z'$ and scalar couplings to be probed. We examine the phenomenology of the sector with a focus on this new process, and determine the limits on the model parameter space from Fermi data on dwarf spheriodal galaxies and other relevant experiments.

We consider a scenario where dark matter (DM) interacts exclusively with Standard Model (SM) leptons at tree level. Due to the absence of tree-level couplings to quarks, the constraints on leptophilic dark matter arising from direct detection and hadron collider experiments are weaker than those for a generic WIMP. We study a simple model in which interactions of DM with SM leptons are mediated by a leptophilic $Z'$ boson, and determine constraints on this scenario arising from relic density, direct detection, and other experiments. We then determine current LHC limits and project the future discovery reach. We show that, despite the absence of direct interactions with quarks, this scenario can be strongly constrained.