On Neumann problems for nonlocal Hamilton-Jacobi equations with dominating gradient terms

We are concerned with the well-posedness of Neumann boundary value problems for nonlocal Hamilton-Jacobi equations related to jump processes in general smooth domains. We consider a nonlocal diffusive term of censored type of order less than 1 and Hamiltonians both in coercive form and in noncoercive Bellman form, whose growth in the gradient make them the leading term in the equation. We prove a comparison principle for bounded sub-and supersolutions in the context of viscosity solutions with generalized boundary conditions, and consequently by Perron's method we get the existence and uniqueness of continuous solutions. We give some applications in the evolutive setting, proving the large time behaviour of the associated evolutive problem under suitable assumptions on the data.


Introduction
The aim of this work is the analysis of the well-posedness of Neumann boundary value problems for partial-integro differential equations (PIDEs in short) of Hamilton-Jacobi type, where the nonlocal terms are singular integrals related to the infinitesimal generator of discontinuous jump processes. To be more specific, we consider the following is an integro-differential operator of censored type and of order strictly less than 1 (see (1.2) for the definition).
In the probabilistic approach to PDEs, Neumann boundary conditions are associated to stochastic processes being reflected on the boundary. The underlying idea is to force the stochastic process to remain inside the domain of the equation. Classically, this is obtained essentially by a reflection on the boundary (see the method developed by Lions and Sznitman [25] in the continuous setting). A key result in the classical setting is that, for a PDE with Neumann boundary conditions, there is a unique underlying reflection process and any consistent approximation will converge to it (see [25] and Barles, Lions [10]).
When dealing with discontinuous jump processes, the underlying idea is the same but the situation is different. This is essentially due to the fact that the jump process may exit the domain without having first hit the boundary. The This work was partially supported by the ERC advanced grant 668998 (OCLOC) under the EU's H2020 research programme.
consequence is that Neumann boundary conditions can be obtained in many ways, depending on the kind of reflection we impose on the outside jumps. Moreover, the choice of a reflection on the boundary changes the equation inside the domain.
The starting point of our work is the paper [6] where Barles, Chasseigne, Georgeline and Jakobsen studied problems as in (1.1) in the case of linear equations (that is, without the Hamiltonian H) and when the domain Ω is the half space. In [6] different models of reflection are presented, among which two types of reflections are particularly relevant for possible extensions in a more general setting. The first is the normal projection, close to the approach of Lions-Sznitman in [25], where outside jumps are immediately projected to the boundary by killing their normal component. This model has been thoroughly investigated in the paper [7] for fully non-linear equations set in general domains.
The second, the censored model, is the one we consider in our paper. In this case, any outside jump of the underlying process is cancelled (censored) and the process is restarted (resurrected) at the origin of that jump.
In particular, in the present work, we consider the boundary value problem (1.1) where I[u] is an integro-differential operator of censored type and of order strictly less than 1 defined as where µ x is a singular non-negative Radon measure representing the intensity of the jumps from x to x + z and satisfying the following integrability condition |z| ∧ 1dµ x (z) < +∞, and j(x, z) is a jump function (see assumptions (M), (J0), (J1) in the following section for details). A meaningful example is where g, f are bounded and Lipschitz functions. Note that I has to be interpreted as a principal value (P.V.) integral. We remark that the domain of integration is restricted to the z such that x + j(x, z) ∈Ω, avoiding thus any outside jump. Note also that, as a consequence of the fact that censored type processes are not allowed to jump outsideΩ, we don't need any conditions on Ω c in the boundary value problem (1.1). We follow the PIDE analytical approach developed in [6], in the sense that we directly work with the infinitesimal generator and not yet with the processes themselves. For more details and probabilistic references on censored processes, we refer to e.g. [14], [21], [24], [22] and to the introduction of [6]. We just mention that the underlying processes in this paper are related to the censored stable processes of Bogdan [14] and the reflected σ-stable process of Guan and Ma [22].
We stress that the boundary value problem (1.1) is interpreted in the sense of viscosity solutions, meaning that the Neumann boundary condition could also be not attained (and in this case the equation holds up to the boundary).
In the case of linear PIDES, as considered in [6], the kind of singularity of µ influences the nature of the boundary value problem (1.1), in the sense that the Neumann boundary condition is attained only if the measure is singular enough.
In particular in [6] it is shown that, when the singularity is of order strictly less than 1 as in (1.3), the equation holds up to the boundary and the process never reaches the boundary. On the other hand, when the singularity of the measure is strong, i.e. when µ is of the type (1.3) with σ ∈ [1, 2), the situation is far more complicated, mainly due to the "ugly" dependence in x of the operator in (1.2) and to the interplay between the singularity of the measure and the geometry of the boundary. In [6] this difficulty is tackled by considering solutions which are in some sense Hölder continuous up to the boundary and the comparison principle is established only in this class. Though the result could not be optimal, it is consistent with the "natural" Neumann boundary condition for the reflected σstable process (proved by Guann and Ma [22] through the variational formulation and Green type formulas) which in the case of the half space reads This allows the normal derivative to growth less than |x N | σ−2 and then suggests that it is appropriate to look for solutions which are β-Hölder continuous, with β > σ − 1, as assumed in [6]. We remark that the previous argument suggests also that, on the contrary, in the case σ < 1 there is no need to assume any further regularity.
The situation is different when dealing with nonlinear equations as (1.1), which we consider in this paper. Indeed, the presence of the Hamiltonian term H in (1.1) entails further difficulties even in the case of measures of order strictly less than 1 (e.g. as in (1.3)). This is due to the fact that the nonlinear term H could force the process to hit the boundary and, consequently, the Neumann boundary condition to be attained.
In order to deal with this difficulty, we consider a class of Hamiltonians with a gradient growth stronger than the diffusive term in the nonlocal operator (1.2). The first example are Hamiltonians H with superfractional coercive growth in the gradient variable, namely (1.5) H(x, p) = a(x)|p| m − f (x), where m > σ, σ ∈ (0, 1), a, f :Ω → R are bounded and continuous functions and a(x) ≥ a 0 > 0 for some fixed constant a 0 . We remark that the positivity of a and the condition m > σ make the first-order term the leading term in the equation. We also observe that we have no other additional restriction to m (in particular, we can deal with Hamiltonians as in (1.5) with m < 1), allowing the study of Hamiltonians which are concave in Du.
The second main example are Hamiltonians H of Bellman type, which arises in the study of Hamilton-Jacobi equations associated to optimal exit time problems, such as (1.6) H(x, p) = sup where A is a compact metric space (the control space) and b, l are continuous and bounded functions (we refer the reader to [2] and [20] for some connections between this type of equations and control problems). Note that the diffusive term of I defined in (1.2) is of weaker order than the first-order term when we assume σ < 1. We also observe that, as in [13] and [28], the well-posedness of (1.1) with Hamiltonians as in (1.6) is based on a careful study of the effects of the drift b at each point of ∂Ω × (0, +∞). Note that the drift b is not allowed to be parallel to the boundary for all controls, differently from the local case. This assumption on b makes the Hamiltonian the leading order term in the equation and thus allows us to control the growth of the nonlocal term.
The main result of our paper is the comparison principle between bounded sub and super-viscosity solutions to (1.1), see Theorem 3.1. We remark that the proof of this result is not standard even in the case σ < 1 in the half space. The difficulties are mainly due to the fact that operators as in (1.2) behave badly in x. The main idea which is behind the proof is to localize the argument on points which have the same distance from the boundary and this is carried out through the use of a non-standard non regular test function. The main assumption which allows us to localize on equidistant points is the superfractional growth of the Hamiltonian term, see in particular the proof of Proposition 5.1 and Proposition 5.6 (more precisely, Lemma 5.2 and Lemma 5.7) for Bellman and coercive Hamiltonians respectively.
After the localization procedure, the rest of the proof in the case of the half space is simple, whereas in the case of general domains, further technical difficulties arise form the way the x-depending set of integration of I interferes with the geometry of the boundary. To face these extra technical difficulties, we rectify the boundary relying on the smoothness of Ω. This is done in Lemma 4.1 which is a key result used in the proof of Theorem 3.1, which we prove before Theorem 3.1 in Section 4.
The first main application of our result is the proof of existence and uniqueness for (1.1), by standard Perron's method (Corollary 3.2).
Then, in Section 6, we present some applications of our results to the evolutive setting. In particular, we prove the well-posedness of the Cauchy problem associated to (1.1) and we study two different kind of asymptotic behaviour under suitable assumptions on the data. We refer to Section 6 for precise assumptions, statement of the results and proofs.
Finally, we remark that in [13] Dirichlet boundary value problems associated to nonlocal Hamilton-Jacobi parabolic equations are studied. The main interest of [13] is to prove the well-posedness in the context of loss of the boundary condition, that is, in the case of Hamiltonians with dominating gradient terms, as the ones considered in the present paper. However, the nonlocal operators treated in [13] are very different from the censored operators (1.2) that we study here, since they are defined in all the space (and can therefore be seen as generalizations of the fractional laplacian). Moreover, differently from (1.1), in [13] also boundary conditions outside the domain are requested. This is strongly connected to the fact that here we consider Neumann boundary value problems with censored operators, that is, roughly speaking, we do not allow the process to jump outside the domain and therefore, no conditions outside is needed.
1.1. Organization of the paper. In Section 2 we state the assumptions on the nonlocal operator and the Hamiltonian and we give the definition of solution to problem (1.1). In Section 3 we state the main results, that is the uniqueness and existence for problem (1.1) for a Hamiltonian either coercive or of Bellman type (Theorem 3.1 and Corollary 3.2). In Section 4 we prove Lemma 4.1 and in Section 5 we prove Theorem 3.1. In Section 6 we treat the associated evolutive problem, studying uniqueness, existence and asymptotic behaviour of the associated evolutive problem for large time. Finally, in the Appendix we prove some lemmas used in the proof of Theorem 3.1.

