Varying constants driven baryogenesis

We study the spontaneous baryogenesis scenario in the early universe for three different frameworks of varying constants theories. We replace the constants by dynamical scalar fields playing the role of thermions. We first obtain the results for baryogenesis driven by the varying gravitational constant, $G$, as in the previous literature, then challenge the problem for varying fine structure constant $\alpha$ models as well as for varying speed of light $c$ models. We show that in each of these frameworks the current observational value of the baryon to entropy ratio, $\eta_B \sim 8.6 \cdot 10^{-11}$, can be obtained for large set of parameters of dynamical constants models as well as the decoupling temperature, and the characteristic cut-off length scale.


I. INTRODUCTION
The problem of an excess of the matter over the antimatter in the universe we observe, is one of the biggest mysteries of contemporary cosmology. Why do we live in the particle-filled universe and not in the anti-particlefilled universe is not an obvious issue, especially taking into account that the anti-particles are observed in the particle-filled universe, too. There are series of explanations to the mystery appealing to the earliest stages of the universe evolution such as to the Planck scale quantum gravity era [1] or even before that, as suggested for example in the context of the multiverse concept [2,3]. However, despite the fact, that it is generally agreed, that quantum gravity does not preserve any global quantum numbers (as is evident from the lost of baryon number in the process of a star collapse forming a black hole), inflationary expansion is considered to dilute any such matter-antimatter asymmetry. Then, one should look for the solution of the problem in the subsequent stages of the evolution of the universe. The first attempt to explain the problem was given in the renowned paper by Sakharov [4].
As it is often referred, he suggested three necessary conditions for the matter-antimatter asymmetry to occur in the universe: the baryon number B violation interactions have to appear; charge C and charge-parity CP violating particle processes have to be possible; de- * katarzyna.leszczynska@usz.edu.pl † mariusz.dabrowski@usz.edu.pl ‡ tomasz.denkiewicz@usz.edu.pl parture from thermal equilibrium to shift the densities of particles with respect to antiparticles have to be present. This idea of Sakharov was developed in many ways [5]. An idea of the spontaneous baryogenesis was later introduced by Cohen and Kaplan [6]. In fact, they challenged the third Sakharov's condition, i.e. the departure from thermal equilibrium, postulating instead a spontaneous breaking of CPT symmetry already at the thermal equilibrium which generates the shift of the energy of the baryons with respect to the energy of anti-baryons in the universe which is responsible for the baryon asymmetry. An exit from thermal equilibrium takes place at some decoupling temperature, but once generated, baryon asymmetry is frozen-in in this scenario. Baryogenesis is driven by a scalar field-the thermion-which decays after the baryon asymmetry is established [7,8]. The role of the thermion can also be played by some gravitationally motivated scalar such as Ricci or Gauss-Bonnet scalars and their combinations-such scenarios are called gravitational baryogeneses [9][10][11]. Baryogenesis in the context of other theories such as in Lorentz symmetry violating models has also been studied [12].
The baryon asymmetry problem is usually referred to the observational number η B , which is the ratio of the baryon number density, n B , to the entropy density, s (or the photon number density, n γ ). According to the latest measurement by the Planck satellite [13], the dimensionless baryon density Ω B h 2 = 0.02225 ± 0.00016 gives the baryon asymmetry equal to η B = n B s = (8.678 ± 0.062) · 10 −11 , (I. 1) where η B = 3.9 · 10 −9 Ω B h 2 . The entropy density and the photon number density are related by s ≈ 7.04n γ .
In this paper we concentrate on the spontaneous baryogenesis approach with the baryon asymmetry generating fields being motivated by the dynamical physical constants. The paper is organised as follows. In Section II we briefly sketch the idea of spontaneous baryogenesis. In Section III, which is the main body of this work, we discuss how to generate baryon asymmetry in spontaneous baryogenesis scenario where the role of a thermion is played by dynamical constants such as the varying gravitational constant G, varying fine structure constant α, and the varying speed of light c. In Section IV we summarise our results and give conclusions.

