Moduli induced cogenesis of baryon asymmetry and dark matter

We study a cogenesis mechanism in which the observed baryon asymmetry of the universe and the dark matter abundance can be produced simultaneously at low reheating temperature without violating baryon number in the fundamental vertex. In particular, we consider a model which could be realized in the context of type IIB large volume string compactifications. The matter superfields in this model include additional pairs of color triplet and singlet superfields in addition to the Minimal Supersymmetric Standard Model (MSSM) superfields. Assuming that the mass of the additional singlet fermions is O(GeV) and color triplet fermions is O(TeV), we show that the modulus dominantly decays into the additional color triplet superfields. After soft supersymmetry (SUSY) breaking, the lightest eigenstate of scalar component of color triplet superfield further decays into fermionic component of singlet superfield and quarks without violating baryon number. Assuming R-parity conservation, it follows that the singlet superfield will not further decay into the SM particles and therefore it can be considered as a stable asymmetric dark matter (ADM) component. We find that the decay of the lightest eigenstate of scalar component of color triplet superfield gives the observed baryon asymmetry in the visible sector, an asymmetric dark matter component with the right abundance and naturally explains cosmic coincidence.


I. INTRODUCTION
In cosmology some of the important puzzles are related to the origin of baryon asymmetry of the universe and the nature of dark matter. The comparable values of dark matter density and baryon density [1] Ω DM h 2 0 ∼ 5 Ω B h 2 0 , points to the possibility that they might have a common origin. However, the standard paradigm adopts completely different mechanisms to explain observable baryon asymmetry of the universe and dark matter abundance. The baryon asymmetry is generated from an initially baryon-antibaryon symmetric universe by considering baryon number (B), C and CP violating processes that went out of equilibrium in the early universe, while the dark matter density is produced by considering weakly interacting massive particles (WIMPs) (with mass around O(100) GeV) with the relic density being determined by the freeze out condition. The fact that they have a comparable abundance is often referred to as the "cosmic coincidence" puzzle. Recently, the CDMS collaboration has reported an excess in the dark matter events [2] which sets an upper limit of O(10 −41 ) cm 2 on the value of spin-independent (SI) dark matter-nucleon cross section for dark matter mass around 10 GeV at 3.1σ significance level. The excess reported by the CoGeNT collaboration [3] also hints at a light dark matter mass, almost in the same region of parameter space. The data taken by the XENON100 experiment [4] also gives a very stringent constraint on SI dark matter-nucleon cross section which points towards a dark matter mass around O( GeV). The light dark matter is often also motivated * Electronic address: mansi@prl.res.in † Electronic address: chandan@prl.res.in ‡ Electronic address: utpal@prl.res.in due to the possibility of explaining 3.5 KeV X-ray line by radiative decay of O( GeV) neutral dark matter particle [5]. However, for an O( GeV) mass the thermal WIMPs give over-abundance of dark matter particle for annihilation cross-section less than 10 −26 cm 2 , and thus the alternative schemes where an O( GeV) mass dark matter can be accommodated have gained significant attention. To this end, the cogenesis scenarios are particularly interesting because they have an attractive feature of explaining the observed baryon asymmetry of the universe together with an asymmetric dark matter component which can naturally satisfy the criterion for O( GeV) mass dark matter. Furthermore, the apparent coincidence of the baryon and dark matter densities can also be addressed in such a framework using the underlying connection between the baryogenesis scenarios and dark matter production. There exists several different mechanisms in the literature [6], which address simultaneous generation of baryon (or lepton) asymmetry and the asymmetric dark matter abundance. The cogenesis of both without violating B or B − L is discussed in Refs. [7].
From a top-down model building perspective, an UV completed supersymmetric model may entail the existence of WIMPs, known as the modulus (moduli). In N = 1 supergravity limit of string theory, the moduli appear while compactification of the extra dimensions takes place [8]. The decay of moduli fields have significant implications for the cosmological history of the universe [9]. The entropy released due to late decay of the lightest modulus dilutes the existing baryon asymmetry of the universe as well as the relic abundance of dark matter produced at high scale. However, the correct amount of dark matter can be produced non-thermally from the decay of modulus into the lightest supersymmetric particle. The non-thermal realization of dark matter is discussed in Refs. [10][11][12][13] in the context of Minimal Supersymmetric Standard Model (MSSM) and string-motivated mod-arXiv:1508.04144v1 [hep-ph] 17 Aug 2015 els. Given that the decay of heavy modulus leads to very low reheating temperature, it renders electroweak baryogenesis and leptogenesis impossible. However, it is possible to accommodate direct baryogenesis and correct dark matter abundance by considering late-decaying moduli in the schemes implementing presence of additional color triplet superfields along with MSSM superfields [14][15][16], and with implementation of other mechanisms [17,18]. The coincidence problem has also been addressed by considering Affleck Dine (AD) baryogenesis in the presence of moduli in Refs. [19,20]. Interestingly, the decay of lightest string modulus into its superpartner axion can also explains the source of dark radiation [21][22][23].
In this work, we propose a model for moduli induced cogenesis which simultaneously generates the baryon asymmetry of the universe and an asymmetric dark matter (ADM) component with a dark matter mass around 5 GeV. In this model, the particle content include two additional iso-singlet color triplet superfields χ andχ with hypercharges −4/3 and 4/3 respectively and two singlet superfields N andN 1 . Due to the presence of the pair of the additional superfields with opposite hypercharges, the terms analogous to the Giudice-Masiero term [24] in the Kähler potential dictate the decay width of modulus into both colored and singlet superfields. In N = 1 supergravity, the effective supersymmetric mass terms as well as soft SUSY breaking terms depend on coupling strength of the hidden sector field (modulus) to the visible sector fields [25]. Interestingly, the effective masses of additional colored(singlet) superfields are also governed by the same Giudice-Masiero like term(s) considered in the kähler potential. Therefore, the coefficient of interaction term responsible for the decay of modulus into the pair of colored(singlet) superfields i.e. the coefficient of new Giudice-Masiero like term(s) can be constrained based on given masses of the superfields. In this model the mass of color triplet superfields being heavier as compared to the mass of singlet superfields, the modulus would preferably decay into pair of color triplet superfields. Now the scalar component of color triplet superfields further decay into quarks and additional singlet fermions and the baryon number of the color triplet superfield gets distributed between quark and additional singlet fermion. The conservation of R-parity ensures that the singlet fermion will further not decay into the Standard Model (SM) particles and therefore can be considered as dark matter component. The decay process is baryon number conserving at tree level as well as at oneloop level, however, it is CP asymmetric at one-loop level due to the presence of soft SUSY breaking terms. Consequently, an asymmetry is generated in both the visible and the dark matter sector. We find that the symmet-ric component of dark matter gets annihilated for a dark matter mass O( GeV), and the required order of baryon asymmetry and dark matter relic abundance can be successfully generated in this mechanism for certain values of Yukawa couplings.
The outline of the rest of the paper is as follows. In section II, we briefly describe a phenomenological model that can be obtained as a low energy limit of the large volume scenario (LVS) proposed in the context of type IIB string compactifications. In section IIA, we discuss all possible decay modes of the modulus and give the corresponding decay widths showing that the modulus dominantly decays into pair of Higgs, axions and additional color triplet and singlet superfields introduced in the low energy spectrum, depending on the coefficients of the interaction terms. In section II B, we show that due to the heavier mass of color triplet superfields as compared to singlet superfields, the modulus preferably decays into pair of color triplet superfields. In section III, we propose a new mechanism which can simultaneously generate the observed baryon asymmetry and the dark matter abundance without violating baryon number in the fundamental vertex. In section III A, we discuss the possible decay mode of scalar component of color triplet superfield into quarks and additional singlet fermions, generating a baryon asymmetry. Next we argue that in a R-parity conserving scenario, the singlet fermion will further not decay into the (MS)SM particles and therefore can be considered as a stable dark matter candidate. In section III B, we discuss a mechanism to annihilate the symmetric component of dark matter, leaving only the asymmetric component accounting for the dark matter relic abundance. In section IV, we summarize our results and conclude.

