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

In a model for leptogenesis based on spontaneous breaking of Lorentz and CPT symmetry [1-3], we examine the consistency of using the approximation of plane-wave solutions for a free spin-(1/2) Dirac (or Majorana) fermion field propagating in a Friedmann-Lemaitre-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).

In a model for leptogenesis based on spontaneous breaking of Lorentz and CPT symmetry [1][2][3], we examine the consistency of using the approximation of plane-wave solutions for a free spin- 1 2 Dirac (or Majorana) fermion field propagating in a Friedmann-Lemaî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).

I. 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 string-inspired 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 right-handed 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 , 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 string-inspired 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 T D ∼ m N ∼ 10 5 GeV, B 0 (T D ) ∼ O(1 keV), (5) 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 d 2 dt 2 a (t) . 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 space-time 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 space-time 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 curvature-induced 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 plane-wave 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] in the SM sector. Sphaleron processes preserve the difference B − L between baryon (B) and lepton (L) numbers [11], This is the route to baryogenesis in the conventional leptogenesis scenario [12]. However we need to check that the presence of the KR background B 0 (known to play the rôle of totally antisymmetric torsion [13] in string theories) does not affect [14,15] the anomaly equations [16] for the baryon and lepton numbers needed in the route [11] to baryogengesis.
The structure of our article is the following: In section II 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 [17,18].
In section III, we obtain systematic corrections to plane-wave solutions of the Dirac equation (in subsection III A) and of the Majorana equation (in subsection III B) on a spatially-flat FLRW space time in the presence of the KR background (4). The results are similar in the two cases.
In subsection III C, for the parameter range (5), we demonstrate that any space-time curvature corrections to the flat space-time 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 [19].
In section IV 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 section V. 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 space-times, in the presence of axial backgrounds of relevance to the leptogenesis scenario of [1][2][3]. In Appendix B we show, following [19], 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 section IV by incorporating high temperature effects properly. Finally, in Appendix E, 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 section (IV).

II. 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 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 four-momenta k µ (i) ≡ k (i) , i = 1, . . . 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 2ǧ denotes the total number of internal degrees of freedom and should not confused with the metric.
where √ −g δ (4) ( i k (i) ) is the curved-space-time-momentum delta-function δ (4) g (k). 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 a † i and annihilation a i operators of the respective quantum fields participating in the processes [17]: 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 [17] ( 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 space-time) 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 (C.11). 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 ϕ n−th 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 amplitiudes 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 timedependent space-time geometries [19] 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 4 Our results differ somewhat from those given in ref. [18], 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 [18] 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 [17]. 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.
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).
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

III. 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 space-time 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

A. Dirac Spinors in FLRW and KR Axial Backgrounds
The spatially flat FLRW space-time is described by the line-element (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 (WKB-like) 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. (C.3), (C.48), (C.49), 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: 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 B0 > 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 (C.4).
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 plane-wave 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.

B. Extension to Majorana-Fermion Case
Although the RHN in (19) is a right-handed field NR, 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 NR 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 NR to the axial KR background now takes the form 7 In our case the phase factor exp(−i t ω 0,λ ) differs from exp(−iω 0,λ t), because of the a(t) dependence of the integrand. However, because these phase factors are irrelevant, as already mentioned, the zeroth order approximation will lead to the results for lepton asymmetry derived in [1][2][3].
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), Bµ = ∂µb = B0 δ0µ, with B0 > 0 given in (4). 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.) C. 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 = mN 10 5 GeV, and T mN TD B0. 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 theoriesǧ ∼ 100, while for supersymmetric extensions this number is larger, but a natural range is 10 2 ǧ 10 3 , 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 ∼ TD ∼ mN 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 section II. 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 S-matrix) the appropriate spinor polarisation, in an expanding universe, should be In our context, this can be justified on noting [19] 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 [19], and reviewed briefly in Appendix B, due to diffeomorphism invariance in general relativity, there are no time-independent 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.: in tangent space notation, where ψ original (x) ≡ ψ original (t, x) 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 B0 background). Its solution is the spinor (40), u λ (E, k, a(t)) 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 S-matrix 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.