II. SPONTANEOUSLY GENERATED BARYON ASYMMETRY IN THE UNIVERSE
As it was mentioned in the Introduction, unlike the Sakharov baryogenesis, the spontaneous baryogenesis is based on two assumptions [6]: 1) baryon number violating interactions appear in thermal equilibrium; 2) CPT is not an exact symmetry of the early universe since its expansion violates Lorentz symmetry and a time-reversal. There is a relation between the Hubble parameter evolution and the size of CPT violation which can be tighten to the "effective" baryon number violating interactions. These being initially large, after the universe cools down and can be approximated by zero temperature, become gradually negligible so that baryon number violating interactions become CPT invariant and then allow the Lorentz invariant vacuum as we observe now.
The key point is to consider a scalar field ϕ-in the original approach called thermion [6]-which spontaneously breaks the baryon symmetry by a term in the action where λ is a characteristic cut-off length scale of the spontaneous baryogenesis model (l pl < λ < l GU T ; l pl is the Planck length and l GU T is the Grand Unified Theory length scale), J µ B is the baryon current, and the Greek indices run from 0 to 3. After integrating (II.2) by parts, one obtains which means that the baryon current J µ B cannot be conserved (or otherwise, ∂ µ J µ B cannot be zero). If it was conserved, the baryon number would be preserved and so there was no baryon asymmetry in the universe. The underlying idea here is to replace the term ∂ µ J µ B by some operator which violates the baryon number and additionally can also give rise to a decay of thermion field at late time to finally reach baryon conservation at the late universe.
Considering a homogeneous and isotropic Friedmann Universe ds 2 = −c 2 dt 2 + a 2 (t) dr 2 1 − kr 2 + r 2 dθ 2 + r 2 sin 2 θdϕ 2 , (II.4) (x ν = (x 0 , x 1 , x 2 , x 3 ) = (ct, r, θ, ϕ), d/dx 0 = d/(cdt), k = 0, ±1), and assuming that the baryon number density is we can write down (II.2) as with µ = (1/c)λ 2φ being the "effective" chemical potential, and (. . .) · = d(. . .)/dt. In fact, the term (II.6) describes a CPT violating interaction which leads to different spectra for baryons and antibaryons, but these spectra remain in thermal equilibrium with different thermal distributions. More precisely, the term (II.6) breaks first, the CP symmetry and then the time symmetry due to having a nonzero vev <φ > = 0, which finally leads to a CPT violation [10]. Now, having these spectra shifted, J µ B represents the baryon current of a particle/an antiparticle, where its null component describes the difference in the number density of particles and antiparticles: (II.7) The CPT symmetry ensures that particles and antiparticles equilibrate with the same thermal distribution, which is not a case when the symmetry is broken. Therefore, the interaction (II.2) or (II.3) contributes to the Einstein equation by the energy-momentum tensor made out of these baryon number violating terms and shifts the energy of baryons with respect to the energy of antibaryons of about 2µ B . This shift is then interpreted as a chemical potential, which enters the particle/antiparticle Hamiltonian through the term [7] where E B is the energy of a baryon and EB of an antibaryon. For the antiparticles the chemical potential is µB = −µ B . The thermodynamical quantities of some species "i"-the number density n i , the energy density, ε i = −ρ i c 2 , and the pressure density P i -are specified by the integrals over their distribution functions: where E is the energy, p is the momentum, g i is a number of the internal degrees of freedom, i.e. g i = 2 for a photon, and 1/(2π ) 3 is a unit size of the phase space.
The distribution function reads as [7,14]: and due to the homogeneity and isotropy of the Friedmann universe it does not depend on the spatial coordinates and the momentum direction, so f ( x, p) → f (p). The sign "+" stands here for fermions (Fermi-Dirac statistics), "−" for bosons (Bose-Einstein statistics), and k B is the Boltzmann constant. Given this, an excess of a baryon number over an antibaryon number (II.7) can be written as which by using (II.7) and (II.8) gives an approximate result for the particle-antiparticle excess as (II.14) The entropy density, s, for bosons (here: the radiation) is given by: and g * s is the effective number of degrees of freedom, which differs from g * present in the solution of the integral (II.10) for the energy density: When all the species have the same temperature and the equation of state may be approximated by p 1/3ρ i c 2 , these quantities appear equal, g * s = g * . Since above the temperature T ∼ 200 GeV all the particles are relativistic, we can find the value of g * = 106.75 by summing up their internal degrees of freedom [7,14]. By combining equations (II.8) and (II.14) we can write the final expression for the baryon asymmetry parameter: which has a dimension of (K/J) in SI units. Another parameter describing the preference of matter over antimatter is the baryon to photon number ratio, η Bγ = ∆n B /n γ . However, until the photon decoupling (T ∼ 0.3 eV) the photon density number n γ vary significantly throughout the epochs of the evolution of the Universe. For this reason, the entropy density s, which remains more or less constant at all energies, seems to be a better quantifier of the baryon asymmetry.

