Arxiv:2012.14213V1 [Math.AP]

Home , Exponential stability

arXiv:2012.14213v1 [math.AP] 28 Dec 2020 in oeEnti itiuin em-ia distribut Fermi-Dirac distribution, Bose-Einstein tion, u n nry h rcs entosof definitions precise conser The the energy. implies and function tum Dirac-delta 4-dimensional The particles. r dnia,w denote we identical, are where (1.3) ttime at (1.2) as tvsi unu otmn qain n[1,Bs ta.poie se provided heavy-ion al. interactions, s et pion-nucleus Buss announced as [11], such [53] interaction In heavy-ion al. quantum equation. between et Boltzmann interaction Li quantum the ativistic about engineering. quant results and relativistic theoretical physics the and of in fields interest various increasing in been has there theory, where (1.1) 1.1. where iitcqatmBlzanequation: Boltzmann quantum tivistic e od n phrases. and words Key h yaiso igeseisfs-oigqatmprilsi des is particles quantum fast-moving single-species of dynamics The H EAIITCQATMBLZANEUTO NEAR EQUATION BOLTZMANN QUANTUM RELATIVISTIC THE oiain n re history. brief a and Motivations τ F σ C × t ( Abstract. otmn qainfrbt oosadfrin hnteini the when classic equilibrium. fermions unique global and a a bosons of both existence best for global equation the the Boltzmann to prove we equation on paper, Boltzmann theory this quantum mathematical relativistic no the been has for el there of importance, simulations such the a of and physics experiments in collision used quantum widely quant the been relativistic has the equation of Boltzmann development quantum recent the With single-spec of cles. dynamics the describes equation) Uhlenbeck ( 1frbsn,and bosons, for +1 = ( ,θ g, ,p t p, x, ∈ F F 1 [0 1 F , eoe h ieeta rs-eto hc ecie h collisio the describes which cross-section differential the denotes ) ( , p ∞ 2 sammnu itiuinfnto ntepaepit( point phase the on function distribution momentum a is ) ′ F , ) F .Terltvsi unu olso prtri ie by given is operator collision quantum relativistic The ). 2 3 h eaiitcqatmBlzaneuto o h relati the (or equation Boltzmann quantum relativistic The ( F , ICA A,JNWOJN,ADSO-A YUN SEOK-BAE AND JANG, WOO JIN BAE, GI-CHAN q W ′ 4 ( + )(1 = ) ( eaiitcqatmBlzaneuto,Rltvsi Ue Relativistic equation, Boltzmann quantum Relativistic ,q p, p C µ ( ∂ Z ,F ,F F, F, F, | τF µ p R F ′ 3 q , 3 τ dq q ( = ′ 0 p = ) = )1+ ))(1 Z p 0 R − ∂ 3 EQUILIBRIUM 1. as ) 2 c o emos hntefunctions the When fermions. for 1 t dp p F sσ ′ Introduction 0 τF ′ + ( C ,θ g, Z ytercn dacsi h eaiitcquantum relativistic the in advances recent the By p 4 R ( s ( F 3 ∇ · q ) 1 and dq δ .Tetasto rate transition The ). )) q o,Nniereeg methods. energy Nonlinear ion, (4) ′ x 0 ′ − F W ( σ p 1+ (1 = ( µ ( ,θ g, ,q p, + C ( τF r ie nScin1.2. Section in given are ) q | ,F ,F F, F, F, p µ e atmvn unu parti- quantum fast-moving ies ′ mmcais h relativistic the mechanics, um 1 q , − crn ngahn.I spite In graphene. in ectrons lslto oterelativistic the to solution al ( p ′ p ) fatos nweg.In knowledge. authors’ of ′ ildsrbto snearby is distribution tial ′ )1+ ))(1 h xsec fsolutions of existence the degneigsc sin as such engineering nd µ − n eto sn h rel- the using neutron and ecin,electron-nucleus reactions, ) q , mBlzanequation Boltzmann um ′ µ W τF ainlw fmomen- of laws vation ) , itcUehling- vistic ( ea xeiet of experiments veral 2 vrlexperimental everal ,q p, ( rbdb h rela- the by cribed ln-hebc equa- hling-Uhlenbeck q F ′ 1 )) | p ,p x, , F ′ F q , 3 2 ( ) ′ , p sbetween ns sdefined is ) ∈ F ) F 3 T 4 and ( 3 q × ) , R F 4 3 2 GI-CHANBAE,JINWOOJANG,ANDSEOK-BAEYUN collisions and neutrino-nucleus interactions using an equation called Giessen Boltzmann- Uehling-Uhlenbeck (GiBUU) transport model. Since the speed of the electron in graphene is close to that of light, Lapitski developed the lattice Boltzmann method to simulate the stream of the electrons in graphene in [51]. We also introduce some other experiments on the relativistic quantum Boltzmann model [49, 54, 78]. Despite such modern advances in physics and engineering, there have been quite few mathematical studies on the relativistic quantum Boltzmann equation. Akama [1] considered several types of relativistic quantum kinetic models such as the Boltzmann, the Fokker-Planck, the Landau equations, and this was the first appearance of the relativistic quantum Boltzmann equation. The unique determination of the relativistic quantum equilibrium satisfying H-theorem and conservation laws has been considered by Escobedo et al. in [25, 26]. To the best of authors’ knowledge, the mathematical theory on the existence of solutions to the equation has not been studied yet. Regarding the Newtonian non-quantum Boltzmann equation, a brief list of the classic literature includes [2, 4, 8, 9, 12–14, 17, 19, 20, 28, 29, 33, 34, 37–39, 41, 43, 50, 69, 70, 81]. The mathematical studies on the relativistic non-quantum Boltzmann equation are also not rich. In 1940, Lichnerowicz and Marrot then suggested the first relativistic Boltzmann equation in [56]. In [23], Dudyński and Ekiel-Je˙zewska studied the linearized relativistic Boltzmann equation and obtained the estimates for the linear term in the hard and soft potential cases. The authors also constructed the global existence of mild solution in [22]. Glassey and Strauss obtained some estimates on the derivatives of the collision map between pre-post collisional momenta including the upperbound on the average of ∂p′/∂p and ∂q′/∂p in [32]. Very recently, Chapman et al. [16] provided an analytic and numerical evidence that the lower bound for the Jacobian determinant for the collision map p p′ or q′ in the center- of-momentum frame is zero. In the nearby-equilibrium regime, Glassey7→ and Strauss [30] established a unique global solution in T3 for the hard potential case, and it was extended to the whole space case via the fourteen-moment compensating function in [31]. For the soft potential case, the global existence and the asymptotic behavior were constructed in torus in [71] and extended to the whole space in [77]. In [40], Guo and Strain obtained the stability of the relativistic Vlasov-Maxwell-Boltzmann system via the use of two different coordinate represenatations for the post-collisional momenta in order to resolve the issue of the momentum singularity. In [21], Duan and Yu provided the global existence of solutions in the weighted L∞ framework. Recently, Wang [82] showed the global wellposedness of the relativistic Boltzmann equation with large amplitude initial data. In the case of Coulombic interaction, Strain and Guo [75] constructed a unique global-in-time classical solution for the relativistic Landau-Maxwell system. The reduction of the relativistic collision operator using the center-of-momentum frame can be found in [26, 72]. Andréasson et al. [3] proved the finite-time blowup of the gain term when the initial data starts with the characteristic ball. The uniform L1 stability of mild solution is established in [42] when the initial data is sufficiently small and it decays exponentially fast. In [46], Jang and Yun proved the regularizing estimate of the gain term. For the spatially homogeneous case, Strain and Yun obtained useful inequalities about the relativistic pre-post collisional velocity and proved the existence of solution in [76]. The uniform L∞ bounds of the solution are established in [45] and the Lp bounds are established in [47]. For the Newtonian limit of the relativistic particles, we refer to [73]. The mathematical studies on the non-relativistic quantum Boltzmann equation are also not rich. We would like to mention that Nordheim in 1928 and Uehling - Uhlenbeck in 1933 discussed the non-relativistic quantum Boltzmann equation in [48] and [79], respectively. THE RELATIVISTIC QUANTUM BOLTZMANN EQUATION NEAR EQUILIBRIUM 3

Benedetto et al. [6,7] showed the rigorous validity of the quantum Boltzmann equation from the N-body Schr¨odinger equation in the weak coupling regime. In spatially homogeneous case for the Fermi-Dirac particles, Lu classiﬁed quantum equilibrium in [61]. In [67], Lu and Wennberg established the strong stability in L1 and proved the convergence of the solution to the equilibrium. The weak solution to the equation is constructed in the soft potential case [62, 63]. For the spatially homogeneous case, several studies on bosons can be found in [10, 27, 55, 58–60, 64–66, 68]. For general mathematical and physical reviews, we refer to [4, 13–15, 17, 19, 29, 37–39, 74, 81].

1.2. Notations. Before introducing our main results, we deﬁne several notations on the relativistic quantities. First of all, we remark that we normalize all physical constants to 1 throughout the paper including the speed of light c. We denote the energy-momentum 4-vector as pµ and usually write it as pµ = (p0,p1,p2,p3). The energy-momentum 4-vector ν with the lower index is written as a product in the Minkowski metric pµ = ηµν p , where the Minkowski metric is given by ηµν = diag( 1, 1, 1, 1). The inner product of energy- µ − momentum 4-vectors p and qµ is deﬁned via the Minkowski metric:

3 pµq = pµη qµ = p0q0 + piqi. µ µµ − i=1 X The energy of a relativistic particle is given by p0 = 1+ p 2. Then we can see that the | | µ inner product of an energy-momentum 4-vector with itself is p pµ = 1. We note that the p µ − µ inner product of energy-momentum 4-vectors is Lorentz invariant p qµ =Λp Λqµ, where Λ is a Lorentz transform which will be deﬁned below. We deﬁne the relative energy s between two energy-momentum 4-vectors pµ and qµ as (1.4) s(pµ, qµ)= (pµ + qµ)(p + q )= 2pµq +2, − µ µ − µ and the relative momentum g between them as

(1.5) g(pµ, qµ)= (pµ qµ)(p q )= 2( pµq 1). − µ − µ − µ − q q Note that s = g2 + 4. In the same sense, we deﬁne the relative momentum between pµ and µ µ µ µ µ µ another momentum p′ and q′ asg ¯ = g(p ,p′ ) andg ˜ = g(p , q′ ), respectively. Also, we µ µ µ µ 2 2 deﬁnes ¯ = s(p ,p′ ) ands ˜ = s(p , q′ ), so thats ¯ =g ¯ +4 ands ˜ =˜g + 4.

1.3. Conservation laws and the Boltzmann H-theorem. The conservation law of the pre-post collisional energy-momentum 4-vectors is given by µ µ µ µ (1.6) p + q = p′ + q′ . The inner product with itself gives µ µ µ µ (p + q )(pµ + qµ) = (p′ + q′ )(pµ′ + qµ′ ). µ µ µ µ µ µ Using p pµ = p′ pµ′ = q qµ = q′ qµ′ = 1, we have p qµ = p′ qµ′ . Similarly, we can have µ µ µ µ − p pµ′ = q qµ′ and p qµ′ = p′ qµ, which implies µ µ µ µ g = g(p , q )= g(p′ , q′ ), µ µ µ µ (1.7) g¯ = g(p ,p′ )= g(q , q′ ), µ µ µ µ g˜ = g(p , q′ )= g(p′ , q ). 4 GI-CHANBAE,JINWOOJANG,ANDSEOK-BAEYUN

Then it turns out [47, Proposition 2.7] that (1.6) can further imply the Pythagorean theorem as (1.8) g2 =¯g2 +˜g2. The angle of the scattering kernel σ(g,θ) in (1.3) is given by µ µ (p q )(p′ q′ ) (1.9) cos θ = − µ − µ . g2 By (1.7) and (1.8), cos θ can also be written in the following form [18, page 277, (A.11)]: g¯2 (1.10) cos θ =1 2 . − g2 Now we introduce the center-of-momentum framework. The relativistic pre-post collisional momenta p, q and p′, q′ satisfying the conservation law (1.6) can be expressed as p + q g (p + q) w p′ = + w (γ 1)(p + q) · , 2 2 − − p + q 2 (1.11) | | p + q g (p + q) w q′ = w (γ 1)(p + q) · , 2 − 2 − − p + q 2 | | where γ = (p0 + q0)/√s, and w denotes S2 component of the unit sphere: (1.12) w = (sin θ cos φ, sin θ sin φ, cos θ),

0 0 on θ [0, π] and φ [0, 2π]. The pre-post collisional energy p′ and q′ is written by ∈ ∈ 0 0 0 p + q g p′ = + (p + q) w, 2 2√s · 0 0 0 p + q g q′ = (p + q) w. 2 − 2√s · Once we divide each side of (1.1) by p0, then we can rewrite the relativistic quantum Boltzmann equation as follows:

∂tF +ˆp xF = Q(F,F,F,F ), (1.13) · ∇ F (x, p, 0) = F0(x, p). 1 Note that Q(F,F,F,F ) = p0 C(F ). We generally denote the normalized momentum asp ˆ where p p pˆ = = . p0 1+ p 2 | | In this framework, the relativistic quantum collisionp operator is reduced to the following form [72]:

µ µ (1.14) Q(F1, F2, F3, F4)= dq dw vø(p , q ) σ(g,θ) R3 S2 Z Z F (p′)F (q′)(1 + τF (p))(1 + τF (q)) (1 + τF (p′))(1 + τF (q′))F (p)F (q) , × 1 2 3 4 − 1 2 3 4 where the Møller velocity vø is given by p q 2 p q 2 g√s v (pµ, qµ)= = . ø p0 − q0 − p0 × q0 p0q0 s

THE RELATIVISTIC QUANTUM BOLTZMANN EQUATION NEAR EQUILIBRIUM 5

One of the properties that the relativistic quantum collision operator satisﬁes is the following identity 1 dp Q(F,F,F,F ) µ =0, R3 p Z which implies the conservation laws of the total mass, momentum and energy: d 1 (1.15) dx dp F µ =0. dt T3 R3 p Z Z The H-theorem for the relativistic quantum Boltzmann equation is established in [26] as

d 1 dx dp F ln F τ − (1 + τF ) ln(1 + τF ) 0. dt T3 R3 − ≤ Z Z 1.4. Global equilibria. It has been shown in [25] that the global equilibria for the relativistic quantum Boltzmann equation (1.1) has the form of 1 (p)= , for ν(p)= ap0 + b p + c, F eν(p) τ · − for some constant a, b, and c. In this paper, we rescale the problem and consider the situation that the macroscopic mean velocity b = 0 is zero. Then we can define the global equilibrium m(p) in the form of 1 (1.16) m(p)= , eap0+c τ − where a> 0 and c a in the case of bosons and a> 0 and c R in the case of fermions. In this paper, we only≥− consider a> 0 and c> a in the case of∈ bosons in order to exclude a possible blow-up of the equilibrium. We choose− and fix the constant a and c and consider the initial distribution F0 nearby the equilibrium m(p). Then we will prove that the particle distribution F whose initial distribution F0 is sufficiently close to the relativistic quantum N 2 global equilibrium m(p) will converge to m(p) in Hx Lv sense for N 3. We also denote the non-quantum relativistic equilibrium J(p) as ≥ 0 ap0 (1.17) J(p )= e− . 1.5. Spaces. Now we define some notations on norms and inner products which are frequently used throughout this paper. The constant C is generically used, whose value can be changed from line to line. Especially, when we want to indicate the dependency of a, we 2 indicate it as Ca. We define the standard L norm as 1 1 2 2 2 2 2 2 f Lp = dp f(p) , f Lx,p = dx dp f(x, p) , k k R3 | | k k T3 R3 | | Z Z Z and we define the weighted L2 norm as 1 1 2 2 2 2 f ν = dp ν(p) f(p) , f x,ν = dx dp ν(p) f(x, p) . k k R3 | | k k T3 R3 | | Z Z Z The standard L2 inner product is given by

2 2 f,g Lp = dp f(p)g(p), f,g Lx,p = dx dp f(x, p)g(x, p). h i R3 h i T3 R3 Z Z Z We use the multi-index notation

α = (α0, α1, α2, α3), 6 GI-CHANBAE,JINWOOJANG,ANDSEOK-BAEYUN to simplify the diﬀerential operator:

