The Dark Z' Portal: Direct, Indirect and Collider Searches

We perform a detailed study of the dark Z' portal using a generic parametrization of the Z'-quarks couplings, both for light (8-15)GeV and heavy (130-1000)GeV dark matter scenarios. We present a comprehensive study of the collider phenomenology including jet clustering, hadronization, and detector artifacts, which allows us to derive accurate bounds from the search for new resonances in dijet events and from mono-jet events in the LHC 7TeV, LHC 8TeV, and Tevatron 1.96 TeV data. We also compute the dark matter relic abundance, the relevant scattering cross sections and pair-annihilation spectrum, and compare our results with the current PLANCK, Fermi-LAT and XENON100/LUX bounds. Lastly, we highlight the importance of complementary searches for dark matter, and outline the excluded versus still viable parameter space regions of the dark Z' portal.


Introduction
The nature of the Dark Matter (DM) that permeates our Universe is one of the most mysterious puzzles in current science. A variety of dark matter candidates have been proposed and investigated. One of the most compelling dark matter candidates is the so-called WIMP (Weakly Interacting Massive Particle): WIMPs naturally arise in many independently well-motivated theories for physics beyond the Standard Model, and produce a thermal relic density often in accordance with the observed dark matter density.
Indirect searches for DM pair-annihilation in our Galaxy have also resulted in tentative evidence for light dark matter particles. In particular, the ∼ 1 − 3 GeV gamma-ray excess observed at the Galactic center in the Fermi-LAT data can be plausibly explained either by a ∼ 8 GeV DM particle annihilating into τ + τ − or a ∼ 25 GeV one annihilating mostly into bb [21][22][23][24][25][26]. This excess has been reported by different groups [21][22][23][24][25][26]. Astrophysical uncertainties, such as the gas distribution surrounding the Galactic Center and unaccounted-for astrophysical sources such as unidentified pulsars, however, blur the significance of this signal, which thus warrants additional critical scrutiny. A few analyses which use AMS-02 and Fermi-LAT data also seem to disfavor a light dark matter scenario [27][28][29]. Besides this controversial gamma-ray excess, the positron excess observed in the AMS-02 data [30] and before that in Fermi-LAT and Pamela data, which might also be explained by a ∼ 500 GeV DM, seems to be more plausibly due to a local pulsar population [31]. Furthermore, an exciting 130 GeV gamma-ray line observed in the Fermi-LAT data and reported by different groups [32][33][34], might potentially be associated with the two-photon annihilation of a 130 GeV DM particle, but might also reflect an instrumental effect. Be it as it may, those different indirect detection signals are quite interesting and promising and are expected to be settled in either direction in the near future. In this work we will use the current Fermi-LAT dwarfs bounds [29] to constrain the annihilation cross section of the WIMPs, but will refer to some of the mentioned tentative signals in our choice of dark matter particle masses.
Regarding collider searches, the basic signature of DM production is the presence of missing energy, because WIMP DM particles generically escape the detector. For this reason collider searches for DM production typically involve jets + E T data. Colliders provide important complementary bounds mostly in the light dark matter window where direct detection experiments are threshold limited, [35][36][37][38]. We will see that collider bounds are moderately sensitive to dark sector features such as the mass of the DM particle and couplings. However, the coupling strength of the Z with the quarks is very important in setting the size of the production cross section and the hardness of the jets. In our analysis, we will derive constraints using the LHC 7 and 8 TeV data plus the Tevatron 1.96 TeV Refs. [39][40][41] Our goal in the present study is to outline the viable parameter space of the Z portal by taking into account the complementarity of direct, indirect and collider searches [42][43][44][45][46]. Our work goes beyond and differs from previous studies in several ways, including the following: • We use an effective Lagrangian approach adding together spin dependent and independent couplings.
• We include indirect detection bounds.
• We perform a comprehensive collider study by including jet clustering and hadronization, and by simulating detector effects.
• We outline the viable parameter space in the Z mass Z − DM − DM plane, after plugging into the most relevant colliders (LHC 8,LHC 7, Tevatron 1.96 TeV), direct (LUX,XENON100) and indirect (Fermi Dwarfs) limits.