III. DYNAMICAL CONSTANTS DRIVEN BARYOGENESIS
An idea of varying physical constants is in a way analogous to the idea of running coupling constants in quantum field theory, i.e. that there is some interaction due to perhaps unknown physics, which causes these constants to vary in time and possibly in space. In practice, what one does is that one replaces the constants of nature by some physical fields, which have their own dynamics. The first fully quantitative framework for this was developed for varying gravitational constant (as a coupling constant of gravitational interaction) by Jordan [15] and Brans-Dicke [16]. They were motivated by the earlier Large Number Hypothesis of Dirac [17] being the consequence of even earlier ideas of Weyl [18] and Eddington [19]. Among the rich set of fundamental constants (for a review see Refs. [20][21][22][23]) the series of them are subject to dynamical studies. These are the gravitational constant G [16], the proton to electron mass ratio µ = m p /m e [24], the fine structure constant α = e 2 / c [25] ( is the Planck constant and c is the speed of light) and related to this charge of an electron e [26] or permittivity of vacuum 0 [27], and the velocity of light c [28].
Our main idea here is to have the scalar fields which are responsible for the dynamics of the constants such as G, α, and c to play the role of a thermion in the spontaneous baryogenesis scenario as described in Section II. In the following we will discuss all these three varying constants scenarios in that context.