α α0 α1 α2 α3 ∂ = ∂t ∂x1 ∂x2 ∂x3 . For a brevity of discussion on the frequently used function m + τm2, we often write the variable only at the end of the function as (m + τm2)(p)= m(p)+ τ(m(p))2. We define the higher-order energy norm as follows: t 1 α 2 α 2 (f(t)) = ∂ f(t) 2 + ∂ f(s) ds. E 2 k kLx,p k kx,ν α N Z0 α N | X|≤ | X|≤ 1.6. Hypothesis on the collision kernel. We assume that the differential cross section σ(g,θ) satisfies the hard potential assumption with an angular cut-off as in [23, 24]: (1.18) σ(g,θ)= g sin θ. This assumption is analogous to the standard hard-sphere assumption for the Newtonian Boltzmann equation; if p + q or p q is sufficiently small or if p q is much larger than p0 q0 , then it behaves| as the| Newtonian| − | hard sphere kernel. | − | | − | 1.7. Main results. In this paper, we prove the existence of a unique global-in-time classical solution of the relativistic quantum Boltzmann equation nearby the global equilibrium. Before we state our main theorem, we would like to introduce the reformulation of the relativistic quantum Boltzmann equation via a special linearization that is relevant to the relativistic and quantum case. In the non-quantum case, the standard decomposition of the perturbed solution f near p 2/2 p0 the global equilibrium µ is F = µ + √µf for µ = e−| | or µ = e− in the Newtonian and the relativistic cases, respectively. But in the quantum case, the previous decomposition F = µ + √µf does not guarantee the non-negativity of Lf,f 2 as in Lemma 2.2. Thus h iLv inspired by the previous work related to the quantum kinetic models in [5, 52, 57, 80], we choose the following decomposition of F : F (x,p,t)= m(p)+ m(p)+ τm2(p)f(x,p,t), where the global equilibrium m(p) is definedp as in (1.16): 1 m(p)= . eap0+c τ − We plug the decomposition into (1.13) and divide each side of the equation by √m + τm2 to have

∂tf +p ˆ xf + Lf = Γ(f)+ T (f), (1.19) · ∇ f(x, p, 0) = f0(x, p). The linear term Lf is given by Lf = ν(p)f + K f K f, 1 − 2 where ν(p) is the collision frequency of a relativistic quantum particle, and K1 and K2 are compact operators. The right-hand side of (1.19) consists of nonlinear terms Γ(f) and T (f) where Γ(f) consists of all the second-order nonlinear terms and T (f) consists of all the third- order nonlinear terms. We can easily check that the fourth-order nonlinear terms disappear by cancellation. In the linearization process, we observe that the collision operator does not THE RELATIVISTIC QUANTUM BOLTZMANN EQUATION NEAR EQUILIBRIUM 7 satisfy the quad-linearity Q(k + h,f,f,f) = Q(k,f,f,f)+ Q(h,f,f,f). For more of the detailed linearization process, see Section 2. Now we are ready to state our main theorem. Theorem 1.1. Let N 3. Suppose that the initial data F satisﬁes ≥ 0 0 F0(x, p) 1 for fermions, 0 ≤ F (x, p) ≤ for bosons, ≤ 0 and the global equilibrium m(p) shares the same total mass, momentum and energy with the initial data: 1 1 (1.20) dxdp F0(x, p) µ = dxdp m(p) µ . T3 R3 p T3 R3 p Z × Z × Then there exist δ > 0 and C > 0 such that if (f0) δ then there exists a unique global- in-time solution of (1.19) such that E ≤ (1) The distribution function F (x,p,t) has the following bounds: 0 F (x,p,t) 1 for fermions, 0 ≤ F (x,p,t) ≤ for bosons. ≤ (2) The energy norm is bounded globally in time:

sup (f(t)) C (f0). t R+ E ≤ E ∈ (3) There exists a uniform constant ǫ > 0 such that the perturbation decays exponentially: α 2 ǫt ∂ f(t) 2 Ce− . k kLx,p ≤ α N | X|≤ (4) Let f and f¯ be the solutions with the initial data f0 and f¯0, respectively. Then there exists a positive constant δ > 0 such that δt f f¯ 2 e− f0 f¯0 2 . k − kLx,p ≤ k − kLx,p To prove the main theorem, we need a coercivity estimate of the linear operator L. The linear operator L satisﬁes the following dissipation property: 2 Lf,f 2 δ (I P )f , h iLv ≥ k − kν for some positive δ > 0. The macroscopic projection Pf denotes the orthonormal projection 2 onto Lp space with respect to the following 5-basis:

2 2 2 2 0 2 (1.21) m + τm ,p1 m + τm ,p2 m + τm ,p3 m + τm ,p m + τm . The 5-dimensionalnp basis abovep constitutesp the kernel of Lp. Then it is crucialp to obtaino some upper-bound estimates of the nonlinear terms Γ and T in order to construct the solution. Some parts of the nonlinear terms can be estimated by a simple change of variables and the Hölder inequality. However, regarding some second-order nonlinear terms whose integrands consist of the product of f(p′) (or f(q′)) and f(p), we need to take a change of variables of either p p′ (or p q′) and there occur some non-trivial difficulties. Different from the non-relativistic7→ case,7→ a uniform positive lower bound for the Jacobian for the change of variables ∂p′/∂p (or ∂q′/∂p ) does not exist as shown in [16], and hence we need another way to deal| with| this| difficulty.| One way is to proceed the estimates of the second-order nonlinear terms by lifting the dq integral to energy-momentum four vector integral dqµ imposing additional Dirac-delta function. This technique has been introduced in [18,44,45]. Then we reduce the integral by computing the Dirac-delta function. In this process, there 8 GI-CHANBAE,JINWOOJANG,ANDSEOK-BAEYUN appears a singularity of 1/g¯ and the exponential growth with respect to the p variable from 0 0 exp(p p′ ). We will explain this more in detail in Section 1.8.2. Once− we obtain the estimates of the linear terms and the nonlinear terms, we define an iteration scheme to construct the local-in-time solution as (1.22) (∂ +ˆp )F n+1 = Q(F n, F n, F n+1, F n). t · ∇x Regarding the case of fermions, the solution F has to be bounded by 1 from above. There- fore, we also need to prove that the function F n+1 is bounded in the closed interval [0, 1] in each iteration scheme based on the induction hypothesis that F n is in [0, 1]. Then we obtain that the collision operator Q(F n, F n, F n+1, F n) is well-defined. On the right-hand side of (1.22), note that we place F n+1 in the p variable position of the collision operator Q (i.e., at the third input of Q( , , , )). This is by the specific structure of the nonlinear operator Q from (1.14). Once we· place· · · F n+1 as a third input of Q then we can rewrite (1.22) in a more clear structure (1.23) ∂ +ˆp τG(F n)+ R(F n) F n+1 = G(F n), { t · ∇x − } where G and R is now defined as

G(F1, F2, F4)= dq dw vøσ(g,θ)F1(p′)F2(q′)(1 + τF4(q)), R3 S2 Z Z R(F1, F2, F4)= dq dw vøσ(g,θ)(1 + τF1(p′))(1 + τF2(q′))F4(q). R3 S2 Z Z Also, we can show that the boundedness 0 F n 1 implies that 0 G(F n) and 0 R(F n), which guarantees the boundedness of≤F n+1≤. Then the additional≤ linearization of≤ F n+1 = m + √m + τm2f n+1 results in creating the dissipation term ν(p)f n+1 in the left- hand side of (1.23), which would then result in the exponential decay of the perturbation f n+1 in the L2 sense. By the induction argument, we obtain the uniform boundedness of the energy locally in time. Then the standard way of extending the local existence to the global one is to eliminate the dissipation of the linear part L. Similarly to the Newtonian case, we substitute f = (I P )f + Pf on each side of (1.19) where Pf is defined by the orthonormal projection of − (1.21). Then the expansion of the linear term (∂t +p ˆ x)Pf yields a linear combination with respect to the 14-basis. Thus we can achieve the following· ∇ coercivity estimate: α α α 2 L∂ f, ∂ f 2 δ ∂ f , h iLx,p ≥ k kx,ν α N α N | X|≤ | X|≤ for some positive constant δ > 0. This full coercivity estimate enables the extension of the local-in-time solution to a global-in-time solution. 1.8. Main difficulties and our strategy. In this subsection, we present the difficulties that arise from the estimates of the nonlinear terms. In the relativistic quantum case, there appear new types of nonlinear terms involving both pre- and post-collisional momenta at the same time as f(p)f(p′) and f(p)f(q′). In general, this kind of terms has been expected to appear only in the non-cutoff Boltzmann theory in the non-quantum case. In the non-cutoff Boltzmann theory, it has been considered very crucial to understand the change of variables q p′ or q′ and the cancellation lemma to understand the fractional diffusive behavior [36,7→ 44]. In the non-relativistic case, these new nonlinear terms can be handled because ∂p′/∂p and ∂q′/∂q are bounded from below, as in [2,19]. However, in the relativistic case under| the| center-of-momentum| | frame, such a positive uniform lower-bound of the Jacobian of the collision map does not exist [16]. THE RELATIVISTIC QUANTUM BOLTZMANN EQUATION NEAR EQUILIBRIUM 9

The remedy of the issue has been introduced in [18,31,71] in the center-of-momentum frame (1.11) where the linear term K2 of the difficulty has been calculated. Namely, the authors lift the dq integral to the energy-momentum four vector integral dqµ by imposing an additional Dirac-delta function, and they take a change of variables via a suitable Lorentz transform. In this paper, we also follow a similar technique to deal with the nonlinearity that occurs in dealing with the terms Γ(f) and T (f). In this direction, we however still encounter several other additional difficulties on the nonlinear term involving f(p)f(p′) as below. Let us first denote the nonlinear term Γ2,1 involving f(p)f(p′) as

dp dq dp′ dq′ (1.24) Γ (f,h), η 2 sσ(g,θ) 2,1 Lp 0 0 0 0 h i ≤ R3 p R3 q R3 p R3 q Z Z Z ′ Z ′ (4) µ µ µ µ 0 0 δ (p + q p′ q′ )J(q )J(p′ /2) f(p) h(p′) η(p) , × − − | || || | and the integral part with respect to the measure dqdq′ as

dq dq′ (4) µ µ µ µ 0 (1.25) B = 0 0 sσ(g,θ)δ (p + q p′ q′ )J(q ). R3 q R3 q − − Z Z ′ Then we present the three main diﬃculties that arise regarding the nonlinear terms.

dp 1.8.1. The Lorentz-invariant measure p0 and the nonlinear structure of sσ(g,θ). The first dp difficulty arises from the Lorentz-invaraint measure p0 including the fraction of the energy and the nonlinear structure of the differential cross section sσ(g,θ) with respect to the collision variables p and q. A good way to remove the fraction of the energy is to lift dq and µ µ dq′ integrals to the energy momentum 4-vector dq and dq′ integrals by considering extra Dirac-delta and unit step functions as

µ µ (4) µ µ µ µ 0 0 0 µ µ B = dq dq′ sσ(g,θ)δ (p +q p′ q′ )J(q )u(q )u(q′ )δ(q qµ+1)δ(q′ qµ′ +1). R4 R4 − − Z Z µ (4) µ If we simply eliminate the dq′ integral by computing the Dirac-delta function of δ (p + µ µ µ q p′ q′ ), the rest of the integral becomes highly complicated to deal with. Thus, motivated− − by de Groot et al. in [18], we instead apply a symmetric change of variables µ µ µ µ µ µ q¯ = q + q′ andq ¯′ = q q′ . Despite the nonlinear structure of the s, g and cos θ in (1.4), (1.5) and (1.9), respectively,− we can also represent them as the terms that depend µ µ µ µ only on the variables p , p′ ,q ¯ andq ¯′ (See Lemma 3.5) as follows: 2 2 2 1 µ µ 2 g¯ gc =¯g (p + p′ )(¯qµ pµ pµ′ ), sc = gc +4, cos θc =1 2 2 . − 2 − − − gc Then we can have B in the following form 0 0 1 µ µ (4) µ µ µ q¯ +¯q′ B = dΘ(¯q , q¯′ )scσ(gc,θc)δ (p p′ +¯q′ )J , 4 R4 R4 − 2 Z × where µ µ µ µ 0 µ µ µ dΘ(¯q , q¯′ )= dq¯ dq¯′ u(¯q )u(¯s 4)δ((¯q q¯ +¯q′ q¯′ )+4)δ(¯q q¯′ ). − µ µ µ µ µ µ By substitutingq ¯′ = p′ p , we now reduce the 4-dimensional Dirac-delta function with − (4) µ µ µ µ the energy-momentum 4-vector δ (p + q p′ q′ ) and obtain a more intuitive form of B as − −

1 µ 0 µ µ µ (1.26) B = dq¯ u(¯q )u(¯s 4)δ(¯q q¯µ + (p′ p )(pµ′ pµ)+4) 4 R4 − − − Z 10 GI-CHANBAE,JINWOOJANG,ANDSEOK-BAEYUN

0 0 0 µ q¯ + p′ p δ(¯q (p′ p ))s σ(g ,θ )J − . × µ − µ c c c 2 1.8.2. Singularity with respect to g¯. Then another difficulty arises from the singularity in the relative momentumg ¯, which occurs from the reduction of the second Dirac-delta function. We first remark that the second Dirac-delta function of (1.26) consists of an inner product between energy momentum 4-vectors. Motivated by the explicit form of the Lorentz transform by Strain [71], we apply the Lorentz trasnform which converts p′ p to (0, 0, 0, g¯) µ − µ (i.e. Λ(p′ p )=(0, 0, 0, g¯)). Then we obtain µ − µ µ µ 3 1 3 δ(¯q (p′ p )) = δ(Λ¯q Λ(p′ p )) = δ(¯q g¯)= δ(¯q ), µ − µ µ − µ g¯ where we used that the Lorentz transform is invariant under the inner-product and that 1 δ(ax) = a δ(x). We can now see that oneg ¯ in the Dirac-delta function comes out as 1/g¯ and creates an additional singularity. However, motivated by (1.10), we use the half angle formula [18, page 277, (A.11)] of the cos θ and observe that 1 cos θ g¯2 sin2(θ/2) = − = . 2 g2 Combining with the assumption of the differential cross section, we obtain

θ θ 1 cos θ θ θ σ(g,θ)= g sin θ =2g sin cos =2g − cos = 2¯g cos . 2 2 r 2 2 2 This allows us to eliminate the singularity ofg ¯.

1.8.3. Exponential growth. The last diﬃculty is regarding the last multiplier on the right- hand side of the equation (1.26). Since the function J(¯q0/2) is contained in dq¯0 integral, 0 0 0 0 we have an exponential decay in p′ and an exponential growth in p from J (p′ p )/2 − at the same time. The decay for p 0 is beneﬁcial for the upper-bound estimate, but there ′ is a problematic term of the exponential growth of p0 even in (1.24). However, we prove that the remaining part of the right-hand side of (1.26) includes the following exponential decaying factor

0 0 2 2 (p + p′ ) p p′ exp | ×2 | . −s 4 − g¯ ! Then, by the estimates in [30, Lemma 3.1, (iii) and (iv)] of

0 0 2 2 2 2 (p + p′ ) p p′ 2 g¯ +4 g¯ 1 2 | × | = p p′ max +1, p p′ , 4 − g¯2 | − | 4¯g2 ≥ 4 4| − | we can have an exponential decay of exp( p p′ /2). Since the diﬀerence of the energy 0 0 −| − | p′ p can further be bounded by the diﬀerence of the momentum p′ p , the exponential growth| − can| be absorbed by the exponential decay as | − |

0 0 0 0 2 2 ′ ′ p′ p (p + p′ ) p p′ 1 (p 0 p0) 1 p p J − exp | × | e− 2 − e− 2 | − | 1. 2 − 4 − g¯2 ≤ ≤ s ! Then we can have properly weighted L2 bounds for the nonlinear terms. THE RELATIVISTIC QUANTUM BOLTZMANN EQUATION NEAR EQUILIBRIUM 11

1.9. Outline of the paper. This paper is organized as follows. In Section 2, we linearize the collision operator of the relativistic quantum Boltzmann equation nearby a global equilibrium. In Section 3, we establish several estimates on the linear and the nonlinear terms. Section 4 is devoted to constructing the unique local-in-time classical solution. In the last section, we prove the coercivity estimate and establish the global-in-time classical solution.