A Z Portal Dark Matter Model
New heavy neutral gauge bosons, Z , appear in many gauge extensions of the Standard Model (SM). In particular, in certain models these additional gauge bosons are responsible for linking the dark and visible sector [53][54][55][56][57][58][59][60][61] producing a so-called Z "portal". In this work we aim to investigate this setup under three different perspectives: direct, indirect and collider searches. To do so, we use a general Lagrangian that reads, where q i 's are denoting the SM quarks, a and b are constant factors, which are equal to unity if one assume that the Z couples equivalently to the SM Z boson, whereas χ and g χ are the dark matter particle and the DM − DM − Z coupling, respectively. Here g i V and g i A are the vector and axial Z-quarks couplings, which read [65], Q i is the charge of the quark i,t 3L is the weak isospin of quarks, with t 3L = +1/2 (−1/2) for up (down) quarks. From Eq.(2.1) we see that: • The Z is Leptophobic, i.e. it just doesn't couple to leptons by construction. This is the most constrained scenario by collider searches concerning a new neutral gauge boson once the Z can only decay to jets and DM pairs. On the other hand, DMnucleon interactions do not depend upon the details of the leptonic sector.
• The constant factors a and b are equal to unity in the SM. Here, we will explore different particle physics models by varying these constants. We will assume that a = b and investigate two scenarios: The latter scenario corresponds to the case where the Z -quarks couplings are suppressed by 50% in comparison with Z-quarks ones.
• As described in Table 2, the Dark Z portal might give rise to four different operators when χ is a Dirac Fermion [45][46][47][48][49][50]. If χ is a Majorana fermion, the vector current is zero. Therefore, in the latter case we are left with the O 2 and O 4 operators. In this work we will not study these operators individually. We consider them all at once. Therefore we assume that χ is a Dirac fermion with vector and axial coupling to fermions, as described in Eq.(2.1).
• Z gauge bosons are predicted to exist in extended gauge theories, such as U (1) X . In principle, the setup investigated here is supposed to include extra fermions to cancel the anomalies induced by the U (1) X . In our scenario the mass of the extra fermion  Table 1. Effective operators for DM-Nucleon scattering. We have classified the Spin Independent (SI) and Spin Dependent (SD) operators. v is the velocity of DM in the lab frame, M n is the nucleon mass and ∆Σ is defined as N | qq γ µ γ 5 q|N = ∆ΣŪ N γ µ γ 5 U N , with U N as the wave function of the nucleon [45]. The g V and g A couplings are determined according to Eq.(2.1).
should be M f ≤ 64π 2 /(a * g 3 V )M Z [66]. However, this bound may be circumvented if extra fermions are introduced in the model. Those extra fermions might not be related to dark matter observables. Since we are studying the Z portal from a general perspective, we will not include extra fermions in our analysis and we will assume that they have negligible impact on the phenomenology associated with the dark matter particle.

Direct Detection
In direct detection, the relevant observables are the dark matter-nucleon scattering cross section and the particle dark matter mass. In general, in the low-velocity limit, WIMP-Nucleus scattering can be either spin-independent (SI) or spin-dependent (SD), depending on what sort of effective couplings are involved. As shown in Table 2, the O 3 and O 4 operators induce spin-dependent interactions, whereas O 1 and O 2 spin independent ones. Therefore the Z portal is subject to both spin-dependent and spin-independent bounds.
Spin-dependent bounds are generically weaker than Spin Independent, because Spinindependent scattering is proportional to A 2 unless a destructive interference happens, which is not the case here [19,20]. As far as direct detection is concerned only the pure vector and vector-axial operators are relevant, because a mixing between these two described by the O 2 and O 3 operators are velocity suppressed, as shown in Table 2. The scattering cross section for the pure vector and vector-axial coupling are given in Table  2. Since spin-dependent bounds provide weaker bounds on the parameter space of the Z portal sort of models we will focus our analysis on the spin-independent case only [42][43][44][45][46]. Note that the argument above fails in theories where the vector coupling is extremely suppressed compared to the vector-axial coupling. The vector and axial couplings might be different, as it happens in many models, but they would presumably not differ by orders of magnitude [53][54][55][56][57][58][59][60][61].
In any case, hereafter our results regarding collider, indirect detection and direct detection searches are obtained using the Lagrangian given in Eq. (2.1) without assuming any particular operator. It is important to notice that despite being a leptophobic theory, the inclusion of leptonic channels would have no impact on the direct detection bounds derived here. The lepton channels are rather relevant for indirect and collider purposes, though [67]. Throughout this work we use the Micromegas package to compute direct detection observables [68][69][70].