Assumptions and definition of solutions
This means that for anyŝ ∈ ∂Ω there exists r = r(ŝ) and a W 2,∞ -diffeomorphism ψ : B r (ŝ) → R N satisfying ψ N (s) = d(s) for any s ∈ B r (ŝ), where d is the signed distance from the boundary of Ω.

Remark 1.
By assumption (O), there exists a neighbourhood of the boundary of Ω where the distance from the boundary d is smooth. Unless otherwise specified, throughout the paper we denote by d a function which coincides with the signed distance from the boundary of ∂Ω in this neighbourhood and is bounded in all the domain. We recall that by n(x) we denote the exterior unit normal vector to ∂Ω and we write n(x) = −Dd(x) in the neighbourhood of the boundary where d is smooth.
We consider non-negative Radon measures with density dµx dz satisfying (M) there exists C µ , D µ > 0, σ ∈ (0, 1) such that for any x, y ∈Ω, For example, (M) is satisfied by where σ ∈ (0, 1), g : R N × R N → R is a non-negative bounded function such that g(·, z) is Lipschitz uniformly with respect to z.
Concerning the jump function j we assume (J0) for any x ∈Ω, j(x, ·) ∈ C 1 (R N ), j(x, ·) is invertible, j −1 (x, ·) ∈ C 1 (R N ) and there exists a constant A j such that (J1) there existC j , C j , D j > 0 such that for any x, y ∈Ω, z ∈ R N , it holds For example (J0), (J1) are satisfied for and satisfies the assumptions below. We assume also: (C) Uniform continuity of the cost l: There exists a modulus of continuity ω l such that |l(x, α) − l(y, α)| ≤ ω l (|x − y|) ∀α ∈ A, ∀x, y ∈Ω; (L) Uniform Lipschitz continuity of the drift b: We introduce the following notation Roughly speaking, Γ in and Γ out can be respectively understood as the set of points where the drift term pushes inside and outside Ω the trajectories. In order to avoid two completely different drifts behaviour for arbitrarily closed points, we assume that each of these subsets is uniformly away from the others, as encoded in the following assumption (B). For example, if ∂Ω is connected, then it consists in one piece belonging to one of Γ in , Γ out and Γ; otherwise, we are able to deal with boundary with several components of different types, precisely each one belonging to one between Γ in , Γ out , Γ.
The assumptions we do on these subsets are the following Remark 2. Note that the strict sign in the definition of Γ in , Γ out and Γ is fundamental, since it makes the Hamiltonian the leading order term in the equation, allowing us to control the growth of the nonlocal term, which is of order strictly less than 1.

Remark 3.
In order to treat the points of Γ in , we use the existence of a blowup supersolution exploding on the boundary. We follow the same approach of [6], where the existence of a blow-up supersolution is proved for censored type operators (of order strictly less than 1) when the measure of integration satisfies specific assumptions (in particular does not depend on x and there exists at least one point where it is strictly positive). In this particular case it is shown in [6] that the integral term computed on the blow-up supersolution does not explode at the boundary. This is not true anymore when considering more general measures as we consider in (M). In order to solve this difficulty, we assume the strict sign in the behaviour of the drift term on Γ in , which allows us to control the growth on the boundary of the integral term computed on this blow-up supersolution. We refer to the proof of Proposition 5.1 and in particular to Lemma 5.5 for further details.

Coercive Hamiltonian and Examples.
We consider superfractional coercive Hamiltonians: (H1) Let σ be as in (M). There exists m > σ, c 0 > 0, D > 0 such that for all We distinguish the case of sub or superlinear coercivity: Sublinear coercivity: We say that H is sublinearly coercive if it satisfies (H1) for m ≤ 1 and the following continuity condition holds: (Ha) There exists a constant C > 0 and modulus of continuity ω 1 such that, for all x, y, q, p ∈ R N , we have Superlinear coercivity: We say that H is superlinearly coercive if: (Hb) There exists m > 1,A,C > 0 such that for all µ ∈ (0, 1), x, y, p ∈ R N , we have (Hc) If m is as in assumption (Hb), there exist C > 0 and a modulus of continuity ω 1 such that, for all x, y, q, p ∈ R N Remark 4. Note that condition (Hb) implies (H1) for m > 1.
As it is classical in viscosity solution's theory, the comparison principle allows the application of Perron's method to conclude the existence of solutions. To this end, we introduce the following assumption, which will allows us to build sub and supersolutions: (E) There exists H R > 0 such that for any p ∈ R N , |p| ≤ R ||H(·, p)|| ∞ ≤ H R .
As a model example for sublinearly coercive Hamiltonians, we consider with m ≤ 1, a 1 ≥ a 0 > 0 for all x ∈ Ω, l < m and a 1 , a 2 , f :Ω → R are continuous and bounded functions and a 1 , a 2 are also Lipschitz continuous.
As a model example for superlinearly coercive Hamiltonian, we consider with m > 1, l ≤ m, b bounded and continuous and a 1 , a 2 , f as before. These Hamiltonians are coercive in p and in the case m > 1 we can include transport terms with a Lipschitz continuous vector field b :Ω → R N . The above assumptions are easily checkable in both cases.

