Physics Letters B 786 (2018) 337–341
Contents lists available at ScienceDirect
Physics Letters B www.elsevier.com/locate/physletb
XENON1T takes a razor to a dark E 6 -inspired model Daniel A. Camargo a , Yann Mambrini b , Farinaldo S. Queiroz a,∗ a b
International Institute of Physics, Universidade Federal do Rio Grande do Norte, Campus Universitário, Lagoa Nova, Natal-RN 59078-970, Brazil Laboratoire de Physique Théorique (UMR8627), CNRS, Univ. Paris-Sud, Université Paris-Saclay, 91405 Orsay, France
a r t i c l e
i n f o
Article history: Received 21 June 2018 Received in revised form 6 August 2018 Accepted 7 September 2018 Available online 11 October 2018 Editor: A. Ringwald
a b s t r a c t Motivated by the recent result of XENON1T collaboration with full exposure, 279 life days, that sets the most stringent limit on the spin-independent dark matter-nucleon scattering cross section we discuss a dark E 6 -inspired model that features the presence of a U (1)d−u gauge symmetry. The dark matter candidate is a Dirac fermion that interacts with Standard Model fermions via a massive Z that preserves the quantum number assignments of this symmetry. We compute the spin-independent scattering cross section off xenon nucleus and compare with the XENON1T limit; find the LHC bound on the Z mass as well as the projection sensitivity of high-energy and luminosity LHC; and derive the Fermi-LAT bounds on the dark matter annihilation cross section based on the observation of gamma-rays in the direction of Dwarf Spheroidal galaxies. We exploit the complementarity between these datasets to conclude that the new bound from XENON1T severely constrain the model, which combined with the LHC upgrade sensitivity rules out this WIMP realization setup below 5 TeV. © 2018 Published by Elsevier B.V. This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/4.0/). Funded by SCOAP3 .
1. Introduction The presence of dark matter (DM) has been established by numerous observations stemming from different scales and epochs in the history of the universe [1]. Apart from gravitation effects, we have not been able to detect the interactions involving dark matter particles and for this reason its nature remains unknown [2,3]. There are several dark matter candidates in the literature with masses spanning over several orders of magnitude [4,5]. However, massive particles that are thermally in the universe are theoretically well motivated the most popular (massive particles that are thermally produced stays the most popular) [4,6,7]. They are commonly referred to as WIMPs (Weakly Interacting Massive Particles) [8]. WIMPs are subject of intense searches that involved different techniques and observables which are classified as direct, indirect and collider probes [9,10]. Direct detection refers to the measurement of the recoil energy deposited in nuclear scatterings by WIMP-like particles on their scattering off of nuclei target used in underground experiments [11]. Using discriminating variable to distinguish nuclear recoil from electron recoil and a good control over background from various sources, direct detection experiments become discover machines since the scattering rate depends
*
Corresponding author. E-mail address:
[email protected] (F.S. Queiroz).
on the properties of dark matter particle, such as the scattering cross section and mass [12–15], and the local dark matter density [16–18]. Therefore, null results from direct detection experiments can be interpreted as upper limits on the pair WIMP-nucleon scattering cross section and dark matter mass [19]. We have several experiments searching for WIMP events using different readout techniques and targets that make them complementary to each other and sensitive to different regions of the scattering cross section vs dark matter mass parameter space [20,21]. Indirect detection covers to the observation of gamma-rays or other possible final states resulted from DM annihilation or decay [22–32]. Fermi-LAT has been continuously observing the gammaray emission in the direction of over a dozen of Dwarf Spheroidal galaxies which are dominated by dark matter, and since no significant excess has been observed, restrictive limits have been placed on the gamma-ray flux resulted from dark matter annihilations [33,34]. The gamma-ray flux originated from dark matter annihilation depends on the dark matter density distribution known as density profile; the gamma-ray energy spectrum that accounts for how many gamma-rays are produced per annihilation; on the annihilation cross section and dark matter mass [35–38]. The gamma-ray flux produced by dark matter annihilation is very important especially for dark matter models that rely entirely on the thermal production of dark matter and the dark matter annihilation cross section is velocity independent as it happens to be the case for our dark E 6 -inspired model [39–43].
https://doi.org/10.1016/j.physletb.2018.09.057 0370-2693/© 2018 Published by Elsevier B.V. This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/4.0/). Funded by SCOAP3 .
338
D.A. Camargo et al. / Physics Letters B 786 (2018) 337–341
Collider searches for WIMPs are characterized by events featuring large missing energy accompanied by visible counter-parts [44–47]. Sometimes the most efficient method to probe a dark sector is not by searching for the dark matter itself at colliders, especially when it has mediators involved have sizable couplings to SM fermions. In the model we consider, the massive Z does have sizable couplings to leptons and thus it is subject to stringent bounds stemming from the LHC [48–50]. With these dark matter searches in mind, a multitude of phenomenological studies have been carried out to illustrate the viable parameter based on simplified dark matter models without specifying a gauge group. One of the most attractive is the Dirac fermion charged under an U (1) gauge symmetry. In this context, the Z gauge boson is a natural mediator between the dark and the visible sector. We advocate that it is more interesting to fix a gauge setup which is motivated by theoretical constructions and assess whether a viable dark matter candidate is possible without worrying about the further complications arising from a UV completion [51–62]. That said, we choose to explore dark matter complementarity in the context of collider, direct and indirect detection searches in a specific setup where the dark matter particle is a Dirac fermion that interacts with a Z boson that possesses interactions with the SM particle as dictated by the U (1)d−u gauge structure which appears in grand unification studies of the E 6 group [63]. The paper is organized as follows: In section 1 we introduce the model; in section 2 we discuss the dark matter observables; in section 3 we exhibit our findings before concluding in section 4. 2. Model The origin of a Z gauge boson is behind a new U (1) gauge symmetry. The charges of the SM fermions under the U (1) gauge symmetry are in principle arbitrary but should respect SU (2) invariance with the left-handed up and down quarks having the same quantum number [64,65]. The same logic applied to the lefthanded charged leptons and their corresponding neutrino flavors. Since the new interaction is not asymptotically free in a UV completion, the model should be well-behaved at high energies and thus anomaly-free [66,67]. The anomaly cancellation can be highly non-trivial and for this reason it is theoretically interesting to consider gauge groups that are generation-independent such as the U (1)d−u group that arises after spontaneous symmetry breaking in some E 6 constructions [63]. The E 6 symmetry is theoretically motivated because it offers a hospital environment for unification of forces and incorporates interesting features of SU (5) or S O (10) groups concerning the generation of fermion masses and flavor symmetries [68–74], super heavy dark matter [75], and collider physics [72]. In order to incorporate dark matter we introduce a Dirac particle we write down the simplified Lagrangian that describes the dark matter interactions with the Z and SM fermion as,
μ ¯f γ μ g f v + g f a γ 5 f , ¯ L ⊃ Zμ χγ gχ v + gχ a γ 5 χ + f ∈SM
(1) where the sum is over all quarks and leptons. The Z -SM couplings in Eq. (1) are determined by the SM fermion charges under the new gauge symmetry as follows,
gZ zu R + z Q L , g ua = zu R − z Q L 2 2 gZ gZ = zd R + z Q L , gda = zd R − z Q L 2 2
g uv = gdv
gZ
Table 1 Standard Model fermions U (1)d−u charges.
glv =
QL
uR
dR
lL
eR
0
−1/3
1/3
0
1/3
gZ 2
=
gν v
ze R + zl L
gZ 2
, gla =
gZ 2
, gνa = −
zl L
ze R − zl L
gZ 2
zl L .
(2)
The z’s charges of the fields as shown in Table 1. In our model (Eq. (1)) we are assuming that the dark matter candidate is a chiral field. Therefore the gauge anomaly cancellation becomes nontrivial and the mechanism behind the anomaly cancellation and Z mass can change our numerical findings [64,67,76–87]. We will assume that Eq. (1) offers a good description of the dark matter phenomenology and to simplify our analysis we adopt g χ v = g χ a . The scenarios with g χ v = g χ a and the g χ a = 0 share a similar dark matter phenomenology and this setup with g χ a = 0 is equivalent to the one of vector-like dark matter fermion. In other words, the dark matter observables and collider bounds derived here are similar to the one present in vector-like fermion dark matter setups. Moreover, we adopt g Z = 1 throughout. That said, Eq. (1) will dictate the entire phenomenology of our work as we describe hereafter. We point out that E 6 group breaks into SU (5) ⊗ U (1) ⊗ U (1). Depending on how the spontaneous symmetry is broken you can generate different U (1) models. The most common breaking E 6 → SU (5) ⊗ U (1)χ ⊗ U (1)η , as adopted by the CMS and ATLAS collab cos θ + Z sin θ , orations, where the physical Z is given by Z = Z ψ χ where θ is the mixing angle which can be in principle suppressed. Our model U (1)d−u resembles the Z boson of the U (1)ψ model because the magnitude of Z couplings to right-handed fermions is the same, but the U (1)d−u model does not have couplings to lepton-handed fields, and for this reason the U (1)d−u model is also known as the right-handed U (1). Anyways, the dark matter phenomenology is dominated by the Z neutral current, thus we do not expect significant differences between our model and the U (1)ψ model. 2.1. Direct detection The Lagrangian with a dark matter field coupling chirally via a Z such as the one presented in Eq. (1) leads to a dark matter-nucleus scattering that is both spin-independent and spindependent. The spin-independent process is coherent, so it scales with the atomic mass of the nucleus. On the other hand, the spindependent scattering goes with the net spin of the nucleus which is not much far from one. Therefore, bounds on spin-independent scattering are much more stringent for heavy nuclei such as xenon that has an atomic mass A = 131. For this reason, we will concentrate our direct detection discussion on spin-independent scattering. The WIMP-nucleon spin-independent scattering cross section in this case is,
σnS I =
μ2χ n π
Z f p + ( A − Z ) fn A
2 ,
(3)
where,
f p =≡
gχ v m2Z
(2g uv + gdv )
(4)
and
fn ≡
gχ v m2Z
( g uv + 2gdv ) ,
(5)
D.A. Camargo et al. / Physics Letters B 786 (2018) 337–341
Fig. 1. Dark matter complementarity plot for mχ = 50 GeV. In blue we show the excluded region based on the recent XENON1T bound on the spin-independent scattering cross section. In green, the region ruled out by Fermi-LAT through the observation of dwarf spheroidal galaxies. The red horizontal curve is the LHC limit on the Z mass using dilepton data at 13 TeV. The black dashed curve delimits the region of parameter space that yields the correct relic density.
339
Fig. 3. Same than Fig. 1 but for mχ = 5 TeV.
2.2. Indirect detection The relevant observable as far as indirect dark matter detection is the differential gamma-ray flux which depends on the dark matter annihilation cross section, the dark matter mass, and the final states produced in the annihilation. There are two possible annihilation modes. One where the Dirac fermion annihilates into SM fermion pairs and another where Z pairs are produced, the latter being relevant only when mχ ≥ m Z . The dark matter particles are slow-moving, thus this annihilation process is non-relativistic. That said, annihilation cross sections are
σ v χ χ¯ → f ¯f
nc
≈
1 2
https://indico.cern.ch/event/726877/. There are competitive limits from other experiments [20,89–91].
m2f m2χ
2π m4Z 4m2χ − m2Z
2
2 × g 2f a 2g χ2 vm4Z m2χ − m2f + g χ2 a m2f 4m2χ − m2Z
, (6) + g χ2 v g 2f vm4Z 2m2χ + m2f
Fig. 2. Same than Fig. 1 but for mχ = 1 TeV.
with μχ n being the WIMP-nucleon reduced mass, Z the atomic number of the target nucleus. With the theoretical prediction at hand, we can compare our finding with the current limit from XENON1T collaboration presented at MPIK on May 28, 20181 which was later followed by the XENON1T publication [88]. The experiment limit is based on 278.8 life days exposure with 1.3 tons fiducial mass. The expected background events are 0.75 neutrons, 5.36 surface events, 0.02 from neutrino-nucleus coherent background among other that total about 8.4 ± 0.63 events. The collaboration observed about 11 events in the signal region thus with a small excess of events. Because of the presence of such excess of events the scientific goal of the collaboration of achieving a sensitivity to the spinindependent scattering cross section down to 1 × 10−47 cm2 was not achieved but the collaboration successfully carried their analysis and set the world leading upper limit on the WIMP-nucleon spin-independent scattering cross section of 4 × 10−47 cm2 [88]. These results severely constrain our model as shown in Figs. 1–3. The blue hashed region represents the XENON1T exclusion region.2
1−
and
σ v χ χ¯ → Z Z ≈
1 16π m2χ m2Z
1−
m2Z m2χ
3/2 1−
m2Z
− 2
2m2χ
× 8g χ2 v g χ2 a m2χ + g χ4 v + g χ4 a − 6g χ2 v g χ2 a m2Z ,
(7)
where v is the relative velocity of the incoming dark matter particles, nc is the color number of the final states SM fermions. The Eq. (6) is not valid near the Z resonance, i.e. when m Z ∼ 2mχ . In this case, the Z width should be properly included. We handled this procedure numerically within the MicrOmegas package [92–94]. With the results we found the region of parameter space that yields that correct dark matter relic density of h2 = 0.12 according to Planck [95]. The region of parameter that furnishes
h2 = 0.12 is displayed in the Figs. 1–3 with a dashed black curve. Concerning indirect detection, there is another quantity to be determined which is the gamma-ray yield. To compute it we implemented the model in Pythia 8 [96] and took into account all relevant annihilation channels. With the annihilation cross section
340
D.A. Camargo et al. / Physics Letters B 786 (2018) 337–341
and energy spectrum at hand, we can compute the dark matter differential flux and compare with upper limits provided by Fermi-LAT [33]. Using this procedure we drew Fermi-LAT exclusion regions displayed in green in Figs. 1–3. There are other competitive limits on the dark matter annihilation cross section but we adopted the ones stemming from dwarf spheroidal galaxies for being more robust [24,97,98]. Since the annihilation cross section is s-wave, the indirect detection limits from Fermi-LAT trace the relic density curve as seen in the figures. We highlight that there are bounds resulted from the production of gamma-ray lines at one-loop level but these bounds are not competitive for vector mediators because of the loop suppression and the relative weak experimental sensitivity to gamma-ray lines [99,100].
FSQ and DC acknowledge support from MEC and UFRN. FSQ also thanks the financial support from ICTP-SAIFR FAPESP grant 2016/01343-7. Y.M. acknowledges partial support from the European Union Horizon 2020 research and innovation program under the Marie Sklodowska-Curie: RISE InvisiblesPlus (grant agreement No. 690575), the ITN Elusives (grant agreement No. 674896), and the France-US PICS MicroDark (grant agreement No. 690575).
3. Collider
References
The collider bounds we derived are based on the dilepton data provided ATLAS collaboration using 13 TeV center-of-mass energy with about 36 fb−1 of data [101]. Dilepton resonance searches are effective laboratories to constrain the properties of Z gauge bosons that possess sizable coupling to charged leptons. The couplings to quarks are also sizable but dilepton searches provide stronger bounds simply because the signal suffers from a smaller SM background. We implemented the model in Madgraph and followed the receipt explained in detail in [102] to find the lower mass bound m Z > 4068 GeV. This bound is applicable when the invisible Z decay into dark matter is not accessible. When m Z > 2mχ this invisible decay opens up diminishing the Z branching ratio into charged lepton and thus weakening the LHC sensitivity. This feature is visible in the Figs. 1 and 2 where the limit on Z mass from LHC searches becomes weaker once g χ = g χ v = g χ a reaches a reasonable value ( 0.1 − 1). On the other hand, this is clearly not visible in Fig. 3 where the dark matter mass (5 TeV) always lies above the dilepton limit on Z mass. 4. Discussion and conclusions In Figs. 1–3 we collect all these bounds to find the allowed region of parameter space. It is quite clear how important the recent XENON1T bound is to probe this model. XENON1T excludes alone most of the parameter of the model that yields the correct relic density even for dark matter masses up to 5 TeV. The indirect detection limits are subdominant in the case of vector mediator simply because of the stringent bounds stemming from LHC on the gauge boson mass. The model only survives for dark matter masses about ∼5 TeV. This happens because the Z resonance that drives the relic density for m Z ∼ 2mχ takes place outside the LHC sensitivity. In summary, the combination of current LHC and XENON1T data pushes the dark E 6 model to live at the multi TeV window. Moreover, we also derived the high-energy LHC (HE-LHC) sensitivity using the code provided in [103] with high-luminosity. The HE-LHC reach is based on a 27 TeV proton-proton collision with an integrated luminosity of 3000 fb−1 . In case of null results, this configuration can probe Z masses up to 12 TeV. This forecast is exhibited with a dashed red line in Figs. 1–3. Considering the prospects for the HE-LHC, we conclude from Figs. 1–3 that the LHC upgrade is paramount to nearly close the window for a Dirac fermion dark matter in the E 6 -inspired model described here for dark matter masses below 5 TeV. Notice that there is tiny region in Fig. 3 that still survives after the HE-LHC prospect is included. We emphasize that this conclusion relies on the fact that the simplified Lagrangian in Eq. (1) provides a good description of the dark matter phenomenology of the dark E 6 model and that g Z = 1. A smaller g Z coupling would certainly weaken the collider and the direct detection bounds, however a
larger g χ would be needed to still reproduce the correct dark matter relic density. At the end, the overall conclusions would remain the same. Acknowledgements
[1] G. Bertone, D. Hooper, J. Silk, Phys. Rep. 405 (2005) 279, arXiv:hep-ph/ 0404175 [hep-ph]. [2] F.S. Queiroz, W. Rodejohann, C.E. Yaguna, Phys. Rev. D 95 (2017) 095010, arXiv:1610.06581 [hep-ph]. [3] B.J. Kavanagh, F.S. Queiroz, W. Rodejohann, C.E. Yaguna, arXiv:1706.07819 [hep-ph], 2017. [4] L. Bergstrom, New J. Phys. 11 (2009) 105006, arXiv:0903.4849 [hep-ph]. [5] M. Klasen, M. Pohl, G. Sigl, Prog. Part. Nucl. Phys. 85 (2015) 1, arXiv:1507. 03800 [hep-ph]. [6] G. Jungman, M. Kamionkowski, K. Griest, Phys. Rep. 267 (1996) 195, arXiv: hep-ph/9506380 [hep-ph]. [7] H. Baer, K.-Y. Choi, J.E. Kim, L. Roszkowski, Phys. Rep. 555 (2015) 1, arXiv: 1407.0017 [hep-ph]. [8] G. Arcadi, M. Dutra, P. Ghosh, M. Lindner, Y. Mambrini, M. Pierre, S. Profumo, F.S. Queiroz, arXiv:1703.07364 [hep-ph], 2017. [9] P. Cushman, et al., in: Community Summer Study 2013: Snowmass on the Mississippi, CSS2013, Minneapolis, MN, USA, July 29-August 6, 2013, 2013, arXiv:1310.8327 [hep-ex]. [10] F.S. Queiroz, in: Proceedings, 51st Rencontres de Moriond on Electroweak Interactions and Unified Theories: La Thuile Italy, March 12–19, 2016, ARISF, ARISF 2016, 2016, pp. 427–436, arXiv:1605.08788 [hep-ph]. [11] T. Marrodan Undagoitia, L. Rauch, J. Phys. G 43 (2016) 013001, arXiv:1509. 08767 [physics.ins-det]. [12] P. Gondolo, in: Dark Matter in Cosmology, Quantum Measurements, Experimental Gravitation. Proceedings, 31st Rencontres de Moriond, 16th Moriond Workshop, Les Arcs, France, January 2–27, 1996, 1996, pp. 41–51, arXiv: hep-ph/9605290 [hep-ph]. [13] M. Pato, L. Baudis, G. Bertone, R. Ruiz de Austri, L.E. Strigari, R. Trotta, Phys. Rev. D 83 (2011) 083505, arXiv:1012.3458 [astro-ph.CO]. [14] K. Freese, M. Lisanti, C. Savage, Rev. Mod. Phys. 85 (2013) 1561, arXiv:1209. 3339 [astro-ph.CO]. [15] E. Del Nobile, G.B. Gelmini, P. Gondolo, J.-H. Huh, in: Proceedings, 13th International Conference on Topics in Astroparticle and Underground Physics, TAUP 2013, Asilomar, California, September 8–13, 2013, Phys. Proc. 61 (2015) 45, arXiv:1405.5582 [hep-ph]. [16] R. Catena, P. Ullio, J. Cosmol. Astropart. Phys. 1008 (2010) 004, arXiv:0907. 0018 [astro-ph.CO]. [17] M. Pato, O. Agertz, G. Bertone, B. Moore, R. Teyssier, Phys. Rev. D 82 (2010) 023531, arXiv:1006.1322 [astro-ph.HE]. [18] F. Iocco, M. Pato, G. Bertone, P. Jetzer, J. Cosmol. Astropart. Phys. 1111 (2011) 029, arXiv:1107.5810 [astro-ph.GA]. [19] C. Kelso, D. Hooper, M.R. Buckley, Phys. Rev. D 85 (2012) 043515, arXiv:1110. 5338 [astro-ph.CO]. [20] D.S. Akerib, et al., LUX, Phys. Rev. Lett. 112 (2014) 091303, arXiv:1310.8214 [astro-ph.CO]. [21] R. Agnese, et al., SuperCDMS, Phys. Rev. Lett. 112 (2014) 041302, arXiv:1309. 3259 [physics.ins-det]. [22] P. Scott, J. Conrad, J. Edsjo, L. Bergstrom, C. Farnier, Y. Akrami, J. Cosmol. Astropart. Phys. 1001 (2010) 031, arXiv:0909.3300 [astro-ph.CO]. [23] A.A. Abdo, et al., Fermi-LAT, Astrophys. J. 712 (2010) 147, arXiv:1001.4531 [astro-ph.CO]. [24] D. Hooper, C. Kelso, F.S. Queiroz, Astropart. Phys. 46 (2013) 55, arXiv:1209. 3015 [astro-ph.HE]. [25] A.X. Gonzalez-Morales, S. Profumo, F.S. Queiroz, Phys. Rev. D 90 (2014) 103508, arXiv:1406.2424 [astro-ph.HE]. [26] Y. Mambrini, S. Profumo, F.S. Queiroz, Phys. Lett. B 760 (2016) 807, arXiv: 1508.06635 [hep-ph]. [27] M.G. Baring, T. Ghosh, F.S. Queiroz, K. Sinha, arXiv:1510.00389 [hep-ph], 2015. [28] S. Li, Y.-F. Liang, K.-K. Duan, Z.-Q. Shen, X. Huang, X. Li, Y.-Z. Fan, N.-H. Liao, L. Feng, J. Chang, Phys. Rev. D 93 (2016) 043518, arXiv:1511.09252 [astro-ph. HE].
D.A. Camargo et al. / Physics Letters B 786 (2018) 337–341
[29] A. Drlica-Wagner, et al., DES, Fermi-LAT, Astrophys. J. 809 (2015) L4, arXiv: 1503.02632 [astro-ph.HE]. [30] Y. Zhao, X.-J. Bi, H.-Y. Jia, P.-F. Yin, F.-R. Zhu, Phys. Rev. D 93 (2016) 083513, arXiv:1601.02181 [astro-ph.HE]. [31] C. Garcia-Cely, J. Heeck, J. Cosmol. Astropart. Phys. 1608 (2016) 023, arXiv: 1605.08049 [hep-ph]. [32] S. Li, et al., arXiv:1805.06612 [astro-ph.HE], 2018. [33] M. Ackermann, et al., Fermi-LAT, Phys. Rev. Lett. 115 (2015) 231301, arXiv: 1503.02641 [astro-ph.HE]. [34] M.L. Ahnen, et al., Fermi-LAT, MAGIC, J. Cosmol. Astropart. Phys. 1602 (2016) 039, arXiv:1601.06590 [astro-ph.HE]. [35] S. Profumo, in: Proceedings, Theoretical Advanced Study Institute in Elementary Particle Physics: Searching for New Physics at Small and Large Scales, TASI 2012, Boulder, Colorado, June 4–29, 2012, 2013, pp. 143–189, arXiv: 1301.0952 [hep-ph]. [36] A. Lopez, C. Savage, D. Spolyar, D.Q. Adams, in: Proceedings, Meeting of the APS Division of Particles and Fields, DPF 2015, Ann Arbor, Michigan, USA, 4–8 Aug. 2015, J. Cosmol. Astropart. Phys. 1603 (2016) 033, arXiv:1501.01618 [astro-ph.CO]. [37] M. Tavani, et al., e-ASTROGAM, arXiv:1711.01265 [astro-ph.HE], 2017. [38] B.S. Acharya, et al., Cherenkov Telescope Array Consortium, arXiv:1709.07997 [astro-ph.IM], 2017. [39] T.R. Slatyer, Phys. Rev. D 93 (2016) 023521, arXiv:1506.03812 [astro-ph.CO]. [40] T.R. Slatyer, C.-L. Wu, Phys. Rev. D 95 (2017) 023010, arXiv:1610.06933 [astroph.CO]. [41] H. Liu, T.R. Slatyer, J. Zavala, Phys. Rev. D 94 (2016) 063507, arXiv:1604.02457 [astro-ph.CO]. [42] N. Bernal, L. Necib, T.R. Slatyer, J. Cosmol. Astropart. Phys. 1612 (2016) 030, arXiv:1606.00433 [astro-ph.CO]. [43] M. Dutra, M. Lindner, S. Profumo, F.S. Queiroz, W. Rodejohann, C. Siqueira, J. Cosmol. Astropart. Phys. 1803 (2018) 037, arXiv:1801.05447 [hep-ph]. [44] J. Goodman, M. Ibe, A. Rajaraman, W. Shepherd, T.M. Tait, et al., Phys. Rev. D 82 (2010) 116010, arXiv:1008.1783 [hep-ph]. [45] A. Rajaraman, W. Shepherd, T.M.P. Tait, A.M. Wijangco, Phys. Rev. D 84 (2011) 095013, arXiv:1108.1196 [hep-ph]. [46] D. Bauer, et al., Snowmass 2013 Cosmic Frontier Working Groups 1-4, Phys. Dark Universe 7–8 (16) (2015), arXiv:1305.1605 [hep-ph]. [47] B. Szczerbinska, R. Allahverdi, K. Babu, B. Balantekin, B. Dutta, T. Kamon, J. Kumar, F. Queiroz, L. Strigari, R. Surman (Eds.), Proceedings, Workshop on Neutrino Physics: Session of CETUP* 2015 and 9th International Conference on Interconnections between Particle Physics and Cosmology, vol. 1743, PPC2015, 2016. [48] G. Aad, et al., ATLAS, Phys. Rev. D 90 (2014) 052005, arXiv:1405.4123 [hepex]. [49] V. Khachatryan, et al., CMS, Phys. Lett. B 768 (2017) 57, arXiv:1609.05391 [hep-ex]. [50] A.M. Sirunyan, et al., CMS, arXiv:1803.06292 [hep-ex], 2018. [51] J.L. Lopez, D.V. Nanopoulos, K.-j. Yuan, Nucl. Phys. B 370 (1992) 445. [52] H. Baer, S. Kraml, S. Sekmen, H. Summy, J. High Energy Phys. 03 (2008) 056, arXiv:0801.1831 [hep-ph]. [53] H. Baer, H. Summy, Phys. Lett. B 666 (2008) 5, arXiv:0803.0510 [hep-ph]. [54] K. Agashe, K. Blum, S.J. Lee, G. Perez, Phys. Rev. D 81 (2010) 075012, arXiv: 0912.3070 [hep-ph]. [55] J. Chakrabortty, S. Mohanty, S. Rao, J. High Energy Phys. 02 (2014) 074, arXiv: 1310.3620 [hep-ph]. [56] A. Anandakrishnan, K. Sinha, Phys. Rev. D 89 (2014) 055015, arXiv:1310.7579 [hep-ph]. [57] P. Ko, Y. Omura, C. Yu, J. High Energy Phys. 06 (2015) 034, arXiv:1502.00262 [hep-ph]. [58] C. Arbelaez, R. Longas, D. Restrepo, O. Zapata, Phys. Rev. D 93 (2016) 013012, arXiv:1509.06313 [hep-ph]. [59] H.-Y. Chen, I. Gogoladze, S. Hu, T. Li, L. Wu, Eur. Phys. J. C 78 (2018) 26, arXiv: 1703.07542 [hep-ph]. [60] G. Arcadi, M. Lindner, Y. Mambrini, M. Pierre, F.S. Queiroz, Phys. Lett. B 771 (2017) 508, arXiv:1704.02328 [hep-ph]. [61] C. Coriano, P.H. Frampton, arXiv:1712.03865 [hep-ph], 2017. [62] S.J.D. King, S.F. King, S. Moretti, arXiv:1712.01279 [hep-ph], 2017. [63] C. Patrignani, et al., Particle Data Group, Chin. Phys. C 40 (2016) 100001. [64] N.F. Bell, G. Busoni, I.W. Sanderson, J. Cosmol. Astropart. Phys. 1703 (2017) 015, arXiv:1612.03475 [hep-ph]. [65] N. Bell, G. Busoni, A. Kobakhidze, D.M. Long, M.A. Schmidt, J. High Energy Phys. 08 (2016) 125, arXiv:1606.02722 [hep-ph].
341
[66] A. Ismail, W.-Y. Keung, K.-H. Tsao, J. Unwin, Nucl. Phys. B 918 (2017) 220, arXiv:1609.02188 [hep-ph]. [67] A. Ismail, A. Katz, D. Racco, J. High Energy Phys. 10 (2017) 165, arXiv:1707. 00709 [hep-ph]. [68] G.-C. Cho, K. Hagiwara, Y. Umeda, Nucl. Phys. B 531 (1998) 65, Erratum: Nucl. Phys. B 565 (2000) 483, arXiv:hep-ph/9805448 [hep-ph]. [69] S. Nie, M. Sher, Phys. Rev. D 64 (2001) 073015, arXiv:hep-ph/0102139 [hepph]. [70] R. Howl, S.F. King, J. High Energy Phys. 01 (2008) 030, arXiv:0708.1451 [hepph]. [71] B. Stech, Z. Tavartkiladze, Phys. Rev. D 77 (2008) 076009, arXiv:0802.0894 [hep-ph]. [72] E. Rojas, J. Erler, J. High Energy Phys. 10 (2015) 063, arXiv:1505.03208 [hepph]. [73] S. Benli, T. Dereli, https://doi.org/10.1007/s10773-018-3757-8, arXiv:1707. 03144 [hep-ph], 2017. [74] B. Dutta, S. Ghosh, I. Gogoladze, T. Li, arXiv:1805.01866 [hep-ph], 2018. [75] J. Schwichtenberg, J. High Energy Phys. 02 (2018) 016, arXiv:1704.04219 [hepph]. [76] F. Kahlhoefer, K. Schmidt-Hoberg, T. Schwetz, S. Vogl, J. High Energy Phys. 02 (2016) 016, arXiv:1510.02110 [hep-ph]. [77] A. Celis, W.-Z. Feng, M. Vollmann, Phys. Rev. D 95 (2017) 035018, arXiv:1608. 03894 [hep-ph]. [78] M. Duerr, F. Kahlhoefer, K. Schmidt-Hoberg, T. Schwetz, S. Vogl, J. High Energy Phys. 09 (2016) 042, arXiv:1606.07609 [hep-ph]. [79] T. Jacques, A. Katz, E. Morgante, D. Racco, M. Rameez, A. Riotto, J. High Energy Phys. 10 (2016) 071, arXiv:1605.06513 [hep-ph]. [80] N.F. Bell, Y. Cai, R.K. Leane, J. Cosmol. Astropart. Phys. 1701 (2017) 039, arXiv: 1610.03063 [hep-ph]. [81] F. D’Eramo, B.J. Kavanagh, P. Panci, J. High Energy Phys. 08 (2016) 111, arXiv: 1605.04917 [hep-ph]. [82] C. Englert, M. McCullough, M. Spannowsky, Phys. Dark Universe 14 (2016) 48, arXiv:1604.07975 [hep-ph]. [83] A.J. Brennan, M.F. McDonald, J. Gramling, T.D. Jacques, J. High Energy Phys. 05 (2016) 112, arXiv:1603.01366 [hep-ph]. [84] V. De Romeri, E. Fernandez-Martinez, J. Gehrlein, P.A.N. Machado, V. Niro, J. High Energy Phys. 10 (2017) 169, arXiv:1707.08606 [hep-ph]. [85] Y. Cui, F. D’Eramo, Phys. Rev. D 96 (2017) 095006, arXiv:1705.03897 [hep-ph]. [86] S.F. Mantilla, R. Martinez, Phys. Rev. D 96 (2017) 095027, arXiv:1704.04869 [hep-ph]. [87] M. Bauer, S. Diefenbacher, T. Plehn, M. Russell, D.A. Camargo, arXiv:1805. 01904 [hep-ph], 2018. [88] E. Aprile, et al., XENON, arXiv:1805.12562 [astro-ph.CO], 2018. [89] D.S. Akerib, et al., arXiv:1608.07648 [astro-ph.CO], 2016. [90] A. Tan, et al., PandaX-II, Phys. Rev. Lett. 117 (2016) 121303, arXiv:1607.07400 [hep-ex]. [91] X. Cui, et al., PandaX-II, Phys. Rev. Lett. 119 (2017) 181302, arXiv:1708.06917 [astro-ph.CO]. [92] G. Belanger, F. Boudjema, A. Pukhov, A. Semenov, Comput. Phys. Commun. 176 (2007) 367, arXiv:hep-ph/0607059 [hep-ph]. [93] G. Belanger, F. Boudjema, A. Pukhov, A. Semenov, Comput. Phys. Commun. 180 (2009) 747, arXiv:0803.2360 [hep-ph]. [94] G. Belanger, F. Boudjema, A. Pukhov, A. Semenov, Comput. Phys. Commun. 185 (2014) 960, arXiv:1305.0237 [hep-ph]. [95] P.A.R. Ade, et al., Planck, Astron. Astrophys. 594 (2016) A13, arXiv:1502.01589 [astro-ph.CO]. [96] T. Sjostrand, S. Mrenna, P.Z. Skands, Comput. Phys. Commun. 178 (2008) 852, arXiv:0710.3820 [hep-ph]. [97] K. Kong, J.-C. Park, Nucl. Phys. B 888 (2014) 154, arXiv:1404.3741 [hep-ph]. [98] B. Dutta, Y. Gao, T. Ghosh, L.E. Strigari, Phys. Rev. D 92 (2015) 075019, arXiv: 1508.05989 [hep-ph]. [99] M. Duerr, P. Fileviez Perez, J. Smirnov, arXiv:1506.05107 [hep-ph], 2015. [100] S. Profumo, F.S. Queiroz, C.E. Yaguna, Mon. Not. R. Astron. Soc. 461 (2016) 3976, arXiv:1602.08501 [astro-ph.HE]. [101] M. Aaboud, et al., ATLAS, J. High Energy Phys. 10 (2017) 182, arXiv:1707.02424 [hep-ex]. [102] D.A. Camargo, L. Delle Rose, S. Moretti, F.S. Queiroz, arXiv:1805.08231 [hepph], 2018. [103] M. Papucci, K. Sakurai, A. Weiler, L. Zeune, Eur. Phys. J. C 74 (2014) 3163, arXiv:1402.0492 [hep-ph].