A review on BGK models for gas mixtures of mono and polyatomic molecules

We consider the socalled Bathnagar-Gross-Krook (BGK) model, an approximation of the Boltzmann equation, describing the time evolution of a single momoatomic rarefied gas and satisfying the same two main properties (conservation properties and entropy inequality). However, in practical applications one often has to deal with two additional physical issues. First, a gas often does not consist of only one species, but it consists of mixture a mixture of different species. Second, the particles can store energy not only in translational degrees of freedom but also in internal degrees of freedom like rotations or vibrations (polyatomic molecules). Therefore, here, we will present recent BGK models for gas mixtures for mono- and polyatomic particles and the existing mathematical theory for these models.


Introduction
In this paper we will concern ourselves with a kinetic description of gas mixtures. In the case of mono atomic molecules and two species this is usually done with the Boltzmann equation for the distribution functions f 1 = f 1 (x, v, t), f 2 = f 2 (x, v, t), see for example [18,19]. Here, x ∈ R d and v ∈ R d are the phase space variables, position and velocity of the particles, and t ≥ 0 denotes the time. The Boltzmann equation for gas mixtures is of the form where the collision operators Q 11 (f 1 , f 1 ) and Q 22 (f 2 , f 2 ) satisfy for k = 1, 2, and Q 12 (f 1 , f 2 ) and Q 21 (f 2 , f 1 ) satisfy Q 21 (f 2 , f 1 )dv = 0, (m 1 vQ 12 (f 1 , f 2 ) + m 2 vQ 21 (f 2 , f 1 ))dv = 0, (m 1 |v| 2 Q 12 (f 1 , f 2 ) + m 2 |v| 2 Q 21 (f 2 , f 1 ))dv = 0. (1. 2) The properties (1.1) ensure conservation of the number of particles, momentum and energy in interactions of one species with itself. The properties (1.2) ensure conservation of the number of particles, total momentum and total energy in interactions of one species with the other species. For the proof see for example [27]. In addition, the collision operator satisfy the inequalities [4] Q kk (f k , f k ) ln f k dv ≤ 0, k = 1, 2 (1.3) Q 12 (f 1 , f 2 ) ln f 1 dv + Q 21 (f 2 , f 1 ) ln f 2 dv ≤ 0 (1. 4) In the first inequality (1.3) we have equality if and only if f k is equal to a Maxwell distribution given by Note that in this paper we shall write T k instead of k B T k , where k B is Boltzmann's constant. Here, for any f 1 , f 2 , f 1 , f 2 ≥ 0, we relate the distribution functions to macroscopic quantities by mean-values of f k , where n k is the number density, u k the mean velocity and T k the mean temperature of species k (k = 1, 2). Additionally to a Maxwell distribution in the second inequality (1.4) we have equality if and only if u 1 = u 2 and T 1 = T 2 [4]. In the following we refer to these inequalities including the characterization of equality as entropy inequalities of H-Theorem.
If we are close to equilibrium [13,49], the complicated interaction terms of the Boltzmann equation can be simplified by a BGK approximation. This consists of a collision frequency ν ij n j multiplied by the deviation of the distribution functions from a local Maxwell distribution.
∂ t f 1 + v · ∇ x f 1 = ν 11 n 1 (M 1 − f 1 ) + ν 12 n 2 (M 12 − f 1 ), The collision frequencies per density ν kj are assumed to be dependent only on x and t and not on the microscopic velocity v. For references taking into account also a dependency on the microscopic velocity v see [54] for the one species case and [26] for the gas mixture case. Here, M k is given by (1.5) and M 12 , M 21 are given by for suitable u 12 , u 21 , T 12 , T 21 such that the conservation properties (1.2) are satisfied. In the literature, there is another type of approximation of the form with only one relaxation term taking into account both type of interactions, inter-and intra-species interactions, with one Maxwell distribution for suitable n (k) , u (k) , T (k) such that the conservation properties (1.2) are satisfied. BGK models give rise to efficient numerical computation, which are asymptotic preserving. This means that they remain efficient even approaching the hydrodynamic regime [7,8,20,22,24,48]. However, the BGK approximation has one drawback. This is its incapability of reproducing the correct Boltzmann hydrodynamic regime in the asymptotic continuum limit. Therefore, a modified version called the ES-BGK approximation was suggested by Holway for one species [34]. Then the H-Theorem of this model was shown in [47] and existence and uniqueness of mild solutions in [60]. Alternatively, the Shakov model [52] and a BGK model with velocity dependent collision frequency [54]was suggested to achieve the correct Prandtl number. For the BGK model with velocity dependent collision frequency, it is shown that a power law for the collision frequency can also lead to the proper Prandtl number. This BGK model with velocity dependent collision frequency should be also constructed in a way such that it satisfies the conservation properties. That this works one has to replace the Maxwell distribution by a different function, for details see [54]. For this model also an H-Theorem can be proven. The existence of these modified functions is proven in [26]. However, since BGK models can be the basis to build extended models as ES-BGK models, Shakov models and BGK models with velocity dependent collision frequency, we will mainly review BGK models for gas mixtures in this paper.
In this paper, we are interested in extensions of a BGK model to gas mixtures since in applications one often has to deal with mixtures instead of a single gas. Moreover, evolution of a polyatomic gas is very important in applications, for instance air consists of a gas mixture of polyatomic molecules. But most kinetic models modelling air deal with the case of a mono atomic gas consisting of only one species.
The outline of the paper is as follows: In section 2 we will present typical ansatzes of modelling gas mixtures with the BGK model and a review on recent results concerning existence of solutions and large-time behaviour. In section 3 we will give a summary over existing BGK models for gas mixtures of polyatomic molecules and a review on existing theoretical results.

