We present a scenario where dark matter is in the form of dark atoms that can accommodate the experimentally observed excess of positrons in PAMELA and AMS-02 while being compatible with the constraints imposed on the gamma-ray ux from Fermi/LAT. This scenario assumes that the dominant component of dark matter is in the form of a bound state between a helium nucleus and a particle and a small component is in the form of a WIMP-like dark atom compatible with direct searches in underground detectors. One of the constituents of this WIMP-like state is a metastable particle with a mass of 1 TeV or slightly below that by decaying to , and produces the observed positron excess. These decays can naturally take place via GUT interactions. If it exists, such a metastable particle can be found in the next run of LHC. The model predicts also the ratio of leptons over baryons in the universe to be close to .

1. Introduction

The possibility of dark matter being in the form of “dark atoms” has been studied extensively [121]. In this scenario, new stable particles are bound by new dark forces (like mirror partners of ordinary particles bound by mirror electromagnetism [2226]). However, it turns out that even stable electrically charged particles can exist hidden in dark atoms, bound by ordinary Coulomb interactions (see [2730] and references therein). Stable particles with charge −1 (and corresponding antiparticles as tera-particles [31]) are excluded due to overproduction of anomalous isotopes. However, negatively doubly charged particles are not constrained by anomalous isotope searches as much as −1 charged particles [32]. There exist several types of particle models where heavy stable −2 charged species, , are predicted:(a)AC-leptons, predicted as an extension of the Standard Model, based on the approach of almost-commutative geometry [3336];(b)technileptons and antitechnibaryons in the framework of Walking Technicolor (WTC) [3743]. All these models also predict corresponding charge particles. If these positively charged particles remain free in the early universe, they can recombine with ordinary electrons in anomalous helium, which is strongly constrained in terrestrial matter. Therefore a cosmological scenario should provide a mechanism which suppresses anomalous helium. There are two possible mechanisms that can provide a suppression.(i)The abundance of anomalous helium in the galaxy may be significant, but in terrestrial matter a recombination mechanism could suppress this abundance below experimental upper limits [33, 35]. The existence of a new gauge symmetry, causing new Coulomb-like long range interactions between charged dark matter particles, is crucial for this mechanism. This leads inevitably to the existence of dark radiation in the form of hidden photons.(ii)Free positively charged particles are already suppressed in the early universe and the abundance of anomalous helium in the galaxy is negligible [29, 44].These two possibilities correspond to two different cosmological scenarios of dark atoms. The first one is realized in the scenario with AC leptons, forming neutral AC atoms [35]. The second assumes a charge asymmetry of the which forms the atom-like states with primordial helium [29, 44].

If new stable species belong to nontrivial representations of the electroweak group, sphaleron transitions at high temperatures can provide the relation between baryon asymmetry and excess of −2 charge stable species, as it was demonstrated in the case of WTC [37, 4547].

After formation in the Big Bang Nucleosynthesis (BBN), screens the charged particles in composite () “atoms” [44]. In all the models of , behaves either as a lepton or as a specific “heavy quark cluster" with strongly suppressed hadronic interactions. Therefore interactions with matter are determined by the nuclear interactions of He. These neutral primordial nuclear interacting objects can explain the modern dark matter density and represent a nontrivial form of strongly interacting dark matter [4856].

The cosmological scenario of the universe can explain many results of experimental searches for dark matter [29]. Such a scenario is insensitive to the properties of , since the main features of the dark atoms are determined by their nuclear interacting helium shell. In terrestrial matter such dark matter species are slowed down and cannot cause significant nuclear recoil in the underground detectors, making them elusive in direct WIMP search experiments (where detection is based on nuclear recoil) such as CDMS, XENON100, and LUX [5761]. The positive results of DAMA and possibly CRESST and CoGeNT experiments [6266] can find in this scenario a nontrivial explanation due to a low energy radiative capture of by intermediate mass nuclei [29, 30].

It has been also shown [37, 4547] that a two-component dark atom scenario is also possible. Along with the dominant abundance, a much smaller excess of positively doubly charged techniparticles can be created. These positively charged particles are hidden in WIMP-like atoms, being bound to . In the framework of WTC such positively charged techniparticles can be metastable, with a dominant decay channel to a pair of positively charged leptons. In this paper we show that even a fraction of such positively charged techniparticles with a mass of 1 TeV or less and a lifetime of  s, decaying to , , and , can explain the observed excess of cosmic ray positrons, being compatible with the observed gamma-ray background.

