Curvature and thermal corrections in tree-level CPT-Violating Leptogenesis

In a model for leptogenesis based on spontaneous breaking of Lorentz and CPT\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {CPT}$$\end{document} symmetry [1–3], we examine the consistency of using the approximation of plane-wave solutions for a free spin-12\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\frac{1}{2}$$\end{document} Dirac (or Majorana) fermion field propagating in a Friedmann–Lemai^\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hat{\mathrm{i}}$$\end{document}tre–Robertson–Walker space time augmented with a cosmic time-dependent (or, equivalently, a temperature-dependent) Kalb–Ramond (KR) background. For the range of parameters relevant for leptogenesis, our analysis fully justifies the use of plane-wave solutions in our study of leptogenesis with Boltzmann equations; any corrections induced by space-time-curvature are negligible. We also elaborate further on how the lepton asymmetry is communicated to the Baryon sector. We demonstrate that the KR background (KRB) does not contribute to the anomaly equations that determine the baryon asymmetry (a) through an explicit evaluation of a triangle Feynman graph and (b) indirectly, on topological grounds, by identifying the KRB as torsion (in the effective string-inspired low energy gravitational field theory).


Motivation and summary
In Refs. [1][2][3] we proposed and discussed a new scenario for leptogenesis induced by an axial background vector field that violates spontaneously Lorentz and CPT (C(charge conjugation), P(parity) and T (time)) symmetry [4]. In stringinspired models such backgrounds might be provided by the spin-one antisymmetric tensor Kalb-Ramond (KR) field [5], part of the massless gravitational string multiplet [6]. Our model for leptogenesis involves heavy sterile Majorana righthanded neutrinos (RHN), which have tree-level decays into lepton and Higgs particles (of the Standard Model (SM)) and their antiparticles to produce a lepton asymmetry L. When the universe is at a temperature T , a e-mail: nikolaos.mavromatos@kcl.ac.uk (corresponding author) where s is the entropy density of the universe and s ∝ T 3 [7]; m N is the RHN mass; z ≡ m N T ; z D ≡ m N /T D 1, with T D the decoupling temperature of RHN; q is a numerical coefficient of order O(1) [2,3]. 1 The constant has mass dimension +1, which equals the temporal component of the Lorentz-(LV) and CPT Violating (CPTV) axial background B 0 evaluated at a decoupling temperature T = T D . In stringinspired cosmological models of [1][2][3] for four-dimensional space time , B 0 is given by the gradient form where b(t) is the massless KR axion field and t is the cosmic time. The analysis of [1][2][3] assumes that, at temperatures near decoupling, one has so that the lepton asymmetry is evaluated to leading order in an expansion in powers of B 0 /m N . For f (z) = 1, B 0 is constant in the local Friedmann-Lemaître-Robertson-Walker (FLRW) frame [1,2]. In [3] we discussed microscopic models for CPTV leptogenesis for which f (z) = z −3 [3] and B 0 varies slowly with T as We shall concentrate on this scaling with temperature in this work. In the model of [3], we took in order for the lepton asymmetry in (1) to have the phenomenologically required [7,8] value L/s ∼ 8 × 10 −11 . This is consistent with (3).
The result (1) for the lepton asymmetry is obtained on using the standard formalism of Boltzmann equations [7] for leptogenesis. The quantum field-theoretic scattering amplitudes in the collision integral in Boltzmann equations were evaluated approximately, ignoring both space-time curvature effects and variation (4) of B 0 with T (or, equivalently cosmic time) [1][2][3]. Consequently we used plane-wave solutions for spinors in evaluating the amplitudes corresponding to the decay of RHN into SM particles. However, at a space-time point x μ in a curved manifold, plane-wave solutions of Dirac or Majorana equations exist only on the tangent space at that point. The use of plane-wave solutions and dispersion relations is thus an approximation, which ignores effects of curvature. Motivated by the current cosmological data [8], we have taken [1][2][3], the manifold to be that of an expanding universe, with a spatially-flat FLRW metric, corresponding to the line element: with x i , i = 1, 2, 3, Cartesian spatial coordinates, t the FLRW time coordinate, and a(t) the scale factor of the universe in units of today's scale factor a 0 . The curvature of the manifold has components proportional to ( d dt a(t)) 2 and Hence, because of the explicit time-dependence in the Dirac (or Majorana) equation, it is important to check that curved space-time effects and the variation of B 0 with t have been consistently accounted for in arriving at (1). During the radiation era of the early universe, when leptogenesis takes place in the scenario of [1][2][3], the scale factor of the universe scales as a(t) rad ∼ t 1/2 ∼ 1/T, (7) and thus, for a spatially-flat FLRW universe, the scalar spacetime curvature (R = 6 ä a + ȧ a 2 , at high temperatures of relevance to leptogenesis [3], exhibits a scaling with T (∼ T 4 ) comparable to that from B 0 (T ) (4). It is necessary to examine in detail whether such temperature scaling affects significantly the Boltzmann analysis of [3] which leads to the lepton asymmetry (1). In this work we shall demonstrate that the expansion of the FLRW universe does not affect the results of [3] for the lepton asymmetry. Our model for leptogenesis requires us to take into account only tree-level decays of RHN to SM particles for the generation of the lepton asymmetry (1); curvature effects will enter through the solution for the spinors, which will be modified compared to the flat spacetime case by terms proportional to powers of the Hubble parameter. Energy-momentum dispersion relations for the various modes will also receive such corrections.
We will present a systematic derivation of curvatureinduced corrections to plane-wave solutions of the Dirac (and Majorana) equation in an axial vector background given in (2). For the range (5) of the parameters of the model, we will show that the corresponding corrections to the planewave solutions of the Dirac (and Majorana) equations for the spinors are negligible . Our derivation extends the analysis of [9] to the standard Dirac equation in both a curved space time and an axial vector background. Such a perturbative analysis is applicable to space times which vary slowly in time, as is the case for spatially flat FLRW space time in the KR background (4) during the era of radiation domination.
Once we have leptogenesis, we use it to induce baryogenesis [1]. The lepton asymmetry generated by the KR background (4) is communicated to the baryon sector via sphaleron processes [10,11] in the SM sector. Sphaleron processes preserve the difference B − L between baryon (B) and lepton (L) numbers [12]. This is the route to baryogenesis in the conventional leptogenesis scenario [13]. However we need to check that the presence of the KR background B 0 (known to play the rôle of totally antisymmetric torsion [14][15][16] in string theories) does not affect [17,18] the anomaly equations [19] for the baryon and lepton numbers needed in the route [12] to baryogengesis.
The structure of our article is the following: In Sect.2 we discuss how the expansion of the universe and the KR background affects the collision terms of the Boltzmann equations used in the leptogenesis scenario of [1][2][3]. We also compare our study with recent results on Boltzmann equations in curved space-times [20,21].
In Sect.3, we obtain systematic corrections to plane-wave solutions of the Dirac equation (in Sect. 3.1) and of the Majorana equation (in Sect. 3.2) on a spatially-flat FLRW space time in the presence of the KR background (4). The results are similar in the two cases.
In Sect.3.3, for the parameter range (5), we demonstrate that any space-time curvature corrections to the flat spacetime result for the Boltzmann collision term are negligible; hence the conclusions in [3] remain unaffected. We provide a further check on the consistency of our calculation by relating the scattering amplitudes in the Boltzmann collision term, to the proper polarisation spinors for the Hermitian Hamiltonian, associated with the relativistic equation of motion for fermions in time-dependent metrics [22,23].
In Sect.4 we discuss in detail how the lepton asymmetry generated in our CPT violating leptogenesis scenario communicates to the baryon sector, via sphaleron processes in the Standard Model sector of the theory. Special attention is paid to discussing some properties of the KR background that are crucial to this effect, namely its non contribution to the baryon-and lepton-number anomaly equations.
Conclusions and outlook are given in Sect.5. Technical aspects of our approach are given in several Appendices.
Specifically, in Appendix A we set up our notation and conventions, and discuss some formal properties of the Dirac equation in (spatially flat) FLRW expanding Universe spacetimes, in the presence of axial backgrounds of relevance to the leptogenesis scenario of [1][2][3]. In Appendix B we show, following [22], that Hermiticity of the associated Hamiltonian is ensured upon taking proper account of (space-time curvature) effects, proportional to time derivatives of the metric. This procedure defines the appropriate polarisation spinors to enter the Boltzmann collision term, and justifies the mathematical self consistency of our model for leptogenesis. In Appendix C, we describe the details of the derivation of the (adiabatic) space-time curvature corrections to the plane-wave solutions of the Dirac equation in an expanding universe, expressed in a perturbative expansion in powers of the Hubble parameter H . In Appendix D, we discuss some thermodynamical aspects of sphaleron-induced baryogenesis, which completes our discussion in Sect.4 by incorporating high temperature effects properly. Finally, in AppendixE, we discuss a topological approach to demonstrating the noncontribution of the Kalb-Ramond torsion to the anomalies, which is of relevance to our baryogenesis considerations in Sect. 4.