Notion of viscosity solutions.
We recall now the definition of solution to problem (1.1). We use the following notation: The I ξ -term is well-defined for any bounded function φ. The I ξ -term is well-defined for φ ∈ C 1 thanks to (M0).
We also denote Following the approach of [6], we give the definition of viscosity solution to (1.1). Let C j be defined as in (J1).
(i) A bounded usc function u is a viscosity subsolution to (1.1) if, for any test-function φ ∈ C 1 (R N ) and maximum point (ii) A bounded lsc function v is a viscosity supersolution to (1.1) if, for any test-function φ ∈ C 1 (R N ) and minimum point (iii) A viscosity solution is both a sub-and a supersolution.

Main results
The main result of this part is the following comparison principle for the problem (1.1). Let u be a bounded usc subsolution of (1.1) and v a bounded lsc supersolution of (1.1). Then u ≤ v inΩ.
Once the comparison holds, we use the Perron's method for integro-differential equations (see [1], [8], [26] and [19], [23] for an introduction on the method) to get as a corollary existence and uniqueness for the problem (1.1) either when H is of Bellman type either for H coercive.

A preliminary key lemma
We prove the following Lemma 4.1, which is a key result used in the proof of Theorem 3.1. Roughly speaking, it deals with the difficulties arising from the way the geometry of the boundary interferes with the singularity of the nonlocal terms. The scope is to estimate the nonlocal terms defined in (4.1) on points near the boundary and equidistant from it. The approach of the proof is essentially based on a rectification of the boundary, relying on its regularity.

Remark 5.
In the case of domains with flat boundary, we do not need Lemma 4.1 in the proof of Theorem 3.1 since the estimation of the nonlocal terms can be carried out more easily. We refer to Remark 7, step 4 of the proof of Theorem 3.1.
We observe that by assumption (O) and (4.2), we have We proceed performing a change of variable in order to write the set of integration in terms of ψ N (s 1 ). In other words, we write Moreover, by (J1), (O) and since ψ ∈ W 2,∞ , we get (4.10) Similarly one can prove an analogous upper bound for |w|, getting (4.11)Ĉ|z| ≤ |w| ≤C|z|, for some constantsĈ,C depending on the Lipschitz constants of ψ, ψ −1 and oñ C j , C j of (J1). In the following step, we rewrite the set of integration of I[J s1 /J s2 ] in a different way using the change of variable (4.8) and the previous estimates.
Step 2. -Rewriting the set of integration. By (4.8) and if ψ N (s 2 + j(s 2 , z)) ≤ 0, we have For convenience of notation, let for the moment . Then, since ψ ∈ W 2,∞ and by (4.12) we write From now on we denote by C any positive constant which may change from line to line. By definition of ζ(t) Then, by (4.14), (J1) and (4.10) we get Now we analyse A 2 . Note that by (4.13) and (4.7) Moreover, since ψ ∈ W 2,∞ , by (J1) and (4.10) Then we have A 2 ≤ C|w||s 1 − s 2 |. We denote a = ψ N (s 1 ) and observe a ≥ 0. By (4.9) and the previous arguments, we get Then, we perform the change of variable in I[J s1 /J s2 ] by (J0), (J1) and using that ψ ∈ W 2,∞ and by (4.15) and (4.11), we get andC,C, C are positive constants, which we supposeC =C = C = 1, by no loss of generality and for simplicity of exposition.
Step 3. -Estimate onD. We introduce the following notation: Note that by the second assumption in (4.2), |s 1 −s 2 | ≤ r. Without loss of generality we can suppose r ≤ 1 2 , so that we have |s 1 − s 2 | ≤ 1/2. Then Note that, if w ∈D, then We identify two cases, depending on the sign of −a + |s 1 − s 2 ||w ′ | and we denote By all the previous observations, we write For F 1 , we use that 1 |w ′ | 2 +|wN | 2 ≤ 1 |w ′ | 2 and by Fubini's Theorem, we integrate in the N -variable and we get By the first of (4.17) and (4.18) and since da ≥ 0, we have −da+ From now on we denote by C any positive constant which may change from line to line. Note that, since w ′ ∈ R N −1 and σ < 1, we have Then by the previous observations, we get (4.24) Now we analyse F 2 . For simplicity of notation, we denote and then, by Fubini's Theorem, we have We split the domain as follows We estimate the first term by (4.27) where in the first inequality we used that −βa + β|s 1 − s 2 ||w ′ | + a ≤ 2|w ′ ||s 1 − s 2 |, since β ≤ 1 and a ≤ |w ′ |, and in the second inequality we used (4.23). Take now the second term in (4.26). Note that, if a > |w ′ |, by (4.17) and is increasing on the negative half line, we have and coupling (4.28) and (4.29), we get Then coupling (4.25), (4.26), (4.27) and (4.30), we obtain and we conclude the proof by coupling (4.16), (4.20), (4.24) and (4.31).

Proof of the comparison principle
We prove Theorem 3.1 and we split the proof into two parts, depending whether H is of Bellman type or coercive.

