Abstract: Although weakly interacting massive particles (WIMPs) have long been among the most studied and theoretically attractive classes of candidates for the dark matter of our universe, the lack of their detection in direct detection and collider experiments has begun to dampen enthusiasm for this paradigm. In this study, we set out to appraise the status of the WIMP paradigm, focusing on the case of dark matter candidates that interact with the Standard Model through a new gauge boson. After considering a wide range of Z' mediated dark matter models, we quantitatively evaluate the fraction of the parameter space that has been excluded by existing experiments, and that is projected to fall within the reach of future direct detection experiments. Despite the existence of stringent constraints, we find that a sizable fraction of this parameter space remains viable. More specifically, if the dark matter is a Majorana fermion, we find that an order one fraction of the parameter space is in many cases untested by current experiments. Future direct detection experiments with sensitivity near the irreducible neutrino floor will be able to test a significant fraction of the currently viable parameter space, providing considerable motivation for the next generation of direct detection experiments.

Abstract: We present halo-independent methods to analyze the results of dark matter direct detection experiments assuming inelastic scattering. We focus on the annual modulation signal reported by DAMA/LIBRA and present three different halo-independent tests. First, we compare it to the upper limit on the unmodulated rate from XENON100 using (a) the trivial requirement that the amplitude of the annual modulation has to be smaller than the bound on the unmodulated rate, and (b) a bound on the annual modulation amplitude based on an expansion in the Earth's velocity. The third test uses the special predictions of the signal shape for inelastic scattering and allows for an internal consistency check of the data without referring to any astrophysics. We conclude that a strong conflict between DAMA/LIBRA and XENON100 in the framework of spin-independent inelastic scattering can be established independently of the local properties of the dark matter halo.

Abstract: In the present work, we show that the Galactic Center Excess (GCE) emission, as recently updated by the Fermi-LAT Collaboration, could be explained by a mixture of Fermi bubbles-like emission plus dark matter (DM) annihilation, in the context of a scalar-singlet Higgs portal scenario (SHP). In fact, the standard SHP, where the DM particle, S, only has renormalizable interactions with the Higgs, is non-operational due to strong constraints, especially from DM direct detection limits. Thus we consider the most economical extension, called ESHP (for extended SHP), which consists solely in the addition of a second (more massive) scalar singlet in the dark sector. The second scalar can be integrated-out, leaving a standard SHP plus a dimension-6 operator. Mainly, this model has only two relevant parameters (the DM mass and the coupling of the dim-6 operator). DM annihilation occurs mainly into two Higgs bosons, SS -> hh. We demonstrate that, despite its economy, the ESHP model provides an excellent fit to the GCE (with p-value similar to 0.6-0.7) for very reasonable values of the parameters, in particular, ms similar or equal to 130 GeV. This agreement of the DM candidate to the GCE properties does not clash with other observables and keep the S – particle relic density at the accepted value for the DM content in the universe.

Abstract: Utilizing an exhaustive set of simplified models, we revisit dark matter scenarios potentially capable of generating the observed Galactic Center gamma-ray excess, updating constraints from the LUX and PandaX- II experiments, as well as from the LHC and other colliders. We identify a variety of pseudoscalar mediated models that remain consistent with all constraints. In contrast, dark matter candidates which annihilate through a spin-1 mediator are ruled out by direct detection constraints unless the mass of the mediator is near an annihilation resonance, or the mediator has a purely vector coupling to the dark matter and a purely axial coupling to Standard Model fermions. All scenarios in which the dark matter annihilates throught-channel processes are now ruled out by a combination of the constraints from LUX/ PandaX-II and the LHC.

Abstract: Using the Fenchel-Eggleston theorem for convex hulls (an extension of the Caratheodory theorem), we prove that any likelihood can be maximized by either a dark matter 1-speed distribution F(v) in Earth's frame or 2-Galactic velocity distribution f(gal) ((u) over right arrow), consisting of a sum of delta functions. The former case applies only to time-averaged rate measurements and the maximum number of delta functions is (N-1), where N is the total number of data entries. The second case applies to any harmonic expansion coefficient of the time-dependent rate and the maximum number of terms is N. Using time-averaged rates, the aforementioned form of F(v) results in a piecewise constant unmodulated halo function (eta) over tilde (BF)-B-0 (v(min)) (which is an integral of the speed distribution) with at most (N-1) downward steps. The authors had previously proven this result for likelihoods comprised of at least one extended likelihood, and found the best-fit halo function to be unique. This uniqueness, however, cannot be guaranteed in the more general analysis applied to arbitrary likelihoods. Thus we introduce a method for determining whether there exists a unique best-fit halo function, and provide a procedure for constructing either a pointwise con fi dence band, if the best-fit halo function is unique, or a degeneracy band, if it is not. Using measurements of modulation amplitudes, the aforementioned form of f(gal) ((u) over right arrow), which is a sum of Galactic streams, yields a periodic time-dependent halo function (eta) over right arrow BF (v(min); t) which at any fixed time is a piecewise constant function of v(min) with at most N downward steps. In this case, we explain how to construct pointwise confidence and degeneracy bands from the time-averaged halo function. Finally, we show that requiring an isotropic Galactic velocity distribution leads to a Galactic speed distribution F(u)that is once again a sum of delta functions, and produces a time-dependent (eta) over tilde BF (v(min); t) function (and a time-averaged (eta) over tilde (0) BF (v(min))) that is piecewise linear, di ff ering significantly from best-fit halo functions obtained without the assumption of isotropy.