II. A PHENOMENOLOGICAL MODEL BASED ON LARGE VOLUME SCENARIO AND MODULUS DECAY
The presence of gravitationally coupled moduli fields can have significant impact on the standard cosmology. During inflation i.e. when the Hubble expansion rate H inf >> m Φ , modulus (Φ) gets significantly displaced from the minimum of its potential [26]. Thus, if one takes into consideration the presence of modulus and high scale inflation, it is a rather generic consequence to expect the modulus to be displaced from the low-energy minimum by an amount |∆Φ| = | Φ inf − Φ 0 | ≈ M P . Since the energy density of these oscillations dilutes in the same way as non-relativistic matter, they will come to dominate the expansion of the universe. This will continue until the modulus decays at a time t ∼ Γ −1 Φ , transferring the remaining oscillation energy into radiation, hence reheating the universe at a late time. The reheating temperature after the modulus decay is given by [9], where g * is the total number of relativistic degrees of freedom and Γ Φ is the decay width of the modulus. The number density of any particle X produced from the decay of the modulus is given by It has been a challenging task to obtain a realistic low energy spectrum in string compactifications. The foremost step while constructing reals in string compactifications is the issue of moduli stabilization. A realistic model should be able to realize de-Sitter minima and also avoid Cosmological Moduli Problem (CMP) [27]. The Large Volume Scenario (LVS) has been considered as an ideal framework to build consistent MSSM-(like) chiral model in which soft terms are calculated explicitly. Let us consider the Large Volume type IIB compactification scheme initially proposed in [28,29]. The volume of the Calabi-Yau (CY) manifold is of swiss-cheese type, given by where τ B denotes the big divisor volume modulus which mainly controls the size of the CY volume and τ α 's correspond to small divisor or blow-up moduli fields. The volume moduli are complexified by associating them with four-form axions. The spectrum also takes into consideration dilaton (S), complex structure moduli (U) and two-form axions. The Kähler potential of the effective theory includes kähler moduli and perturbative α corrections. The superpotential includes nonperturbative contribution effects on the small blow-up mode. Interestingly, the (non)perturbative corrections in the effective potential lead to the nonsupersymmetric anti-di Sitter minima when the volume of the Calabi-Yau manifold is very large [28]. The visible/MSSM sector in this scenario is realized by including D3/D7 branes on blow-up modes [30]. The soft SUSY breaking terms are calculated from the nonvanishing F-terms corresponding to the hidden-sector moduli. The pattern of the softterms depend on the location of the D-brane in the bulk geometry. If the MSSM-like divisor is placed in the proximity of main source of SUSY breaking sector, it will give soft term masses of the order of gravitino mass. If location of D-brane is geometrically separated from dominant SUSY breaking sector, it will give smaller mass of sparticles as compared to the gravitino mass, known as sequestered string models.
In this work, we follow the sequestered Large Volume Models discussed in Refs. [31,32]. The relevant scales in this model are given by The soft terms as well as supersymmetric mass terms are evaluated by expanding the Kähler potential and super-potential as a power series expansion in the matter superfields respectively. The analysis of soft terms in Ref. [32] is given in two limits, namely (i) local limit (ii) ultralocal limit. The classification is based on the precise form of the Kähler metric which is used to obtain Yukawa couplings independent of the compactification volume. Since the analysis of the soft terms also depends on that particular form of Kähler metric, it generates different pattern of soft SUSY breaking terms in different limits. In this paper, we consider ultra local limit in which The requirement of m soft ∼ O( TeV) constrains the value of CY volume V ∼ O(1) × 10 7 in string length units. The choice of V ∼ 10 7 also provides 60 e-folds of inflation, generating right amount of density perturbations in this model [33].
We assume that the blow-up mode upon which visible sector is realized by wrapping D3/D7-brane, includes aforementioned new set of (s)particles in addition to the MSSM spectrum. Therefore, the matter Kähler potential for our model includes the soft SUSY breaking terms corresponding to additional superfields also.