IV. 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 [11]. 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. [12] 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 B0 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 baryogenesis mechanisms, and then the baryogenesis through leptogenesis approach. We will emphasise those features that will be essential for our approach.

A. Review of Basic Features of Electroweak Baryogenesis: Sphalerons & 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]; 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 mW/g 2 , where mW 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 sphaleron-mediated (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 particleantiparticle 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 [11]. 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 ≡ NF is the total fermion number in the SM sector.
From the detailed strudies of [11], 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 [11]: ∼ 10 12 GeV T Tew ∼ 100 GeV, and Tew denotes the temperature of the electroweak phase transition. The detailed computation of [11] 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 fermionnumber nonconservation of [11] 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 tini 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 Appendix D we discuss a more physical way [11] 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 TD ∼ mN ∼ 100 TeV in the scenario of [1][2][3], is well within the range of active sphaleron processes in the SM. If TD is identified with a freeze-out time tF , then we can take tini = tF. In the scenario of [1][2][3], ∆(B(tini) = 0, and hence, at the sphaleron-freezout time t sph , which is later than tF , (t sph > tF ), the sphaleron-induced baryon asymmetry is of the same order as the lepton asymmetry generated at tF : 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 B0(T) ∼ B(T0) ( T T 0 ) 3 background (where T0 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 [12] and is standard in scenarios of baryogenesis through leptogenesis.

B. 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µ = B0 δ 0 µ the fermion propagator SF 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 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 .
For the high-temperature regime of interest, the momenta | p| ∼ T, and hence such propagators can be expanded in powers of the weak background B0 T. Hence, where the . . . denote higher powers of γ 0 i B0γ 5 i / p . This expansion in terms of γ 0 i B0γ 5 i / p is actually a general way of using the diagrammatic analysis to prove that the contribution from the (constant) B0 background to the anomaly vanishes: one may consider switching on the torsion B0 background adiabatically, starting from an infinitesimal value. To first order in the expansion in γ 0 i B0γ 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 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 B0-dependent part of the anomaly with Y1(k, p, q) = (k − p) 2 k 0 (k + q) 2 + k 0 (k 2 + p · q + p · k) + q 0 (k 2 − k · p) + p 0 (k 2 + k · q) , 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), (63) and (64) that all Yi = 0, i = 1,8 Hence the B0-dependent terms do not contribute to the triangle anomaly. It should be also remarked that a generic nonconstant B0-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 B0 on the same footing as the background photon field used for the computation of the triangle anomaly. It can be shown that the B0 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 [20]. On a 4-dimensional closed Euclidean manifold X with flat metric, 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 [21]. Explicit computations [14,15] show that (42) is independent of the KR H-torsion. One naively finds KR, H-torsion contributions to the integrand of the expression of the index (66), which, however, conspire to yield total derivatives and thus do not contribute [15]. This cancellation has its roots in the renormalisation-group properties of the lowenergy 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 B0(T )-background contributions to the baryon-asymmetry 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 [14]. More details are given in Appendix E.
This concludes our demonstration of the noncontribution of the KR background to the triangle anomaly, and thus to the rates for baryon asymmetry during the electroweak baryogenesis, based on it. In Appendix D we present yet another derivation of this result based on thermal-equilibrium aspects of sphaleron processes.

V. 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 WKB-plane-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 [19]) associated with the Dirac equation.
Finally, we have explained in some detail how the lepton asymmetry generated by the CPT violating decays of heavy righthanded 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 [17,18] 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. 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 [22]: (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 I2×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 [19], 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 (Aa) and axial vector (Ba) potentials are given by and Ba = dbc a e µ d ω µ bc .
The vector potential Aa 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 appearing in the correct "Schrödinger equation" which the Dirac equation (A5) is mapped to [19]. 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 10 ω abµ = g βα ea α ∇µe b β and ∇µ is the gravitational covariant derivative. 11 0123 = +1 and the other components of abcd are determined by antisymmetry.
with the symbol [. . . ] denoting complete antisymmetrisation of the respective indices. The three-form acts as torsion in the effective gravitational field theory [13]. On account of (A6), and for backgrounds with only a temporal component non trivial, B0 = 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 Aa in (A12) and its presence in the Hamiltonian formalism is discussed in detail in ref. [19], and will be reviewed in Appendix B below.
we (naively) identify as the Hamiltonian. As noted in [19], 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 [19] (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: 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). 13 Due to the signature of our metric, we use the inner product which differs in sign from the definition given in [19].
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 ∂t = ∂ ∂ t , 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 [19]: From inspection, it can be seen that a solution [19] of (B9) is where the basis | x is time independent. Consequently, in view of (B7) and (B9), for time-dependent metric backgrounds (in general relativity), it was postulated [19] 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: where the matrix elements of the operator H are defined so as to satisfy locality in space time: 15 t, x| H|t, with the correct Hermitian Hamiltonian operator [19] 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 Di term, corresponding to vector potentials Ai (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: 14 It is useful to note that γ 0γ0 (x) . 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). 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 ψ original (t, x) of (B1), the naive Schrödinger equation (with the non-hermitian "Hamiltonian" (B2)) by [19]: 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) satrisfying 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 [19]: where we used (A16) and for the Fock-Ivanenko coefficients. The extension of the above results to the case with a non-trivial KR axial background with a non-trivial temporal component B0 = 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 B0, 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 (C.1) are evaluated with the spinor ansatz (C.3): Putting these terms together (including the mass term) gives These two equations can be written compactly as: and and The quantities α and β are both real. We should note that the machinery, that we will develop, is not needed for the case m = 0 since F λ −1 is diagonal. From the Dirac orthogonality condition If this holds at some t = t0 , 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 B0 (t) in the radiation era, as in [3]: (where t0 is present time) and with the integer n ≥ 3 and a0 and b0 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 universeȧ > 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: is a real symmetric matrix. We note that 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 F 2,λ = D 2,λ − iU † 2,λ ∂tU 2,λ , where 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 U2 is a real symmetric matrix.The spinor solution is obtained with the help of Equivalently h−1 = U0U1U2h2.
In the expression for ξ1, B0 will be ignored since it is very small in comparison to the other terms. Hence 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. (C.13)) ω1 ω0 1 + (∂tα) 2 m 2 8ω 6 0 − m 4 (∂tα) 4 128ω0 12 . (C.31) We also have from (C.30) and .