One should note that, as it was shown in [35, 37, 44, 45] (for a review, see [29, 33] and references therein), the case of −2 charged stable particles is significantly different from the case of stable or metastable particles with charge −1, avoiding severe constraints on charged particles from anomalous isotope searches and BBN due to their catalytic effects (see, e.g., [6769]). In essence this difference comes from the fact that primordial He formed in BBN captures −2 charged particles in neutral states, while −1 charged particles are captured by in charged ions, which either (if stable) form anomalous isotopes of hydrogen or (if long-lived, but metastable) catalyze processes of light element production and influence their abundance. Nuclear physics of is in the course of development, but a qualitative analysis has shown [46] that the interactions with matter should not lead to overproduction of anomalous isotopes, while catalytic effects in BBN can lead to primordial heavy element production, but not to overproduction of light elements.

The paper is organized as follows. In Section 2 we give a brief review of dark atoms made of stable charged techniparticles. In Section 3 we present the constraints and the predictions of the scenario with respect to the parameters of the Technicolor model we use, as well as how the ratio of lepton over baryon number is deduced. In Section 4 we show what GUT operators can implement the decay of the doubly charged particle to leptons. In Section 5, we show how the scenario of decaying dark matter can be realized, and how it can explain the PAMELA and AMS-02 results while satisfying the Fermi/LAT constraints. We conclude in Section 6.

2. Dark Atoms from Techniparticles

Technicolor theories that do not violate the electroweak precision tests, while not introducing large flavor changing currents, have been extensively studied lately (see [70] and references therein). Old models where fermions transformed under the fundamental representation of the gauge group required a large number of flavors (for a given number of colors) in order to be close to the conformal window and thus to suppress the flavor changing neutral currents. The need for many flavors coupled to the electroweak sector (that violates the electroweak precision measurements) disfavored Technicolor in the past. However, it has been demonstrated that once one allows fermions to transform under higher representations of the gauge group, quasi-conformality can be achieved even with a small number of colors and flavors [3840]. This means that there is a set of Technicolor models that evade the strict constraints of the electroweak tests, making Technicolor a viable candidate for the TeV energy scale. Apart from the perturbative calculation of the oblique parameters [41] in this type of models, nonperturbative calculations based on holographic descriptions [7173] showed that indeed the oblique parameter can be small. Note that the oblique parameters (e.g., , , and ) measure the modifications of the Standard Model gauge boson vacuum polarization amplitudes caused by contributions of new physics. These parameters are severely constrained by electroweak precision tests. Extra flavors that couple with the electroweak sector contribute to these parameters and can potentially exclude a model.

One of the simplest models that possesses the features described above, is the so-called Minimal Walking Technicolor [38, 42, 74]. The theory consists of two techniquarks transforming under the adjoint representation of an gauge group, and an extra family of leptons and coupled to the electroweak in order to cancel the global Witten anomaly. The hypercharge assignment can be chosen consistently (without introducing gauge anomalies) such that one of the techniquarks has zero electric charge. Such a simple theory can have a variety of dark matter candidates, ranging from dark matter particles that are Goldstone bosons of the theory (with nonzero technibaryon number) [43, 75, 76] or Majorana WIMPs [7782]. Apart from these possibilities, there is another intriguing scenario that is of an electromagnetic bound state between a charged helium nucleus and a charged techniparticle [37, 45]. More specifically in [37], we examined the possibility where the dark matter bound state is or . Recall that and are the two techniquarks of the theory and and are the extra leptons. There is a gauge anomalous free hypercharge assignment where the charges of , , , and are, respectively, , 0, , and . We should also emphasize that, due to the fact that techniquarks transform under the adjoint representation of the gauge group, some of the Goldstone bosons are colorless diquarks (carrying technibaryon number). Apparently and have charges . This candidate (with being or ) is a Strongly Interacting Massive Particle (SIMP) rather than a WIMP due to the large geometric cross section of the helium component. Despite the large cross section, this candidate has not been ruled out by any experiment so far. Amazingly enough, although such a candidate interacts strongly with matter, it cannot be detected in earth based detectors (based on measuring the recoil energy) like CDMS, Xenon, or LUX. By the time such a particle reaches the detector, and it has lost most of its kinetic energy making it impossible to produce recoil energies above the detection threshold. In [45], we examined a generalized version of the aforementioned scenario, where although the majority of dark matter is (or ), a small component can be of the WIMP form (or ). Such a WIMP component must be small since it is constrained by direct detection experiments.

