Bertrand Lods

Mustapha Mokhtar-Kharroubi
email: mustapha.mokhtar-kharroubi@univ-fcomte.fr

Ryszard Rudnicki
email: rudnicki@us.edu.pl

Invariant density and time asymptotics for collisionless kinetic equations with partly diffuse boundary operators

Keywords: 2010 Mathematics Subject Classification. Primary: 82C40; Secondary: 35F15, 47D06 Kinetic equation, Stochastic semigroup, convergence to equilibrium

INTRODUCTION

Kinetic transport equations in bounded geometry is an important field of investigation which can be traced back to the seminal work [START_REF] Bardos | Problèmes aux limites pour les équations aux dérivées partielles du premier ordre à coefficients réels; théorèmes d'approximation; application à l'équation de transport[END_REF] where absorbing boundary conditions have been considered. For more general boundary conditions, relating the incoming and outgoing fluxes at the boundary of the physical domain, the well-posedness of associated transport equations with general force terms -including Vlasov-like equations -have been considered in [START_REF] Beals | Abstract time-dependent transport equations[END_REF][START_REF] Arlotti | A new approach to transport equations associated to a regular field: trace results and well-posedness[END_REF][START_REF] Arlotti | On general transport equations with abstract boundary conditions. The case of divergence free force field[END_REF] while a thorough analysis of the free transport equation with abstract boundary conditions on general domains have been performed in [START_REF] Voigt | Functional analytic treatment of the initial boundary value problem for collisionless gases[END_REF] (see also [16, Appendix of § 2, p. 249]). Notice that, for a nonlinear and collisional kinetic equation such as Boltzmann equation, taking into account general boundary conditions induces notoriously additional difficulties; we just mention here the works [START_REF] Guo | Decay and continuity of the Boltzmann equation in bounded domains[END_REF] (dealing with close-to-equilibrium solutions) and [START_REF] Mischler | Kinetic equations with Maxwell boundary conditions[END_REF] (for renormalized solutions) and the references therein.

The object of this paper is to build a general theory of time asymptotics (t → ∞) for multidimensional collisionless kinetic semigroups with partly diffuse boundary operators. Our construction is twofold:

(1) On the one hand, we continue previous functional analytic works [START_REF] Arlotti | Substochastic semigroups for transport equations with conservative boundary conditions[END_REF][START_REF] Arlotti | Transport semigroup associated to positive boundary conditions of unit norm: a Dyson-Phillips approach[END_REF][START_REF] Mokhtar-Kharroubi | On collisionless transport semigroups with boundary operators of norm one[END_REF][START_REF] Voigt | Functional analytic treatment of the initial boundary value problem for collisionless gases[END_REF] on substochastic semigroups governing collisionless transport equations with conservative boundary operators in L 1 -spaces and combine them to recent developments on the asymptotics of stochastic partially integral semigroups in L 1 -spaces motivated by piecewise deterministic processes [START_REF] Rudnicki | Piecewise Deterministic Processes in Biological Models[END_REF]. [START_REF] Arkeryd | Boltzmann asymptotics with diffuse reflection boundary conditions[END_REF] On the other hand, we investigate the problem of the existence of invariant densities for collisionless transport equations. Such existence theory depends heavily on our understanding of compactness properties induced by the diffuse parts of the boundary operators. These compactness properties rely on the fine knowledge of smoothness and transversality properties of the dynamical flow induced by the semigroup. More precisely, we consider transport equations of the form

∂ t ψ(x, v, t) + v • ∇ x ψ(x, v, t) = 0, (x, v) ∈ Ω × V, t 0 (1.1a)
with initial data ψ(x, v, 0) = ψ 0 (x, v), (x, v) ∈ Ω × V, (1.1b) under abstract (conservative) boundary conditions

ψ |Γ -= H(ψ |Γ + ), (1.1c)
where

Γ ± = {(x, v) ∈ ∂Ω × V ; ±v • n(x) > 0}
(n(x) being the outward unit normal at v ∈ ∂Ω, see Figure 1) and H is a linear boundary operator relating the outgoing and incoming fluxes ψ |Γ + and ψ |Γ -and is bounded on the trace spaces

L 1 ± = L 1 (Γ ± ; |v • n(x)|π(dx) ⊗ m(dv)) = L 1 (Γ ± , dµ ± (x, v
)) where π denotes the Lebesgue surface measure on ∂Ω. We will focus our attention to the case of nonnegative and conservative boundary conditions, i.e.

Hψ 0 and

Hψ L 1 -= ψ L 1 + , for any nonnegative ψ ∈ L 1 + .

(1.2)

Here

Ω ⊂ R d (d 2
) is an open subset with C 1 boundary ∂Ω and our analysis takes place in the functional space

X = L 1 (Ω × V ; dx ⊗ m(dv))
where V ⊂ R d is the support of a nonnegative Borel measure m which is orthogonally invariant (i.e. invariant under the action of the orthogonal group of matrices in R d ). Such a measure covers the Lebesgue measure on R d , the surface Lebesgue measure on spheres (one speed or multi-group models) or even combinations of them. Very precise one-dimensional results corresponding to slab geometry have been obtained in [START_REF] Mokhtar-Kharroubi | On asymptotic stability and sweeping of collisionless kinetic equations[END_REF]. Their extension to multi-dimensional geometries (d 2) is far from being elementary and is completely open. It is the main concern of the present work to provide such a generalization.

Let

W = ϕ ∈ X; v • ∇ x ϕ ∈ X, ϕ |Γ ± ∈ L 1 (Γ ± ) FIGURE 1.
x ∈ ∂Ω; Γ 0 (x) -the tangent space to ∂Ω at x; Γ + (x) -outward velocities; Γ -(x) -inward velocities.

where v • ∇ x ϕ is meant in a distributional sense, (see Section 2 below for a reminder of the trace theory) and let T H : D(T H ) ⊂ X → X be defined by

T H ϕ = -v • ∇ x ϕ, D(T H )= {ϕ ∈ W ; ϕ |Γ -= H(ϕ |Γ + )}.
In contrast to the one-dimensional case [START_REF] Mokhtar-Kharroubi | On asymptotic stability and sweeping of collisionless kinetic equations[END_REF], in general, T H needs not be a generator. However, there exists a unique extension A ⊃ T H which generates a positive contraction C 0 -semigroup (U H (t)) t 0 , see [START_REF] Arlotti | Substochastic semigroups for transport equations with conservative boundary conditions[END_REF][START_REF] Mokhtar-Kharroubi | On collisionless transport semigroups with boundary operators of norm one[END_REF][START_REF] Voigt | Functional analytic treatment of the initial boundary value problem for collisionless gases[END_REF]. Notice that (U H (t)) t 0 needs not be stochastic, i.e. mass-preserving on the positive cone X + of X. Actually (U H (t)) t 0 is stochastic if and only if

A = T H (1.3)
and different characterizations of this property are also available [START_REF] Arlotti | Substochastic semigroups for transport equations with conservative boundary conditions[END_REF][START_REF] Mokhtar-Kharroubi | On collisionless transport semigroups with boundary operators of norm one[END_REF]. A general sufficient condition for (U H (t)) t 0 to be stochastic is given in Proposition 4.1 below. Let us briefly describe the main contributions of this paper. We restrict ourselves to the stochastic case (1.3). A very important role is played here by the irreducibility of (U H (t)) t 0 (see Definition 4.3 below). When T H is not a generator, it is not possible to handle easily its closure A = T H . Despite this fact, the resolvent of A is given by an "explicit" series converging strongly, see (2.1) below. By exploiting this series one can derive a very general sufficient criterion of irreduciblity of (U H (t)) t 0 in terms of properties of the stochastic boundary operator H, see Proposition 4.6 below. It is well known (see [START_REF] Davies | One-parameter Semigroups[END_REF]) that if the kernel of the generator of an irreducible stochastic C 0 -semigroup is not trivial (and consequently one-dimensional) then the semigroup is ergodic and converges strongly in Cesarò means to its one-dimensional (positive) ergodic projection (as t → +∞). Thus the existence of an invariant density of (U H (t)) t 0 is a cornerstone of this construction and is a fundamental problem for the understanding of the long-time behaviour of (1.1).

We mainly consider (local in space) stochastic boundary operators H : L 1 + → L 1 -which are (locally in space) convex combinations of reflection and diffuse operators of the form Hϕ(x, v) = α(x)Rϕ(x, v) + (1 -α(x))Kϕ(x, v) = α(x)ϕ(x, V(x, v)) + (1 -α(x))

Γ + (x) k(x, v, v )ϕ(x, v )µ x (dv )
where

Γ -(x, v) -→ (x, V(x, v)) ∈ Γ +
is a general µ-preserving reflection law, µ x (dv) = |v • n(x)| m(dv) and

Γ -(x) k(x, v, v )µ x (dv) = 1, (x, v ) ∈ Γ +
where α : x ∈ ∂Ω -→ α(x) ∈ [0, 1] is a measurable function.

Regarding the long-time behaviour of the solution to (1.1), when

sup x∈∂Ω α(x) < 1, (1.4)
we show under quite general assumptions on the kernel k(x, v, v ) that (U H (t)) t 0 is partially integral (i.e. for each t > 0, U H (t) dominates a non trivial integral operator). It follows that if (U H (t)) t 0 has an invariant density Ψ H then (U H (t)) t 0 is asymptotically stable, i.e.

lim t→+∞ U H (t)f -Ψ H = 0
for any density f ; see Theorem 7.5 for a precise statement. This result provides us with a much more precise result than the mere Cesarò convergence given by the general theory. Converse results are also given; indeed we show that if (U H (t)) t 0 has no invariant density then (U H (t)) t 0 is sweeping with respect to suitable sets. In a more precise way, the total mass of any trajectory of (1.1) t 0 -→ U H (t)ψ 0 concentrates for large time t → ∞ near small (or large) velocities or near the boundary ∂Ω × V , see Theorem 8.3 for a precise statement. Such asymptotics follow from general results on partially integral stochastic semigroups [START_REF] Ór | Continuous Markov semigroups and stability of transport equations[END_REF][START_REF] Ór | Asymptotic decomposition of substochastic operators and semigroups[END_REF][START_REF] Rudnicki | On asymptotic stability and sweeping for Markov operators[END_REF] which we recall in Appendix B of the paper. These general theorems on asymptotic stability or sweeping of stochastic collisionless kinetic semigroups (U H (t)) t 0 (and also some related results) are the first object of this paper. Our second object is to deal with the existence of an invariant density for stochastic collisionless kinetic semigroups (U H (t)) t 0 . As far as we know, the existence of an invariant density is known only for the classical Maxwell diffuse model (see Example 6.3 below) for which it is known that (U H (t)) t 0 is asymptotically stable [START_REF] Arkeryd | Boltzmann asymptotics with diffuse reflection boundary conditions[END_REF].

Thus our second object is to provide a general existence theory of invariant density for such kinetic models. We show first, for general stochastic boundary operators H, that 0 is an eigenvalue of T H associated to a nonnegative eigenfunction if and only if there exists a nonnegative solution ϕ ∈ L 1 + to the eigenvalue problem

M 0 Hϕ = ϕ, (1.5)
which satisfies the additional condition

Γ + ϕ(x, v)|v| -1 dµ + (x, v) < +∞ (1.6)
where

M 0 : L 1 -→ L 1 +
is the stochastic operator defined by

(M 0 ϕ) (x, v) = ϕ(x -τ -(x, v)v, v); (x, v) ∈ Γ + , ϕ ∈ L 1 - where τ -(x, v)
is the exit time function (see the definition in Section 2 below).

To study the existence of an invariant density, we introduce the sub-class of regular partly diffuse boundary operators such that the diffuse part is "weakly compact with respect to velocities" (see Definition 3.5 below) which enjoys nice approximation properties. The part of the paper concerned with the existence of an invariant density is very involved and is based on a series of highly technical results culminating in a key spectral stability result

r ess (M 0 H) = r ess (M 0 (αR)) (1.7)
(see Theorem 5.6) where r ess refers to the essential spectral radius. Since

r ess (M 0 (αR)) sup x∈∂Ω α(x)
then the spectral problem (1.5) has a solution under (1.4) (i.e. when the diffuse reflection is active everywhere on ∂Ω). If the corresponding eigenfunction satisfies the additional condition (1.6) then (U H (t)) t 0 is asymptotically stable. If not we show a more precise sweeping behaviour: the total mass of any trajectory t 0 -→ U H (t)ψ 0 of (1.1) concentrates near the zero velocity as t → +∞, see Theorem 8.5.

The above spectral stability result is a consequence of a key weak compactness theorem namely: for any integers k, 1

K(M 0 R) k M 0 K(M 0 R) M 0 K : L 1 + → L 1 - is weakly compact.
The proof of this important result (Theorem 5.1), using the Dunford-Pettis criterion, is highly technical and is given in numerous steps. Roughly speaking, the main difficulty lies in the fact that K induces compactness only in the velocity variables and several iterations and changes of variables are necessary to produce the missing compactness in the space variable x ∈ ∂Ω. Such changes of variables are non trivial and have to be carefully justified. To do this, we take advantage of the stochastic character of the various operators involved and we show (see Lemma A.11), up to µ-null sets, smoothness and transversality properties of the µ-preserving iterates

U -k • ξ -1 : Γ -→ Γ + (k ∈ N)
where ξ is the ballistic flow

ξ : (x, v) ∈ Γ + -→ ξ(x, v) = (x -τ -(x, v)v, v) and U : (x, v) ∈ Γ + -→ U(x, v) = (x -τ -(x, v)v, V (ξ(x, v))) ∈ Γ -.
A thorough analysis of the ballistic flow and of U is performed in Appendix A where several non trivial smoothness and non degeneracy results are given culminating in Lemma A.11 and involving intrinsic tools from differential geometry. These results are postponed in Appendix A for the simplicity of reading but we wish to point out that our analysis of the flow induced by (U H (t)) t 0 is new (even if results similar to some of ours appear e.g. in [START_REF] Guo | Decay and continuity of the Boltzmann equation in bounded domains[END_REF], see Remark A.6) and has its own interest independently of the main motivation of this paper.

As far as we know, most of our results are new and appear here for the first time. Finally, we note that the assumption that ∂Ω is of class C 1 plays a role only for the results on smoothness and transversality of the flow stated in Appendix A; it is likely that the results stated there remain valid for ∂Ω which is only piecewise of class C 1 .

The paper is organized as follows: in Section 2, we introduce the mathematical framework and notations used in the rest of the paper and establish several properties of the various operators involved in our subsequent analysis. In Section 3 we introduce and analyse the general class of boundary operators we investigate in the rest of the paper. Section 4 is devoted to general criteria for the ergodic convergence of the semigroup (U H (t)) t 0 (see Theorem 4.7) which is related to the study of the eigenvalue problem (1.5) as well as the irreducibility property of (U H (t)) t 0 . In Section 5 we establish the main technical result of the paper (Theorem 5.1) as well of its consequence on the stability of the essential radius (1.7), see Theorem 5.6. Section 6 is devoted to the main existence result for an invariant density, Theorem 6.6. The question of the asymptotic stability of (U H (t)) t 0 is then discussed in Section 7 while the sweeping properties of (U H (t)) t 0 , when no invariant density exists, are given in Section 8. As already mentioned, the paper ends with two Appendices. A first one, Appendix A contains all the technical results regarding the smoothness and transversality of the ballistic flow while Appendix B recall several important results about partially integral semigroup and sweeping properties used in Section 7 and 8.

We end this Introduction by mentioning that a related work dealing with rates of convergence to equilibrium is now in preparation [START_REF] Lods | Algebraic convergence to equilibrium for the transport equation with partly diffuse boundary conditions[END_REF] extending the results of [START_REF] Mokhtar-Kharroubi | Rates of convergence to equilibrium for collisionless kinetic equations in slab geometry[END_REF] devoted to slab geometry. Moreover, we hope also to take advantage of the tools developed here to revisit some important works (see e.g. [START_REF] Comets | Billiards in a General Domain with Random Reflections[END_REF][START_REF] Evans | Stochastic billiards on general tables[END_REF] and references therein) on stochastic billiards [START_REF] Lods | On stochastic billiards[END_REF].

MATHEMATICAL SETTING AND USEFUL FORMULAE

2.1. Functional setting. We introduce the partial Sobolev space

W 1 = {ψ ∈ X ; v • ∇ x ψ ∈ X}.
It is known [START_REF] Cessenat | Théorèmes de traces Lp pour les espaces de fonctions de la neutronique[END_REF][START_REF] Cessenat | Théorèmes de traces pour les espaces de fonctions de la neutronique[END_REF][START_REF] Dautray | Mathematical analysis and numerical methods for science and technology[END_REF]] that any ψ ∈ W 1 admits traces ψ |Γ ± on Γ ± such that