A. Modulus decay
There are multiple moduli present in the LVS model described in Refs. [31,32], however only the lightest modulus couples to the (MS)SM as well as the additional fields. Below we give the kinematically possible decay modes of the lightest "big" divisor volume modulus into different modes, studied extensively in Ref. [21].
The Kähler potential involving the "big" divisor volume is given by with T B = τ B +ia B , where τ B is the real volume modulus and a B is four-form axion. This leads to the interaction terms Canonically normalizing the volume modulus Φ = 3 2 ln τ B , one obtains (i) Decay into axions: Utilizing the above, the decay width for the modulus decaying into axions is given by (ii) Decay into gauge bosons: The coupling of the modulus to the gauge bosons is obtained through gauge kinetic function. In LVS model, the SM arises from wrapping D-branes on the blown-up modulus, so there is no direct tree-level coupling of the volume modulus τ B to gauge bosons (A µ ). The effective interaction term appearing at one-loop level is given by This leads to decay width given by (iii) Decay into MSSM Scalars: The couplings to matter scalars are given by After canonical normalization, the interaction term is given by The decay width for the matter scalars is given by (iv) Decay into matter fermions and gauginos: Using equation (10), the coupling of matter fermions and gauginos to the volume modulus is given by It leads to a decay width (v) Decay into the Higgs: The coupling of the modulus to the Higgs fields is dominated by the Guidice-Maisero coupling in the Kähler potential After canonical normalization, the decay width is given by If the Higgs possess a shift symmetry [34], then z = 1 and decay width is proportional to m 3 Φ /M 2 P .
(vi) Decay into color triplets: We consider an analogue of Guidice-Maisero term for additional color triplets superfields under consideration. So, the couplings of the modulus to the pair of color triplets is given by where z χ is an undetermined constant. After canonical normalization of the volume modulus, the interaction term is given by The decay width for color triplets is given by (vii) Decay into singlets: Similar to the case of color triplets, the coupling of the modulus to the pair of singlets is given by where z N is an undetermined constant. After canonical normalization of the volume modulus, the interaction term is given by The decay width for the modulus decay into singlets is given by Since m soft , m gaugino << m Φ , it is clear from equation (12) and (14) that couplings of the modulus into (s)particles as well as gauginos are suppressed. The modulus can dominantly decay into any of the pair of Higgs, axions, color triplets χ and singlets N depending on the value of the coefficients of the interaction coupling of the modulus to the same. In the next subsection, we show that the coefficient of the interaction term coupling the modulus to the pair of color triplets χ and singlets N , e.g. z χ and z N , can be constrained based on the masses of these particles.