2. Linearization of the relativistic quantum Boltzmann equation In this section, we introduce the reformulation of the equation (1.1) via the linearization of the relativistic quantum Boltzmann equation nearby the global equilibrium: 1 m(p)= . eap0+c τ − Proposition 2.1. If we substitute F = m + √m + τm2f in (1.13), then we have ∂ f +ˆp f + Lf = Γ(f)+ T (f) t · ∇x where the linear term Lf is decomposed as following form: Lf = ν(p)f + K f K f, 1 − 2 where the collision frequency ν(p) is given by 1 (2.1) ν(p)= dq dw vøσ(g,θ)m(q)(1 + τm(p′))(1 + τm(q′)), 1+ τm(p) R3 S2 Z Z and the compact operator K1 and K2 are deﬁned by

2 2 K1f(p)= dq dw vøσ(g,θ) m + τm (p′) m + τm (q′)f(q), R3 S2 (2.2) Z Z p p 2 2 K2f(p)=2 dq dw vøσ(g,θ) m + τm (q) m + τm (q′)f(p′). R3 S2 Z Z The nonlinear term Γ(f) is represented as follows:p p 6 4 Γ(f)= Γi(f,f),T (f)= Ti(f,f,f). i=1 i=1 X X We denote the precise deﬁnition of nonlinear terms at the end of this proof. Proof. We substitute F = m + √m + τm2f in (1.13) to have m + τm2∂ f + m + τm2pˆ f = Q(m + m + τm2f). t · ∇x Dividing √m +pτm2 on each sidep gives an equation for the perturbationp f:

1 2 (2.3) ∂tf +ˆp xf = Q(m + m + τm f). · ∇ √m + τm2 To decompose the right-hand side into linear and nonlinearp terms, we ﬁrst deﬁne the zeroth- order-in-f term Q0:

Q0 = Q(m,m,m,m). As we can see in (1.14), the quantum collision operator includes (1 + τF ) terms. Because of these terms, we cannot have the quad-linear property of Q. In other words, we have Q(k + h,f,f,f) = Q(k,f,f,f)+ Q(h,f,f,f). So we deﬁne the following four ﬁrst-order-in- 6 f terms Q1, Q2, Q3 and Q4: Q = Q(m + m + τm2f,m,m,m) Q , 1 − 0 p 12 GI-CHANBAE,JINWOOJANG,ANDSEOK-BAEYUN

Q = Q(m,m + m + τm2f,m,m) Q , 2 − 0 Q = Q(m,m,mp+ m + τm2f,m) Q , 3 − 0 Q = Q(m,m,m,mp+ m + τm2f) Q . 4 − 0 In view of this notation, the collection of the zerothp and the first-order terms in f can be written as 4 1 Qi + Q0 . √ 2 m + τm i=1 ! X Thus we divide the right-hand side of (2.3) into the collection of zeroth, first-order terms and other terms as 1 Q(m + m + τm2f)= Lf + Γ(f), √m + τm2 − where p 1 4 1 Lf = Qi Q0, −√ 2 − √ 2 m + τm i=1 m + τm X and 1 Γ(f)= Q(m + m + τm2f)+ Lf. √m + τm2 p We first calculate the linear part Lf. By definition of Q0, we have

Q0 = dq dw vøσ(g,θ) m(p′)m(q′)(1 + τm(p))(1 + τm(q)) R3 S2 Z Z

(1 + τm(p′))(1 + τm(q′))m(p)m(q) . − We observe from the conservation of energy 0 0 0 0 p + q = p′ + q′ , that

(2.4) m(p′)m(q′)(1 + τm(p))(1 + τm(q)) = (1 + τm(p′))(1 + τm(q′))m(p)m(q), which gives Q0 =0. For the ﬁrst-order linear term, an explicit computation gives

2 Q1 = dq dw vøσ(g,θ) m + τm (p′)m(q′)(1 + τm(p))(1 + τm(q)) 3 2 R S Z Z p 2 τ m + τm (p )(1 + τm(q′))m(p)m(q) f(p′), − ′ p By (2.4) on the ﬁrst-order term Q1, we have

Q1 = dq dw vøσ(g,θ)m(p)m(q)(1 + τm(p′))(1 + τm(q′)) R3 S2 Z Z 2 2 m + τm (p′) m + τm (p′) τ f(p′), × p m(p′) − p1+ τm(p′) ! THE RELATIVISTIC QUANTUM BOLTZMANN EQUATION NEAR EQUILIBRIUM 13 which is equal to

f(p′) Q1 = dq dw vøσ(g,θ)m(p)m(q)(1 + τm(p′))(1 + τm(q′)) . R3 S2 m + τm2(p ) Z Z ′ With similar computations, combining with Q2,Q3 and Q4 yields p 1 Lf = dq dw vøσ(g,θ)m(p)m(q)(1 + τm(p′))(1 + τm(q′)) m + τm2(p) R3 S2 (2.5) Z Z f(p) f(q) f(p ) f(q ) p + ′ ′ . × 2 2 − 2 − 2 m + τm (p) m + τm (q) m + τm (p′) m + τm (q′) We can easily seep that the ﬁrst termp of Lf is equalp to νf. We deﬁnep the second term of Lf as K1f: m(p) m(q) (2.6) K1f = dq dw vøσ(g,θ) m + τm2(p) R3 S2 m + τm2(q) Z Z (1 + τm(p′))(1 + τm(q′))f(q). p p × We observe

m(p) 1 (ap0+c) (2.7) = e− 2 . m + τm2(p) Combining with the energy conservationp law gives m(p) m(q) m(p ) m(q ) (2.8) = ′ ′ . 2 2 2 2 m + τm (p) m + τm (q) m + τm (p′) m + τm (q′) Applying it top (2.6), we havep p p

2 2 K1f = dq dw vøσ(g,θ) m + τm (p′) m + τm (q′)f(q). R3 S2 Z Z We deﬁne the collection of the third and thep fourth term ofpLf in (2.5) as K f: − 2 m(p) 1+ τm(p′) (2.9) K2f = dq dw vøσ(g,θ)m(q) (1 + τm(q′))f(p′) m + τm2(p) R3 S2 m + τm2(p ) Z Z ′ m(p) 1+ τm(q′) + p dq dw vøσ(g,θ)m(q)(1p + τm(p′)) f(q′). m + τm2(p) R3 S2 m + τm2(q ) Z Z ′ We would callp the ﬁrst line of the K2f in (2.9) as K2,1f and thep second line of (2.9) as K2,2f. We then write the dw integral of K2,2 in the spherical coordinate w (φ, θ) as in (1.12) as follows: 7→ m(p) 2π π K2,2f = dq dφ sin θ dθ vøσ(g,θ) m + τm2(p) R3 Z Z0 Z0 1+ τm(q′) p m(q)(1 + τm(p′)) f(q′). 2 × m + τm (q′) Then we apply the change of variables θ π θ and φ π + φ. The change of variables → − → p would result in the exchanged roles of p′ and q′, since w in (1.12) changes into w. Thus we have − m(p) 3π 0 K2,2f = dq dφ sin θ ( dθ) m + τm2(p) R3 − Z Zπ Zπ p 14 GI-CHANBAE,JINWOOJANG,ANDSEOK-BAEYUN

1+ τm(p′) vøσ(g, π θ)m(q)(1 + τm(q′)) f(p′). 2 × − m + τm (p′) p By the assumption of the diﬀerential cross section σ(g,θ) in (1.18), we have σ(g,θ) = g sin θ = σ(g, π θ). This further gives −

m(p) 1+ τm(p′) K2,2f = dq dw vøσ(g,θ)m(q)(1 + τm(q′)) f(p′). m + τm2(p) R3 S2 m + τm2(p ) Z Z ′ p p This shows that K2,2 = K2,1. Thus we have

(2.10) K2f =2K2,1f

m(p) 1+ τm(p′) =2 dq dw vøσ(g,θ)m(q) (1 + τm(q′))f(p′). m + τm2(p) R3 S2 m + τm2(p ) Z Z ′ p p Similarly, we observe that

1+ τm(p) 1 ( 0+ ) (2.11) = e 2 ap c . m + τm2(p) p Combining with (2.7), we have

m(p) 1+ τm(p ) 1+ τm(q) m(q ) (2.12) ′ = ′ . 2 2 2 2 m + τm (p) m + τm (p′) m + τm (q) m + τm (q′) p p p p Substituting it in (2.10), we have

2 2 K2f =2 dq dw vøσ(g,θ) m + τm (q) m + τm (q′)f(p′). R3 S2 Z Z p p This completes the derivation of the linear term νf, K1f and K2f. Now we consider the nonlinear terms. Since the quantum collision operator is not quad- linear (Recall that Q(k+h,f,f,f) = Q(k,f,f,f)+Q(h,f,f,f)), the second-order nonlinear term is highly complicated to be represented.6 Thus, we ﬁrst observe one of the second-order nonlinear terms involving f(p′) and f(q′):

v σ(g,θ) dq dw ø R3 S2 m + τm2(p) Z Z 2 2 pm + τm (p′)f(p′) m + τm (q′)f(q′)(1 + τm(p))(1 + τm(q)) × 2 2 p p m + τm (p )f(p′) m + τm (q )f(q′)m(p)m(q) . − ′ ′ p p There are 6 second-order nonlinear terms of the same kind similar to the term above where the number 6 is coming from the number of choices for choosing 2 variables among the four variables p, q, p′, and q′. We represent all of the second-order nonlinear terms (6 nonlinear THE RELATIVISTIC QUANTUM BOLTZMANN EQUATION NEAR EQUILIBRIUM 15 terms) as follows: (2.13) v σ(g,θ) Γ(f,h)= dq dw ø R3 S2 m + τm2(p) Z Z 2 2 (m(p′)m(q′)p (1 + τm(p′))(1 + τm(q′))) m + τm (p)f(p) m + τm (q)h(q) × − p 2 p 2 + τ(m(q′)(1 + τm(q)) m(q)(1 + τm(q′))) m + τm (p)f(p) m + τm (p )h(p′) − ′ 2 2 + τ(m(q′)(1 + τm(p)) m(p)(1 + τm(q′)))pm + τm (q)f(q)pm + τm (p )h(p′) − ′ 2 2 + τ(m(p′)(1 + τm(q)) m(q)(1 + τm(p′)))pm + τm (p)f(p)pm + τm (q )h(q′) − ′ 2 2 + τ(m(p′)(1 + τm(p)) m(p)(1 + τm(p′)))pm + τm (q)f(q)pm + τm (q )h(q′) − ′ p 2 p 2 +(1+ τm(p))((1 + τm(q)) m(p)m(q)) m + τm (p )f(p′) m + τm (q )h(q′) − ′ ′ p p =Γ1(f,h)+Γ2(f,h)+Γ3(f,h)+Γ4(f,h)+Γ5(f,h)+Γ6(f,h). We also denote as

Γ(f)= Γi(f,f), 1 i 6 ≤X≤ when f = h. Lastly, we consider the following third-order nonlinear terms involving f(p′), f(q′) and f(p):

2 2 dq dw vøσ(g,θ) τ m + τm (p′)f(p′) m + τm (q′)f(q′)f(p)(1 + τm(q)) R3 S2 Z Z p p2 2 m + τm (p )f(p′) m + τm (q )f(q′)f(p)m(q) . − ′ ′ Similarly, we represent all of the third-orderp nonlinear termsp (4 nonlinear terms) as

2 2 T1(f,h,η)= τ dq dw vøσ(g,θ)f(p) m + τm (q)h(q) m + τm (p′)η(p′), − R3 S2 Z Z p p 2 2 T2(f,h,η)= τ dq dw vøσ(g,θ)f(p) m + τm (q)h(q) m + τm (q′)η(q′), − R3 S2 Z Z p p (2.14) 2 2 T3(f,h,η)= τ dq dw vøσ(g,θ)f(p) m + τm (p′)h(p′) m + τm (q′)η(q′), R3 S2 Z Z p p 1 2 T4(f,h,η)= τ dq dw vøσ(g,θ) m + τm (q)f(q) m + τm2(p) R3 S2 Z Z 2 2 p p m + τm (p′)h(p′) m + τm (q′)η(q′). × Similarly, we deﬁne p p

T (f)= Ti(f,f,f). 1 i 4 ≤X≤ We can easily check that the fourth-order nonlinear term is cancelled.

By the linearization proposition above and substituting F = m + √m + τm2f in (1.13), we obtain the linearized equation for the relativistic quantum Boltzmann model (1.13) as 16 GI-CHANBAE,JINWOOJANG,ANDSEOK-BAEYUN follows:

∂tf +p ˆ xf + Lf = Γ(f)+ T (f), (2.15) · ∇ f(x, p, 0) = f0(x, p).

2 where f0(x, p) = (F0(x, p) m)/√m + τm . Then the conservation laws (1.15) can be written as following form: −

2 2 dx dp f(x,p,t) m + τm = dx dp f0(x, p) m + τm , T3 R3 T3 R3 Z Z p Z Z p 2 2 (2.16) dx dp f(x,p,t)p m + τm = dx dp f0(x, p)p m + τm , T3 R3 T3 R3 Z Z p Z Z p 0 2 0 2 dx dp f(x,p,t)p m + τm = dx dp f0(x, p)p m + τm . T3 R3 T3 R3 Z Z p Z Z p Now we state several useful properties for the linear term Lf.

Lemma 2.1. For any smooth function φ, we have

dp dq dw vøσ(g,θ)m(p)m(q)(1 + τm(p′))(1 + τm(q′)) R3 R3 S2 Z Z Z f(p) f(q) f(p ) f(q ) + ′ ′ φ(p) × 2 2 − 2 − 2 m + τm (p) m + τm (q) m + τm (p′) m + τm (q′)

= pdp dq dwp vøσ(g,θ)m(p)m(qp)(1 + τm(p′))(1 +pτm(q′)) R3 R3 S2 Z Z Z f(p) f(q) f(p ) f(q ) + ′ ′ φ(q) × 2 2 − 2 − 2 m + τm (p) m + τm (q) m + τm (p′) m + τm (q′)

= pdp dq dwp vøσ(g,θ)m(p)m(qp)(1 + τm(p′))(1 +pτm(q′)) R3 R3 S2 Z Z Z f(p) f(q) f(p′) f(q′) + ( φ(p′)) × 2 2 − 2 − 2 − m + τm (p) m + τm (q) m + τm (p′) m + τm (q′)

= pdp dq dwp vøσ(g,θ)m(p)m(qp)(1 + τm(p′))(1 +pτm(q′)) R3 R3 S2 Z Z Z f(p) f(q) f(p′) f(q′) + ( φ(q′)). × 2 2 − 2 − 2 − m + τm (p) m + τm (q) m + τm (p′) m + τm (q′)

Proof. Notep that m(p)m(qp)(1+τm(p′))(1+τmp (q′)) is invariantp under the change of variables (p, q) (p′, q′) from (2.4). Thus, by the change of variables p q and (p, q) (p′, q′) introduced↔ in [29], we have the desired results. ↔ ↔

Lemma 2.2. We have the following properties for the linear operator L:

(1) L is a symmetric operator: Lf,g 2 = f,Lg 2 . h iLp h iLp (2) Lf =0 if and only if f = Pf. 2 where Pf is deﬁned as the orthonormal projection to the Lp space which is spanned by following 5-dimensional basis:

2 2 2 2 0 2 m + τm ,p1 m + τm ,p2 m + τm ,p3 m + τm ,p m + τm . np p p p p o THE RELATIVISTIC QUANTUM BOLTZMANN EQUATION NEAR EQUILIBRIUM 17

Proof. (1) We take an inner product of g with Lf of (2.5) to have g(p) dp gLf = dp dq dw vøσ(g,θ)m(p)m(q)(1 + τm(p′))(1 + τm(q′)) R3 R3 m + τm2(p) R3 S2 Z Z Z Z f(p) f(q) f(p ) f(q ) p + ′ ′ . × 2 2 − 2 − 2 m + τm (p) m + τm (q) m + τm (p′) m + τm (q′) By Lemma 2.1,p we substitute p p p 1 φ(p)= g(p) m + τm2(p) and obtain p 1 dp gLf = dp dq dw vøσ(g,θ)m(p)m(q)(1 + τm(p′))(1 + τm(q′)) R3 4 R3 R3 S2 Z Z Z Z f(p) f(q) f(p ) f(q ) + ′ ′ 2 2 2 2 × m + τm (p) m + τm (q) − m + τm (p′) − m + τm (q′)! p g(p) p g(q) p g(p ) p g(q ) + ′ ′ . 2 2 2 2 × m + τm (p) m + τm (q) − m + τm (p′) − m + τm (q′)! This proves the symmetricityp of L. p p p (2) Once we substitute g = f in the equation above, then we have 1 dp fLf = dp dq dw vøσ(g,θ)m(p)m(q)(1 + τm(p′))(1 + τm(q′)) R3 4 R3 R3 S2 Z Z Z Z f(p) f(q) f(p ) f(q ) 2 + ′ ′ × 2 2 − 2 − 2 m + τm (p) m + τm (q) m + τm (p′) m + τm (q′) 0. ≥ p p p p Therefore Lf = 0 implies f(p) f(q) f(p ) f(q ) + = ′ + ′ , 2 2 2 2 m + τm (p) m + τm (q) m + τm (p′) m + τm (q′) which is satisﬁedp if and onlyp if f = Pf. Moreover,p L(Pf) = 0 givesp the desired results.