A. Dynamical gravitational constant G driven baryogenesis
The action for (varying G) Brans-Dicke theory reads as [16]: where: where S BD is the standard Brans-Dicke action, S m is the matter action, φ is the Brans-Dicke scalar field, ω is a constant Brans-Dicke parameter, and we have added an extra term S B which describes the interaction responsible for baryogenesis [7]. It is worth to mention that we have used the notation for the action (III. 19) with the speed of light being c 3 rather than c 4 in front of the integral. Here we follow the notation of Refs. [29,30] compensating one c to be kept in the definition of the null coordinate x 0 = ct rather than x 0 = t as in most of the textbooks (e.g. [31]). We will come back to this problem in section III C where the models with varying speed of light are considered (cf. also the detailed discussion of Ref. [32]). The scalar field φ is related to the varying gravitational constant G as The action (III.18) varied with respect to the metric yields the field equations: where the tensors T i ν µ are given by and the equation of motion of the field φ takes the form: where Γ µ µγ are the Chrisstoffel connection coefficients. Assuming that the field (III.22) is homogeneous and isotropic we can write the Friedmann equation for the flat Universe as follows where H denotes the Hubble parameter, H φ =φ/φ is the rate of variation of G(t), ρ m is the matter mass density, and λ G is the characteristic cut-off scale for G-varying models.
In order to calculate (II.17), we need to solve the equation of motion (III.26), which for Friedmann metric takes the form: where w is an index of the barotropic equation of state p = wρ m c 2 , and p is the pressure. The value of the chemical potential (II.8) now reads as: Given this, we find (II.14) for the varying G case, which reads as: and in the next step, by inserting (III.30) into (III.28) we obtain:φ where β is constant with the dimension of We simplify (III.31) by making an assumption that during baryogenesis the temperature does not change significantly, i.e.Ṫ ≈ 0. We also consider the length scale λ G in the range l P l ∼ 10 −35 m < λ G < l GU T ∼ 10 −31 m. This allows us to check, that βT 2 1 for considered temperatures: 10 13 GeV < T < 10 15 GeV, and for 500 < ω < 40000. Consequently, the right hand side of (III.31) appears to be negligible and the equation of motion simplifies toφ + 3Hφ = 0 . (III.33) We describe the dynamics of the field φ by making an ansatz where a dimensionless parameter q is a measure of the variation of G(t) ≡ G, and a(t) ≡ a, a(t 0 ) ≡ a 0 are the scale factors at times t and t 0 . For q = 0, G(t) is equal to a currently measured value of the gravitational constant, which has the following solution: and allows us to find the Hubble parameter H: which then gives the values of φ and its derivativeφ as the functions of time: FIG. 1. The decoupling temperature, TD, as function of the parameter q, in the model with varying G using the currently measured value of the asymmetry, ηB 8.6 · 10 −11 . The thin line corresponds to ω = 1 000, the middle line corresponds to ω = 10 000, and the thick line to ω = 40 000. The plots were made for λG = 10 −19 Gev −1 , ain = 10 −25 , and Hin = 10 8 GeV.
The indices "in" denote the initial values of the scale factor a in , the Hubble parameter H in =ȧ in /a in and the field φ in = 1/G in at the beginning of baryogenesis, when t = t in . In order to define the temperature dependence of φ we combine equations (III.27), (II.16), (III.38), (III.39), and (III.40). Due to the fact that the right hand side of (III.28) is negligible, we feel excused to neglect the term cλ 2 G /6H φ J 0 B in the Friedmann equation. This gives a time-temperature relation in our model as which can be implemented into (III.40) and together with eq. (II.17) and (III.29) gives the final expression for the baryon asymmetry: The quadratic equation in (III.42) relates ω with the parameter q and gives a bound on the allowed values of the field φ. For ω > 0 the bound is given by (III.44) The limits for q for some specific ω has been listed in TA-BLE I. In order to be consistent with the current measurements, we take ω = 40 000 and the corresponding TABLE I. Limits for the parameter q for some specific values of ω (ω = −3/2 for conformal relativity, ω = −1 for superstring theory, ω = ∞ for Einstein gravity limit).  limit for q [33][34][35]. Despite the fact that ω could have been smaller in the early Universe, taking ω = 40 000 seems to be justified by the fact that it has small impact on the decoupling temperature T D in (III.42).
We limit ourselves to consider the negative values of q, only (G decreases during the evolution of the universe). The positive values of q would result in negative η B , hence a universe with an excess of antimatter on matter. The currently measured value of ω imposes even stronger bound onto the value of q and narrows the limit to the range (0; 0.0122) (see TABLE I). Nevertheless, the parameter ω does not have a strong influence onto T D for the values of q taken from this range. Its influence becomes more pronounced for smaller q (see Fig. 1). An interesting observation from (III.44) is that for the well-known from the literature case-the conformally invariant gravity [36], ω = −3/2-an allowed value of q is positive, so this theory seems to contradict observed baryon asymmetry. On the other hand, the limit of q for the low-energy superstring gravity [37], ω = −1, allows creation of the observed baryon asymmetry in a G-varying universe. Note, that in the general relativistic limit, ω → ∞, the parameter q vanishes, and according to (III.34), the value of G remains constant as it should be.
We have found that the Brans-Dicke baryogenesis gives a currently measured value of the baryon asymmetry, η B , for large range of the parameters q and T D . In order to be consistent with the results of BBN (−0.10 < ∆G/G < 0.13) [21]), we have taken q ∼ 10 −3 . This results in the decoupling temperature T D ∼ 10 14 GeV. The measurements of CMB indicate that q should be rather of the order of 10 −2 (−0.083 < ∆G/G < 0.095), which corresponds to a change of G between the recombination (z ≈ 10 3 ), and today (z = 0) [21]. In fact the parameter q characterising the dynamics of G can be calculated from the formula: where z G is a corresponding value of the redshift for which ∆G/G was measured.
We have performed calculations and plotted the results with the initial condition for the scale factor a in = 10 −25 and the corresponding Hubble parameter H in = 10 8 (calculated from ΛCDM). However, shifting the beginning of baryogenesis even from a in = 10 −25 to a in = 10 −30 only slightly changes η B . Therefore, the second most relevant parameter to drive baryogenesis is the fundamental length, λ G . In this case, a small change in the value of λ G results in a big change in the decoupling temperature, T D (see Fig. 2). The baryon to entropy ratio η B (III.42) as a function of the parameter q for three possible values of the temperature T D (see the plot on the left), as well as a function of temperature T D for three values of q (see the plot on the right) has been presented in Fig. 3.