B. Constraints on the modulus interaction to color triplet (singlet) superfields
In N = 1 supergravity limit of any superstring model, the effective mass term of any matter superfield obtains contribution from both superpotential as well as Kähler potential, and the values of both supersymmetric mass term and soft SUSY breaking parameter corresponding to any matter superfield depend on the interaction of the hidden sector field with the matter superfields [25]. Hence, one can naively expect that the coefficient of the interaction coupling of the hidden sector field to the matter superfields can be constrained depending on the mass of the matter superfields. We explicitly describe this situation in the context of sequestered large volume compactification model below.
The soft terms as well as supersymmetric mass terms for the MSSM have already been calculated in the context of sequestered LVS model in Refs. [31,32]. Similar to the MSSM superfields, the soft SUSY mass terms of additional superfields under consideration will also depend on their interaction with moduli (hidden sector fields). Including additional color triplet superfields (χ,χ) and singlet superfields (N ,N ), the matter superpotential is given by The Kähler potential corresponding to matter superfields is given by where and C i , i = 1, 2, 3 correspond to the three generations of MSSM superfields. The general expression for supersymmetric mass term associated with superfield χ is given by where the index I corresponds to the number of moduli. The first term in the above expression appears from the superpotential while the second and third terms appear from the Kähler potential. Similar to the higgsino mass parameter µ in the MSSM, M χ is just an input parameter which can be made to vanish under general assumptions [35]. Following Refs. [31,32], mass term M χ similar to the higgsino mass parameter µ in the MSSM will be proportional to 2 The soft SUSY breaking Lagrangian for additional superfields is given by The soft scalar mixing term associated with superfield χ is given by Following Refs. [31,32], the value of B χ M χ (similar to the soft SUSY breaking term corresponding to the higgsino mass parameter (Bµ) in the MSSM) is given by The mass matrix in the (χ,χ) basis can be diagonalized by the transformations In the diagonalized basis, the mass eigenvalues for the scalar component of χ are given by Due to universality of the soft terms massesm 2 χ =m 2 χ = m 2 soft , the above equation reduces to For z χ = 1, one obtains B χ M χ =m 2 soft . Hence, the lightest mass eigenvalue of scalar component of χ is the same as the fermion mass M χ . Thus both scalar and fermionic components of χ will be of O( TeV) if z χ = 1.
Similarly, we have The mass eigenvalues for the scalar component of N is given by Form 2 N ∼ m 2 soft and B N M N ∼ 10 −3 m 2 soft , the eigenvalues of the scalar component of N will be O( TeV) even for z N ∼ 10 −3 . By using equations (19) and (22), it implies that which clearly shows that the modulus will dominantly decay into the pair of color triplet superfield (χ,χ) as compared to the pair of singlet superfields (N ,N ).

