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 α\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha $$\end{document} 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, ηB∼8.6·10-11\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\eta _B \sim 8.6 \cdot 10^{-11}$$\end{document}, 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.


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 particlefilled universe and not in the anti-particle-filled 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 a e-mail: katarzyna.leszczynska@usz.edu.pl b e-mail: mariusz.dabrowski@usz.edu.pl c e-mail: tomasz.denkiewicz@usz.edu.pl 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 C P violating particle processes have to be possible; departure 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,7]. 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 [8,9]. 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 [10][11][12][13][14][15][16]. Baryogenesis in the context of other theories such as in Lorentz symmetry violating models has also been studied [17].
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 [18], 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 , (1.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 Sect. 2 we briefly sketch the idea of spontaneous baryogenesis. In Sect. 3, 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 Sect. 4 we summarise our results and give conclusions.

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,7]: (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,7] -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 (2.1) 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
where n B describes the difference in the number density of particles and antiparticles: In fact, the term (2.4) describes a CPT violating interaction which leads to different spectra for baryons and antibaryons. More precisely, the term (2.4) breaks first, the CP symmetry and then the time symmetry due to having a nonzero vev <φ > = 0, which finally leads to a CPT violation [11][12][13].
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 (2.1) or (2.2) 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 [8] 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πh) 3 is a unit size of the phase space. The distribution function reads as [8,19]: 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 (2.5) can be written as which by using (2.5) and (2.6) gives an approximate result for the particle-antiparticle excess as 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 (2.8) 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 [8,19].
By combining Eqs. (2.6) and (2.12) 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.

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 [20] and Brans-Dicke [21]. They were motivated by the earlier Large Number Hypothesis of Dirac [22,23] being the consequence of even earlier ideas of Weyl [24] and Eddington [25]. Among the rich set of fundamental constants (for a review see Refs. [26][27][28][29]) the series of them are subject to dynamical studies. These are the gravitational constant G [21], the proton to electron mass ratio μ = m p /m e [30], the fine structure constant α = e 2 /hc [31,32] (h is the Planck constant and c is the speed of light) and related to this charge of an electron e [33] or permittivity of vacuum 0 [34], and the velocity of light c [35]. Though one usually considers the dynamics of the constants separately, the models in which two of the constants vary instantaneously have also been considered. Out of them the most natural are modified varying both G and c models [36][37][38] since these constants show up together in the Einstein-Hilbert action for gravity and in the Einstein field equations. In fact, they can be classified as an extension of Brans-Dicke models into a varying c case. Another extension of this type which is based on Brans-Dicke model are varying both G and α models [39]. On the other hand, varying both α and c models would not perhaps make so reasonable because α and c are related via the definition of the fine structure constant and the effects of changes of these constants would have to be indicated separately in the construction.
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 Sect. 2. In the following we will discuss all these three varying constants scenarios in that context.