ψ |Γ ± ∈ L 1 loc (Γ ± ; dµ ± (x, v)) where dµ ± (x, v) = |v • n(x)|π(dx) ⊗ m(dv),
denotes the "natural" measure on Γ ± . Notice that, since dµ + and dµ -share the same expression, we will often simply denote them by

dµ(x, v) = |v • n(x)|π(dx) ⊗ m(dv),
the fact that it acts on Γ -or Γ + being clear from the context. Note that

∂Ω × V := Γ -∪ Γ + ∪ Γ 0 ,
where

Γ 0 := {(x, v) ∈ ∂Ω × V ; v • n(x) = 0}.
We introduce the space

W = ψ ∈ W 1 ; ψ |Γ ± ∈ L 1 ± .
One can show [START_REF] Cessenat | Théorèmes de traces Lp pour les espaces de fonctions de la neutronique[END_REF][START_REF] Cessenat | Théorèmes de traces pour les espaces de fonctions de la neutronique[END_REF] that

W = ψ ∈ W 1 ; ψ |Γ + ∈ L 1 + = ψ ∈ W 1 ; ψ |Γ -∈ L 1 -.
Then, the trace operators B ± :

B ± : W 1 ⊂ X → L 1 loc (Γ ± ; dµ ± ) ψ -→ B ± ψ = ψ |Γ ± , are such that B ± (W ) ⊆ L 1 ± .
Let us define the maximal transport operator T max as follows:

T max : D(T max ) ⊂ X → X ψ → T max ψ(x, v) = -v • ∇ x ψ(x, v), with domain D(T max ) = W 1 . Now, for any bounded boundary operator H ∈ B(L 1 + , L 1 -), define T H as T H ϕ = T max ϕ for any ϕ ∈ D(T H ),
where

D(T H ) = {ψ ∈ W ; ψ |Γ -= H(ψ |Γ + )}.
In particular, the transport operator with absorbing conditions (i.e. corresponding to H = 0) will be denoted by T 0 . We recall here that there exists a unique minimal extension (A, D(A)) of (T H , D(T H )) which generates a nonnegative C 0 -semigroup (U H (t)) t 0 in X. We note that D(A) ⊂ W 1 and Aϕ = -v • ∇ x ϕ = T max ϕ for any ϕ ∈ D(A) but the traces B ± ϕ need not to belong to L 1 (Γ ± , dµ ± ). The resolvent of A is given by

R(λ, A)f = R λ f + ∞ n=0 Ξ λ H (M λ H) n G λ f, ∀f ∈ X , λ > 0 (2.1)
where the series is strongly converging in X. See [4, Theorem 2.8] for details. Moreover, (U

H (t)) t 0 is a stochastic C 0 -semigroup, i.e. U H (t)f X = f X ∀f ∈ X + ; t 0 if and only if A = T H .
Actually, under suitable assumptions on H (see Prop. 4.1), A = T H so that (U H (t)) t 0 is stochastic.

2.2.

Exit time and integration formula. Let us now introduce the exit time of particles in Ω (with the notations of [START_REF] Arlotti | A new approach to transport equations associated to a regular field: trace results and well-posedness[END_REF]), defined as:

Definition 2.1. For any (x, v) ∈ Ω × V, define t ± (x, v) = inf{ s > 0 ; x ± sv / ∈ Ω}.
To avoid confusion, we will set τ ± (x, v)

:= t ± (x, v) if (x, v) ∈ ∂Ω × V.
With the notations of [START_REF] Guo | Decay and continuity of the Boltzmann equation in bounded domains[END_REF], t -is the backward exit time t b . From a heuristic viewpoint, t -(x, v) is the time needed by a particle having the position x ∈ Ω and the velocity -v ∈ V to reach the boundary ∂Ω. One can prove [START_REF] Voigt | Functional analytic treatment of the initial boundary value problem for collisionless gases[END_REF]Lemma 1.5]

that t ± (•, •) is measurable on Ω × V . Moreover τ ± (x, v) = 0 for any (x, v) ∈ Γ ± whereas τ ∓ (x, v) > 0 on Γ ± . It holds (x, v) ∈ Γ ± ⇐⇒ ∃y ∈ Ω with t ± (y, v) < ∞ and x = y ± t ± (y, v)v.
In that case, τ ∓ (x, v) = t ± (y, v). Notice also that,

t ± (x, v)|v| = t ± (x, ω) , ∀(x, v) ∈ Ω × V, v = 0, ω = |v| -1 v ∈ S d-1 . (2.2)
We have the following integration formulae from [START_REF] Arlotti | A new approach to transport equations associated to a regular field: trace results and well-posedness[END_REF].

Proposition 2.2. For any h ∈ X, it holds

Ω×V h(x, v)dx ⊗ m(dv) = Γ ± dµ ± (z, v) τ ∓ (z,v) 0 h (z ∓ sv, v) ds, (2.3)
and for any ψ ∈ L 1 (Γ -, dµ -),

Γ - ψ(z, v)dµ -(z, v) = Γ + ψ(x -τ -(x, v)v, v)dµ + (x, v).
(2.4)

Remark 2.3. Notice that with the notations introduced in [START_REF] Arlotti | A new approach to transport equations associated to a regular field: trace results and well-posedness[END_REF],

Γ ±∞ = {(x, v) ∈ Γ ± ; τ ∓ (x, v) = ∞} = {(x, v) ∈ Γ ± ; v = 0}
so that µ ± (Γ ±∞ ) = 0. This explains why the above integration formulae do not involve the sets Γ ±∞ . Moreover, because µ -(Γ 0 ) = µ + (Γ 0 ) = 0, we can extend the above identity (2.4) as follows:

for any ψ ∈ L 1 (Γ -∪ Γ 0 , dµ -) it holds Γ -∪Γ 0 ψ(z, v)dµ -(z, v) = Γ + ∪Γ 0 ψ(x -τ -(x, v)v, v)dµ + (x, v). (2.5)
2.3. About the resolvent of T H . For any λ ∈ C such that Reλ > 0, define

M λ : L 1 --→ L 1 + u -→ M λ u(x, v) = u(x -τ -(x, v)v, v)e -λτ -(x,v) , (x, v) ∈ Γ + ; Ξ λ : L 1 --→ X u -→ Ξ λ u(x, v) = u(x -t -(x, v)v, v)e -λt -(x,v) 1 {t -(x,v)<∞} , (x, v) ∈ Ω × V ;    G λ : X -→ L 1 + ϕ -→ G λ ϕ(x, v) = τ -(x,v) 0 ϕ(x -sv, v)e -λs ds, (x, v) ∈ Γ + ; and    R λ : X -→ X ϕ -→ R λ ϕ(x, v) = t -(x,v) 0 ϕ(x -tv, v)e -λt dt, (x, v) ∈ Ω × V ;
where 1 E denotes the charateristic function of the measurable set E. All these operators are bounded on their respective spaces. More precisely, for any Reλ > 0

M λ 1, Ξ λ (Reλ) -1 , G λ ( Reλ) -1 , R λ (Reλ) -1 .
The interest of these operators is related to the resolution of the boundary value problem:

(λ -T max )f = g, B -f = u, (2.6)
where λ > 0, g ∈ X and u is a given function over Γ -. Such a boundary value problem, with u ∈ L 1 -can be uniquely solved (see [START_REF] Arlotti | A new approach to transport equations associated to a regular field: trace results and well-posedness[END_REF])

Theorem 2.4. Given λ > 0, u ∈ L 1 -and g ∈ X, the function f = R λ g + Ξ λ u
is the unique solution f ∈ D(T max ) of the boundary value problem (2.6).

Remark 2.5. Notice that Ξ λ is a lifting operator which, to a given u ∈ L 1 -, associates a function f = Ξ λ u ∈ D(T max ) whose trace on Γ -is exactly u. More precisely,

T max Ξ λ u = λΞ λ u, B -Ξ λ u = u, B + Ξ λ u = M λ u, ∀u ∈ L 1 -. (2.7)
We can complement the above result with the following whose proof can be extracted from [7, Theorem 4.2]:

Proposition 2.6. If r σ (M λ H) < 1 (λ > 0), then A = T H and R(λ, T H ) = R λ + Ξ λ HR(1, M λ H)G λ
where the series converges in B(X).

2.4. Some auxiliary operators. For λ = 0, we can extend the definition of these operators in an obvious way but not all the resulting operators are bounded in their respective spaces. However, we see from the above integration formula (2.4), that

M 0 ∈ B(L 1 -, L 1 + ) with M 0 u L 1 + = u L 1 -, ∀u ∈ L 1 -.
In the same way, one deduces from (2.3) that for any nonnegative ϕ ∈ X:

Γ + G 0 ϕ(x, v)dµ + (x, v) = Γ + dµ + (x, v) τ -(x,v) 0 ϕ(x -sv, v)ds = Ω×V ϕ(x, v)dx ⊗ m(dv) (2.8) which proves that G 0 ∈ B(X, L 1 + ) with G 0 ϕ L 1 + = ϕ X , ∀ϕ ∈ X.
To be able to provide a rigorous definition of the operators Ξ 0 and R 0 we need the following Definition 2.7. Introduce the function spaces

Y ± 1 = L 1 (Γ ± , |v| -1 dµ ± ) with its associated L 1 -norm • Y ± 1 and X τ = L 1 (Ω × V, t + (x, v)dx ⊗ m(dv)) with the associated L 1 -norm • τ .
The interest of the above boundary spaces lies in the following:

Lemma 2.8. For any u ∈ Y - 1 one has Ξ 0 u ∈ X with Ξ 0 u X = Γ - u(x, v)τ + (x, v)dµ + (x, v) D u Y - 1 , ∀u ∈ Y - 1 (2.9)
where D is the diameter of Ω,

D = sup x,y∈∂Ω |x -y|. Moreover, if u ∈ Y - 1 then M 0 u ∈ Y + 1 with M 0 u Y + 1 = u Y - 1 .
(2.10)

If f ∈ X τ then G 0 f ∈ Y - 1 and R 0 f ∈ D(T 0 ) ⊂ X and T 0 R 0 f = -f . Proof. From (2.3), for nonnegative u ∈ L 1 -: Ω×V Ξ 0 u(x, v)dx ⊗ m(dv) = Ω×V Ξ 0 u(x, v)dx ⊗ m(dv) = Γ + dµ + (z, v) τ -(z,v) 0 u(z -sv -t -(z -sv, v)v, v)1 {t -(z-sv,v)<∞} ds = Γ + u(z -τ -(z, v)v, v)τ -(z, v)dµ + (z, v)
which, using now (2.4) yields (2.9). If now u ∈ Y - 1 , then

Γ + M 0 u(x, v) |v| -1 dµ + (x, v) = Γ + u(x -τ -(x, v)v, v)|v| -1 dµ + (x, v)
and we deduce from (2.4) that

Γ + M 0 u(x, v)|v| -1 dµ + (x, v) = Γ - u(z, v)|v| -1 dµ -(z, v)
which is (2.10). If now f ∈ X τ is nonnegative, one has directly from (2.3) that

Γ + |v| -1 G 0 f (x, v)dµ + (x, v) = |v| -1 f X D f τ < ∞ which proves that G 0 f ∈ X. Moreover, using (2.3), Ω×V R 0 f (x, v)dx ⊗ m(dv) = Γ + dµ + (z, v) τ -(z,v) 0 [R 0 f ](z -sv, v)ds = Γ + dµ + (z, v) τ -(z,v) 0 ds t -(z-sv,v) 0 f (z -sv -tv, v)dt and, since t -(z -sv, v) = τ -(z, v) -s for all (z, v) ∈ Γ + and all 0 < s < τ -(z, v) we get Ω×V R 0 f (x, v)dx ⊗ m(dv) = Γ + dµ + (z, v) τ -(z,v) 0 ds τ -(z,v) s f (z -tv, v)dt = Γ + dµ + (z, v) τ -(z,v) 0 t f (z -tv, v)dt Now, since t + (z -tv, v) = t for any (z, v) ∈ Γ + , the above reads Ω×V R 0 f (x, v)dx ⊗ m(dv) = Γ + dµ + (z, v) τ -(z,v) 0 t + (z -tv, v) f (z -tv, v)dt
and, using again (2.3), one gets

Ω×V R 0 f (x, v)dx ⊗ m(dv) = Ω×V t + (x, v)f (x, v)dx ⊗ m(dv).
This proves that R 0 f ∈ X. Now, it is easy to see that actually g = R 0 f satisfies T max g = -f and B -g = 0, i.e. g ∈ D(T 0 ) with T 0 g = -f. Remark 2.9. Notice that, for any nonnegative u ∈ L 1 + ,

Γ + M 0 u(x, v)τ -(x, v)dµ + (x, v) = Γ + u(x -τ -(x, v)v, v)τ -(x, v)dµ + (x, v) and, since τ + (x -τ -(x, v)v, v) = τ -(x, v) for any (x, v) ∈ Γ -, we deduce from (2.4) that Γ + M 0 u(x, v)τ -(x, v)dµ + (x, v) = Γ - u(z, v)τ + (z, v)dµ -(z, v)
This shows that, in (2.10), we can replace

Y ± 1 with L 1 (Γ ± , τ ∓ (x, v)dµ ± (x, v)).
In the same way, one see that, for g ∈ X τ it holds G 0 g L 1 (Γ -,τ + dµ -) = g τ .

One has the following result: [START_REF] Bogachev | Measure Theory[END_REF]. To prove the uniqueness, it suffices to assume that g = u = 0 but then (2.11) reads T 0 f = 0 which admits the unique solution f = 0.

Proposition 2.10. Let g ∈ X τ be given and u ∈ L 1 (Γ -, τ + (x, v)dµ -). The boundary value problem -T max f = g B -f = u (2.11) admits a unique solution f ∈ X given by f = R 0 g + Ξ 0 u. Proof. Let u ∈ L 1 (Γ -, τ + (x, v)dµ -) and g ∈ X τ , since Ξ 0 u ∈ D(T max ) with T max Ξ 0 u = 0, B + Ξ 0 u = M 0 u and B -Ξ 0 u = u one sees that f ∈ D(T max ) with T max f = T max R 0 g + T max Ξ 0 u = T 0 R 0 g = -g while B -f = B -R 0 g + B -Ξ 0 u = B -Ξ 0 u = u. This shows that f = R 0 g + Ξ 0 u is a solution to (2.

GENERAL STOCHASTIC PARTLY DIFFUSE BOUNDARY CONDITIONS

Let us explicit here the general class of boundary conditions we aim to deal with. Typical boundary operators arising in the kinetic theory of gases are local with respect to x ∈ ∂Ω. In order to exploit this local nature of the boundary conditions, we introduce the following notations. For any x ∈ ∂Ω, we define

Γ ± (x) = {v ∈ V ; ±v • n(x) > 0}, Γ 0 (x) = {v ∈ V ; v • n(x) = 0}
and we define the measure µ x (dv) on Γ ± (x) given by

µ x (dv) = |v • n(x)|m(dv).
This allows to define the L 1 -space L 1 (Γ ± (x), dµ x ) in an obvious way. We shall denote the

L 1 (Γ ± (x), µ x ) norm by • L 1 (Γ ± (x)) . Since, for any ϕ ∈ L 1 (Γ ± , µ ± ) one has ϕ L 1 ± = ∂Ω Γ ± (x) |ϕ(x, v)|µ x (dv) π(dx) = ∂Ω ϕ(x, •) L 1 (Γ ± (x)) π(dx)
we can identify isometrically any ϕ ∈ L 1 ± to the field

x ∈ ∂Ω -→ ϕ(x, •) ∈ L 1 (Γ ± (x)).
(3.1)

3.1. Reflection boundary operators. We begin with the following definition of pure reflection boundary conditions (see [34, Definition 6.1, p.104]):

Definition 3.1. One says that R ∈ B(L 1 + , L 1 -) is a pure reflection boundary operator if R(ϕ)(x, v) = ϕ(x, V(x, v)) ∀(x, v) ∈ Γ -, ϕ ∈ L 1 + where V : x ∈ ∂Ω → V(x, •
) is a field of bijective bi-measurable and µ x -preserving mappings

V(x, •) : Γ -(x) ∪ Γ 0 (x) → Γ + (x) ∪ Γ 0 (x) such that i) |V(x, v)| = |v| for any (x, v) ∈ Γ -. ii) If (x, v) ∈ Γ 0 then (x, V(x, v)) ∈ Γ 0 , i.e. V(x, •) maps Γ 0 (x) in Γ 0 (x).
iii) The mapping