III. COGENESIS MECHANISM
In this section, we discuss the cogenesis mechanism in which both the baryon asymmetry as well as the dark matter abundance are produced via baryon number conserving decay of a singlet superfield N . As mentioned earlier, we are consider the mass of the fermionic component of (N ,N ) to be of the order of GeV, while the masses of fermionic components of (χ,χ) to be of the order of TeV. The modulus decaying into any of these particles as well as the other particles will dilute the preexisting baryon asymmetry as well as the dark matter density produced at high temperature. As a result, the baryon asymmetry and the dark matter density have to be created via some other mechanism at low temperature. It follows from the results for the decay of the modulus into different species (discussed in the previous section) that the modulus can dominantly decay into the pair of Higgs, axions and color triplets. We are interested in a scenario where the baryon asymmetry can be produced at a low temperature through the decay of χ and in the process, one ends up with an asymmetric component of dark matter with a mass O( GeV) giving the correct relic density 3 .

A. Baryon asymmetry and asymmetric dark matter
In this subsection, we describe a mechanism of baryogenesis where the fundamental interactions conserve baryon number, and the baryogenesis happens by generating certain amount of asymmetry in both the visible and the dark sector. We show that the decay products of the lightest eigenstate of χ − (χ − ) can simultaneously explain the observed baryon asymmetry and give rise to an asymmetric dark matter component, if N is light (mass O( GeV)). The ratio of dark matter abundance to baryon abundance is given by The cosmic coincidence Ω ADM /Ω B ∼ 5 is satisfied if Y ADM ∼ Y B ∼ 10 −10 for a dark matter mass M N around 5 GeV. In the previous section, we have seen that modulus decays dominantly into pair of χ andχ superfields. A cartoon showing the decay of the modulus into color triplet superfields which further decay into a singlet and a quark is given in Figure 1. Using the interaction term and B = −1 respectively. Therefore, the interaction of χ with N andū conserves baryon number. The scalar and fermionic components of χ andχ have R-parity assignments +1 and -1 respectively, while the scalar and fermionic components of N andN have R-parity assignments -1 and +1 respectively. Now R-parity conservation implies that the fermionic component of N (N ) will further not decay into the SM particles. Consequently, the fermionic component of N superfield produced during decay ofÑ − can be considered as a stable asymmetric dark matter particle. Thus, it follows that the decay of a scalar component of N generates an equal and opposite amount of baryon asymmetry in the visible sector and the dark sector i.e. Let us first calculate the baryon asymmetry generated in the visible sector. We begin by considering a single generation of singlet superfield N . The CP violation arises due to the soft SUSY breaking terms. The degeneracy between the two states belonging to the supermultiplet of a same generation can be removed by including the SUSY breaking effects and the CP violation occurs due to the interference between the two states of a single generation [36,37], as compared to conventional baryogenesis mechanism where at least two generations are required for CP violation. The SUSY breaking Lagrangian involving χ andχ superfields is given by where indices i, j, k correspond to the different generations of the particles. The evolution of system governing χ − χ † mixing is given by [36,37] χ|H|χ † = M χ(12) − i 2 Γ χ (12) .
This induces a mass difference ∆M χ and a decay width difference ∆Γ between two states given by The tree level and the one-loop diagrams responsible for generating CP violation are shown in Figure 2. The interaction Lagrangian in the mass basis (χ + ,χ − ) is given by, (43) Now, if we start with equal densities ofχ andχ † at t = 0, then after some time t the states evolve intõ where f + (t) = e −iMχt e −Γχt/2 cos(∆M χ t/2), Now the excess of N ,ū over N c ,ū c is given by The decay width for the decay modes ofχ(t) andχ(t) † is given by 4 where X 1 is the normalization factor and M uN is the amplitude of the decay mode considered. Incorporating these expressions into equation (46), the asymmetry parameter is given by [36,37] .
The decay width for χ − decaying into Nū is given by Γ χ = κκ † 4π M χ . The asymmetry is given by Using equation (II) we have B χ M χ = M χ ∼ O( TeV), and κ 1, implying Γ χ < B χ M χ . In the small κ limit, the above equation reduces to The baryon asymmetry generated from the modulus decay is given by Since the Yukawa couplings do not depend on the compactification volume in this model, we choose κ ∼ 10 −2 . For A ijk ≡ κ ijk A ijk ∼ O( TeV), we get B = 10 −3 . Now using equation (2), the value of Y χ is given by For Br χ = 1, T R ∼ O( GeV) and m Φ ∼ 5 × 10 6 GeV, we get Y χ ∼ 10 −7 . Finally, from equation (53) the baryon asymmetry is given by Y B ∼ 10 −3−7 = 10 −10 , which is the observed baryon asymmetry of the universe. From equation (39), the asymmetry in the dark matter sector is also given by Y DM ≡ Y N ∼ 10 −10 .