Indirect Detection
Indirect searches for dark matter probe different and complementary dark matter observables to direct detection, namely: the annihilation cross section, the dark matter particle mass, and the dark matter halo model. Different halo models have tend to converge to similar trends at distances far enough from the center of given astrophysical objects such as the Galactic center, but they might vary from cuspy to core types at small radial distances [71][72][73]. It is still open to debate whether or not the halo profile of our own Milky Way is cuspy or cored. The Fermi-LAT collaboration [29], as well as other independent studies [26][27][28], have set stringent limits on the annihilation cross section of dark matter particles under different halo profile and annihilation channels assumptions. Here we will use as reference the current Fermi-LAT ones, which are focused on the bb annihilation channel.
It is important to emphasize that, because we are in the leptophobic regime, the dark matter particle χ will always annihilate into quarks, which all feature similar γ-ray spectra, resulting from the neutral pion two-photon decay. We will sum the annihilation modes into all quarks as if they produced the same gamma-ray spectrum of bb case. This approximation is quite reasonable because in fact all quarks produce basically the same gamma-ray spectrum (the possible exception being the top quark near threshold). Once we have summed up the total annihilation cross section into quarks, we compare the result with the current Fermi-LAT bounds, which was obtained from stacked 4-years of observations of local dwarf spheroidal galaxies, assuming an annihilation 100% into bb. The Fermi-LAT constraints can be thus straightforwardly applied to our leptophobic setup.

Bounds from the search for new resonances in dijet events
Several extensions of the SM predict the existence of new heavy particles that are within the reach of hadron colliders. A new heavy neutral gauge boson as the Z , for example, may lead to a resonance in jj and + − invariant masses, m jj and m , respectively.
The Tevatron and LHC collaborations have been searching for resonances in dijet events with null results until now [39][40][41]. These null results allow to place strong constraints on any new model predicting resonances in dijet events, as in our case, where the dark Z has a large branching fraction into quarks due to leptophobia.
Concerning the way a dijet search might constrain a Z dark matter model, the effect on the mass and couplings of the DM particle is indirect. For a fixed Z mass, the branching Also, decreasing the coupling g χ between the DM and Z increases the branching faction into quarks. We show in Fig. 1 the branching ratio into quarks of a Z of fixed mass as a function of the DM mass M χ and the g χ coupling.
From Fig. 1 we see that a dijet search for a Z is moderately sensitive to the DM model details. For not too large g χ couplings, the branching ratio is never smaller than ∼ 70% independent of the DM mass. On the other hand, the Z mass and couplings to quarks are crucial, not only because of the size of the production cross sections, but also because heavier resonances give rise to harder yields, in this case harder jets which are likely to pass selection cuts.
We now discuss in detail the constraints from Tevatron, LHC 7 and LHC 8 data.

Tevatron 1.96 TeV
To evaluate the impact of the Tevatron search on dijet resonances after 1.13 fb −1 of integrated luminosity [39], we simulated the process plus up to two extra jets using MadGraph5 [74]-FeynRules [75], clustering and hadronizing jets with Pythia [76], and simulating detector effects with PGS4 [77]. Soft and collinear jets from QCD radiation generated by Pythia are consistently merged with the hard radiation calculated from matrix elements in MLM scheme [78] at appropriate matching scales. We adopted the CTEQ6L parton distribution functions computed at µ F = µ R = M Z . All the collider events simulated in this work were obtained using these packages. The signal events were then selected with the same criteria adopted in Ref. [39] which, by the way, are rather inclusive. The only requirement is on the jets rapidities y j of an event, As in the Tevatron analysis of Ref. [39], we also multiplied the cross sections by a K-factor of 1.3 to account for higher order QCD corrections. We generated events for Z of masses from 300 GeV to 1.4 TeV, the mass region for which the Tevatron data have sensibility, and applied the 95% C.L. upper limits, quoted in Ref. [39], on the production cross section times branching ratio for a Z after imposing Eq. (5.3).
The exclusion regions, in the M Z versus g χ , from Tevatron searches, can be seen in Figures 2 and 4. The wavy aspect of the curves just reflects the shape of the experimental exclusion regions [39].