3. Estimates of the linear and the nonlinear terms This section is devoted to proving several estimates of the linear and the nonlinear terms. Especially, we emphasize that the estimates of the nonlinear term Γ2 involve the main diﬃculties as we mentioned in Section 1.8. Before we move onto it, we present some useful properties that are commonly used throughout this paper.

Lemma 3.1. There exist positive constants C1 > 0 and C2 > 0 such that C J(p0) m(p) C J(p0). 1 ≤ ≤ 2 c c Proof. In the fermionic case, C1 =1/(1 + e ) and C2 = e− give the desired estimates: 1 1 1 1 . ec +1 eap0 ≤ eap0+c +1 ≤ eap0+c 18 GI-CHANBAE,JINWOOJANG,ANDSEOK-BAEYUN

c c a In the case of bosons, we choose C = e− and C = 1/(e e− ). Since c > a, we have 1 2 − − C2 < , and we get ∞ 1 1 1 1 ap0+c ap0+c c a ap0 . e ≤ e 1 ≤ e e− e − −

By this Lemma, the relativistic quantum equilibrium m(p) can be treated as the non- quantum relativistic equilibrium J(p0). Now we present the estimates for collision frequency ν(p) and operator K1 and K2. Lemma 3.2. There exists a positive constant C > 0 such that

1 0 1 0 1 (p ) 2 ν(p) C(p ) 2 . C ≤ ≤ 2 The integral operator Ki(f) is a compact operator in Lp for i =1, 2. Proof. By Lemma 3.1,

0 eap +c ecm(p) (3.1) 1+ τm(p)= = C. eap0+c τ J(p0) ≤ − Then by (3.1), we have the upper bound of ν as

ν(p) C dq dw vøσ(g,θ)m(q). ≤ R3 S2 Z Z Similarly we have

1 0 2 (ap +c) e 1 0 0 (3.2) m + τm2(p)= Ce 2 ap = CJ(p /2), eap0+c τ ≤ − which yields p

0 0 K1f(p) C dq dw vøσ(g,θ)J(p′ /2)J(q′ /2)f(q), and ≤ R3 S2 Z Z 0 0 K2f(p) C dq dw vøσ(g,θ)J(q /2)J(q′ /2)f(p′). ≤ R3 S2 Z Z We note that the upper bounds of ν and Ki are identical to those of non-quantum relativistic linear terms in [23]. Therefore we obtain the desired results.

Lemma 3.3. The linear operator L satisﬁes following dissipation property for some positive constant δ: 2 Lf,f 2 δ (I P )f . h iLv ≥ k − kν Proof. Since K1 and K2 are compact and ν satisﬁes Lemma 3.1, we have the desired results by following the same proof of Lemma 3.4 of [71].

3.1. Estimates of the second-order nonlinear terms. In this subsection, we establish the estimates on the second-order nonlinear terms. Lemma 3.4. We have

Γ(f,h), η 2 C f 2 h ν + f ν h 2 η ν. h iLp ≤ k kLp k k k k k kLp k k

THE RELATIVISTIC QUANTUM BOLTZMANN EQUATION NEAR EQUILIBRIUM 19

Proof. As we can see in (2.13), there are six second-order nonlinear terms. By the change of variables p′ q′, we can see that Γ (f,h)=Γ (f,h) and Γ (f,h)=Γ (f,h). Thus we ↔ 2 4 3 5 write the proof in three parts. Firstly, we prove the estimates of Γ1(f,h) and Γ6(f,h), since they are similar to the non-quantum relativistic case. Secondly, we present the estimate of Γ3(f,h). Lastly, we prove the most difficult part Γ2(f,f). For this we need several techniques to deal with the difficulties that arise from the relativistic integrals. (1) Estimates of Γ1 and Γ6: We first consider Γ1(f,h) term. By (3.1), (3.2), and Lemma 3.1, we have

0 2 (3.3) Γ1(f,h), η Lp C dp dq dw vøσ(g,θ)J(q /2) f(p) h(q) η(p) . h i ≤ R3 R3 S2 | || || | Z Z Z

By the H¨older inequality, we have 1 1 2 2 2 2 0 2 2 Γ1(f,h), η Lp C dp dw dq vøσ (g,θ)J(q ) dq h(q) f(p) η(p) . h i ≤ R3 S2 R3 R3 | | | || | Z Z Z Z

For the ﬁrst dq integral, we recall the deﬁnition of vø and σ(g,θ): g√s v = , σ(g,θ)= g sin θ. ø p0q0 By an explicit computation, we have s =2+2p0q0 2p q 2p0q0 + 2(1 cos θ) p q 4p0q0, − · ≤ − | || | ≤ and g = √s 4 √s 2 p0q0. − ≤ ≤ Applying above inequalities, we have p

2 2 0 0 dq vøσ (g,θ)J(q ) Cp . R3 ≤ Z We apply the H¨older inequality again to obtain the following estimate:

0 1 2 2 2 2 Γ1(f,h), η Lp C h Lp dp dw (p ) f(p) η(p) C f ν h Lp η ν. h i ≤ k k R3 S2 | || |≤ k k k k k k Z Z 0 0 0 0 Now we consider Γ6(f,h ) term. We ﬁrst use the energy conservation law p + q = p′ + q′ to have 2 2 0 0 0 0 m + τm (p ) m + τm (q ) CJ(p′ /2)J(q′ /2) = CJ(p /2)J(q /2), ′ ′ ≤ which givesp p

0 2 (3.4) Γ6(f,h), η Lp C dp dq dw vøσ(g,θ)J(q /2) f(p′) h(q′) η(p) . h i ≤ R3 R3 S2 | || || | Z Z Z

Applying the H¨older inequality yields 1 2 2 0 2 Γ6(f,h), η Lp C dp dw dq vøσ (g,θ)J(q ) h i ≤ R3 S2 R3 Z Z Z 1 2 2 2 dq vø f(p′) h(q′) η(p) . × R3 | | | | | | Z 0 0 0 1 0 1 Similarly we use s 4p q and g 2(p ) 2 (q ) 2 to have ≤ ≤ 1 2 0 1 2 2 2 2 Γ6(f,h), η Lp C dp dw (p ) dq vø f(p′) h(q′) η(p) . h i ≤ R3 S2 R3 | | | | | | Z Z Z

20 GI-CHANBAE,JINWOOJANG,ANDSEOK-BAEYUN

By the H¨older inequality for the dpdw integral, we have 1 2 0 1 2 2 2 2 Γ6(f,h), η Lp C η ν dp dw (p ) dqvø f(p′) h(q′) . h i ≤ k k R3 S2 R3 | | | | Z Z Z

The energy conservation law implies 0 1/2 0 0 1/2 0 1/2 0 1/2 (p ) (p′ + q′ ) (p′ ) + (q′ ) , ≤ ≤ which yields

0 1 0 1 2 2 2 2 2 Γ6(f,h), η Lp η ν dw dp dq vø (p′ ) + (q′ ) f(p′) h(q′) . h i ≤k k S2 R3 R3 | | | | Z Z Z Then we consider the pre-post change of variables (p, q) (p′, q′) with the Jacobian in [32]: ↔ dpdq dp′dq′ (3.5) 0 0 = 0 0 , p q p′ q′ and recall from (1.7) that we already have µ µ µ µ µ µ µ µ s(p , q )= s(p′ , q′ ), g(p , q )= g(p′ , q′ ). Thus we have

0 1 0 1 2 2 dw dp dq vø (p′ ) 2 + (q′ ) 2 f(p′) h(q′) . S2 R3 R3 | | | | Z Z Z µ µ µ µ g(p′ , q′ ) s(p , q ) 1 1 ′ ′ 0 2 0 2 2 2 = dw dp′ dq′ 0 0 (p′ ) + (q′ ) f(p′) h(q′) . S2 R3 R3 p q | | | | Z Z Z ′p′ µ µ Finally by v (p′ , q′ ) C and the H¨older inequality, we obtain the desired results: ø ≤ Γ6(f,h), η 2 C f 2 h ν + f ν h 2 η ν . h iLp ≤ k kLp k k k k k kLp k k

(2) Estimates of Γ3: We split the Γ3 term into two parts: τ Γ3(f,h)= dq dw vøσ(g,θ) m(q′)(1 + τm(p)) m(p)(1 + τm(q′)) m + τm2(p) R3 S2 { − } Z Z 2 2 p m + τm (q)f(q) m + τm (p′)h(p′) × =Γ3p,1(f,h)+Γ3,2(f,h)p.

By 1+ τm(p) C and 1+ τm(q′) C, we obtain ≤ ≤ 0 0 0 J(q′ )J(q /2)J(p′ /2) Γ3,1(f,h) C dq dw vøσ(g,θ) 0 f(q) h(p′) | |≤ R3 S2 J(p /2) | || | Z Z 0 0 C dq dw vøσ(g,θ)J(q )J(q′ /2) f(q) h(p′) , ≤ R3 S2 | || | Z Z and

0 0 0 Γ3,2(f,h) C dq dw vøσ(g,θ)J(p /2)J(q /2)J(p′ /2) f(q) h(p′) , | |≤ R3 S2 | || | Z Z 0 0 0 0 where we used J(p′ /2)J(q′ /2) = J(p /2)J(q /2) for the exponential term. In the case of 0 0 Γ3,2, there are decays of two precollisional momentum variables J(p /2)J(q /2). Thus we can have the exponential decay with respect to all the pre-post momentum variables by the energy conservation: 0 0 0 0 0 0 J(p /2)J(q /2) = J(p /4)J(q /4)J(p′ /4)J(q′ /4). THE RELATIVISTIC QUANTUM BOLTZMANN EQUATION NEAR EQUILIBRIUM 21

Thus, the estimate of Γ3,2 can be absorbed in that of Γ3,1. Thus, without loss of generality, we only consider the estimates of Γ3,1:

0 0 2 Γ3,1(f,h), η Lp C dp dq dw vøσ(g,θ)J(q )J(q′ /2) f(q) h(p′) η(p) . h i ≤ R3 R3 S2 | || || | Z Z Z

To separate pre-collisional variable and post-collisional variable, we apply the H¨older inequality as follows: 1 2 0 0 2 2 2 Γ3,1(f,h), η Lp C dp dq dw vøσ(g,θ)J(q )J(q′ /2) f(q) η(p) h i ≤ R3 R3 S2 | | | | Z Z Z 1 2 0 0 2 dp dq dw vøσ(g,θ)J(q )J(q′ /2) h(p′) × R3 R3 S2 | | Z Z Z = II II . 1 × 2 0 0 0 1 0 1 Using s 4p q and g 2(p ) 2 (q ) 2 , we have ≤ ≤ 2 0 1 0 1 0 0 2 2 (II1) C dp dq dw (p ) 2 (q ) 2 J(q )J(q′ /2) f(q) η(p) . ≤ R3 R3 S2 | | | | Z Z Z 0 1 0 0 Since (q ) 2 J(q /2) and J(q′ /2) are bounded by constants, we obtain

1 2 2 0 2 2 2 2 (II1) C dq f(q) dp (p ) η(p) C f L2 η ν. ≤ R3 | | R3 | | ≤ k k p k k Z Z For II2, we apply the pre-post momentum change of variables (p, q) (p′, q′) similarly to g√s ↔ the estimate of Γ6(f,h). Note that vø = p0q0 is invariant under the change of variables, and we obtain µ µ µ µ 2 g(p′ , q′ ) s(p′ , q′ ) µ µ 0 0 2 (II2) = dp′ dq′ dw 0 0 σ(g(p′ , q′ ),θ)J(q )J(q′ /2) h(p′) . R3 R3 S2 p q | | Z Z Z ′p′ µ µ 0 1/2 0 1/2 µ µ 0 0 Similarly, we use g(p′ , q′ ) 2(p′ ) (q′ ) and s(p′ , q′ ) 4p′ q′ to have ≤ ≤ 2 0 1 0 1 0 0 2 (II2) dp′ dq′ dw (p′ ) 2 (q′ ) 2 J(q )J(q′ /2) h(p′) ≤ R3 R3 S2 | | Z Z Z 0 1 0 0 1 2 dw dq′ (q′ ) 2 J(q′ /2) dp′ (p′ ) 2 h(p′) ≤ S2 R3 R3 | | Z Z Z C h 2. ≤ k kν Combining the estimates of II1 and II2 yields

Γ3,1(f,h), η 2 C f 2 h ν η ν. h iLp ≤ k kLp k k k k

(3) Estimates of Γ2: We ﬁrst separate Γ2 by two terms:

Γ2(f,h)= τ dq dw vøσ(g,θ) m(q)(1 + τm(q′)) m(q′)(1 + τm(q)) − R3 S2 − Z Z 2 m + τm (p )f(p)h(p′) × ′ =Γ (f,h) Γ (f,h). 2p,1 − 2,2 Using 1 + τm(p) C, we estimate each term as follows: ≤ 0 0 Γ2,1(f,h) C dq dw vøσ(g,θ)J(q )J(p′ /2) f(p) h(p′) , ≤ R3 S2 | || | Z Z

22 GI-CHANBAE,JINWOOJANG,ANDSEOK-BAEYUN and

0 0 Γ2,2(f,h) C dq dw vøσ(g,θ)J(q′ )J(p′ /2) f(p) h(p′) . ≤ R3 S2 | || | Z Z

Since Γ2,2 has the exponential decays with respect to two post-collisional variables p′ and q′, we can have the decays for all pre- and post-collisional momentum variables by using the energy conservation: 0 0 0 0 0 0 J(p /2)J(q /2) = J(p /4)J(q /4)J(p′ /4)J(q′ /4).

Thus the estimate of Γ2,2 can be absorbed into that of Γ2,1. Thus, we only consider Γ2,1 term. However, different from the previous cases, it is difficult to separate f,h and η using the Hölder inequality. Thus, we proceed the estimate via lifting the dw integral to dp′dq′ integral imposing a 4-dimensional Dirac-delta function:

dp dq dp′ dq′ (3.6) Γ (f,h), η 2 C sσ(g,θ) 2,1 Lp 0 0 0 0 h i ≤ R3 p R3 q R3 p R3 q Z Z Z ′ Z ′ (4) µ µ µ µ 0 0 δ (p + q p′ q′ )J(q )J(p′ /2) f(p) h(p′) η(p) . × − − | || || | Then we focus on the following estimate:

dq dq′ (4) µ µ µ µ 0 (3.7) B = 0 0 sσ(g,θ)δ (p + q p′ q′ )J(q ). R3 q R3 q − − Z Z ′ We also lift the dq and dq′ integrals to the relativistic energy-momentum 4-vector integral µ µ dq and dq′ imposing Dirac-delta function and unit step function: (3.8) µ µ (4) µ µ µ µ 0 0 0 µ µ B = dq dq′ sσ(g,θ)δ (p +q p′ q′ )J(q )u(q )u(q′ )δ(q qµ+1)δ(q′ qµ′ +1), R4 R4 − − Z Z where the unit step function u(x) is defined by 1 when x 0 and 0 when x< 0. Note that µ µ ≥ µ µ B has the integration over dq and dq′ , but s and g are functions depending on p and q . µ µ Thus it is convenient to split the q and the q′ parts from g as in the lemma below. Before we proceed furtherly, we first provide a preliminary lemma on the properties of g, g¯ andg ˜. Lemma 3.5 ( [44, 47, 71]). Define g, g¯ and g˜ as (1.7). Then we have (1) g2 =¯g2 +˜g2,

2 1 µ µ (2)g ¯ = (p + q′ )(q + p′ p q′ ), −2 µ µ − µ − µ 2 1 µ µ (3)g ˜ = (p + p′ )(q + q′ p p′ ). −2 µ µ − µ − µ Proof. The proofs can be found in [47,71], but for the readers’ convenience, we present the proof. (1) By deﬁnition of g in (1.7), we have 2 2 2 µ µ µ g +¯g +˜g = 2+2p q 2p p′ 2p q′ . − − µ − µ − µ Since pµ is energy momentum 4-vector, we have pµp = 1, which gives µ − 2 2 2 µ µ µ µ g +¯g +˜g =2p p +2p q 2p p′ 2p q′ − µ µ − µ − µ µ =2p (p + q p′ q′ ). µ µ − µ − µ THE RELATIVISTIC QUANTUM BOLTZMANN EQUATION NEAR EQUILIBRIUM 23