5.1.
Hamiltonians of Bellman type. The proof of Theorem 3.1 follows mainly by the following proposition, which we prove first. At the end of the proof of Proposition 5.1, we will prove Theorem 3.1.
satisfies, in the viscosity sense, the equation where B is a positive constant depending on the data.
Proof. Let x 0 ∈Ω and φ ∈ C 1 (R N ) such that ω − φ has a maximum point at x 0 inB Cj ξ (x 0 ) ∩Ω for some 0 < ξ < 1. Note that the restriction ξ < 1 is only for simplicity of exposition. We can suppose that x 0 is a strict maximum point by consideringφ = φ(x) + k|x − x 0 | 2 for k sufficiently large. In particular, note that we have the following property From now on, we consider the test functionφ and we still denote it by φ, with some abuse of notation. We observe that if x 0 ∈ Ω the proof is rather standard, since in this case the maximum points (x, y) of u − v − φ converge as ε → 0 to (x 0 , x 0 ) and hence they are bounded away from the boundary for ε small enough. This last property implies that we can directly use the equations and then proceed as in the following case. Let Γ in , Γ out , Γ be defined respectively in (2.3), (2.4) and (2.5) and recall they satisfy (B). We suppose x 0 ∈ ∂Ω and we split the proof depending if In case (a) we use the existence of the blow-up supersolution which explodes at the boundary and allows us to keep the maximum points far from the boundary. Since the proof in this case is inspired by a similar approach used in [6], we give the details at the end of the proof of (b) and (c) in Remark 9. Now we treat case (b) and (c). Since the proofs are similar, we treat them at the same time.
We suppose that and prove that the F-viscosity inequality of Definition (2.1) hold for ω.
Step 1. Localising on equidistant points (that is, d(x) = d(y)). Let ε > 0 and d be a function as in Remark 1. We double the variable by introducing the function for ε, δ > 0 where χ ε : R → R (and similarly χ δ ) is defined as follows and K > 0 is a constant large enough such that where γ, C 2 > 0 depend on x 0 and are precisely defined in Lemma 6.10 in the Appendix (forŝ = x 0 ). Let We observe that (x,ȳ) depends now also on δ and we omit the dependence. Now consider where d is the signed distance from the boundary (see Remark 1), χ ε is defined as in (5.4) and K is as in (5.5). Note that the test function in (5.8) is not differentiable on the points such that d(x) = d(y). By upper-continuity, Ψ in (5.7) attains its maximum over . By classical arguments in viscosity solution theory, we get as ε → 0 . Note also that by the second of (5.9), it directly follows We prove the following key lemma.
Lemma 5.2. Under the above notation, we have Proof. Note that (i) follows by classical argument in viscosity solution theory. We remark that the proof of (ii) is slightly different in case (b) and case (c). We argue by contradiction and we suppose that d(x) = d(y). First we prove that the Fviscosity inequalities for u and v of Definition (2.1) hold. Suppose that x ∈ ∂Ω, then d(x) = 0 and d(y) = 0. We denote Note that by (5.9), we can suppose that x, y are close to the boundary, by taking ε small enough. By the Taylor's formula for the distance function, we have for ε small enough and then By (5.11), (5.13), (5.14) and (5.9), we have . Note that, from (5.2), for ε small enough we have also By (5.12), (5.16), (5.15) and since K ≥ 0, we conclude for ε small enough Then, since u is a viscosity subsolution and the function u(·) − v(y) − ψ(·, y) has a local maximum at x, the F -viscosity inequality of Definition (2.1)(i) holds. A similar argument can be carried out for v. From now on, we treat separately Case Then, by Definition (2.1)(i) and by (5.17), we have We apply Lemma 6.10, (6.17) withŝ = x 0 , p = ε −1 χ ′ ε (|x − y|)p + q and λ = ε −1 K and by the definition (5.5) of K, we get for ε small where by C, here and in the following, we denote any positive constant independent of ε which may change from line to line. To estimate the nonlocal terms we use the following lemma, which we prove in the Appendix.
Lemma 5.3. Let I ξ , I ξ be as in (2.7), (2.6) and assume the first of (M ) and (J1). Under the above notation, for any ξ > 0, there exist a positive constants C 1 independents of ε such that Then, by (5.20), by Lemma 5.3 (i) with ξ ′ = ε and by the boundedness of u, we write (5.18) as follows and we reach a contradiction for ε small enough, since C is independent of ε and σ < 1. Now suppose d(x) < d(y). In this case we use the following F -viscosity inequality wherep is defined in (5.11) and q is defined in (5.19). Then , for ε small enough, we apply Lemma 6.10, (6.18) We proceed as in the previous case, we apply Lemma 5.3 (ii) with ξ = ε and by (5.22) and the boundedness of v , we get and we reach a contradiction for ε small enough as above.
Step 2. Writing the viscosity inequalities. By (5.9) and Lemma 5.2 (i), from now on we consider δ, ε small enough so that Now we prove that the F -viscosity inequalities for u and v hold. We takex ∈ ∂Ω and we show that the boundary conditions do not hold, so the F -viscosity inequalities hold as in Definition (2.1). We proceed exactly as in Step 1, Lemma 5.2, so we omit the details. We recall that for all δ > 0 for all x, y ∈Ω, and we note only that since d(x) = 0, we have for ε, δ small enough Sinceφ ∈ C 1 , by (J1) and the first of (M), we have where C j is as in (J1) and o ξ ′ (1) is independent of δ. The same holds for For δ, ε small enough, we suppose thatx,ȳ belong to the neighbourhood of the boundary where the distance is smooth. By (5.24) and the smoothness of the distance function we have and since where C > 0 is a constant depending on B, K and ||Dφ|| L ∞ (BC j (x0)) . By coupling (5.30), (5.32), (5.10) and (i) of Lemma 5.2, we get where o δ,ε (1) means lim δ→0 o δ,ε (1) = o ε (1). Plugging (5.33) and (5.26) into (5.25), we get Step 3. Sending δ → 0. We want to send first δ → 0 in (5.34) and we observe that the nonlocal terms are uniformly bounded in δ. Consider I ξ ′ [u](x), observing that the same argument works similarly for I ξ ′ [v](ȳ). Note that by (5.23) and ≤φ(x + j(x, z),ȳ) −φ(x,ȳ). (5.35) Note that χ δ is Lipschitz with Lipschitz constant independent of δ thanks to (5.24). Then, by the definition ofφ, since χ ε , χ δ , φ are Lipschitz and by (J1), we have which, by the first of (M), gives the uniform boundedness in δ of I ξ ′ [u](x) when |z| < 1. When |z| ≥ 1, the claim simply follows by the boundedness of u and the first of (M). Then, we send δ → 0 in (5.34) and, by the semicontinuity and boundedness of u and v and Lemma 5.2 (i), apply Fatou's Lemma and we get In the next step, we estimate the nonlocal terms by bounds independent of the parameter ξ ′ small. Indeed, at the end of the proof, we will send first ξ ′ → 0 and then ε → 0. Note also that now, thanks to Lemma 5.2 (ii), we have that d(x) = d(y).
Step 4. Estimate of the nonlocal terms. We prove the following lemma.
Lemma 5.4. Under the above notation, we have where K ξ , P ξ are as in (5.47) and (5.48) and C > 0 is independent of all the parameters.
Proof. For simplicity of exposition, we first conclude the proof when the measure µ in the nonlocal terms has no dependence on x, i.e. µ x ≡ µ. We refer to Remark 8 for details in the case of x-dependence. We write where J x = {z ∈ R n | x + j(x, z) ∈Ω} and Letting δ → 0 in (5.42), by (i) of Lemma 5.2, the definition ofφ and the semicontinuity of u, v, we get for any If |z| < 1, then by (5.23) and (J1), x + j(x, z), y + j(y, z) ∈ B 2Cj ξ ′ (x 0 ). Then we write (5.43) for x ′ = x + j(x, z), y ′ = y + j(y, z) and we have Note that by the Lipschitz continuity of χ ε , (J1) and (5.9), we have where D j is defined in (J1) and then Then by (5.45) and the first of (M), we get by (J1), the first of (M) and the Lipschitz character of φ. Now we consider the term I ξ ′ [J x /J y ], defined in (5.40), observing that the same argument works similarly for I ξ ′ [J y /J x ]. Take 0 < δ 0 < 1 enough small (note that δ 0 will be defined more precisely at the end of the proof of Lemma 5.38). We split the domain of integration in {z : |z| ≥ δ 0 } and {z : ξ ′ ≤ |z| ≤ δ 0 }. We write By the boundedness of u, we have and since |J x /J y | → 0 as ε → 0, by the first of (M) and the Dominated Convergence theorem, we get ] we use again the maximum point inequality (5.43) with x ′ = x + j(x, z), y ′ = y and since φ ∈ C 1 and by (J1), the first of (M), we get where we remark that C > 0 is independent of all the parameters. We couple (5.39), (5.46), (5.49), (5.50), (5.51) and (5.52) with (5.39) and we get (5.53) where for all x, y ∈ R N , we denote Then the claim of the lemma follows by plugging (5.55) into (5.53).
Note that Lemma 4.1 is not necessary when dealing with domains with flat boundary. In the following remark we consider the case when Ω is the half space and we show how the estimate of the nonlocal terms can be carried out more easily without Lemma 4.1.
Consider the nonlocal terms in (5.25) and restrict ourselves to a subsequence such thatx N ≥ȳ N (ifx N ≤ȳ N the argument is similar). Then we can write where in the last line we used the same notation as in the previous step, see in particular (5.40), (5.41). The term T ξ ′ [Jx ∩ Jȳ] is treated exactly as in the non flat case (see the previous step). On the contrary, note that in this case the estimate of the term I ξ ′ [Jx/Jȳ] is easier, since by (5.56) |Jx/Jȳ| → 0 as δ → 0 and then by the Dominated Convergence Theorem, we have I ξ ′ [Jx/Jȳ] → 0 as δ → 0.
Step 5. -Sending the other parameters to their limits. We couple (5.37) with (5.38) and we get where C > 0 is independent of the parameters. Then, we first send ξ ′ → 0 by the Dominated Convergence Theorem and we get where C is a constant independent of ξ. Moreover, since φ is C 1 , by the first of (M), the Dominated Convergence Theorem and since x, y → x 0 as ε → 0, we have lim sup and by the boundedness and semicontinuity of u, v and Lemma 5.2 (i), we apply Fatou's lemma for each ξ > 0 fixed and we get lim sup and, by the previous estimates, we conclude by sending ε → 0 in (5.57).
Remark 8. We give some details of the analysis of the nonlocal terms in step 4 when the measure µ depends on x. We write (5.39) with For I ξ ′ [J x /J y ] and I ξ ′ [J y /J x ] we proceed as above (Step 4), noting that the xdependence plays no role by the first of (M). For the T -term, we write For T ξ ′ 1 [J x ∩ J y ], we proceed as above (in Step 3, for T ξ ′ [J x ∩ J y ] defined in (5.41)) and we prove (5.46). Now consider For T ξ ′ 2 [B ξ ] we use the maximum point inequality (5.43) and we write for |z| ≤ ξ Then by the Lipschitz continuity of χ ε and φ, (J1), (M) and (5.9) we get where we observe o ε (1) is independent of ξ ′ and from now on may change from line to line in the following. For T ξ ′ 2 [B c ξ ], by the boundedness of u, (M), (5.9), we get Then, by (5.59) and (5.60), we get From now on the proof is the same as above.
Remark 9. We give the details of the proof of Proposition 5.1 in case (a), when The strategy of the proof relies on the existence of a blow-up supersolution exploding on the boundary, which allows us to keep the maximum points away from the boundary. The existence of such a supersolution is stated in the following lemma, whose proof is given in the Appendix.
Lemma 5.5. For anyx ∈ Γ in , there exists r = r(x) > 0 and a positive function U r ∈ C 2 (B r (x) ∩ Ω) satisfying for any ξ small enough (with respect to r, that is, in B r (x)∩Ω, for some function ω r which is non-negative, continuous, strictly increasing in a neighbourhood of 0 and satisfies ω r (0) = 0.
Proof of case (a). Let r = r(x 0 ) be defined in Lemma 5.5 forx = x 0 . We localize the argument in a ball of radius r around x 0 and we use the existence of the blowup function U r defined in Lemma 5.5 forx = x 0 . Let 0 < ξ ′ < min{ξ, C −1 j r 4 } and ε > 0. We double the variable and we consider (x, y) maximum point on Note that, by (ii) of Lemma 5.5, we have that (x, y) ∈B Cj ξ ′ (x 0 )∩Ω×B Cj ξ ′ (x 0 )∩Ω; moreover, again by (ii) of Lemma 5.5, we have for k small enough Note that the existence of the blow-up function plays its mayor role here to get (5.62). This estimate tells us, roughly speaking, that the maximum points are away from the boundary. For fixed k, a standard argument shows that By the previous estimate on x, y and extracting subsequences if necessary, we can assume, without loss of generality, that as ε, k → 0 Thanks to (5.64), we can take ε, k small enough so that x, y ∈ B r 4 (x 0 ) ∩ Ω. We proceed as in Step 2 in the above proof, we write the viscosity inequalities (5.25) and, using thatφ ∈ C 1 , (J1) and the first of (M), we get . For simplicity of exposition, we conclude the proof in the case the measure µ in the nonlocal terms has no dependence on x, i.e. µ x ≡ µ. The result can be easily extended in the case of x-dependence analogously as already shown in Remark 8 for case (b) and (c). We use the same notation of Proposition 5.1, see (5.40), (5.41) and we write As in the proof of Proposition 5.1, the term T ξ ′ [J x ∩ J y ] can be estimated as follows where we use the notation (5.47), (5.48) of Proposition 5.1. Note that o ε (1) is independent of ξ ′ . Since U r is Lipschitz, by (J1) and the first of (M), we have where o ε (1) is independent of ξ ′ . Now we estimate the terms I ξ ′ [J x /J y ] and I ξ ′ [J y /J x ]. Thanks to (5.62), in this case the estimate is easier than in the cases (b) and (c) treated above. Take for example I ξ ′ [J x /J y ] (the argument being analogous for I ξ ′ [J y /J x ]) and note that by (5.62) the integral is independent of ξ ′ as soon as ξ ′ <δ whereδ is defined in (5.62). Then by the boundedness of u, we have and since |J x /J y | → 0 as ε → 0, by the first of (M), the Dominated Convergence theorem, we get where o ε (1) is independent of ξ ′ . Then plugging (5.68) and (5.67) into (5.66) and then coupling it with (5.65), we get where o ε (1) is independent of ξ ′ . Now, by (i) of Lemma 5.5, we estimate the integral terms of the left hand side of (5.69) together with the first order terms involving U r in H(y, −D[φ(x, ·)](y)) − H(x, D[φ(·, y)](x) and we get Then, plugging (5.70) into (5.69), we get where o ε (1) is independent of ξ ′ . The rest of the proof is the same as in the previous cases, by sending first ξ ′ → 0, then ε → 0 and finally k → 0. For the details we refer to the end of the proof of Proposition 5.1. for 0 ≤ s < 1 2 , increasing for 1 2 ≤ s < 1, ||u|| ∞ + ||v|| ∞ + 1 for s ≥ 1. and d is the signed distance from the boundary (see Remark 1). Note that sup Φ → M as R → ∞ and ν → 0. Since Φ ≤ −1/2 for |x| large and ν small enough and M > 0, the function Φ achieves its positive maximum sup Φ > M 2 at a point x for R big and ν small enough. We give the details in the case where all maximum points x are located on the boundary. We have We use φ(·) := ψ(R −1 |·|)−νd(·) as a test function at x. Note that, if x ∈ ∂Ω and for ν > R −1 ||ψ ′ || L ∞ , we have ∂φ ∂n ≥ −R −1 ||ψ ′ || L ∞ + ν > 0. Then, by Proposition 5.1, we get

Coercive Hamiltonians.
We proceed analogously as for Hamiltonians of Bellman type and we prove Proposition 5.6. Once proved Proposition 5.6, the proof of Theorem 3.1 for H coercive follows by standard arguments as already showed for Hamiltonians of Bellman type. We sketch first the proof of Theorem 3.1 and then we prove Proposition 5.6.
Proof of Theorem 3.1. We just observe that we proceed again by contradiction, supposing that M = sup Ω {u − v} > 0. We fix 0 < µ < 1 and define for x ∈ Ω ω µ (x) = µu(x) − v(x). We proceed as in the Bellman case and we use Proposition 5.6 to get µ). Then, by letting R → ∞, ν → 0 and finally µ → 1, we get a contradiction since M > 0 and we conclude the proof.
where A, C m,µ are positive constants which depend on the data. Precisely, if µ = 1, Proof. We give the details when H has superlinear form i.e. when m > 1 > σ, since the proof in the sublinear case is similar with easier computations. Since the proof is similar to that of Proposition 5.1, we focus only on the main differences.
We start by noting that if u is a subsolution of (1.1), thenū = µu is a viscosity subsolution to Let x 0 ∈Ω and φ a smooth function such that ω − φ has a strict maximum point inB Cjξ (x 0 ) ∩Ω at x 0 for some 0 < ξ < 1. We suppose that x 0 ∈ ∂Ω, since the other case being similar and even simpler.
Step 1. -Localising on equidistant points (i.e. d(x) = d(y)). We double the variable and consider the function whereφ is as in (5.3) with K = 2. Let (x,ȳ) be a maximum point of Φ over the set , where ψ is defined as in (5.8) with K = 2 and let (x, y) be a point of maximum of Ψ over A. By classical argument, we have as ε → 0 We have the following lemma.
Lemma 5.7. Under the above notation, we have Proof. The proof is very similar to that of Lemma 5.2. We give a sketch of the proof of (ii). We suppose that ∂φ ∂n (x 0 ) > 0, then for ε small enough ∂ψ ∂n ((x + y)/2) > 1 2 ∂φ ∂n (x 0 ) > 0. We assume that d(x) > d(y). Then for 0 < µ < 1, we have where C > 0 is a constant independent of ε. For the integral terms in (5.78) we proceed as in Proposition 5.1 by using Lemma 5.3. For the Hamiltonian terms we use assumption (H1) together with (5.79) and we get which is a contradiction for ε small, since σ < m.
Step 2. -Writing the viscosity inequalities. We write the viscosity inequalities for u and v where we estimated the I ξ ′ [φ] terms by o ξ ′ (1) (independent of δ) as in (5.26). Denote We recall that By the smoothness of the distance function Thanks to (Hb) and (Hc), we get for ε small Note that, by (5.77) and (iii) of Lemma 5.7, we can take ǫ = ǫ(µ), δ = δ(µ) small enough so that and we can write where o δ,ε (1) means that lim δ→0 o δ,ε (1) = o ε (1). Note that the infimum in the previous expression is attained and therefore Then we couple (5.86), (5.81) and (5.80), we let δ → 0 and we get Step 3. Estimate of the nonlocal terms. In order to estimate the nonlocal terms we use the following lemma. We omit the proof since it is exactly the same as that of Lemma 5.4. We remark that in the proof we deeply rely on the assumption σ ∈ (0, 1).
Lemma 5.8. Under the above notation, we have where C > 0 is independent of all the parameters.
Then the rest of the proof easily follows as in Proposition 5.1, step 5.
6. Applications to evolutive problems: existence, uniqueness and asymptotic behavior In this section we present some applications of our results proved in the previous sections for the stationary case to the evolutive setting and we consider the associated Cauchy problem where Ω ⊂ R N satisfies assumption (O), u 0 ∈ C(Ω), I[u] is an integro-partial differential operator of censored type and of order strictly less than 1, defined as where µ x is a singular non-negative Radon measure satisfying (M), and j(x, z) is a jump function satisfying (J0), (J1) (see Section 2) . The main example are measures µ x with density dµx dz = g(x, z)|z| −(N +σ) with σ < 1 and g a non-negative bounded function Lipschitz in x uniformly with respect to z. Note that the operator (6.2) is the natural extension to the evolutive case of the nonlocal operator considered in the stationary case and defined in (1.2). Moreover H :Ω×[0, +∞)×R×R N → R is a continuous function whose growth in the gradient makes it the leading order term in the equation, and can be coercive or of Bellman type. We refer to the following section for the precise assumptions on the Hamiltonian. The well-posedness of problems as (6.1) follows from analogous arguments used for the stationary problem with some standard adaptations. We refer to Theorem 6.1 and Theorem 6.3 for further details and proofs of uniqueness and existence.
We study two different kind of asymptotic behaviour of the solutions of (6.1). First we consider a Hamiltonian either of Bellman type either coercive under assumptions ensuring uniqueness of the solution of the associated stationary problem. We prove, by classical methods based on the weak-relaxed semilimits, the convergence as t → +∞ of the solution of (6.1) to the unique solution of the associated stationary problem. On the other hand, when the associated stationary problem has not unique solution, we consider a Hamiltonian with superfractional coercive growth and we study the so-called ergodic large time behaviour, proving that the solution of (6.1) approaches a solution of the so-called ergodic problem as t → +∞. We follow the methods of [9], which rely on the Hölder regularity in Ω of the subsolutions of the associated ergodic problem and on a Strong Maximum Principle. We recall also [5], where analogous methods were applied to uniformly elliptic integro-differential equations. We refer to subsection 6.4 for more details.
6.1. Assumptions. We are going to consider the finite time horizon problem associated to (6.1) The definition of viscosity solutions to (6.3) (and then to (6.1)) is the natural extension of Definition 2.1 to the corresponding Cauchy problem. We assume the following condition.
(H') For all R > 0, there exists γ R ≥ −K such that for all x ∈Ω, u, v ∈ R, |u|, |v| ≤ R, 0 ≤ t ≤ R and p ∈ R n , we have Remark 10. Note that in the comparison principle (Theorem 6.1), assumption (H') is assumed mainly for simplicity of exposition and can be relaxed assuming only γ R ≥ 1. Indeed if γ R < 1 we perform the changeũ = ue −(γR−1)t (analogously for the supersolution) and prove Theorem 6.1 forũ andṽ.
As it is classical in viscosity solutions theory, the comparison principle allows the application of Perron's method to conclude the existence. To this end, we introduce the following assumption, which will allows us to build sub and supersolutions: We consider Hamiltonian either of Bellman type either coercive.
We say that the Hamiltonian H is of Bellman type if for t ∈ [0, +∞), x ∈Ω, p ∈ R N , H(x, t, r, p) can be written as (C') Uniform continuity of l and λ: There exist modulus of continuity ω l , ω λ such that such that ∀α ∈ A, ∀x, y ∈ Ω, t, s ∈ [0, +∞) |l(x, t, α) − l(y, s, α)| ≤ ω l (|x − y| + |t − s|); |λ(x, t, α) − λ(y, s, α)| ≤ ω λ (|x − y| + |t − s|); (L') Uniform Lipschitz continuity of the drift b: There exists C > 0 such that Remark 11. Note that assumption (L') may seem unusual since it requires also the uniform Lipschitz continuity of the drift b in the time variable. This is due to the fact that in the proof of the comparison principle (Theorem 6.1) we need to double also the time variable in the test function, which as a consequence will have the same dependence on x and t, which we estimate by assumption (L'). We refer to Remark 12.
We assume that the components of the boundary are of three different types whose definition is analogous to that of Γ in , Γ out , Γ defined in (2.3), (2.4), (2.5) in the stationary case, with the only difference that now they are part of the parabolic boundary. For simplicity of exposition, we do not repeat the definition. For the sake of simplicity and for the rest of the paper, we adopt the following notation. We assume that the Hamiltonian satisfies in x uniformly in t and uniformly on compact sets in r the assumptions made in the stationary case, namely, (H1), (Ha) (sublinear coercivity) or (H1), (Hb), (Hc) (superlinear coercivity). For the sake of brevity and simplicity of exposition, we omit to repeat them and for the rest of the paper we adopt the following notation. or (He"). Assume also (H'). Let u, v ∈ L ∞ (Ω × [0, T ]) for all T > 0 be respectively a bounded usc sub and bounded lsc supersolution of (6.1). Then The proof is based on the following proposition, which is the analogous for the evolutive case to Propositions 5.1 and 5.6. The proof is also similar with some adaptation to the evolutive setting. Note that Proposition 6.2 is used also later in the proof of Proposition 6.5.

Remark 12.
We omit the proof of Theorem 6.1, since it follows closely the arguments presented in Proposition 5.1 and Proposition 5.6 for the stationary case adapted by standard arguments to the evolution setting. We only remark that, since the strategy of the proof of Proposition 5.1 relies on an asymmetric use of the viscosity inequalities satisfied by the sub-and supersolution (see in particular in Step 1), the standard approach by the Ishi's Lemma for evolutive equation is not applicable. Then, we need to double also the time variable in the test function, which as a consequence will have the same dependence on x and t. Because of this, in particular in order to estimate the Hamiltonian terms in the viscosity inequalities satisfied by u and v (in the Bellman case), we need the uniform continuity also with respect to time as stated in assumption (L').
For both the coercive and Bellman case, the application of Perron's method on a sequence of finite-time horizon problems with the form (6.3) with T → +∞ and the strong comparison principle allows us to get the existence of a solution which is defined for all time. Note that, in order to apply Perron's method, we ask the initial datum to be bounded, i.e. u 0 ∈ C b (Ω).
Moreover, in order to have the uniform boundedness in T of the solutions of (6.1), we suppose the following assumption.
(H") There exists γ 0 > 0 such that γ R ≥ γ 0 for all R > 0 where γ R is defined in (H'). Moreover, the unique viscosity solution u of (6.1) converges uniformly on compact sets inΩ to u ∞ , the unique viscosity solution of the problem (6.5).
For the existence and uniqueness for the problem (6.5) we refer to Theorem 3.1 and Corollary 3.2. We omit the proof of Theorem 6.4 since it is rather classical and follows the same arguments used in [13], where the same kind of results have been given in the case of the Dirichlet problem for nonlocal equations (fractional Laplacian) with Hamiltonians both coercive both of Bellman type. 6.4. Large time behavior II: convergence to the ergodic problem. In this subsection we prove large time behaviour for the problem where Ω is a bounded open subset of R N satisfying (O), u 0 ∈ C(Ω) and H is a Hamiltonian in superfractional coercive form, that is, m > 1 in (H1'). Existence and uniqueness for the problem (6.6) follow from Theorem 6.3. Note that H does not depend on u, so it does not satisfy (H').
The main result of this section is Theorem 6.9, namely the convergence as t → +∞ of the solution of (6.6) to a solution of the ergodic problem, which we solve in Proposition 6.8. The proof of Theorem 6.9 strongly relies on the Hölder regularity up to the boundary for subsolutions of (6.6) and on the control of their oscillation. This result was proved in the stationary setting by Barles, Ley, Koike, Topp in [9], Corollary 2.14 (see also Barles and Topp [12]), in the case of censored operators and coercive Hamiltonians with m > 1. We recall this result in Proposition 6.7.
For the evolutive problem we refer to the proof of Theorem 4.5 of [9], where first Lipschitz regularity in time is proved and then Corollary 2.14 is used to conclude the regularity in time and space.
We remark that, differently from [11], where the Lipschitz regularity of the solutions is used to linearise the equations in order to apply the Strong Maximum Principle, our proof relies mainly on the use of a Strong Maximum Principleà la Coville [17], [18] (see also Ciomaga [16]). This means that it relies mainly on a topological property of the support of the measure defining the nonlocal operator. Note that in this final part we assume Ω bounded for technical reasons related to the proof of the Strong maximum principle, we refer to the proof of Proposition 6.5.
6.4.1. A strong maximum principle. We need some notation for the statement of the Strong Maximum Principle. Let µ, j be as in the definition of the nonlocal operator I, that is, satisfying (M), (J0), (J1) and denote by suppµ the support of the measure µ. For x ∈ R n we define inductively The Strong Maximum Principle presented in this paragraph relies in the non locality of the operator under the "iterative covering property" (6.7) X (x) = Ω, for all x ∈ Ω.
The most basic example is the case where j(x, z) = z and there exists r > 0 such that B r ⊂ supp{µ}. For further details and examples we refer to [9]. The following proposition states the Strong Maximum Principle. . Let u, v be respectively a bounded sub and supersolution of (6.6), such that there exists (x 0 , t 0 ) ∈Ω × (0, +∞) satisfying Then, the function u − v is constant inΩ × [0, t 0 ]. Moreover we have The proof of Proposition 6.5 uses the following lemma, which is a consequence of the comparison principle, see [11], Theorem 4.1.
Lemma 6.6. Let assumptions of Proposition 6.5 hold. Let u, v be bounded sub and supersolution to equation (6.6) and for t ∈ [0, +∞) define Then, for all 0 ≤ s ≤ t, we have k(t) ≤ k(s).
We prove the strong maximum principle. We essentially follow the argument of [9], Proposition 4.1, with some changes due to presence of the Neumann boundary condition. We give a sketch of the proof, focusing on the main differences.
Proof (Proof of Proposition 6.5). We divide the proof into several steps.
Step 1. -Preliminaries. We want to prove that (u − v)(x, t) = k(0) for each (x, t) ∈ Ω × [0, t 0 ]. Since k(t 0 ) is a global maximum value of k in [0, +∞), by Lemma 6.6 we have k(t) = k(0) for all t ∈ [0, t 0 ]. Then, we have just to prove that for each τ ∈ (0, t 0 ). By upper-semicontinuity, we derive the result up to τ = 0 and τ = t 0 . Fix τ ∈ (0, t 0 ) and define the set We observe that by the upper-semicontinuity of u − v, B τ is non-empty. Then, the claim of the proposition follows once proved that B τ =Ω.
Step 4. -Writing the viscosity inequality for ω µ . Denote We test ω µ with the function φ in (x µ , t µ ). We suppose x µ ∈ ∂Ω, since the other case being analogous and even simpler. By the Taylor's formula for the distance function, we have and then Take µ, ε small enough so that By (6.11), by the definition of K and (6.12), we have Then, the rest of the proof follows exactly as in [9], Proposition 4.1, by writing the viscosity inequality for ω µ , letting first µ → 1, then ε → 0 and using the iterative covering property (6.7). We give some details for completeness of exposition. By Proposition 6.2, we get for ξ > 0 where C µ,m and A R are defined in Proposition 6.2. Since We recall that t µ → τ as µ → 1 and we observe that by the smoothness of ψ ε,α the term in parenthesis in (6.13) remain bounded as µ → 1. Moreover, by the Dominated Convergence Theorem, we get where we recall thatx is the limit of x µ as µ → 1. Then and letting ε → 0 and recalling thatx → x τ as ε → 0 and (u − v)(x τ , τ ) = k(τ ) we finally conclude Since ξ > 0 is arbitrary, we get that Therefore we can proceed in the same way as above, and conclude by induction that Then we conclude the proof by the upper-semicontinuity of u − v and applying the iterative convering property (6.7).
We solve the ergodic problem in Proposition 6.8. Note that we need the compactness of the family of solutions {u δ }, which relies mainly on the regularity result for subsolutions of (6.14) which we recall in the following proposition. We remark that the Hölder regularity for subsolution of the first equation of (6.14) up to the boundary was proved in [9], Theorem 5.5. Note that here the Neumann boundary conditions play no role, since the regularity inside the domain follows by [9], Theorem 5, and the extension up to the boundary is carried out similarly to [9], applying the method used by Barles in [3] (see also [15]). Proposition 6.7. Let assumptions of Proposition 6.5 hold. Then any bounded viscosity subsolution u : R N → R to (6.14) is Hölder continuous inΩ with Hölder exponent γ 0 = m−σ m and Hölder seminorm depending on Ω, the data and osc Ω (u) and not on δ. Moreover, there exists K > 0 such that for any bounded viscosity subsolution of (6.14) we have (6.15) osc Ω (u) ≤ K.

Proof.
A key ingredient is Proposition 6.7, which gives the compactness of the family of solutions of the approximating equation (6.14). Once we have the compactness, the proof follows standard arguments which we do not repeat. The uniqueness follows by the comparison principle for (6.6) and the application of the strong maximum principle for the problem (6.6) (Proposition 6.5). 6.4.3. Convergence as t → +∞. Theorem 6.9. Let assumptions of Proposition 6.5 hold. Let u be the unique solution to problem (6.6). Then, there exists a pair (w, λ) solution to (6.16) such that u(x, t) − λt − w(x) → 0 as t → +∞, uniformly onΩ.
We omit the details of the proof since we follow closely the arguments given in [9] (see also [4], [13] for the local framework and [27] for the nonlocal one).
We just observe that again a crucial ingredient of the proof is the Hölder regularity of the solutions of (6.6). The proof is carried out as in Theorem 4.5 of [9], by first proving Lipschitz regularity in time and then applying Proposition 6.7 to conclude the regularity in time and space. Once established the regularity, the proof follows by the application of the strong maximum principle proved in Proposition 6.5.
We state Lemma 6.11 and we omit the proof since it follows by standard arguments.
Let d denote the distance from the boundary of Ω in a neighbourhood V of the boundary and extend it C 1 and bounded in all the domain. Let I as in (1.2) and assume (M) and (J1). Then the function φ = ψ R + νd satisfies 6.6. Blow-up supersolution. In this section we prove Lemma 5.5, used in the proof of Theorem 3.1 (see in particular Proposition 5.1, Remark 9). We construct the function U r as showed in the following and then we prove Lemma 6.12. Note that Lemma 5.5 follows as a consequence of Lemma 6.12 and we prove it after the statement of Lemma 6.12.
Letx ∈ Γ in and r = r(x) be given as in assumption (O). We recall that by (O), there exists a W 2,∞ -diffeomorphism ψ : B r (x) → R N , satisfying (6.27) ψ N (s) = d(s) for any s ∈ B r (x), where d is the signed distance from the boundary of Ω. We define U r in a suitable neighbourhood of the pointx, where we rectify the boundary and carry on the computations. In particular we define U r as follows (6.28) U r (x) =Ū r (d(x)) for x ∈ B r (x) ∩ Ω, whereŪ r (s) = − log(s) + 3 2 log r if 0 < s ≤ r.
Lemma 6.12. For anyx ∈ Γ in , let r = r(x) be defined as in assumption (O). Let U r be defined as in (6.28). Then we have for ξ small enough (with respect to r) Remark 13. Note that the strict positivity of the drift term on the points of Γ in is essential here to prove (6.30), since the drift term controls the integral term which explodes on the boundary, as (6.29) shows.
As a consequence of Lemma 6.12, we prove Lemma 5.5.
Proof of Lemma 5.5. Taker as defined in Lemma 6.12 and let Ur be defined as in (6.28) for r =r. Then Ur is a non-negative decreasing function which trivially satisfies (ii) of Lemma 5.5 with ωr(s) = 1 Ur (s) . Moreover, (i) Lemma 5.5 follows as a direct application of (6.30) of Lemma 6.12. Now we prove Lemma 6.12.
We proceed performing a change of variable in order to write the set of integration in terms of ψ N (x). In other words, we write ψ(x + j(x, z)) − ψ(x) = w.
Then by (J0), (J1), the first of (M) and since ψ is W 2,∞ , (6.31) becomes for someC, C > 0. By the change of variable y = w ψN (x) , we get (6.32) Note that the integral in the right hand side is finite and does not depend on ξ. For convenience of notation we denote A :=C 0≥yN ≥−1 | ln (1 + y N ) | dy |y| N +σ . Then (6.32) becomes which is exactly (6.29). Now we prove (6.30). First note that, by the definition of U r , we have (6.33) DU r (x) = d(x) −1 n in B Cjr (x) ∩ Ω.
Since we are in a neighbourhood of Γ in , σ < 1 and A is compact, there exists 0 <r < r (depending only on A, σ and r), such that if x ∈ B r 2 (x) ∩ Ω b(x, α) · DU r (x) + I ξ [U r ](x) ≤ 0 ∀α ∈ A.
Then (6.30) follows and we conclude the proof of the Lemma.

Acknowledgement
This work has been carried out during the permanence of the author at the Laboratoire de Mathématiques de Tours, which the author wants to thank for the warm welcome. In particular the author wants to express her deepest gratitude to Prof. Guy Barles, for proposing the problem and for the important support given. Without him, this work would not have been carried out. The author kindly thanks also the anonymous referees, whose comments and suggestions were precious in the improving of the final presentation of the paper.