(x, v) ∈ Γ -→ (x, V(x, v)) ∈ Γ + is a C 1 diffeomorphism.
Remark 3.2. This last regularity property on V may require additional regularity of ∂Ω as seen in Example 3.3.

Note that

[Rϕ](x, v) = [R(x)ϕ(x, •)](v), ∀ϕ ∈ L 1 + , (x, v) ∈ Γ - where we identify (isometrically) ϕ ∈ L 1 + to the integrable field (3.1) and x ∈ ∂Ω → R(x) ∈ B(L 1 (Γ + (x), Γ -(x))
is the field of operators defined by

R(x)ψ(v) = ψ(V(x, v)), ψ ∈ L 1 (Γ + (x)), v ∈ Γ -(x). It holds R(x)ψ L 1 (Γ -(x)) = ψ L 1 (Γ + (x)) , ∀x ∈ ∂Ω, ψ ∈ L 1 (Γ + (x)), therefore, R ∈ B(L 1 + , L 1 -) is stochastic since Rϕ L 1 -= ∂Ω R(x)ϕ(x, •) L 1 (Γ -(x)) π(dx) = ∂Ω ϕ(x, •) L 1 (Γ + (x)) π(dx) = ϕ L 1 + .
Notice that this last identity is equivalent to the property that the mapping

(x, v) ∈ Γ --→ (x, V(x, v)) ∈ Γ + is µ-preserving.
Example 3.3. In practical situations, the most frequently used pure reflection conditions are (a) the specular reflection boundary conditions which corresponds to the case in which V and m are invariant under the orthogonal group and

V(x, v) = v -2(v • n(x)) n(x) (x, v) ∈ Γ -.
Notice that, for V to be a C 1 diffeormorphism, we need ∂Ω to be of class C 2 . (b) The bounce-back reflection conditions for which V(x, v) = -v, (x, v) ∈ Γ -.

Diffuse boundary operators. We introduce the following definition

Definition 3.4. One says that K ∈ B(L 1 + , L 1 -) is a stochastic diffuse boundary operator if K ψ(x, v) = Γ + (x) k(x, v, v )ψ(x, v )µ x (dv ), (x, v) ∈ Γ -, ψ ∈ L 1 + (3.2)
where the kernel k(x, v, v ) induces a field of nonnegative measurable functions

x ∈ ∂Ω → k(x, •, •) where k(x, •, •) : Γ -(x) × Γ + (x) → R + is such that Γ -(x) k(x, v, v )µ x (dv) = 1, ∀(x, v ) ∈ Γ + . FIGURE 2.
Regular and diffuse reflection: v -an outward vector, V(x, v) -the specular reflection, thin vectors -diffuse reflection.

As we did for reflection operators, we identify

K ∈ B(L 1 + , L 1 -) to a field of integral operators x ∈ ∂Ω -→ K(x) ∈ B(L 1 (Γ + (x), Γ -(x))
by the formula

[Kψ](x, v) = [K(x)ψ(x, •)] (v)
where, for any x ∈ ∂Ω

K(x) : ψ ∈ L 1 (Γ + (x)) -→ [K(x)ψ] (v) = Γ + (x) k(x, v, v )ψ(v )µ x (dv ) ∈ L 1 (Γ -(x)). Note that K(x) : L 1 (Γ + (x)) → L 1 (Γ -(x)
) is stochastic for any x ∈ ∂Ω and therefore so is

K ∈ B(L 1 + , L 1 -), i.e. Kψ L 1 -= ψ L 1 + ∀ψ ∈ L 1 + .
We introduce now a useful class of diffuse boundary operators. Before giving the formal definition, let us recall that, if

K ∈ B(L 1 + , L 1 -) given by (3.2) is such that K(x) ∈ B(L 1 (Γ + (x)), L 1 (Γ -(x))) is weakly compact for any x ∈ ∂Ω (3.3)
then, according to the Dunford-Pettis criterion (see [START_REF] Ézis | Functional analysis, Sobolev spaces and partial differential equations[END_REF]Theorem 4.30,p. 115 & Exercise 4.36, p. 129]), for any x ∈ ∂Ω and any ε > 0, there is δ > 0 such that

sup v ∈Γ + (x) A k(x, v, v )µ x (dv) < ε ∀A ⊂ Γ -(x) such that µ x (A) < δ and lim m→∞ sup v ∈Γ + (x) Γ -(x)\Am k(x, v, v )µ x (dv) = 0 for any sequence (A m ) m ⊂ Γ -(x) with A m ⊂ A m+1 , µ x (A m ) < ∞ and ∪ m A m = Γ -(x).
In particular, for any x ∈ ∂Ω,

lim m→∞ sup v ∈Γ + (x) {v∈Γ -(x) ; |v| m} k(x, v, v ) µ x (dv) = 0.
Moreover, since

1 = Γ -(x) k(x, v, v )µ x (dv) {v∈Γ -(x) ; k(x,v,v ) m} k(x, v, v )µ x (dv) m µ x {v ∈ Γ -(x) ; k(x, v, v ) m} , ∀m ∈ N, (x, v ) ∈ Γ + ,
we have lim

m→∞ sup v ∈Γ + (x) µ x {v ∈ Γ -(x) ; k(x, v, v ) m} = 0.
In other words, for any x ∈ ∂Ω, the following holds

lim m→∞ sup v ∈Γ + (x) Sm(x,v ) k(x, v, v ) µ x (dv) = 0 (3.4)
where, for any m ∈ N and any

(x, v ) ∈ Γ + S m (x, v ) = {v ∈ Γ -(x) ; |v| m} ∪ {v ∈ Γ -(x) ; k(x, v, v ) m}.
We introduce then the following class of diffuse boundary operators: Definition 3.5. We say that a diffuse boundary operator

K ∈ B(L 1 + , L 1 -) is regular if the family of operators K(x) ∈ B(L 1 (Γ + (x)), L 1 (Γ -(x))),
x ∈ ∂Ω is collectively weakly compact in the sense that (3.3) holds true for any x ∈ ∂Ω and the convergence in (3.4) is uniform with respect to x ∈ ∂Ω. Remark 3.6. A diffuse boundary operator K is regular for instance whenever there exists h :

V → R + such that V h(v) |v| m(dv) < +∞ and k(x, v, v ) h(v) ∀x ∈ ∂Ω, v ∈ Γ + (x), v ∈ Γ -(x).
In particular, the classical Maxwell boundary operator (see Example 6.3 below) is a regular diffuse boundary operator.

We have then the following approximation result.

Lemma 3.7. Assume that K ∈ B(L 1 + , L 1
-) is a regular diffuse boundary operator in the sense of the above definition. Then, there exists a sequence

(K m ) m ⊂ B(L 1 + , L 1 -) such that (1) 0 K m K for any m ∈ N; (2) lim m→∞ K -K m B(L 1 + ,L 1 -) = 0; (3) For any m ∈ N and any nonnegative f ∈ L 1 + it holds K m f (x, v) ψ m (v) Γ + (x) f (x, v ) |v • n(x)| m(dv ), (x, v) ∈ Γ - (3.5) with ψ m = m1 Bm where B m = {v ∈ R d ; |v| m}. Proof. Let k(x, v, v ) be the kernel associated to K through (3.2). Introduce then k m (x, v, v ) = inf{k(x, v, v ) ; m1 Bm (v)} for any m ∈ N
, where B m is the ball of R d centered in 0 and with radius m, and set

K m ϕ(x, v) = Γ + (x) k m (x, v, v )ϕ(x, v ) |v • n(x)| m(dv ), ϕ ∈ L 1 + , (x, v) ∈ Γ -.
Clearly, K m ∈ B(L 1 + , L 1 -) is a diffuse boundary operator with 0 K m K and (3.5) holds. Moreover, for any x ∈ ∂Ω and any ϕ ∈ L 1 (Γ + (x)), it is easy to check that

K(x)ϕ -K m (x)ϕ L 1 (Γ -(x)) ϕ L 1 (Γ + (x)) × sup v ∈Γ + (x) {v∈Γ -(x) ; k(x,v,v ) m1 Bm (v)} k(x, v, v )µ x (dv), i.e. K(x) -K m (x) B(L 1 (Γ + (x)),L 1 (Γ -(x))) sup v ∈Γ + (x) {v∈Γ -(x) ; k(x,v,v ) m1 Bm (v)} k(x, v, v )µ x (dv) sup v ∈Γ + (x) Sm(x,v ) k(x, v, v )µ x (dv).
One sees then that

K -K m B(L 1 + ,L 1 -) = sup x∈∂Ω K(x) -K m (x) B(L 1 (Γ + (x)),L 1 (Γ -(x)))
goes to zero as m → ∞ since the convergence in (3.4) is uniform with respect to x ∈ ∂Ω.

We complement the above result with a different kind of approximation which will turn useful in Section 8: Lemma 3.8. Let K be a regular stochastic diffuse boundary operator with kernel k(x, v, v ). Let

β n (x, v ) = Γ -(x)∩{|v|> 1 n } k(x, v, v )µ x (dv), (x, v ) ∈ Γ + and k n (x, v, v ) = k(x, v, v ) β n (x, v ) 1 {|v|> 1 n } , x ∈ ∂Ω, v ∈ Γ + (x), v ∈ Γ -(x)
Then, denoting by K n the regular stochastic diffuse boundary operator with kernel k n , it holds (i)

lim n→+∞ β n (x, v ) = 1 uniformly in (x, v ) ∈ Γ + . (ii) lim n→∞ K n -K B(L 1 + ,L 1 -) = 0. Proof. (i) For any x ∈ ∂Ω, set A n (x) = v ∈ Γ -(x) ; |v| < 1 n . One has µ x (A n (x)) = Γ -(x)∩{|v|<n -1 } |v • n(x)|m(dv) m(B 1 ) n
where B 1 is the unit ball of R d . Thus,

lim n→∞ sup x∈∂Ω µ x (A n (x)) = 0. (3.6)
Now, since K is regular, (3.4) holds uniformly with respect to x ∈ ∂Ω, i.e. for any ε > 0, there is m ∈ N large enough so that sup

(x,v )∈Γ + Sm(x,v ) k(x, v, v )µ x (dv) < ε.
Then, for any n ∈ N and any (x, v ) ∈ Γ + ,

An(x) k(x, v, v )µ x (dv) = An(x)∩Sm(x,v ) k(x, v, v )µ x (dv) + An(x)\Sm(x,v ) k(x, v, v )µ x (dv) ε + m µ x (A n (x)) since k(x, v, v ) m on Γ -(x) \ S m (x, v ). Using then (3.6) we get lim n→∞ sup (x,v )∈Γ -Γ -(x)∩{|v|<n -1 } k(x, v, v )µ x (dv) = 0 (3.7) which shows (i) since Γ -(x) k(x, v, v )µ x (dv) = 1 for any x ∈ ∂Ω.
(ii) Set K n the boundary operator with kernel

1 {|v|>n -1 } k(x, v, v ). One checks easily that K -K n B(L 1 + ,L 1 -) sup (x,v )∈Γ -Γ -(x)∩{|v|<n -1 } k(x, v, v )µ x (dv) so that lim n K -K n B(L 1 + ,L 1 -)
= 0 from (3.7). Since moreover

K n -K n B(L 1 + ,L 1 -) sup (x,v )∈Γ + 1 -β n (x, v ) K B(L 1 + ,L 1 -)
which goes to zero from point (i), we get the desired result.

3.3. Stochastic partly diffuse boundary operators. We introduce now the general class of boundary operator we aim at investigating. Definition 3.9. We shall say that a boundary operator

H ∈ B(L 1 + , L 1 -) is stochastic partly diffuse if it writes Hψ(x, v) = α(x) Rψ(x, v) + (1 -α(x)) Kψ(x, v), (x, v) ∈ Γ -, ψ ∈ L 1 + (3.8)
where α(•)

: ∂Ω → [0, 1] is measurable, R is a reflection operator, and K ∈ B(L 1 + , L 1 -) is a stochastic diffuse boundary operator given by (3.2).
If the diffuse part K is regular we say that H is a regular stochastic partly diffuse boundary operator.

Remark 3.10. Notice that, being a convex combination of stochastic operators, a stochastic partly diffuse operator H is stochastic.

GENERAL RESULTS FOR ABSTRACT STOCHASTIC BOUNDARY OPERATORS

We begin with the following which is a direct consequence of [START_REF] Mokhtar-Kharroubi | On collisionless transport semigroups with boundary operators of norm one[END_REF]Theorem 21]:

Proposition 4.1. Let H ∈ B(L 1 + , L 1
-) be a stochastic boundary operator. Let there exist ϕ ∈ L 1 + such that ϕ > 0 µ-a.e. and ϕ M 0 Hϕ. Then, A = T H and (U H (t)) t 0 is a stochastic C 0semigroup, where we recall that (A, D(A)) is the generator of (U H (t)) t 0 .

We give a general result about the spectrum of T H :

Proposition 4.2. Let H ∈ B(L 1 + , L 1 -) be a stochastic boundary operator. Then, there is a nonnega- tive Ψ ∈ D(T H ) with T H ψ = 0 if and only if 1 is an eigenvalue of M 0 H associated to a nonnegative eigenfunction ϕ ∈ Y + 1 such that Hϕ ∈ Y - 1 . Proof. Assume first that there exists ϕ ∈ Y + 1 such that M 0 Hϕ = ϕ,
and Hϕ ∈ Y - 1 . Then, as already seen (see (2.9)), Ξ 0 Hϕ ∈ X. Let Ψ = Ξ 0 Hϕ. One has

T max Ψ = 0, B -Ψ = B -Ξ 0 Hϕ = Hϕ, i.e. B -Ψ = HB + Ψ. This means that Ψ ∈ D(T H ) with T H Ψ = 0. Assume now that 0 ∈ S p (T H ) is associated to a nonnegative eigenfunction Ψ ∈ D(T H ). Let ϕ = B + Ψ and u = B -Ψ. It holds u = Hϕ and, solving the boundary value problem (2.11) (see Proposition 2.10) yields Ψ = Ξ 0 u. It is easy to check then that ϕ = M 0 Hϕ. Since u = 0 =⇒ ϕ = 0, we get 1 ∈ S p (M 0 H).
We recall the definition of irreducible operators or semigroups in L 1 -spaces and refer to [START_REF] Arendt | One-parameter Semigroups of Positive Operators[END_REF] for more details. ii) B is positivity improving if for any h ∈ L 1 (E, Σ, m) non identically zero

h(s) 0 for m-a. e. s ∈ E =⇒ Bh(s) > 0 for m-a. e. s ∈ E.
iii) B is irreducible if for any non trivial and nonnegative h ∈ L 1 (E, Σ, m) and any non trivial nonnegative g ∈ L ∞ (E, Σ, m), there exists n ∈ N such that

B n h , g 1,∞ > 0 where •, • 1,∞ is the duality bracket between L 1 (E, Σ, m) and L ∞ (E, Σ, m). iv) A positive C 0 -semigroup (S(t)) t 0 on L 1 (E, Σ, m) with generator (G, D(G)) is irreducible
if, for any non trivial nonnegative h ∈ L 1 (E, Σ, m) and any non trivial nonnegative g ∈ L ∞ (E, Σ, m), there exists t 0 such that S(t)h, g 1,∞ > 0. This property is equivalent to the fact that R(λ, G) is positivity improving for λ > 0 large enough.

We introduce the following:

Assumption 4.4. H ∈ B(L 1 + , L 1 -) is a stochastic operator such that M 0 H ∈ B(L 1 +
) is irreducible and there exists ψ 0 ∈ L 1 + such that Hψ 0 > 0 µ-a.e. on Γ + . Remark 4.5. If H is stochastic partly diffuse operator of the form (3.8)

and if M 0 (1 -α)K is irreducible then so is M 0 H. This is the case for instance if α ∞ < 1 and if, for π-a.e. x ∈ ∂Ω k(x, v, v ) > 0 for µ x -a.e. v ∈ Γ -(x), v ∈ Γ + (x). (4.1)
In addition, if α(x) > 0 π-a.e. x ∈ ∂Ω, the second condition Hψ 0 > 0 µ-a.e. is satisfied by any ψ 0 > 0 µ-a.e. Otherwise, we assume that, for any x ∈ ∂Ω such that α(x) = 0 there is a subset

γ + (x) ⊂ Γ + (x) of positive µ x -measure and such that k(x, v, v ) > 0 for µ x -a.e. v ∈ Γ + (x), v ∈ Γ + (x).
Then, again the second condition is satisfied by any ψ 0 > 0 µ-a.e. In particular, the second assumption in Assumptions 4.4 is always satisfied under (4.1).