Then the conservation law of energy momentum 4-vector lead to the desired results. (2) We denote the right-hand side as R and expand as follows:

1 µ µ µ µ 1 µ µ µ R = (p q + p p′ + q′ q + q′ p′ )+ (p p +2p q′ + q′ q′ ). −2 µ µ µ µ 2 µ µ µ µ µ µ µ We recall from (1.7) that p qµ = p′ qµ′ and p pµ′ = q qµ′ to have

1 µ µ 1 µ µ µ µ R = (2p q +2p p′ )+ ( 2+2p q′ )= p q + p q′ p p′ 1. −2 µ µ 2 − µ − µ µ − µ − Using pµp = 1, we have µ − µ R = p ( q + q′ p′ + p ). − µ µ − µ µ Then the conservation laws of energy momentum 4-vector yields µ µ R =2p ( p′ + p )= 2p p′ 2. − µ µ − µ − By deﬁnition ofg ¯, we derived desired results. (3) The proof can be obtained by the same way as that of (2). So we omit it. By Lemma 3.5 (1) and (3) above, we use the following representation of g:

2 2 1 µ µ g =g ¯ (p + p′ )(q + q′ p p′ ). − 2 µ µ − µ − µ Now we go back to the estimates of B and apply the following change of variables µ µ µ µ µ µ (3.9) q¯ = q + q′ , q¯′ = q q′ . − Then the reverse relation can be written by

µ 1 µ µ µ 1 µ µ q = (¯q +¯q′ ), q′ = (¯q q¯′ ), 2 2 − and the Jacobian is given by µ µ ∂(q , q′ ) 1 µ µ = . ∂(¯q , q¯′ ) 16

Then the gc, sc and θc are now expressed as 2 2 2 1 µ µ 2 g¯ (3.10) gc =¯g (p + p′ )(¯qµ pµ pµ′ ), sc = gc +4, cos θc =1 2 2 . − 2 − − − gc To calculate the Dirac-delta function and the unit step function in (3.8), we use followings. First we use δ(x)δ(y)=2δ(x + y)δ(x y) to have − µ µ µ µ µ µ δ(q q + 1)δ(q′ q′ +1)=2δ(q q + q′ q′ + 2)δ(q q q′ q′ ) µ µ µ µ µ − µ µ µ µ =4δ((¯q q¯µ +¯q′ q¯µ′ )+4)δ(¯q q¯µ′ ). 0 0 0 0 0 0 0 0 Note that q 0 and q′ 0 are equivalent to q + q′ 0 and q q′ 0. Also, q q′ 0 ≥ 2 0 ≥0 ≥ ≥ µ ≥ is equivalent tog ¯ = 2q q′ 2q q′ 2 0 under the assumption that q qµ +1=0 and µ − · − ≥ q′ qµ′ + 1 = 0. Thus we have 0 0 µ µ 0 µ µ µ u(q )u(q′ )δ(q q + 1)δ(q′ q′ +1)= u(¯q )u(¯s 4)4δ((¯q q¯ +¯q′ q¯′ )+4)δ(¯q q¯′ ). µ µ − µ µ µ Thus we have 0 0 1 µ µ (4) µ µ µ q¯ +¯q′ B = dΘ(¯q , q¯′ )scσ(gc,θc)δ (p p′ +¯q′ )J , 4 R4 R4 − 2 Z × where µ µ µ µ 0 µ µ µ dΘ(¯q , q¯′ )= dq¯ dq¯′ u(¯q )u(¯s 4)δ((¯q q¯ +¯q′ q¯′ )+4)δ(¯q q¯′ ). − µ µ µ 24 GI-CHANBAE,JINWOOJANG,ANDSEOK-BAEYUN

µ µ µ µ Now we substituteq ¯′ = p′ p to reduce the 4-dimensional integration dq¯′ by reducing the 4-dimensional Dirac-delta− function as follows: 0 0 0 1 µ q¯ + p′ p B = dΘ(¯q )scσ(gc,θc)J − , 4 R4 2 Z where µ µ 0 µ µ µ µ dΘ(¯q )= dq¯ u(¯q )u(¯s 4)δ(¯q q¯ + (p′ p )(p′ p )+4)δ(¯q (p′ p )). − µ − µ − µ µ − µ To remove one more Dirac-delta function, we follow that 0 µ µ µ 0 µ µ u(¯q )δ(¯q q¯ + (p′ p )(p′ p )+4)= u(¯q )δ(¯q q¯ 2p p′ + 2) µ − µ − µ µ − µ = u(¯q0)δ( (¯q0)2 + q¯ 2 +¯s) − | | δ(¯q0 q¯ 2 +¯s) = − | | . 2 q¯ 2 +¯s |p| Then the dq¯0 integral with the delta function above is reducedp as follows: 0 0 0 1 p′ p dq¯ q¯ µ B = J − 0 scσ(gc,θc)J u(¯s 4)δ(¯q (pµ′ pµ)), 4 2 R3 q¯ 2 − − Z whereq ¯0 is deﬁned by q¯ 2 +¯s. Sinces ¯ 4=¯g2 0, the last unit step function is always equal to 1. In order to make| | the ﬁnal Dirac-delta− ≥ function even simpler, we consider the p Lorentz transform satisfying

Λ(p + p′ ) = (√s,¯ 0, 0, 0), Λ(p p′ )=(0, 0, 0, g¯). µ µ µ − µ − We can see that in [71] that the Lorentz transform satisfying the relation above is uniquely determined by following form: 0 ′0 ′ ′ ′ p +p p1+p1 p2+p2 p3+p3 √s¯ − √s¯ − √s¯ − √s¯ Λ01 Λ11 Λ21 Λ31  ′ ′ ′  (3.11) Λ= (p p )1 (p p )2 (p p )3 , 0 p× p′ p× p′ p× p′  0 ′0 | × |′ | × |′ | × |′   p p p1 p1 p2 p2 p3 p3   − − − −   g¯ − g¯ − g¯ − g¯  where   0 0 µ 0 0 µ 2(pi p + p′ p p′ + p′ p′ + p p p′ ) Λi1 = { µ} i{ µ} . g¯√s¯ p p | × ′| Then the exponential part can also be expressed carrying out the Lorentz transform as follows: J q¯0/2 = exp( q¯ /2) = exp( q¯µU /2) = exp( Λ¯qµΛU /2) , − 0 − µ − µ where we used following simple four-vectors: U µ = (1, 0, 0, 0),U = ( 1, 0, 0, 0). µ − Applying this Lorentz transform yields 0 0 1 p′ p dq¯ µ µ B = J − 0 sΛσ(gΛ,θΛ)exp( Λ¯q ΛUµ/2) δ(Λ¯q Λ(pµ′ pµ)), 4 2 R3 q¯ − − Z where gΛ is written via the Lorentz transform:

2 2 1 µ µ g =¯g Λ(p + p′ )Λ(¯q p p′ ) Λ − 2 µ − µ − µ THE RELATIVISTIC QUANTUM BOLTZMANN EQUATION NEAR EQUILIBRIUM 25

1 =¯g2 + (√s,¯ 0, 0, 0)(Λ¯q (√s,¯ 0, 0, 0)). 2 µ − We can also represent sΛ and cos θΛ in terms ofg ¯ and gΛ as 2 2 g¯ sΛ = gΛ +4 and cos θΛ =1 2 2 . − gΛ We now apply the change of variables Λ¯qµ =q ¯µ. Since dq/¯ q¯0 is Lorentz invariant, we have 0 0 1 p′ p dq¯ µ 3 B = J − 0 sλσ(gλ,θλ)exp( q¯ ΛUµ/2) δ(¯q g¯), 4 2 R3 q¯ − Z where 1 g¯2 2 2 √ 0 √ 2 (3.12) gλ =g ¯ + s¯(¯q s¯), sλ = gλ +4, cos θλ =1 2 2 . 2 − − gλ Now we use the spherical coordinates for the variableq ¯: q¯ = q¯ (sin ψ cos φ, sin ψ sin φ, cos ψ). | | Then B is equal to 0 0 2π π 2 1 p′ p ∞ q¯ d q¯ J − dφ dψ sin ψ | | | |s σ(g ,θ )exp( q¯µΛU /2) δ( q¯ cos ψg¯). 4 2 q¯0 λ λ λ − µ | | Z0 Z0 Z0 Using δ(ax) =(1/a)δ(x) and substituting ψ = π/2 reduce the dψ integral with δ(cos ψ) as follows: 0 0 2π 1 p′ p ∞ q¯ d q¯ (3.13) B = J − dφ | | | |s σ(g ,θ )exp( q¯µΛU /2) . 4 2 g¯q¯0 λ λ λ − µ |ψ=π/2 Z0 Z0 We ﬁrst consider the scattering kernel sλσ(gλ,θλ). The half-angle formula from (3.12) gives 2 2 θλ g¯ cos θλ =1 2 sin =1 2 2 , − 2 − gλ which implies θ g¯ (3.14) sin λ = , 2 gλ for 0 θ 2π. Thus we have ≤ λ ≤ θ θ θ σ(g ,θ )= g sin θ =2g sin λ cos λ = 2¯g cos λ 2¯g. λ λ λ λ λ 2 2 2 ≤ From (3.12) and the deﬁnition ofq ¯0, we have 1 1 s = g2 +4=¯g2 +4+ √s¯(¯q0 √s¯)=¯s + √s¯( q¯ 2 +¯s √s¯). λ λ 2 − 2 | | − Then we consider the exponential part of (3.13). Using the Lorentp z transform (3.11), we have 0 0 0 0 p + p′ 2 p p′ p p′ ΛUµ = , | × |, 0, − . √s¯ g¯√s¯ g¯ Combining with the spherically expression ofq ¯, we calculate 0 0 0 0 µ 2 p + p′ 2 p p′ p p′ q¯ ΛUµ ψ=π/2 = ( q¯ +¯s, q¯ cos φ, q¯ sin φ, 0) , | × |, 0, − | | | | | | | √s¯ g¯√s¯ g¯ p 0 0 p + p′ 2 p p′ = q¯ 2 +¯s + q¯ cos φ | × |. − | | √s¯ | | g¯√s¯ p 26 GI-CHANBAE,JINWOOJANG,ANDSEOK-BAEYUN

3/2 0 0 2 s¯ p′ p ∞ ydy 1+ y +1 B J − ≤ 2 2 2 2 Z0 y +1 p 0 0 2π p + p′ p p′ pexp y2 +1 dφ exp y cos φ| × | , × − 2 0 g¯ p Z For the notational simplicity, we denote 0 0 p + p′ p p′ (3.15) R = , r = | × |, 2 g¯ then we can write 3/2 0 0 s¯ π p′ p ∞ y 2 B J − dy + y exp R y +1 I0(ry), ≤ 2 2 0 y2 +1 ! − Z p where I0 denotes the modified Bessel functionp of the first kind: 1 π I (y)= ey cos φdφ. 0 π Z0 The formula including above Bessel function can be found in [35] and [30] when R > r 0: ≥ √R2 r2 ∞ y R√y2+1 e− − I = e− I0 (ry) dy = , 2 2 2 0 y +1 √R r (3.16) Z − ∞ p R√y2+1 R 1 √R2 r2 II = ye− I0 (ry) dy = 1+ e− − . 2 2 2 2 0 R r √R r Z − − Therefore, we have 3/2 0 0 s¯ π p′ p 1 R R √R2 r2 (3.17) B J − + + e− − . ≤ 2 2 √R2 r2 R2 r2 (R2 r2)3/2 − − − We can also find the useful estimates of R2 r2 in [30] as − 2 2 2 2 2 g¯ +4 g¯ 1 2 (3.18) R r = p p′ max +1, p p′ , − | − | 4¯g2 ≥ 4 4| − | which lead to

0 0 ′ ′ p′ p √R2 r2 1 (p 0 p0) 1 p p J − e− − e− 2 − e− 2 | − | 1. 2 ≤ ≤ THE RELATIVISTIC QUANTUM BOLTZMANN EQUATION NEAR EQUILIBRIUM 27

In the last inequality, we used

0 0 2 2 ( p′ p )( p′ + p ) p′ p = 1+ p 1+ p = | | | − | | | | | | | p′ p . ′ 2 2 | − | | | − | | 1+ p′ + 1+ p ≤ | − | p p | | | | The estimates (3.18) also implies that p p 1 4 R2 r2 1, . − ≥ R2 r2 ≤ s¯ − Thus the three terms inside thep large bracket of (3.17) are bounded as follows: 1 R 2R p0 + p 0 + 4 ′ . √R2 r2 (R2 r2)3/2 ≤ R2 r2 ≤ s¯ − − − Combining the estimates above yields 1/2 0 0 (3.19) B Cs¯ (p + p′ ). ≤ Now substituting it in (3.6), we turn back to estimate of Γ2,1 as

dp dp′ 1 0 0 0 Γ (f,h), η 2 C s¯2 (p + p′ )J(p′ /2) f(p) h(p′) η(p) . 2,1 Lp 0 0 h i ≤ R3 p R3 p | || || | Z Z ′ 0 0 We uses ¯ 4p p′ to have ≤ 0 0 p p′ 0 Γ (f,h), η 2 C dp dp′ + J(p′ /2) f(p) h(p′) η(p) . 2,1 Lp 0 0 h i ≤ R3 R3 p p | || || | Z Z s ′ s !

0 0 0 Using the boundedness 1 p and p J(p′ /4) C yields ≤ ′ ≤ 0 1 0 2 p 2 Γ2,1(f,h), η Lp C dp dp′(p ) J(p′ /4) f(p) h(p′) η(p) . h i ≤ R3 R3 | || || | Z Z

Finally, we apply the H¨older inequality and obtain 0 1 2 2 2 Γ2,1(f,h), η Lp C h Lp dp (p ) f(p) η(p) h i ≤ k k R3 | || | Z C f ν h 2 η ν. ≤ k k k kLp k k

Regarding the linear terms K1f and K2f in (2.2), we have the following estimates.

Lemma 3.6. For i =1, 2, the compact operator Ki satisﬁes the following estimate.

Kif,h 2 C f ν h ν. h iLp ≤ k k k k Proof. Recall the deﬁnition of K1 in (2.2). Using the energy conservation law, we can replace post-collisional variables by pre-collisional variables on the exponential term:

0 0 K1f C dq dw vøσ(g,θ)J(p /2)J(q /2) f(q) . ≤ R3 S2 | | Z Z Then by v 1 and g (p0)1/2(q0)1/2 and the H¨older inequality, we have ø ≤ ≤ K1f,h 2 C f 2 h 2 . h iLp ≤ k kLp k kLp For K , by the simple boundedness m + τm2(p) CJ(p0/2), we have 2 ≤ p 0 0 K2f C dq dw vøσ(g,θ)J(q /2)J(q′ /2) f(p′) . ≤ R3 S2 | | Z Z 28 GI-CHANBAE,JINWOOJANG,ANDSEOK-BAEYUN

Then we take the inner product with h and write it as

0 0 2 K2f,h Lp C dp dq dw vøσ(g,θ)J(q /2)J(q′ /2) f(p′) h(p) . h i ≤ R3 R3 S2 | || | Z Z Z 0 If we substitute h(q′)= J(q′ /2) in Γ6 in (3.4), then we get

0 0 2 Γ6(f,J),h Lp C dp dq dw vøσ(g,θ)J(q /2)J(q′ /2) f(p′) h(p) . h i ≤ R3 R3 S2 | || | Z Z Z

Thus, the estimates of Γ6(f,h) gives

K2f,h 2 C f ν h ν. h iLp ≤ k k k k

3.2. Estimates of the third-order nonlinear terms. In this subsection, we provide the upper-bound estimates for the third-order nonlinear terms. Lemma 3.7. We have