LHC 7 TeV
As in the case of the Tevatron, we generated signal events for the process at the LHC at 7 and 8 TeV center-of-mass energy, plus one and two additional hard jets. The ATLAS collaboration performed an early resonance search in dijet invariant masses using 1 fb −1 of data, placing 95% C.L. upper limits on σ×BR [40]. Background suppression was achieved imposing the following cuts on dijets events which we also applied to our signal events, where y * = (y 1 −y 2 )/2 is the rapidity of the two highest p T jets in their mutual CM system. The constraints on the dark Z , in this case, could not be straightforwardly taken from the upper limits quoted in Ref. [40], once a Gaussian model template, where events are normally distributed in m jj , was used to derive the limits, and the jets invariant mass of the heavy Z resonances are skewed distributions with considerable asymmetries.
In order to obtain the attainable limits from the LHC 7 TeV data we fit the observed dijet invariant mass to the functional form in the 1 to 4 TeV range where x = m jj / √ S and the p i are fit parameters.
Invariant mass distributions for signal events corresponding to Z masses from 900 GeV to 3 TeV were generated and the same Bayesian method used in the experimental study to obtain the new 95% C.L. limits on σ×BR. We have assumed a flat prior probability density for the number of signal events and marginalized over a nuisance parameter to account for systematic uncertainties on the signal acceptance. This nuisance parameter was first tuned to reproduce the limits quoted in Ref. [40] for a Gaussian m jj distribution for each M Z . After tunning the systematic uncertainty, we calculated the likelihood function and the posterior probability density function upon which we obtained the new upper limits.
Compared to Tevatron, the LHC 7 TeV exclusion regions are narrower concerning the Z masses, but it reaches bigger g χ values as can be seen in Figures 2. This is consequence of the much more restrictive selection cuts of LHC 7 TeV analysis necessary to suppress the QCD backgrounds, and larger production cross sections compared to Tevatron. Next we comment the dijet production at the LHC 8 TeV.

LHC 8 TeV
The CMS collaboration performed a search for narrow resonances using dijets with 19.6 fb −1 of data at the LHC 8 TeV [41]. Contrary to the ATLAS search [40], 95% C.L. limits on a Z model were derived in the 1 to 5 TeV range. In this case, all jets in the analysis were requested to have and if one of two highest p T (leading) jets fails to pass these cuts, the event is discarded.
To reduce the sensitivity to gluon radiation, the remaining jets are combined into wide jets [79], J, which are constructed from the leading jets, j 1,2 , in an event by adding the four-vectors of all other jets to the closest leading jet if ∆R ij = (∆R ij ) 2 + (∆φ ij ) 2 < R w , i = 1, 2 and j denotes a non-leading jet. As in the experimental study we set R w = 1.1 to form wide jets in signal events. Now, the dijet system is composed of two wide jets upon which we impose the following additional cuts Compared to the Tevatron and LHC 7 TeV searches, the CMS LHC 8 TeV analysis is able to exclude bigger g χ couplings for a given DM mass in the ∼ 1.2 to 1.7 TeV mass range as can be seen in Fig. 2 for DM masses from 8 to 500 GeV. The LHC 7 TeV has a deeper reach compared to the Tevatron, but in a narrower mass range. This is consequence of a more selective set of kinematic cuts, as we discussed. The LHC 8 TeV, on the other hand, excludes a larger region of the M Z versus g χ plane not only because of the larger production cross sections, but also due to the more efficient selection of signal events and the much larger integrated luminosity as compared to previous studies. Note, however, that there exist small gaps between the M Z regions covered by the three sets of data, mainly between LHC 7 TeV and LHC 8 TeV. All those gaps could be closed using the whole of data accumulated by the Tevatron and 7 TeV run of the LHC. This notwithstanding, there is an almost complete complementarity between the 3 experiments for small g x couplings. It should also be pointed out the complementarity between dijet and monojet searches as we are going to discuss next.

