We have numerically explored the scalar field condensate baryogenesis model for numerous sets of model’s parameters, within their natural range of values. We have investigated the evolution of the baryon charge carrying field, the evolution of the baryon charge contained in the scalar field condensate, and the final value of the generated baryon charge on the model’s parameters: the gauge coupling constant α, the Hubble constant at the inflationary stage HI, the mass m, and the self-coupling constants λi.
1. Introduction
There exists baryon asymmetry β≠0 in the neighborhood of our Galaxy, within 20 Mpc. β is usually parameterized as β=(nb-nb-)/nγ, where nb is the number density of baryons, nb- is the number density of antibaryons, and nγ is the number density of photons. Contemporary observational knowledge on the baryon density of the Universe is based mainly on the following sets of precise observational data: data based on Big Bang Nucleosynthesis (BBN), that is, the determination of the baryon density from the requirement of consistency between theoretically predicted abundance and observationally measured abundance of the primordially produced light elements [1]; measurements of Deuterium towards low metallicity distant quasars compared with BBN predicted D [2]; and CMB anisotropy measurements (see WMAP [3] and Planck [4]), allowing precise determination of the main Universe characteristics, including the baryon density.
Namely, the consistency between theoretically obtained abundance and observationally measured abundance of the light elements produced in BBN [1] requires that the baryon-to-photon density is in the range:(1)5.1×10-10≤nbnγBBN≤6.5×10-10at 95% CL.
The information from measurements of Deuterium towards low metallicity quasars combined with BBN data [2] points to(2)nbnγD=6±0.3×10-10at 95% CL.
The most precise determination is provided by the measurements of the CMB anisotropy (z~1000). Recent results by WMAP9 [3] point to(3)nbnγWMAP=6.19±0.14×10-10at 68% CL.
The up-to-date data from Planck project [4] point to(4)nbnγPlanck=6.108±0.038×10-10at 68% CL.
Observational constraints exist on the presence of the antimatter in our local vicinity (within tens of Mpc), mainly based on Cosmic Ray data [5–9] and Gamma Ray data [10–14]. No antimatter in astronomically considerable amounts has been observed/detected (still, domains of antimatter are not absolutely ruled out [15–17]). Hence, the baryon asymmetry in our neighborhood is β~nb/nγ. Observational evidence for matter-antimatter asymmetry in the Universe has been recently reviewed in [17, 18].
In case this locally observed asymmetry is a global characteristic of the Universe, that is, baryon asymmetry of the Universe (BAU), it may be due to the generation of a baryon excess at some early stage of the Universe. This excess eventually diluted during the Universe further evolution determined the BAU observed today.
Sakharov [19] defined the conditions for the generation of predominance of matter over antimatter from initially symmetric state of the early Universe. Namely, these are nonconservation of baryons, C and CP-violation, and deviation from thermal equilibrium. None of these conditions is obligatory. Different baryogenesis scenarios, where some of the Sakharov’s requirements are not fulfilled, have been discussed in [20].
The exact (nature chosen) baryogenesis mechanism is not known yet. It is known that baryon asymmetry may not be postulated as an initial condition, in case of inflationary early stage of the Universe evolution, and it should have been generated in the period after inflation and before BBN epoch [21].
Numerous baryogenesis scenarios exist today, whose aim is to explain the observed baryon asymmetry, its sign, and its value. For a review see [22–25]. In [24] different mechanisms for generation both of the baryon asymmetry and the dark matter of the Universe and their dependence on the reheating temperature have been discussed. The most studied among the baryogenesis scenarios are Grand Unified Theories (GUT) baryogenesis [19], electroweak (EW) baryogenesis [26–30], baryogenesis-through-leptogenesis (often called leptogenesis) [31–33], Affleck-Dine (AD) baryogenesis [34], and so forth.
GUT baryogenesis is the earliest baryogenesis scenario, proceeding at GUT unification scale MGUT. However, most inflationary models predict reheating temperature below this scale. Besides, successful unification requires supersymmetry. SUSY implies the existence of gravitinos, which are too numerously produced unless the reheating temperature is well below MGUT [35, 36].
EW baryogenesis is theoretically attractive because it relies only upon weak scale physics and is experimentally testable scenario. For a review see [28, 37–39]. However, the simplest and most appealing version of this scenario cannot generate within the Standard Model the observed value of the baryon asymmetry, because of the insufficient CP-violation [40–42] induced by CKM phase and the requirement of first-order electroweak transition, possible only for Higgs boson mass considerably smaller than the detected one by ATLAS and CMS collaborations.
EW baryogenesis in Minimal Supersymmetric Standard Model was considered [43]. Now MSSM window for EW baryogenesis is substantially narrowed by the experimental data from the Large Hadron Collider (LHC) [44–46] and recent electron dipole measurements. Constraints on EW baryogenesis in case of a minimal extension of the Standard Model from current data from LHC have been discussed as well [47]. The viable parameter space is considerably reduced. Baryogenesis in next-to-minimal SM is being discussed now [48].
Baryogenesis-through-leptogenesis is a plausible possibility. Baryon asymmetry in this scenario is created before the electroweak phase transition, which then gets converted to the baryon asymmetry in the presence of (B+L) violating anomalous processes. For a review see [23, 49–51]. It has become especially attractive after the discovery of nonzero neutrino masses in neutrino oscillations experiments. Baryogenesis through leptogenesis mechanisms in different extensions of the SM is studied.
Possibilities for falsifying concrete realizations of high scale leptogenesis from recent LHC data have been proposed [52, 53].
Neutrino Minimal Standard Model (NMSM) can potentially account simultaneously for baryon and dark matter generation and neutrino oscillations [54]. NMSM is testable at colliders and in astrophysical observations. For a recent review and constraints from collider experiments, astrophysics, and cosmology see [55].
AD baryogenesis scenario [22, 34] is one of the most promising today baryogenesis scenarios, compatible with inflation. Nice reviews of AD baryogenesis contemporary status can be found in [22, 56]. AD baryogenesis has numerous attractive features. A short list of these is as follows: (i) it is extremely efficient: it can produce equal or much bigger baryon asymmetry than the observed one; (ii) it can be realized at lower energy, that is, relatively late in the Universe evolution. That is, it is consistent with the low energy scales after inflation; (iii) AD condensate can be generated generically in different cosmological models; (iv) it can explain simultaneously the generation of the baryon and the dark matter in the Universe and explain their surprisingly close values; (v) AD model, due to its high efficiency, can be successful even in case of significant production of entropy at late times, predicted by some particle physics models.
AD scenario is based on SUSY. In supersymmetric models scalar superpartner of baryons and leptons φ exists. The potential U(φ) of such scalar field may have flat directions, along which the field can have a nonzero vacuum expectation value due to quantum fluctuations during inflation. After inflation φ evolves down to the equilibrium point φ=0 and if the potential U(φ) is not symmetric with respect to the phase rotation it acquires nonvanishing and typically large baryon charge. Subsequent B-conserving decay of φ into quarks and leptons transforms baryon asymmetry into the quarks sector. In contrast to other scenarios of baryogenesis, in which the generated asymmetry usually is insufficient, the original Affleck-Dine scenario leads to higher value of β and additional mechanisms are needed to dilute it down to the observed value.
AD mechanism was reexamined in [57]. It was realized that finite energy density of the early Universe breaks SUSY and induces soft parameters in the soft potential along flat directions, which are of the order of the Hubble parameter. Then, contrary to the original AD mechanism the observed value of the baryon asymmetry may be generated, without the requirement of subsequent entropy release. Different issues on AD baryogenesis were presented in [58–60]. AD baryogenesis mechanism was used in numerous SM extensions and different inflationary scenario. A short list of several of the more recent studies includes: AD in effective supergravity [61], AD in anomaly mediated SUSY breaking models [62], AD in SUSY with R-parity violation [63], AD in D-term inflation [64], and so forth. Most of AD baryogenesis models can naturally explain the origin of the dark matter in the Universe. Constraints on the subclass of AD models were obtained from current CMB data, based on the backreaction of the flat direction on the inflationary potential [65].
Here we discuss the scalar field condensate baryogenesis model (SFC baryogenesis), which is among the preferred today baryogenesis scenarios, compatible with inflation. It is based on the Affleck-Dine scenario.
SFC baryogenesis model was first discussed and studied analytically in [66, 67]. There it was shown that the account of particle creation by the time varying scalar field during postinflationary period will lead to strong reduction of the produced baryon excess in the Affleck-Dine scenario. Namely, it was proven that fast oscillations of φ result in particle creation due to the coupling of the scalar field to fermions gφf1f2, where g2/4π=α. For λi3/4>g the rate of particle creation Γ exceeds the ordinary decay rate of φ at the stage of baryon nonconservation and, therefore, its amplitude is damped. Hence, the baryon charge, contained in the condensate, is reduced due to particle creation at this stage with considerable baryon violation (when Γ = constant, the baryon charge of the condensate is reduced exponentially and does not survive till φ decays to quarks and leptons. However, when Γ is a decreasing function of time, the damping process may be slow enough for the baryon charge in the condensate to survive until its decay; i.e., this case is more promising).
The importance of a precise numerical account for the particle creation processes was explored further in [68, 69]. Different possibilities of SFC baryogenesis models were discussed. The possibility to generate simultaneously, within inhomogeneous SCF baryogenesis model, the observed baryon asymmetry and the observed large scale structure quasiperiodicity of the baryonic matter was studied in [68, 70, 71]. On the basis of inhomogeneous SCF baryogenesis elegant mechanisms were proposed for achieving sufficient separation between domains of matter and antimatter (to inhibit the contact and evade annihilation of matter and antimatter regions with big density) that allow the production of considerable antimatter domains with different size in the Universe and their observational signatures were analyzed [15, 17, 71–74].
In series of papers [69, 75, 76], we explored numerically SFC baryogenesis model. Here we present the results of our extended numerical analysis of the evolution of the baryon excess in SFC baryogenesis model and its dependence on the model parameters.
In the next section we briefly describe the SCF baryogenesis model and the numerical approach we have used. The last section presents the results; that is, we present the value of the produced baryon density for numerous sets of model’s parameters.
2. SFC Baryogenesis Model2.1. Description
The essential ingredient of the model is a baryon charged complex scalar field φ, present together with the inflaton. A condensate φ≠0 with a nonzero baryon charge is formed during the inflationary period as a result of the rise of quantum fluctuations of the φ field [77–80]: φ2=H3t/4π2 until the limiting value φ2~H2/λ in case that λφ4 terms dominate in the potential energy of φ.
The baryon charge of the field is not conserved at large field amplitude due to the presence of the B nonconserving self-interaction terms in its potential.
We choose the form of the potential as follows:(5)Uφ=m2φ2+λ12φ4+λ24φ4+φ∗4+λ34φ2φ2+φ∗2.
The mass parameters of the potential are assumed to be small in comparison with the Hubble constant during inflation m≪HI. In supersymmetric theories the self-coupling constants λi are of the order of the gauge coupling constant α. A natural range of m is 102–104 GeV.
We examine the case when after inflation there exist two scalar fields, the inflaton ψ and the scalar field φ, and the inflaton density dominates prior to the decay of φ: ρψ>ρφ. Hence, at the end of inflation the Hubble parameter is H=2/(3t).
In the expanding Universe, spatially homogeneous field φ satisfies the equation of motion:(6)φ¨+3Hφ˙+14Γφ˙+Uφ′=0,where a(t) is the scale factor and H=a˙/a, Γ accounts for the particle creation processes.
The initial values for the field variables are derived from the natural assumption that the energy density of φ at the inflationary stage is of the order HI4; then(7)φomax~HIλ-1/4,φ˙o=HI2.
After inflation φ oscillates around its equilibrium point and its amplitude decreases due to the Universe expansion and the particle creation by the oscillating scalar field. In case Γ is a decreasing function of time the damping process may be slow enough for the baryon charge contained in φ to survive until the B-conservation epoch [67].
At low φ baryon violation (BV) becomes negligible. At the B conserving stage the baryon charge contained in the field is transferred to that of quarks during the decay of the field φ→qq-lγ at tb. As a result, in case φ has not reached the equilibrium point at tb, the baryogenesis makes a snapshot of φ(tb) and a baryon asymmetric plasma appears. This asymmetry, eventually further diluted during the following evolution of the Universe, gives the present observed baryon asymmetry of the Universe.
2.2. Evolution of the Baryon Charge Carrying Field
We have solved the system of ordinary differential equations, corresponding to the equation of motion for the real and imaginary components of φ=x+iy:(8)x¨+3Hx˙+14Γxx˙+λ+λ3x3+λ′xy2=0,y¨+3Hy˙+14Γyy˙+λ-λ3y3+λ′yx2=0,where λ=λ1+λ2, λ′=λ1-3λ2.
It is convenient to make the substitutions x=HI(ti/t)2/3u(η), y=HI(ti/t)2/3v(η), where η=2(t/ti)1/3. Then the functions u(η) and v(η) satisfy the equations:(9)u′′+0.75αΩuu′-2uη-1+uλ+λ3u2+λ′v2-2η-2+m2Hη4=0,v′′+0.75αΩvv′-2vη-1+vλ-λ3v2+λ′u2-2η-2+m2Hη4=0.
The baryon charge in the comoving volume V=Vi(t/ti)2 is given by(10)B=NB·V=2u′v-v′u.
We have solved numerically the system of ordinary differential equations (9), corresponding to the equation of motion for the real and imaginary part of φ and B contained in it, using Runge-Kutta fourth-order method and Fortran 77. The Runge-Kutta fourth-order routine from [81] is used.
We studied numerically the evolution of φ(η) and B(η) in the period after inflation until the BC epoch. The typical range of energies discussed was 1015–100 GeV. Therefore, serious computational resources were used. A single calculation took between several hours and three weeks, depending on the concrete parameters.
We analyzed φ and B evolution for natural ranges of values of the model’s parameters: λ=10-2–5×10-2, α=10-3–5×10-2, HI=107–1016 GeV, and m=100–1000 GeV. The numerical analysis was provided for around seventy sets of parameters.
We have accounted numerically for the particle creation processes by varying φ, which allowed us to describe more precisely the evolution of B and determine its final value which was transferred to quarks (antiquarks) at tb epoch and defined the baryon asymmetry. In this work we calculated Γ numerically in contrast to our previous papers, where for the rate of particle creation Γ the analytical estimation was used Γ=αΩ, where Ω~λ11/4φ. In the program Ωu and Ωv were calculated at each step in separate routine procedures.
The results of our numerical study are presented in the next section.
3. The Generated Baryon Charge for Different Parameters Values: Numerical Results
We have calculated B for different sets of values of model’s parameters: gauge coupling constant α, Hubble constant during inflation Hi, mass of the condensate m, and self-coupling constants λi. We have not made calculations for all the possible values of the parameters, because each single point requires days or weeks of CP time. Our main aim was to find the dependence of the final B on the parameters and choose the more promising ranges of the parameters for successful baryogenesis, rather than providing full systematic numerical study. Therefore, some entries in Tables 1, and 3 are missing.
The baryon charge B contained in the SFC at the time of its decay for different HI and m and fixed set of λ1=10-2, α=10-2, λ2=λ3=10-3, φo=HIλ-1/4, and φ˙o=HI2. The particle creation processes are accounted for numerically.
m, GeV
H, GeV
1.00E+14
1.00E+12
1.00E+11
1.00E+10
1.00E+09
1.00E+08
1.00E+07
350
-2.75E-07
2.10E-06
1.18E-04
1.61E-03
2.89E-02
-4.08E-01
-2.31E-01
500
7.44E-05
-7.72E-04
-6.72E-04
-1.00E-02
2.69E-01
1.58E-01
800
-1.07E-04
-1.40E-03
2.27E-03
2.86E-02
1.11E-01
9.04E-01
3.1. Dependence on Hubble Constant during Inflation HI
We have followed the evolution B(η) varying HI for fixed values of the other parameters. The results of our preliminary analysis of this dependence have been first discussed in [69, 75].
In this work we have performed an extended analysis of this dependence, studying wider range of models parameters.
In the rows in Tables 1, 2, and 3 we present the results of the generated baryon charge in the SFC baryogenesis model for different values of the Hubble constant HI and for different fixed values of the other parameters of the model. The first table presents B(HI) for λ1=10-2, α=10-2, and λ2=λ3=10-3, the second one presents it for λ1=5×10-2, α=10-2, and λ2=λ3=10-3, and the third one presents it for λ1=5×10-2, λ2=λ3=10-3, m=350 GeV, and two different values of α, namely, 10-2 and 5·10-2.
The baryon charge B contained in the SFC at the time of its decay for different HI and m and fixed set of λ1=5×10-2, α=10-2, λ2=λ3=10-3, φo=HIλ-1/4, and φ˙o=HI2.
m, GeV
H, GeV
1.00E+12
1.00E+11
1.00E+10
1.00E+09
1.00E+08
1.00E+07
350
1.00E-05
1.00E-04
-5.00E-05
-1.00E-04
7.00E-03
-4.00E-03
500
2.57E-06
1.00E-06
-2.00E-04
-1.34E-02
7.87E-02
5.00E-01
The baryon charge B contained in the SFC at the time of its decay for different α and HI and λ1=5×10-2, λ2=λ3=10-3, m=350 GeV, φo=HIλ-1/4, and φ˙o=HI2.
α
H, GeV
1.00E+12
1.00E+11
1.00E+10
1.00E+09
1.00E+08
1.00E+07
1.00E-02
1.00E-05
1.00E-04
-5.00E-05
-1.00E-04
7.00E-03
-4.00E-03
5.00E-02
8.01E-07
5.50E-05
2.40E-03
8.10E-03
6.90E-02
Figure 1 presents the dependence of the evolution of the baryon charge on HI at fixed values of the other parameters.
The evolution of the baryon charge B(η) for λ1=10-2, λ2=λ3=10-3, α=10-2, m=500 GeV, φo=HIλ-1/4, and φ˙o=HI2. (a) For HI=109 GeV, (b) for HI=1010 GeV, (c) for HI=1011 GeV, and (d) for HI=1012 GeV.
Our detailed analysis for numerous different parameters of the SCF model shows that B evolution becomes longer and the final B value decreases with the increase of HI. This result is in agreement with previous numerical and analytical studies. It is an expected result because particle creation, which reduces β, is proportional to φ, Γ~Ω~φ, and the initial value of φ is proportional to HI. Thus, the bigger HI is, the more efficient is the decrease of β due to particle creation.
3.2. Dependence on Gauge Coupling Constant α
Using the numerical account for Γ we have calculated B(η) for α varying in the range 10-3–10-2 and fixed other parameters; see also [69]. The dependence of B on α is very strong, as can be expected, knowing that particle creation processes play essential role for the evolution of the field and the baryon charge, contained in it, and keeping in mind that the analytical estimation is Γ=αΩ.
With increasing α, B evolution becomes shorter and the final B decreases. An illustration of this dependence B(α) is given in Figure 2.
The evolution of the baryon charge B(η) for λ1=5×10-2, λ2=λ3=10-3, H=1010 GeV, m=350 GeV, φo=HIλ-1/4, and φ˙o=HI2. (a) For α=10-3, (b) for α=10-2, and (c) for α=5×10-2.
In the columns of Table 3 we present the results for the baryon charge contained in the SFC at the time of its decay for different values of α and for λ1=5×10-2, λ2=λ3=10-3, and m=350 GeV.
3.3. Dependence on the Mass m of the Condensate
The dependence of the final baryon charge on m for fixed λ1, λ2, λ3, α, and HI has been first discussed in [69, 75]. Here we present the results of a more detailed study. In Table 2 the dependence of generated baryon asymmetry values on the mass of the field is given for different fixed sets of the other parameters of the model.
It has been found that in general B decreases with the increase of the value of m. This behavior is more clearly and more strongly expressed for big values of HI and then corresponds to the expected one from analytical estimations (namely, as far as m defines the onset of BC epoch: ts~1/αm, and hence for lower values of m, B evolution is longer and the final B value is smaller). The dependence is illustrated in Figure 3.
The evolution of the baryon charge B(η) for λ1=5×10-2, λ2=λ3=10-2, α=5×10-2, H=1011 GeV, φo=HIλ-1/4, and φ˙o=HI2. (a) Curve is for m=100 GeV, (b) for m=200 GeV, and (c) for m=350 GeV.
For smaller values of HI the dependence is weaker and not so straightforward and clear. In Table 2 we present the generated baryon asymmetry values for different mass of the field and fixed values of the other parameters of the model.
In the columns of Tables 1 and 2 we present the results for the baryon charge contained in the SFC at the time of its decay for different values of the mass of the field m and for different fixed values of the other parameters of the model.
3.4. Dependence on the Self-Coupling Constants
The dependence of the baryon charge, at the B-conservation epoch, on the value of the coupling constants λi was first discussed in [76]. Our extended analysis confirms that the final B value decreases when increasing λ1 and B evolution becomes shorter. The effect provides a difference in the final B value of an order of magnitude.
The results of the extended analysis showed that the final value of B may be much sensitive to λ2 and λ3 compared to that previously estimated. Namely, we have found that the final values of B may differ up to 3 orders of magnitude even for small changes of these parameters. For example, for m=350 GeV, HI=1010, α=10-3, and λ1=5×10-2 the final B value for λ2=λ3=5×10-4 is B=3.88×10-5, while for λ2=λ3=10-3 it is B=2.36×10-2.
Thus, we have found that the dependence of baryon generation on the self-coupling constants is important for determination of the parameters range for the successful baryogenesis model.
4. Estimation of the Generated Baryon Asymmetry
In order to estimate the baryon asymmetry on the basis of the obtained results for the produced baryon density it is necessary to know the temperature of the relativistic plasma after the decay of φ and the decay of the inflaton.
In case the inflaton energy density dominates until the decay of φ, that is, prior to reheating, ρψ>ρφ, the entropy is mainly defined by the relativistic particles from inflaton decay. Thus, the temperature after the decay of ψ at tψ will be approximately (11)TR~ρψ1/4=ρψ01/4η0ηψ3/2.
Then the baryon asymmetry will be given by(12)β~NBTR3~BTRHI,where TR is the reheating temperature after the decay of the inflaton.
Hence, from these estimations it is seen that the lower the reheating temperature after inflaton decay is the lower the produced baryon asymmetry will be. Also, the later the inflaton decays the smaller the produced β will be. Knowing that the reheating temperature should be sufficiently low to avoid gravitino problem, that is, in our model it should be several orders or more lower than the value of HI, and having the results for B, it is easy to obtain the value of the observed baryon asymmetry for different sets of parameters in this model.
Of course HI, the decay time of ψ, and the value of the reheating temperature may be different in different inflationary scenarios. We would like only to note here that the results of the numerical analysis of the SCF baryogenesis model are encouraging. The analysis points that this model provides an opportunity to produce baryon asymmetry β, consistent with its observed value for natural values of the model’s parameters. Therefore, this model deserves further considerations.
5. Conclusions
It was found that the analytical estimations of the baryon charge evolution and its final value in SFC baryogenesis model may considerably differ from the exact numerically calculated ones. Therefore, in this work we have numerically explored the SFC baryogenesis model for numerous sets of model’s parameters.
We have investigated the dependence of the evolution of the field and the evolution of the baryon charge contained in it, as well as the final value of the baryon charge contained in it on the model’s parameters: the gauge coupling constant α, the Hubble parameter at the inflationary stage HI, the mass m, and the self-coupling constants λi. Qualitative dependence of the final B on these parameters has been found. Namely, it was shown that the produced baryon excess is a strongly decreasing function of α, and it also is a decreasing function of HI. The dependence on m is not so straightforward. For small m values B decreases with m increase; however, for larger m the dependence is more complicated.
The analysis may be used to indicate the values of the model’s parameters for which baryon asymmetry β, consistent with its observed value, may be produced in a given inflationary scenario.
The results of this analysis may be used for constructing realistic SCF baryogenesis models. Moreover, assuming SCM baryogenesis and assuming a concrete inflationary scenario, from the observed value of the baryon asymmetry it is possible to put cosmological constraints on the model’s parameters, provided by physics theories, that is, constrain physics beyond Standard Model.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
The authors would like to thank M. V. Chizhov for useful conversations and for the technical support. D. Kirilova expresses her gratitude to O. V. Teryaev for the overall help during her visiting position at BLTP, JINR. The authors are grateful to the unknown referee for the useful suggestions and criticism.
OliveK. A.GolwalaS.VogelP.Review of particle physics201438909000110.1088/1674-1137/38/9/090001PettiniM.CooceR.A new, precise measurement of the primordial abundance of Deuterium201242542477248610.1111/j.1365-2966.2012.21665.xBennettC. L.LarsonD.WeilandJ. L.Nine-year Wilkinson Microwave Anisotropy Probe (WMAP) observations: final maps and results201220822010.1088/0067-0049/208/2/20AdeP. A. R.AghanimN.ArnaudM.Planck 2015 results. XIII. Cosmological parametershttp://arxiv.org/abs/1502.01589AbeK.FukeH.HainoS.Search for antihelium with the BESS-polar spectrometer201210813130110.1103/physrevlett.108.131301AdrianiO.BarbarinoG. C.BazilevskayaG. A.PAMELA results on the cosmic-ray antiproton flux from 60 MeV to 180 GeV in kinetic energy201010512512110110.1103/PhysRevLett.105.121101MayorovA. G.GalperA. M.AdrianiO.Upper limit on the antihelium flux in primary cosmic rays2011931162863110.1134/s0021364011110087AlcarazJ.AlvisiD.AlpatB.Search for antihelium in cosmic rays1999461438739610.1016/S0370-2693(99)00874-6KapplR.WinklerM. W.The cosmic ray antiproton background for AMS-0220142014905110.1088/1475-7516/2014/09/051SteigmanG.Observational tests of antimatter cosmologies19761433937210.1146/annurev.aa.14.090176.002011SteigmanG.When clusters collide: constraints on antimatter on the largest scales2008200810, article 00110.1088/1475-7516/2008/10/001SteckerF. W.On the nature of the baryon asymmetry1985252253610.1016/0550-3213(85)90423-72-s2.0-0002991852CohenA. G.De RújulaA.GlashowS. L.A matter-antimatter universe?1998495253954910.1086/3053282-s2.0-22044438743von BallmoosP.Antimatter in the Universe: constraints from gamma-ray astronomy20142281–39110010.1007/s10751-014-1024-92-s2.0-84912025502DolgovA.SilkJ.Baryon isocurvature fluctuations at small scales and baryonic dark matter199347104244425510.1103/physrevd.47.4244KhlopovM. Y.RubinS. G.SakharovA. S.Possible origin of antimatter regions in the baryon dominated universe20006208350510.1103/physrevd.62.083505DolgovA.Antimatter in the universe and laboratory20159503007CanettiL.DrewesM.ShaposhnikovM.Matter and antimatter in the universe20121409501210.1088/1367-2630/14/9/0950122-s2.0-84867018997SakharovA.Violation of CP invariance, С asymmetry, and baryon asymmetry of the universe1967524DolgovA. D.Non-GUT baryogenesis1992222630938610.1016/0370-1573(92)90107-bDolgovA. D.ZeljdovichY.SazhinM.1988DineM.KusenkoA.Origin of the matter-antimatter asymmetry200476110.1103/revmodphys.76.12-s2.0-1542309866BuchmullerW.PecceiR. D.YanagidaT.Leptogenesis as the origin of matter20055531135510.1146/annurev.nucl.55.090704.151558BuchmüllerW.Baryogenesis, dark matter and the maximal temperature of the early universe201243122153218510.5506/aphyspolb.43.2153DolgovA. D.ZeljdovichYa. B.Cosmology and elementary particles198153114110.1103/revmodphys.53.1KuzminV. A.RubakovV. A.ShaposhnikovM. E.Anomalous electronweak baryon number non-conservation and GUT mechanism for baryogenesis19871911-217117310.1016/0370-2693(87)91340-2KuzminV. A.RubakovV. A.ShaposhnikovM. E.On anomalous electroweak baryon-number non-conservation in the early universe19851551-2364210.1016/0370-2693(85)91028-72-s2.0-0003017602ShaposhnikovM. E.Possible appearance of the baryon asymmetry of the universe in an electroweak theory1986444651CohenA. G.KaplanD. B.NelsonA. E.Baryogenesis at the weak phase transition1991349372774210.1016/0550-3213(91)90395-eCohenA.KaplanD.NelsonA.Progress in electroweak baryogenesis199343277010.1146/annurev.ns.43.120193.000331FukugitaM.YanagidaT.Barygenesis without grand unification19861741454710.1016/0370-2693(86)91126-3AsakaT.BlanchetS.ShaposhnikovM.The νmSM, dark matter and neutrino masses2005631415115610.1016/j.physletb.2005.09.0702-s2.0-27944475229AsakaT.BlanchetS.ShaposhnikovM.The νMSM, dark matter and baryon asymmetry of the universe20056201-2172610.1016/j.physletb.2005.06.020AffleckI.DineM.A new mechanism for baryogenesis1985249236138010.1016/0550-3213(85)90021-5RubakovV. A.SazhinM. V.VerjaskinA. V.Graviton creation in the inflationary universe and the grand unification scale1982115318919210.1016/0370-2693(82)90641-4WeinbergS.Cosmological constraints on the scale of supersymmetry breaking198248191303130610.1103/physrevlett.48.1303ClineJ.Status of electroweak phase transition and baryogenesis2000551-2334210.1007/s12043-000-0081-6RubakovV. A.ShaposhnikovM. E.Electroweak baryon number nonconservation in the early universe and in high-energy collisions1996166493537MorrisseyD. E.Ramsey-MusolfM. J.Electroweak baryogenesis20121412500310.1088/1367-2630/14/12/1250032-s2.0-84871910908GavelaM. B.HernándezP.OrloffJ.PèneO.Standard model CP-violation and baryon asymmetry19949979580910.1142/S0217732394000629GavelaM.HernandezP.OrloffJ.PéneO.QuimbayC.Standard model CP-violation and baryon asymmetry (II). Finite temperature1994430238242610.1016/0550-3213(94)00410-2HuetP.SatherE.Electroweak baryogenesis and standard model CP violation199551237939410.1103/physrevd.51.3792-s2.0-33744543821HuetP.NelsonA. E.Electroweak baryogenesis in supersymmetric models19965384578459710.1103/physrevd.53.45782-s2.0-0000855405CurtinD.JaiswalP.MeadeP.Excluding electroweak baryogenesis in the MSSM201220128, article 00510.1007/JHEP08(2012)005CarenaM.NardiniG.QuirosM.WagnerC. E. W.MSSM electroweak baryogenesis and LHC data201313021CirglianoV.LiY.ProfumoS.Ramsey-MusolfM. J.MSSM baryogenesis and electric dipole moments: an update on the phenomenology201020101, article 00212310.1007/jhep01(2010)002DamgaardP. H.O'ConnellD.PetersenT. C.TranbergA.Constraints on new physics from baryogenesis and large hadron collider data20131112222180410.1103/physrevlett.111.221804BalázsC.MazumdarA.PukartasE.WhiteG.Baryogenesis, dark matter and inflation in the next-to-minimal supersymmetric standard model20142014107310.1007/jhep01(2014)073BlanchetS.Di BariP.The minimal scenario of leptogenesis20121412501210.1088/1367-2630/14/12/125012DavidsonS.NardiE.NirY.Leptogenesis20084664-510517710.1016/j.physrep.2008.06.0022-s2.0-50149119949FongC. S.NardiE.RiottoA.Leptogenesis in the universe201220125915830310.1155/2012/1583032-s2.0-84872763388FrèreJ.-M.HambuyeT.VertonenG.Is leptogenesis falsifiable at LHC?200920091, article 05110.1088/1126-6708/2009/01/051DhuriaM.HatiC.RangarajanR.SarkarU.Falsifying leptogenesis for a TeV scale WR± at the LHChttp://arxiv.org/abs/1503.07198CanetiL.ShaposhnikovM.Baryon asymmetry of the Universe in the NuMSM201020109, article 001CanettiL.DrewesM.FrossardT.ShaposhnikovM.Dark matter, baryogenesis and neutrino oscillations from right-handed neutrinos201387909300610.1103/physrevd.87.093006EnqvistK.MazumdarA.Cosmological consequences of MSSM flat directions20033803-49923410.1016/S0370-1573(03)00119-4ZBL1027.835332-s2.0-0038239848DineM.RandallL.ThomasS.Baryogenesis from flat directions of the supersymmetric standard model19964581-229132310.1016/0550-3213(95)00538-2DolgovA.KirilovaD.On the temperature of the boson condensate evaporation and the baryon asymmetry of the universe in the Affleck-Dine scenario19895016211629AnisimovA.DineM.Some issues in flat direction baryogenesis20016191–372974010.1016/s0550-3213(01)00550-82-s2.0-0035945045AllaverdiR.MazumdarA.Affleck-Dine condensate, late thermalization, and the gravitino problem20087804351110.1103/physrevd.78.043511DuttaB.SinhaK.Affleck-Dine baryogenesis in effective supergravity20108209500310.1103/physrevd.82.095003KawasakiM.NakayamaK.Affleck-Dine baryogenesis in anomaly-mediated supersymmetry breaking200720072, article 00210.1088/1475-7516/2007/02/002HigakiT.NakayamaK.SaikawaK.TakahashiT.YamaguchiM.Affleck-Dine baryogenesis with R-parity violation20149004500110.1103/physrevd.90.045001KasuyaS.KawassakiM.Towards the robustness of the Affleck-Dine baryogenesis2006746606350710.1103/PhysRevD.74.063507MarshD.Towards constraining Affleck-Dine baryogenesis20122012504110.1007/jhep05(2012)041DolgovA. D.KirilovaD.Production of particles by a variable scalar field199051172177DolgovA. D.KirilovaD.Baryon charge condensate and baryogenesis19911217ChizhovM. V.KirilovaD. P.Generation of the 128 h-1 Mpc periodicity of the universe in the scalar field condensate baryogenesis scenario1996101697510.1080/10556799608203247KirilovaD.PanayotovaM.The role of particle creation processes in the scalar condensate baryogenesis model200734s2330335KirilovaD. P.ChizhovM. V.Non-GUT baryogenesis and large-scale structure of the Universe2000314225626210.1046/j.1365-8711.2000.03193.x2-s2.0-0001079476ChizhovM.DolgovA.Baryogenesis and large-scale structure of the universe19923721-252152910.1016/0550-3213(92)90330-eDolgovA. D.KawasakiM.KevlishviliN.Inhomogeneous baryogenesis, cosmic antimatter, and dark matter20098071-222925010.1016/j.nuclphysb.2008.08.029ZBL1192.830762-s2.0-53949107377KirilovaD.Baryogenesis model predicting antimatter in the Universe200312240440810.1016/s0920-5632(03)80429-2KirilovaD.PanayotovaM.ValchanovT.Vast antimatter regions and scalar condensate baryogenesisProceedings of the 14th Rencontres de Blois ‘Matter-Antimatter Asymmetry’June 2003439KirilovaD.PanayotovaM.Inhomogeneous baryogenesis model and antimatter in the UniverseProceedings of the 8th Serbian-Bulgarian Astronomical Conference (VIII SBGAC '12)May 2012Leskovac, SerbiaKirilovaD.PanayotovaM.The dependence of the baryon asymmetry generation on the couplings of the baryon charge carrying scalar field20142045VilenkinA.FordL. H.Gravitational effects upon cosmological phase transitions19822661231124110.1103/physrevd.26.12312-s2.0-0000264088LindeA. D.Coleman-Weinberg theory and the new inflationary universe scenario1982114643143510.1016/0370-2693(82)90086-7MR669167BunchT. S.DaviesP. C. W.Quantum field theory in de Sitter space: renormalization by point-splitting1978360170011713410.1098/rspa.1978.0060StarobinskyA. A.Dynamics of phase transition in the new inflationary universe scenario and generation of perturbations19821173-417517810.1016/0370-2693(82)90541-xPressW. H.20012ndCambridge University Press