B. Dynamical fine structure α driven baryogenesis
In this section we examine the application of the mechanism of spontaneous baryogenesis to the Bekenstein-Sandvik-Barrow-Magueijo (BSBM) [38] model of the varying fine structure constant α. Such models were first proposed by Teller [39], and later by Gamow [40], following the original path of the Large Number Hypothesis by Dirac [17]. A fully quantitative framework was developed by Bekenstein [26] in which a change in the fine structure constant α was fully identified with a variation of the constant electric charge, e 0 (cf. also Ref. [41]). By assuming that α can vary, we also assume that the electric charge become space-time dependent. This gives a path to a charge conservation, but maintains the Lorentz invariance, which is usually violated in the theories of varying α, where e and are kept constant, and c varies. The electric charge variability was introduced by defining a dimensionless scalar field, (x µ ), and as a consequence, e 0 was replaced by e = e 0 (x µ ). The electromagnetic tensor was then redefined to the form where the standard form of it can be restored for the constant . For simplicity, in [38,41] an auxiliary gauge potential, a µ = A µ , and the electromagnetic field strength tensor, f µν = F µν , were introduced, as well as a variable change: → ψ ≡ ln was performed. The field ψ in this model couples only to the electromagnetic energy, disturbing neither the strong, nor the electroweak charges, nor the particle masses.
The BSBM baryogenesis action is composed of where S g is the gravitational action, S em is the electromagnetic part of the theory with the kinetic term S ψ and S B is the baryogenesis term with the field ψ derivatively coupled to the baryon current J µ B . Similar to the original BSBM theory, the coupling constant Ω = c/λ, is a constant introduced for dimensional reason (J/m), where λ is considered the length scale of the electromagnetic part of the theory. The constant λ α is a cut-off length scale of the spontaneous baryogenesis model and is taken to be λ P l < λ α < λ GU T . The field ψ is given by: and is dimensionless. The field equations read as where the tensors T i ν µ are given by: (III.55) and the equation of motion of the field ψ is: The Friedmann equation for the Friedmann metric (II.4) and the homogeneous field ansatz ψ = ψ(t) reads as and the energy densities of ψ field, the electromagnetic field, and baryons respectively, are expressed by Similarly to the section III A, where the varying G baryogenesis was discussed, the interaction (III.50) violates the CPT symmetry. This results in a different thermal distributions for particles and antiparticles and contributes to the stress-energy tensor by (III.60), which can be understood as a chemical potential: This together with (II.17) leads to the baryon to entropy ratio in the form: (III.62) In order to calculate (III.62), we need to solve the equation of motion of the field ψ (III.56). For the pure radiation L em vanishes, so we can safely neglect this term. Since we assumed, that the field ψ is homogeneous and isotropic and only the null component of J µ B gives a contribution to the difference in the number densities, we can reformulate (III.56) to: When we insert (II.7), (II.14), and (III.61) into (III.63), we can try to estimate the value of the right hand side of this equation and its impact onto the evolution of ψ as: whereβ is a constant of unit K −2 : (III.65) We assume that the temperature of the Universe in a short period of baryogenesis did not change significantly (Ṫ ≈ 0). We also evaluate the value ofβT 2 , which for λ α = λ P l and ωλ of the order of few tens of MeV, is much smaller than one. For this reason we are excused to simplify (III.64) to: In order to solve (III.66) we need to describe the dynamics of α by an explicit dependence on the scale factor a(t).
We make the following ansatz for α: where a constant parameter m measures a change in α and the index "0" denotes the current values of the fine structure constant, α(t 0 ) = α 0 , and the scale factor, a(t 0 ) = a 0 . A scenario with no variation of α can be restored for q = 0. By inserting (III.67) into (III.66) we find that:ä The above equation can be integrated from t to t in , where time t in stands for the onset of baryogenesis. This gives where a(t in ) = a in , and H(t in ) = H in are, respectively, the initial value of the scale factor, and the corresponding Hubble parameter. The evolution of the scale factor (III.69) is not sensitive to the parameter m, nevertheless the baryon to entropy ratio exhibits a strong dependence on m. By using (III.51), (III.67) and (III.69) we can find that In order to write η B as a function of the temperature T , we combine the Friedmann equation (III.57) and the energy density (II.15) to yield and finally express (III.62) in terms of temperature as: We have found that it is possible to achieve the currently measured value of the baryon asymmetry, η B , in the BSBM model of baryogenesis, as well as in the model of varying G discussed in III A. A possible parameter space is presented in Fig. 4 for three different values of λ α . We have compared the ansatz (III.67) with the measurements of time variation of α to find the order of magnitude of the parameter m. We have found m to be in the range: where z α is a corresponding value of redshift for which ∆α/α has been measured. Using the bound from Ref. [25] we have decided to restrict m to be of the order of 10 −6 . This corresponds to the decoupling temperature T D ∼ 10 16 GeV. In our model m takes positive values only, which stands for the smaller α in the past. However, positive m can also be admitted according to the so-called α-dipole measurement [25]. Similarly to the model of Gdriven baryogenesis, any small change in the initial value of the scale factor a(t) at the moment of baryogenesis does not have any strong impact on η B . The sensitivity of η B increases with the growth of m. Again, the second most significant parameter is the length λ α , which was chosen to be of the order of 10 −19 GeV (see Fig. 4). The baryon to entropy ratio η B (III.73) as a function of the parameter m for three values of the decoupling temperature T D (see the plot on the left), as well as a function of the temperature T D for three values of m (see the plot on the right) has been shown in Fig. 5.