BGK models for gas mixtures
In the following, we will present existing BGK models for gas mixtures. Then, we will give an overview over existing theoretical results (existence of solutions, large-time behaviour).

2.1.
Overview on existing BGK models for gas mixtures in the literature. Here we will focus on gas mixtures modelled via a BGK approach. In the literature one can find two types of BGK models for gas mixtures. Just like the Boltzmann equation for gas mixtures has a sum of collision terms on the right-hand side, one type of BKG models also contains a sum of BGK-type relaxation terms on the right-hand side (1.7). Examples are the models of Asinari [5], Cercignani [19], Garzo, Santos, Brey [25], Greene [28], Gross and Krook [31], Hamel [32], Sofena [53], and recent models by Bobylev, Bisi, Groppi, Spiga, Potapenko [12]; Haack, Hauck, Murillo [33] and by Klingenberg, Pirner, Puppo [38]. The other type of models contains only one collision term on the right-hand side (1.9). Examples for this are Andries, Aoki and Perthame [2], and the models in [15,29].A comparison of these models concerning their hydrodynamic limit can be found in [13] There are also many results concerning the hydrodynamic limit via the Chapman Enskog expansion, see for example [2,8,15,20,33] and extensions to ES-BGK models, Shakov models and BGK models with velocity dependent collision frequency [14,26,29,55].
2.1.1. BGK models for gas mixtures with one collision term. BGK models for gas mixtures with one interaction term [2,15,29] on the right-hand side have the form (1.9). Now, the interspecies velocities u (k) and temperatures T (k) in (1.10) have to be determined such that the conservation of total momentum and total energy as in (1.2) is satisfied. This gives d + 1 constraints for the 2(d + 1) quantities u (k) , T (k) , k = 1, 2. Therefore there is additional freedom to choose these quantities. Examples in the literature are given in [2,15,29]. In the first case [2] the quantities u (k) and T (k) are chosen such that the exchange terms of momentum and energy coincide with the exchange terms of momentum and energy of the Boltzmann equation for Maxwell molecules. For the details, see [2]. This leads to the choice where χ 12 , χ 21 are parameters which are related to the differential cross section. For the detailed expressions see [2]. The model also satisfies the conservation properties (1.2) and the H-theorem (1.4) . In the Htheorem one has equality if and only if the distribution functions are Maxwell distributions with the same mean velocity and temperature. The model of Andries, Aoki and Perthame has another property (see proposition 3.2 in [2]). It is called the indifferentiability principle. This means the following. When the masses m k , k = 1, 2 and the collision frequencies ν kj , k, j = 1, 2 are the same for each species, the total distribution function f = f 1 + f 2 satisfies a single species BGK equation.
A derivation of the Navier-Stokes system in the compressible regime and the corresponding transport coefficients can be found in section 4 of [2].
Another model in the literature with shape (1.9) is the model [29]. Here u (k) , T (k) are chosen such that all interspecies velocities u (k) and temperatures T (k) are equal ν s n s In [30], it is proved that the positivity of the temperatureT is guaranteed and the H-Theorem holds for the space-homogeneous case. The hydrodynamic limit and corresponding transport coefficients of this models can be found in section 5 in [13].
Another model with shape (1.9) is the model [15]. Here the aim was to derive the BGK model for gas mixtures from an entropy minimization principle ensuring that the model satisfies the exact Fick and Newtons laws in the hydrodynamical limit to the Navier-Stokes equations. This leads to a choice of different values for u (k) and equal values T (k) = T * for all k = 1, 2 for the temperatures. For the detailed expressions, see [15]. The transport equations of the hydrodynamic regime for this model can be found in section 5 of [15].
2.1.2. BGK models for gas mixtures with two collision terms. Now, we review BGK models for gas mixtures with two collision terms. In case of a gas mixture, if we assume that we only have binary interactions, there are two possibilities. The particles of one species can interact with themselves or with particles of the other species. One can take this into account by writing two interaction terms in the equations (1.7). This means that the right-hand side of the equations now consists of a sum of two relaxation operators. This structure is also described in [18,19]. This leads to two different types of equilibrium distributions. Due to an interaction of a species k with itself, we expect a relaxation to an equilibrium distribution M k . And due to the interaction of a species with the other species, we expect a relaxation towards a different mixture equilibrium distribution M kj .
The quantities ν kk n k are the one species collision frequencies, while the collision frequencies ν kj n j are related to interspecies interactions. To be flexible in choosing the relationship between the interspecies collision frequencies, we assume the following relationship The restriction on ε is without loss of generality. If ε > 1, we can exchange the notation 1 and 2 and choose 1 ε instead. Let us provide an example. We consider a plasma with electrons and ions. Let us first denote the electrons with the index e and ions with the index i. Then a common relationship found in the literature [6] is ν ie = me mi ν ei or equivalent ν ei = mi me ν ie . Now, if we want to use the notation 1 and 2, we have two possibilities. The first one is to choose the notation 1 for electrons and 2 for ions. In this case the mass ratio of the two particles is m2 m1 >> 1, and we get ν 12 = m2 m1 ν 21 . So we have ε = m2 m1 > 1. The other possibility is to choose the notation 1 for ions and 2 for electrons. In this case the mass ratio of the two kinds of particles is m2 m1 << 1, and we have ν 12 = m2 m1 ν 21 , therefore ε = m2 m1 < 1. So in this case, we would use the second choice for the notation. The condition (2.2) will enter in the proof of the H-Theorem.
Additionally, we assume that all collision frequencies are strictly positive. The Maxwell distribution M k in (1.8) has the same density, mean velocity and temperature as f k . With this choice, it can be guaranteed that we have conservation of the number of particles, momentum and energy in interactions of a species with itself (see section 2.2 in [38]). The remaining parameters u 12 , u 21 , T 12 , T 21 will now be determined using conservation of total momentum and total energy. Due to the choice of the densities, one can prove conservation of the number of particles, see Theorem 2.1 in [38]. We further assume that u 12 is a linear combination of u 1 and u 2 then we have conservation of total momentum provided that [38]. If we further assume that T 12 is of the following form then we have conservation of total energy provided that see Theorem 2.3 in [38]. In order to ensure the positivity of all temperatures, we need to restrict δ and γ to see Theorem 2.5 in [38]. For this model, one can prove an H-theorem as in (1.4) with equality if and only if f k , k = 1, 2 are Maxwell distributions with equal mean velocity and temperature, see [38].
The second last [12] presents an additional motivation how it can be derived formally from the Boltzmann equation. The last one [33] presents a Chapman-Enskog expansion with transport coefficients in section 5, a comparison with other BGK models for gas mixtures in section 6 and a numerical implementation in section 7.