T (f,h,η), ξ 2 C f 2 h 2 η ν + f 2 h ν η 2 + f ν h 2 η 2 ξ ν . h iLp ≤ k kLp k kLp k k k kLp k k k kLp k k k kLp k kLp k k Proof. Recall the deﬁnition of third-order nonlinear terms in (2.14). Applying the change of variables φ π φ and θ π + θ, we have T (f,h,η)= T (f,h,η). Thus we only consider → − → 1 2 the proof for T1,T3 and T4. (1) Estimates of T : Using the boundedness m + τm2(p) CJ(p0/2), we have 1 ≤ p 0 0 T1(f,h,η) C dq dw vøσ(g,θ)J(p′ /2)J(q /2) f(p) h(q) η(p′) . | |≤ R3 S2 | || || | Z Z We take the inner product with ξ to have

T1(f,h,η), ξ 2 h iLp 0 0 C dp dq dw vøσ(g,θ)J(p′ /2)J(q /2) f(p) h(q) η(p′) ξ(p) . ≤ R3 R3 S2 | || || || | Z Z Z Using the H¨older inequality, we disunite the integrand as follows: 1 2 0 0 2 2 2 T1(f,h,η), ξ Lp C dp dq dw vøσ(g,θ)J(p′ /2)J(q /2) f(p) h(q) h i ≤ R3 R3 S2 | | | | Z Z Z 1 2 0 0 2 2 dp dq dw vøσ(g,θ)J(p′ /2)J(q /2) η(p′) ξ(p) × R3 R3 S2 | | | | Z Z Z = T T . 11 × 12 We ﬁrst consider T . Using v C and g 2 p0q0, we have 11 ø ≤ ≤ 2 0 1 p0 1 0 0 2 2 (T11) dp dq dw (p ) 2 (q ) 2 J(p′ /2)J(q /2) h(q) f(p) ≤ R3 R3 S2 | | | | Z Z Z 2 0 1 2 C dq h(q) dp(p ) 2 f(p) ≤ R3 | | R3 | | Z 2 2 Z C f ν h 2 , ≤ k k k kLp 0 1/2 0 2 where we used (q ) J(q /2) C. To estimate T12, we rewrite (T12) as in the form of (3.6): ≤ THE RELATIVISTIC QUANTUM BOLTZMANN EQUATION NEAR EQUILIBRIUM 29

2 dp dq dp′ dq′ 0 0 (T12) = 0 0 0 0 sσ(g,θ)J(p′ /2)J(q /2) R3 p R3 q R3 p R3 q Z Z Z ′ Z ′ 2 2 (4) µ µ µ µ η(p′) ξ(p) δ (p + q q′ q′ ). × | | | | − − Then we can ﬁnd that the dqdq′ integral has a similar form as of (3.7):

dq dq′ 0 (4) µ µ µ µ B′ = 0 0 sσ(g,θ)J(q /2)δ (p + q q′ q′ ). R3 q R3 q − − Z Z ′ The only diﬀerence is that J(q0) is modiﬁed by J(q0/2). Thus the exponential growth 0 0 0 0 J((p′ p )/2) in (3.17) is changed by J((p′ p )/4). But since the R and r in (3.15) are also changed− by R/2 and r/2, we can have −

0 0 ′ ′ p′ p √R2 r2/2 1 (p 0 p0) 1 p p J − e− − e− 4 − e− 4 | − | 1. 4 ≤ ≤ Thus we can apply the estimate in (3.19): 1/2 0 0 B Cs¯ (p + p′ ), ≤ which gives

2 dp dp′ 1/2 0 0 0 2 2 (3.20) (T12) C 0 0 s¯ (p + p′ )J(p′ /2) η(p′) ξ(p) . ≤ R3 p R3 p | | | | Z Z ′ 0 0 0 1/2 0 Usings ¯ Cp p′ and (p′ ) J(p′ /2) C, we have ≤ ≤ 2 0 1/2 2 2 2 2 (T12) C dp(p ) ξ(p) dp′ η(p′) C η L2 ξ ν. ≤ R3 | | R3 | | ≤ k k p k k Z Z Combining the estimates of T11 and T12 yields

T1(f,h,η), ξ 2 C f ν h 2 η 2 ξ ν . h iLp ≤ k k k kLp k kLp k k (2) Estimates of T : Using m + τm 2(p) J(p0/2), we estimate T as follows: 3 ≤ 3 p 0 0 T3(f,h,η) dq dw vøσ(g,θ)J(p′ /2)J(q′ /2) f(p) h(p′) η(q′) . | |≤ R3 S2 | || || | Z Z Via the H¨older inequality, we can split the variable as follows: 1 2 0 0 2 2 2 T3(f,h,η), ξ Lp C dp dq dw vøσ(g,θ)J(p′ /2)J(q′ /2) f(p) h(p′) h i ≤ R3 R3 S2 | | | | Z Z Z 1 2 0 0 2 2 dp dq dw vøσ(g,θ)J(p′ /2)J(q′ /2) ξ(p) η(q′) × R3 R3 S2 | | | | Z Z Z = T T . 31 × 32 Via the change of variables p′ q′, we can see that T31 and T32 are indeed identical. Thus ↔ 2 we only consider T31 part here. Similarly to the estimates of Γ2, we rewrite (T31) as follows:

dp dq dp′ dq′ 0 0 2 2 (4) µ µ µ µ 0 0 0 0 sσ(g,θ)J(p′ /2)J(q′ /2) f(p) h(p′) δ (p +q q′ q′ ). R3 p R3 q R3 p R3 q | | | | − − Z Z Z ′ Z ′ Then we note that the energy conservation law and the exponential decays with respect to the two post-collisional momenta can together imply the decays with respect to all pre-post collisional variables: 0 0 0 0 0 0 J(p′ /2)J(q′ /2) = J(p /4)J(q /4)J(p′ /4)J(q′ /4). 30 GI-CHANBAE,JINWOOJANG,ANDSEOK-BAEYUN

Then, we use the boundedness of B in (3.19) as

dq dq′ 0 (4) µ µ µ µ 1/2 0 0 B = 0 0 sσ(g,θ)J(q /4)δ (p + q q′ q′ ) Cs¯ (p + p′ ), R3 q R3 q − − ≤ Z Z ′ which yields

2 dp dp′ 1/2 0 0 0 0 2 2 (T31) C 0 0 s¯ (p + p′ )J(p /4)J(p′ /4) f(p) h(p′) . ≤ R3 p R3 p | | | | Z Z ′ 0 0 0 1/2 0 Usings ¯ Cp p′ and (p ) J(p /4) C, we have ≤ ≤ 2 2 2 2 2 (T31) C dp f(p) dp′ h(p′) C f L2 h L2 . ≤ R3 | | R3 | | ≤ k k p k k p Z Z Since T32 has the same form with T31, we have

T3(f,h,η), ξ 2 C f 2 h 2 η 2 ξ 2 . h iLp ≤ k kLp k kLp k kLp k kLp (3) Estimates of T : Recall the deﬁnition of T in (2.14). Using Lemma 3.1, we can 4 4 have CJ(p0/2) m + τm2(p) CJ(p0/2), which yields ≤ ≤ 0 0 0 p J(q /2)J(p′ /2)J(q′ /2) T4(f,h,η) C dq dw vøσ(g,θ) 0 f(q) h(p′) η(q′) . | |≤ R3 S2 J(p /2) | || || | Z Z 0 0 0 0 Using the energy conservation law J(p′ /2)J(q′ /2) = J(p /2)J(q /2), we have

0 T4(f,h,η) dq dw vøσ(g,θ)J(q ) f(q) h(p′) η(q′) . | |≤ R3 S2 | || || | Z Z Then we take the inner product with ξ to have

0 2 T4(f,h,η), ξ Lp C dp dq dw vøσ(g,θ)J(q ) f(q) h(p′) η(q′) ξ(p) . h i ≤ R3 R3 S2 | || || || | Z Z Z

Via the H¨older inequality, we divide the variable into two pre-collisional variables and two post-collisional variables as follows:

1 2 0 2 2 2 T4(f,h,η), ξ Lp C dp dq dw vøσ(g,θ)J(q ) f(q) ξ(p) h i ≤ R3 R3 S2 | | | | Z Z Z 1 2 0 2 2 dp dq dw vøσ(g,θ)J(q ) h(p′) η(q′) × R3 R3 S2 | | | | Z Z Z = T T . 41 × 42 For the estimates of T , we use v σ(g,θ) C(p0)1/2(q0)1/2 to have 41 ø ≤ 1 1 2 0 2 0 2 0 2 2 2 2 (T41) C dp dq dw (p ) (q ) J(q ) f(q) ξ(p) C f L2 ξ ν, ≤ R3 R3 S2 | | | | ≤ k k p k k Z Z Z 0 1 0 where we used (q ) 2 J(q ) C. Similarly, we have ≤ 2 0 1 2 2 (T42) C dp dq dw (p ) 2 h(p′) η(q′) . ≤ R3 R3 S2 | | | | Z Z Z 0 1/2 0 1/2 0 1/2 We use (p ) (p′ ) + (q′ ) to have ≤ 2 0 1 0 1 2 2 (T42) C dp dq dw (p′ ) 2 + (q′ ) 2 h(p′) η(q′) . ≤ R3 R3 S2 | | | | Z Z Z THE RELATIVISTIC QUANTUM BOLTZMANN EQUATION NEAR EQUILIBRIUM 31

Then applying the pre-post collisional change of variables (p, q) (p′, q′) as in (3.5) yields ↔ T42 C h 2 η ν + h ν η 2 . ≤ k kLp k k k k k kLp Combining the estimates of T41 andT42 yields

T4(f,h,η), ξ 2 C h 2 η ν + h ν η 2 f 2 ξ ν . h iLp ≤ k kLp k k k k k kLp k kLp k k

4. Local existence In this section, we construct the local-in-time classical solution. Here we briefly introduce a main difference of the quantum Boltzmann theory from the Newtonian one. From the statistical description of quantum Boltzmann equation in [17], we denote that the term (1 F (p)) in the collision operator for fermions is the probability that a fermion is being placed− in the p momentum place after a collision. Thus the ratio (1 F (p)) has to be non- negative. Then, in each iteration scheme in the proof of the local− existence, the collision operator depending on F n+1 is well-defined only when F n+1 is in the interval [0, 1] in the case of fermions. Thus we need to prove the boundedness of the solution F n+1 in each iteration step as well. Theorem 4.1. Let N 3. Suppose that the initial data F satisfies ≥ 0 0 F (x, p)= m(p)+ m(p) m2(p)f (x, p) 1, for fermion ≤ 0 − 0 ≤ 0 F (x, p)= m(p)+ m(p)+ m2(p)f (x, p), for boson, ≤ 0 p 0 M0 Then there exist M0 > 0 and T >p0 such that if (f0) 2 then there exists a unique local-in-time solution f(x,p,t) of ∗(2.15) satisfying E ≤ (1) The higher order energy is uniformly bounded:

sup (f(t)) M0. 0 t T∗ E ≤ ≤ ≤ (2) The distribution function is bounded in t [0,T ]: ∈ ∗ 0 F (x,p,t)= m(p)+ m(p) m2(p)f(x,p,t) 1, for fermions ≤ − ≤ 0 F (x,p,t)= m(p)+ m(p)+ m2(p)f(x,p,t), for bosons, ≤ p (3) The higher order energy normp (f(t)) is continuous in t [0,T ]. E ∈ ∗ Proof. We first take an iteration scheme as follows: (4.1) (∂ +ˆp )F n+1 = Q(F n, F n, F n+1, F n), t · ∇x 0 n+1 with F (x,p,t) = F0(x, p). Note that F in the operator Q is placed in the p variable location. We proceed the proof by applying the induction argument. Let us assume that n 0 F (x,p,t) 1 for fermions, n (4.2) ≤ n ≤ and sup (f (t)) M0. 0 F (x,p,t) forbosons, 0 t T∗ E ≤ ≤ ≤ ≤ Note that the collision operator Q(F n, F n, F n+1, F n) is well-defined when F n+1 is in the interval [0, 1]. We define

G(F )= dq dw vøσ(g,θ)F (p′)F (q′)(1 + τF (q)), R3 S2 Z Z R(F )= dq dw vøσ(g,θ)(1 + τF (p′))(1 + τF (q′))F (q), R3 S2 Z Z 32 GI-CHANBAE,JINWOOJANG,ANDSEOK-BAEYUN to have Q(F n, F n, F n+1, F n)= G(F n)(1 + τF n+1(p)) R(F n)F n+1(p). − Then we rewrite (4.1) as follows: (4.3) (∂ +ˆp τG(F n)+ R(F n))F n+1 = G(F n). t · ∇x − In the case of fermions, we additionally consider the upper bound of F n+1. We observe that the induction hypothesis (4.2)1 implies 0 G(F n), 0 R(F n). ≤ ≤ Since R(F n) is positive, we have (∂ +ˆp + G(F n)+ R(F n))F n+1 G(F n)+ R(F n). t · ∇x ≤ The associated ODE for the particle characteristic trajectory is given by dX(s)/ds =p ˆ(s) where X(s)= X(s; t,x,p). We integrate over the particle path to obtain

n+1 R t(G(F n)+R(F n))dτ F (X(t),p,t) e− 0 F (X(0),p) ≤ 0 t R t(G(F n)+R(F n))dτ n n + e− s (G(F )+ R(F ))(X(s),p,s)ds. Z0 We observe from

d R t(G(F n)+R(F n))dτ R t(G(F n)+R(F n))dτ n n e− s = e− s (G(F )+ R(F ))(X(s),p,s), ds that n o n+1 R t(G(F n)+R(F n))dτ R t(G(F n)+R(F n))dτ F (x,p,t) e− 0 F (x pt,pˆ )+1 e− 0 ≤ 0 − − R t(G(F n)+R(F n))dτ = e− 0 (F (x pt,pˆ ) 1)+1 0 − − 1, ≤ where we used the boundedness of the initial data F0 1. We now consider the lower bound of F n+1. By G(F≤n) 0 on (4.3), we have ≥ (∂ +ˆp τG(F n)+ R(F n))F n+1 0, t · ∇x − ≥ which yields n+1 R t( τG(F n)+R(F n))dt F (x,p,t) e− 0 − F (x pt,pˆ ) 0, ≥ 0 − ≥ where we used the boundedness of the initial data F0 0. Now the collision operator of (4.1) is well-deﬁned. ≥ Substituting F n+1 = m + √m + τm2f n+1 on (4.1) gives (∂ +ˆp + ν)f n+1 = Kf n + Γ(f n,f n+1), t · ∇x where K = K K and 2 − 1 n n+1 n+1 n n n Γ(f ,f )= Γi(f ,f )+ Γi(f ,f ) i=1,2,4 i=3,5,6 X X n+1 n n n n n + Ti(f ,f ,f )+ T4(f ,f ,f ). i=1,2,3 X We take ∂α on each side to have ∂ ∂αf n+1 +ˆp ∂αf n+1 + ∂α(νf n+1)= ∂αKf n + ∂αΓ(f n,f n+1). t · ∇x THE RELATIVISTIC QUANTUM BOLTZMANN EQUATION NEAR EQUILIBRIUM 33

2 α n+1 We take the Lx,p inner product with ∂ f , then the nonlinear estimates Lemma 3.4 and Lemma 3.7 with the induction hypothesis (4.2)2 yield

2 2 3 M0 (1 C M0 CM0 CT M0 CM0 CM0 ) sup n+1(t) − − − ∗ − − 0 t T∗ E ≤ 2 p ≤ ≤ For suﬃciently small T and M0, we have ∗ sup n+1(t) M0. 0 t T∗ E ≤ ≤ ≤ Taking the limit as n gives a local-in-time classical solution. The remaining proof is standard as in [38, 39]→ and ∞ we omit it.

5. Global existence In this section, we extend the local solution constructed in Theorem 4.1 to a global solution. For this we ﬁrst recover the full coercivity estimates of the linear operator L. 5.1. Coercivity estimate. Recall the deﬁnition of Pf in Lemma 2.2. Since Pf is the 2 orthonormal projection to the Lp space with the following basis,

2 2 2 2 0 2 m + τm ,p1 m + τm ,p2 m + τm ,p3 m + τm ,p m + τm , Pf can benp written as followsp by the Gram-Schmidtp process:p p o Pf = m + τm2 + p m + τm2 + p0 m + τm2, A B· C where , and are givenp by p p A B C 1 2 λ0 1 0 2 λ0 2 = dp f m + τm 2 dp fp m + τm dp f m + τm , λ R3 λ λ0 R3 λ R3 A − λ00 − Z p − λ Z p Z p 1 i 2 i = dp fp m + τm , B λi R3 Z p 1 0 2 λ0 2 = 2 dp fp m + τm dp f m + τm , λ0 R3 λ R3 C λ00 − − λ Z p Z p and

