The ANtarctic Impulsive Transient Antenna (ANITA) collaboration has observed two EeV-energy, upward going events originating from below the horizon. As no standard model (SM) particles, propagating through the earth at such energy and exit angles, can give rise to the required survival probability for the observed events, several beyond standard model (BSM) proposals have come up. We propose a scenario where a $Z_2$ odd sector is responsible for such anomalous events. The next to lightest $Z_2$ odd particle or the next to lightest stable particle (NLSP), created from ultra high energy neutrino interactions with nucleons, can pass through the earth and then decay into the lightest $Z_2$ odd particle or the lightest stable particle (LSP) and a tau lepton. The long-lived nature of the NLSP requires its coupling with the LSP to be very small, keeping the LSP out of thermal equilibrium in the early universe. The LSP can then be a non-thermal dark matter while the tau leptons produced from NLSP decay after passing through the earth can explain the ANITA events. We first show that a purely non-thermal dark matter scenario can not give rise to the required event rates while a hybrid scenario where dark matter can have both thermal as a well non-thermal contribution to its relic abundance, serves the purpose.

We show that the weakly interacting massive particle (WIMP) scenario of proton-philic spin-dependent inelastic dark matter can still provide a viable explanation of the observed DAMA effect in compliance with the constraints from other experiments. We also show that, although the COSINE-100 collaboration has recently tested the DAMA effect using the same target material, for the time being the comparison between DAMA and COSINE-100 still depends on the particle-physics model.

We analyze the DAMA/Libra-phase2 modulation result using a basis of 16 effective operators describing the WIMP interaction with photons, gluons and quarks up to mass dimension seven. For each operator we fix the effective theory at the scale of 2 GeV and parametrize WIMP-quark interactions in terms of two independent couplings for up-type and down-type quarks. We discuss the connection with the non-relativistic limit of the effective theory in terms of operators invariant by Galilean transformations, and the impact of the ensuing momentum-dependent Wilson coefficients due to long-range interactions or light meson poles. Most relativistic operators yield a good fit of the DAMA modulation effect, although for parameters excluded by the constraints from XENON1T and PICO60.

Assuming for Weakly Interacting Massive Particles (WIMPs) a Maxwellian velocity distribution in the Galaxy we provide an assessment of the sensitivity of existing Dark Matter (DM) direct detection (DD) experiments to operators up to dimension 7 of the relativistic effective field theory describing dark matter interactions with quarks and gluons . In particular we focus on a systematic approach, including an extensive set of experiments and large number of couplings, both exceeding for completeness similar analyses in the literature. The relativistic effective theory requires to fix one coupling for each quark flavor, so in principle for each different combination the bounds should be recalculated starting from direct detection experimental data. To address this problem we propose an approximate model-independent procedure that allows to directly calculate the bounds for any combination of couplings in terms of model-independent limits on the Wilson coefficients of the non-relativistic theory expressed in terms of the WIMP mass and of the neutron-to-proton coupling ratio $c^n/c^p$. We test the result of the approximate procedure against that of a full calculation, and discuss its possible pitfalls and limitations. We also provide a simple interpolating interface in Python that allows to apply our method quantitatively.

Assuming a standard Maxwellian for the WIMP velocity distribution, we obtain the bounds from null WIMP search results of 59.5 days of COSINE-100 data on the DAMA/LIBRA-phase2 modulation effect within the context of the non-relativistic effective theory of WIMP-nucleus scattering. Here, we systematically assume that one of the effective operators allowed by Galilean invariance dominates in the effective Hamiltonian of a spin-1/2 dark matter (DM) particle. We find that, although DAMA/LIBRA and COSINE-100 use the same sodium-iodide target, the comparison of the two results still depends on the particle-physics model. This is mainly due to two reasons: i) the WIMP signal spectral shape; ii) the expected modulation fractions, when the upper bound on the time-averaged rate in COSINE-100 is converted into a constraint on the annual modulation component in DAMA/LIBRA. We find that the latter effect is the dominant one. For several effective operators the expected modulation fractions are larger than in the standard spin-independent or spin-dependent interaction cases. As a consequence, compatibility between the modulation effect observed in DAMA/LIBRA and the null result from COSINE-100 is still possible for several non-relativistic operators. At low WIMP masses such relatively high values of the modulation fractions arise because COSINE-100 is mainly sensitive to WIMP-sodium scattering events, due to the higher threshold compared to DAMA/LIBRA. A next COSINE analysis is expected to have a full sensitivity for the 5$\sigma$ region of DAMA/LIBRA.