In [37, 45], we had assumed that techniparticles are stable. In particular with respect to the technibaryons, the symmetry associated with the technibaryon number protected the lightest diquark Goldstone boson from decaying. Here we reexamine the scenario of [45] allowing decays of the techniparticles. It has been demonstrated that decaying dark matter can provide a possible explanation of the unexpected positron excess seen in PAMELA [83, 84]. Decaying of dark matter particles through a dimension-6 operator gives a lifetime where is the mass of the dark matter particle. Note that we have normalized the lifetime with respect to a GUT scale by an order of magnitude lower than the typical value of suggested by supersymmetry. As we are going to argue a small component of dark matter with a mass of ~TeV or less and a lifetime of  s can accommodate nicely the positron excess seen in PAMELA and AMS-02 data. In addition such a lifetime is sufficiently large in order not to deplete the density of this component of dark matter by today since it is a few orders of magnitude larger than the age of the universe. As it was stressed in [83], dimension-6 operators are very natural objects in Technicolor, and therefore such a framework becomes very appealing.

3. Techniparticle Excess

We already mentioned that the MWT has two techniquarks and in the adjoint representation of the Technicolor with charges and 0 and two new leptons and with charges and , respectively. The theory possesses a global symmetry that breaks spontaneously to an . Out of the 9 Goldstone bosons, three of them (with the quantum numbers of the usual pions) are eaten by the and bosons, while the rest 6 are the colorless diquarks , , and and their antiparticles [43].

We are going to consider two possibilities. The first one is to have an excess of charge and a little of . The main component of dark matter is the SIMP . There is also a small WIMP component of . The second scenario is to have an excess of and a little of , in such a way that the main SIMP component of dark matter is and the small WIMP one is . In both cases we have assumed that is the lightest among the technibaryons and similarly is the lightest of the new leptons. The calculation of the relic density of the technibaryons taking into account sphaleron violating processes, weak equilibration, and overall charge neutrality gives similarly to [43] where , , , and are the technibaryon, baryon, lepton, and new lepton family number, respectively. are statistical factors for the specific particle given by where the functions and are defined as follows: is the freeze-out temperature for the sphaleron process, usually taken somewhere between 150 and 250 . In the first aforementioned possibility, the dark matter density is where , , and are the masses of , , and proton, respectively. We have taken the ratio of dark matter to baryonic matter to be . If denotes the fraction of the WIMP component () of dark matter, then the ratio of leptons over baryons is given by In the second scenario (that of and ), where is the mass of . The ratio is here is again the fraction of the WIMP-like component of dark matter. There are two points we would like to emphasize. The first one is that both possibilities give a ratio of lepton over baryon numbers very close to if the masses of and remain around 1 TeV. In fact the first scenario gives a ratio slightly above and the second gives a ratio slightly below. starts deviating (exponentially) as a function of the mass of and/or once we go to masses much higher than 1.5 to 2 TeV (see Figure 1). The second point we would like to stress is that is constrained by earth based direct detection search experiments. In [45] we found that the WIMP component of this dark matter scenario cannot be more than 1% (or ). Since then, the constraint from the CDMS and Xenon experiments has improved significantly and more severe constraints from LUX appeared. The cross section of (or its antiparticle) with a proton is [85] where , that is, the sum of the hypercharge of left and right components. It is easy to check that in our case . This is because has and and has and . The total sum is . In addition since both and are much heavier than the proton, the reduced mass is approximately the mass of the proton. The form factor depends on the target nucleus and the recoil energy. For example, for Ge detector with recoil energies between 20 and 50 keV, the form factor ranges from 0.43 to 0.72 [43]. Here in this estimate of the WIMP-proton cross section we have set . The results of the LUX experiment [61] exclude WIMPs with a cross section for a typical WIMP mass of 1 TeV. This means that WIMPs with the cross section of can make up only a component of or smaller of the total dark matter. Here we are going to use a typical value of .

4. Decaying Dark Matter