C. Dynamical speed of light c baryogenesis
Early ideas about varying speed of light c were even distributed by Einstein [42] and then many years later recalled by Petit [43] and Moffat [28,44]. Moffat developed a fully consistent theory which was designed to alternatively solve all the problems of standard cosmology which were originally resolved by the inflationary scenario [45]. Different types of varying speed of light models were also suggested by Albrecht and Magueijo [46], Barrow and Magueijo [47,48], and further developed by Magueijo [49,50]. These models are also useful to solve the standard cosmological problems such as the horizon problem, the flatness problem, the Λ−problem, and has recently been proposed to solve the singularity problem [51]. Another different class of varying speed of light models was given by Avelino and Martins [52]. All the above models have recently been subject to statistical evaluation against observational data [53] showing the preference of Moffat's models, which we have selected to study in the context of baryogenesis.
Here we combine the most recent Moffat's approach [54] with the theory of the spontaneous baryogenesis.
The appropriate action is made up of three terms: The action (III.76) is the gravitational action with the field Φ coupled to the curvature and the kinetic term with a constant κ. We also introduce the matter term S m , since at the moment of baryogenesis the Universe was filled-in with radiation. The baryon asymmetry is produced by the interaction term (III.78). As in the previous chapters, the length λ c is the cut-off length of the applicability of the theory, J µ B is the baryon current, and its null component describes a difference in the particle and the antiparticle number densities (cf. eq. (II.7)). Similarly to the section III A, we follow the notation of Refs. [29,30,32] for the Einstein-Hilbert action.
This means that we take x ν = (x 0 , x 1 , x 2 , x 3 ) and so (II.4) is replaced by (III.79) where x 0 = c(t)t (cf. the discussion of Appendix A in Ref. [46]). As a consequence, the dynamics of the speed of light field Φ is given by: (III.80) which differs from the Moffat's definition of the Φ field in [28] and also in other references which take Φ = c 4 [46,48]. However, both formulations are equivalent. In fact, the original Moffat's theory consists of the action representing the dynamics of four scalar fields. In a later paper [54] a vector field was driving the spontaneous violation of SO(3, 1) Lorentz invariance, while a dimensionless scalar field, minimally coupled to gravity, was responsible for quantum primordial fluctuations. Nevertheless, unlike in [28], we are dealing with a small Lorentz violation, and consequently with a small change in the speed of light c. In this paper we do not intend to explain the fast exponential expansion of the early universe to make it alternative to inflation and for this reason we have dropped the part of the Moffat's theory, which exhibits the strong Lorentz symmetry breaking and therefore, the large change in c (c ≈ 10 28 c 0 , where c 0 is the current value of the speed of light). The Lagrangian for the quantum primordial fluctuations has not been included, either. At the moment of baryon asymmetry generation, this term is not relevant anymore, and can safely be neglected.
The variation of (III.75) with respect to the metric g µν leads to the field equations: where the tensors T i ν µ are given by: The equation of motion of the field Φ takes the form: where Γ µ µγ are the Christoffel connection coefficients, and T m is the trace of the radiation energy-momentum tensor. We assume the barothropic equation of state of the fluid p = wρc 2 , which gives the trace T m = −ρc 2 (1− 3w) and vanishes for pure radiation field, w = 1/3. The equation (III.85) reads then as: where ( ) stands for the derivative with respect to the coordinate x 0 ,H = a /a is the Hubble parameter, andH Φ = Φ /Φ. Both,H andH Φ are of the unit m −2 , instead the usual s −2 . This is a consequence of a chosen definition of the action (III.75) and of coordinates in (III.79). The spontaneous baryogenesis occurs when CPT symmetry is broken in the Universe, which is in thermal equilibrium. This leads to a conclusion that particles, as well as the antiparticles are in thermodynamical equilibrium, but possess different energies. This is what we call the energy shift and can find it by investigating the contribution of (III.78) to the total energy density. The chemical potential takes the form: Unlike in the previous sections III A and III B of varying G and varying α, here not only the derivative of the field Φ enters the chemical potential µ B , but also the field Φ itself. The ratio of the baryon asymmetry to the entropy density is given by: which together with (III.87) yields: (III.89) In order to find the value of (III.89), we need to solve the equation of motion (III.86). First, we find J 0 B , which is the matter-antimatter excess: and then insert (III.90) into (III.86). This gives a relation, which connects Φ and its derivatives with the temperature T : where χ is a constant of the unit K 2 ms −1 : It has been checked, that the right hand side of (III.91) is small, and so it can safely be neglected. A resulting simplified equation of motion is then: In order to solve (III.93), we make an ansatz for the field Φ as follows: where x 0 0 = c 0 t 0 and n is a parameter, which indicates the variation in c (presumably small since we deal with approximate Lorentz symmetry). In the limit n → 0, a currently measured value of the speed o light c 0 is restored and the field is just equal to Φ = c 3 0 . We will denote a(x 0 ) ≡ a and a(x 0 0 ) ≡ a 0 , later on. The solution of (III.93) is: where Φ in = c 3 0 (a in /a 0 ) 3n is the initial value of the field at the beginning of baryogenesis, x 0 in = c in t in , , andH is the initial value of the Hubble parameter. Taking the derivative of (III.95), we find the final expression for the baryon asymmetry η B : .
The relation (III.96) is a function of the temperature T and the x 0 -coordinate. In order to express η B as a function of temperature only, we use the temperature dependent expression for the energy density of relativistic particles (II.16), which depends on temperature, as well as on the speed of light c at a given moment. For this reason, we find it reasonable to replace c 3 in the denominator by the field Φ. By combining the modified equation (II.16) with the Friedmann equation for a flat universe: (III.97) where: we can find an approximate time-temperature relation in the varying c models for the radiation dominated epoch. In order to do so, we have neglected ε B in (III.97) treating this term as a small perturbation on the background of the main fluid, which is radiation (similarly as we did in (III.93)). This leads to: The limits for the parameter n are as follows for κ > 0, and n ∈ −∞, 3 + √ 9 + 6κ 3κ ∪ 3 − √ 9 + 6κ 3κ , +∞ , (III.103) for κ < 0. For κ = 0, n > −1/3. Finally, we find that the baryon asymmetry reads as In summary, for the Moffat's model of the varying speed of light, c, we have found the parameter space, for which the desirable asymmetry is possible for a given range of the parameter n and the temperature T D (see the Fig. 6). We have found a relation between n and the constant κ from (III.101). The limits on the n values has been shown in the Table II. We have noticed, that the limits corresponding to the negative and positive κ partially overlap and have decided to proceed the calculation only for positive κ, even though there are no observational or experimental bounds on its value. Nevertheless, for the chosen order of magnitude of the parameter n, the influence of κ onto the final result is negligible. However, κ becomes more relevant, when bigger n is taken into account. In order to estimate an order of magnitude of n, we have assumed that any possible c variation would find its manifestation in the variation of the fine structure constant and thus it would become visible in the measurement of the α-time variation. By comparison of the value of ∆α/α with the ansatz (III.67) we have achieved the following expression for n: n = log 1+zα 1 + ∆α α . (III.105) Calculated limits on the parameter n have turned to be of the order of ∼ 10 −6 . The corresponding temperature for baryogenesis is about 10 16 GeV. However, the value of n does not need to be necessarily compared with the results for time variation of α, and one could consider even a bigger change of the speed of light. Our model favours the positive values of n, and thus the increase of the speed of light. This is not in the spirit of the varying speed of light models, which solve the basic cosmological problems and stand as an alternative to the inflation theories. In Fig. 7 we have shown the baryon to entropy ratio η B (III.104) as a function of the parameter n for three possible values of the temperature T D (see the plot on the left) and as a function of the temperature T D for three possible values of n (see the plot on the right).

