Producing an acceptable thermal relic abundance of dark matter with masses $\gg 10^2$ TeV is a challenge. We propose a novel mechanism where GeV-scale states establish a tiny thermal relic abundance for dark matter, which is later promoted to ultra massive status by a very light scalar. We refer to this dark matter as a THermal Ultra Massive Particle (THUMP). Direct detection of THUMPs can be naturally expected due to large scattering cross sections mediated by low mass states that couple THUMPs to the Standard Model. Our model generically leads to signals for the associated GeV-scale states at accelerator experiments.

Measurements from the Gaia satellite have greatly increased our knowledge of the dark matter velocity distributions in the Solar neighborhood. There is evidence for multiple cold structures nearby, including a high-velocity stream counterrotating relative to the Sun. This stream could significantly alter the spectrum of recoil energies and increase the annual modulation of dark matter in direct detection experiments such as DAMA/Libra. We reanalyze the experimental limits from Xenon1T, CDMSlite, PICO-60 and COSINE-100, and compare them to the results of the DAMA/Libra experiment. While we find that this new component of the dark matter velocity distribution can greatly improve the fit to the DAMA/Libra data, both spin-independent and spin-dependent interpretations of the DAMA/Libra signal with elastic and inelastic scattering continue to be ruled out by the null results of other experiments, in particular Xenon1T.

A massive $U(1)^{\prime}$ gauge boson known as a "dark photon" or $A^{\prime}$, has long been proposed as a potential explanation for the discrepancy observed between the experimental measurement and theoretical determination of the anomalous magnetic moment of the muon, ($g_{\mu} - 2$) anomaly. Recently, experimental results have excluded this possibility for a dark photon exhibiting exclusively visible or invisible decays. In this work, we revisit this idea and consider a model where $A^{\prime}$ couples inelastically to dark matter and an excited dark sector state, leading to a more exotic decay topology we refer to as a semi-visible decay. We show that for large mass splittings between the dark sector states this decay mode is enhanced, weakening the previous invisibly decaying dark photon bounds. As a consequence, $A^{\prime}$ resolves the $g_{\mu} - 2$ anomaly in a region of parameter space the thermal dark matter component of the Universe is readily explained. Interestingly, it is possible that the semi-visible events we discuss may have been vetoed by experiments searching for invisible dark photon decays. A re-analysis of the data and future searches may be crucial in uncovering this exotic decay mode or closing the window on the dark photon explanation of the $g_{\mu} - 2$ anomaly.

We present MadDM v.3.0, a numerical tool to compute particle dark matter observables in generic new physics models. The new version features a comprehensive and automated framework for dark matter searches at the interface of collider physics, astrophysics and cosmology and is deployed as a plugin of the MadGraph5_aMC@NLO platform, inheriting most of its features. With respect to the previous version, MadDM v.3.0 can now provide predictions for indirect dark matter signatures in astrophysical environments, such as the annihilation cross section at present time and the energy spectra of prompt photons, cosmic rays and neutrinos resulting from dark matter annihilation. MadDM indirect detection features support both $2\to2$ and $2 \to n$ dark matter annihilation processes. In addition, the ability to compare theoretical predictions with experimental constraints is extended by including the Fermi-LAT likelihood for gamma-ray constraints from dwarf spheroidal galaxies and by providing an interface with the nested sampling algorithm PyMultinNest to perform high dimensional parameter scans efficiently. We validate the code for a wide set of dark matter models by comparing the results from MadDM v.3.0 to existing tools and results in the literature.

We propose that neutrino masses can be zero $in ~vacuo$ and may be generated by the local distribution of dark matter through a feeble long range scalar force. We discuss potential phenomenological constraints and implications of this framework. Our model typically implies that the relic neutrino background left over from the Big Bang is mostly absent in our Galactic neighborhood. Hence, a positive detection signal from future proposed experiments, such as PTOLEMY, could in principle falsify our scenario.