2 i 2 2 λ = dp m(p)+ τm (p), λi = dp (p ) (m(p)+ τm (p)), R3 R3 Z Z 0 2 0 2 2 λ0 = dp p (m(p)+ τm (p)), λ00 = dp (p ) (m(p)+ τm (p)), R3 R3 Z Z for i =1, 2, 3. Since Pf has the exponential decay √m + τm2, we can have α α α α (5.1) ∂ Pf x,ν ∂ 2 + ∂ 2 + ∂ 2 . k k ≤ k AkLx k BkLx k CkLx α N α N | X|≤ | X|≤ Now we substitute f = Pf + (I P )f in the perturbation equation (2.15) to have − (∂ +ˆp )(Pf)= (∂ +ˆp )((I P )f) L(I P )f + Γ(f)+ T (f), t · ∇x − t · ∇x − − − where we used L(Pf) = 0 by Lemma 2.2. We expand the left-hand side as follows:

3 p 3 p p ∂ + ∂ i + (∂ + ∂ )p + ∂ i j + ∂ p0 m + τm2,  tA xi Ap0 tBi xi C i xi Bj p0 tC  i=1 i=1 1 i,j 3 X X ≤X≤ p   34 GI-CHANBAE,JINWOOJANG,ANDSEOK-BAEYUN which is a linear combination of following 14-basis: (5.2) p p p m + τm2, i m + τm2, p m + τm2, i j m + τm2, p0 m + τm2 , p0 i p0 for 1p i, j 3. To denotep the right-handp side, we define p p ≤ ≤ (5.3) l = (∂ +ˆp + L)((I P )f), h = Γ(f)+ T (f). − t · ∇x − By expanding l and h with respect to the 14-basis elements in (5.2), we obtain the following macro-micro system: ∂ = l + h , tA a a ∂ = l + h , xi A i i (5.4) ∂t i + ∂xi = lbci + hbci, B C ∂ + ∂ = l + h , xi Bj xj Bi ij ij ∂ = l + h , tC c c for i, j = 1, , 3. On the right-hand side, la,li,lbci,lij ,lc and ha,hi,hbci,hij ,hc are the coefficients of··· expansion of l and h with respect to the basis in (5.2), respectively. We define the summation of the coefficients of l and h as follows:

˜l = la + lc + (li + lbci)+ lij , 1 i 3 1 i,j 3 ≤X≤ ≤X≤ h˜ = ha + hc + (hi + hbci)+ hij . 1 i 3 1 i,j 3 ≤X≤ ≤X≤ Note that the macro-micro system of (5.4) is identical to the macro-micro system of the relativistic Landau-Maxwell model (98)-(102) in [75] when the electromagnetic ﬁeld is zero B = E = 0. However, we provide the full proof here for the readers’ convenience. We ﬁrst observe following property of , and . The conservation laws of f in (2.16) can be written in the following form. A B C Lemma 5.1. We have (x, t)dx = (x, t)dx = (x, t)dx =0. T3 A T3 B T3 C Z Z Z Proof. The conservation laws in (1.15) implies 1 dxdp (F F0) µ =0. T3 R3 − p Z × Combining with the assumption (1.20), we have

2 1 dxdp f m + τm µ =0. T3 R3 p Z × p We derived desired results. We now establish the estimates of , , . The proof of the following lemma is motivated by the proof of the estimate of (108)A inB [75].C Lemma 5.2. Let N 3. We have ≥ α α α α α ∂ 2 + ∂ 2 + ∂ 2 C ∂ l˜ 2 + ∂ h˜ 2 . k AkLx k BkLx k CkLx ≤ k kLx k kLx α N α N 1 | X|≤ | |≤X− THE RELATIVISTIC QUANTUM BOLTZMANN EQUATION NEAR EQUILIBRIUM 35

Proof. We first prove the estimate of . Taking ∂α on the first equation of (5.4) gives A ∂α∂ = ∂αl + ∂αh . tA a a we take the innerproduct with ∂α∂ and apply the Hölder inequality to have tA α 2 α α α α ∂ ∂t 2 C ∂ la L2 ∂ ∂t L2 + ∂ ha L2 ∂ ∂t L2 . k AkLx ≤ k k x k Ak x k k x k Ak x α Dividing each side by ∂ ∂t 2 , we have k AkLx α α α ∂ ∂t 2 C( ∂ la 2 + ∂ ha 2 ). k AkLx ≤ k kLx k kLx Similarly we take ∂α on the second equation of (5.4) to have ∂α∂ = ∂αl + ∂αh . xi A i i Multiplying ∂α∂ on each side and applying the Hölder inequality, we have xi A α α α ∂ x 2 C( ∂ li 2 + ∂ hi 2 ). k ∇ AkLx ≤ k kLx k kLx We apply the Poincaréinequality on to have A 1 2 (x)dx C 2 . Lx 3 x Lx kAk − T T3 A 2 ≤ k∇ Ak Z Lx | | Then the conservation law for in Lemma 5.1 implies A 2 C x 2 . kAkLx ≤ k∇ AkLx This completes the estimate of . Since and are connected by the third equation of (5.4), we consider the fourth andA the fifth equationsB C of (5.4) first. For the notational brevity, we use ∂j to denote ∂xj . Using the fourth equation of (5.4), we calculate = ∂ = ∂ + ∂ △Bi jj Bi jj Bi iiBi 1 j 3 j=i ≤X≤ X6 1 = (∂ l + ∂ h ∂ )+ (∂ l + ∂ h ). j ij j ij − jiBj 2 i ii i ii j=i X6 We substitute the fourth equation for i = j case to have 1 1 = ∂ l + ∂ h (∂ l + ∂ h ) + (∂ l + ∂ h ). △Bi j ij j ij − 2 i jj i jj 2 i ii i ii j=i X6 Taking ∂α and multiplying ∂α on each side yield Bi α α α (5.5) ∂ x 2 C ∂ lij 2 + ∂ hij 2 . k ∇ BkLx ≤ k kLx k kLx 1 i,j 3 ≤X≤ We take ∂α on the fifth equation of (5.4) to have α α α (5.6) ∂ ∂t 2 C ∂ lc 2 + ∂ hc 2 . k CkLx ≤ k kLx k kLx Now we use the third equation of (5.4). We consider the estimate of when there are more B n than one temporal derivatives. For this, we take the temporal derivative ∂t for 1 n N 1 on the third equation of (5.4) to have ≤ ≤ − ∂n∂ = ∂nl + ∂nh ∂n∂ . t tBi t bci t bci − t xi C Then the estimate (5.6) gives n+1 n n n 1 n 1 (5.7) ∂ i 2 ∂ lbci 2 + ∂ hbci 2 + C ∂x ∂ − lc 2 + ∂x ∂ − hc 2 . k t B kLx ≤k t kLx k t kLx k i t kLx k i t kLx 36 GI-CHANBAE,JINWOOJANG,ANDSEOK-BAEYUN

For the estimate of when there are less than two temporal derivatives, we apply the Poincar´einequality onB and ∂ to have B tB n 1 n n ∂ 2 ∂ dx ∂ 2 , t i Lx 3 t i x t i Lx k B k − T T3 B ≤ k∇ B k Z | | for n =0, 1. Then the conservation law in Lemma 5.1 yields

n n ∂ i 2 x∂ i 2 . k t B kLx ≤ k∇ t B kLx Combining with (5.5) and (5.7), we conclude that

α α α (5.8) ∂ 2 C ∂ l˜ 2 + ∂ h˜ 2 . k BkLx ≤ k kLx k kLx α N α N 1 | X|≤ | |≤X− For the estimate of , we take ∂α on third equation of (5.4) to have C ∂α∂ = ∂αl + ∂αh ∂α∂ . xi C bci bci − tBi Applying (5.8) yields

α β β (5.9) ∂ ∂x 2 C ∂ l˜ 2 + ∂ h˜ 2 . k i CkLx ≤ k kLx k kLx β = α | X| | | Then the Poincar´einequality with Lemma 5.1 gives

2 C x 2 . kCkLx ≤ k∇ CkLx Combining with the estimates (5.6) and (5.9), we derive the desired result. Lemma 5.3. Let N 3. Suppose that ≥ α 2 ∂ f 2 M0. k kLx,p ≤ α N | X|≤ Then we have α α (1) ∂ ˜l 2 C (I P )∂ f x,ν, k kLx ≤ k − k α N 1 α N | |≤X− | X|≤ α α (2) ∂ h˜ 2 C( M0 + M0) ∂ f x,ν. k kLx ≤ k k α N α N | X|≤ p | X|≤ Proof. For the Newtonian Boltzmann equation and the relativistic Landau-Maxwell system, analogous proofs can be found in [38] and [75], respectively. The diﬀerence now is that the estimate of the nonlinear term h includes the third-order nonlinear terms. (1) We denote the 14-basis of (5.2) as ei 1 i 14, and let ei∗ 1 i 14 be the corresponding orthonormal basis. Then the orthonormal{ } basis≤ ≤ can be written{ } ≤ by≤ a linear combination of the original basis ei 1 i 14 as follows: { } ≤ ≤ 14

ei∗ = Cij ej , j=1 X for j =1, , 14. We consider the orthonormal expansion of l as follows: ··· 14 14 14 14 l = l,e∗ 2 e∗ = l, Cij ej Cikek. h i iLp i i=1 i=1 * j=1 + 2 k=1 X X X Lp X THE RELATIVISTIC QUANTUM BOLTZMANN EQUATION NEAR EQUILIBRIUM 37

Then the coeﬃcient of ek can be read as follows:

Cij Cik l,ej 2 , h iLp 1 i,j 14 ≤X≤ which correspond to la,li,lbci,lij , and lc. By the deﬁnition of l in (5.3) and the linear operator L in Proposition 2.1, we can write l as l = (∂ +ˆp + ν + K K )((I P )f). − t · ∇x 1 − 2 − For α = N 1, we have | | − α dp ∂ l ei(p) R3 · 2 Z Lx

β α dp (1+p ˆ)(I P )∂ f + (ν + K1 K2)(I P )∂ f ei(p) . ≤ R3 − − − · 2 Z β =N Lx | X|

We use the H¨older inequality on the (1 +p ˆ)(I P )∂βf term and the ν(I P )∂αf term. Then we have − −

β α α 2 dp (1+p ˆ)(I P )∂ f + ν(I P )∂ f ei(p) (I P )∂ f Lx,p , R3 − − · ≤ k − k Z β =N α N | X| | X|≤ 2 where we used the exponential decay √m + τm of the basis ei(p) in (5.2) that implies

2 2 (1+p ˆ) + ν(p) ei (p)dp C. R3 ≤ Z From the estimate of the compact operator Kj from Lemma 3.6, we also have

Kj f,ei 2 C f ν, h iLp ≤ k k for j =1, 2 and i =1, , 14. The estimates above together yield ··· α α α 2 ∂ l ei(p)dp C (I P )∂ f Lx,p + (I P )∂ f x,ν R3 · 2 ≤ k − k k − k Z Lx α =N | X| α C (I P )∂ f x,ν. ≤ k − k α N | X|≤ (2) In the same manner, the coeﬃcient of ek of expansion of l can be written as

Cij Cik h,ei 2 . h iLp 1 i,j 14 ≤X≤ For the second-order nonlinear terms, by Lemma 3.4, we have

α α1 α2 α1 ∂ Γ(f,f),ei 2 C ∂ f 2 ∂ f x,ν C M0 ∂ f x,ν. h iLx,p ≤ k kLx,p k k ≤ k k α1 + α2 α α1 α | | X| |≤| | p | X|≤| | For the third-order nonlinear term, Lemma 3.7 yields

α α1 α2 α3 ∂ T (f,f,f),ei C ∂ f 2 ∂ f 2 ∂ f ν. h i ≤ k kLp k kLp k k α1 + α2 + α3 α | | | X| | |≤| | 2 3 3 Then the Sobolev embedding H (T ) L∞(T ) implies ⊂⊂

α α1 ∂ T (f,f,f) ei(p)dp CM0 ∂ f x,ν. R3 · L2 ≤ k k Z x α1 α | X|≤| |

38 GI-CHANBAE,JINWOOJANG,ANDSEOK-BAEYUN

So we obtain the desired results. We now have all the estimates to recover the full coercivity. We combine (5.1) with Lemma 5.2 and Lemma 5.3 to have α α α ∂ Pf x,ν ∂ l 2 + ∂ h 2 k k ≤ k kLx k kLx α N α N 1 | X|≤ | |≤X− C (I P )∂αf + C M ∂αf . ≤ k − kx,ν 0k kx,ν α N | X|≤ p Thus, for a suﬃciently small M0, there exists δ > 0 such that α α α 2 (5.10) L∂ f, ∂ f 2 δ ∂ f , h iLx,p ≥ k kx,ν α N α N | X|≤ | X|≤ by Lemma 3.3. 5.2. Global existence. We now have all the ingredients for the proof of Theorem 1.1. Extending the local solution constructed in Theorem 4.1 to a global solution is standard as in [37, 38]. We only sketch the proof here. Recall that substituting F = m + √m + τm2f on (1.13) yields the following linearized equation: ∂ f +p ˆ f + Lf = Γ(f)+ T (f), t · ∇x f(x, p, 0) = f0(x, p). We take ∂α on each side to have ∂ ∂αf +p ˆ ∂αf + L∂αf = ∂αΓ(f)+ ∂αT (f). t · ∇x Taking the inner product with ∂αf yields

1 d α 2 α α α α α α ∂ f 2 + L∂ f, ∂ f L2 = ∂ f, ∂ Γ(f) L2 + ∂ f, ∂ T (f) L2 . 2 dtk kLx,p h i x,p h i x,p h i x,p Applying the full coercivity estimate (5.10), we have

1 d α 2 α 2 α α α α ∂ f 2 + δ ∂ f ∂ f, ∂ Γ(f) L2 + ∂ f, ∂ T (f) L2 . 2 dtk kLx,p k kx,ν ≤ h i x,p h i x,p The estimate of the second-order nonlinear terms in Lemma 3.4 gives

α α α1 α2 α1 α2 α 2 2 2 ∂ f, ∂ Γ(f) Lx,p C dx ∂ f Lp ∂ f ν + ∂ f ν ∂ f Lp ∂ f ν . h i ≤ T3 k k k k k k k k k k α1 + α2 N Z | | X| |≤ Without loss of generality, we assume that α1 is less than or equal to α2. Since N 3, we 2 3 3 ≥ have α +2 N. Thus the Sobolev embedding H (T ) L∞(T ) implies | 1| ≤ ⊂⊂ α α α α α ∂ f, ∂ Γ(f) 2 C ∂ f 2 ∂ f x,ν ∂ f x,ν h iLx,p ≤ k kLx,p k k k k α N α N α N α N | X|≤ | X|≤ | X|≤ | X|≤ C (t) ∂αf 2 . ≤ E k kx,ν α N p | X|≤ For the estimate of the third-order nonlinear terms, we apply Lemma 3.7 to have

α α α1 α2 α3 2 2 2 ∂ f, ∂ T (f) Lx,p C dx ∂ f Lp ∂ f Lp ∂ f ν h i ≤ T3 k k k k k k α1 + α2 + α3 N Z | | | X| | |≤ α1 α2 α3 α1 α2 α3 α + ∂ f 2 ∂ f ν ∂ f 2 + ∂ f ν ∂ f 2 ∂ f 2 ∂ f ν. k kLp k k k kLp k k k kLp k kLp k k THE RELATIVISTIC QUANTUM BOLTZMANN EQUATION NEAR EQUILIBRIUM 39

Similarly we assume that α1 and α2 are less than or equal to α3. Since N 3, we have 2 3 ≥3 α +2 N and α +2 N. By the Sobolev embedding H (T ) L∞(T ), we have | 1| ≤ | 2| ≤ ⊂⊂ α α α α α α ∂ f, ∂ T (f) 2 C ∂ f 2 ∂ f 2 ∂ f x,ν ∂ f x,ν h iLx,p ≤ k kLx,p k kLx,p k k k k α N α N α N α N α N | X|≤ | X|≤ | X|≤ | X|≤ | X|≤ C (t) ∂αf 2 . ≤ E k kx,ν α N | X|≤ Combining these estimates, we conclude that

1 d α 2 α 2 α 2 ∂ f 2 + δ ∂ f x,ν C (t)+ (t) ∂ f x,ν. 2 dtk kLx,p k k ≤ E E k k α N α N | X|≤ p | X|≤ The remaining proof can be established by the standard continuity argument as in [37, 38]. This completes the proof. Acknowledgement: J. W. Jang is supported by CRC 1060 The mathematics of emergent eﬀects at the University of Bonn funded through the German Science Foundation (DFG). S.-B. Yun is supported by Samsung Science and Technology Foundation under Project Number SSTF-BA1801-02.