2.2.
Theoretical results of BGK models for gas mixtures. In this section, we present theoretical results for the models presented in section 2.1. We start with reviewing some existing theoretical results for the one-species BGK model. Concerning the existence of solutions the first result was proven by Perthame in [46]. It is a result on global weak solutions for general initial data. This result was inspired by Diperna and Lion from a result on the Boltzmann equation [23]. In [47], the authors consider mild solutions and also obtain the uniqueness in the periodic bounded domain. There are also results of stationary solutions on a 1-dimensional finite interval with inflow boundary conditions in [58]. In a regime near a global Maxwell distribution, the global existence in the whole space R 3 was established in [59]. Concerning convergence to equilibrium Desvillettes proved strong convergence to equilibrium considering the thermalizing effect of the wall for reverse and specular reflection boundary condition in a periodic box [21]. In [51], the fluid limit of the BGK model is considered.
In the following we will present theoretical results for BGK models for gas mixtures.
2.2.1. Existence of solutions. First, we will present an existence result of mild solutions under the following assumptions for both type of models.
(1) We assume periodic boundary conditions in x. Equivalently, we can construct solutions satisfying for all i = 1, ..., d and a suitable {a i } ∈ R d with positive components, for k = 1, 2.
With these assumptions, we can show the following Theorem, existence of mild solutions in the following sense.
under the conditions of the collision frequencies (2.9) iff f 1 , f 2 satisfy
The proof can be found in [39].The main idea consists in proving Lipschitz continuity of the Maxwell distribution M kj and bounds on the macroscopic quantities needed for this.  .5) and (2.6), and to the initial value problem to (1.9) with (2.1) . Moreover, for all t > 0 the following bounds hold: for k, j = 1, 2, k = j and some constants A(t), B(t).