As we mentioned in the previous section, we might have a (or less) WIMP component in our dark matter framework. This comes in the form of (first scenario) or (second scenario). Our goal is to consider decay processes that can produce the excess of positrons seen in PAMELA and AMS-02. For this, it is generically better if the objects decay accordingly.

In the first considered scenario we assume that is stable, and therefore the SIMP component which consists the overwhelming part of dark matter is unaffected. On the other hand, we assume that can decay to leptons. By construction since and belong to the same electroweak doublet, couples to and . Since is a lepton with an electric charge , it can in principle slightly mix with the usual leptons, that is, electrons, muons, and taus. The tiny WIMP component of dark matter made of decays due to the fact that can decay to a and (via ) to positrons, antimuons, and antitaus. We assume that is heavier than , and therefore the decay is suppressed. In order not to get very fast decays of , the mixing of with positrons and so forth has to be extremely small. However, this is something expected due to experimental constraints as well as due to the fact that is much heavier than the leptons. It is also expected that the mixing between and would be larger than and or and positrons. The decay in this scenario can be accommodated via a dimension-5 operator. However, decays of to positrons, , or can lead to unwanted production of hadrons via decays of . Therefore we focus on the second case.

In the second scenario the small WIMP component is made of . In this case we assume that is stable (and no mixing with other leptons exists), but the Goldstone boson decays via a GUT interaction. A natural dimension-6 operator that can accommodate the decay can be of the form where is an electron, muon, or tau. Notice that, due to the transpose instead of the bar, such an operator violates both the lepton and the technibaryon number. It allows a possible decay of to two positrons (or two antimuons or antitaus (in principle we can have an even more general operator where decays to different species of antileptons, i.e., a positron, an antimuon, etc)): It is understood that although and can be generic Dirac matrices, has to be the charge conjugate matrix in order for to be the pseudo-Goldstone boson . If we require that parity is not violated by the interaction, must also be the charge conjugate matrix. In case parity is violated, can be (as it is a well-known fact that is a scalar). Of course nothing forbids a similar decay of to two quarks or even a quark and a lepton, as it would depend on the details of the GUT interaction. However, here we do not want to speculate regarding the GUT interactions but simply to demonstrate that such a realization can in fact produce the positron spectrum seen by experiments. As we already mentioned, a dimension-6 operator of the above form would give according to (1) a lifetime of the order of for a mass of of the order of TeV. If does not decay to hadrons, this scenario is more appropriate for explaining the positron excess compared to the first scenario we mentioned because in the first scenario the decay of will always be accompanied by hadronic decays that are not seen by PAMELA.

5. Positron Excess and Fit to the PAMELA and AMS-02 Data

Here we show the impact of decaying particles on the cosmic positron flux and diffuse gamma radiation. The so-called “PAMELA anomaly” in the cosmic positron spectrum [86] has been recently confirmed also by AMS-02 [87]. This anomaly cannot be explained by positrons of only secondary origin, and therefore primary positron sources are needed to explain the data. There are attempts to realize it based on decaying or annihilating dark matter models. Any scenario that provides positron excess is constrained by other observational data mainly from the data on cosmic antiprotons, gamma-radiation from our halo (diffuse gamma-background), and other galaxies and clusters [8895]. If dark matter does not produce antiprotons, then the diffuse gamma-ray background gives the most stringent and model-independent constraints.

In our scenario the component of a tiny WIMP component of dark matter decays as in principle with different branching ratios. All decay modes give directly or through intermediate particle decays the positrons and photons, which are hereafter referred as final state radiation (FSR). In Figure 2 we show the produced positron and gamma spectra for each decay mode individually. Note that, unlike PAMELA, the AMS-02 disfavors decays purely to and (although it does not exclude them).