References

[1] Akama, H.: Relativistic Boltzmann equation for plasmas. Journal of the Physical Society of Japan, 28 (1970), no. 2, 478-488. [2] Alexandre, R., Villani, C.: On the Boltzmann equation for long-range interactions. Comm. Pure Appl. Math. 55 (2002), no. 1, 30-70. [3] Andréasson, H., Calogero, S., Illner, R.: On blowup for gain-term-only classical and relativistic Boltz- mann equations. Math. Methods Appl. Sci. 27 (2004), no. 18, 2231-2240. [4] Arkeryd, L., Esposito, R., Marra, R. and Nouri, A.: Exponential stability of the solutions to the Boltzmann equation for the Benard problem. Kinet. Relat. Models 5 (2012), no. 4, 673-695. [5] Bae, G.-C., Yun, S.-B.: Quantum BGK model near a global Fermi-Dirac distribution. SIAM J. Math. Anal. 52 (2020), no. 3, 2313-2352. [6] Benedetto, D., Castella, F., Esposito, R. and Pulvirenti, M.: A short review on the derivation of the nonlinear quantum Boltzmann equations. Commun. Math. Sci. 5 (2007), 55-71. [7] Benedetto, D., Castella, F., Esposito, R. and Pulvirenti, M.: From the N-body Schrödinger equation to the quantum Boltzmann equation: a term-by-term convergence result in the weak coupling regime. Comm. Math. Phys. 277 (2008), no. 1, 1-44. [8] Bobylev, A. V., Cercignani, C. and Gamba, I. M.: On the self-similar asymptotics for generalized nonlinear kinetic Maxwell models. Comm. Math. Phys. 291 (2009), no. 3, 599-644. [9] Boblylev, A. V., Pulvirenti, M. and Saffirio, C.: From particle systems to the Landau equation: a consistency result. Comm. Math. Phys. 319 (2013), no. 3, 683-702. [10] Briant, M., Einav, A.: On the Cauchy problem for the homogeneous Boltzmann-Nordheim equation for bosons: local existence, uniqueness and creation of moments. J. Statist. Phys. 163 (2016), no. 5, 1108-1156. [11] Buss, O., Gaitanos, T., Gallmeister, K., Van Hees, H., Kaskulov, M., Lalakulich, O., Larionov, A. B., Leitner, T., Weil, J. and Mosel, U.: Transport-theoretical description of nuclear reactions. Physics Reports, 512 (2012), no. 1-2, 1-124. [12] Cercignani, C.: Slow rarefied flows. Theory and application to micro-electro-mechanical systems. Progress in Mathematical Physics, 41. Birkhäuser Verlag, Basel, 2006. [13] Cercignani, C.: The Boltzmann equation and its applications. Applied Mathematical Sciences, 67. Springer-Verlag, New York, 1988. [14] Cercignani, C., Illner, R. and Pulvirenti, M.: The mathematical theory of dilute gases. Applied Math- ematical Sciences, 106. Springer-Verlag, New York, 1994. [15] Cercignani, C., Kremer, G. M.: The relativistic Boltzmann equation: theory and applications. Progress in Mathematical Physics, 22. Birkhäuser Verlag, Basel, 2002. 40 GI-CHANBAE,JINWOOJANG,ANDSEOK-BAEYUN

[16] Chapman, J., Jang, J. W., and Strain, R. M.: On the Determinant Problem for the Relativistic Boltzmann Equation. arXiv preprint arXiv:2006.02540. (2020). [17] Chapman, S., Cowling, T. G.: The mathematical theory of non-uniform gases. An account of the kinetic theory of viscosity, thermal conduction and diffusion in gases. Third edition, prepared in co-operation with D. Burnett Cambridge University Press, London 1970. [18] De Groot, S. R., Van Leeuwen, W. A., Van Weert, C. G.: Relativistic kinetic theory. Principles and Applications. North-Holland Publishing Co., Amsterdam-New York, 1980. [19] Desvillettes, L., Mouhot, C.: Stability and uniqueness for the spatially homogeneous Boltzmann equation with long-range interactions. Arch. Ration. Mech. Anal. 193 (2009), no. 2, 227-253. [20] DiPerna, R. J., Lions, P.-L.: On the Cauchy problem for Boltzmann equations: global existence and weak stability. Ann. of Math. 130 (1989), no. 2, 321-366. [21] Duan, R., Yu, H.: The relativistic Boltzmann equation for soft potentials. Adv. Math. 312 (2017), 315-373. [22] Dudyński, M., Ekiel-Je˙zewska, M. L.: Global existence proof for relativistic Boltzmann equation. J. Statist. Phys. 66 (1992), no. 3-4, 991-1001. [23] Dudyński, M., Ekiel-Je˙zewska, M. L.: On the linearized relativistic Boltzmann equation. I. Existence of solutions. Comm. Math. Phys. 115 (1988), no. 4, 607-629. [24] Dudyński, M., Ekiel-Je˙zewska, M. L.: The relativistic Boltzmann equation-mathematical and physical aspects. J. Tech. Phys, 48 (2007), no. 1, 39-47. [25] Escobedo, M., Mischler, S. and Valle, M. A.: Entropy maximisation problem for quantum relativistic particles. Bull. Soc. Math. France 133 (2005), no. 1, 87-120. [26] Escobedo, M., Mischler, S. and Valle, M. A.: Homogeneous Boltzmann equation in quantum relativistic kinetic theory. Electronic Journal of Differential Equations. Monograph, 4. Southwest Texas State University, San Marcos, TX, 2003. [27] Escobedo, M., Velázquez, J. J. L.: On the blow up and condensation of supercritical solutions of the Nordheim equation for bosons. Comm. Math. Phys. 330 (2014), no. 1, 331-365. [28] Gallagher, I., Saint-Raymond, L. and Texier, B.: From Newton to Boltzmann: hard spheres and short- range potentials. Zurich Lectures in Advanced Mathematics. European Mathematical Society (EMS), Zürich, 2013. [29] Glassey, R. T.: The Cauchy problem in kinetic theory. Society for Industrial and Applied Mathematics (SIAM), Philadelphia, PA, 1996. [30] Glassey, R. T., Strauss, W. A.: Asymptotic stability of the relativistic Maxwellian. Publ. Res. Inst. Math. Sci. 29 (1993), no. 2, 301-347. [31] Glassey, R. T., Strauss, W. A.: Asymptotic stability of the relativistic Maxwellian via fourteen moments. Transport Theory Statist. Phys. 24 (1995), no. 4-5, 657-678. [32] Glassey, R. T., Strauss, W. A.: On the derivatives of the collision map of relativistic particles. Transport Theory Statist. Phys. 20 (1991), no. 1, 55-68. [33] Grad, H.: On the kinetic theory of rarefied gases. Comm. Pure Appl. Math. 2 (1949), 331-407. [34] Grad, H.: Principles of the kinetic theory of gases. Handbuch der Physik (herausgegeben von S. Flügge), Bd. 12, Thermodynamik der Gase pp. 205-294 Springer-Verlag, Berlin-Göttingen-Heidelberg, 1958. [35] Gradshteyn, I. S., Ryzhik, I. M.: Table of integrals, series, and products. Translated from the Russian. Translation edited and with a preface by Daniel Zwillinger and Victor Moll. Eighth edition. Revised from the seventh edition. Elsevier/Academic Press, Amsterdam, 2015. [36] Gressman, P., Strain, R. M.: Global classical solutions of the Boltzmann equation without angular cut-off. J. Amer. Math. Soc. 24 (2011), no. 3, 771-847. [37] Guo, Y.: The Boltzmann equation in the whole space. Indiana Univ. Math. J. 53 (2004), no. 4, 1081- 1094. [38] Guo, Y.: The Vlasov-Maxwell-Boltzmann system near Maxwellians. Invent. Math. 153 (2003), no. 3, 593-630. [39] Guo, Y.: The Vlasov-Poisson-Boltzmann system near Maxwellians. Comm. Pure Appl. Math. 55 (2002), no. 9, 1104-1135. [40] Guo, Y., Strain, R. M.: Momentum regularity and stability of the relativistic Vlasov-Maxwell- Boltzmann system. Comm. Math. Phys. 310 (2012), no. 3, 649-673. [41] Harris, S.: An introduction to the theory of the Boltzmann equation. Courier Corporation. 2004. [42] Ha, S.-Y., Jeong, E. and Strain, R. M.: Uniform L1-stability of the relativistic Boltzmann equation near vacuum. Commun. Pure Appl. Anal. 12 (2013), no. 2, 1141-1161. [43] Huang, K.: Statistical Mechanics, John Wiley and Sons. New York. 1963. THE RELATIVISTIC QUANTUM BOLTZMANN EQUATION NEAR EQUILIBRIUM 41

[44] Jang, J. W.: Global classical solutions to the relativistic Boltzmann equation without angular cut-off, Ph.D. thesis, University of Pennsylvania, 2016. [45] Jang, J. W., Strain, R. M. and Yun, S.-B.: Propagation of uniform upper bounds for the spatially homogeneous relativistic Boltzmann equation, arXiv preprint arXiv:1907.05784. (2019). [46] Jang, J. W., Yun, S.-B.: Gain of regularity for the relativistic collision operator. Appl. Math. Lett. 90 (2019), 162-169. [47] Jang, J. W., Yun, S.-B.: Propagation of Lp estimates for the Spatially Homogeneous Relativistic Boltzmann Equation. arXiv preprint arXiv:2001.11672. (2020). [48] Kikuchi, S., Nordheim, L.: Uber¨ die kinetische Fundamentalgleichung in der Quantenstatistik. Zeitschrift für Physik A Hadrons and nuclei, 60 (1930), no. 9-10, 652-662. [49] Kim, M., Lee, C., Kim, Y. and Jeon, S.: Introduction to the DaeJeon Boltzmann-Uehling-Uhlenbeck (DJBUU) Project. New Physics, 66 (2016), no. 12, 1563-1570. [50] Lanford, O. E.: Time evolution of large classical systems. Dynamical systems, theory and applications (Rencontres, Battelle Res. Inst., Seattle, Wash., 1974), pp. 1-111. Lecture Notes in Phys., Vol. 38, Springer, Berlin, 1975. [51] Lapitski, D.: Development of the Quantum Lattice Boltzmann method for simulation of quantum elec- trodynamics with applications to graphene. Thesis (D.Phil.)–University of Oxford (United Kingdom). 2014. [52] Lemou, M.: Linearized quantum and relativistic Fokker-Planck-Landau equations. Math. Methods Appl. Sci. 23 (2000), no. 12, 1093-1119. [53] Li, B. A., Chen, L. W. and Ko, C. M.: Recent progress and new challenges in isospin physics with heavy-ion reactions. Physics Reports, 464 (2008), no. 4-6, 113-281. [54] Li, B. A., Ko, C. M., and Ren, Z.: Equation of state of asymmetric nuclear matter and collisions of neutron-rich nuclei. Physical Review Letters, 78 (1997), no. 9, 1644. [55] Li, W., Lu, X.: Global existence of solutions of the Boltzmann equation for Bose-Einstein particles with anisotropic initial data. J. Funct. Anal. 276 (2019), no. 1, 231-283. [56] Lichnerowicz, A., Marrot, R.: Propriétés statistiques des ensembles de particules en relativitérestreinte. C. R. Acad. Sci. Paris 210 (1940), 759-761. [57] Liu, S., Ma, X., Yu, H.: Optimal time decay of the quantum Landau equation in the whole space. J. Differential Equations. 252 (2012), no. 10, 5414-5452. [58] Lu, X.: A modified Boltzmann equation for Bose-Einstein particles: isotropic solutions and long-time behavior. J. Statist. Phys. 98 (2000), no. 5-6, 1335-1394. [59] Lu, X.: Long time convergence of the Bose-Einstein condensation. J. Statist. Phys. 162 (2016), no. 3, 652-670. [60] Lu, X.: On isotropic distributional solutions to the Boltzmann equation for Bose-Einstein particles. J. Statist. Phys. 116 (2004), no. 5-6, 1597-1649. [61] Lu, X.: On spatially homogeneous solutions of a modified Boltzmann equation for Fermi-Dirac particles. J. Statist. Phys. 105 (2001), no. 1-2, 353-388. [62] Lu, X.: On the Boltzmann equation for Fermi-Dirac particles with very soft potentials: averaging compactness of weak solutions. J. Statist. Phys. 124 (2006), no. 2-4, 517-547. [63] Lu, X.: On the Boltzmann equation for Fermi-Dirac particles with very soft potentials: global existence of weak solutions. J. Differential Equations. 245 (2008), no. 7, 1705-1761. [64] Lu, X.: The Boltzmann equation for Bose-Einstein particles: condensation in finite time. J. Statist. Phys. 150 (2013), no. 6, 1138-1176. [65] Lu, X.: The Boltzmann equation for Bose-Einstein particles: regularity and condensation. J. Statist. Phys. 156 (2014), no. 3, 493-545. [66] Lu, X.: The Boltzmann equation for Bose-Einstein particles: velocity concentration and convergence to equilibrium. J. Statist. Phys. 119 (2005), no. 5-6, 1027-1067. [67] Lu, X., Wennberg, B.: On stability and strong convergence for the spatially homogeneous Boltzmann equation for Fermi-Dirac particles. Arch. Ration. Mech. Anal. 168 (2003), no. 1, 1-34. [68] Lu, X., Zhang, X.: On the Boltzmann equation for 2D Bose-Einstein particles. J. Statist. Phys. 143 (2011), no. 5, 990-1019. [69] Saint-Raymond, L.: Hydrodynamic limits of the Boltzmann equation. Lecture Notes in Mathematics, 1971. Springer-Verlag, Berlin, 2009. [70] Spohn, H.: Large scale dynamics of interacting particles. Springer Science and Business Media. 2012. [71] Strain, R. M.: Asymptotic stability of the relativistic Boltzmann equation for the soft potentials. Comm. Math. Phys. 300 (2010), no. 2, 529-597. 42 GI-CHANBAE,JINWOOJANG,ANDSEOK-BAEYUN

[72] Strain, R. M.: Coordinates in the relativistic Boltzmann theory. Kinet. Relat. Models 4 (2011), no. 1, 345-359. [73] Strain, R. M.: Global Newtonian limit for the relativistic Boltzmann equation near vacuum. SIAM J. Math. Anal. 42 (2010), no. 4, 1568-1601. [74] Strain, R. M., Guo, Y.: Exponential decay for soft potentials near Maxwellian. Arch. Ration. Mech. Anal. 187 (2008), no. 2, 287-339. [75] Strain, R. M., Guo, Y.: Stability of the relativistic Maxwellian in a collisional plasma. Comm. Math. Phys. 251 (2004), no. 2, 263-320. [76] Strain, R. M., Yun, S.-B.: Spatially homogeneous Boltzmann equation for relativistic particles. SIAM J. Math. Anal. 46 (2014), no. 1, 917-938. R3 [77] Strain, R. M., Zhu, K.: Large-time decay of the soft potential relativistic Boltzmann equation in x. Kinet. Relat. Models. 5 (2012), no. 2, 383-415. [78] Succi, S.: Lattice Boltzmann equation for relativistic quantum mechanics. Philosophical Transactions of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences, 360 (2002), no. 1792, 429-436. [79] Uehling, E. A., Uhlenbeck, G. E.: Transport phenomena in einstein-bose and fermi-dirac gases. i. Physical Review, 43 (1933), no. 7, 552. [80] Van Weert, C. G.: Generalized hydrodynamics from relativistic kinetic theory. Physica A: Statistical Mechanics and its Applications, 111 (1982), no. 3, 537-552. [81] Villani, C.: A review of mathematical topics in collisional kinetic theory. Handbook of mathematical ﬂuid dynamics, Vol. I, 71-305, North-Holland, Amsterdam, 2002. [82] Wang, Y.: Global well-posedness of the relativistic Boltzmann equation. SIAM J. Math. Anal. 50 (2018), no. 5, 5637-5694.

Department of mathematics, Sungkyunkwan University, Suwon 16419, Republic of Korea Email address: [email protected]

University of Bonn, Institute for Applied Mathematics, Endenicher Allee 60, 53115 Bonn, Germany Email address: [email protected]

Department of mathematics, Sungkyunkwan University, Suwon 16419, Republic of Korea Email address: [email protected]