(C.33)
We solve the Dirac equation using an adiabatic procedure, which assumes that the time derivatives of α satisfy | 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 U0 ∼ O 0 , U1 ∼ O 2 and U2 ∼ O 4 . 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 (C.13) and (C.14)) 16 . From C.33 we can deduce that approximations In view of (C.13) and (C.14), 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 h−1 = U0U1U2h2.
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): Using these results, we find that the zeroth order ((. . . ) (0) ) terms in an expansion in powers of H for h ↑,↓ −1 (λ), obtained on setting ω 0,λ = ω 1,λ (= ω 2,λ ), are with ω 0,λ > 0 given by (C.11) and α λ by (C.8), (and we assume [1-3] a fixed sign for B0 > 0. Up to the energy-dependent normalisation factors 1 √ 2 ω 0,λ , and irrelevant phase factors, this result coincides with the corresponding expressions for the spinors of [1], for helicity λ, which provides a self-consistency check of our approach. The energy ω0 satisfies the dispersion relation (C.11), which is the same as the dispersion relation of [1] upon the correspondence of the momentum with the physical momentum (12), k = k/a(t), here. At the next order, we obtain from (C.44) and (C.45): are given in (C.47), and we used (C.12), (C.13). The energies (frequencies) ω0 > 0 are taken to be positive. The reader should notice that one passes from (C.48) to (C.49) upon flipping the sign of m, m → −m, and changing ↑ to ↓, and vice versa., where appropriate.
The expanding Universe corrections (proportional to powers of the Hubble parameter) enter the spinor solutions (C.3), which in turn participate in the expression for the scattering amplitudes (19) that enter the interaction terms in the Boltzman equations for leptogenesis in the scenario of [1][2][3]. We estimate the order of such corrections for the range (5) of the parameters of the model of [3] in section III, and show that they are negligible, thus justifying the plane-wave approximation for leptogenesis used in that work.
Appendix D: Thermal-equilibrium treatment of electroweak baryogenesis, and connection with CPT -violating leptogenesis In this Appendix we review some thermal-equilibrium aspects of sphaleron processes, which are relevant for our discussion of baryogenesis in section IV. In our leptogenesis model [1][2][3] we have not needed chemical potentials for the generation of tree-level lepton asymmetries induced by the KR background B0. A detailed discussion of the communication of the lepton asymmetry to the baryon sector, requires us, however, to introduce chemical potentials, that would implement the pertinent conservation laws (47) in a path integral formulation of the effective action. Following the second reference of [11] (see section 11.2.1), in the regime of (high) temperatures, we may consider the number densities of leptons and quarks in the SM sector to be in thermal equilibrium, as is the case where the sphaleron processes are active. Assuming for simplicitly T mW, with mW the electroweak scale, and using the standard finite-temperature distribution function, we see that the difference ∆n between particles and antiparticles in the equlibrium number density of Bosons (Fermions), with a chemical potential µB(µF) T , behaves as In the case of a SM sector with N f fermionic generators and Ns Higgs doublets, 17 considered in the second reference in [11], 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 (D.1), we obtain for the asymmetry of Higgs-like particles and for fermions of all N f generations ∆nH ∼ Ns µH T 2 3 , ∆nF ∼ N f µF T 2 6 . (D. 3) The requirement of "neutrality" of the plasma of particles under the U(1)Y hypercharge, is expressed as [11]: 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: 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 (D.7) can be finally given as [11]: 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 (D.10) 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 (D.10) by ∆(B − L)(tini) 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, (D.9) 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, Ns = 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 B0(T). Its presence could in principle modify the previous derivation leading to (D.10), 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 B0, one may expand the integrand in powers of B0/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 ∆nF in (D.3) is to add to their hand side terms of the form [1][2][3] ∆n KR,λ F = ∆nF(B0 = 0) + c1 g F λ B0 T 2 , (D. 12) with ∆nF(B0 = 0) given by (D.3), g F denotes the number of degrees of freedom of the fermion in question, and the constant coefficient c1 −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 (D.11).
The B0-dependence of (D.12) disappears though once we average over quark helicities λ = ±1, as becomes necessary in order to evaluate the Baryon asymmetry using (D.6): given that λ λ = 0. This is in agreement with our arguments in section (IV B) 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 B0(T ) comes through the chemical potential term ∆(B−L) which is conserved, and thus can be evaluated at the leptogenesis heavy-right-handed-neutrino-decay freezeout point in the scenario of [1][2][3], cf. Eq. (52).
we obtain from (E.10): Thus, under a change in the torsion part of the generalised spin connection the relevant Pontryagin index changes by an exact form d X, with X(K, ω) transforming covariantly, as follows from the covariant transformation properties of both K and R under tangent-space rotations.
Since the homotopically extended Dirac genus (integrand of (E.1), after homotopic extension through (E.6)) is an infinite series of various polynomials in Pontryagin indices of various orders (cf. (E.2), (E.3)), which are linear combinations of terms of the form (E.8), it will also change under infinitesimal t-changes by the exterior derivative d of a covariantly transformed form. Indeed, consider a generic polynomial f(Pt i) of homotopically extended Pontryagin indices Pt i, of order n, f(Pt i) = n k=1 α k P k t i . Under the action of the t-homotopic exterior derivative dt, we have: 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 (E.13) d dt Pt i dt = d X , and the functionals F k (Pt i) contain generalised curvature two-forms only, it follows immediately, on account of the Bianchi identity (E.11), 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 [21] are independent of the torsion, and so independent of our KR background field B0, in the case of [1][2][3].