2.2.2.
Large-time behaviour. In this section, we will give an overview over existing results on the largetime behaviour for BGK models for gas mixtures. We denote by H(f ) = f ln f dv the entropy of a function f and by H(f |g) = f ln f g dv the relative entropy of f and g. Then, one can prove the following theorems. The proofs are given in [20].  Therefore, this theorem can also be proven in a similar way for the model (1.9) with (2.1), since a corresponding inequality for the model (1.9) with (2.1) of the form is proven in [2]. The next two theorems can also be easily extended to the model (1.9) with (2. where the constants are defined by C 1 = (1 − α)ν 12 (n 2 + n 1 ) , The proofs can be found in [20]. There are also results in the space-inhomogeneous case for linearized collision operator of the model (1.7) with (1.6), (2.3), (2.4), (2.5) and (2.6) for two species, see [40]. For this, we consider a solution (f 1 , f 2 ) to (1.7) which is close to the equilibrium ( (2.13) Now, we do Taylor expansion of the terms M 1 , M 2 , M 12 , M 21 with respect to σ 1 , σ 2 , µ 1 , µ 2 , τ 1 and τ 2 around 0 and only take first order terms. Moreover, one neglects quadratic terms of the form σ k σ l , σ k µ l and σ k τ l . Then one obtains the linearized system (2.14) and the following hypocoercivity result with some positive constants c d , C d that depend on L. Moreover, assume that ν 11 n ∞,1 + ν 12 n ∞,2 = 1 and ν 22 n ∞,2 + ν 21 n ∞,1 = 1, then any solution (h 1 (t), h 2 (t)) to (2.14) in dimension d = 1 with e(h 1 (0), h 2 (0)) < ∞, normalized according to then satisfies e(h1(t), h2(t)) ≤ e −Ct e(h1(0), h2(0)), whereC is given bỹ Here, µ d (L) is a one species decay rate developed in theorem 1.1 in [1].

BGK models for gas mixtures of polyatomic molecules
In this section, we review recent models for gas mixtures of polyatomic molecules. This means that we take into account that the particles can also store energy in degrees of freedom in internal energy like rotations and vibrations. We start with an overview on the one species models to introduce typical modeling aspects in this context.
First, we introduce a dependency on the degrees of freedom in internal energy. In the literature, this is done in different ways: • Discrete dependency on the degrees of freedom in internal energy Let us consider a system with l > 0 internal energy states E l . Then, we introduce l distribution functions f l (x, v, t, E l ), one distribution function for each internal energy state E l . This is for example considered in [42].

• A continuous scalar dependency on the degrees of freedom in internal energy
In this description we take into account degrees of freedom in internal energy by introducing an internal energy parameter I which takes into account all degrees of freedom with a continuous scalar variable. Then, we introduce a distribution function f (2) = f (2) (x, v, t, I). This is done for example in [3].
• A discrete and continuous dependency of the degrees of freedom in internal energy In this description, the rotations and vibrations of a diatomic gas are described in different ways. Here we define the distribution function f (3) = f (3) (x, v, t, I, i) where I ∈ R + 0 describes the internal energy in a continuous way, whereas i denotes the ith vibrational energy level of the corresponding vibrational energy iRhν kB (h is the Planck constant, R the fundamental gas constant, while ν is the fundamental vibrational frequency of the molecule). This is done for example in [41].
• A vector-valued continuous dependency of the degrees of freedom in internal energy Last, we consider a distribution function f (4) = f (4) (x, v, t, η) where η ∈ R l is a vector-valued continuous variable for the internal degrees of freedom, one component for each degree of freedom. Then |η| 2 has the meaning of microscopic energy in the internal degrees of freedom. Especially, one can also describe rotations and vibrations in a separate way. This is introduced in [9].
In this article, we focus on the treatment of the additional continuous variable. For this we consider a distribution function f (x, v, t, E) where E represents the dependency on the internal degrees of freedom and can be either a scalar I or a vector η. So f can be either f (2) , f (3) for a fixed i or f (4) . Then we define the macroscopic quantities as where we have e(I) = I 2/l in the scalar case and e(η) = |η −η| 2 in the vector-valued case. In many casesĒ is assumed to be zero, but we keep it here to be most general. Here, we note that in contrast to the monoatomic case (1.6), we have an additional temperature related to the degrees of freedom in internal energy T int . From the physical principle equipartition of internal energy, one expects that in equilibrium these two temperatures are the same. To achieve this, there are different strategies in modelling • Relaxation to an equilibrium distribution with equal temperature We define the equilibrium temperature as Then, we consider the equilibrium distribution with Λ l being a constant ensuring that the integral of M equ with respect to v and E is equal to n. and the BGK model If we multiply this equation with respect to e(E), we obtain the macroscopic equation in the spacehomogeneous case.
so T int relaxes towards T tr with a rate dependent on the collision frequency ν.
• Relaxation of the temperature due to a convex combination of temperatures in the Maxwell distribution We define T rel = θT equ + (1 − θ)T int , with 0 < θ ≤ 1. Then, we consider the following Maxwell distribution with the temperature T = (1 − θ)T tr + θT equ . Λ l is a constant ensuring that the integral of G[f ] with respect to v and I is equal to the density n. Then the model is given by with the collision frequency A ν . If we choose f (2) as distribution function, this corresponds to the model in [3]. For this model one can show conservation of the number of particles, momentum and total energy. Moreover, one can prove an entropy inequality. Here, the equilibrium is characterized by a Maxwell distribution with equal temperatures T equ = T tr = T int , for details see section 3 in [3]. One can also show that there exists a unique mild solution to this model. This is proven in [45].
With the convex combination in θ one takes into account that T tr and T int relax to the common value T equ . In the space-homogeneous case one can compute the following macroscopic equations These macroscopic equations describe a relaxation of T tr and T int towards T equ with a speed depending on the additional parameter θ. We see that the model captures the regime where this relaxation of the temperatures is slower than the relaxation of the distribution function to a Maxwell distribution since θ satisfies θ ≤ 1, so it reduces the speed of relaxation from A ν to A ν θ. This model satisfies the following assymptotic behaviour proven in [44] in the space-homogeneous case.
Theorem 3.0.1. Let 0 < θ ≤ 1. The distribution function for the spatially homogeneous case converges to equilibrium with the following rate: with the relative entropy H(f |g) = f ln f g dvdI for two functions f and g, and the Maxwell distribution M equ given by (3.2).
• Relaxation of the temperatures with an additional kinetic equation This concept was introduced in [9] for the distribution function f (4) . Here, we discribe the time evolution in the following way where νn is the collision frequency. Here, there appear to additional artificial temperatures Λ and Θ. In order to describe the time evolution of these two temperatures, we couple this kinetic equation with an algebraic equation for conservation of internal energy and a relaxation equation ensuring that the two temperatures Λ and Θ relax to the same value in equilibrium where Z r is a given parameter corresponding to the different rates of decay of translational and rotational/vibrational degrees of freedom. Here M is given by Note that we have The second equality follows from (3.6). The equation (3.7) is used to involve the temperature Θ. If we multiply (3.7) by e(E), integrate with respect to v and E and use (3.9), we obtain ∂ t (nΘ) + ∇ x · (nΘu) = νn Z r n(Λ − Θ) (3.10) a relaxation to a common value with a speed Z r not restricted to a slower speed as in (3.3). Therefore, the relaxation of the two temperatures to a common value can be slower or faster than the relaxation of the distribution function to the Maxwell distribution. This depends on the choice of Z r . The initial data of Λ is determined using (3.6).We see that in this model the term νn Zr d+l d (M equ − M ) plays the role to produce the relaxation of the two temperatures Λ and Θ to the same value. So in this model the effect of the relaxation to equal temperatures is done by coupling the BGK equation with an additional kinetic equation. If we choose f (4) as distribution function, this corresponds to the model in [9]. For this model one can prove conservation of the number of particles, momentum and total energy, and also an entropy inequality. Additionally, the equilibrium can be characterized by a Maxwell distribution with equal temperatures T equ = T tr = T int , for details see [9,50]. The existence of a unique mild solution can be proven similar to the existence in the momatomic case [47].
This model satisfies the following asymptotic behaviour in the space-homogeneous case proven in [50] for f = f (4) .
Assume that (f (4) , M ) is a solution of (3.4) coupled with (3.7) and (3.6). Then, in the space homogeneous case, we have the following convergence rate of the distribution functions f : where C is given by and the index 0 denotes the value at t = 0.
• Relaxation of the temperatures with an additional relaxation term This concept was introduced in [56]. Here, we add an additional relaxation term into the right-hand side where τ is the relaxation time of f towards a Maxwell distribution with the temperatures T tr and T int given by and τ Z ν with Z ν > 1 the relaxation time of this Maxwell distribution to the equilibrium distribution M equ with equal temperatures given by (3.8). So the relaxation to equilibrium is divided into two parts, first a relaxation towards an intermediate equilibrium distribution where the temperatures T tr and T int are separate, then to the final equilibrium with common temperatures.
3.1. Summary of existing BGK models for gas mixtures of polyatomic molecules in the literature. Now, for the gas mixture case, we will present different models [11,50,56] combining different ansatzes from the one species polyatomic case and the mixture modelling.
3.1.1. A BGK model for mixtures of polyatomic gases with one relaxation term. The BGK model we consider in this section was introduced by Bisi, Monaco and Soares in [11]. We introduce here two distribution functions with scalar continuous dependency on the degrees of freedom of internal energy f 1 (x, v, t, I) and f 2 (x, v, t, I). Then the time evolution of these distribution functions is described by two kinetic equations with one relaxation term on the right-hand side to the equilibrium distribution with common temperatures with the partition function q k (T ) = ∞ 0 φ k (I) exp(− I T )dI. Then, the parametersñ k ,ũ andT will be determined to have conservation of mass, total momentum and total energy. For the computation and the detailed expression see [11]. For this model also an entropy inequality can be proven in the space-homogeneous case, see [11]. Transport coefficients in the hydrodynamical limit of this model can be find in section 5 of [11].
3.1.2. A BGK model for mixtures of polyatomic gases with two relaxation terms. In this section, we present the model developed in [50]. This model has a vector-valued dependency on the internal energy. For this we introduce two numbers related to the degrees of freedom in internal energy. One is the total number of different rotational and vibrational degrees of freedom M ; and the other is l k , the number of internal degrees of freedom of species k, k = 1, 2. Moreover, η ∈ R M is the variable for the internal energy degrees of freedom, whereas η l k ∈ R M coincides with η in the components corresponding to the internal degrees of freedom of species k and is zero in all the other components. In this way, it is possible that the two species can have a different number of degrees of freedom in internal energy. Then, we have distribution functions f 1 (x, v, t, η l1 ) and f 2 (x, v, t, η l2 ). Their time evolution is described by , for j, k = 1, 2, j = k with the conditions (3.13) ν 12 = εν 21 , 0 < l 1 l 1 + l 2 ε ≤ 1.
The equation is coupled with conservation of internal energy (3.6) for each species, and an additional relaxation equation for j, k = 1, 2, j = k. M equ,k is given by (3.8) for each species. The additional M kj is defined by where T kj is given by For a certain choice of Λ kj and Θ kj one can prove conservation of mass, total momentum and total energy. For details see [50]. The existence of solutions for this model can be proven in the same way as it is done in [38] for the monoatomic case. In [50] they also prove an entropy inequality and the following decay to equilibrium Theorem 3.1.1. Assume that (f 1 , f 2 , M 1 , M 2 ) is a solution of (3.11) coupled with (3.14) and (3.6). Then, in the space homogeneous case, we have the following convergence rate of the distribution functions f 1 and f 2 : where C is given by C = min ν 11 n 1 + ν 12 n 2 , ν 22 n 2 + ν 21 n 1 , ν 11 n 1 z 1 + ν 12 n 2 , ν 22 n 2 z 2 + ν 21 n 1 , and the index 0 denotes the value at time t = 0.
There are also numerical results for this model in [57].

BGK model for mixtures of polyatomic gases with intermediate relaxation terms.
The model in [56] extends the idea of additional relaxation terms with intermediate equilibrium distributions from the one species case to gas mixtures. The model is of the form with Z r , Z ν > 1 and intermediate equilibrium distributions m s1 and m s2 . The detailed expressions of the intermediate equilibrium distributions can be found in [56] with a proof of the conservation properties. With standard methods one can also prove an entropy inequality. Transport coefficients of the hydrodynamic regime of this model can be found in section 4 of [56].

Conclusions
This paper reviews various existing BGK models for gas mixtures of mono and polyatomic molecules from the literature. In the case of monoatomic particles two types of models are presented.
One contains only one relaxation term on the right-hand side taking into account all types of interactions in one relaxation term. The other type of model separates the inter-and intra-species interactions by writing a sum of relaxation terms on the right-hand side. For both types of models a review on theoretical results concerning existence of solutions and convergence to equilibrium are given. The results on convergence to equilibrium consider both the space-homogeneous case for the full non-linear model and the space-inhomogeneous case for a linearized model.
In the polyatomic case first different ansatzes of modelling the degrees of freedom in internal energy for one species are considered: discrete or continuous; and scalar-or vector-valued.
Next, different ansatzes for modelling the relaxation of the temperature related to translational degrees of freedom and the temperature related to internal degrees of freedom to the same value for one species are presented. Here, also theoretical results on the convergence to equilibrium are presented.
Finally, three existing BGK models in the literature concerning gas mixtures of polyatomic molecules using the different ansatzes of polyatomic and gas mixture modelling are presented with existing theoretical results on the convergence to equilibrium.
However, BGK-type models often lack on correct parameters in the continuum limit like the Prandtl number. Therefore these models can be used as a basis for more extended models like ES-BGK models or velocity dependent collision frequency. As a future work the Chapman-Enskog expansion of the missing models can be computed and then the transport coefficients of all these models can be compared and eventually extended to match all parameters in the macroscopic equations. Here, the free parameters in the BGK model for monoatomic molecules with a sum of interaction terms might be useful.
funding: This research was funded by the Alexander von Humboldt foundation.