Assuming for Weakly Interacting Massive Particles (WIMPs) a Maxwellian velocity distribution in the Galaxy we explore in a systematic way the relative sensitivity of an extensive set of existing and projected Dark Matter (DM) direct detection experiments to each of the 14 couplings that parameterize the most general non-relativistic (NR) effective Hamiltonian allowed by Galilean invariance for a contact interaction driving the elastic scattering off nuclei of WIMPs of spin 1/2. We perform our analysis in terms of two free parameters: the WIMP mass $m_{\chi}$ and the ratio between the WIMP-neutron and the WIMP-proton couplings $c^n/c^p$. We include the modified signal spectral shape due to non-standard interactions when it is needed in the determination of the bound, such as in the case of background subtraction or of the application of the optimal-interval method. For each coupling, in the $m_{\chi}$-$c^n/c^p$ plane we provide contour plots of the most stringent 90 % C.L. bound on the WIMP-nucleon cross section and show the experiment providing it. We also introduce NRDD_constraints, a simple interpolating code written in Python that allows to obtain the numerical value of the bound as a function of the WIMP mass $m_{\chi}$ and of the coupling ratio $c^n/c^p$ for each NR coupling. We find that 9 experiments out of the 14 present Dark Matter searches considered in our analysis provide the most stringent bound on some of the effective couplings for a given choice of $(m_{\chi},c^n/c^p)$: this is evidence of the complementarity of different target nuclei and/or different combinations of count-rates and energy thresholds when the search of DM is extended to a wide range of possible interactions.

We show that the Weakly Interacting Massive Particle scenario of proton-philic spin-dependent inelastic Dark Matter (pSIDM) can still provide a viable explanation of the observed DAMA modulation amplitude in compliance with the constraints from other experiments after the release of the DAMA/LIBRA-phase2 data and including the recent bound from COSINE-100, that uses the same $NaI$ target of DAMA. The pSIDM scenario provided a viable explanation of DAMA/LIBRA--phase1 both for a Maxwellian WIMP velocity distribution and in a halo-independent approach. At variance with DAMA/LIBRA-phase1, for which the modulation amplitudes showed an isolated maximum at low energy, the DAMA/LIBRA-phase2 spectrum is compatible to a monotonically decreasing one. Moreover, due to its lower threshold, it is sensitive to WIMP-iodine interactions at low WIMP masses. Due to the combination of these two effects pSIDM can now explain the yearly modulation observed by DAMA/LIBRA only when the WIMP velocity distribution departs from a standard Maxwellian. In this case the WIMP mass $m_{\chi}$ and mass splitting $\delta$ fall in the approximate ranges 7 GeV $\lesssim m_{\chi}\lesssim$ 17 GeV and 18 keV$\lesssim\delta\lesssim$29 keV. The recent COSINE-100 bound is naturally evaded in the pSDIM scenario due to its large expected modulation fractions.

We re-examine the case of anapole dark matter as an explanation for the DAMA annual modulation in light of the DAMA/LIBRA-phase2 results and improved upper limits from other DM searches. If the WIMP velocity distribution is assumed to be a Maxwellian, anapole dark matter is unable to provide an explanation of the DAMA modulation compatible with the other searches. Nevertheless, anapole dark matter provides a better fit to the DAMA-phase2 modulation data than an isoscalar spin-independent interaction, due to its magnetic coupling with sodium targets. A halo-independent analysis shows that explaining the DAMA modulation above 2 keVee in terms of anapole dark matter is basically impossible in face of the other null results, while the DAMA/LIBRA-phase2 modulation measurements below 2 keVee are marginally allowed. We conclude that in light of current measurements, anapole dark matter does not seem to be a viable explanation for the totality of the DAMA modulation.