Bounds from DM searches in the monojet channel
The associated production of DM and jets, photons or gauge bosons, has been extensively studied both from the theoretical and experimental sides in the search for DM in colliders.
In models with a dark mediator, as the Z , the process we are interested in is where a jet from QCD radiation is irradiated from the initial state partons alongside Z which, then, decays to a pair of DM particles giving rise to a hard jet and large missing energy signal. As in the case of dijets, we simulated events with one additional jet and matching the hard matrix element contributions from MadGraph5 to soft and collinear jets from Pythia at an appropriate matching scale in the MLM scheme. The renormalization and factorization scales were chosen dynamically in this case as µ R = µ F = p 2 T + 1 2 M 2 Z where p T is the transverse momentum of the hardest jet in the event.
Combining dijet and monojet searches for DM is interesting once they cover complementary regions of the M Z versus g χ plane. As we discussed, dijets have a larger sensitivity to small g χ couplings, as a larger Z χχ coupling increases the branching ratio do DM. On the other hand, monojets have an increased production rate for larger g χ once Z decays to DM in this case. Also, the sensitivity to the DM mass is complementary between the two processes. As the DM gets heavier, the BR(Z → χχ) decreases for a fixed M Z , which enhances the dijet rates, but suppresses the number of monojet events. This effect can be seen in Figure 3 where we show the production rates as a function of the DM mass.
The CMS collaboration has performed a search for new physics in monojet events with 19.5 fb −1 of data at the 8 TeV LHC [80], placing 95% C.L. upper limits on new physics events passing the following selection criteria where p T j 1 is the transverse momentum of the leading jet of the event. A second jet is allowed provided it is not too far away from the leading jet in azimuthal direction, ∆φ(j 1 , j 2 ) < 2.5. Events with more than two jets with p T > 30 GeV and |η| < 4.5 are vetoed.
Seven E miss T regions were used to select signal events: 250, 300, 350, 400, 450, 500, and 550 GeV. Requiring the signal events to satisfy these criteria excludes adjacent and intersecting regions of the M Z versus g χ plane. The resulting 95% C.L. exclusion region is shown in Fig. 2 in the case of an 8, 50, 130 and 500 GeV DM mass.
It should be pointed out how the complementarity between dijet and monojet searches excludes entire regions of the parameters space. Except for small gaps, all Z masses from 300 GeV to 1.7 TeV are excluded for all Z couplings to DM at 95% confidence level in the case of an 8 GeV DM mass, for example, as is shown in the upper left panel of Figure 2. In fact, very similar regions are excluded for DM masses up to 130 GeV, at least, as can be seen in Figure 2. Even for a 500 GeV DM, Z masses fro 1.2 to 1.6 TeV are excluded at 95% C.L. for all g χ as we see in the lower right panel of Fig. 2. Some models predicting a new heavy gauge boson, such as 331 models, and other models with extended gauge groups, might present reduced gauge couplings between quarks and the new gauge bosons. Reducing the Z − q − q couplings, g Z qq , shrink all exclusion regions collectively, and both dijet and monojet data are more easily evaded. We show the effect of reducing g Z qq by 90% and 50% for a 130 GeV DM in Fig. 4.
We now discuss the complementarity among direct, indirect and collider searches of DM and highlight their interplay to constrain a dark Z model such as the one considered in this work.

Complementary Results: Direct, Indirect Detection and Collider
So far we have discussed direct, indirect and collider searches for DM in the particular Z portal scenario under consideration. We here investigate how the parameter space of the Z portal is constrained by these bounds. We will point out how crucial complementarity is when it comes to dark matter searches. Direct detection bounds such as the XENON100/LUX and the LHC 8 TeV jet+E miss T limits, are both more sensitive to larger g χ (DM − DM − Z ) couplings. Indirect detection constraints coming from Fermi-LAT Dwarf Spheroidals data rule out light Z masses and a wide range of g χ couplings. Other collider bounds coming from the LHC 7 TeV and Tevatron 1.96 TeV dijet data give complementary bounds on small g χ couplings for larger Z masses. We select WIMP masses based on the tantalizing signals reported in both direct and indirect searches, although we will not attempt to fit for those signals with the model under investigation.

