Baryogenesis, Dark Matter, and Flavor Structure in Non-thermal Moduli Cosmology

The appearance of scalar/moduli fields in the early universe, as motivated by string theory, naturally leads to non-thermal"moduli cosmology". Such cosmology provides a consistent framework where the generation of radiation, baryons, and dark matter can occur while maintaining successful Big Bang Nucleosynthesis and avoiding the cosmological moduli problem. We present a relatively economical construction with moduli cosmology, building on a variety of string-inspired components (e.g. supersymmetry, discrete symmetries, Green-Schwarz anomaly cancellation). We address a range of outstanding problems of particle physics and cosmology simultaneously, including the fermion mass hierarchy and flavor puzzle, the smallness of neutrino masses, baryogenesis and dark matter. Our setup, based on discrete $\mathrm{Z}_{12}^{R}$ symmetry and anomalous $\mathrm{U}(1)_A$, is void of the usual issues plaguing the Minimal Supersymmetric Standard Model, i.e. the $\mu$-problem and the overly-rapid proton decay due to dimension-4,-5 operators. The model is compatible with $\mathrm{SU}(5)$ Grand Unification. The smallness of Dirac neutrino masses is automatically established by requiring the cancellation of mixed gravitational-gauge anomalies. The decay of the moduli field provides a common origin for the baryon number and dark matter abundance, explaining the observed cosmic coincidences, $\Omega_{B} \sim \Omega_{DM}$.