The DAMA/LIBRA collaboration has recently released updated results from their search for the annual modulation signal expected from Dark Matter (DM) scattering in their NaI detectors. We have fitted the updated DAMA result for the modulation amplitudes in terms of a Weakly Interacting Massive Particle (WIMP) signal, parameterizing the interaction with nuclei in terms of the most general effective Lagrangian for a WIMP particle spin up to 1/2, systematically assuming dominance of one of the 14 possible interaction terms, and assuming for the WIMP velocity distribution a standard Maxwellian. We find that most of the couplings of the non-relativistic effective Hamiltonian can provide a better fit compared to the standard Spin Independent interaction case, and with a reduced fine-tuning of the three parameters (WIMP mass, WIMP-nucleon effective cross-section and ratio between the WIMP-neutron and the WIMP-proton couplings). Moreover, effective models for which the cross section depends explicitly on the WIMP incoming velocity can provide a better fit of the DAMA data at large values of $m_{\chi}$ compared to the standard velocity-independent cross-section due to a different phase of the modulation amplitudes. All the best fit solutions are in tension with exclusion plots of both XENON1T and PICO60.

We consider a simple extension of the minimal left-right symmetric model (LRSM) in order to explain the PeV neutrino events seen at the IceCube experiment from a heavy decaying dark matter. The dark matter sector is composed of two fermions: one at PeV scale and the other at TeV scale such that the heavier one can decay into the lighter one and two neutrinos. The gauge annihilation cross sections of PeV dark matter are not large enough to generate its relic abundance within the observed limit. We include a pair of real scalar triplets $\Omega_{L,R}$ which can bring the thermally overproduced PeV dark matter abundance into the observed range through late time decay and consequent entropy release thereby providing a consistent way to obtain the correct relic abundance without violating the unitarity bound on dark matter mass. Another scalar field, a bitriplet under left-right gauge group is added to assist the heavier dark matter decay. The presence of an approximate global $U(1)_X$ symmetry can naturally explain the origin of tiny couplings required for long-lived nature of these decaying particles. We also show, how such an extended LRSM can be incorporated within a non-supersymmetric $SO(10)$ model where the gauge coupling unification at a very high scale naturally accommodate a PeV scale intermediate symmetry, required to explain the PeV events at IceCube.

We consider scenarios in which the annihilation of self-conjugate spin-1 dark matter to a Standard Model fermion-antifermion final state is chirality suppressed, but where this suppression can be lifted by the emission of an additional photon via internal bremsstrahlung. We find that this scenario can only arise if the initial dark matter state is polarized, which can occur in the context of self-interacting dark matter. In particular, this is possible if the dark matter pair forms a bound state that decays to its ground state before the constituents annihilate. We show that the shape of the resulting photon spectrum is the same as for self-conjugate spin-0 and spin-1/2 dark matter, but the normalization is less heavily suppressed in the limit of heavy mediators.

In a Dark left-right gauge model, the neutral component of right-handed lepton doublet is odd under generalized R-parity and thus the lightest one serves as the dark matter (DM) candidate. The coannihilation of the dark matter with the singly charged Higgs triplet produces the correct relic abundance. We explain AMS-02 positron excess by the annihilation of 800 GeV dark matter into $\mu^+\mu^-\gamma$, through a t-channel exchange of the additional charged triplet Higgs boson. The DM is leptophilic which is useful for explaining the non-observation of any antiproton excess which would generically be expected from DM annihilation. The large cross-section needed to explain AMS-02 also requires an astrophysical boost. In addition, we show that the muon $(g-2)$ receives required contribution from singly and doubly charged triplet Higgs in the loops.