One has the following (see also [START_REF] Mokhtar-Kharroubi | On collisionless transport semigroups with boundary operators of norm one[END_REF]Remark 20]):

Proposition 4.6. Let H ∈ B(L 1 + , L 1
-) be a stochastic operator and let Assumptions 4.4 be satisfied. Then, the C 0 -semigroup (U H (t)) t 0 is irreducible.

Proof. Let h ∈ X and g ∈ L ∞ (Ω × V, dx ⊗ m(dv)) be nonnegative and non trivial. Denoting the duality parity between X and its dual simply by •, • , we have for any λ > 0

R(λ, A)h, g = R λ h, g + ∞ k=0 Ξ λ H (M λ H) k G λ g, g = R λ h, g + ∞ k=0 (M λ H) k G λ h, H Ξ λ g
where H and Ξ λ denote the dual operator of H and Ξ λ respectively. Notice that g = Ξ λ g is nonnegative and nontrivial and the same holds for H Ξ λ = H g since, under Assumption 4.4 ψ 0 , H g = Hψ 0 , g > 0.

Now, since the irreducibility of M 0 H is equivalent to that of M λ H for any λ > 0, we deduce that R(λ, A)h, g > 0 for any nonnegative and nontrivial g ∈ L ∞ (Ω × V, dx ⊗ m(dv)) which proves that R(λ, A)h is positive a. e. on Ω × V and R(λ, A) is positivity improving.

The main result of this section is then the following:

Theorem 4.7. Let H ∈ B(L 1 + , L 1
-) be a stochastic boundary operator and let Assumptions 4.4 be satisfied. Assume there exists Ψ H ∈ D(T H ) (with unit norm) such that

T H Ψ H = 0
and Ψ H 0 a.e. on Ω × V.

Then, T H generates a irreducible and stochastic C 0 -semigroup (U H (t)) t 0 on X and Ψ H is the unique invariant density of (U H (t)) t 0 . Moreover, (U H (t)) t 0 is ergodic with ergodic projection

Pf = f Ψ H , where f = Ω×R d f (x, v)dx ⊗ m(dv), f ∈ X, i.e. lim t→∞ 1 t t 0 U H (s)f ds -f Ψ H X = 0, ∀f ∈ X (4.2)
and X = Ker(T H ) ⊕ Range(T H ).

Proof. According to Proposition 4.1, there exists ϕ ∈ L 1 + such that ϕ 0 µ-a.e. and ϕ = M 0 Hϕ. Since M 0 H is irreducible, ϕ > 0 µ-a.e. and we deduce from Proposition 4.1 that A = T H . Moreover, Lemma 4.6 ensures that (U H (t)) t 0 is irreducible. Since Ker(T H ) = 0, the ergodicity of (U H (t)) t 0 follows from [17, Theorem 7.3, p. 174 and Theorem 5.1 p. 123].

MAIN COMPACTNESS RESULT AND EXISTENCE OF AN INVARIANT DENSITY

5.1. Weak compactness result. We prove here the main compactness result of the paper. The proof of the result is based on a series of important geometrical results regarding regularity and transversality of the ballistic flow

ξ : (x, v) ∈ Γ + → ξ(x, v) = (x -τ -(x, v)v, v) ∈ ∂Ω × V.
Such highly technical results have been postponed to Appendix A for the clarity of the reading and will be used repeatedly in the proof of the following.

Theorem 5.1. Let K ∈ B(L 1 + , L 1
-) be regular diffuse boundary operator and let R ∈ B(L 1 + , L 1 -) be a pure reflection operator. Then, for any integers k, 1, one has

K(M 0 R) k M 0 K(M 0 R) M 0 K ∈ B(L 1 + , L 1 -) is weakly compact. Proof. Let k, 1 be fixed. Let (K m ) m ⊂ B(L 1 + , L 1
-) be the sequence of approximation obtained from Lemma 3.7, which is such that

lim m K -K m B(L 1 + ,L 1 -) = 0. It is then enough to prove the weak-compactness of K n (M 0 R) k M 0 K m (M 0 R) M 0 K p for any n, m, p ∈ N. Still using the notations of Lemma 3.7, introduce K m f (x, v) = ψ m (v) v •n(x)>0 f (x, v )µ x (dv ), f ∈ L 1 + , (x, v) ∈ Γ -, m ∈ N
where ψ m (v) = m1 Bm (v). Given n, m, p ∈ N, using (3.5) and a domination argument, the weakcompactness of K n (M 0 R) k M 0 K m (M 0 R) M 0 K p would imply the result. To avoid too heavy notations, and setting for instance F (v) = max(ψ n , ψ m , ψ p ), it suffices to prove that

K(M 0 R) k M 0 K(M 0 R) M 0 K is weakly-compact for Kϕ(x, v) = F (v) Γ + (x) ϕ(x, v )µ x (dv ). (5.1)
Since F is compactly supported and bounded we can assume without loss of generality that

F = 1 B 1 (5.2)
which amounts to consider only velocities

|v| 1. Write Rϕ(x, v) = ϕ(x, V(x, v)), (x, v) ∈ Γ -, ϕ ∈ L 1 +
, where V satisfies Definition 3.1. Recalling the transformation (see Appendix A)

U : (x, v) ∈ Γ + → U(x, v) = (ξ s (x, v); V(ξ(x, v))) ∈ Γ +
which preserves the measure µ where here, as well as in Appendix A, ξ is the ballistic flow and we adopt the following notations: for any element z = (x, v) ∈ Ω × V , we write x = z s and v = z v .

Consequently, for any j ∈ N, the iterated U j is also µ-preserving. Notice that, for any ϕ ∈ L 1 + and any j ∈ N, (M 0 R) j ϕ(x, v) = ϕ(U j (x, v)), and

(M 0 R) j M 0 ϕ(x, v) = ϕ(ξ • U j (x, v)), (x, v) ∈ Γ + . (5.3)
For simplicity, we write

Z j = (M 0 R) j M 0 ∈ B(L 1 -, L 1 + ). It is clear that Z j is stochastic for any j ∈ N and K(M 0 R) k M 0 K(M 0 R) M 0 K = KZ k KZ K.
Using the notations introduced in the proof of Proposition A.10,

U -• ξ -1 : Γ -\ γ → Γ ( +1) + ⊂ Γ + is a C 1 diffeomorphism from Γ -\ γ onto its image where Γ ( +1) + = Γ + \ +1
j=1 Λ j and γ are of zero µ-measure. Denote here for simplicity

Γ ( ) -= Γ -\ γ .
We may make the identification

L 1 + = L 1 (Γ ( +1) + , dµ), L 1 -= L 1 (Γ ( ) -, dµ) (5.4)
so that we only have to prove the weak-compactness of

KZ k KZ K : L 1 (Γ ( +1) + , dµ) -→ L 1 (Γ ( ) -, dµ).
Notice that, by (5.1) and (5.2), the range of KZ k KZ K can be rather considered as

L 1 (Γ ( ) -, dμ) where dμ(x, v) = F (v)dµ(x, v) = 1 B 1 (v)|v • n(x)|π(dx) ⊗ m(dv) (5.5)
is nothing but the restriction of µ -to ∂Ω × B 1 . In particular, μ is a finite measure. From a simple consequence of the Dunford-Pettis criterion (see [START_REF] Bogachev | Measure Theory[END_REF]Corollary 4.7.21, p. 288]), we need to prove that, for any nonincreasing sequence of measurable subsets (A j ) j ⊂ Γ -with j A j = ∅, it holds

lim j→∞ sup ϕ L 1 + 1 A j |KZ k KZ Kϕ(x, v)| dµ(x, v) = 0. (5.6)
Since K and Z j are nonnegative operators, it suffices of course to consider nonnegative ϕ ∈ L 1 + in (5.6). Let us fix a sequence (A j ) j ⊂ Γ -with j A j = ∅ and consider a nonnegative ϕ ∈ L 1 + . From now on, the identification (5.4) will be assumed and -for notation simplicity -we shall always omit the superscript in Γ ± , which means that, in the all the sequel, Γ ± = Γ ( ) ± and also

Γ ± (x) = Γ ( ) ± (x) = {v ∈ R d ; (x, v) ∈ Γ ( )
± } for any x ∈ ∂Ω. We set

I j (ϕ) = A KZ k KZ Kϕ(x, v)dµ(x, v).
Introduce then, for n ∈ N, the µ-preserving change of variables

(x, v) ∈ Γ + -→(y, w) = ξ • U n (x, v ) = (Y n (x, v), V n (x, v )) ∈ Γ -,
where Y n and V n are respectively the "position" and "velocity" components, i.e. Y n (x, v) ∈ ∂Ω with V n (x, v)•n(Y n (x, v)) < 0. We denote then the position and velocity components of the inverse

(ξ • U n ) -1 : Γ -→ Γ + by Y n , V n , i.e. given (y, w) ∈ Γ -, U -n • ξ -1 (y, w) = (ξ • U n ) -1 (y, w) = (Y n (y, w), V n (y, w)) ∈ Γ + .
Defining, for any x, y ∈ ∂Ω,

B j,k (y) = Γ -(y) G j (Y k (y, w))F (w)µ y (dw), G j (x) = Γ -(x) 1 A j (x, v)F (v)µ x (dv),
one can show (see the Lemma 5.5 hereafter) that

I j (ϕ) = ∂Ω π(dz) Γ + (z) F (u)B j,k (Y (z, u))µ z (du) Γ + (z) ϕ(z, w 2 )µ z (dw 2 ).
It is clear then that (5.6) will hold true if

lim j→∞ sup z∈∂Ω Γ + (z) F (u)B j,k (Y (z, u))µ z (du) = 0. (5.7)
The proof of this property is given in the next Lemma yielding the desired weak-compactness.

Lemma 5.2. With the notations of the proof of Theorem 5.1, given z ∈ ∂Ω, introduce

P j (z) = Γ + (z) F (u)B j,k (Y (z, u))µ z (du).
Then, lim

ψ ∈ L 1 (R d , m) it holds R d ψ(v)m(dv) = 1 |S d-1 | ∞ 0 m 0 (d ) S d-1 ψ( ω)σ(dω)
where dσ denotes the Lebesgue measure on S d-1 with total mass |S d-1 |.

Remark 5.4. We shall use in the proof of Lemma 5.2 that, with the notations of Proposition A.8, for any y ∈ ∂Ω we can construct an orthonormal basis {e 1 (y), . . . , e d (y)} of R d depending continuously on y ∈ ∂Ω with e d (y) = -n(y)

and such that, in this basis, any ∈ S d-1 can be written as

= d i=1 ω i e i (y)
where ω = ω(θ) = (ω 1 , . . . , ω d ) is given by (A.14) in terms of the polar coordinates

θ = (θ 1 , . . . , θ d-1 ) ∈ U = [0, 2π] × [0, π] d-3 × 0, π 2 .
In this case, ω is independent of y ∈ ∂Ω. We also recall that, according to Remark A.9, for any ε > 0, one can define S ε (y) as those ∈ S d-1 for which

θ ∈ U ε := (θ 1 , . . . , θ d-1 ) ∈ U ; sin d-2 θ d-1 sin d-3 θ d-2 . . . sin θ 2 ε .
and prove that lim ε→0 + sup y∈∂Ω σ(S ε (y)) = 0. See Remark A.9 for more details.

Proof of Lemma 5.2. Recall that (5.4) is in force and (A j ) j ⊂ Γ -= Γ -is non-increasing with j A j = ∅. Recall also that

P j (z) = Γ + (z) F (u)µ z (du) Γ -(Y (z,u)) G j (Y k (Y (z, u), w))F (w)µ Y (z,u) (dw) = Γ + (z) F (u)µ z (du) Γ -(Y (z,u)) F (w)µ Y (z,u) (dw)× × Γ -(Y k (Y (z,u),w)) 1 A j (Y k (Y (z, u), w), v)F (v)µ Y k (Y (z,u),w) (dv)
where F = 1 B 1 . Therefore,

P j (z) = Γ + (z) F (u)µ z (du) Γ -(Y (z,u)) F (w)G j (Y k (Y (z, u), w))µ Y (z,u) (dw).
Introducing the polar coordinates w = , u = rν with r, > 0 and , ν ∈ S d-1 and using Lemma 5.3 (recall that µ x (dv) = |v • n(x)|m(dv)) we get

P j (z) = |S d-1 | -2 1 0 rm 0 (dr) Γ + (z) |ν • n(z)|σ(dν) × 1 0 m 0 (d ) Γ -(Y (z,rν)) G j (Y k (Y (z, rν), ))| • n(Y (z, rν))|σ(d ).
Since {0} is not an atom for the measure m 0 (d ), according to the dominated convergence theorem, it is enough to prove that, for any given , r ∈ (0, 1) and ν ∈ S d-1 ,

lim j→∞ sup z∈∂Ω Γ -(Y (z,rν)) G j (Y k (Y (z, rν), ))| • n(Y (z, rν))|σ(d ) = 0.
It suffices then to prove that

lim j→∞ sup y∈∂Ω Γ -(y) G j (Y k (y, ))| • n(y)|σ(d ) = 0 (5.8)
and, since > 0, there is no loss of generality in proving the result only for = 1. Notice that, for any ε > 0 and y ∈ ∂Ω,

Γ -(y) G j (Y k (y, ))| • n(y)|σ(d ) = Γ -(y)∩Sε(y) G j (Y k (y, ))| • n(y)|σ(d ) + Γ -(y)\Sε(y) G j (Y k (y, ))| • n(y)|σ(d )
where S ε (y) has been introduced in the above Remark 5.4. Clearly, since

G j ∞ R d F (u)|u|m(du), there is C > 0 independent of j such that sup y∈∂Ω Γ -(y)∩Sε(y) G j (Y k (y, ))| • n(y)|σ(d ) C sup y∈∂Ω σ (S ε (y))
which goes to 0 as ε → 0 according to Remark A.9. Therefore, to show (5.8), we only have to prove that, for any ε > 0,

lim j→∞ sup y∈∂Ω Γ -(y)\Sε(y) G j (Y k (y, ))σ(d ) = 0.
Recall that, with the notations of Section A.2,

U -k • ξ -1 : Γ (k) -= Γ -\ γ k -→ Γ (k+1) + , Y k (y, ) = U -k s • ξ -1 (y, ) ∀y ∈ ∂Ω, ∈ Γ -(y),
and

Γ -(y) = { ∈ S d-1 ; (y, ) ∈ Γ (k)
-}. Thus, we have to show that, for any ε > 0,

lim j→∞ sup y∈∂Ω Γ -(y)\Sε(y) G j U -k s • ξ -1 (y, ) σ(d ) = 0.
(5.9)

It follows from Lemma A.11 that, for any y ∈ ∂Ω, the mapping

∈ Γ -(y) -→ U -k s • ξ -1 (y, ) ∈ ∂Ω is of class C 1 with differential of rank d -1.
(5.10) On the other hand, with the notations and parametrization used in Proposition A.8 and recalled in Remark 5.4, the mapping (θ, y) ∈ U × ∂Ω -→ = (θ, y) is continuous while

θ ∈ U -→ = (θ, y) is of class C 1 with a continuous derivative (θ, y) ∈ A × ∂Ω -→ ∂ θ (θ, y). Since the mapping θ ∈ U ε → (θ, y) is a C 1 parametrization of S ε (y)
, by virtue of (5.10) we have that, for any

y ∈ ∂Ω, θ ∈ U \ U ε -→ U -k s • ξ -1 (y, (θ, y)) ∈ ∂Ω is a regular parametrization of E y := U -k s • ξ -1 (y, ) ; ∈ Γ -(y) \ S ε (y) ⊂ ∂Ω. Then,
J Y k (y, )dθ 1 . . . dθ d-1 = J Y k (y, )dθ where J Y k (y, ) = det (∂ θ i Y k (y, ) , ∂ θ Y k (y, )) 1 i, d-1 1/2 > 0 on U \ U ε . Since the mapping (θ, y) ∈ U × ∂Ω -→ ∂ θ i Y k (y, (θ, y)) = d Y k (y, (θ, y))∂ θ i (θ, y))
is continuous for any i = 1, . . . , d -1, then so is the mapping

(θ, y) ∈ U × ∂Ω -→ J Y k (y, (θ, y))
and there exists C ε > 0 such that

J Y k (y, (θ, y)) C ε > 0, ∀(θ, y) ∈ (U \ U ε ) × ∂Ω.
Hence, for any y ∈ ∂Ω