Dynamical gravitational constant G driven baryogenesis
The action for (varying G) Brans-Dicke theory reads as [21]: where: 3) where S B D is the standard Brans-Dicke action, S m is the matter action, φ is the Brans-Dicke scalar field, ω is a constant Brans-Dicke parameter. We have added an extra term, S B , which describes the interaction responsible for baryogenesis [8]. The quantity λ G is the characteristic cut-off scale for G-varying baryogenesis models. It is worth to mention that we have used the notation for the action (3.2) with the speed of light being c 3 rather than c 4 in front of the integral. Here we follow the notation of Refs. [40,41] 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. [42]). We will come back to this problem in Sect. 3.3 where the models with varying speed of light are considered (cf. also the detailed discussion of Ref. [43]). The scalar field φ is related to the varying gravitational constant G as .
The action (3.1) varied with respect to the metric yields the field equations: where the tensors T i ν μ are given by 8) and the equation of motion of the field φ takes the form: where μ μγ are the Christoffel connection coefficients.
Assuming that the field (3.5) 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) and ρ m is the matter energy density with corresponding pressure, p m . The acceleration equation is given bÿ In order to calculate (2.15), we need to solve the equation of motion (3.9), 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 main problem with the set of Eqs. (3.10)-(3.12) (compare Ref. [44] for example) is that in radiation dominated universe the first term on the right-hand side of Eq. (3.12) vanishes. This, after additionally neglecting the last term of baryogenesis which is in fact small, leads to a pure scalar field (or stiff-fluid) domination with a simple integral (3.13) Despite this solves easily, still the solutions of the whole set of equations for a(t) and φ(t) are non-trivial. In fact, one can postulate the power-law solutions which would include both an early universe "scalar field domination" and late universe radiation domination together. However, as it has been shown in Ref. [45] that the unique power-law solutions which allow current acceleration of the universe a ∝ t 4/3 , φ ∝ t −2 are possible for the dust ( p = 0) models if one also adds a specific scalar field potential. Because the matter applied is dust, then the relation (3.13) is modified accordingly.
In fact, if one neglected radiation also in the Friedmann equation (3.10), then one would get a unique solution which would pick up some specific powers in the scale factor and scalar field power-law time dependence. However, in our case we deal with the early universe and so radiation is the crucial component. Because of that we need to rely on the relation (3.13), though modified slightly by the baryogenesis term which is pretty small. The solutions which in fact keep relation (3.13) valid but also include radiation have been studied in Ref. [46]. In particular, it was found that there exist two regimes in which simple power-law solutions (with radiation and the scalar field present) exist. One of them applies close to a big-bang a → 0, where the scalar field is dominating (behaving as a stiff-fluid), and another to the late time evolution a → ∞, when the radiation comes to dominate. The first solution gives a simple power law for the scale factor a ∝ t 1/3 , while the second gives standard radiationdominated power law behaviour a ∝ t 1/2 .
What is crucial here is that there exist solutions fulfilling set of Eqs. (3.10)-(3.12) which include radiation and asymptote from the power-law solution a ∝ t 1/3 to possibly another power law solution a ∝ t 1/2 . They can be parametrised by the values of some extra parameter which takes some specific value for the asymptotic a ∝ t 1/3 solution. In other words, we can consider the solutions which are in stiff-fluid regime, but which are slightly modified by the presence of radiation. Such an approach have been applied to Brans-Dicke theory already in Refs. [47][48][49] though only some simple examples of baryogenesis out of the whole set of admissible values of the extra parameter (in Ref. [47] parameter n) have been studied.
In the following we will explain the above approach step by step.
We start with the value of the chemical potential (2.6) which now reads as: (3.14) Given this, we find (2.12) for the varying G case, which is: and in the next step, by inserting (3.15) into (3.12) we obtain the modified with baryogenesis term equation of motion for the scalar field: The Eq. (3.16) can be solved, but as it can be proven numerically, the contribution from the right-hand side is small because during baryogenesis the temperature does not change significantly, i.e.Ṫ ≈ 0. Besides, for the length scale λ G in the range l Pl ∼ 10 −35 m < λ G < l GU T ∼ 10 −31 m, βT 2 1 for considered temperatures 10 13 GeV < T < 10 15 GeV, and for 500 < ω < 40,000. Consequently, the right hand side of (3.16) can be neglected and the equation of motion simplifies tö which immediately gives the solution (3.13). Without any loss of generality in looking for the power law solutions we will describe the dynamics of the field φ by making an ansatz for the variability of the gravitational constant G: 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, G(t 0 ) = G 0 . The parameter q (which is equivalent to a parameter n = 1/(3 − q) of Ref. [47]) also measures the deviation from a ∝ t 1/3 as mentioned earlier. Applying the 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: 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 . It is worth mentioning that in the limit q → 0 the above solutions give an asymptotic early time behaviour a ∝ t 1/3 as it should be following the work of Ref. [46]. In fact, the scale factor (3.21) is a superposition of a solution for pure radiation and for pure stiff fluid which mimics the scalar field. In the above mentioned limit q → 0, the solutions (3.21) and (3.24) solve simultaneously the whole set of the field Eqs. (3.10), (3.11), (3.18), when satisfied: for the radiation energy density taking form in the Brans-Dicke theory: Notice that in order to get (3.18) the contribution (3.4) from baryogenesis term has been neglected in the Friedmann equation for the same reason for which we neglected it in the equation of motion (3.18). In fact, if we have also dropped radiation contribution in the Friedmann equation so ρ m0 = 0, we would obtain the condition linking the values of q-parameter and ω: which would restrict the freedom of choice of q in a similar way as in the dust case considered in Ref. [45].
In order to define the temperature dependence of φ we combine Eqs. (3.10), (2.14), (3.23), (3.24), and (3.25). Due to the fact that the right hand side of (3.12) 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 (3.25) and together with eq. (2.15) and (3.14) gives the final expression for the baryon asymmetry: The quadratic equation in (3.28) relates ω with the parameter q and gives a bound on the allowed values of the field φ.
For ω > 0 the bound is given by while for ω < 0 it reads as The limits for q for some specific ω has been listed in Table 1.
In order to be consistent with the current measurements, we take ω = 40,000 and the corresponding limit for q [50][51][52].
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 (3.29).
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 Table 1 Limits for the parameter q for some specific values of ω (ω = −3/2 for conformal relativity, ω = −1 for superstring theory, ω = ∞ for Einstein gravity limit) The decoupling temperature, T D , 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 ω = 1000, the middle line corresponds to ω = 10,000, and the thick line to ω = 40,000. The plots were made for λ G = 10 −19 Gev −1 , a in = 10 −25 , and H in = 10 8 GeV value of q and narrows the limit to the range (0; 0.0122) (see Table 1). 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 (3.31) is that for the well-known from the literature case -the conformally invariant gravity [53], ω = −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 [54], ω = −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 (3.19), 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) [27]), 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) [27]. 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 (3.29) 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.

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) [55] model of the varying fine structure constant α. Such models were first proposed by Teller [56], and later by Gamow [57], following the original path of the Large Number Hypothesis by Dirac [22,23]. A fully quantitative framework was developed by Bekenstein [33] 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. [58]). By assuming that α can vary, we also assume that the electric charge become spacetime 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 andh 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 [55,58] 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 S = S g + S ψ + S em + S B (3.33) and 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 . Similar to the original BSBM theory, the coupling constant =hc/λ, 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 λ Pl < λ α < λ GU T . The field ψ is given by: and is dimensionless. The field equations read as where the tensors T i ν μ are given by: 42) and the equation of motion of the field ψ is: The Friedmann equation for the flat Friedmann metric (2.3) and the homogeneous field ansatz ψ = ψ(t) reads as where ρ em is the electromagnetic field energy density which will be, later on, re-scaled as follows: The acceleration equation is given by: where p em is the electromagnetic pressure which we re-scale as: Similarly to the Sect. 3.1, where the varying G baryogenesis was discussed, the interaction (3.37) violates the CPT symmetry. This results in a different thermal distributions for particles and antiparticles and contributes to the stress-energy tensor. this contribution may be understood as a chemical potential: This together with (2.15) leads to the baryon to entropy ratio in the form: . (3.48) In order to calculate (3.48), we need to solve the equation of motion of the field ψ (3.43). 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 (3.43) to: When we insert (2.5), (2.12), and (3.47) into (3.49), 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 : (3.51) 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 λ α = λ Pl and ωλ of the order of few tens of MeV, is much smaller than one. For this reason we are excused to simplify (3.50) to: In order to solve (3.52) 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 m = 0. By inserting (3.53) into (3.52) 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 scale factor (3.55) differs from the one which is expected for the radiation dominated Universe. The presence of the scalar field ψ(t) shifts its value from a ∼ t 1/2 to a ∼ t 1/3 , this means that the solution (3.55) scales like a solution for the stiff fluid. By using (3.38), (3.53) and (3.55) we can find that In the limit m → 0 the solutions (3.55) and (3.56) solve simultaneously all the field equations when: , (3.58) for ρ em0 being a positive constant in the expression for the energy density of the stiff fluid: (3.59) In order to write η B as a function of the temperature T , we combine the Friedmann equation (3.44) and the energy density (2.13) to yield 60) and finally express (3.48) in terms of temperature as: (3.61) 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 3.1. A possible parameter space is presented in Fig. 4 for three different values of λ α . We have compared the ansatz (3.53) with the measurements of time variation of α where z α is a corresponding value of redshift for which α/α has been measured. Using the bound from Ref. [31,32] 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 [31,32]. Similarly to the model of G-driven 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 (3.61) 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.