B. Annihilation of the symmetric dark matter component
The requirement that an overall dark matter abundance is asymmetric follows from the condition that the symmetric component of dark matter should get depleted. For this to happen, the annihilation cross section of the symmetric component should be higher than the value of the cross section at freeze out temperature [12]. The freeze out temperature can be estimated by the rule of thumb given by Γ = n σ |v| = H.
(55) similar to the case of WIMP dark matter. However, in this case the reference temperature is the reheating temperature after the decay of the modulus. The Hubble expansion rate is given by H = The decay of the modulus produces a baryon number symmetric component of dark matter with a density given by n = Br χ 3T R 4mΦ × s, where s is the entropy density. Using equation (55) for T R = O( GeV) and s = 2π 2 45 g * T 3 R , the thermally averaged cross section is given by For the modulus mass m Φ = 5 × 10 6 GeV and g * = O(100), we get σ |v| freeze out ∼ O(10 −14 ) GeV −2 .
The annihilation of N and N c can be mediated through electroweak neutral Z boson. The cross section for N N c → Z → ff is given by σ |v| ∼ 1 4π where g corresponds to the gauge coupling of neutral Z boson to the pair of singlet superfields (N ,N ), g corresponds to the gauge coupling of Z boson to a pair of fermions and M Z is the mass of Z boson. For g ∼ 1, M N = 5 GeV, we get σ |v| ∼ 10 −8 GeV −2 , which is higher than the freeze out cross section calculated above. Thus, the symmetric component of the abundance of singlet fermion N gets annihilated and the left-over asymmetric relic abundance gives the correct dark matter relic abundance.

IV. CONCLUDING REMARKS
We have discussed a cogenesis mechanism unifying the generation of both baryon asymmetry of the universe and dark matter abundance in a model which can be obtained from a low energy limit of type IIB LVS string model. The decay of the lightest modulus (generically present in string models) dilutes the pre-existing baryon asymmetry and dark matter abundance of the universe. To avoid the cosmological moduli problem (CMP) one requires a very heavy modulus, which decays post inflation to give a low reheating temperature. Therefore one needs to consider a post reheating mechanism for generating the baryon asymmetry as well as the dark matter abundance. In this work, we show that both the baryon asymmetry and the non-thermal dark matter abundance can be generated simultaneously from the decay of a pair of color triplets produced after reheating. We consider the mass of the fermionic component of the pair of color triplets and singlets to be 1 TeV and 5 GeV respectively. We demonstrate that in the context of N = 1 supergravity, the interaction coupling of the modulus to the pair of singlet superfields as well as the colored superfields can be constrained depending upon the masses of the same. We find that the branching ratio of the modulus decaying into the pair of additional singlets is suppressed by a factor of 10 −6 as compared to the decay of the modulus into the pair of color triplets. Therefore, we conclude that the modulus will dominantly decay into pair of color triplets. The lightest eigenstate of scalar component of the color triplet further decays into singlet fermion and up type quark. Imposition of R-parity conservation ensures that the singlet fermion does not further decay into the SM particles and therefore it can be considered as a stable dark matter candidate. The CP asymmetry is generated via the interference of tree level and one loop diagrams for the decay of color triplets in the presence of soft SUSY breaking terms. We find that it is possible to obtain the observed baryon asymmetry of the universe and the asymmetric dark matter abundance by considering dark matter mass around 5 GeV, and the cosmic coincidence is natural in this scenario. Thus, if TeV scale colored fields are found at the LHC, it will have a very profound consequences for explaining the observed baryon asymmetry of the universe, the dark matter abundance and the cosmic coincidence.