Γ -(y)\Sε(y) G j (U -k s • ξ -1 (y, ))σ(d ) 1 C ε U \Uε G j U -k s • ξ -1 (y, (θ, y)) J Y k (y, (θ, y))dθ = 1 C ε Ey G j (Y )π Ey (dY ).
Clearly, recalling the definition of G j -and because the measures π and π Ey coincide on E y -we get

Ey G j (Y )π Ey (dY ) = Ey Γ -(Y ) F (v)1 A j (Y, v)µ Y (dv) π Ey (dY ) ∂Ω Γ -(Y ) F (v)1 A j (Y, v)µ Y (dv) π(dY ) = μ(A j ),
where μ is given by (5.5). Thus,

sup y∈∂Ω Γ -(y)\Sε(y) G j (U -k s • ξ -1 (y, ))σ(d ) μ(A j ) C ε ∀j ∈ N.
Since (A j ) j is non-increasing with j A j = ∅ and μ is a finite measure, we have lim j μ(A j ) = 0 which implies (5.9) and proves the Lemma.

Lemma 5.5. With the notations of the proof of Theorem 5.1, it holds, for any j ∈ N,

I j (ϕ) = ∂Ω π(dz) Γ + (z) F (u)B j,k (Y (z, u))µ z (du) Γ + (z) ϕ(z, w 2 )µ z (dw 2 ) (5.11)
where, for any x, y ∈ ∂Ω

B j,k (y) = Γ -(y) G j (Y k (y, w))F (w)µ y (dw), G j (x) = Γ -(x) 1 A j (x, v)F (v)µ x (dv).
Proof. We use the notations of the above proof and recall that (5.4) is in force. In particular, we assume K to be given by (5.1). Notice that, for any nonnegative ψ ∈ L 1 -,

A j KZ k ψ(x, v)dµ -(x, v) = A j F (v) Γ + (x) Z k ψ(x, v )µ x (dv ) dµ -(x, v) =: A j F (v)Ψ k (x)dµ -(x, v)
where, using (5.3), one has

Ψ k (x) = Γ + (x) Z k ψ(x, v )µ x (dv ) = Γ + (x) ψ(ξ • U k (x, v ))µ x (dv ), x ∈ ∂Ω.
Simple use of Fubini's theorem yields

A j F (v)Ψ k (x)dµ -(x, v) = ∂Ω G j (x)Ψ k (x)π(dx) = Γ + G j (x)ψ(ξ • U k (x, v ))dµ + (x, v ).
Introduce then the µ-preserving change of variables

(x, v ) ∈ Γ + →(y, w) = ξ • U k (x, v ) = (Y k (x, v ), V k (x, v )) ∈ Γ - (5.12)
and noticing that x = Y k (y, ω) we have,

A j KZ k ψ(x, v)dµ -(x, v) = Γ - G j (Y k (y, w))ψ(y, w)dµ -(y, w).
Applying this with ψ = Kg for some nonnegative g ∈ L 1 + , we get

A j KZ k Kg(x, v)dµ -(x, v) = Γ - G j (Y k (y, w))Kg(y, w)dµ -(y, w) Γ - G j (Y k (y, w))F (w)dµ -(y, w) Γ + (y)
g(y, w 1 )µ y (dw 1 )

and, with g = Z Kϕ for some nonnegative ϕ ∈ L 1 + , noticing that g(y, w 1 ) = F (V (y, w 1 ))

Γ + (Y (y,w 1 )) ϕ(Y (y, w 1 ), w 2 )µ Y (y,w 1 ) (dw 2 )
where we used the notations introduced in (5.12), we get

I j (ϕ) := A j KZ k KZ Kϕ(x, v)dµ -(x, v) = A j KZ k Kg(x, v)dµ -(x, v) = Γ - G j (Y k (y, w))F (w)dµ -(y, w)× × Γ + (y) Γ + (Y (y,w 1 ))
ϕ(Y (y, w 1 ), w 2 )µ Y (y,w 1 ) (dw 2 ) F (V (y, w 1 ))µ y (dw 1 ).

Thanks to Fubini's Theorem, this reads

I j (ϕ) = Γ + F (V (y, w 1 ))B j,k (y)dµ + (y, w 1 ) Γ + (Y (y,w 1 )) ϕ(Y (y, w 1 ), w 2 )µ Y (y,w 1 ) (dw 2 ) with B j,k (y) = Γ -(y)
G j (Y k (y, w))F (w)µ y (dw), y ∈ ∂Ω.

We use the µ-preserving change of variables in (5.12), namely

(y, w 1 ) ∈ Γ + → (z, u) = (Y (y, w 1 ), V (y, w 1 )) = ξ • U (y, w 1 ) ∈ Γ - so that y = Y (z, u), w 1 = V (z, u)
and we obtain

I j (ϕ) = Γ - F (u)B j,k (Y (z, u))dµ -(z, u) Γ + (z) ϕ(z, w 2 )µ z (dw 2 )
which is exactly (5.11).

5.2. About the essential spectral radius of M 0 H. A consequence of the above is the following stability result for the essential spectral radius of M 0 H:

Theorem 5.6. Let H ∈ B(L 1 + , L 1
-) be a stochastic regular partly diffuse boundary operator given by (3.8). Then,

r ess (M 0 H) = r ess (M 0 α(•)R). In particular, r ess (M 0 H) sup x∈∂Ω α(x).
Proof. For notation simplicity, we simply denote by αR the operator α(•)R and by β K the operator

(1 -α(•))K (i.e. β(•) = 1 -α(•)). Let λ ∈ C with |λ| > 1. Define K = (M 0 β K)R(λ, M 0 αR). One has K = λ -1 ∞ j=0 M 0 β K λ -1 M 0 αR j and K 2 (M 0 β K) = λ -2 k, 0 λ -k-M 0 β K (M 0 αR) k M 0 β K (M 0 αR) M 0 β K.
Notice that, for any k, 0, the operator

M 0 β K (M 0 αR) k M 0 β K (M 0 αR) M 0 β K is dominated by M 0 K (M 0 R) k M 0 K (M 0 R) M 0 K
and, from Theorem 5.1, this last operator is weaklycompact. This shows that

K 2 (M 0 β K) : L 1 + → L 1 +
is a weakly-compact operator. This results in the weak-compactness of the operator K 3 and, as well-known, (M 0 β KR(λ, M 0 αR)) 6 is compact. Using the terminology of [START_REF] Voigt | A perturbation theorem for the essential spectral radius of strongly continuous semigroups[END_REF], we have M 0 β K is (M 0 αR)-power compact and this is enough to ensure that r ess (M 0 H) = r ess (M 0 α(•)R) (see [START_REF] Voigt | A perturbation theorem for the essential spectral radius of strongly continuous semigroups[END_REF]Corollary 1.4]). The last estimate is obvious since R is stochastic and

r ess (M 0 α(•)R) M 0 α(•)R B(L 1 + ) .
Remark 5.7. For regular diffuse boundary operator (i.e. H = K or, equivalently, α = 0), we have r ess (M 0 H) = 0 and M 0 H is power-compact.

KINETIC SEMIGROUP FOR REGULAR PARTLY DIFFUSE BOUNDARY OPERATORS

We introduce the following set of Assumptions:

Assumption 6.1. The regular stochastic partly diffuse boundary operator

H = α(•)R + (1 -α(•))K is such that A1) Range(K) ⊂ Y - 1 ; A2) R(Y + 1 ) ⊂ Y - 1 ; A3) sup x∈∂Ω α(x) < 1.
Remark 6.2. In the above set of Assumptions, it is possible to replace Y ± 1 with L 1 (Γ ± , τ ∓ dµ ± ). However, in this case, Assumption A2) is not necessarily satisfied for practical examples of boundary conditions (see Example 6.4). Example 6.3. Consider the classical Maxwell diffuse boundary condition for which

[Hf ](x, v) = M(v) γ(x) v •n(x)>0 f (x, v )|v • n(x)|dv , ∀(x, v) ∈ Γ - with M(v) = 1 (2πθ) d/2 exp - |v| 2 2θ
, and γ(x)

= u•n(x)<0 M(u)|u•n(x)|du, ∀x ∈ ∂Ω
for some θ > 0. Notice that, actually, γ is independent of x and

γ(x) = γ d := C d R d |v|M(v)dv, ∀x ∈ ∂Ω for some universal constant C d = |S d-2 | |S d-1 | 1 0 t(1 -t 2 ) d-3 2 dt. One has then H(L 1 + ) ⊂ Y - 1 . Indeed, for f ∈ L 1 + nonnegative, one has Hf Y - 1 = 1 γ d ∂Ω π(dx) v•n(x)<0 |v| -1 M(v)|v • n(x)|dv × v •n(x)>0 f (x, v )|v • n(x)|dv 1 γ d ∂Ω π(dx) v •n(x)>0 f (x, v )|v • n(x)|dv = 1 2γ d Γ + f dµ
where we used that R d M(v)dv = 1 for the first inequality. This shows that Range H ⊂ Y - 1 . This result extends easily to the case in which the temperature θ depends on x ∈ ∂Ω, i.e. θ = θ(x) with θ(x) θ 0 > 0 for any x ∈ ∂Ω. Example 6.4. Recalling that both V and the measure m are invariant under the orthogonal group let us consider the pure reflection boundary operator

Rϕ(x, v) = ϕ(x, v -2(v • n(x))n(x)), (x, v) ∈ Γ -, ϕ ∈ L 1 + . and let ϕ ∈ Y + 1 . Then, with the change of variables w = v -2(v • n(x))n(x) such that |w| = |v| and v = w -2(w • n(x))n(x) (which preserves the measure dµ ± ) we get Rϕ Y - 1 = Γ - ϕ(x, v -2(v • n(x))n(x)) |v| -1 |v • n(x)|m(dv)π(dx) = Γ + ϕ(x, w)|w| -1 |w • n(x)|π(dx)m(dw), i.e. R(Y + 1 ) ⊂ Y - 1 .
Notice that, in this example, in full generality, we can not replace Y -

1 with L 1 (Γ -, τ + dµ -). Actually, requiring that R(L 1 (Γ + , τ -dµ + )) ⊂ L 1 (Γ -, τ + dµ -) amounts to assume that there exists c > 0 such that τ + (x, V(x, v)) cτ -(x, v) for any (x, v) ∈ Γ + which is a geometrical condition not satisfied if Ω is not strictly convex.
A key point is that, under Assumptions 6.1, the following holds:

Lemma 6.5. Assume H ∈ B(L 1 + , L 1 -) satisfies Assumptions 6.1. Then, for any ϕ ∈ L 1 + ψ = R(1, M 0 (αR))M 0 ((1 -α)K)ϕ ∈ Y + 1 so that Hψ ∈ Y - 1 . Proof. Notice that, since sup x∈∂Ω α(x) = α 0 < 1, one has M 0 (αR) B(L 1 + ) α 0 < 1 and R(1, M 0 (αR)) is well-defined. From Assumption 6.1 A1), (1-α(•))Kϕ ∈ Y - 1 . Then, from (2.10), g = M 0 (1 -α)Kϕ ∈ Y + 1 . From Assumption 6.1 A2), αRg ∈ Y - 1 and, from (2.10), M 0 αRg ∈ Y + 1 . More precisely, M 0 α R B(Y + 1 ) α 0 < 1 so that ψ = ∞ n=0 (M 0 (αR)) n g ∈ Y + 1 . Now, it is clear that Hψ ∈ Y - 1 since H(Y + 1 ) ⊂ Y - 1 (notice that K maps any function in Y - 1 while R(Y + 1 ) ⊂ Y - 1 )
). We can now state our main existence and uniqueness result about invariant density:

Theorem 6.6. Let H ∈ B(L 1 + , L 1
-) be a regular stochastic partly diffuse boundary operator and let Assumptions 6.1 and 4.4 be satisfied. Then,(T H ,D(T H )) is the generator of a stochastic C 0semigroup (U H (t)) t 0 . Moreover, (U H (t)) t 0 is irreducible and has a unique invariant density

Ψ H ∈ D(T H ) with Ψ H (x, v) > 0 for a. e. (x, v) ∈ Ω × R d , Ψ H X = 1
and Ker(T H ) = Span(Ψ H ). Moreover, (U H (t)) t 0 is ergodic, Eq. (4.2) holds and X = Ker(T H )⊕ Range(T H ).

Proof. We begin with proving that, under Assumptions 6.1, 1 ∈ S p (M 0 H). Indeed, being both M 0 and H stochastic, the spectral radius of M 0 H is r σ (M 0 H) = 1 ∈ S(M 0 H). According to Theorem 5.6, one also has r ess (M 0 H) = r ess (M 0 (αR)) sup x∈∂Ω α(x) < 1, i.e.

r ess (M 0 H) < 1 = r σ (M 0 H).
As well-known (see [

in Y + 1 . Being M 0 Hϕ = ϕ, we have ϕ = M 0 (αR)ϕ + M 0 (1 -α)Kϕ so that, since 1 -M 0 (αR) is invertible, ϕ = R(1, M 0 (αR))M 0 ((1 -α)K)ϕ.
From Lemma 6.5, we get that ϕ ∈ Y + 1 . We deduce then from Proposition 4.2 that there exists Ψ H ∈ D(T H ) nonnegative and such that T H Ψ H = 0. We conclude with Theorem 4.7.

Remark 6.7. The fact that T H is the generator of (U H (t)) t 0 does not depend on A1) and A2) in Assumptions 6.1.

ASYMPTOTIC STABILITY OF COLLISIONLESS KINETIC SEMIGROUPS

The object of this section is to complement Theorem 4.7 and Theorem 6.6 where a convergence in Cesarò means of (U H (t)) t 0 to its ergodic projection is given. Indeed, under a quite weak additional assumption on the kernel of K we will show that (U H (t)) t 0 is asymptotically stable, i.e. U H (t)f converges in norm as t → +∞. In particular

lim t→∞ U H (t)f -Ψ H X = 0 (7.1)
for any density f ∈ X, i.e. any nonnegative f with f X = 1. For the sake of simplicity, we restrict ourselves to the case in which m(dv) = dv is the Lebesgue measure over

V = v ∈ R d ; m |v| M
where 0 m < M ∞, although the surface Lebesgue measure on the unit sphere can also be dealt with, see Remark 8.2 below.

In order to prove asymptotic stability of (U H (t)) t 0 we first describe the movement of particles as a piecewise deterministic Markov process. Then we explain how the stochastic semigroup (U H (t)) t 0 can be defined by this process and finally we prove the asymptotic stability of this semigroup.

7.1. Piecewise deterministic Markov process. Consider the following stochastic process which describes the movement of particles. A particle is moving in the space Ω with a constant velocity and when it strikes the boundary ∂Ω a new direction is chosen randomly from the directions that point back into the interior of Ω and the motion continues. We recall that if x ∈ ∂Ω and v ∈ Γ + (x) then the distribution of velocity v after reflection is given by a probability measure P (x, v , •) defined on Borel subsets B of Γ -(x) by

P (x, v , B) = α(x)δ V(x,v ) (B) + (1 -α(x)) B h(x, v, v ) m(dv).
where V(x, v) is the regular reflection law. From the Assumptions 6.1 A3) it follows that there exists ε 0 = 1 -sup x∈∂Ω α(x) > 0 such that 1 -α(x) ε 0 for all x ∈ ∂Ω. This implies that

P (x, v , dv) ε 0 h(x, v , v) m(dv). (7.2)
Let a particle starts at time t = 0 from some point x ∈ Ω with some initial velocity v ∈ V \ {0} or from x ∈ ∂Ω with velocity v ∈ Γ -(x). Let x(t) be the position and v(t) be the velocity of the particle at time t and let

t 1 = t + (x, v). Then x(t) = x + vt and v(t) = v for t ∈ [0, t 1 ). Let 0 < t 1 < t 2 < .
. . be a sequence of times when a particle hits the boundary ∂Ω. Then

Prob(v(t n ) ∈ B | x(t - n ) = x, v(t - n ) = v ) = P (x, v , B) for every Borel subset B of Γ + (x)
, where x(t - n ) and v(t - n ) are the left-hand side limits of x(t) and v(t), respectively, at the point t n . Moreover

x(t) = x(t n ) + v(t n )(t -t n ), v(t) = v(t n ) for t ∈ [t n , t n+1 ), x(t n ) = x(t - n ) and t n+1 = t n + t + (x(t n ), v(t n )) for n 1.
It is easy to observe that ξ(t) = (x(t), v(t)), t 0, defines a piecewise deterministic Markov process [START_REF] Rudnicki | Piecewise Deterministic Processes in Biological Models[END_REF] with values in the space