Motivated by the recently reported diboson and dijet excesses in Run 1 data at ATLAS and CMS, we explore models of mixed dark matter in left-right symmetric theories. In this study, we calculate the relic abundance and the elastic scattering cross section with nuclei for a number of dark matter candidates that appear within the fermionic multiplets of left-right symmetric models. In contrast to the case of pure multiplets, WIMP-nucleon scattering proceeds at tree-level, and hence the projected reach of future direct detection experiments such as LUX-ZEPLIN and XENON1T will cover large regions of parameter space for TeV-scale thermal dark matter. Decays of the heavy charged W' boson to particles in the dark sector can potentially shift the right-handed gauge coupling to larger values when fixed to the rate of the Run 1 excesses, moving towards the theoretically attractive scenario, gR = gL. This region of parameter space may be probed by future collider searches for new Higgs bosons or electroweak fermions.

It has been proposed that a large population of unresolved millisecond pulsars (MSPs) could potentially account for the excess of GeV-scale gamma-rays observed from the region surrounding the Galactic Center. The viability of this scenario depends critically on the gamma-ray luminosity function of this source population, which determines how many MSPs Fermi should have already detected as resolved point sources. In this paper, we revisit the gamma-ray luminosity function of MSPs, without relying on uncertain distance measurements. Our determination, based on a comparison of models with the observed characteristics of the MSP population, suggests that Fermi should have already detected a significant number of sources associated with such a hypothesized Inner Galaxy population. We cannot rule out a scenario in which the MSPs residing near the Galactic Center are systematically less luminous than those present in the Galactic Plane or within globular clusters.

We explore detection prospects of a non-standard dark sector in the context of boosted dark matter. We focus on a scenario with two dark matter particles of a large mass difference, where the heavier candidate is secluded and interacts with the standard model particles only at loops, escaping existing direct and indirect detection bounds. Yet its pair annihilation in the galactic center or in the Sun may produce boosted stable particles, which could be detected as visible Cherenkov light in large volume neutrino detectors. In such models with multiple candidates, self-interaction of dark matter particles is naturally utilized in the {\it assisted freeze-out} mechanism and is corroborated by various cosmological studies such as N-body simulations of structure formation, observations of dwarf galaxies, and the small scale problem. We show that self-interaction of the secluded (heavier) dark matter greatly enhances the capture rate in the Sun and results in promising signals at current and future experiments. We perform a detailed analysis of the boosted dark matter events for Super-Kamiokande, Hyper-Kamiokande and PINGU, including notable effects such as evaporation due to self-interaction and energy loss in the Sun.

Type 1a supernova magnitudes are used to fit cosmological parameters under the assumption the model will fit the observed redshift dependence. We test this assumption with the Union 2.1 compilation of 580 sources. Several independent tests find the existing model fails to account for a significant correlation of supernova color and redshift. The correlation of magnitude residuals relative to the $\Lambda CDM$ model and $color \times redshift$ has a significance equivalent to 13 standard deviations, as evaluated by randomly shuffling the data. Extending the existing $B-V$ color correction to a relation linear in redshift improves the goodness of fit $\chi^{2}$ by more than 50 units, an equivalent 7-$\sigma$ significance, while adding only one parameter. The $color-redshift$ correlation is quite robust, cannot be attributed to outliers, and passes several tests of consistency. We review previous hints of redshift dependence in color parameters found in bin-by-bin fits interpreted as parameter bias. We show that neither the bias nor the change $\Delta \chi^{2}$ of our study can be explained by those effects. The previously known relation that bluer supernovae have larger absolute luminosity tends to empirically flatten out with increasing redshift. The best-fit cosmological dark energy density parameter is revised from $ \Omega_{\Lambda} =0.71 \pm 0.02$ to $ \Omega_{\Lambda} = 0.74 \pm 0.02$ assuming a flat universe. One possible physical interpretation is that supernovae or their environments evolve significantly with increasing redshift.