8 GeV WIMP in Light of DAMA,CDMS-Si,CoGeNT and CRESST
An 8 GeV WIMP is a well motivated dark matter candidate because of recent positive signals coming from direct detection. In particular, the CoGeNT collaboration has claimed to observe a ∼ 2σ modulation consistent with a 8 GeV WIMP scattering off nuclei with a spin independent cross section of 3−4×10 −5 pb [12,13]. Moreover, the DAMA experiment has observed with ∼ 9σ an annual modulation also consistent with a WIMP scattering [11]. Furthermore, the three excess events reported by CDMS-Si could also be plausibly explained by a similar WIMP [15]. For this reason, it is interesting to ask whether a 8 GeV WIMP via the Z portal is a feasible explanation to these signals, bearing in mind the mentioned LUX and XENON100 constraints.
First, we show in Fig.5 the spin independent cross section as a function of the Z mass with g χ varying randomly from zero to one, where the blue (green) points correspond to the underabundant (overabundant) regimes. In the left panel we have assumed the coupling factors a = b = 1. Thus we are in the regime that the Z boson couples equivalently to the SM Z. While in the right panel we have used a 50% suppression in the Z -quarks couplings, i.e, with a = b = 0.5. We thus conclude that in both cases the Z portal with a 8 GeV DM particle is excluded by the XENON and LUX constraints in the sense that the region of the parameter space that sets the right abundance (thin line between green and blue points) is ruled out by orders of magnitude. In other words, Z mediated processes do not offer the necessary "XENONPHOBIA" described in [19,20], and therefore do not provide a viable explanation for these modulation signals, while still being consistent with other searches. This conclusion might be different if the Z portal is not responsible for setting the WIMP abundance.
With this in mind, in Fig.6 we show the importance of the complementary search for bounds on the Z portal. Assuming that the connection between the visible and the dark sector is given exclusively through the Z portal with the coupling factors a = b = 1 we have computed the abundance and the black line in Fig.6 reproduces Ωh 2 = 0.12. The region beneath the black line provides Ωh 2 < 0.12, whereas the one above the line Ωh 2 > 0.12.
In Fig.6 we also exhibit the parameter space M Z × g χ along with the collider, direct and indirect detection bounds. The white gaps are not ruled out by any constraint. Pink and Grey regions are ruled out by the LHC 8 TeV and 7 TeV using the dijet data. Yellow is M Z' (GeV) 5 Fig.6 that with the recent LUX results the light mediator region is excluded only the mass range 1200 < M Z < 1100 with g χ < 0.25 we can evade all constraints. Apart from these, only heavy mediators M Z > 1.7 TeV, are allowed by the current data.
In order to probe different Z portal particle physics models we changed the Z -quarks coupling by 50%, i.e, with the coupling factors determined in Eq.2.1 as a = b = 0.5. Again we have calculated the abundance in this regime and drawn a black line in Fig.7 which reproduces Ωh 2 = 0.12. The region under (above) the black line provides Ωh 2 < 0.12 (Ωh 2 > 0.12). It is important to emphasize that we have assumed that the connection between the visible and the dark sector is given exclusively through the Z portal with the coupling factors a = b = 0.5. Furthermore, in Fig.7 we exhibit the parameter space M Z × g χ along with the collider, direct and indirect detection bounds. The white gaps are not ruled out by any constraint. Pink region is ruled out by the LHC 8 TeV dijet data. Yellow is excluded by Tevatron (1.96 TeV) with the dijet data. The dashed blue region is excluded by LHC at 8 TeV using the jet + E miss T data. The red (blue) dashed region is the XENON100 (LUX 2013) excluded region, whereas the green region is ruled out be Fermi-LAT Dwarfs. We note how significant the impact of suppressing the Z -quarks couplings by 50% is. It opens up a large region of the parameter space consistent with the current bounds. Therefore light and heavy Z bosons with suppressed couplings with quarks are totally a viable annihilation channel as long as they do not set the abundance of the dark matter particle. As we can see in Fig.7, the black line that delimits the right abundance parameter space lies in the very light mediators region only, but such light mediators are excluded by both current XENON100/LUX and Fermi-LAT Dwarfs bounds.