E = (Ω × V ) ∪ Γ -.
The process {ξ(t)} t 0 has càdlàg sample paths, i.e., they are right-continuous with left limits. Let P(t, x, v, B) be the transition probability function for this process, i.e.

P(t, x, v, B) = Prob(ξ(t) ∈ B | ξ(0) = (x, v)),
where B are Borel subsets of E. The semigroup (U H (t)) t 0 can be uniquely determined by the transition probability function P(t, x, v, B) because the following relation holds

B U H (t)f (y, w) dy ⊗ m(dw) = Ω×V P(t, x, v, B)f (x, v)dx ⊗ m(dv)
for all f ∈ X, Borel subsets B of Ω × V and t 0.

Remark 7.1. It should be noted that we do not assume here that Ω is a strictly convex set and it can happen that at some boundary points x some outward or inward vectors belong to the tangent space Γ 0 (x). In such cases trajectories can be tangent to the boundary ∂Ω, especially in the case when we consider the specular reflection (see Fig. 3). But there is no need to consider such pathological trajectories because the set Γ 0 has zero measure for µ and does not play any role in the definition of the boundary operator H .

7.2. Asymptotic stability. Now we check that the semigroup (U H (t)) t 0 is partially integral (see Appendix B.1 for precise definition), i.e. that for some t > 0 there there exists an integrable function q

: Ω × V × Ω × V → [0, ∞), q ≡ 0, such that P(t, x, v, B) B q(x, v, y, w) dy ⊗ m(dw). (7.3)
In order to prove this property we need a rather weak assumption concerning function h(x, v, v ).

Definition 7.2. Let H ∈ B(L 1 + , L 1
-) be a stochastic partly diffuse boundary operator of the form (3.8). We say that the boundary operator H is weakly locally diffuse (WLD) if for each point x ∈ ∂Ω and v 0 ∈ Γ + (x) there exists a v 0 ∈ Γ -(x) and δ > 0 such that

k(x, v, v ) > 0 for µ x -a. e. v ∈ Γ -(x) ∩ B(v 0 , δ), v ∈ Γ + (x) ∩ B(v 0 , δ). (7.4)
If we replace condition (7.4) by a stronger one:

k(x, v, v ) δ for all v ∈ Γ -(x) ∩ B(v 0 , δ), v ∈ Γ + (x) ∩ B(v 0 , δ). (7.5)
then the boundary operator H will be called strongly locally diffuse (SLD).

Lemma 7.3. Assume that the operator H is weakly locally diffuse and satisfies Assumptions 4.4 and 6.1. Then the semigroup (U H (t)) t 0 is partially integral.

Proof. Let (x, v) ∈ E 0 = Ω × (V \ {0}) be the initial position and velocity of a particle. At time t 1 = t + (x, v) it hits the point x 1 = x + t + (x, v)v on the boundary ∂Ω. Then we choose a new velocity v1 ∈ Γ -(x 1 ) and at time t 1 + t + (x 1 , v1 ) the particle hits the boundary for the second time at the point x2 =

x 1 + t + (x 1 , v1 )v 1 . We choose a new velocity v2 ∈ Γ -(x 2 ). Let t > 0 satisfies inequalities t 1 + t + (x 1 , v1 ) < t < t 1 + t + (x 1 , v1 ) + t + (x 2 , v2 ) (7.6)
and let τ = t -t 1 . We will find an neighborhood and v(t) = v 2 is the position at time t of the particle if it starts from x with velocity v, and after hitting the boundary we choose velocities v 1 and v 2 (see Fig. 4). We define the function F : U → V by

U of (v 1 , v2 ) ∈ V such that for (v 1 , v 2 ) ∈ U we have v 1 ∈ Γ -(x 1 ), v 2 ∈ Γ -(x 1 + t + (x 1 , v 1 )v 1 ) and (7.6) is satisfied for (v 1 , v 2 ) ∈ U . Then x(t) = x 2 + (τ -t + (x 1 , v 1 ))v 2 = x 1 + t + (x 1 , v 1 )v 1 + τ -t + (x 1 , v 1 ) v 2
F (v 1 , v 2 ) = x 1 + t + (x 1 , v 1 )v 1 + τ -t + (x 1 , v 1 ) v 2 , v 2 = (F 1 (v 1 , v 2 ), F 2 (v 1 , v 2 )) .

Now we check that if

∇ v1 t + (x 1 , v1 ) • v2 = t + (x 1 , v1 ) + ∇ v1 t + (x 1 , v1 ) • v1 (7.7)
then the function F is a local diffeomorphism in some neighborhood of (v 1 , v2 ). Indeed, let us denote the Jacobian matrix of F by J F = ∂F ∂v 1 , ∂F ∂v 2 . One checks easily that

det J F (v 1 , v 2 ) = det ∂F 1 ∂v 1 = det ∂ ∂v 1 t + (x 1 , v 1 )(v 1 -v 2 ) =: det M
where the matrix M is given by

M = c Id + A, A = a ⊗ u = a t u,
where a is the vector

a = ∂ v 1 t + (x 1 , v 1 ), u = v 1 -v 2 and c = t + (x 1 , v 1 ). We check then that * det M = c d-1 (c + a • u) = t d + (x 1 , v 1 ) + t d-1 + (x 1 , v 1 )∇ v 1 t + (x 1 , v 1 ) • (v 1 -v 2 )
. Consequently, if (7.7) holds then the Jacobian matrix J F (v 1 , v2 ) is non singular and F is a diffeomorphism in some neighborhood of (v 1 , v2 ). Observe that condition (7.7) does not hold only * Indeed, let p(z) be the characteristic polynomial of A, p(z) = det(A -zId). Since the rank of A is less or equal to 1, z = 0 is a root of p with multiplicity at least d -1 while the trace tr(A) = u • a should also be a root of p. Therefore,

p(z) = (-1) d z d-1 (z -u • a)
and, taking z = -c gives the result. on a (2d -1)-dimensional differentiable manifold in R 2d and we can change equality in (7.7) to inequality after a small perturbation of the vector v2 . We have

P(t, x, v, B) = Prob((x(t), v(t)) ∈ B) Prob((v 1 , v 2 ) ∈ U : F (v 1 , v 2 ) ∈ B) F -1 (B) ε 2 0 k(x 1 , v 1 , v)k(x 1 + v 1 t + (x 1 , v 1 ), v 2 , v 1 ) m(dv 1 ) m(dv 2 ).
Since

x 1 = x + vt + (x, v) we can define κ (x,v) (v 1 , v 2 ) := ε 2 0 k(x 1 , v 1 , v)k(x 1 + v 1 t + (x 1 , v 1 ), v 2 , v 1 ).
Since k satisfies WLD, for each (x 0 , v 0 ) ∈ E 0 there exist δ > 0 and v1 ∈ Γ

-(x 1 ), v2 ∈ Γ -(x 1 + t + (x 1 , v1 )v 1 ) such that κ (x,v) (v 1 , v 2 ) > 0 a. e. for (x, v) ∈ B((x 0 , v 0 ), δ ) and (v 1 , v 2 ) ∈ B((v 1 , v2 ), δ ).
Without lost of generality we can assume that condition (7.7) holds and F is a diffeomorphism from U 0 = B((v 1 , v2 ), δ ) onto F (U 0 ). Then P(t, x, v, B)

F -1 (B) κ (x,v) (v 1 , v 2 ) m(dv 1 ) m(dv 2 ) = B∩F (U 0 ) κ (x,v) (F -1 (y, w)) |det J F -1 (y, w)| dy ⊗ m(dw)
where J F -1 (y, w) = ∂F -1 ∂y , ∂F -1 ∂w is the Jacobian matrix of F -1 . From the last inequality it follows that (7.3) holds for q(x, v, y, w) = 1 F (U 0 ) (y, w)κ (x,v) (F -1 (y, w)) |det J F -1 (y, w)| and the semigroup (U H (t)) t 0 is partially integral. Remark 7.4. It is very likely that an analytical proof based upon the Dyson-Phillips-like representation of the semigroup (U H (t)) t 0 obtained in [START_REF] Arlotti | Explicit transport semigroup associated to abstract boundary conditions[END_REF][START_REF] Arlotti | Transport semigroup associated to positive boundary conditions of unit norm: a Dyson-Phillips approach[END_REF] may replace the adopted probabilistic proof. Such a proof seems more involved than the probabilistic one given here and we did not investigate further on this point.

Combining Theorem 4.7, Theorem B.2 and Lemma 7.3 we obtain: Theorem 7.5. Let the assumptions of Theorem 4.7 be satisfied. Assume moreover that H is weakly locally diffuse, then the semigroup (U H (t)) t 0 is asymptotically stable. Remark 7.6. In particular, under the conditions of Theorem 6.6, the semigroup (U H (t)) t 0 is asymptotically stable.

SWEEPING PROPERTIES OF COLLISIONLESS KINETIC SEMIGROUPS

The asymptotic stability of the semigroup (V H (t)) t 0 is strictly connected with the existence of an invariant density which was assumed in Theorem 4.7 and proved in Theorem 6.6. We investigate here the behaviour of (U H (t)) t 0 when this semigroup has no invariant density. A crucial role is played by sweeping property (see Appendix B.2). We first establish the following which complements Lemma 7.3: Lemma 8.1. If H is strongly locally diffuse (in the sense of Definition 7.2), then, defining E 0 = Ω × (V \ {0}), for every (x 0 , v 0 ) ∈ E 0 there exist ε > 0, t > 0, and a measurable function

η 0 : E 0 → [0, ∞) such that E 0 η(x, v) dx ⊗ m(dv) > 0 and q(x, v, y, w) η(y, w)1 B 0 (ε) (x, v) for (y, w) ∈ Ω × V , (8.1)
where B 0 (ε) is the open ball in E 0 centered at (x 0 , v 0 ) with radius ε.

Proof. The proof uses the notations introduced in the proof of Lemma 7.3. Recall that (U H (t)) t 0 is partially integral with

q(x, v, y, w) = 1 F (U 0 ) (y, w)κ (x,v) (F -1 (y, w)) |det J F -1 (y, w)| .
If the operator H is strongly locally diffuse then there there exist δ > 0 and v1 ∈ Γ

-(x 1 ), v2 ∈ Γ -(x 1 + t + (x 1 , v1 )v 1 ) such that κ (x,v) (v 1 , v 2 ) δ for all (x, v) ∈ B((x 0 , v 0 ), δ ) and (v 1 , v 2 ) ∈ B((v 1 , v2 ), δ ). Now setting η(y, w) = 1 F (U 0 ) (y, w) |det J F -1 (y, w)|
we check that η satisfies the desired properties.

Remark 8.2. We note that Lemma 8.1 and Lemma 7.3 are also true when m(dv) is a surface Lebesgue measure on a sphere but the proofs are slightly more technical. Indeed, instead of two reflections at the boundary (see Figure 4) we need one more reflection to achieve the property that the semigroup is partially integral. Proof. Since ((Ω × V ) \ E 0 ) is of zero measure for the measure dx ⊗ m(dv), we can assume that the semigroup (U H (t)) t 0 is defined on the space L 1 (E 0 , B(E 0 ), dx ⊗ m(dv)). Then, on this space, Lemma 8.1 exactly means that (U H (t)) t 0 satisfies property (K) of Theorem B.4.

Remark 8.4. For any ε > 0 and M > ε we define the set

F ε,M = {(x, v) ∈ Ω × V : ε |v| M, dist(x, ∂Ω) ε},
where dist(x, ∂Ω) = inf{|x -y| : y ∈ ∂Ω}. Since the set F ε,M is compact in the space E 0 , for every f ∈ X we have

lim t→∞ F ε,M U H (t) f (x, v) dx ⊗ m(dv) = 0. (8.2)
This result has the following probabilistic interpretation. If the semigroup (U H (t)) t 0 has no invariant density, the velocity of almost all particles converges to 0 or to ∞ or particles get close to the boundary ∂Ω when time goes to infinity.

We complement Theorem 8.3 by a more precise sweeping result:

Theorem 8.5. Assume that H ∈ B(L 1 + , L 1
-) is a regular stochastic partly diffuse operator given by (3.8) and satisfying Assumptions 4.4. Assume moreover that H is weakly locally diffuse (WLD), sup x∈∂Ω α(x) < 1 and

µ + (x, v ) ∈ Γ + ; Γ -(x) k(x, v, v )τ + (x, v)µ x (dv) = +∞ > 0. (8.3) Then lim t→∞ Ω×V 1 {|v| ε} U H (t)f (x, v) dx ⊗ m(dv) = 0, ∀ε > 0, ∀f ∈ X. (8.4)
Proof. Note first that T H is the generator of (U H (t)) t 0 (see Remark 6.7). By virtue of Theorem B.5, the proof simply consists in showing that (U H (t)) t 0 has no invariant density and in constructing a function Ψ = Ψ(x, v) such that

0 < Ψ(x, v) < ∞ a. e. on Ω × V, Ω×V Ψ(x, v)dx ⊗ m(dv) = +∞, Ω×V 1 {|v| ε} Ψ(x, v)dx ⊗ m(dv) < +∞ (ε > 0) (8.5)
and

U H (t)Ψ Ψ, ∀t 0. (8.6)
The proof will be given in several steps. First of all, according to Theorem 5.6, there exists ϕ ∈ L 1 + such that

M 0 Hϕ = ϕ, ϕ L 1 + = 1. (8.7)
Since M 0 K is irreducible then this ϕ is unique.

• First step: The function Ψ = Ξ 0 Hϕ satisfies (8.5). Indeed, one first notices that

Γ - [Kϕ](x, v)τ + (x, v)dµ -(x, v) = ∂Ω dπ(x) Γ -(x) Γ + (x) k(x, v, v )ϕ(x, v )µ x (dv ) τ + (x, v)µ x (dv) = ∂Ω dπ(x) Γ + (x) Γ -(x) k(x, v, v )τ + (x, v)µ x (dv) ϕ(x, v )µ x (dv ) = Γ + Γ -(x) k(x, v, v )τ + (x, v)µ x (dv) ϕ(x, v )dµ + (x, v ).
Therefore, under assumption (8.3) and because ϕ(x, v ) > 0 a.e. on Γ + , we have

Γ - [Kϕ] (x, v)τ + (x, v)dµ -(x, v) = +∞.
Using Lemma 2.8 -identity (2.9) -and with

ε 0 = 1 -sup x∈∂Ω α(x) > 0 Ω×V Ψ(x, v)dx ⊗ m(dv) = Γ - [Hϕ] (x, v)τ + (x, v)dµ -(x, v) ε 0 Γ - [Kϕ] (x, v)τ + (x, v)dµ -(x, v) = +∞.
Hence,

Γ - [Hϕ] (x, v) |v| -1 dµ -(x, v) = ∞ since τ -(x, v) |v| -1
D (where we recall that D is the diameter of Ω). Thus, Hϕ / ∈ Y - 1 and 0 is not an eigenvalue of T H (associated to a nonnegative eigenvalue) according to Proposition 4.2. Since T H is the generator of (U H (t)) t 0 , this means that (U H (t)) t 0 has no invariant density. Moreover,

Ω×V 1 {|v| ε} Ψ(x, v)dx ⊗ m(dv) = Γ - 1 {|v| ε} [Hϕ] (x, v)τ + (x, v)dµ -(x, v) D ε Hϕ L 1 - = D ε ϕ L 1 + .
Using that M 0 Hϕ = ϕ, one has ϕ(x, v) = [Hϕ] (x -τ -(x, v)v, v), for any (x, v) ∈ Γ + and, from the irreducibility of M 0 H, we get 0 < ϕ(x, v) < +∞ a.e. on Γ + which in turns implies that

0 < Ψ(x, v) = Ξ 0 Hϕ(x, v) = [Hϕ] (x -t -(x, v)v, v) < +∞ a.e. on Ω × V.
This proves that Ψ satisfies (8.5).

In order to prove that Ψ satisfies also (8.6) we shall resort to Lemma 3.8 and for any n ∈ N, introduce the regular diffuse operator given by H n = αR + (1 -α)K n with K n is defined as in Lemma 3.8. As before, for any n ∈ N, there exists

ϕ n ∈ L 1 + such that M 0 H n ϕ n = ϕ n , ϕ n L 1 + = 1.
• Second step: lim n ϕ n -ϕ L 1 + = 0. To prove this, we notice that ϕ n -M 0 αRϕ n = M 0 [(1 -α)K n ] ϕ n so that

ϕ n = R(1, M 0 αR) (M 0 (1 -α)K n ) ϕ n =: L n ϕ n , ∀n ∈ N.
According to Theorem 5.1 (see also the proof of Theorem 5.6), for any n ∈ N, the operator L 3 n is weakly compact so that L 6 n is compact. By virtue of Lemma 3.8, it holds that lim