Boltzmann Equations for tree-level CPT -violating Leptogenesis
In our study of leptogenesis [1][2][3], we considered the Boltzmann equation for the number density n r of a fermion with helicity λ r = (−1) r −1 (r = 1, 2), in a homogenous and isotropic spatially flat FLRW space time [8]. The Boltzmann equation reads where f r (E, t) is the phase space density associated with n r (n r =ǧ 8π 3 d 3 k f r (E, t)), 2 H is the Hubble parameter and g is the determinant of the metric tensor; it is assumed [1][2][3] that B 0 min(T, m N ). On summing over the helicity λ r of the fermion [1][2][3], makes the second term on the left-hand side of (8) vanish. The term on the right-hand side of (8) is the collision integral C[ f r ]. In general for a species χ the collision integral describes the process 2ǧ denotes the total number of internal degrees of freedom and should not confused with the metric.
In curved space time, the collision integral is proportional to the square of the modulus of the amplitude of the scattering operator M for the decay processes relevant to leptogenesis: The delta function in (9) ensures conservation of the fourmomenta k . . N , the number of scattered particles at the interaction point, for both incoming and outgoing particles. In curved space time, we have used the covariant momentum integration element where √ −g δ (4) . For the spatially flat FLRW metric (6), we have that √ −g ∼ a 3 (t), and so where k [7] is "physical" spatial momentum, 3 Energy does not change under the redefinition of k to include the scale factor a(t) of the expanding universe. As standard, the scattering amplitudes for the appropriate interaction processes can be expressed, in terms of creationâ † i and annihilation a i operators of the respective quantum fields participating in the processes [20]: with T denoting time-ordered product; a (generic) quantum field operator φ(x) can be expanded in terms of the functions f i (x) which are solutions to the classical equations of motion for the (free) field φ(x) in curved space time: In a curved space time with metric g μν (x), the inner product between two functions f i (x) and g j (x) is defined as [20] ( The normalised solutions f i satisfy From this, it becomes evident that the functions f i in curved space time will be proportional to a normalisation factor that depends on the square root of the covariant volume V ∝ √ −g(x) at a given space time point x: On account of (14), (15) and (16), one obtains the relationŝ which implies that the creation and annhiliation operators are independent of √ −g. Hence, on account of (18), such volume normalisation factors will cancel out in the expression for the squared amplitude for the heavy-neutrino decay processes (13). However there remain space-time curvature corrections in the scattering amplitudes per se, as a result of modifications of the polarisation tensor and spinors entering such amplitudes, and, in loop cases, due to the curved space-time modifications of the dispersion relations of the fields circulating in the loops.
In our scenario of CPTV-induced leptogenesis [1][2][3], due to the non trivial background B 0 = 0, the dominant amplitudes of relevance to our discussion are the ones describing tree level decays of a right-handed neutrino N to standard model Higgs (h = h ± , h 0 )) and lepton ( = ( ± , ν L ) fields (all to be considered massless at the high temperatures of interest). In a plane-wave (i.e. Minkowski spacetime) approximation, a generic amplitude has the structure [1] iM Minkowski where the factor (1 ± γ 5 )/2 depends on the particular products of the decay; Y is the Yukawa coupling that appears in the so-called Higgs portal interaction of the model and connects the right-handed neutrino sector to the Standard Model sector; p i , i = , N are the relevant field momenta; u s ( p), u r ( p ) are the Dirac polarisation spinors with helicities λ r,s = (−1) (r,s)−1 , s, r = 1, 2. (The (Higgs) scalar polarisation is 1, independent of the space-time metric). There are restrictions in the various decay channels, as discussed in detail in [1][2][3]. These details will not be relevant for our discussion here, as we shall only restrict our attention to the potential effects on the amplitude of the slowly varying time dependence of a(t) and the KR field, through the relevant modifications of the spinor polarisation and the modified energy-momentum dispersion relations. This t-dependence implies that the spinors have also an explicit t-dependence in addition to the four-momentum dependence: u s ( p, t) and v r ( p , t) are solutions of the free Dirac equation in a spatially flat FLRW and KR backgrounds (4) [3].
In Appendix C we will discuss in detail, an n-th-order expansion in powers of H [9] for the spinor solutions of the Dirac equation in our time-dependent backgrounds. The respective spinor polarisation (of a given helicity λ) assumes the form (in the standard helicity basis ξ λ ): where ϕ n is a phase with corrections up to and including order n [9], and u ↑, ↓(0) λ (E (0) , k a(t) ) has the form of the corresponding polarization spinor in Minkowski space time, but with the spatial momenta being replaced by the "physical" momenta (12), while the energy E (0) is given by the Minkowski-form of the dispersion relation, but with the replacement (12) and contribution form the KR field(C11). The total energy E receives corrections from the expansion of the universe and the time-dependence of the KR field. (As shown in Appendix C, the phase ϕ nth coincides with the total energy E to this order. In our case, such phase factors are not relevant, since we are interested only in the collision terms (9) of the Boltzmann equation (8), which involve the square of the modulus of the scattering amplitudes and so phase contributions cancel out.) We note that in (20) the presence of the volume factors V ∼ √ −g ∼ a 3 (t). However, as we shall discuss in this article, it is important to note that the quantities which appear in the scattering amplitudes should have the volume factors removed. This will be linked with the hermiticity of the proper form of the Dirac Hamiltonian in time-dependent space-time geometries [22] and will result in the elimination of any potential dependence of the scattering amplitudes from such factors, although the space-time curvature-dependent corrections will remain. 4 The above corrections are assumed adiabatic, as appropriate for a slowly-expanding universe, and a background B 0 (4), which also exhibits comparable mild cosmic-time dependence, as appropriate for the conditions of leptogenesis in the model of [3]. As we shall show in this work, such corrections are proportional to powers of the Hubble parameter and the background B 0 . For the conditions of leptogenesis described in [3], the dominant corrections are of order H , and turn out to be negligible for the relevant range of the model parameters (5). Therefore, upon integrating over the redefined spatial momenta (12), one obtains the same Boltzmann equations as in [1][2][3], proving that, for spatially flat Robertson-Walker Universes, the flat space-time formalism to solve the Boltzmann suffices to produce results that are both qualitatively and quantitatively correct.
Before closing this section, we would also like to remark that the scaling (4), is found in [3] by computing in a flat space-time background the thermal condensate of the axial current for the fermions, summed over helicities λ, and showing that such a condensate vanishes: The temperature-dependent background (4) emerges in that case as a consistent solution of the equations of motion of the KR field [3]. In fact our analysis in [3] also implies that the result (21) remains valid in our expanding universe case with curved metric (6), despite the presence of the scale factor in the "physical" momenta (12). 4 Our results differ somewhat from those given in Ref. [21], which were based on a detailed derivation of the Boltzmann equation from the Kadanoff-Baym formalism in the context of a scalar field theory. In [21] it was claimed that the only effect of the curved space-time appears on the left-hand-side of the Boltzmann equation, describing the dilution of the particle number density due to the Hubble expansion H . These authors assume that the collision terms in the Boltzmann equation for the scalar-field scattering amplitudes are the same as for the case of flat space time, upon redefining their momenta to the "physical ones" (12). For us this is not the case. There are adiabatic corrections to the spinor polarisations; when loop contributions to the leptogenesis scenario are considered, there will be corrections as well to field propagators (including, in our case, the Higgs-scalar propagator). Such corrections for scattering amplitudes have been demonstrated clearly in [20]. Using Riemann normal coordinates (RNC), the corrections were shown to be proportional to positive powers of the space-time curvature. In the case of a spatially-flat Robertson-Walker universe, such corrections are expected to be encoded in the higher-order terms of the adiabatic (WKB-like) expansion (20) discussed here and in [9]. The explicit connection between the two works, via appropriate coordinate transformations that link the RNC expansion to the adiabatic expansion is still lacking. Nonetheless, for our purposes in Ref. [1][2][3] all such corrections turn out to be negligible.
Since the scaling of B 0 is not affected, compared to the case studied in [3] this will yield the same value for B 0 today as the one determined in that work. To an excellent approximation (for the parameter range (5)) the entire phenomenology of the flat space-time analysis of our earlier work [1][2][3] carries over to the full curved space time case,.
We now proceed to evaluate the space-time curvature corrections to the spinors due to the expansion of the universe. Although the RHN in the model [1][2][3] are Majorana, nonetheless our analysis is valid for both Dirac and Majorana spinors. 5

Spinors in spatially-flat expanding universe space-times with axial Kalb-Ramond (KR) backgrounds
In our model of leptogenesis, particle interactions occur on a background of a string gravitational multiplet which consists of graviton, Kalb-Ramond and dilaton 6 fields. The graviton background will be that of flat FLRW cosmological spacetime and the Kalb-Ramond field varies inversely as a power of temperature (2). Since in our leptogenesis scenarios both type of spinors, Dirac and Majorana, are involved in general, we cover here both case. We commence our discussion with the Dirac case

Dirac spinors in FLRW and KR axial backgrounds
The spatially flat FLRW space-time is described by the lineelement (6). The Dirac equation reads (for notations and conventions see Appendix A): where the Dirac matrices are tangent space ones, γ 5 , γ 0 , γ j , j = 1, 2, 3, satisfying the Clifford algebra (A4), and we adopt the chiral representation (A3). In Appendix C we solve (22) using an adiabatic (WKBlike) perturbative method, appropriate for slowly varying a(t), and B 0 (t) which is of relevance to our leptogenesis scenario [3]. The method we shall follow is developed in [9]. The corrections can be expanded in appropriate powers of the Hubble parameter H ; it follows from the parameter range (5) of the leptogenesis model [3] that |B 0 | H and that we are in the high temperature regime T T D ∼ m N . As shown in Appendix C, (cf. (C3), (C48), (C49), up to and including second order terms in an expansion in powers of H , we find the Dirac spinor for a fermion of mass m (of mass m) and helicity λ to be: and where n = 3 for the model of [3]; we will restrict our attention to this case. The quantities h and we assume [1-3] a fixed sign for B 0 > 0, the energies (frequencies) and ω 0 are taken to be positive. The reader should notice that for m = 0, one passes from (24) to (25), upon flipping the sign of m, m → −m and changing ↑ to ↓, and vice versa, where appropriate. Moreover, the expanding universe corrections vanish for massless fermions m → 0, as is the case of the SM leptons in the decay channels (19). Hence such spinors remain unaffected by the inclusion of curvature effects, apart from the overall factor a −3 (t) which appears as a result of their normalisation (C4).
The adiabatic corrections in (24), (25), will enter the expression for the modulus squared of the scattering amplitudes (19) that appears in the interaction terms in the Boltzman equations for leptogenesis in the scenario of [3]. The phase factors in these expressions are irrelevant as they cancel out in the Boltzmann collision term (9). The zero-th order term in the expansion coincides formally with the planewave solutions discussed in [1], provided one uses physical momenta (12). 7 As we shall demonstrate below, for the range of parameters (5), the curvature corrections in (24), (25), that take proper account of the Universe expansion, are negligible. Hence, the plane-wave approximation used in [1][2][3] to calculate the lepton asymmetry is fully justified in this case.

Extension to Majorana-Fermion case
Although the RHN in (19) is a right-handed field N R , with a Majorana mass M term, the results remain the same as in the Dirac case, apart from a relative normalisation factor of 1 2 in the kinetic terms of Majorana spinors in the Lagrangian. Indeed, if N R is the right-handed Neutrino spinor, then the Majorana mass term in the Lagrangian can be written as where N C R is the Dirac-charge-conjugate field, and N denotes the corresponding Majorana field defined as On the other hand, the kinetic term is also expressed (up to total derivative terms) in terms of the Majorana field N as Compared to the corresponding term in the Dirac case there is a factor of a 1 2 . (Majorana spinors, unlike Dirac fermions, do not couple to gauge fields, as they cannot be charged. They couple but only to gravity and so only ∇ μ the gravitational covariant derivative appears in their kinetic term.) The coupling of N R to the axial KR background now takes the form In addition, the model of [1,3] involves the Higgs-portal interactions which give rise to the decays (19). The Higgs field is viewed as an excitation from the standard vacuum, since in the leptogenesis scenario of [3] we are in the unbroken electroweak symmetry breaking phase. From (27), (29), (30), we therefore obtain the analogue of (22) for Majorana N spinors (28) in the model of [1][2][3]: where the axial background is of the form (2), Thus, apart from the relative factors of 1 2 , the analysis of the Majorana case would proceed in the same way as the Dirac case (24), (25), and will not be repeated here. (Such factors can be absorbed into the definition of the axial background field.)

Eastimates of curvature effects and connection with the plane-wave approximation for leptogenesis
We will now estimate the order of magnitude of the leading correction, proportional to H in (24) (or, equivalently, (25)).
For the leptogenesis scenario of [3], we have ( (5)): m = m N 10 5 GeV, and T m N T D B 0 . Also, during the radiation era of the universe, we have a(t) rad ∼ 1/T , and the Hubble parameter where M Pl ∼ 2.4 × 10 18 GeV is the reduced Planck mass, andǧ is the number of effective degrees of freedom of the system under consideration. For Standard Model like theorieš g ∼ 100, while for supersymmetric extensions this number is larger, but a natural range is which we assume for our purposes here (and in [1][2][3]).
The decays (19) preserve the helicity [1]. As follows from (26), for massless fermions such as the SM leptons in these decays, the zeroth order solution vanishes for one of the helicities [1][2][3], e.g.: For massive spinors, on the other hand, the leading O(H ) effects are easily estimated from from (24), (25). However, in view of the integration over momenta k ≡ k/a(t) in the collision term of the Boltzmann equation, we shall treat k as an integration variable, independent of a(t), and discuss the order of both quantities: at various k regimes. The temperature T (and, hence, H ((32)) is kept fixed, assuming that the universe in the radiation era behaves as a black body, and we are interested in the RHN decoupling temperature region T ∼ T D ∼ m N for the regime of parameters of the model of [3], (5), (33). We have: for the regime (5), (33).
We will now remark on the dependence of the polarisation spinors (20) on a(t) 3/2 . Such volume factors, if present, would be inconsistent with the general properties of the scattering amplitudes (13), discussed in Sect. 2. Any dependence of the scattering amplitudes on such factors is absent, due to the fact that the creation and annihilation operators of states that define the scattering (S-matrix) amplitudes are defined through appropriate inner products for curved space time (18), (15).
In the case of our spinors, therefore, a state a † i |0 = |i entering the corresponding scattering amplitude (19) should correspond to a spinor polarisation (20) without the a −3/2 (t) factors. This would imply that (for the evaluation of the Smatrix) the appropriate spinor polarisation, in an expanding universe, should be (2) .
In our context, this can be justified on noting [22] that in the case of time-dependent space-time metrics there are some subtleties in demonstrating Hermiticity of the Hamiltonian associated with the Dirac equation in curved space-time. The naive expression for the Hamiltonian, obtained by rewriting the Dirac equation as a Schrödinger equation, is not Hermitian, as explained in Appendix B,. One needs to appropriately redefine the Hamiltonian, in order to have a Hermitian Hamiltonian operator (B14). As discussed in detail in [22], and reviewed briefly in Appendix B, due to diffeomorphism invariance in general relativity, there are no timeindependent state-basis vectors (in contrast to the case of nonrelativistic quantum mechanics). If one uses the appropriate time-dependent basis (B10), then the correct generally covariant, Schrödinger equation with Hermitian Hamiltonian emerges from the original Dirac equation; in the case of the FLRW universe with axial KR background, the Dirac equation assumes the form (B22), i.e.: is the solution of the original Dirac equation (22). We note that equation (41), apart from the a(t) factors in the spatial derivative parts, looks like a Minkowski-space-time Dirac equation (in a B 0 background). Its solution is the spinor (40), S−matrix , which is independent of the covariant volume factor a 3/2 . The spinor u λ (E, k, a(t)) (2) S−matrix is used in the S-matrix amplitude. It is natural for the unitary Smatrix operator S, to be related to a Hermitian Hamiltonian operator, via S ∼ exp(−i H t). Thus, the scattering amplitude of the collision term (9) in the Boltzmann equation (8), is independent of any volume factors √ −g, and so in the limit where the adiabatic corrections to the spinors (24), (25) are ignored, one obtains exactly the flat Minkowski space-time results of leptogenesis of [1][2][3].
The above results demonstrate, therefore, that the adiabatic effects of the expansion of the universe in the presence of KR torsion on the Boltzmann collision term are negligible compared to the zeroth-order terms for the regime of parameters (5), (33), for the leptogenesis model of [3]. Thus the plane-wave approximation for the estimation of the lepton number in [1][2][3] is a very good one.

Generation of baryon asymmetry through the CPT -violating leptogenesis
In our earlier works [1][2][3] we simply stated that baryogenesis can proceed through Baryon (B)-minus-lepton (L)-number (B-L)-conserving sphaleron processes in the SM sector of the theory, following the seminal works of [12]. Sphaleron processes may lead directly to Electroweak Baryogenesis which, in its original form, however is not currently considered to be phenomenologically viable. In the spirit of the pioneering contribution of Ref. [13] we combine these processes with our Beyond-the-Standard-Model (BSM) leptogenesis mechanism, so as to obtain a baryon asymmetry through leptogenesis. In our context there are some subtleties and non-trivial mathematical features, due to the presence of the Kalb-Ramond background field B 0 in sphaleron processes. For the viability of our scenario for baryogengesis, we will need to show that the implications for the baryon sector remains unaltered from our previous work [1][2][3]. It will be instructive to first review briefly the electroweak baryoge-nesis mechanisms, and then the baryogenesis through leptogenesis approach. We will emphasise those features that will be essential for our approach.

Review of basic features of electroweak baryogenesis: sphalerons and triangle anomalies
Triangle anomalies lie behind the nonconservation of B and L numbers at a quantum level in the field theory of the SM. In Minkowski space time, for chiral (left-handed) fermion currents, pertaining to quarks and leptons, one has the anomaly equations where the corresponding currents J B(L) μ are defined over chiral (left-handed ( )) fermions, either quarks (B) or leptons (L) repsectively; J μ = species ψ γ μ ψ , where the sum is over the appropriate set of species of fermion. For our purposes here, this compact notation suffices. We do not give the detailed form of the currents. N f is the number of fermion families/generations (f). L f , denotes the lepton number for each family, with the total lepton number being defined as the sum L = f L f . We will restrict ourselves to SM where N f = 3; f = e, μ, τ for leptons; F a μν is the field strength of the weak SU(2) L gauge bosons, with a = 1, 2, 3 the SU(2) adjoint-representation index; g is the weak SU(2) L coupling; the hypercharge (Y) U(1) Y has anomalous gauge field contributions which are Abelian but are similar in form to the weak SU(2) L contribution and have not been given explicitly. The standard notation F a μν = 1 2 μνρσ F a ρσ denotes the dual tensor with μνρσ the (totally antisymmetric) contravariant Levi-Civita tensor.
Since the combinations F a μν F a μν = ∂ μ K μ are total derivatives, the integral is an integer, and a topological winding number. For perturbative gauge field configurations N = 0, but there are nonperturbative configurations for which this number is nonzero, and such configurations for the SM theory are instantons, and sphalerons [10,11]; the latter are unstable saddle-point (local maxima) solutions of the electroweak theory, for which the potential exhibits a periodic form, with a height separating the minima (at zero) of order m W /g 2 , where m W is the electroweak scale. This is the barrier that has to be overcome for B+L violation to occur. At zero temperatures, the instantons lead to tunneling through the periodic vacua, which leads to a very strong suppression of the baryon and lepton (B+L) number violation. For high temperatures, however, of relevance to the early Universe, the unstable sphaleron configurations can climb up the potential barrier ("thermal jump" on the saddle point), leading to relatively unsuppressed sphaleronmediated (B+L)-violating processes. By integrating the equations (42) over three space, and defining the corresponding charges of d 3 x J B(L) 0 as particle-antiparticle asymmetries: in the B(L) numbers, 8 we obtain the important relations: which imply the following conservation laws, that are respected by the sphaleron processes in the SM: The notation refers to particle-antiparticle asymmetry. In short-hand notation, since the antiparticles carry B and L numbers of opposite sign but equal in magnitude with the particle, the conservation laws (45) are expressed as the set of the following quantities being conserved by the (B+L)-violating sphaleron processes during the electroweak baryogenesis in the SM sector [12]. For our purposes here we concentrate on the B − L conservation law, (46). Adding the two equations (42), and using (44) and the B − L conservation (46), we readily obtain where B + L ≡ N F is the total fermion number in the SM sector.
From the detailed strudies of [12], we know that the rate where τ is the rate of the anomalous sphaleron-mediated processes for temperatures T, in the range where the sphaleron proicesses are active [12]: ∼ 10 12 GeV T T ew ∼ 100GeV, and T ew denotes the temperature of the electroweak phase transition. The detailed computation of [12] indicated that τ −1 = CT, where C is a function depending on the coupling constants of the SM. The temperature dependence of C can be inferred from the detailed studies of the anomalous fermion-number nonconservation of [12] but C has not been calculated analytically. Due to the nonperturbative gauge dynamics, C can be calculated using lattice gauge theories. Fortunately, we will not need the precise form of τ −1 (T).
From (49), we infer where t ini denotes some initial time within the temperature range that the sphaleron processes are active and in thermal equilibrium. Integrating over the time t (48) and using (50), we readily obtain for the Baryon asymmetry at time t: where we took into account that for the range of temperatures for which the sphaleron processes are active, the second (exponential) term on the right-hand-side of the first equality in (51) is heavily suppressed due to the large absolute value of the exponent. The above result was based only on the anomaly equation and the generic relation (49) but not on any detailed thermal behaviour of the sphaleron processes. In AppendixD we discuss a more physical way [12] of deriving (51), which makes use of the thermal equilibrium properties of the system in the range of temperatures where sphaleron processes are active. However, as we shall see, the two separate derivations of the baryon antisymmetry agree in order of magnitude. When the more detailed thermal properties are considered the form of the relation(51) remains unchanged, but the proportionality coefficient between B and B − L changes from 1/2 in (51) to 28/79 0.354 .
It should be noted that the above result is not affected by an extension to curved space-times, present in the early universe, since the triangle gauge anomaly (42), on which it is based is topological and as such is independent of the metric. For generic space times in addition to the gauge terms in (42), there are also gravitational anomaly terms, proportional to R μνρσ R μνρσ , where (. . . ) again denotes the corresponding dual in curved space time. For a FLRW universe, however, the latter terms vanish.
The temperature T D ∼ m N ∼ 100 TeV in the scenario of [1][2][3], is well within the range of active sphaleron processes in the SM. If T D is identified with a freeze-out time t F , then we can take t ini = t F . In the scenario of [1][2][3], (B(t ini ) = 0, and hence, at the sphaleron-freezout time t sph , which is later than t F , (t sph > t F ), the sphaleron-induced baryon asymmetry is of the same order as the lepton asymmetry generated at t F : as asserted in [1][2][3]. The numerical factor q ∼ O(1) (cf. (1)) has been estimated in [1][2][3] and remains approximately unchanged in the case of a slowly varying KR background where T 0 is the CMB temperature in the current-epoch). The reader should notice the opposite signs between lepton and baryon asymmetries, but this is not of concern, given that such a relative minus sign can be absorbed in the definitions of the baryon and lepton current in (42). The conventions are such that matter dominates antimatter in both baryon and lepton sectors. A similar relative sign difference between baryon and lepton asymmetries also appears in the approach of [13] and is standard in scenarios of baryogenesis through leptogenesis.

Independence of the anomaly equation from the KR background: two arguments
We shall check if the axial anomaly (42) is affected by the presence of our CPT -Violating KR background in two ways. The first uses an explicit calculation of the triangle graph and the second uses a topological argument . Both methods show that the KR background does not affect the generic result (42), and thus the mechanism of baryogenesis through leptogenesis survives. The arguments used are instructive and nontrivial and so are worth discussing.
• I. Diagrammatic argument We will follow the standard procedure and evaluate the one-loop triangle graph between two vector and one axial-vector vertices (see Fig. 1). In the presence of a constant KR background B μ = B 0 δ 0 μ the fermion propagator S F for the internal lines of the graph is where we have used the standard notation / A = γ μ A μ . Matter fermions, in the triangle anomaly calculation, can be considered to be massless at high temperatures. For Fig. 1 Generic triangle anomaly diagrams, with one axial vector (γ μ γ 5 ) and two vector (γ α,β ) vertices. The wavy lines indicate external Abelian (or non-Abelian) gauge bosons the case of the U(1) chiral anomaly 9 where J A(V ) is the axial (vector) fermion current, the · in the exponent of the exponential denotes the inner product between two four-vectors, and the Fourier-space quantity i μαβ ( p, q) is determined by applying the appropriate Feynman rules (for the U(1) gauge theory): The last terms in the parenthesis on the right-hand side of above indicates the Bose symmetry of the graph i μαβ (p, q) = i μβα (q, p).
The anomaly equation is obtained by evaluating the quantity by contracting it with the polarisatrion tensors for the external gauge bosons, and by passing into configuration space time. The external gauge bosons satisfy the on-shell conditions since they are massless (at temperatures above the electroweak phase transition). Gauge invariance requires: p α i μαβ (p, q) = 0, and q β i μαβ (p, q) = 0. (59) 9 Extension to the non-Abelian triangle anomaly, of interest for (42), is straightforward.
For the high-temperature regime of interest, the momenta | p| ∼ T, and hence such propagators can be expanded in powers of the weak background B 0 T. Hence, where the · · · denote higher powers of This expansion in terms of γ 0 i B 0 γ 5 i / p is actually a general way of using the diagrammatic analysis to prove that the contribution from the (constant) B 0 background to the anomaly vanishes: one may consider switching on the torsion B 0 background adiabatically, starting from an infinitesimal value.
To first order in the expansion in γ 0 i B 0 γ 5 i / p , a straightforward computation of the graphs of Fig. 1 can be performed, using the following identity for the trace of a product of n-even Dirac matrices Tr γ 1 γ 2 . . . γ n = Tr 1 2 {γ 1 , γ 2 · · · γ n } = n k=2 (−1) k g 1 k T r γ 2 · · · ((γ k )) · · · γ n , where g αβ is the metric tensor, and the notation . . . ((γ k )) . . . indicates that this particular Dirac matrix is absent from the respective product. Using some straightforward manipulations for the momentum integrals over k, we find that we need to evaluate the trace (61) for n = 6. This yields the following structure for the B 0 -dependent part of the anomaly where with Taking into account the symmetry of the graph under α ↔ β, and the conditions (59) for (on-shell) gauge invariance, it can then be seen immediately from (62), It should be also remarked that a generic nonconstant B 0 -torsion, also yield zero contributions to the triangle anomaly. This follows from the topological argument given below. Within the diagrammatic approach the method of using the expanded propagators (60) leading to (62), does not apply. One has to treat the field B 0 on the same footing as the background photon field used for the computation of the triangle anomaly. It can be shown that the B 0 contributions to the anomaly vanish on account of the Bianchi identity for the field strength (A13) of the background antisymmetric tensor KR field B μν : ∂ [μ H νρσ ] = 0 (with [. . . ] denoting total antisymmetrisation of the indices). • II. Topological argument: There is a topological method for understanding anomalies which is in terms of the Atiyah-Singer index theorem [24]. On a 4dimensional closed Euclidean manifold X with flat met-ric, the index theorem requires that where n ± denotes the number of ± chiral zero modes of the Dirac operator. This framework can be generalised to a curved manifold and applied to our case on noting that the KR-background-dependent terms in an effective low energy string action, can be interpreted in terms of generalised curvature and gravitational covariant derivative terms with torsion ("KR H-torsion") [1][2][3].
The pertinent Atiyah-Singer index theorem, associated with the zero modes of the generalised Dirac operator corresponding to a space-time manifold (M 4 ) with contorted spin-connection (ω + 1 2 H ), is given by : (omitting, for brevity, the gauge terms (. . . ), see Appendix E); as for the case of the flat manifold, the index theorem is related to the triangle anomalies appearing in (42) in the path-integral method of Fujikawa [25].
Explicit computations [17,18] show that (42) is independent of the KR H-torsion. One naively finds KR, Htorsion contributions to the integrand of the expression of the index (66), which, however, conspire to yield total derivatives and thus do not contribute [18]. This cancellation has its roots in the renormalisation-group properties of the low-energy field theory stemming from the underlying microscopic string theory. Indeed, at the level of the effective action, such H-torsion terms, and hence the potential B 0 (T )-background contributions to the baryonasymmetry rates, are renormalisation-scheme dependent; consequently these contributions, can always be removed by a judicious choice of renormalisation-group counterterms between the gauge and metric sectors of the theory [17]. More details are given in Appendix E.
This concludes our demonstration of the non contribution of the KR background to the triangle anomaly, and thus to the rates for baryon asymmetry during the electroweak baryogenesis, based on it. In AppendixD we present yet another derivation of this result based on thermal-equilibrium aspects of sphaleron processes.

Conclusions
In this paper, we have given a careful treatment of the Boltzmann equation used in the CPT violating tree-level leptogenesis scenario of [1][2][3] in the presence of time dependence from the expansion of the universe and the Kalb-Ramond background field. We have explained quite rigorously why the flat space-time analysis of the collision term leads to accurate results.
Following [9] we have explained why the zeroth-order WKB-expanded (plane-wave) solutions to the equation (22) (equivalent to a flat space-time analysis, as far as the collision terms in the Boltzmann equation are concerned) suffice to produce qualitatively and quantitatively correct results for the lepton asymmetry. In the specific parameter range (5) of [1][2][3], which is phenomenologically relevant, all the space-time curvature effects that characterise the higher-order WKB corrections are negligible. It must be stressed though, that for generic parameters, such curvature effects might lead to physically relevant corrections in the pertinent Boltzmann equations.
As an interesting by-product of our analysis of the WKBplane-wave solutions of the Dirac equation over space times with time-dependent metrics, we have related aspects of the solution to a properly defined Schrödinger equation with a Hermitian Hamiltonian (for a particular inner product [22]) associated with the Dirac equation.
Finally, we have explained in some detail how the lepton asymmetry generated by the CPT violating decays of heavy right-handed neutrinos in the scenarios of [1][2][3], can be transmitted to the baryon sector by means of sphaleron processes in the standard model sector of the theory. Some interesting properties of the KR background, namely its noncontribution to the anomaly equations relevant for lepton and baryon number violation, have been highlighted in that discussion.
The results presented here go beyond the particular example of the leptogenesis model of [1][2][3] and are nontrivial. They pertain to properties of the Dirac equation in curved space-time and KR backgrounds and attempt to examine in detail the influence of these backgrounds on the Boltzmann collision terms. Only a few studies pay attention to these important issues [20,21] and are not complete.We therefore hope that, in view of the above results for our particular model for leptogenesis [1][2][3], the discussion in this article will also make a useful contribution to the literature on quantum field theories in curved space times and the corresponding Boltzmann equations and their generalisations.

Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors' comment: This is a theoretical, formal work, not involving use of any kind of data.] Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/. Funded by SCOAP 3 .

Appendix A: Dirac equation in curved space times with time-dependent metrics: notation and some mathematical properties
In curved four-dimensional space-time with metric g μν (x) ≡ g μν (t, x), whose signature is (+, −, −, −) and μ, ν = 0, . . . 3, the motion of a free spin- 1 2 fermion of mass m, is determined by the Dirac equation [26]: (A1) (We use the notation tr for matrix transposition, ψ(x) for a spinor, and x for spatial coordinate vectors and the summation convention for a repeated index.) D μ is the spinorial covariant derivative andγ μ is a curved-space-time Dirac 4×4 matrix. Theγ μ satisfy the Dirac algebra where σ i , i = 1, 2, 3 are the Hermitian 2 × 2 Pauli matrices and I 2×2 is the unit matrix. The Dirac matrices γ a satisfy the Clifford algebra In terms of the spin connection 10 ω ab μ , where [ , ] denotes a commutator; the Latin indices are raised or lowered by the Minkowski metric η ab . The spin connection is is related to the vielbeins and the Christoffel symbol μ αβ via: The quantities μ in (A5), the Fock-Ivanenko coefficients [22], can be expressed as: On using the identity with abcd the totally antisymmetric Levi-Civita symbol 11 , γ 5 = iγ 0 . . . γ 3 , and {γ 5 , γ d } = 0, in the space of spinors, the Dirac operatorγ μ D μ is where the vector (A a ) and axial vector (B a ) potentials are given by and The vector potential A a may lead to non-Hermitian terms, which need to be interpreted through a modified inner product so as to preserve the Hermiticity of the Hamiltonian operator 10 ω abμ = g βα e a α ∇ μ e b β and ∇ μ is the gravitational covariant derivative. 11 0123 = +1 and the other components of abcd are determined by antisymmetry. appearing in the correct "Schrödinger equation" which the Dirac equation (A5) is mapped to [22].
In our work we will consider Lagrangians rather than Hamiltonians; the A α vector potential drops out of the Lagrangian: At this point we remark that in the leptogenesis scenario of [1][2][3] the quantity B μ is associated with an axial background stemming from the KR antisymmetric tensor B μν = −B νμ , with field strength with the symbol [. . . ] denoting complete antisymmetrisation of the respective indices. The three-form acts as torsion in the effective gravitational field theory [14][15][16]. On account of (A6), and for backgrounds with only a temporal component non trivial, B 0 = 0, it follows from (A12) the following Dirac equation: where we expressed the equation in tangent space, with γ a , γ 5 the tangent space Dirac and chirality matrices, respectively. 12 For a flat FLRW space-time we have the relations: In Appendix C we shall solve this equation so as to determine the effects of the expansion of the Universe on the spinor solutions of (A14). Before embarking on that task, however, we remark that the equivalence of the absence of A a in (A12) and its presence in the Hamiltonian formalism is discussed in detail in Ref. [22], and will be reviewed in Appendix B below. 12 Notice that the chirality matrix γ 5 in spatially flat Robertson-Walker space equals its flat Minkowski counterpart, since it is defined by where ε μνρσ = √ −g μνρσ , with μνρσ the flat Levi-Civita totally antisymmetric symbol, 0123 = +1, etc., and we used √ −g = a 3 (t) and (A16).

Appendix B: Mapping the Dirac equation to a "Schrödinger" equation with a Hermitian Hamiltonian
By writing the Dirac equation (A1) as a Schrödinger equation for the wave function of the fermion ψ(x) = ψ(t, x): we (naively) identify as the Hamiltonian. As noted in [22], this Hamiltonian is not hermitian for time dependent metrics g μν (t, x), in the usual inner product 13 between two wave functions φ 1 ( x) and φ 2 ( x) defined as [22] (where, in standard bra-ket notation, it is assumed that φ i ( x) = x|φ i for a bra x| and a ket |φ i , i = 1, 2) since However, in general relativity, with space and time dependent metrics, the complete basis to define φ(x) is not the time independent | x but rather |t, x . Consequently we write the Dirac fermion wave function ψ(x)(= ψ(t, x)) as The completeness relation leads to a modified inner product for wave functions: In view of the time dependence of the basis vectors one has:

x|ψ(t, x)
= i ∂ ∂t t, x| |ψ + i t, x| ∂ ∂t ψ 13 Due to the signature of our metric, we use the inner product which differs in sign from the definition given in [22].
where in the last equality we have used (B1).
To evaluate the first term in the left hand side of (B7), we take the time derivative of (B6). Hence, on using the notation Observing that the term in the second line is obtained from the first line by simply taking the Hermitean conjugate, we conclude that the equation (B8) is satisfied if the coefficients of the |t, x and t, x| vanish independently, 14 which leads to the relations [22]: From inspection, it can be seen that a solution [22] of (B9) is where the basis | x is time independent. Consequently, in view of (B7) and (B9), for timedependent metric backgrounds (in general relativity), it was postulated [22] that the correct quantum-mechanical Schrödinger equation is By making use of the completeness relation (B6), we may express the right-hand side of (B11) as: 14 It is useful to note that γ 0γ 0 (x) where the matrix elements of the operator H are defined so as to satisfy locality in space time: 15 t, x| H|t, y = − H δ (3) with the correct Hermitian Hamiltonian operator [22] in spinor space; its Hermiticity follows from the fact that the first term on the right-hand side of (B14) cancels the non-Hermitian part on the right-hand side of (B4), leading to for matrix elements of H on Dirac-spinor wave functions. From (B14), we note that the first term on the right-hand side, contains non-Hermitian parts for time-dependent metrics, which are cancelled by the the corresponding non-Hermitian parts of of the covariant derivative D i term, corresponding to vector potentials A i (cf. (A9), (A10)). Notice that in view of (B10), we have that i ∂ ∂ t t, x|ψ = i t, x| ∂ ∂t ψ and thus the non-Hermitian operator H in (B2) is not the proper Hamiltonian of the system. By writing the left-hand side of (B11) as: and using (B9), (B12), (B13) and (B14), we readily observe from (B11) that, if one accepts (B14), then the wave function t, x|ψ satisfies the equation: where H is given by (B2). This equation is identical to Eq. (B1), and, upon multiplication byγ 0 , to the original Dirac equation (A1). From (B10) and (B5) we observe that the solution ψ(t, x) of the Schrödinger equation (B11), with the hermitian Hamtiltonian (B14), is formally related to the solution 15 Notice the factor − γ 0γ 0 (t, y) −1 on the right-hand side of (B13), whose presence is a necessary consequence of the form of the completeness relation (B6) corresponding to the inner product (B3). ψ original (t, x) of (B1), the naive Schrödinger equation (with the non-hermitian "Hamiltonian" (B2)) by [22]: where we have used (B10).
For the case of a spatially-flat Robertson-Walker cosmological space-time, of interest to us, (B17) becomes with ψ herm (x) satisfying the equation: Were it not for the factor a(t) −1 this would be the Dirac equation in Minkowski space-time. However, the effect of the expansion of a spatially flat universe on the dynamics of spinors is encoded in that factor. The factor a −3/2 (t) in (B18) is related to the standard normalization factor 1/ √ V of a quantum field in a covariant volume V ∝ √ −g(x) for FLRW. The hermitian "Schrödinger" Hamiltonian (B14) corresponding to (B19) reads [22]: where we used (A16) and for the Fock-Ivanenko coefficients. The extension of the above results to the case with a nontrivial KR axial background with a non-trivial temporal component B 0 = 0, as is the case in [1][2][3], is straightforward. In that case the analogue of (B19) is It is possible to calculate modifications of plane-wave solutions of (B19), (B22), in a systematic adiabatic approximation following [9]. We do this in Appendix (C) for (B22) (or, equivalently for (A14)), of interest to us here, in order to determine the effects of the expansion of the Universe on the collision term of the Boltzmann equation (8), (19), used in the leptogenesis scenario of [3].

Appendix C: Dirac spinors in an expanding universe and a Kalb-Ramond background
In this Appendix we discuss the corrections to the form of the Dirac spinors induced by an expanding Universe, up to quadratic order in a perturbative adiabatic expansion in the Hubble parameter H . This is only required for massive fermions, since the massless case can be solved easily. Our analysis follows that of [9].
We shall be concerned with solutions of the Dirac equation in the presence of an axial constant background B 0 , given in (A14), which we give here again for convenience of the reader (we use tangent-space γ 0 , γ j , j = 1, 2, 3, Dirac matrices): Our representation of the Dirac game matrices is the chiral one (A3). We use the following notation for the helicity basis spinors: ξ r − → p , σ i p i p ξ r = λ r ξ r , λ r = ±1, p = | p|. In terms of creation and annihilation operators, spinors, and antispinors, the corresponding quantum field ψ(t, x) reads: where the Dirac polarisation spinor in the above helicity basis in a FRW universe with scale factor a (t) and Hubble parameter H =ȧ a are given by: with denoting the operation of complex conjugation and C the charge conjugation operator. The spinors u k,λ satisfy the orthonormality condition in the expanding universe: with similar relations for the antispinors v k,λ . The various terms in (C1) are evaluated with the spinor ansatz (C3): Putting these terms together (including the mass term) gives These two equations can be written compactly as: and The quantities α and β are both real. We should note that the machinery, that we will develop, is not needed for the case If this holds at some t = t 0 , it will hold at all t owing to unitary evolution. We will derive the adiabatic method to quadratic order in the Hubble parameter, which suffices for our purposes within the frameworkl of the leptogenesis scenarios of [1][2][3].
In this work we consider models for the time-dependence of a (t) and B 0 (t) in the radiation era, as in [3]: (where t 0 is present time) and with the integer n ≥ 3 and a 0 and b 0 positive. In the leptogenesis scenario of [3], we have n = 3 (cf. (4)), but in this Appendix we keep the general scaling power n. Thus, we have and where H =ȧ (t) a(t) = 1 2t is the Hubble parameter during the radiation era we are interested in. For an expanding universė a > 0. Thus, a perturbative expansion of the spinors in powers of ∂ t α is equivalent to an expansion in powers of H 1. Notice that the perturbative expansion measure the deviation of the scale factor of the Universe from constancy, and as such is independent of whether the mass of the fermions is zero or not. Thus, the expansion can equally apply to the Standard Model leptons, which are approximately massless at high temperatures in our leptogenesis scenarios, and the massive right-handed neutrinos.
For clarity, we shall keep explicit the λ dependence in the expressions below. .
We start the sequence of approximations: This leads to In contrast to F −1 , F 0,λ is complex but Hermitian but somewhat similar in structure. We follow similar steps to the previous steps otherwise. We diagonalise F 0,λ . Let where It is easy to check that .
We define and so The structure of F 1,λ resembles F −1 and so we can proceed as before.
We are going to do the iteration one more time.
and then and We have emphasised the dependence on λ in this formalism since it plays an important role in our theory. However, in order to lessen the burden on our notation, the dependence on λ will no longer be indicated; any λ dependence can be found in earlier formulae. Note that U 2 is a real symmetric matrix.The spinor solution is obtained with the help of We have explicit expressions for U 0 , U 1 and U 2 except for ξ 1 which we need to determine. Since it is clear from (C21) that and we need to determine sgn −ω 1 ∂ t ω 0 +ω 0 ∂ t ω 1 2ω 1 ω 2 1 −ω 2 0 on using (C12), (C13) and (C14). We have from (C17) So, from (C25) and (C24), Since Hence To make further progress we use (C12), (C13) and (C14).
In the expression for ξ 1 , B 0 will be ignored since it is very small in comparison to the other terms. Hence For the early Universe regime we are interested in, in the scenario of [1][2][3], we have that sgn 2 k a 2 − m 2 = +1, hence we set from now on We shall now summarise the key formulae for our analysis: where we keep the positive sign when taking the square root, due to the positivity of the energy ω 1 (assuming ω 0 > 0 to lowest order). Since in our perturbative expansions in this work we shall not consider terms higher than H 2 , we may truncate the above expression to (cf. (C13)) We also have from (C30) and We solve the Dirac equation using an adiabatic procedure, which assumes that the time derivatives of α satisfy 1 for j = 1, 2, . . .. As a bookkeeping device (for the order of adiabaticity) we will introduce the parameter in front of ∂ t . In the context of this notation we can say that Although, it will be seen to be, a posteriori, negligible for our application to leptogenesis, we shall retain expressions up to second order in the Hubble parameter H 2 (see (C13) and (C14)) 16 . From C33 we can deduce that approximations and In view of (C13) and (C14), we observe that the ratio ω 2 ω 1 differs from 1 by terms of order H 4 , which we ignore in our analysis.
Let us return to In this approximation and so the phase in ϕ 2 of (20) can be identified with ω 2 (on suppressing the λ dependence). Forβ = −1, we have (from now on we denote t ω 2 = ω 2 for brevity): (C36) 16 In our expressions corrections to O 4 will be computed, and then subsequently truncated to O 2 . and Forβ = 1, we have: In the case of a SM sector with N f fermionic generators and N s Higgs doublets, 17 considered in the second reference in [12], one has two conservation laws, for which one introduces chemical potentials: the B−L conservation, corresponding to a chemical potential μ and the hypercharge U(1) Y comnservation, corresponding to a chemical potential μ Y . Thus for the "I-th" particle we introduce the chemical potential with the chemical potential for the corresponding antiparticle being μ I = −μ I . From (D1), we obtain for the asymmetry of Higgs-like particles and for fermions of all N f generations n H ∼ N s μ H T 2 3 , The requirement of "neutrality" of the plasma of particles under the U(1) Y hypercharge, is expressed as [12]: I Y I n I = 0. (D4) In our case with two chemical potentials, this relation allows expression of, say, the (B − L) chemical potential μ in terms of the hypercharge chemical potential μ Y : The baryon number asymmetry is then determined by computing the quantity (n B(B) denote number densities of baryons (antibaryons): 3 n le f t−handed−quarks + n right−handed−quarks .
(D6) 17 As we consider here only SM fields (and their antiparticles), which participate in the sphaleron processes, there is no right-handed neutrino. In our case the RHN has decayed long before the sphaleron processes freeze out, and baryon asymmetry is generated. Moreover, as we are in a regime above the electroweak symmetry breaking, the Higgs particle (antiparticle) spectra contain both charged h ± and neutral Higgs h 0 .
Using then (D3), and the conventional quantum-number assignments of for the SM particles/antiparticles, we easily arrive at [12]: in the high temperature regime of interest,for which μ I T . The lepton asymmetry at this temperature L = n le f t−handed−lepton−doublets + n right−handed−lepton−singlets where, as we mntioned prefiously, we do not consider the heavy right-handed neutrino, which had already decoupled at the temperatures of the creation of the baryon asymmetry we are interested in (i.e. the freeze out of the sphaleron processes, which is slightly above the electroweak symmetry breaking temperature, of O(100)GeV).
We also have which allows μ Y to be expressed in terms of B − L, implying that the baryon asymmetry (D7) can be finally given as [12]: which is to be evaluated at the temperature at which the sphaleron processes decouple, which is of the order of the electroweak transition, slightly above it. The relation (D10) needs to be compared with (51). To this end, the reader should first recall the conservation by the sphaleron processes of (B−L), which implies that the latter quantity can be replaced in (D10) by (B − L)(t ini ) at some initial time value, in the scenario of [1][2][3] this can bne taken as the freezeout point of the heavy right-handed neutrino decays. In this respect, (D9) should be understood as being valid for a fixed equilibrium temperature, which does not change with time. Second, for a SM we have N f = 3, N s = 1, which implies B 28 79 (B − L) for temperatures above the electroweak phase transition. We thus observe that in order of magnitude our (51) was in excellent agreement with the more detailed derivation above.
Before closing we would like to discuss the rôle of the KR background B 0 (T ). Its presence could in principle modify the previous derivation leading to (D10), since the fermion dispersion relations get modified with λ the helicity of the spinor (i.e. the projection of the third component of the spin to the direction of the spatial momentum). When considering the distribution functions that enter the expressions for the particle (antiparticle) equilibrium number densities, taking into account that at the high temperature regime we are interested in T B 0 , one may expand the integrand in powers of B 0 /T 1. Thus, to first order in this small quantity, the reader can readily verify that the effects of the KR backround to the thermal equilibrium n F in (D3) is to add to their hand side terms of the form [1][2][3] n KR,λ with n F (B 0 = 0) given by (D3), g F denotes the number of degrees of freedom of the fermion in question, and the constant coefficient c 1 −0.16 has been computed in [2] (cf. Eq. (157) of that work; the equation refers to leptons, but it can be straightforwardly generalised to quarks). The dependence on the helicity of the KR-background correction term is due to the dispersion relation (D11).
The B 0 -dependence of (D12) disappears though once we average over quark helicities λ = ±1, as becomes necessary in order to evaluate the Baryon asymmetry using (D6): given that λ λ = 0. This is in agreement with our arguments in Sect. 4.2 on the non-contribution of the KR background to the anomalies, which determine the rate of the baryon asymmetry B. As we have seen above, the only dependence of B on B 0 (T ) comes through the chemical potential term (B − L) which is conserved, and thus can be evaluated at the leptogenesis heavy-right-handed-neutrinodecay freezeout point in the scenario of [1][2][3], cf. Eq. (52).
with torsion where R = R(ω) = d ∧ ω + ω ∧ ω is the curvature twoform and D(ω) the covariant derivative one-form, both with respect to the torsion-free (Riemannian) connection and K is a generic contorsion tensor (in our case in (E1)) K ≡ 3 2 H ), we observe that the torsion-dependent parts integrate to zero, being exact forms.
The proof can be extended to include a generic higher order term of order k in the perturbative expansion of the integrand of (E1) in powers of the gerneralised curvature form R, where the symbol P denotes form ordering, given the "antiderivation of degree 1" character of the exterior derivative d, when acting on wedge products of forms (d (a ∧ b) = d a ∧ b + (−1) p a ∧ b, where a is a p-form). Since, as we have seen above (E13) d dt P t i dt = d X , and the functionals F k (P t i ) contain generalised curvature two-forms only, it follows immediately, on account of the Bianchi identity (E11), and Stokes theorem, that the volume integral for a manifold without boundary, as we assume to be the case for the FLRW Universe of interest to us here. Thus the index of the generalised Dirac operator coincides with that in a Riemannian manifold. In particular, this means that by considering the Dirac operator in the presence of external gauge potentials in flat Minkowski space times, the resulting triangle anomalies [25] are independent of the torsion, and so independent of our KR background field B 0 , in the case of [1][2][3].