15 GeV WIMP in Light of CDMS-Si Excess Events
A 15 GeV WIMP is also well motivated because of the three recent excess events observed by CDMS-Si. Despite a likelihood analysis favors a 8 GeV WIMP, a 15 GeV DM particle with a SI cross section of 2 × 10 −6 pb is perfectly capable of explaining the signal. Besides, the excess of gamma-ray emission from the Galactic Center which might as well be partially and plausibly explained by a ∼ 15 GeV WIMP that annihilates mostly in bb [21]. For these reasons we will investigate the Z portal for M χ = 15 GeV.
First, we show in Fig.8 the spin independent cross section as a function of the Z mass with g χ varying randomly from zero to one. In the left panel we have assumed the coupling factors a = b = 1, while in the in the right panel we have used a 50% suppression in the Z -quarks couplings, i.e, with a = b = 0.5. We again conclude that in both cases the Z portal with 15 GeV dark matter particle is excluded by the XENON100 (top) and LUX (bottom) bounds. Therefore a 15 GeV DM particle which scatters off nuclei via the Z boson is not a viable explanation for this excess events while obeying other constraints. One would need more exotic gauge coupling structures to acquire the "XENONPHOBIA" necessary to evade the XENON and LUX bounds.
In Fig.9 draw a black line which delimits the right abundance (Ωh 2 = 0.12) parameter space for a 15 GeV DM particle as a function of the DM-DM-Z coupling (g χ ) assuming that the connection between the visible and the dark sector is given through the Z portal only with the coupling factors a = b = 1, i.e, in the regime that the Z boson has identical couplings to the SM Z. The region on top of the black line provides Ωh 2 > 0.12, whereas the beneath the line one set Ωh 2 < 0.11. Hence Z mediated processes overproduce the dark matter particles unless we are in the regime of light mediator M Z < 500 GeV.
In Fig.9 we show as well the parameter space M Z × g χ along with the collider, direct and indirect detection bounds. The white regions are not ruled out by any constraint. Pink and Grey regions are ruled out by the LHC 15 TeV and 7 TeV using the dijet data. Yellow is excluded by Tevatron (1.96 TeV) with the dijet data. The dashed blue region is excluded by LHC at 15 TeV using the jet + E miss T analysis. The red (blue) dashed region is the XENON100 (LUX 2013) excluded region, whereas the green one is excluded by Fermi-LAT Dwarfs data [29]. We thus conclude from Fig.9 that the light mediator window is now closed and the region with M Z ∼ 1100 for g χ < 0.05 is still viable. Apart from this one, only heavy mediators M Z > 1.7 TeV with g χ < 0.4, are allowed by the current data. Notice that the light mediator regime is still totally ruled out XENON/LUX bounds in comparison with the 8 GeV DM case. This can be explained simply by the energy threshold of XENON/LUX. At sufficiently low energies the XENON/LUX efficiency goes down, and therefore for very light dark matter particle M χ < 8 GeV the XENON/LUX bounds gets weaker. On the order hand, as we increase the mass of the DM particle up to ∼ 15 GeV, XENON/LUX constraints kick in. We clearly see this effect by comparing the XENON100/LUX excluded region in the left panels of Fig.6 and Fig.9. Assuming that the connection between the visible and the dark sector is given exclusively through the Z portal with the coupling factors a = b = 0.5 we have computed the abundance and the black line in Fig.10 reproduces Ωh 2 = 0.12. The region beneath the black line provides Ωh 2 < 0.12, whereas the one above the line Ωh 2 > 0.12. Notice that setting a = b = 0.5 means that we are taking the Z − quarks suppressed in 50% in comparison with the SM Z-quarks ones.
In Fig.10 we plot the parameter space M Z × g χ allowed (white) and excluded (shaded and colorful) by the combined collider, direct and indirect detection bounds. Pink region is ruled out by the LHC 8 TeV dijet data. Yellow is excluded by Tevatron (1.96 TeV) with the dijet data. The dashed blue region is excluded by LHC at 8 TeV using the jet + E miss T data. The red (blue) dashed region is the XENON100 (LUX 2013) excluded region, whereas in green we show Fermi-LAT Dwarfs one. Notice that when we suppress the Z -quarks couplings in 50% new regions of the parameters space show up. In particular, 1 TeV mediators with g χ 0.1 are allowed by the data.