n L 6 n -L 6 B(L 1 + ) = 0 where L := R(1, M 0 αR) (M 0 (1 -α)K) . Therefore, since ϕ n = L 6 n ϕ n , the sequence (ϕ n ) n is relatively compact in L 1 + and a subsequence -still denoted (ϕ n ) n -converges in L 1 + to some g ∈ L 1
+ with unit norm. It is easy to check then that g = L 6 g which implies that M 0 Hg = g. We deduce from this that g = ϕ by uniqueness. This shows finally that the whole sequence (ϕ n ) n converges to ϕ in L 1 + . • Third step: Introducing the semigroup (V Hn (t)) t 0 associated to the boundary operator

H n n ∈ N, it holds lim n V Hn (t)f -U H (t)f X = 0, ∀t 0, f ∈ X. (8.8)
Indeed, for any n ∈ N the resolvent of the generator T Hn is given by

R(λ, T Hn ) = Ξ λ H n (λ -M λ H n ) -1 G λ + R λ
and it is easy to check, using again Lemma 3.8 and Eq. (2.1) that

lim n R(λ, T Hn )f -R(λ, A)f X = 0, ∀λ > 0, f ∈ X
where we recall that (A, D(A)) is the generator of (U H (t)) t 0 . We deduce the second step from the Trotter-Kato approximation Theorem [START_REF] Davies | One-parameter Semigroups[END_REF]Theorem 3.19,p. 83].

• Fourth step. Introduce then Ψ n = M 0 H n ϕ n . According to Theorem 6.6, V Hn (t)Ψ n = Ψ n , ∀n ∈ N, t 0.

On the other hand, since lim n H n ϕ n -Hϕ L 1 -= 0, we have, for any ε > 0, lim

n 1 {|v|>ε} H n ϕ n -1 {|v|>ε} Hϕ L 1 - = 0 and also lim n Ξ 0 (1 {|v|>ε} H n ϕ n ) -Ξ 0 (1 {|v|>ε} Hϕ) L 1 + = 0 or equivalently lim n 1 {|v|>ε} Ψ n -1 {|v|>ε} Ψ L 1 + = 0. Let then Ψ ε n = 1 {|v|>ε} Ψ n , Ψ ε = 1 {|v|>ε} Ψ, n ∈ N, ε > 0, we note that V Hn (t)Ψ ε n Ψ n and 1 {|v|>ε} V Hn (t)Ψ ε n Ψ ε n for any n ∈ N, ε > 0, t > 0.
Using the Third step, we can pass to the limit in norm in this inequality as n → +∞ and get

1 {|v|>ε} U H (t)Ψ ε Ψ ε Ψ.
From the previous part of the proof, we know that lim n x n -τ -(x n , ω n )ω n = x 0 -τ -(x 0 , ω 0 )ω 0 and, since n(•) is continuous, we get

ω 0 • n(x 0 -τ -(x 0 , ω 0 )ω 0 ) = 0,
which contradicts the assumption that (x 0 , ω 0 ) ∈ Γ + . Therefore, no sequence with the above properties can exist and Γ + is open.

Lemma A.4. For any x ∈ ∂Ω, the mapping

ω ∈ Γ + (x) -→ τ -(x, w) ∈ R + is differentiable and ∇ ω τ -: (x, ω) ∈ Γ + -→ τ -(x, ω) ∈ R + (A.3)
is continuous.

Proof. As before, let us fix x 0 ∈ ∂Ω and ω 0 ∈ Γ + (x 0 ). Since the normal vector n(•) is continuous on ∂Ω, we deduce from Lemma A.3 that there exists a radius r > 0 such that

ω • n(x) > 0 and ω • n(ξ s (x, ω)) < 0 ∀(x, ω) ∈ U(x 0 ) × V(ω 0 ), (A.4) where U(x 0 ) = B(x 0 , r) ∩ ∂Ω is an open neighbourhood of x 0 and V(ω 0 ) = B(ω 0 , r) ∩ S d-1
is an open neighbourhood of ω 0 . The continuity of τ -implies that there exists t 0 > 0 such that

τ -(x, ω) t 0 > 0 ∀(x, ω) ∈ U(x 0 ) × V(ω 0 ). (A.5) Since the mapping ξ : (x, ω) ∈ Γ + → ξ(x, ω) = (x -τ -(x, ω)ω, ω) ∈ Γ -is continuous, invertible with inverse ξ -1 : (y, ω) ∈ Γ -→ ξ -1 (y, ω) = (y + τ + (y, ω)ω, ω) ∈ Γ +
and since ξ( Γ + ) ⊂ Γ -, one has ξ -1 continuous. In particular

{(x -τ -(x, ω)ω, ω), (x, ω) ∈ U(x 0 ) × V(ω 0 )}
is an open neighbourhood of ξ(x 0 , ω 0 ) and

W(z 0 ) := {x -τ -(x, ω)ω, (x, ω) ∈ U(x 0 ) × V(ω 0 )} is an open neighbourhood of z 0 = x 0 -τ -(x 0 , ω 0 )ω 0 = ξ s (x 0 , ω 0 ) ∈ ∂Ω. Since ∂Ω is of class C 1 then (up to choosing a smaller neighbourhood W(z 0 ) if necessary), there exists a C 1 bijective mapping Ψ : y ∈ (-1, 1) d-1 -→ Ψ(y) ∈ W(z 0 )
with Ψ(0) = z 0 and such that such that the range of the differential dΨ(y) has dimension d -1 for any y ∈ (-1, 1) d-1 . We introduce open pieces of ∂Ω indexed by x ∈ ∂Ω

S x = {x -τ -(x, ω)ω, ω ∈ V(ω 0 )} ⊂ W(z 0 ). Define then O x = Ψ -1 (S x ), x ∈ U(x 0 ),
one sees that, for any x ∈ ∂Ω, the mapping Ψ : y ∈ O x → Ψ(y) ∈ S x is a parametrization of S x . Namely, given (x, ω) ∈ U(x 0 ) × V(ω 0 ), there is a unique y ∈ O x such that x -τ -(x, ω)ω = Ψ(y). Thus,

τ -(x, ω) = |x -Ψ(y)| and ω = x -Ψ(y) |x -Ψ(y)| . (A.6)
In particular x 0 -τ -(x 0 , ω 0 )ω 0 = z 0 = Ψ(0), and ω 0 = x 0 -Ψ(0) |x 0 -Ψ(0)| . Introduce the C 1 mapping

H : (x, y) ∈ U(x 0 ) × (-1, 1) d-1 → H(x, y) = x -Ψ(y) |x -Ψ(y)| ∈ V(ω 0 ) (A.7)
and, for any z ∈ R d \ {0}, let P z denote the orthogonal projection on the hyperplane orthogonal to z,

P z h = h ⊥ z := h -h, z z, z = z |z| ∈ S d-1 , h ∈ R d . Notice that P z = P z for any z ∈ R d \ {0}. Because the differential of the mapping z ∈ R d \ {0} → z |z| is given by h ∈ R d → -|z| -1 P z h, h ∈ R d it follows that the differential d y H(x, y) of H is given by d y H(x, y) : h ∈ R d -→ - P x-Ψ(y) (dΨ(y)h) |x -Ψ(y)| = - P ω (dΨ(y)h) |x -Ψ(y)| , (A.8)
where ω = H(x, y) = |x -Ψ(y)| -1 (x -Ψ(y)). Note that the differential d y H(x, y) depends continuously on (x, y) ∈ U(x 0 ) × (-1, 1) d-1 . Let us assume for a while that

Rank (d y H(x 0 , 0)) = d -1. (A.9)
Then, the dimension of the range of d y H(x, y) remains of dimension d -1 for x close enough to 0.

Recalling that ω 0 = |x 0 -Ψ(0)| -1 (x 0 -Ψ(0)), we deduce from the local inverse function theorem that, in some open neighbourhood U (x 0 ) × (-δ, δ) d-1 of (x 0 , 0) and a neighbourhood V (ω 0 ) of ω 0 such that the equation

ω = H(x, y) (x, y) ∈ U (x 0 ) × (-δ, δ) d-1 is solved uniquely as y = G(x, ω), (x, ω) ∈ U (x 0 ) × V (x 0 ) where G(x, •) is a C 1 mapping on a neighbourhood V (ω 0 ) of ω 0 and the mapping (x, ω) → G(x, ω) is continuous on U (x 0 ) × V (ω 0 ). It follows that, for x ∈ U (x 0 ) the mapping ω ∈ V (ω 0 ) -→ τ -(x, ω) = |x -Ψ(y)| = |x -Ψ(G(x, ω))| is differentiable with differential d ω τ -(x, ω) given by d ω τ -(x, ω) : h ∈ R d -→ - x -Ψ(G(x, ω)), dΨ(G(x, ω))d ω G(x, ω)h |x -Ψ(G(x, ω))| . (A.10) Since d ω G(x, ω) = (d y H(x, G(x, ω))) -1 and the mapping (x, ω) → G(x, ω) is continuous then so is (x, ω) ∈ U (x 0 ) × V (ω 0 ) -→ d ω τ -(x, ω) where G(ω, •) is a C 1 mapping and (ω, x) ∈ V (ω 0 ) × U (x 0 ) → G(ω, x) is continuous. It follows that the mapping x ∈ U (x 0 ) → τ -(x, ω) = Ψ( G(ω, x)) -x is differentiable for any ω ∈ V (ω 0 ) with differential given by d x τ -(x, ω) : h ∈ T x -→ - Ψ( G(ω, x)) -x, dΨ(G(ω, x))d x G(ω, x)h -h |Ψ( G(ω, x)) -x| (A.11)
where T x is the tangent space of ∂Ω at x ∈ ∂Ω. Let us prove now the continuity of

d x τ -(•, •). Because ω + Ψ( G(ω,x))-x |Ψ( G(ω,x))-x|
= 0, differentiating with respect to x along the direction h tangent at ∂Ω at x yields

P Ψ( G(ω,x))-x dΨ( G(ω, x))d x G(ω, x)h -h |x -Ψ( G(ω, x))| = 0 i.e. P ω dΨ( G(ω, x))d x G(ω, x)h = P ω h, ∀h ∈ T x
where we used that ω is the unit vector in the direction of (Ψ( G(ω, x)) -x). This implies that

dΨ( G(ω, x))d x G(ω, x)h = P ω h + dΨ( G(ω, x)))d x G(ω, x)h, ω ω. (A.12) Since dΨ( G(ω, x))d x G(ω, x)h is a tangent vector to ∂Ω at Ψ( G(ω, x)) = ξ s (x, ω) = x - τ -(x, ω
)ω, taking the inner product of the above identity with the normal unit vector n(ξ s (x, ω)) yields dΨ( G(ω, x))d x G(ω, x)h, ω ω, n(ξ s (x, ω)) = -P ω h, n(ξ s (x, ω)) .

Inserting this into (A.12) and since ω • n(ξ s (x, ω)) = 0 we get

dΨ( G(ω, x))d x G(ω, x)h = P ω h - P ω h, n(ξ s (x, ω)) ω, n(ξ s (x, ω)) ω
which, plugged into (A.11), yields d x τ -(x, ω)h = -ω, P ω h -P ω h , n(ξ s (x, ω)) ω , n(ξ s (x, ω)) ω -h = h ; ω -P ω h ; n(ξ s (x, ω)) ω ; n(ξ s (x, ω))

i.e. The conditions θ d-1 ∈ {0, π} only means θ d-1 = 0 (recall that θ d-1 π/2) which means that ω = (0, . . . , 0, 1). Then, for d 3, the condition θ d-2 ∈ {0, π} means that ω = (0, . . . , 0, ± sin θ d-1 , cos θ d-1 ),

i.e ω belongs to some (half) unit circle of S d-1 . More generally, the condition θ d-j ∈ {0, π} for 2 j d-2 describes a unit (j -1)-dimensional (half)-spheres of S d-1 . This means that S(x) can be written as S(x) = d-2 j=2 C j where C j is a closed set with positive (d -j)-Lebesgue measure † (2 j d -2). Therefore, σ(S(x)) = 0 (where we recall that dσ is the Lebesgue surface measure over S d-1 ) and the conclusion follows.

Remark A.9. For any x ∈ ∂Ω and any ε > 0, introduce the set S ε (x) of all ω ∈ Γ -(x) ⊂ S d-1 whose polar coordinates θ = (θ 1 , . . . , θ d-1 ) in the basis {e 1 (x), . . . , e d (x)} are such that det (J θ (ω)) ε where we recall that the determinant det (J θ (ω)) is given by (A.15) with ω = ω(θ) given by (A. [START_REF] Chacon | Linear modulus of linear operator[END_REF]. Notice that det (J θ (ω)) is actually independent of x ∈ ∂Ω. Then, the surface Lebesgue measure σ(S ε (x)) of S ε (x) is given by We have then the following.

Proposition A.10. Assume that, for π-a. e. x ∈ ∂Ω, V(x, •) : Γ -(x) → Γ + (x) is a field of measurable mappings associated to a pure reflection boundary operator as in Definition 3.1 and let

U : (x, v) ∈ Γ + ∪ Γ 0 -→ U(x, v) = (x -τ -(x, v)v, V(x -τ -(x, v)v, v))
= (ξ s (x, v), V(ξ(x, v))) ∈ Γ + ∪ Γ 0 .

For any k ∈ N there exists a subset γ k ⊂ Γ -such that:

(1) γ k is a closed subset of Γ -with µ(γ k ) = 0.

(2) U -k • ξ -1 (Γ -\ γ k ) is an open subset of Γ + and

U -k • ξ -1 : Γ -\ γ k → Γ + is a C 1 diffeomorphism from Γ -\ γ k onto its image.
Proof. We first notice that, thanks to Corollary A. for each density f . If moreover (E, ρ) is a metric space and Σ = B(E) is the σ-algebra of Borel subsets of E, a partially integral semigroup {P (t)} t 0 with kernel k(t, ξ, z) is said to satisfy the property (K) on E if the following holds:

(K) for every ξ 0 ∈ E there exist ε > 0, t > 0, and a measurable function η : E → [0, ∞) such that E η(ξ)m(dξ) > 0 and k(t, ξ, z) η(z)1 B(ξ 0 ,ε) (ξ), (ξ, z) ∈ E × E, where B(ξ 0 , ε) = {ξ ∈ E : ρ(ξ, ξ 0 ) < ε}.

We have the following which is a simple consequence of a general result concerning asymptotic decomposition of stochastic semigroups (see [START_REF] Ór | Asymptotic decomposition of substochastic operators and semigroups[END_REF]Corollary 2]):

Theorem B.4. Let {P (t)} t 0 be a stochastic semigroup on L 1 (E, Σ, m), where E is a separable metric space, Σ = B(E), and m is a σ-finite measure on (E, Σ). Assume that {P (t)} t 0 has the property (K) and has no invariant density. Then {P (t)} t 0 is sweeping from all compact sets. B.3. Foguel alternative. If a stochastic semigroup has no invariant density but we are able to find a subinvariant function f * > 0, then we can precisely point out all sets having the sweeping property [START_REF] Rudnicki | On asymptotic stability and sweeping for Markov operators[END_REF]. We start with some general description.

Let a stochastic semigroup {P (t)} t 0 be given and assume that this semigroup is partially integral. If the kernel k(t, x, y) satisfies E ∞ 0 k(t, x, y) dt m(dx) > 0 y -a.e., then {P (t)} t 0 is called a pre-Harris semigroup. In particular, if a semigroup is partially integral and irreducible then it is pre-Harris semigroup. The following condition plays a crucial role in studying sweeping.

(KT): There exists a measurable function f * such that: 0 < f * < ∞ a.e., P (t)f * ≤ f * for t 0, f * / ∈ L 1 and A f * (x) m(dx) < ∞.

In (KT) we have written P (t)f * for a non-integrable function. We can use such notation because any substochastic operator P may be extended beyond the space L 1 (see [START_REF] Foguel | The Ergodic Theory of Markov Processes[END_REF] Chap. I). If f is an arbitrary non-negative measurable function, then we define P f as a pointwise limit of the sequence P f n , where (f n ) is any monotonic sequence of non-negative functions from L 1 pointwise convergent to f almost everywhere.

Theorem B.5 ([30], Corollary 3). Let {P (t)} t 0 be a pre-Harris stochastic semigroup which has no invariant density. Assume that the semigroup {P (t)} t 0 and a set A ∈ Σ satisfy condition (KT). Then the semigroup {P (t)} t 0 is sweeping with respect to A.