Dynamical speed of light c baryogenesis
Early ideas about varying speed of light c were even distributed by Einstein [59] and then many years later recalled by Petit [60][61][62] and Moffat [35,63]. 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 [64]. Different types of varying speed of light models were also suggested by Albrecht and Magueijo [65], Barrow and Magueijo [36,66], and further developed by Magueijo [67][68][69][70]. These models are also useful to solve the standard cosmological problems such as the horizon problem, the flatness problem, theproblem, and has recently been proposed to solve the singularity problem [71]. Another different class of varying speed of light models was given by Avelino and Martins [72,73]. All the above models have recently been subject to statistical evaluation against observational data [74] 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 [75] with the theory of the spontaneous baryogenesis. The appropriate action is made up of three terms:

65)
The action (3.64) 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 (3.66). 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. (2.5)). Similarly to the Sect. 3.1, we follow the notation of Refs. [40,41,43] for the Einstein-Hilbert action.
This means that we take x ν = (x 0 , x 1 , x 2 , x 3 ) and so (2.3) is replaced by where x 0 = c(t)t (cf. the discussion of Appendix A in Ref. [65]). As a consequence, the dynamics of the speed of light field is given by: which differs from the Moffat's definition of the field in [35] and also in other references which take = c 4 [36,65]. 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 [75] 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 [35], 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 (3.63) 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: , (3.73) 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 Eq. (3.73) 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 (3.63) and of coordinates in (3.67). 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 (3.66) to the total energy density. The chemical potential takes the form: Unlike in the previous Sects. 3.1 and 3.2 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 (3.75) yields: In order to find the value of (3.77), we need to solve the equation of motion (3.74). First, we find J 0 B , which is the matter-antimatter excess: and then insert (3.78) into (3.74). 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 (3.79) is small, and so it can safely be neglected. A resulting simplified equation of motion is then: In order to solve (3.81), 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 (3.81) 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. The solution (3.83) solve the full set of the field equations in the limit n → 0 when: forρ m0 being a positive constant in the expression for the energy density of the stiff fluid: (3.85) Taking the derivative of (3.83), we find the final expression for the baryon asymmetry η B : (3.86) The relation (3.86) 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 (2.14), 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 Eq. (2.14) with the Friedmann equation for a flat universe: 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 (3.87) treating this term as a small perturbation on the background of the main fluid, which is radiation (similarly as we did in (3.81)). This leads to: (3.91) The limits for the parameter n are as follows for κ > 0, and for κ < 0. For κ = 0, n > −1/3. Finally, we find that the baryon asymmetry reads as (3.94) 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 (3.91). The limits on the n values has been shown in the Table 2. 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 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 (3.94) 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).

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 (3.19), (3.53), and (3.82), 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.
It is advisable to note that there exist models in which two of the three considered in this paper constants vary simultaneously [36][37][38][39]. However, in the most interesting case of varying G and α it has been shown that an overall evolution of the universe is determined by G and follows Brans-Dicke model so and extra influence of α on baryogenesis is not expected to be large. It is then expected that similar small effect of c-variability would remain in both varying G and c models. The detailed quantitative considerations of such models will be considered in some future work. Vincenzo Salzano and Tomáš Husek for discussions and valuable suggestions.

Data Availability Statement
This manuscript has no associated data or the data will not be deposited [Authors' comment: This is a theoretical study and no experimental data has been listed.] Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecomm ons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. Funded by SCOAP 3 .