In the context of indirect dark matter searches from cosmic rays (CR), the leptonic decay modes have been studied extensively (see, e.g., [8895]), using a variety of different approaches in estimating the CR signals. For our estimate, we adopt the following model of positron propagation in the galaxy. Due to energy losses, positrons have a finite diffusion length at given energy where is a typical value for the diffusion coefficient [96], is the rate of energy losses with , and is the initial energy. All energies are measured in . The effect of the diffusion in the propagation can be estimated by assuming a homogeneous distribution of the sources. In fact, the result of diffusion is not sensitive to the effects of inhomogeities, because it depends on the averaged density within the diffusion length. Since we are interested in positron energies above , which corresponds to  kpc (see (12)) over which no essential inhomogeneity effects are expected, this simple approximation we make here is good. At , secondary positrons dominate the spectrum. If exceeds the size of the magnetic halo (MH) ( kpc in height and  kpc in width), the leakage of particles from the halo should be taken into account. We consider this effect by introducing a suppression factor, which is equal to the ratio of the volume of MH contained within the sphere of radius : where is the step function. If is the number of positrons produced in a single decay (see Figure 2), the positron flux near the Earth can be estimated as where is the local number density of particles with . Recall that is the fraction of dark matter in the WIMP component.

The effect of solar modulation becomes important at the less interesting low energy part of the positron spectrum. To account for this effect, we have adopted the forced field model [98] with two different parameters for positrons and electrons. They are easily adjusted so they can fit the data points at low energy. The positron and electron background components were taken from [99]. In Figure 3 we present the positron excess due to decays for two values of the mass of ,  TeV and  TeV. We also show the lifetime of and the branching ratios that fit the experimental data optimally for each choice of . They evade the existing constraints of [8895].

The gamma-ray flux from decays has two main contributions: one from FSR (shown in Figure 2) and another one from Inverse Compton (IC) scattering of positrons on background photons (star light, infrared background, and CMB).

For the FSR photons produced by decays in our galaxy, the flux arriving in the Earth is given by where we use an isothermal profile , and are the distances from the Galactic center and the Earth, respectively. We obtain the averaged flux over the solid angle corresponding to , . For the IC photons from our galaxy, we have estimated the contribution following [100]. In Figure 4 we show both contributions in the gamma-ray flux for the same parameters as in Figure 3.

Decays of , which are outside our Galaxy being homogeneously distributed over the Universe, should also contribute to the observed gamma-ray flux. For FSR photons this contribution can be estimated as where corresponds to the recombination epoch, is the current cosmological number density of , is the inverse value of the Hubble parameter with being the age of the universe, and and are, respectively, the current vacuum and matter relative densities. Note in (16) the transition between distributions at different , . This extragalactic contribution to FSR increases significantly the total gamma-ray flux as shown in Figure 4 by dot-dashed lines.

It is not expected that extragalactic IC photons can contribute significantly to the spectrum. Indeed, mainly only low energetic CMB photons are present in the medium outside the galaxy (or before the galactic stage). After the scattering of electrons with energy off CMB photons with energy  eV, the recoiled photons acquire at redshift energy , which is below 1  in the modern epoch. It makes therefore this contribution indifferent for the energy range of Fermi/LAT.

To conclude, on the basis of Figure 4, one may assert that the considered scenarios of decays satisfy the Fermi/LAT constraints. In addition, although we used the best fit values for the branching ratios, we have found that some small variation of the branching ratios is possible. If one chooses  TeV, a possible satisfaction of the constraints is possible at the expense of the positron spectrum fit.

6. Conclusions

Dark matter can potentially be in the form of neutral dark atoms made of stable heavy doubly charged particles and primordial He nuclei bound by ordinary Coulomb interactions. This scenario sheds new light on the nature of dark matter and offers a nontrivial solution for the puzzles of direct dark matter searches. It can be realized in the framework of Minimal Walking Technicolor, in which an exact relation between the dark matter density and baryon asymmetry can be naturally obtained predicting also the ratio of leptons over baryons in the universe. In the context of this scenario a sparse component of WIMP-like dark atoms of charged techniparticles can also appear. Direct searches for WIMPs put severe constraints on the presence of this component. However, we demonstrated in this paper that the existence of a metastable positively doubly charged techniparticle, forming this tiny subdominant WIMP-like dark atom component and satisfying the direct WIMP searches constraints, can play an important role in the indirect effects of dark matter. We found that decays of such positively charged constituents of WIMP-like dark atoms to the leptons , , and can explain the observed excess of high energy cosmic ray positrons, while being compatible with the observed gamma-ray background. These decays are naturally facilitated by GUT scale interactions. This scenario makes a prediction about the ratio of leptons over baryons in the universe to be close to . The best fit of the data takes place for a mass of this doubly charged particle of 1 TeV or below making it accessible in the next run of LHC.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.


C. Kouvaris is supported by the Danish National Research Foundation, Grant no. DNRF90.