Definition 4 . 3 .

 43 Let (E, Σ, m) be a given σ-finite measure space. Let B ∈ B(L 1 (E, Σ, m)) be given. Then, we say that i) B is positive and write B 0, if B leaves invariant the cone of nonnegative functions of L 1 (E, Σ, m) i.e. for any h ∈ L 1 (E, Σ, m), h(s) 0 for m-a. e. s ∈ E =⇒ Bh(s) 0 for m-a. e. s ∈ E.

Lemma 5 . 3 .

 53 The proof of the above will use the the following polar decomposition theorem (see[START_REF] Voigt | Functional analytic treatment of the initial boundary value problem for collisionless gases[END_REF], Lemma 6.13, p.113]): If m is a orthogonally invariant Borel measure with support V ⊂ R d , introduce m 0 as the image of the measure m under the transformation v ∈ R d → |v| ∈ [0, ∞), i.e. m 0 (I) = m {v ∈ R d ; |v| ∈ I} for any Borel subset I ⊂ R + . Then, for any

FIGURE 3 .

 3 FIGURE 3. Examples of pathological trajectories in the case of the specular reflection v -outward vector, v -inward vector.

FIGURE 4 .

 4 FIGURE 4. Position x(t) of the particle after two collisions with the boundary.

d

  x τ -(x, ω)h = h, n(ξ s (x, ω)) ω, n(ξ s (x, ω)) (A.[START_REF] Cessenat | Théorèmes de traces pour les espaces de fonctions de la neutronique[END_REF] This gives directly the continuity of the mapping (x, ω) → d x τ -(x, ω) since ξ s is continuous on Γ + .with θ 1 ∈ [0, 2π] and θ 2 , . . . , θ d-1 ∈ [0, π]. Notice that the assumption ω ∈ Γ -(x) actually implies that θ d-1 ∈ [0, π/2). We will write θ = (θ 1 , . . . , θ d-2 , θ d-1 ) and U := [0, 2π]×[0, π] d-3 × 0, π 2 , so that ω ∈ Γ -(x) =⇒ θ ∈ U . Notice that the set U is independent of x ∈ ∂Ω. Within this frame and with the above set of coordinates, we writeT (x, θ 1 , . . . , θ d-2 , θ d-1 ) := τ + (x, ω)andG(x, θ 1 , . . . , θ d-2 , θ d-1 ) = x + T (x, θ 1 , . . . , θ d-2 , θ d-1 )ω, ∀ω ∈ Γ -(x).From Theorem A.2, for any x ∈ ∂Ω, the mappingθ ∈ U → T (x, θ) is of class C 1 . It is then clear that the set of ω ∈ S d-1given by (A.14) such that the vectors∂ θ j G(x, θ) j=1,...,d-1are linearly independent coincide with the set at which the differential of the mappingω ∈ Γ -(x) → x + τ + (x, ω)ω ∈ ∂Ω is of full rank d -1. One has ∂ θ j G(x, θ) = ∂ θ j T (x, θ)ω + T (x, θ)∂ θ j ω and, because |G(x, θ) -x| = T (x, θ), it is easy to check that ∂ θ j T (x, θ) = ω • ∂ θ j G(x, θ), i.e. ∂ θ j G(x, θ) -(∂ θ j G(x, θ) • ω)ω = T (x, θ)∂ θ j ω. Notice that ∂ θ j G(x, θ) -(∂ θ j G(x, θ) • ω)ω is nothing but the projection of ∂ θ j G(x, θ) on the hyperplane ω ⊥ . Recalling that T (x, θ) > 0, we see that ∂ θ j G(x, θ) -(∂ θ j G(x, θ) • ω)ω j=1,...,d-1are independent if and only if (∂ θ j ω) j=1,...,d-1 are independent, or equivalently, if the Gram matrixJ θ (ω) = ∂ θ i ω , ∂ θ j ω i,j is not singular. It is well known that det (J θ (ω)) = sin d-2 θ d-1 sin d-3 θ d-2 . . . sin θ 2 . (A.15) In other words, if det (J θ (ω)) = 0, then ∂ θ j G(x, θ) -(∂ θ j G(x, θ) • ω)ω j=1,...,d-1 are independent and one deduces easily that then ∂ θ j G(x, θ) j=1,...,d-1 are also independent. We define then S(x) = {ω ∈ Γ -(x) ; det (J θ (ω)) = 0} the mapping ω ∈ Γ -\ S(x) → x + τ + (x, ω)ω has a differential of rank d -1. It is clear that S(x) is closed. Let us now prove that indeed S(x) has a zero surface Lebesgue measure. Using then (A.15), we get that ω ∈ S(x) if and only if θ j ∈ {0, π} for some j = 2, . . . , d -1.

S d- 1 1

 1 Sε(x) (ω)σ(dω) = Uε det (J θ (ω)) dθ ε π d-2whereU ε = {θ ∈ U ; sin d-2 θ d-1 sin d-3 θ d-2 . . . sin θ 2 ε}. Therefore lim ε→0 + sup x∈∂Ω σ (S ε (x)) = 0.

2 .

 2 7, U : Γ + → Γ + is a C 1 diffeomorphism from Γ + onto its image U( Γ + ) which is an open set of Γ + . Let us introduce ∆ := Γ + \ Γ + = {(x, v) ∈ Γ + ∪ Γ 0 ; ξ(x, v) ∈ Γ 0 } As already noticed in (A.1), µ(∆) = 0, i.e. ∆ is a closed set of Γ + of zero dµ-measure. Since U (Γ + ∪ Γ 0 ) = U(Γ 0 ) ∪ U( Γ + ) ∪ U(∆) † namely Cj = {ω ∈ Γ+(x) ; θ d-j = 0 or θ d-j = π} B.Sweeping property. We define now the sweeping property for stochastic semigroups: Definition B.3. A stochastic semigroup {P (t)} t 0 on the space L 1 (E, Σ, m) is sweeping from a set A if lim t→∞ A P (t)f (x) m(dx) = 0

  according to [32, Lemma 5.2.11 & Theorem 5.2.16, pp. 128-131], the Lebesgue surface measure π Ey (dY ) on E y is given by

  23, Theorem 2.1]), this implies that 1 is an isolated eigenvalue of M 0 H. Moreover, being M 0 H irreducible, we deduce from [23, Theorem 2.2] the uniqueness and the strict positivity (almost everywhere) of a nonnegative eigenfunction ϕ. Let us consider now λ ∈ C with Reλ > 0. Considering the modulus operator |M λ H| (see [14] for a precise definition) one has |M λ H| M 0 H and |M λ H| = M 0 H. We conclude that A = T H thanks to [4, Theorem 4.5]. Let us now show that the eigenfunction ϕ lies

In particular, from

[START_REF] Marek | Frobenius theory of positive operators: Comparison theorems and applications[END_REF] Theorem 4

.3], r σ (|M λ H|) < r σ (M 0 H) = 1. Since moreover, r σ (M λ H) r σ (|M λ H|) according to [14, Theorem 1], this proves that r σ (M λ H) < 1, i.e. 1 /

∈ S(M λ H).

  According to Theorem B.4 and the previous Lemma, we have: Theorem 8.3. Let us assume that (U H (t)) t 0 is stochastic and has no invariant density. If the boundary operator H is strongly local diffuse then (U H (t)) t 0 is sweeping from all compact subsets of E 0 .

This paper was partially supported by the Polish National Science Centre Grant No. 2017/27/B/ST1/00100 (RR). This work was written while R.R. was a visitor to Université de Franche-Comté and the authors thank Université de Franche-Comté for financial support of this visit.

Letting ε → 0, the monotone convergence theorem yields to U H (t)Ψ Ψ, i.e. Ψ satisfies (8.6) and the proof is concluded.

APPENDIX A. ABOUT THE BALLISTIC FLOW

We establish in this appendix several important properties of the so-called ballistic flow

which are fundamental for the proof of our main weak compactness result Theorem 5.1. For the clarity of exposition, we postponed these results in an Appendix but strongly believe that the results stated here have their own mathematical interest. In this Appendix, we will use the following notations: for any element z = (x, v) of the extended phase space Ω × V , we will call x the space component of z and v the velocity component of z, writing x = z s and v = z v .

With the notations of [START_REF] Guo | Decay and continuity of the Boltzmann equation in bounded domains[END_REF], ξ = x b . Notice that, as already observed in [START_REF] Guo | Decay and continuity of the Boltzmann equation in bounded domains[END_REF][START_REF] Voigt | Functional analytic treatment of the initial boundary value problem for collisionless gases[END_REF], in non convex domain this deterministic flow does not avoid the grazing set Γ 0 , i.e. in full generality ξ(x, v) ∈ Γ -∪ Γ 0 and -as far as the regularity of ξ is concerned -the set {(x, v) ∈ Γ + ; ξ(x, v) ∈ Γ 0 } will be particularly relevant. Notice though that

Moreover, according to (2.5) with ψ = 1 Γ 0 we see that

We introduce the following where we focus on velocity which are unit vectors (this is no loss of generality by virtue of (2.2))

and introduce, for any x ∈ ∂Ω the section

A.1. Regularity of the travel time. The main result of this section is the following:

Theorem A.2. The set Γ ± are open subsets of Γ ± and

We will split the proof of the above in a series of Lemma -dealing with τ -but all the results have their counterpart for τ + :

Proof. Let us fix x 0 ∈ ∂Ω and ω 0 ∈ Γ + (x 0 ), i.e. ω 0 • n(x 0 ) > 0 and ω 0 • n(ξ 0 ) < 0 where ξ 0 = x 0 -τ -(x 0 , ω 0 )ω 0 is the space component of ξ(x 0 , ω 0 ). For simplicity, set

Let (x n , v n ) ⊂ Γ + be a given sequence such that lim n (x n , v n ) = (x 0 , ω 0 ). In particular, we can assume that |v n | = 0 for any n ∈ N. Set then ω n = |v n | -1 v n ∈ S d-1 for any n ∈ N and

Taking a subsequence if necessary (recall that ∂Ω is compact), we may assume that (y n ) n converges to some y 0 ∈ ∂Ω. Then, since τ

and, consequently, letting n goes to infinity in the definition of y n yields

This in particular shows that τ 0 τ 1 . To prove that τ 1 = τ 0 , let us argue by contradiction and assume that τ 1 > τ 0 . Since both y 0 = x 0 -τ 1 ω 0 and ξ 0 = x 0 -τ 0 ω 0 belong to ∂Ω and since ω 0 • n(ξ 0 ) < 0, the set

∈ Ω is open and not empty. Therefore, there exists δ > 0 such that τ 1 > δ + τ 0 and

we get that, for n ∈ N large enough,

Letting then n goes to infinity, we obtain x 0 -tω 0 ∈ Ω for any t ∈ (0, τ 0 + δ) which contradicts (A.2). Therefore, τ 1 = τ 0 which proves the continuity of τ -on Γ + . Let us now show that Γ + is open. We keep the previous notations, fixing (x 0 , ω 0 ) ∈ Γ + . Let us assume that there exists a sequence

which proves the Lemma under assumption (A.9). It only remains to prove (A.9). Notice that

is the (d -1)-dimensional tangent space of ∂Ω at z 0 = Ψ(0) with x 0 -Ψ(0) = τ -(x 0 , ω 0 )ω 0 and

is the orthogonal projection of Range(dΨ(0)) on the orthogonal hyperplane to ω 0 . One sees that

and consequently P ω 0 (Range(dΨ(0))) coincides with the orthogonal hyperplane to ω 0 . In particular, it has dimension d -1 which is exactly (A.9).

Lemma A.5. For any ω ∈ S d-1 , the mapping

is differentiable and

Proof. As in the previous Lemma, we fix x 0 ∈ ∂Ω and ω 0 ∈ Γ + (x 0 ) and consider an open neighbourhood U(x 0 ) × V(ω 0 ) of (x 0 , ω 0 ) on which (A.4) and (A.5) hold. For any ω ∈ V(ω 0 ), define

and, with the notation of the previous Lemma, S ω ⊂ W(z 0 ) for any ω ∈ V(ω 0 ) where W(z 0 ) is the image of the C 1 function

with Ψ(0) = z 0 and Rank(dΨ(y)) = d -1 for any y ∈ (-1, 1) d-1 . This allows to introduce, as in the previous Lemma, O ω = Ψ -1 (S ω ) and S ω is parametrized by Ψ (defined now on O ω ), i.e. given (x, ω) ∈ U(x 0 ) × V(ω 0 ), there is a unique y ∈ O ω such that x -τ -(x, ω)ω = Ψ(y) and (A.6) and (A.7) still hold. We have seen in the proof of Lemma A.4 that H is C 1 with differential d y H(x, y) given by (A.8) and depending continuously on (x, y) ∈ U(x 0 ) × (-1, 1) d-1 . In particular, as seen earlier, at (x, y) = (x 0 , 0) the differential d y H is given by

and has (d -1)-dimensional range. As before, from the implicit function theorem, there is a neighbourhood (-δ, δ) d-1 ×U (x 0 ) of (0, x 0 ) and a neighbourhood V (ω 0 ) of ω 0 on which the equation

Remark A.6. Notice that (A.13) allows to recover the expression

which was obtained in [20, Lemma 3] for some special structure of Ω. Moreover, using (A.10) and using the range of the differential of Ψ(G(x, ω)) is the orthogonal of Span(n(G(x, ω)) we can prove

which, again, is a result obtained in a special case in [20, Lemma 3].

Proof of Theorem A.2. The above three Lemmas give directly the proof of Theorem A.2 for τ -and Γ + . The proof for τ + and Γ -is done similarly.

An immediate consequence of Theorem A.2 is the following regularity of the ballistic flow:

Corollary A.7. The ballistic flow:

A.2. Further non degeneracy results. We introduce a local polar parametrization of the boundary which will turn useful later on:

Proposition A.8. For any x ∈ ∂Ω, there is a closed subset S(x) ⊂ Γ -(x) with zero surface Lebesgue measure dσ and such that the mapping ω ∈ Γ -(x) \ S(x) → x + τ + (x, ω)ω ∈ ∂Ω has a differential of rank d -1. As a consequence, the differential of the mapping

Proof. For any x ∈ ∂Ω, we choose an orthonormal basis {e 1 (x), . . . , e d-1 (x), e d (x)} -depending continuously on x ∈ ∂Ω -where e d (x) = -n(x).

Let us write the components of ω ∈ Γ -(x) in this basis using polar coordinates

we have that µ U( Γ + ) = µ (U (Γ + ∪ Γ 0 )) = Γ + ∪ Γ 0 since U is µ-preserving as a composition of the µ-preserving mapping ξ and V. Therefore

Introduce then Γ

(1)

+ . One has

+ onto its image. Since U 2 is µ-preserving, writing the disjoint unions Γ

(1)

By induction, assuming that there is Γ

) which is of full µ-measure, i.e.

. As before, writing the disjoint unions

On the other hand, according to Corollary A.7,

Using the definition (A.16) for k + 1 we see that

. Arguing as before and writing

and since ξ • U k is µ-preserving, we have that

where we used that ξ is also µ-preserving. Therefore

) is an open subset of Γ -then setting

one sees that γ k is a closed subset of Γ -with µ(γ k ) = 0 and

In the next Lemma, we used the notations of Proposition A.8:

Lemma A.11. For any x ∈ ∂Ω, let S(x) ⊂ Γ -(x) be the closed set introduced in Proposition A.8. For any k ∈ N, let

where ∆ k+1 is defined thanks to (A.16). Then, for any ω ∈ O k (x) \ S(x), the differential

Proof. Notice that open subsets of ∂Ω × S d-1 are C 1 manifolds of dimension 2(d -1) and then, so are open subsets of Γ + or Γ -. According to Proposition A.8, the differential

), the differential dU -k (y, ω ) is an isomorphism between the (2d -1) dimensional tangent spaces. Finally, the differential of the mapping

APPENDIX B. REMINDERS ON PARTIALLY INTEGRAL STOCHASTIC SEMIGROUPS

We collect here several results on partially integral stochastic semigroups in L 1 (E, Σ, m) where (E, Σ, m) is a given σ-finite measure space. for some t > 0.

We have then the following (see [START_REF] Ór | Continuous Markov semigroups and stability of transport equations[END_REF])

Theorem B.2. Let {P (t)} t 0 be a partially integral stochastic semigroup. Assume that the semigroup {P (t)} t 0 has a unique invariant probability density f * . If f * > 0 a.e., then the semigroup {P (t)} t 0 is asymptotically stable.

Let P(t, y, B) be a probability transition function for the semigroup {P (t)} t 0 , i.e.