IV. RESULTS AND CONCLUSIONS
In this paper we have investigated the scalar fields for the dynamical constants: the gravitational constant G, the fine structure constant α, and the speed of light c, which, as we have assumed, could drive the baryogenesis in the universe. The spontaneous baryogenesis model was investigated, in which the baryon number violating processes occur in thermal equilibrium, while the Universe allows a period of CPT symmetry breaking. We have formulated and solved the dynamical equations for the scalar fields corresponding to varying G, α, and c, acting as thermions. We have applied some special ansätze for the scale factor for each of the fields, as given in (III.34), (III.67), and (III.94), which related the scalar fields with the evolution of the scale factor, and the resulting parameters q, m, and n which measured the degree of variability of G, α, and c, accordingly. We have calculated the cosmological equations and used them to find the relation between the time t and the temperature T in the radiation dominated epoch. We have used the standard statistical mechanics tools in order to introduce the temperature dependence into the fields φ(G(x 0 )), ψ(α(x 0 )),  Φ(c(x 0 )) (where the coordinate x 0 = ct), and in order to calculate the baryon asymmetry ratio η B of the net number density of baryons and antibaryons to the entropy density of photons. As a result of our calculations, we have obtained similar conclusion as in the previous literature i.e. that varying G can drive baryogenesis in the universe. Our main new results (not yet considered in the literature) are obtained for varying fine structure constant α models, as well as for varying speed of light c models. We have shown that in each of these frameworks the current observational value of the baryon to entropy ratio η B ∼ 8.6 × 10 −11 can be obtained for large set of parameters q, m, n, as well as the decoupling temperature T D , and the characteristic cut-off length scale λ. This means that not only varying-G-driven baryogenesis is possible, but also varying-α-driven and varying-c-driven baryogeneses are admissible.