Introduction
Supersymmetry (SUSY), the maximal possible extension of the Poincare symmetry, remains one of the most appealing solutions to the gauge hierarchy problem of the Standard Model (SM). Among its salient features are the gauge coupling unification and the natural existence of dark matter (DM) candidates, the neutralinos. It is a difficult task to construct consistent cosmological scenarios based on standard thermal production of baryons as well as dark matter, allowing one to address the "cosmic coincidence puzzle" why their observed abundances are similar [1], i.e. Ω DM ∼ Ω B , as well as the matter-antimatter asymmetry. A non-thermal cosmology with an intermediate matterdominated epoch can be a viable setting to resolve these issues. From string theory, which generically requires supersymmetry at some scale for consistency, one expects existence of multiple scalar fields (moduli) with flat potentials. They acquire mass through non-perturbative effects after SUSY breaking. During inflation moduli fields develop large values and come to dominate the energy density of the universe as they coherently oscillate near the bottom of their potential. This naturally leads to a non-thermal "moduli cosmology". However, if the moduli are very light, their late-time decay on-sets radiation-dominated universe at energies below the scale allowed by the Big Bang Nucleosynthesis, leading to the "cosmological moduli problem" [2,3]. Making the moduli massive (i.e. "stabilizing moduli") [4][5][6][7][8] allows to avoid this. On other hand, massive moduli decaying to gravitinos can result in "moduli-induced cosmological gravitino problem" [9,10].
An attractive model with non-thermal moduli cosmology that addresses baryon-DM coincidence as well as baryogenesis, based on discrete Z 9 symmetry, was proposed by Kitano, Murayama and Ratz (KMR) [11] (for earlier related works see [12][13][14][15][16]). Here, a heavy scalar φ field, which is the scalar component of a chiral supermultiplet Φ = (φ,φ, F φ ), dominates the energy density of the early universe during its coherent oscillations and then later decays. The resulting radiation influx from decay washes out any dangerous relics that are present, such as the gravitinos. Enforced R-parity in the model prevents direct φ decays into gravitino pairs, thus avoiding the possible gravitino problem. On the other hand, through the operator ΦU DD present in the model, the decay of the φ field produces baryons as well as superpartners, for example through φ → qqq. The decay transfers any initial φ-number asymmetry into baryon number, explaining the baryon asymmetry. Additionally, the DM candidate, which is a Wino or Higgsino, is also produced by the φ decays in an appropriate amount. Since these processes happen right before the BBN, the problem of overproducing the dangerous relics is avoided. Since DM and baryon asymmetry have a common origin, the "cosmic coincidence puzzle" is naturally addressed in this setup (see [15,[17][18][19][20] for some alternative proposals). While the mechanism put forth by KMR produced a compelling cosmological setting, a variety of significant problems, particularly those related to particle physics, remained unsolved. Among the outstanding issues is the µ problem of the Minimal Supersymmetric Standard Model (MSSM), the flavor puzzle and neutrino mass generation. In this work we aim to construct a more complete model for non-thermal moduli cosmology, while addressing several of the major issues related to particle physics and the Standard Model. In particular, we incorporate the attractive KMR mechanism to produce a common origin for baryogenesis and DM. We address the µ problem of the MSSM and the issue of fermion masses, including the neutrinos. As has been shown in [21], only R symmetries can forbid the µ term. This motivates us to base our construction on an anomaly free (ensured by the 4-dimensional Green-Schwarz mechanism [22]) discrete R-symmetry. Further, we employ anomalous U (1) A flavor symmetry (see e.g. [23][24][25][26][27]) to address the 'origin of fermion mass hierarchy and mixing and establish the desired flavor structure through the Froggatt-Nielsen (FN) mechanism [28]. The interplay between U (1) A and discrete symmetries play an essential role in suppressing dangerous proton decay operators in the MSSM. As we Further, as we will show, the size of the smallness of the Dirac neutrino masses is closely related to the mixed U (1) A -gravitational anomaly cancellation. Our model is compatible with SU(5) Grand Unification (GUT). Since anomalous U(1) A as well as discrete R symmetries are expected to appear in string constructions, our setting is well motivated from a more fundamental theory.
This manuscript is structured as follows. Section 2 describes our SU(5)-compatible model based on a discrete Z R 12 symmetry combined with a anomalous U(1) A flavor symmetry. We discuss Green-Schwarz anomaly cancellation and review the general features of the model. In Section 3, we promote U(1) A to a full flavor symmetry. We comment how modifying the charge assignment of the flavor model proposed in [29], allows to cancel the gravitational anomaly without extra singlets and naturally explain the size of the Dirac neutrino masses in our model, in addition to explaining the top-bottom hierarchy in the Yukawa sector. The textures we obtain are shown to be fully compatible with SU(5) and correctly reproduce the mixing in the quark and lepton sector, as shown by performing a χ 2 fit. In Section 4, we discuss constraints from proton decay and comment on neutron-anti-neutron oscillations. We then demonstrate in Section 5 that the co-genesis of baryon number asymmetry and dark matter abundance is possible in our specific setting. Finally, Section 6 concludes the paper.

Symmetries and particle content
The model is based on a generation-independent Z R 12 R-symmetry in combination of an anomalous U (1) A flavor symmetry. The particle content of the model with the charge assignments of the superfields is given in Table 1. From the charge assignment, it is evident that our model is SU(5)-compatible with 10 = (Q, U , E) and 5 = (D, L) being the SU(5) chiral GUT super-multiplets unifying the matter fields in each SM generation, while Φ and Φ being the SU(5) singlet chiral superfields. W hid is the "hidden" sector superpotential and X, S and σ are the spurion fields parametrizing the SUSY/R breaking, the dilaton and the flavon, respectively. The "→ Z 3 " notation represents the residual Z 3 symmetry of Z R 12 after SUSY/R breaking.
The gauge invariant superpotential, up to order 4, is given by, where Y u, d, e, ν , implicitly depend on powers of the flavon field, σ. After the scalar component of the flavon field, σ, acquires a vacuum expectation value (VEV) breaking the U (1) A flavor symmetry, the effective Yukawa couplings are generated, as discussed in Section 3.1. Note that Φ n and Φ n are forbidden in the superpotential to all orders in n. Additionally, there will be terms in the effective superpotential that appear "non-perturbatively" after SUSY breaking, The order parameter for SUSY/R-symmetry breaking is the VEV of "hidden sector" superpotential W hid ∼ m 3/2 as in gravity mediation, which allows to parametrize the size of the effective terms.
Some key features of the model are described in the following: • µ-term: The µ term is forbidden in the superpotential due to the discrete R-symmetry, and appears only after SUSY and R-symmetry breaking. It can be generated with the correct magnitude through the Giudice-Masiero mechanism [30]. Specifically, the µ term can effectively appear from the Kähler potential once the anti-holomorphic X † field acquires a VEV in the F -term, F X ∼ m 3/2 M pl , with M pl = 2.435 × 10 18 GeV being the reduced Planck mass. This leads to µ ef f to be of the order of the gravitino mass, as in gravity mediation. Note that in the model of [11], based on non-R discrete symmetry, the µ-problem persists. On the other hand, discrete R-symmetries allow to naturally address this issue [21].
• Neutrinos: The Weinberg operator, LH u LH u , which leads to Majorana neutrino masses, is forbidden in both the superpotential and the Kähler-potential to all orders due to the discrete R symmetry. The neutrinos are Dirac fermions. As we show below, the suppression in the Dirac neutrino masses in our model naturally follow from the family U(1) A symmetry, and its magnitude is determined by the requirement of cancellation of mixed U (1) A -gravitational anomaly.
• Flavor structure: The flavor structure of the model is determined by the Froggatt-Nielsen mechanism [28]. This approach has already been extensively studied in the literature (e.g. [23][24][25][26][27]). Here, the Yukawa couplings are assumed to be of order one, with the specific textures and hierarchy determined by the powers of a small expansion parameter ε. More specifically, the SM fermions carry various positive charges under the U(1) A symmetry, while the MSSM singlet flavon field σ is negatively charged. As the flavon field acquires a VEV σ and breaks U(1) A , the Yukawa terms obtain (Y σ /M ) n couplings, where Y O(1), = ( σ /M ) and M is some suppression scale.
• Proton decay: The potentially dangerous proton decay operators, such as QQQL, do not appear in the superpotential, due to the discrete R-symmetry (see e.g. [31]). The "non-perturbative" contribution from the Kähler potential is insignificant.
• Baryogenesis: The chiral superfield Φ = (φ,φ, F φ ) couples to baryon number-carrying operator U DD and appears as a higher-demensional superpotential term This term will be responsible for generating the baryon asymmetry in the modelá la Affleck-Dine mechanism [32], as follows. First, an asymmetry in φ is generated. Coherent oscillations of the φ-field come to dominate the energy density of the universe after inflation, on-setting matter-domination era. As φ later decays, prior to the BBN, the universe is reheated and the φ-asymmetry is transferred to the baryons through the operator in Eq. (2.4). Our treatment of baryogenesis closely follows that of [11] (see [12][13][14][15][16] for earlier similar works).

U(1) A anomalies
The gauge anomalies associated with U(1) A can be canceled via the four-dimensional Green-Schwarz mechanism as follows. The scalar component of the dilaton chiral superfield S, S| θ=0 = s + ia, contains the dilaton field s as well as the string axion field a. In the presence of S, the gauge boson Lagrangian is given by where W α A , W α a are the superfield strengths of vector supermultiplets V A , V a , and k A , k a are the Kac-Moody algebra levels, associated with the gauge groups U(1) A and G a , respectively. For non-Abelian groups, Kac-Moody levels are positive integers, while for U(1) they need not be integers. The gauge couplings are related to the dilaton VEV as k a S = 1/g 2 a . The above leads to the axion-coupled field strength Lagrangian terms where F A , F a denote the gauge field strength of U(1) A , G a and R is the Riemann curvature tensor, with the last term describing the gravity contribution. Under the U(1) A gauge transformations the dilaton S shifts as where Λ is a parameter chiral superfield and δ GS is a real number. The associated vector superfield transforms as Hence, gauge invariance requires the modified Kähler potential of the dilaton to be of the form Chiral anomalies contribute terms of the form FF A, where A is the anomaly coefficient and F is the field strength of some gauge group. The gravitational anomaly contribution is similar. As the dilaton transformations of Eq. (2.7) induce an axion shift in terms of Eq. (2.6), the anomaly contributions, which enter with the opposite sign, are canceled by an appropriate choice of δ GS . The requirement of anomaly cancellation is fulfilled by (e.g. [33]) The anomaly coefficients are given as follows (see, for example, [27]) where we have already accounted for SU(5) charge multiplet assignment and the standard normalization of 1/2 for the fundamental representation of SU(N ). The hypercharge Y L is not quantized in general and the uncertainty in its normalization renders the anomaly coefficients A Y AA not very informative (see discussion in [34]). The additional possible anomaly coefficient A Y AA due to U(1) A 2 × U(1) Y also suffers from the hypercharge normalization uncertainty and is often neglected. For SU(5) matter charge assignment as well as q Hu = q H d , which is the case of our model, it vanishes. The hypercharge normalization can be fixed by an underlying GUT. The pure [U(1) A ] 3 and gravitational anomalies can obtain potential contributions from the hidden sector SM singlet fields and are thus often neglected. However, the existence of additional light singlet states is potentially of phenomenological interest and gravitational anomaly could still be useful. Thus, we shall only focus on A CCA , A W W A , A Y Y A and A GGA as has been typically done in previous similar studies (e.g. [29]). We note that since SM singlet fields do not carry a non-trivial Dynkin index they only contribute to the gravitational anomaly coefficient. Assuming gauge-coupling unification, which can occur without a simple covering group as in string theory, one has [35] where g s is the string coupling constant. Enforcing the observed gauge-coupling unification in the MSSM requires resulting in This is the same as SU (5) (2.14)

Discrete R-symmetry anomalies
For a local U(1) R-symmetry, the superspace variable θ carries a charge of q θ = R and hence the superpotential carries 2q θ , with R being an integer. For a chiral superfield carrying a charge of q x , the associated fermion carries a charge of (q x − q θ ). This is also the case for Higgsinos, while gauginos carry a charge of q θ . The corresponding anomaly coefficients are then [38,39] For a gauged discrete Z N symmetry anomalies should also cancel. While one cannot construct the corresponding triangle diagram, as there is no gauge boson associated with the discrete symmetry, the anomaly coefficient can still be computed via the Fujikawa method using the path integral measure transformation. The anomaly cancellation condition now becomes where ∆ GS is the dilaton shift constant, and (2.17) In the case of discrete R symmetry, Z R N , the computation is similar. For Z R 12 and charges of Table 1 the corresponding anomaly coefficients are Since in our setup both U(1) A as well as discrete gauge symmetries are present, mixed anomaly coefficients involving both of these can appear. These coefficients are usually neglected, as they depend on the normalization of the U(1) A charges [40] and we will ignore them.

Anomalous U(1) A flavor symmetry
The SM Cabibbo-Kobayashi-Maskawa (CKM) quark matrix can be parametrized through the Cabbibo angle sin θ c 0.2. The Yukawa coupling suppression ( σ /M ) naturally provides the desired small expansion parameter for the flavor structure of the model of the size of the Cabbibo angle, assuming that the model is string-inspired and identifying the scale M with the string scale M s . For simplicity, we take M s = M pl throughout this work. Anomalous flavor U(1) A can also be used to forbid certain operators, through supersymmetric zeroes 3 (e.g. [41]).
The modified Kähler potential invariant under the dilaton transformations (2.7) leads to a radiatively-generated string Fayet-Iliolopuolos (FI) U(1) A term from the gravitational anomaly [42][43][44][45], with a coefficient where A GGA is the gravitational anomaly coefficient as before. Inclusion of a non-zero FI term leads to a D-term contribution to the scalar potential of the form where g A is the U(1) A gauge coupling as before and φ i is the scalar component of the Φ i chiral superfield carrying a U(1) A charge of q φ i . Requirement of SU(3) C × SU(2) W × U(1) Y gauge invariance ensures that the scalar components of MSSM-superfields have vanishing VEVs. In order to maintain SUSY at the U(1) A -breaking scale, the D-term must vanish. This requires that at least one remaining chiral superfield has a negative charge to cancel the ξ FI contribution 4 . This role is played by the flavon field σ, setting its VEV to σ ∝ √ ξ FI . The previously obtained gravitational anomaly coefficient in Eq. (2.14) of (16 + 3p)/4, with p ∈ {0, 1, 2}, thus leads to a desired value of the small expansion parameter where we have taken g 2 s /4π 1/24. After U(1) A breaking the corresponding flavor gauge boson and the gaugino become 3 In this case, the flavon coupling is family-independent and matter fields are allowed to carry negative charges under U(1)A. Since holomorphicity of superpotential only allows for σ and not σ † to appear, any combined operator that carries a total negative U(1)A charge will be forbidden. 4 A U(1)A charge normalization resulting in qσ = −1 is assumed.

massive, with a mass of
Between the scale M pl and M A they can contribute to flavor violating processes, with potentially observable consequences [46].

Textures and massive Dirac neutrinos
With the flavon field, the Yukawa terms appear in the superpotential as For the fermion U(1) A charge assignment of Table 1 as well as Yukawa coupling parametrization, neglecting right handed neutrinos, we employed the results of Model 2 (i.e. α = 1) of [29]. These SU(5)-compatible flavor textures are a variation of textures presented in [47]. Analysis of renormalization group equations (RGEs) as well as arbitrary O(1) coefficients in front of the entries have been also provided by [29]. Using REAP/MPT packages [48], we have confirmed that these textures yield approximately correct mixings. From SU(5) compatibility, Y e = Y T d . However, this leads to a known issue of unacceptable fermion mass relations, with m d /m e = m s /m µ = m b /m τ . A possible solution has been proposed by Georgi and Jarlskog [49], which introduced additional Higgs multiplets that don't couple to all the generations. Adopting their approach, we modify M e as follows For the neutrinos, in contrast to [29], we do not rely on the see-saw mechanism for generation of neutrino masses. Neutrinos in our model are thus Dirac fermions. Without the see-saw mechanism in play, the smallness of the neutrino masses begs for an explanation. We achieve this by requiring that the mixed U(1) A -gravitational anomaly (2.14) is canceled by the GS mechanism without the need of invoking any extra singlet fields, in contrast to [29]. This condition sets the right-handed neutrino charges to the values presented in Table 1 and automatically leads to the desired neutrino mass suppression through the Froggatt-Nielson expansion parameter . The resulting modified Dirac neutrino Yukawa texture is

Fitting fermion masses and mixings
We demonstrate that the above Yukawa textures can lead to reasonable fermion mixings and mass hierarchies by explicitly fitting the O(1) coefficients in front of factors. As the σ flavon field obtains a VEV at the U(1) A breaking scale of ∼ 10 15 GeV, the resulting Yukawa matrices should in principle be evolved with RGEs down to the scale M Z , where the parameters are to be compared to the measured experimental values. Since our goal here is only to show that we can obtain approximately correct flavor structure at the GUT scale, we will not discuss the RGE analysis (see e.g. [29] for a possible implementation).
The Yukawa matrices are diagonalized using singular value decomposition (SVD), which explicitly gives the left-and right-handed unitary rotation matrices V i . The respective rotation matrices for the quark and lepton fields are The CKM and PMNS mixing matrices are then 5 The current measured values of the CKM mixing angles (see PDG [50] where we took the higher of the upper/lower uncertainty and have assumed for illustration purposes the normal hierarchy ordering, which is currently slightly favored by experiments. In addition to the above mixing, if the model is to be compatible with SU (5), at the GUT scale or near it 6 quark and fermion masses should be related. As we have adopted the approach by Georgi and Jarlskog [49] and modified the Y e texture, the resulting GUT-scale masses should be related as (3.13) This will ensure that at low energies the observed mass relations approximately hold. Further, compatibility with SU(5) requires that not only do the textures obey Y d = Y T e , but the respective O(1) coefficients are also related. We note that coefficients presented in the models of [29] are unrelated and thus, strictly speaking, their resulting flavor structures are not compatible with SU(5).
To find the relevant O(1) coefficients for the Yukawa matrices, we have generated 500,000 realizations of Y e , Y d , Y u , Y ν . Each realization includes randomized complex O(1) coefficients in front of 's. From each realization we obtain the CKM and PMNS mixings as well as fermion masses with the help of the MPT package [48].
To select the best fitting realization, we define a χ 2 -like statistic 14) 5 The convention for the PMNS and CKM matrices in the weak currents is L quark where the PMNS (denoted "l") and CKM (denoted "q") components are given as Here, superscript "th" represents the theoretical values obtained from one realization, while "ex" denotes the experimentally observed values as well as their respective uncertainties. The mass χ 2 tot component enforces the Georgi-Jarlskog fermion mass relations.
For p = 0, the best fit coefficients are found to be A better fit accuracy can be obtained with increasing the number of simulated realizations.

Nucleon stability
The QQQL and U U DE baryon and lepton number violating operators of Eq. (2.2) can directly mediate proton decay via p → νK + , which is strongly constrained by the experiment [52]. Since these terms appear only non-perturbatively in the effective superpotential, they are suppressed, with coefficients Flavor structure provides additional strong suppression, which depends on the considered family to which the particles belong and we have specified the minimum (family 3) and the maximum (family 1) possibilities. When coupled together, the R-parity violating dimension-4 terms U DD and LQD can again lead to p → νK + . Thus, there is a constraint on the VEV the Φ field coming from the U DDΦ and LQDΦ terms.
We can neglect these conditions for the Φ field, since the minimum of the potential (5.3) is at φ = 0 if the φ 6 term was absent. Higher dimensional terms in the superpotential polynomial F (x) can in principle modify this, but they are highly suppressed. Similarly, there is a constraint on the VEV of ν from the U DDν and LQDν terms, which appear non-pertubatively and thus carry extra suppression. These are where we have assumed third family contribution, which is most restrictive since it has the smallest flavor suppression and thus leads to a lower (B − L) scale. A possible fine-tuning of the ν VEV can be easily avoided by imposing that (B − L) symmetry is unbroken to some lower energy scale or alternatively that one or both of the contributing / B or / L operators is suppressed or absent. It is also possible that ∆B = 2 neutron-anti-neutron oscillations [53] can appear within the model with different mass-scale suppression than expected. Namely, upon integrating out the Φ-field from the effective superpotential after SUSY/R breaking one obtains the (n − n) mediating operator where in the above we have taken flavor structure corresponding to first generation that describes neutrons, which leads to the resulting operator being strongly suppressed. Similarly, after ν acquires a VEV the U DDν operator can lead to (U DD) 2 , which just highlights the usual statement that (B − L) scale can be related to n − n oscillations. However, since these terms appear non-perturbatively they are highly suppressed. The free n − n lifetime is The constraint of Eq. (4.12) is automatically satisfied in our model when cosmological bounds on reheating are taken into account through Eq. (5.11), even for the most stringent case of p = 2.

Baryon asymmetry and dark matter abundance from moduli decay
In our considerations we remain agnostic about the details of the inflationary phase. During inflation, moduli fields can develop large values and then dominate the energy density of the universe as they coherently oscillate. If they are light and take too long to decay, they can on-set radiation dominated universe at energies below the scale allowed by the BBN, resulting in a "cosmological moduli problem" [2,3]. Making the moduli massive ("stabilizing moduli") [4][5][6][7][8] helps to circumvent this issue. On other hand, massive moduli decaying to gravitinos can give rise to "moduli-induced cosmological gravitino" problem [9,10]. In our model both of these problems are avoided, as the moduli field Φ is massive and is forbidden to decay into gravitinos by the matter parity Z 2 . Here, unlike in the thermal cosmology, reheating radiation originates from decays of the Φ field and not the inflation. Since there are spontaneously broken discrete symmetries, resulting domain walls can be problematic. This issue is avoided, if the inflation scale is not extremely high and the domain walls have time to form and inflate away [55].

Baryogenesis
The U(1) A flavor structure implies that the U DDΦ operator that drives baryogenesis also carries a flavor charge of 5 + 2p (family 1 fields), 2 + 2p (family 2 fields) and 2p (family 3 fields), respectively. The contribution from the third family dominates. For p = 0 (tan β = 20) there is no suppression, while for p = 1 (tan β = 10) and p = 2 (tan β = 5) the operator is suppressed by 4 × 10 −2 and 2 × 10 −3 , respectively. Thus, for p = 0 the resulting estimates for baryon asymmetry and DM abundance would need to be adjusted accordingly. For simplicity, we focus on the p = 0 case below.
Following [11], as in the Affleck-Dine mechanism [32], we define the φ-number as which effectively denotes the difference in the number densities n φ and n φ * of φ and φ * , respectively. For a model of baryogenesis to be successful, all three of the Sakharov's conditions [56] must be fulfilled: CP violation, baryon number violation ( / B) as well as out of equilibrium interactions. Here, q φ -violating terms in the potential generate / B and are also responsible for CP is violation. Since q φ -violation decouples, the condition of being out of equilibrium is also satisfied.
The asymmetry in φ is generated by the suppressed φ 6 Kähler potential term where X is the spurion parametrizing SUSY/R breaking. To preserve q φ during coherent oscillations, the φ-violating terms should be sufficiently suppressed. Taken together with other terms, one obtains the following potential where F (x) is a general polynomial function, κ is a coupling that a priori is assumed to be of O(1) and dots represent higher order contributions. Here, SUSY/R breaking has been taken to be of the gravitino mass size, with m X ∼ m 3/2 .
The dynamics of φ field are described by the equation of motion where H is the Hubble parameter and Γ φ is the φ decay rate. Treating φ as constant during H m φ allows to obtain the initial q φ right after the end of inflation. Rewriting Eq. (5.4) as time evolution of the φ-number giveṡ Then, the initial condition for q φ at t 0 = 1/m φ (i.e. when H ∼ m φ ) is found to be where the initial φ in = (m φ M 3 ) 1/4 was taken as approximately M . Assuming that M is large, i.e. O(M GU T − M pl ), the φ-asymmetry stays nearly constant as φ starts to coherently oscillate. Namely, the amplitude of φ in the right-hand side of Eq. (5.5) becomes small relative to M . Hence, the respective φ-asymmetry is given by where n φ and similarly n φ * are given by ρ φ /m φ , with ρ φ m 2 φ |φ| 2 + |φ| 2 .
After oscillations φ decays (e.g. φ → qqq), reheating the universe and producing superpartners. The φ decay rate is given by where the numerical factor is ξ = 27/(256π 3 ) = 3 × 10 −3 [11], assuming all φ couplings to quarks as one. The reheat temperature T r is given by [13] where g * counts the effective number of relativistic degrees of freedom and the whole numerical pre-factor is O(1) around the MeV scale, which corresponds to BBN. At H Γ Φ , T r is found to be T r 20 GeV m φ 5 × 10 4 TeV Assuming that φ dominates the energy density before its decay, T r is related to q φ via radiation density Since the baryon number density n b is set by q φ via n b = q φ , the resulting baryon asymmetry is found to be (see Model A of [11]) (5.13)

Dark matter
With matter parity conserved, the lightest supersymmetric partner (LSP) is stable and thus constitutes the DM candidate, with a possible natural choice being the neutralino. Depending on the details of supersymmetry breaking, the sparticle spectrum will be altered and so will the identity of the neutralino that is the LSP. Since we require a heavy gravitino of m 3/2 50 TeV for the scenario considered, loop corrections to soft masses (gauginos, etc.) are relevant. This is a natural setting for anomaly-mediated supersymmetry breaking [4,5] with Wino DM [16]. If the modulus and anomaly mediation contributions are competitive, Higgsino with mass below O(TeV) becomes the DM candidate [58]. This is easily achieved in the context of "mirage mediation" [59,60], which could be further extended to "deflected mirage mediation" [61], where the gaugino masses are given by (1 + 0.1α m )g 2 2 : (1 + 0.03α m )g 2 1 . (5.14) Here, α m = m 3/2 /M 0 log(M pl /m 3/2 ) parametrizes the relative strength of the anomaly and modulus-mediated contributions, with M 0 denoting the modulus-mediated contribution at the GUT scale, and g 1,2,3 are the gauge coupling constants. In the limit of α m → 0, anomalymediated contribution vanishes and Bino becomes the DM candidate. where m nuc is the nucleon mass. The Ω DM /Ω b factor is required to be of O(1) if one is to address the observed baryon-DM abundance coincidence. We first discuss Wino LSP, arising in the context of anomaly-mediated SUSY breaking. Here, the gaugino masses are proportional to It has been suggested (e.g. [63]) that indirect detection already strongly constrains the nonthermal Wino. These results, however, strongly depend on the assumptions regarding the uncertain astrophysical J-factor associated with the DM halo shape in dwarf spheroidal galaxies [64]. Assuming a more cored DM profile allows to alter the constraints by more than an order of magnitude in the m χ ∼ O(1) TeV parameter range. Indirect detection does provide a relatively robust constraint for m χ < 300, which in turn suggests a higher reheating temperature as well as a heavier modulus in the model. Future gamma-ray observations will allow to further probe the m χ ∼ O(1) TeV region. While indirect detection constraints are subject to large astrophysical uncertainties, disappearing track searches at colliders have proven to be a powerful probe for wino DM. Current Large Hadron Collider (LHC) searches already constrain wino with m χ < 460 GeV and upcoming high luminosity HL-LHC run studies will be sensitive to wino masses of up to 900 GeV [65]. For non-thermal Higgsino, the mass region of m χ is already also restricted [66]. However, there exists an open parameter space region with T r ∼ 20 − 50 GeV and m χ ∼ 600 − 800 GeV, which will be probed by the future Cherenkov Telescope Array (CTA) and 100 TeV collider experiments [67]. For representative values of T r = 35 GeV and m χ = 700 GeV, the reheating (Eq. (5.11), Eq. (5.10), Eq. (5.16)), baryogenesis (Eq. (5.13)), DM abundance (Eq. (5.18)) and DM-baryon coincidence (Eq. (5.19)) constraints are automatically satisfied for m φ 5 × 10 4 TeV and m 3/2 10 3 TeV, assuming M = M pl as before.
We note that the mentioned DM constraints can be further alleviated by reducing the contribution of neutralinos to the overall DM abundance. This can be achieved, for example, within the frameworks where DM is a mix of axion(axino) and Higgsino [68,69]. As has been shown in [70], this is realizable within a setting based on discrete gauge symmetries and can also lead to resolution of the strong CP problem via the Peccei-Quinn axion [70]. A full implementation of above into our model, however, is beyond the scope of the present work. Other DM candidates, such as primordial black holes, have also been suggested in SUSY-based setting with scalar fields whose coherent-oscillations lead to an onset of a matterdominated phase in the early universe [71].

Conclusion
The origin of the cosmological matter-antimatter asymmetry, the gauge hierarchy problem, the identity of the dark matter as well as the pattern of the fermion mass hierarchy and mixing angles remain among the greatest challenges for the Standard Model and Cosmology. While minimal supersymmetric extension of the Standard Model remains the most promising solution for the gauge hierarchy problem, it has difficulty confronting the longevity of the nucleon lifetime and the smallness of the µ term. Further, the standard cosmological scenarios with thermal production of baryons and dark matter often face severe inconsistencies.
In this work, we consider non-thermal moduli cosmology, which provides a consistent framework for generation of radiation, baryons, and dark matter, while maintaining successful BBN and avoiding the cosmological moduli problem. We construct an SU(5)-compatible model based on discrete Z R 12 R-symmetry combined with an anomalous U(1) A flavor symmetry, in which the aforementioned problems are resolved simultaneously. Specifically, the Z R 12 symmetry plays an important role in solving the µ problem, while the U(1) A symme-try addresses the flavor puzzle. The interplay between Z R 12 and U(1) A symmetries leads to the absence or high suppression of proton decay operators as well as the prediction of Dirac neutrinos. The requirement of mixed U(1) A -gravitational anomaly cancellation automatically give rise to the correct size of the Dirac neutrino masses. The decay of the moduli φ field non-thermally produces baryons as well as dark matter, identified as the neutralino LSP. Consequently, the coincidence between the baryon number and dark matter abundance, Ω B ∼ Ω DM , can be explained.