Slac-Pub-15306 Dichromatic Dark Matter
Both the robust INTEGRAL 511 keV gamma-ray line and the recent tentative hint of the 135 GeV gamma-ray line from Fermi-LAT have similar signal morphologies, and may be produced from the same dark matter annihilation. Motivated by this observation, we construct a dark matter model to explain both signals and to accommodate the two required annihilation cross sections that are different by more than six orders of magnitude. In our model, to generate the low-energy positrons for INTEGRAL, dark matter particles annihilate into a complex scalar that couples to photon via a charge-radius operator. The complex scalar contains an excited state decaying into the ground state plus an off-shell photon to generate a pair of positron and electron. Two charged particles with non-degenerate masses are necessary for generating this charge-radius operator. One charged particle is predicted to be long-lived and have a mass around 3.8 TeV to explain the dark matter thermal relic abundance from its late decay. The other charged particle is predicted to have a mass below 1 TeV given the ratio of the two signal cross sections. The 14 TeV LHC will concretely test the main parameter space of this lighter charged particle.
Although dark matter serves as the dominant component of matter in our universe, its various properties remain unknown. From astrophysical evidence, there is no doubt that dark matter can interact with the Standard Model (SM) particles through gravitational interaction. However, whether there are additional interactions between dark matter and SM particles is still a mystery to us. Among several approaches to search for dark matter particles, measuring the cosmic ray spectrum provides the indirect detection of dark matter. Observing a high-energy gamma-ray line has long been believed to be the “smoking gun” of the dark matter detection [1, 2, 3, 4, 5, 6]. Furthermore, the propagation of energetic photons in our Galaxy is less affected by the interstellar gas or Galactic magnetic field. The gamma-ray line signal can even provide the dark matter geometrical profile in our Galaxy.
The detection of celestial gamma-ray line at keV from the inner galaxy, which is believed to be caused by annihilations, was first reported by  and later confirmed by [8, 9, 10, 11]. The total flux of the 511 keV line has been estimated to be around . About annihilations proceed through the intermediate state of a positronium atom, and 25 of these annihilations with opposite spin of and can produce 511 keV line emission [13, 14]. Although this gamma-ray line has been known for decades, the identification of the positron source remains undetermined. Different astrophysical sources have been suggested during the years, but each of the models faces various challenge to explain the observations consistently. The relatively high ratio of the bulge to disk 511 keV emission towards the inner Galaxy seems against its origin from hypernovae and gamma ray bursts, while the constraints on the production rate of high energy positrons also disfavors millisecond pulsars, as well as proton-proton collisions from e.g. microquasars, low luminosity X-ray binary jets, and the central supermassive black hole. Furthermore, pulsars, magnetars, and Galactic cosmic rays are not favored as major sources to the observed 511 keV from the bulge, and stringent constraints on these origin of the 511 keV line was suggested [15, 16].
Besides these astrophysical suggestions, the possibility that dark matter may create the 511 keV gamma-ray line has been widely discussed, mainly motivated by the rather spheroidal, symmetric, bulge-centered morphology. The lack of higher energy gamma ray requires the injection energy of positrons to be below 3 MeV . This motivates studies for both MeV-scale dark matter models [17, 18, 19, 20, 21, 22] and TeV-scale dark matter models with a MeV mass splitting among different dark matter states [23, 24, 25, 26]. Since TeV-scale dark matter with electroweak interaction strength can naturally gives correct thermal relic abundance, those models are more favored. Interestingly, the morphology of the 511 keV signal profile has a peaked structure around the Galactic center, and the sharpness of the peak prefers to have dark matter annihilation rather than decaying as an explanation . Thus we focus on the heavy dark matter scenario, and try to explain the keV INTEGRAL signal via dark matter annihilation.
One popular dark matter model to explain the INTEGRAL signal is the excited dark matter model with an MeV mass splitting . This class of models suffer from the requirement of a large kinematic energy of dark matter to excite the ground state, hence relying on the Boltzmann tail of dark matter velocity distribution. It is under a debate whether the excited dark matter model can generate enough positrons to explain the large gamma-ray flux for the INTEGRAL data [27, 28]. For the 100 GeV dark matter mass region that we will consider in this paper, the situation is worse, because it requires a higher velocity to obtain enough kinematic energy comparing to a TeV mass dark matter. In our paper, we will address this problem and propose a new scenario of the Down-scattering excited Dark Matter (DeDM) to solve the Boltzmann suppression problem of the vanilla excited dark matter models.
More recently, the hint for another gamma-ray line around 130 GeV from the Galactic center has been suggested by analyzing the public data from Fermi Gamma-ray Space Telescope (Fermi-LAT) [29, 30]. The hint becomes even stronger with the template fitting approach, which takes into account the spatial distribution of the LAT events towards the inner Galaxy along with the spectral information . Fermi-LAT Collaboration has confirmed the hint of the peak at GeV using Pass 7 data. The peak shifts to a higher mass at GeV and the significance becomes weaker using the reprocessed Pass 7 data . Such high energy gamma-ray line emission has been considered as a clean signature from dark matter annihilations. Many dark matter models have been constructed to explain the 130(135) GeV gamma-ray line feature (see  and references therein).
The morphology of the INTEGRAL 511 keV and Fermi-LAT 130(135) GeV line shares similar structures: (1) the signal events concentrate at the center of the Galaxy with non-disk like distributions; (2) after smoothing Fermi-LAT signal with INTEGRAL’s angular resolution, they have comparable full widths at half maximum (FWHM) in both the longitudinal and latitude directions. This motivates us to explain both signals using same dark matter particle in our universe. The fittings for both signals prefer annihilation rather than decaying [12, 34]. Having worked out the required annihilation cross sections, we find that the INTEGRAL 511 keV cross section is six to seven orders of magnitude larger than that of the Fermi-LAT 130 GeV line. This large hierarchy of cross sections sets a challenge when constructing a detailed model. However, the order of magnitude is comparable with an electromagnetic loop factor of if the INTEGRAL and Fermi-LAT signals are coming from tree-level and loop-level processes, respectively. This serves as a clue for our model building.
Our paper is organized as following. In 2, we emphasize the similarities of morphologies for both signals and work out the required cross sections. In 3, we propose our model, Down-scattering excited Dark Matter model, to explain both signals. In section 3.1, we first provide a general operator analysis to illustrate the essence of our model and calculate the scales of cutoffs of the effective operators. Then we build up a concrete UV-completion for the operator analysis in 3.2. In 3.3, we discuss the thermal history of our model. One charged particle needs to be long-lived in our UV model, so that we have a semi-natural model to explain the final dark matter relic abundance.
2 Experimental Data
In this section, we discuss the INTEGRAL and Fermi-LAT oberservations in more detail. Since photon is not much affected during its propagation in the Galaxy, photon coming from dark matter annihilation can be used to determine the dark matter distribution in our Galaxy. However, there are subtleties on how to map the INTEGRAL 511 keV gamma-ray line signal profile to the dark matter distribution profile. This is because the low-energy positrons that are generated from dark matter particles can propagate through the interstellar medium and annihilate with electrons to photons away from the production site and bias the inferred dark matter distribution from the 511 keV line morphology. In this paper, we assume that the positron propagation is negligible comparing to the spatial resolution of the INTEGRAL, thus the dark matter profile can be estimated by measuring 511 keV emission morphology. On the other hand, the Fermi-LAT 130(135) GeV photons could directly be generated from dark matter particles, and its morphology can therefore tell us the dark matter profile.
We first compare the morphologies of the INTEGRAL keV and Fermi-LAT 130(135) GeV lines. After smoothing the Fermi-LAT 130(135) GeV line using the angular resolution of INTEGRAL, we find the spatial distributions are comparable to each other. Furthermore, assuming that both signals are generated by dark matter annihilation, we estimate the annihilation cross sections for the two processes. They will serve as inputs for our model building in the rest of the paper.
2.1 Experimental Data
Thanks to a coded mask located m above the detector plane and a specific dithering observational strategy, the spectrometer (SPI) onboard the INTEGRAL observatory can image the sky with a spatial resolution of 2.6 (FWHM). Based on observations recorded from February , 2003 to January, 2 2009, the study in  has obtained the morphology of the 511 keV line towards the inner Galaxy. In Figure 1, we compare the intensity of the 511 keV gamma-ray line as a function of Galactic longitude and latitude with the 130 GeV line profile obtained by fitting 3.7 years Fermi-LAT observations . Especially, the dark green line shows the 130 GeV line profile further smoothed by SPI 2.6 FWHM beam.
Interestingly, both longitudinal and latitude distributions of INTEGRAL are comparable to those of Fermi-LAT after smoothing. Furthermore, both distributions show the tendency to be off-center in the negative longitudinal direction[31, 35]. These similarities motivate the attempts to build one dark matter model to explain both these two signals.
2.2 Dark matter annihilation cross sections for INTEGRAL and Fermi-LAT
As discussed in the previous section, both the 511 keV line and the Fermi-LAT 130(135) GeV line could potentially be explained by dark matter annihilation. In this section, we estimate the required annihilation cross sections for both experiments.
The gamma-ray line intensity in a given direction provided by dark matter annihilation is the line-of-sight integral of the squared dark matter number density along that given direction
with the -factor defined as:
where and are longitude and latitude, and the integral of is along the line of sight. Here, kpc is the distance from the Sun to the galactic center; is the dark matter halo profile; GeV cm is the often-used dark matter density in the Solar system ; the relation between and is ; () is the number of photons (positrons) generated from each dark matter annihilation hard process; is the dark matter mass; and are the annihilation cross sections. We define for self-conjugated dark matter, e.g. a real scalar or a Majorana fermion, and for a complex scalar or a Dirac fermion. is the number of monochromatic photons that the final states could convert to. For Fermi-LAT, , since we assume that only monochromatic photons are produced in the hard process. For INTEGRAL, observations suggest that about of positrons annihilate through positronium formation . Only 1/4 of annihilation takes place in the parapositronium state, which produces two 511 keV photons. So, we have for INTEGRAL.
We consider both the Einasto and the Navarro-Frenk-White (NFW) dark matter profile
with kpc and for Einasto  and for NFW . Using the fitted fluxes for the INTEGRAL signal (the dark matter+disk ones) in Ref. , we obtain the annihilation cross sections as 111Here, we use different parameters for dark matter profiles compared to the ones in Ref. . We simply rescale their signal flux by the ratio of functions, which could bring an uncertainty of .
For the Fermi-LAT 130(135) GeV gamma line, we use the fitted fluxes from Ref.  for both profiles to calculate the cross sections,
To quantify the ratio of the required cross sections for two experimental results, we define . Taking , we have the experimentally measured ratios as
where clearly show a large hierarchy for the two required cross sections. We want to also stress that the astrophysical uncertainties are fairly large and a global fit by combining the INTEGRAL and Fermi-LAT might bring the uncertainty down. The cross section ratio between these two expertiments is . This will be the input for model building in latter sections. Interestingly, this ratio is comparable to the square of the electromagnetic loop factor , which implies these two experimental results may be related by a loop with two electromagnetic vertices. It serves as a clue for model building.
3 Down-scattering excited Dark Matter
There are several interesting features required to construct dark matter models if both signatures are to be explained by the same dark matter particle with a mass at the 100 GeV scale.
The required cross section for the INTEGRAL data is amazingly large. For a simplest estimation on the annihilation rate, we get for around 100 GeV. This estimation is three orders of magnitude smaller than the required one. Additional mechanisms are therefore required to increase the annihilation rate. There are several ways to achieve this and we pay special attention on the resonance enhancement [42, 43, 44, 45, 46, 47].
To explain the INTEGRAL data, primary positron injections from dark matter are required. Since we don’t see any excess for other cosmic rays, the underlying dark matter model should be arranged to treat positron/electron differently from other particles. In principle, this can be achieved either from kinematic constraints or symmetry reasons.
The ratio of the two cross sections is or. The dark matter model should also provide a natural explanation for this hierarchy of two cross sections.
The model should provide correct amount of dark matter relic abundance to be consistent with observation.
In the following, we will provide a particle physics model to incorporate all above four ingredients. Specifically, we will use a resonance particle in the -channel to increase the dark matter annihilation cross section required for 511 keV gamma-ray line. The kinematic constraints from a small mass splitting will be used in this paper to select positron/electron as the signals from dark matter annihilation. Instead of introducing a light mediator, e.g. dark photon, for the dark matter sector connecting to the positron/electron, we use photon as a more natural mediator to achieve this goal. Noticing that a neutral scalar field cannot decay into another neutral scalar field plus one on-shell photon, which is the reason why , a neutral scalar coupling to photon with the charge-radius operator can naturally induce without generating a photon signal in the meanwhile. The mass difference of the two scalar fields is chosen to be small such that the kinematic energy of is small enough to be consistent with observation. To explain the ratio of the cross sections, we will have the cross section for INTEGRAL to be controlled by coefficients of renormalizable operators, while loop-generated higher-dimensional operators for Fermi-LAT. We will first perform an operator analysis and then provide a UV-complete model.
3.1 Operator analysis
We introduce one Dirac fermion and one complex scalar field in the dark matter sector. Both and are stable particles and coexist in our current Universe. In our study, we will assume that the dark matter component occupies the majority of the dark matter energy, but we will come back to discuss the relative relic abundances of them later. The interactions of the dark matter sector to the SM particles are described by the following set of effective operators
where we implicitly assume that the higher-dimensional operators can be generated at one-loop level. The annihilation of ’s is through exchanging the real scalar in the -channel. For the INTEGRAL data, a small mass scale at around 1 MeV is required to generate positrons almost at rest. In our model, we introduce this small mass scale as the mass splitting of and from such that and MeV. Noticing that the parameter explicitly breaks the global , so the smallness of is technically natural. Expanding the last operator in terms of and , we have
Using the equation of motion, one can rewrite the above operator as . This indicates that cannot decay to a mass-on-shell photon. For , we have the leading decay channel of as
Photon, naturally, behaves as a mediator for the dark matter sector to generate positrons.
The processes to generate positrons for INTEGRAL and photons for Fermi-LAT are shown in Fig. 2, where the solid thick points indicate higher-dimensional operators for those vertices. Although it looks like that the relative cross sections for those two processes are unrelated to each other, we will show in a concrete renormalizable model that the overall cross sections could have a relation in §3.2. In order to generate slowly moving positron from dark matter annihilation, as preferred from the INTEGRAL data, there are two conditions required: (1) the mass splitting should be close to ; (2) cannot have a large boost. The first condition can be satisfied by choosing . The second condition can be arranged by choosing .
We first calculate the annihilation cross section for INTEGRAL. Using the interactions of in Eq. (7), one gets the annihilation cross section of at leading order in as
We are interested in the parameter space with . The decay width of is calculated to be with
Here we treat the decay width of to be approximately the same as .
For INTEGRAL, we need to calculate the velocity-averaged annihilation rate, which is given by
where with determining the variance of the Gaussian dark matter velocity distribution. In our numerical calculation, we neglect the upper limit of the integration, which is controlled by the escaping velocity of dark matter in the galaxy and has only a small effect on our final results.
In Fig. 3, we show the contours of the annihilation rates of in terms of and in the left panel, also and in the right panel. To obtain a large annihilation rate around pb to explain the INTEGRAL data, the resonance mass has to be very close to twice of the dark matter mass. The mass splitting should be a few GeV for the parameter GeV. In Fig. 3, we only presented the results for . The case with has an additional contribution to the total width of from , and has similar results. From the right panel of Fig. 3, we can see that the annihilation rate is insensitive to except for the region with . One might think that can be as light as possible. However, a light generated from dark matter annihilation can have a large Lorentz boost. As a consequence, from decays is also boosted and too energetic to explain the INTEGRAL data . Therefore, we restrict the parameter space in our later study to have at least above 50 GeV.
For the Fermi-LAT signal, instead of obtaining the absolute annihilation rate, we calculate the ratio of the required signal strengths for Fermi-LAT and INTEGRAL. The ratio is equal to the branching ratio of the two decay channels of in Eq. (12), assuming that the additional contribution from the process is small. By taking the ratio, the dependence on the resonance propagator is cancelled and we have
From the effective operator analysis in this section, we have seen that it is possible to explain the required annihilation rates for both INTEGRAL and Fermi-LAT. Our model is economical in a sense that only a few operators and a small number of degrees of freedom are required to explain the data. On the other hand, we should also admit that the resonance requirement of is a tuning point of the parameter space of the current model. Additional ingredients are therefore required to explain this delicate mass relation. We leave this direction of exploration to a future study. Here we emphasize that the ratio of INTEGRAL and Fermi-LAT signals are independent on the way we enhance the annihilation cross sections. Thus one can attach the rest of the model to any other ways of enhancement, e.g. a light mediator in the -channel plus the Sommerfeld enhancement . In the next section, we construct a renormalizable model to UV complete the Lagrangian in Eq. (7) and explain the common origin of the last two operators in Eq. (7).
3.2 Renormalizable model
One way to UV complete the effective Lagrangian in the previous section is to introduce electromagnetic charged states to connect the dark matter sector to photon. In order to have the state stable, at least two charged particles are required to preserve the discrete symmetry associated with . As one example, we introduce two charged complex scalar fields, and . One could also study fermionic charged states in the same procedure. Under or after electroweak symmetry breaking, and have charge one. The global symmetries that we introduce contain a symmetry responsible for the stability of the dark matter particles and a protecting the mass degeneracy of and . We show the field content and symmetries in Table 1.
Based on the symmetries in Table 1, we have the following subset of operators allowed by the symmetries,
Here we only list the operators which are relevant to the processes we concern in this paper. Especially the operator is neglected, which is assumed to have a small coefficient. The last operator generically introduces lepton flavor violation processes, so the couplings ( is the flavor index) should be small.
We first note that the vertices could generate the charge radius operator for the field as shown in the last operator in Eq. (7). After a calculation of the triangle diagram with and propagating in the loop, we get the Feynman rule of the following operator
where and are momenta of and with opposite directions towards the vertex. In the limit , we can match to the coefficient of the effective operator in Eq. (7) as
We notice that the above formula vanishes when . This can be understood by the enhanced discrete symmetry, , in the Lagrangian when and have degenerate masses.222Operator does not preserve this symmetry, but this operator could have a very small coefficient and is irrelevant to this calculation. The charge-radius operator violates this discrete symmetry, thus cannot be generated when . Another more intuitive explanation is to think as a composite particle of and . If the mass of is much heavier than , one can treat as a particle rotating around and have a nonzero charge radius. However, for the mass degenerate case, and should be treated with equal foot and rotate around the center with the same radius. As a result, for each orbit the net charge is zero and the charge radius is zero.
Similarly, we can integrate out and to generate the effective operator coupling to two photons. To match the coefficient in Eq. (7), we have
In the limit , we have
Using the values of in Fig. 4, we anticipate at least the charged particle to have a mass below 1 TeV. This charged particle can decay into one lepton plus one neutrino, for example via the higher dimensional operator .
3.3 Dark Matter Relic Abundance
In our DeDM model, we have two stable particles in our spectrum: and . In our previous analysis, we have assumed that the majority of dark matter in our universe is mainly composed of . To justify our assumption, it is important to study the thermal history of and . In this section, we demonstrate that our setup contains enough ingredients to induce a right relic abundance for , thus it could be the dominant part of the dark matter in our current Universe.
where GeV is the electroweak vacuum expectation value. The function is the width of a Higgs boson in the SM with a mass at . For , GeV and GeV, we have pb and . Thus the relic abundance of can be naturally small.
To satisfy the dark matter relic abundance, a non-trivial thermal history of is needed. This is because a large annihilation cross section in Eq. (10) is needed to explain the INTEGRAL data. The thermal relic abundance of is very small compared to the required dark matter energy density. Noticing that the last operator in Eq. (15) can introduce the decay channel, , the late decay of thermally abundant particles can generate enough , and therefore explain why could be the majority of dark matter.
We first calculate the thermal relic abundance of the charged particle before it decays into and a positron/electron. There are two classes of annihilation channels for . The first class has a photon or boson exchanging in the -channel with final states as a pair of the SM fermions, gauge bosons, and . The second class includes the -channel diagrams, interfering with seagull diagrams. Assuming that the mass is far above the SM particle masses and neglecting the SM particle masses, we have
Here, is the electric charge of the SM fermion; is the axi-vector (vector) couplings of the to the SM fermion up to the electric coupling ; for quarks and 1 for leptons. To derive the above formulas, we have only included the leading terms in for each equation. If the charged particle could occupy the total energy of dark matter in our universe, , the required mass is calculated to be GeV. To derive this mass, we have found that the -wave suppressed annihilation cross section or the terms at is subdominant compared to the total annihilation cross section.
For a lifetime of not too long in the cosmological time scale, we should anticipate that has already decayed into its daughter particle and the final dark matter in our current universe is composed of . On the other hand, the lifetime of can not be too short. Otherwise, the produced particles from decays in the early universe can easily annihilate away and do not provide enough dark matter energy density. To calculate the thermal history of the field, one needs to solve for the following coupled Boltzmann equations between and
Here, in the radiation dominated era, , , , and , where is the number of degrees of relativistic freedom and and . It is convenient to rescale the number density by the entropy and to define the quantity with . The coupled equations become
where and . The final relic abundance is given by , where is the critical density corresponding to a flat universe and with being the entropy today.
At the temperature region with , the decaying terms in Eq. (27) and (28) are not important. The number densities of and reach their separate freeze-out values. Since the cross section of is much larger than , the freeze-out number density for is much below the one of . At a later time, only the last terms in Eq. (27) and (28) become important. One can easily show that the quantity is a conserved number. As a result, the final number density of should just match to the number density of at . So, approximately we have the relic abundance of as
So, the charged particle is predicted to be 3.8 TeV and the other charged particle should be below around 1 TeV to explain the ratio of INTEGRAL and Fermi-LAT cross sections in Eq. (18).
We solve the coupled equations in Eq. (28) numerically and show both comoving number densities of and in Fig. 5. We find that if the mass of is 3.8 TeV, it generates the relic abundance for which satisfies the total dark matter energy density, . In the blue solid and the red solid lines, for pb 333The annihilation cross section of is not necessarily related to its annihilation cross section at the current time. This is because its main production here is from the heavy particle decay, and it has a relativistic velocity and hence a smaller cross section. and s we show the evolutions of the and comoving number densities as a function of temperature. As can be seen from Fig. 5 and at , both and have reached ordinary relic abundances according to their respective annihilation cross sections. At , starts to decay and its number density drops rapidly. Meanwhile, the stable particle number density increases and reaches a plateau at around . The final number density of the field is found to be independent on the lifetime , as long as the decay happens late enough so that the annihilation of is not important any more. The actual time for to reach its eventual number density is proportional to .
To satisfy the dark matter relic abundance, the lifetime of the charged particle can be . For such a late decayed particle, we need to worry about its modification on the Big Bang nucleosynthesis (BBN) history. Since the main decaying product of is into leptons plus the stable , the BBN constrains are fairly weak. From Ref. , the ratio constrains the lifetime of to be for if it would have not decayed. As pointed in Ref. [53, 54], the long-lived charged particle, with a lifetime , can form a bound state with nuclei and enhance the Li production. The parameter space in our model can indeed satisfy the BBN constraints.
4 Discussion and conclusions
The charged particle in our model behaves as a heavy stable charged particle (HSCP) at colliders. The current searches from CMS at TeV and 5.0 fb have set a lower limit on the mass of to be 223 GeV at 95% C.L. . For the predicted mass of around 3.8 TeV, the existing studies have shown that the 14 TeV LHC with 100 fb can reach the HSCP up to a mass around 1 TeV. So, unlikely the stable charged particle can be discovered at the 14 TeV LHC. However, for the other charged particle its mass should be below 1 TeV and could be a long-lived particle or decay into SM particles, for instance . The parameter space of the particle will be well explored at the LHC 14 TeV running.
One feature of our model is directly using photon as a mediator to link the dark matter sector to positron/electron. Unfortunately, other than searching for the charged particles responsible for the charge radius operator, in the near future there is no additional observable dark matter direct or indirect signatures for the field, which has interactions with SM particles suppressed by the TeV scale cutoff of the effective operators. The minor component of dark matter, , may have detectable effects. However, that highly relies on the parameters one chooses, thus we do not pursue that in detail here. Another ingredient that we utilize is the -channel resonance particle to increase the annihilation cross section. We want to stress that this option is not a unique one and is introduced just for convenience. One can also introduce a light mediator in the -channel plus the Sommerfeld enhancement to achieve the same goal .
In summary, we have constructed a realistic model to have the same dark matter particle responsible for both the INTEGRAL 511 keV and Fermi-LAT 135 GeV lines. Through an -channel resonance, the dark matter particles annihilate into a complex scalar, which couples to photon via a charge-radius operator. For a few MeV mass splitting between the real and imaginary parts of the complex scalar, two pairs of electron and positron are the main visible particles from dark matter annihilation. We have worked out the parameter space and have found that both the large cross section required for INTEGRAL and the small cross section for Fermi-LAT can be simultaneously accommodated in our model. The thermal relic abundance of dark matter is achieved by the late decay of a charged particle, which also generates the charge-radius operator. The other charged particle responsible for the charge-radius operator is predicted to have a mass below 1 TeV. The 14 TeV LHC will concretely test the scenario presented in this paper.
We would like to thank James Cline, Tim Cohen, Douglas Finkbeiner, JoAnne Hewett, Dan Hooper, Jessie Shelton, Tracy Slatyer, Aaron Vincent and Jay Wacker for useful discussions and comments. YB is supported by start-up funds from the University of Wisconsin, Madison. YB thanks SLAC for their warm hospitality. SLAC is operated by Stanford University for the US Department of Energy under contract DE-AC02-76SF00515. We also thank the Aspen Center for Physics, under NSF Grant No. 1066293, where part of this work was completed. Support for the work of M.S. was provided by NASA through Einstein Postdoctoral Fellowship grant number PF2-130102 awarded by the Chandra X-ray Center, which is operated by the Smithsonian Astrophysical Observatory for NASA under contract NAS8-03060.
-  L. Bergstrom and P. Ullio, Full one loop calculation of neutralino annihilation into two photons, Nucl.Phys. B504 (1997) 27–44, [hep-ph/9706232].
-  Z. Bern, P. Gondolo, and M. Perelstein, Neutralino annihilation into two photons, Phys.Lett. B411 (1997) 86–96, [hep-ph/9706538].
-  L. Bergstrom, P. Ullio, and J. H. Buckley, Observability of gamma-rays from dark matter neutralino annihilations in the Milky Way halo, Astropart.Phys. 9 (1998) 137–162, [astro-ph/9712318].
-  P. Ullio and L. Bergstrom, Neutralino annihilation into a photon and a Z boson, Phys.Rev. D57 (1998) 1962–1971, [hep-ph/9707333].
-  M. Perelstein and A. Spray, Indirect Detection of Little Higgs Dark Matter, Phys.Rev. D75 (2007) 083519, [hep-ph/0610357].
-  G. Bertone, C. Jackson, G. Shaughnessy, T. M. Tait, and A. Vallinotto, Gamma Ray Lines from a Universal Extra Dimension, JCAP 1203 (2012) 020, [arXiv:1009.5107].
-  W. N. Johnson, III, F. R. Harnden, Jr., and R. C. Haymes, The Spectrum of Low-Energy Gamma Radiation from the Galactic-Center Region., Astrophys. J. Lett. 172 (Feb., 1972) L1.
-  W. N. Johnson, III and R. C. Haymes, Detection of a Gamma-Ray Spectral Line from the Galactic-Center Region, Astrophys. J. 184 (Aug., 1973) 103–126.
-  R. C. Haymes, G. D. Walraven, C. A. Meegan, R. D. Hall, F. T. Djuth, and D. H. Shelton, Detection of nuclear gamma rays from the galactic center region, Astrophys. J. 201 (Nov., 1975) 593–602.
-  M. Leventhal, C. J. MacCallum, and P. D. Stang, Detection of 511 keV positron annihilation radiation from the galactic center direction, Astrophys. J. Lett. 225 (Oct., 1978) L11–L14.
-  R. W. Bussard, R. Ramaty, and R. J. Drachman, The annihilation of galactic positrons, Astrophys. J. 228 (Mar., 1979) 928–934.
-  Y. Ascasibar, P. Jean, C. Boehm, and J. Knoedlseder, Constraints on dark matter and the shape of the Milky Way dark halo from the 511-keV line, Mon.Not.Roy.Astron.Soc. 368 (2006) 1695–1705, [astro-ph/0507142].
-  E. Churazov, R. Sunyaev, S. Sazonov, M. Revnivtsev, and D. Varshalovich, Positron annihilation spectrum from the Galactic Center region observed by SPI/INTEGRAL, Mon.Not.Roy.Astron.Soc. 357 (2005) 1377–1386, [astro-ph/0411351].
-  G. Weidenspointner, C. Shrader, J. Knoedlseder, P. Jean, V. Lonjou, et. al., The sky distribution of positronium annihilation continuum emission measured with spi/integral, Astron.Astrophys. (2006) [astro-ph/0601673].
-  J. F. Beacom and H. Yuksel, Stringent constraint on galactic positron production, Phys.Rev.Lett. 97 (2006) 071102, [astro-ph/0512411].
-  N. Prantzos, C. Boehm, A. M. Bykov, R. Diehl, K. Ferriere, N. Guessoum, P. Jean and J. Knoedlseder et al., arXiv:1009.4620 [astro-ph.HE].
-  C. Picciotto and M. Pospelov, Unstable relics as a source of galactic positrons, Phys.Lett. B605 (2005) 15–25, [hep-ph/0402178].
-  D. Hooper and L.-T. Wang, Possible evidence for axino dark matter in the galactic bulge, Phys.Rev. D70 (2004) 063506, [hep-ph/0402220].
-  C. Boehm, D. Hooper, J. Silk, M. Casse, and J. Paul, MeV dark matter: Has it been detected?, Phys.Rev.Lett. 92 (2004) 101301, [astro-ph/0309686].
-  M. Pospelov, A. Ritz, and M. Voloshin, Secluded WIMP dark matter, Physics Letters B 662 (Apr., 2008) 53–61, [arXiv:0711.4866].
-  D. Hooper and K. M. Zurek, Natural supersymmetric model with MeV dark matter, Phys. Rev. D 77 (Apr., 2008) 087302, [arXiv:0801.3686].
-  J.-H. Huh, J. E. Kim, J.-C. Park, and S. C. Park, Galactic 511 keV line from MeV millicharged dark matter, Phys. Rev. D 77 (June, 2008) 123503, [arXiv:0711.3528].
-  M. Pospelov and A. Ritz, The galactic 511 keV line from electroweak scale WIMPs, Phys.Lett. B651 (2007) 208–215, [hep-ph/0703128].
-  D. P. Finkbeiner and N. Weiner, Exciting Dark Matter and the INTEGRAL/SPI 511 keV signal, Phys.Rev. D76 (2007) 083519, [astro-ph/0702587].
-  N. Arkani-Hamed, D. P. Finkbeiner, T. R. Slatyer, and N. Weiner, A theory of dark matter, Phys. Rev. D 79 (Jan., 2009) 015014, [arXiv:0810.0713].
-  F. Chen, J. M. Cline, and A. R. Frey, New twist on excited dark matter: Implications for INTEGRAL, PAMELA/ATIC/PPB-BETS, DAMA, Phys. Rev. D 79 (Mar., 2009) 063530, [arXiv:0901.4327].
-  F. Chen, J. M. Cline, A. Fradette, A. R. Frey, and C. Rabideau, Exciting dark matter in the galactic center, Phys.Rev. D81 (2010) 043523, [arXiv:0911.2222].
-  R. Morris and N. Weiner, Low Energy INTEGRAL Positrons from eXciting Dark Matter, arXiv:1109.3747.
-  C. Weniger, A Tentative Gamma-Ray Line from Dark Matter Annihilation at the Fermi Large Area Telescope, arXiv:1204.2797.
-  T. Bringmann, X. Huang, A. Ibarra, S. Vogl, and C. Weniger, Fermi LAT Search for Internal Bremsstrahlung Signatures from Dark Matter Annihilation, arXiv:1203.1312.
-  M. Su and D. P. Finkbeiner, Strong Evidence for Gamma-ray Line Emission from the Inner Galaxy, arXiv:1206.1616.
-  F.-L. Collaboration, Search for GammSearch Gamma-ray Spectral Lines in the Milky Way Diffuse with the Fermi Large Area Telescope, .
-  T. Bringmann and C. Weniger, Gamma Ray Signals from Dark Matter: Concepts, Status and Prospects, arXiv:1208.5481.
-  W. Buchmuller and M. Garny, Decaying vs Annihilating Dark Matter in Light of a Tentative Gamma-Ray Line, JCAP 1208 (2012) 035, [arXiv:1206.7056].
-  L. Bouchet, J.-P. Roques, and E. Jourdain, On the morphology of the electron-positron annihilation emission as seen by SPI/INTEGRAL, Astrophys.J. 720 (2010) 1772–1780, [arXiv:1007.4753].
-  M. Su and D. P. Finkbeiner, Double Gamma-ray Lines from Unassociated Fermi-LAT Sources, ArXiv e-prints (July, 2012) [arXiv:1207.7060].
-  G. Jungman, M. Kamionkowski, and K. Griest, Supersymmetric dark matter, Phys.Rept. 267 (1996) 195–373, [hep-ph/9506380].
-  P. Jean, J. Knodlseder, W. Gillard, N. Guessoum, K. Ferriere, et. al., Spectral analysis of the galactic e+ e- annihilation emission, Astron.Astrophys. 445 (2006) 579–589, [astro-ph/0509298].
-  J. F. Navarro, E. Hayashi, C. Power, A. Jenkins, C. S. Frenk, et. al., The Inner structure of Lambda-CDM halos 3: Universality and asymptotic slopes, Mon.Not.Roy.Astron.Soc. 349 (2004) 1039, [astro-ph/0311231].
-  J. F. Navarro, C. S. Frenk, and S. D. White, A Universal density profile from hierarchical clustering, Astrophys.J. 490 (1997) 493–508, [astro-ph/9611107].
-  A. C. Vincent, P. Martin, and J. M. Cline, Interacting dark matter contribution to the Galactic 511 keV gamma ray emission: constraining the morphology with INTEGRAL/SPI observations, JCAP 1204 (2012) 022, [arXiv:1201.0997].
-  M. Ibe, H. Murayama, and T. Yanagida, Breit-Wigner Enhancement of Dark Matter Annihilation, Phys.Rev. D79 (2009) 095009, [arXiv:0812.0072].
-  H. M. Lee, M. Park, and W.-I. Park, Fermi Gamma Ray Line at 130 GeV from Axion-Mediated Dark Matter, ArXiv e-prints (May, 2012) [arXiv:1205.4675].
-  M. R. Buckley and D. Hooper, Implications of a 130 GeV Gamma-Ray Line for Dark Matter, ArXiv e-prints (May, 2012) [arXiv:1205.6811].
-  H. M. Lee, M. Park, and W.-I. Park, Axion-mediated dark matter and Higgs diphoton signal, arXiv:1209.1955.
-  Y. Bai and J. Shelton, Gamma Lines without a Continuum: Thermal Models for the Fermi-LAT 130 GeV Gamma Line, arXiv:1208.4100.
-  G. Chalons, M. J. Dolan, and C. McCabe, Neutralino dark matter and the Fermi gamma-ray lines, arXiv:1211.5154.
-  N. Arkani-Hamed, D. P. Finkbeiner, T. R. Slatyer, and N. Weiner, A Theory of Dark Matter, Phys.Rev. D79 (2009) 015014, [arXiv:0810.0713].
-  R. E. Shrock and M. Suzuki, INVISIBLE DECAYS OF HIGGS BOSONS, Phys.Lett. B110 (1982) 250.
-  V. Barger, P. Langacker, M. McCaskey, M. J. Ramsey-Musolf, and G. Shaughnessy, LHC Phenomenology of an Extended Standard Model with a Real Scalar Singlet, Phys.Rev. D77 (2008) 035005, [arXiv:0706.4311].
-  C. Burgess, M. Pospelov, and T. ter Veldhuis, The Minimal model of nonbaryonic dark matter: A Singlet scalar, Nucl.Phys. B619 (2001) 709–728, [hep-ph/0011335].
-  K. Jedamzik, Big bang nucleosynthesis constraints on hadronically and electromagnetically decaying relic neutral particles, Phys.Rev. D74 (2006) 103509, [hep-ph/0604251].
-  M. Pospelov, Particle physics catalysis of thermal Big Bang Nucleosynthesis, Phys.Rev.Lett. 98 (2007) 231301, [hep-ph/0605215].
-  K. Kohri and F. Takayama, Big bang nucleosynthesis with long lived charged massive particles, Phys.Rev. D76 (2007) 063507, [hep-ph/0605243].
-  CMS Collaboration, S. Chatrchyan et. al., Search for heavy long-lived charged particles in pp collisions at sqrt(s)=7 TeV, Phys.Lett. B713 (2012) 408–433, [arXiv:1205.0272].