130 GeV WIMP in light of the Fermi line
The 130 GeV gamma-ray line observed in the Fermi-LAT data [32,33] has appeared as a smoking gun signature for a 130 GeV DM particle which has a fairly large annihilation cross section into γγ. Different groups have found the same line feature in the Fermi data [34]. The Fermi-LAT Collaboration, though, has claimed no evidence for such a line [81]. The new PASS8 software plus additional accumulated data will hopefully shed some light on this interesting possibility [82]. Hereunder, we will be agnostic as to the particular particle mechanism responsible for the production of the line, and simply consider the particular value of 130 GeV as interesting for the dark matter particle mass.
First, we show in Fig.11 the spin independent cross section as a function of the Z mass with g χ varying randomly from zero to one. In the left panel we have assumed the coupling factors a = b = 1, while in the in the right panel we have used a 50% suppression in the Z -quarks couplings, i.e, with a = b = 0.5. We can clearly see in the left and right panels of Fig.11 the region between the blue and green point, which delimits the right abundance regime, is ruled out or on the verge of being excluded. Therefore a particle physics model which has a Z only mediating the annihilation and scattering channels is rather disfavored by the current data.
In Fig.12 we have drawn a black line that sets Ωh 2 = 0.12 for a 130 GeV DM particle as a function of the DM-DM-Z coupling (g χ ) assuming that the connection between the visible and the dark sector is given through the Z portal only with the coupling factors a = b = 1, i.e, in the regime that the Z boson has identical couplings to the SM Z. The area beneath the black curve provides Ωh 2 < 0.12, whereas the one above the line Ωh 2 > 0.12. Hence Z mediated processes may provide the right abundance as long as M Z 1 TeV.
In Fig.12 we have shown as well the parameter space M Z × g χ along with the collider, analysis. The red (blue) dashed region is the XENON100 (LUX 2013) excluded region, whereas the green one is excluded by Fermi-LAT Dwarfs data [29]. We can conclude from Fig.12 that the whole light mediator window is excluded and only the region with M Z ∼ 1100 for g χ < 0.05 circumvents the bounds. Apart from this one, only heavy mediators M Z > 1.7 TeV, are allowed by the current data for g χ < 0.3. As we increase the mass larger couplings are allowed by data. For larger couplings the Z mass is excluded up to 3 TeV. Now in Fig.13 we show the impact of suppressing in 50% the Z -quarks couplings with M χ = 130 GeV. It is clear from Fig.13 that a small window for M Z ∼ 500 GeV with suppressed couplings opens up and a somewhat large one for 700 GeV < M Z < 1.2 TeV with g χ < 0.1 becomes allowed by the joined data. Besides these only heavy mediator are consistent with the current data.
6.4 Results for M χ = 50, 500, 1000 GeV WIMPs In this final section we consider three additional values for the WIMP mass: 50, 500 and 1000 GeV, and study how our results extrapolate to larger masses. The structure of all plots is the same as before: for each value of the mass we show the spin independent cross section as a function of m Z for a = b = 1 (left) and for a = b = 0.5 (right). We then show the over-and under-abundant thermal relic density regions and the bounds from direct, indirect and collider searches both for suppressed and unsuppressed quark couplings.
While nothing qualitatively new emerges for the M χ = 50 GeV case, in the M χ = 500 and 1000 GeV case we notice the new feature associated with the presence of a resonance at M χ = m Z /2. While portions of that resonance are highly constrained by collider searches and by indirect searches (the resonance enhances the pair-annihilation cross section via an M Z' (GeV) 5  s-channel Z exchange), we find interesting portions of parameter space that are viable and that feature the correct thermal relic density.
The feature we highlighted for the 500 and 1000 GeV mass persists for even larger masses, where the portion of parameter space compatible with the thermal relic density requirement and with collider searches continues to be present. The only constraint is given by gamma-ray searches, which benefit from the large annihilation cross sections associated with the resonant annihilation mode.

Conclusions
In this study we have investigated a particular incarnation of the Z dark portal, where the additional gauge boson does not interact with leptons. We carried out a detailed study of direct and indirect searches for the resulting dark matter candidate, as well as an extensive study of the collider phenomenology using Tevatron and LHC results.
We found a high degree of complementarity at several different levels: between direct and indirect dark matter searches, between dark matter searches and collider studies, and between Tevatron and LHC searches. We focused our detailed analyses on specific values of the dark matter particle mass, motivated by tentative signals in both direct and indirect detection. Inspecting the thermal relic density we found that only for large masses can we obtain regions that are not ruled out by dark matter or collider searches and that possess M Z' (GeV) 5   the correct universal abundance. In those regions, the mass of the dark matter is about half the mass of the Z and resonant annihilation produces large cross sections contributing to suppressing the relic density, with other parameters compatible with direct and collider searches.   Figure 19. Result for M χ = 500 GeV in the M Z × DM-DM-Z coupling plane, with Z -quarks couplings suppressed. The gaps evade all bounds. The pink region is ruled out using the dijet data from the LHC 8 TeV. The dashed blue region is excluded by LHC at 8 TeV using the jet + E miss T analysis. The red dashed region is the XENON100 excluded region. The black line reproduces Ωh 2 = 0.12, whereas the regions below and beneath the lines set Ωh 2 > 0.12 and Ωh 2 < 0.12 respectively.