Cutting rules on a cylinder: a bottom-up approach to quantum kinetic theory for the early universe

Nonequilibrium quantum field theory is often used to derive an approximation for the evolution of number densities and asymmetries in astroparticle models when a more precise treatment of quantum thermal effects is required. This work presents an alternative framework using the zero-temperature quantum field theory, $S$-matrix unitarity, and classical Boltzmann equation as starting points leading to a set of rules for calculations of thermal corrections to reaction rates. Statistical factors due to on-shell intermediate states are obtained from the cuts of forward diagrams with multiple spectator lines. It turns out that it is equivalent to cutting closed diagrams on a cylindrical surface.

Nonequilibrium quantum field theory is often used to derive an approximation for the evolution of number densities and asymmetries in astroparticle models when a more precise treatment of quantum thermal effects is required. This work presents an alternative framework using the zerotemperature quantum field theory, S-matrix unitarity, and classical Boltzmann equation as starting points leading to a set of rules for calculations of thermal corrections to reaction rates. Statistical factors due to on-shell intermediate states are obtained from the cuts of forward diagrams with multiple spectator lines. It turns out that it is equivalent to cutting closed diagrams on a cylindrical surface.

I. INTRODUCTION.
Exploring particle physics models' implications in cosmology, we often focus on their predictions for dark matter relic density or CP asymmetries in the standard model sector. Despite sophisticated modern methods being developed, at low temperatures, the classical Boltzmann equation is often used for the sake of simplicity. The particles' interactions are computed in terms of the zero-temperature quantum theory, while the kinetic description of the thermal medium is entirely classical.
Improving the calculation's accuracy demands the inclusion of thermal corrections. Though this may be straightforward at leading order in the coupling constants, higher-order calculations necessary for matter asymmetry generation require more advanced treatment to be consistent with no asymmetry produced in thermal equilibrium [1]. The usual (top-down) procedure to derive the quantum-corrected Boltzmann equation is based on the in-in (closed time path) formalism of nonequilibrium quantum field theory [2][3][4][5][6][7]. Then, the Keldysh-Schwinger contour is stretched from −∞ to ∞ and back, allowing for the asymptotic particle states to be used approximatively. We shall assume the approximation is reasonable for systems of our interest.
This work presents an alternative approach to finitetemperature effects based on an in-out formalism, the Smatrix unitarity, and CP T symmetry. We show that constructing forward diagrams with each propagator wound on a cylindrical surface any number of times is a helpful heuristic tool [8] and correct statistical factors for the on-shell intermediate states are obtained. Moreover, the CP T and unitarity relations for the equilibrium reaction rate asymmetries are respected, as they are in zerotemperature calculations [9][10][11][12][13]. These features of our general framework will be demonstrated within the seesaw type-I leptogenesis, where the results are easily comparable to the broad literature on the subject [14][15][16][17][18][19][20][21][22][23][24] In this work, the primary focus is to simplify the inclusion of thermal effects at higher orders. However, to introduce a new diagrammatic concept, we first consider a simple leading-order example of the N i → lH decay.
In the seesaw type-I model Lagrangian, right-handed neutrino singlets N i couple to the standard model lepton and Higgs doublets via Yukawa interaction where i and α label the right-handed neutrino and lepton family, respectively. Here, the only massive particles are the right-handed neutrinos with Majorana masses denoted M i , while the other particles are massless. How do these interactions affect the N i number densities in the expanding universe? At the O(Y 2 ) order, the only contribution expected within the classical approach comes from the tree-level decays entering the Boltzmann equation aṡ where H denotes the Hubble constant while thermally averaged decay rates are conventionally written in the form ofγ with abbreviated notation and [dp] = d 3 p/((2π) 3 2p 0 ). The small circle in Eq. (3) indicates the reaction rate computed in terms of the classical Maxwell-Boltzmann phase-space densities. We assume kinetic equilibrium for all particle species, including the right-handed neutrinos (see Refs. [25][26][27] (Fig. 1a -1c) and the cylindrical diagram corresponding to Fig. 1c (Fig.  1c ). Diagrams in Figs. 1a, 1b can also be obtained in the representation of Refs. [28][29][30][31][32][33]. The black and blue dots in Figs. 1c and 1c do not represent vertices and only serve for the reader's reference. more detailed discussion). When the temperature drops below their masses, a departure from chemical equilibrium parametrized by µ Ni occurs. The inclusion of general out-of-equilibrium distributions will be discussed in section IV. In Eq. (3), the zero-temperature amplitude squared, summed over spin, isospin, and lepton flavors, reads We would like to point out that even when describing particles as classical point-like objects revealing their quantum nature only through the interactions when they come close to each other, the right-hand side of Eq. (2) is incomplete. Let us continue systematically, considering the quantum description of interactions as a black box that can be asked for a complete list of processes occurring at specific order in the couplings. Each reaction contributing to Eq. (2) should be represented by a square of an amplitude, where certain diagrams interfere. As a consequence of S-matrix unitarity and optical theorem, each such contribution can be related to cuts of forward diagrams. Examples of those appearing at the second order in the Yukawa coupling, and containing one N i in the initial state, are shown in Fig. 1. In Fig. 1a, the contribution leading to the simple N i decay rate in Eq. (3) is shown. Next, in Fig. 1b, the disconnected diagrams involve standalone Higgs lines acting as Lorentz invariant momentum conservation deltafunctions, (2π) 3 2p 0 δ (3) (p − p ), leading to the same amplitude squared as in Fig. 1a though the presence of additional Higgs particle in the initial state brings one extraf H to the decay rate that is now written as Once the 'Pandora's box' of disconnected diagrams has been opened, one can equally well include any number of Higgses, leptons, and right-handed neutrinos, all as initial states from the black box point of view. However, the disconnected amplitudes should be taken with care. If, for example, the standalone Higgs lines in Fig. 1b were joined to each other and not connected to the rest of the diagram, the contribution to the decay rate would be equal toγ Ni→lH from Eq. (3) times the overall number of Higgs particles in the universe and thus already accounted for in Fig. 1a. To avoid such singular and spurious contributions, we should only consider those in which the spectator lines joined together as a single line are connected to the interaction part of the diagram. In other words, we can start drawing forward diagrams such as in Fig. 1a on a cylindrical surface joining the righthanded neutrino lines. Next, we wind the propagators around as shown in Fig. 1c . In the end, the contributions of these diagrams will only differ by the powers of phase-space densities and overall sign when the sum of fermionic winding numbers is odd. The sum of all corresponding reaction rates is denoted γ Ni→lH and equals with newly introduced Starting the sum with w = 0, factors 1 − f eq l and 1 + f eq H come out on the right-hand side of Eq. (7). The 'eq' superscript denotes the equilibrium Fermi-Dirac or Bose-Einstein distribution, respectively. We emphasize that there is no direct manifestation of Pauli's exclusion principle in this picture as we started from the kinetic theory of classical particles. Otherwise, the initial states with identical right-handed neutrinos leading to Eq. (7) could not be considered, and quantum statistics would not emerge. Instead, the alternating sign in Eq. (8) comes from the crossings of legs in Feynman diagrams. This is similar to the t-and u-channel relative sign in 2 → 2 scattering with identical fermions in the final state.
At this place, the appearance of quantum statistics in Eq. (7) may seem a miraculous coincidence. It will be shown in section IV that it is not so and this appearance is general. Now we continue with thermal corrections to the reaction rate asymmetries, where the inclusion of the proper statistical factors is more subtle.
III. CP ASYMMETRIC THERMAL RATES FOR LEPTOGENESIS.
In leptogenesis, the collision terms in the Boltzmann equation for the lepton number generation, to the first order in asymmetric quantities, can be divided into two parts. The wash-out terms are proportional to the net number density, i.e. differences between particle and antiparticle densities, and quickly erase any accidentally produced asymmetry if a thermal equilibrium is maintained. The source terms, on the other hand, generate the asymmetry when the departure from equilibrium occurs (at temperatures lower than the right-handed neutrino mass) and are expressed in terms of the symmetric parts of the densities and asymmetries of reaction rates.
In this section, we focus on thermal corrections to the latter. To that purpose, let us briefly comment on the connection of the presented framework to the common accounting of thermal effects. We focus on the role of the scalarH propagator in the right-handed neutrino selfenergy. In Fig. 2a, its contribution takes a usual form of Whenever the N j self-energy occurs as a subdiagram, depending on the overall kinematics, the positive frequency on-shell part of this propagator may contribute or not. If it does, the diagrams in Figs. 2b and 2c do as well. However, each of them contributes by some number ofH particles to the initial state of the process and thus by certain power offH to the corresponding reaction rate. The sum of all such diagrams may be represented by a single self-energy diagram with the scalar propagator replaced by equivalent to the usual form of the positive frequency part of a thermal propagator. The negative frequencies contribute by the diagrams similar to those in Figs. 2b and 2c withH legs flipped to the opposite side of each diagram. In the N i decay rate, it is the on-shell part, the second term on the right-hand side of Eq. (9), that contributes to the CP violation at zero temperature. Therefore, to take the thermal corrections to the reaction rate asymmetries into account, we should consider all the windings of the propagators involved in the contributing diagrams.
As a next step, to keep the simplicity of diagrammatic notation, the standard approach to CP asymmetries relying on the Cutkosky rules and imaginary parts of loop integrals needs to be modified. Here we follow our previous work devoted to asymmetries at zero-temperature (see Ref. [34]). Following this procedure will allow us later (in Eqs. (19)) to treat the unitary CP -violating cuts in the same way (with the same statistical factors) as for the final states.
To sketch the idea briefly, let us expand the unitarity condition S † = S −1 for S = 1 + iT . Employing we obtain where the T -matrix elements for the i → f reaction can be written as (12) is always truncated at a finite order given by the power of the coupling constants. For the squared amplitude entering the reaction rate, we use that should be divided by V 4 = (2π) 4 δ (4) (0) to obtain (2π) 4 δ (4) (p f − p i )|M f i | 2 . Each term in Eq. (13) may be interpreted as a contribution of a forward diagram with one or more simultaneous on-shell cuts. With iT † from Eq. (12), the CP asymmetry ∆|T f i | 2 = |T f i | 2 − |T if | 2 can be expressed as [34] ∆|T

From Eq. (14), the CP T and unitarity relations
are manifest. Analogous relations hold for the equilibrium reaction rate asymmetries when the classical Boltzmann approach is used [9][10][11][12][13]. Concerning thermalcorrected rates, similar cancelations were studied in Ref. [35] in the scenario of asymmetric dark matter freeze-in. The meaning of Eq. (14) can be understood from Figs. 3a and 3b. In both diagrams, i = |N i , f = |l α H and n = |l βH . After subtracting similar contributions with 3. Lepton number violating contributions to the CP asymmetry in the Ni → lH decay (Figs. 3a and 3b) and the s-channel part of the real-intermediate-state-subtracted lH →lH scattering (Fig. 3c). The latter can be obtained by a cyclic permutation of the diagrams in Fig. 3a using the approach of Refs. [10,34].
the intermediate states in a reversed order (by means of Eq. (14)), the asymmetry of the squared amplitude for the N i → lH decay, before the intermediate state momentum integration is performed, reads Here the two terms correspond to the diagrams in Figs. 3a and 3b, respectively. Within the standard Cutkosky approach, the same contributions to the asymmetry are obtained from the imaginary part of the self-energy and vertex loop integrals. Integrating Eq. (16) over the l, H,l andH phase-space, and dividing by twice the symmetric part of the squared amplitude from Eq. (5), we obtain the usual expression for the lepton number violating part of the zero-temperature CP violating parameter [36] where Moreover, the asymmetries of thelH → N i inverse decay and lH →lH scattering (see Fig. 3c), summed over spin, isospin, and flavor, can be obtained by cyclic permutations of the diagrams in Figs. 3a and 3b, leading to the same ∆|M | 2 as in Eq. (16). It is the advantage of Eq. (14) that it allows treating the asymmetry in scattering, at this order known as the real-intermediatestate-subtracted scattering, in the same way as the asymmetry in the N i decay. In the standard approach, the tree-level lH →lH amplitude is considered with finite right-handed neutrino width [9]. Taking the square and subtracting the part quadratic in the width results in the same asymmetry in lH →lH as can be obtained directly using Eq. (14) and the diagram from Fig. 3c with its t-channel counterpart (see Ref. [34] for the details). This feature is crucial for the inclusion of thermal corrections to the reaction rate asymmetries entering the lepton number source term. Those are obtained by including windings of the l, H,l,H, and N i lines in Fig. 3, summing over all winding numbers, the powers of circled densities. Then for the asymmetry in N i → lH, thel,H intermediate state receives the same quantum statistical factors as were derived in section II for the l, H final states.
Here the phase-space integration is understood as in Eq. (4) with additional (2π) 4 δ (4) (p l + p H − pl − pH ) and thel,H three-momentum integrations. We further parametrize f Ni as f Ni = f eq Ni +δf Ni . Note that again, we use the circled notation in Eq. (19) to emphasize the calculation in terms of the zero-temperature Feynman rules. One of this paper's points is that no 'uncircled' amplitudes are needed to include quantum thermal corrections to the Boltzmann equation.
When the detailed balance condition in the lH ↔lH scattering is taken into account, one can see that the asymmetries in Eqs. (19b) and (19c) are equal. Moreover, replacinglH with lH and vice versa in Fig. 3 changes the order of the two terms in Eq. (14), leading to an extra minus sign for lH → N i asymmetry. Thus we obtain an example of CP T and unitarity relations ∆γ lH→Ni + ∆γ lH→lH = 0 (20) that is well known in the literature [9][10][11][12][13]. Here, however, crucial differences appear, and they belong to the main results of this paper. The reaction rate asymmetries are computed with statistical factors for the final and intermediate states due to the summation of the windings on the cylinder. Thus it is remarkable that employing the cylindrical diagrammatic representation of the previous section allows us to rephrase the thermal-corrected asymmetries in terms of the zero-temperature ones. Relations for thermal-corrected reaction rates have been considered in Ref. [35]. There the standard Cutkosky approach is applied to asymmetry freeze-in in two sectors at different temperatures. The point of the calculation is the extraction of the loop's imaginary part. Sticking to these thermal cutting rules becomes rather complicated when higher-order terms of the expansion in Eq. (14) are needed (see Ref. [34] for the zero-temperature case). Another advantage of Eq. (14) over the common procedure involving the loop's imaginary part extraction is that the unitarity relations, such as in Eq. (20), are always evident and simple to check. The lepton number asymmetry evolution arises from the difference of the Boltzmann equations for lepton and antilepton number densities. Considering the reactions mentioned above, the complete O(Y 4 ) thermal-corrected source term, earmarked at the beginning of this section, equals to Using the detailed balance conditions and Eq. (19), we can rewrite the source term in Eq. (21) as [14,15,37] ∆γ Interestingly enough, the last bracket in Eq. (22) is linear inl,H distributions. Refs. [14,15] pointed out that using a finite-temperature generalization of Cutkosky cutting rules [38,39] leads to the wrong statistical factors with incorrect additional −2f eq l f eq H term [40][41][42]. Instead, the causal n-point functions [37] or direct closed-time-path derivation should be used [14][15][16][17][18][19][20][21][22][23][24]. Remarkably, our result is in exact agreement with the latter, which seems to be the first time the correct form of Eq. (22) is derived using the Boltzmann equation and S-matrix elements [43].

IV. QUANTUM KINETIC THEORY AND GENERAL ONE-PARTICLE DISTRIBUTIONS.
Up until now, it might have seemed that the presented method works for equilibrium densities only. Here we show that it is suitable for nonequilibrium cases too, and that it is not an accident. In general, a quantum state of a many-particle system is given in terms of a density matrixρ. Its hermiticity and positive definiteness allows us to writê with a hermitian operatorF [44]. To simplify the description, a finite volume V 3 and time T 1 are considered. The one-particle states are labeled simply by p, q, . . . suppressing the information on particle species and their degrees of freedom. We proceed as if all particles were bosonic [45]. In the Boltzmann equation, we only deal with oneparticle densities. Therefore, we approximate operator F byF Tracing over the Fock space basis |i 1 , i 2 , . . . labeled by one-particle occupation numbers, we obtain [46] with Z p = exp{F p }/(exp{F p }−1). The mean occupation number of a one-particle state equals In analogy to the equilibrium case, it is natural to intro- In practice,f p is determined in terms of f p through Eq. (26). Usingρ = SρS † to describe the temporal evolution, the S-matrix unitarity leads to [47] where . . . stands for the terms that vanish in Tr[ρa † p a p ] forF in Eq. (24). From Eqs. (27) and (12) we obtain where {i}, {n} denote sequences of one-particle occupation numbers. Dividing Eq. (28) by T 1 V 3 gives a change of the particle density per unit time on the left-hand side, while the reaction rates counting the interactions per unit volume per unit time appear to the right. The T -matrices in Eq. (28) are computed employing Wick's contractions represented by Feynman diagrams. The trace allows to draw them as diagrams on a cylindrical surface with two or more on-shell cuts. The two cuts that are always there are marked by i and n. For each term in Eq. (28), the one-particle states other than p fall into two distinct groups.
In the first group, there are those that do not participate in interactions -their a's or a † 's are not contracted with fields in any of the T -matrices. If q is such a state, at the diagrammatic level, it can be represented by a standalone loop whose winding number is i q = n q . Therefore, in addition to each diagram containing vertices, there are many contributions from cuts of such standalone loops bringing in some powers of circled densities. Finally, they completely factor out by the contribution of the same states to Z.
The second group then contains one-particle states connected to a vertex at least in one diagram. Within the chain of k + 1 diagrams in Eq. (28), for each state labeled by r = p, three cases may occur: 1. the occupation number i r goes up by one unit in some T , then down in another, 2. the same in a reversed order, 3. i r is lowered and increased (the external a's contracted with a † r a r from the fields) in a single T . At the level of Wick's contractions, the first case leads to in a T -matrix with n r = 1 + i r . The opposite process, with a † r after a nr r replaced by a r in Eq. (29), brings in another factor of √ 1 + i r . Thus, in this case, we get ∞ ir=0 (1 + i r )(f r ) ir = Z r ∞ ir=0 (f r ) ir = Z r (1 + f r ) (30) where Z r is canceled by Z in Eq. (28). Similarly, in the second and third case is obtained. Note that we do not consider more contractions within the external states and fields. For V 3 → ∞, such terms can be neglected as they are restricted to a lower-dimensional submanifold of the phase-space. Finally, state p in Eq. (28) requires n p = i p such that one of the above three cases, with n p −i p = 1 (case 1.) or −1 (cases 2. and 3.), takes place in T ni . As in Eqs. (30) and (31), we obtain factors Z p (1 + f p ) or −Z p f p . These correspond to the gain or loss terms, respectively, in the Boltzmann equation, while the presence of Z p ensures that 1/Z in Eq. (28) is completely canceled. Moreover, the sums in Eqs. (30) and (31) imply that occupation numbers higher than one may be represented by windings of the external lines in exact agreement with what we saw in sections II and III.

V. CONCLUSIONS.
A diagrammatic concept connecting the classical Boltzmann equation and quantum kinetic theory has been introduced in this work. Statistical factors due to on-shell intermediate states have been formally represented by cuts of forward diagrams with multiple spectator lines. The resulting reaction rates and their asymmetries have been expressed in terms of zero-temperature amplitudes, while the statistical factors arose from explicit diagrammatic summation over propagator windings on a cylinder. We showed how this relation is connected to the quantum formulation of kinetic theory and that it is not an accident. Furthermore, for the first time the unitarity and CP T relations between thermalcorrected reaction rate asymmetries have been formulated analogously to the zero-temperature case -see Eq. (20) for the simplest example -and less trivial examples will be considered in our follow-up studies.