On the dynamics of vortices in viscous 2D flows

We study the 2D Navier--Stokes solution starting from an initial vorticity mildly concentrated near $N$ distinct points in the plane. We prove quantitative estimates on the propagation of concentration near a system of interacting point vortices introduced by Helmholtz and Kirchhoff. Our work extends the previous results in the literature in three ways: The initial vorticity is concentrated in a weak (Wasserstein) sense, it is merely $L^p$ integrable for some $p>2$, and the estimates we derive are uniform with respect to the viscosity.


Introduction
Coherent flow structures are widely observed in two-dimensional turbulent fluid motions.Isolated regions of concentrated vorticity emerge, for instance, as a result of stirring or solid body interactions, and they persist for rather long time scales while moving with the ambient background flow and interacting with nearby vortex regions.Vortex regions that are sufficiently far from each other can be considered as rigid disc-like objects, whose interactions reduce approximately to the interactions of their centers.
The idea of studying the motion of vortices in two-dimensional fluid flows by means of an idealized point vortex system dates back to the second half of the nineteenth century and the work of Helmholtz [19].Arguing formally, Helmholtz derived an ODE model for interacting point vortices in inviscid incompressible fluids.Kirchhoff later documented the Hamiltonian structure of this system [22], and, as a result, the ODE model is nowadays commonly referred to as the Helmholtz-Kirchhoff point vortex system.The Helmholtz-Kirchhoff system is a reasonable description of the actual vortex dynamics of isolated and highly concentrated vortex regions.In these situations, if the fluid kinematic viscosity is small, it seems quite natural to expect that the vortices will exhibit a behavior similar to the Helmholtz-Kirchhoff idealization.Likewise, it is reasonable to expect that this description will become invalid in the event of collisions and merging of vortices or if vortex regions dissolve in the ambient background flow as a result of viscous friction.The goal of the present paper is to derive estimates that describe to what extent the point vortex system captures the motion of vortex regions in incompressible viscous flows.
For inviscid fluids, the first rigorous work in this direction is due to Marchioro and Pulvirenti [27], who proved that bounded vortex patch solutions to the Euler equation remain close to the point vortices during the evolution.Their strategy exploits the regularity of the velocity field sufficiently far from the vortex patches and the symmetry properties of the two-dimensional Biot-Savart kernel.The original work was gradually improved in many subsequent articles, see for example [28,26,3].Other methods have been employed to study the connection of the point vortex system with the Euler equation.We mention the approach by Turkington [31], which relies on the conservation of the energy.This has the advantage that it works also in situations where the Biot-Savart operator does not possess certain symmetry features, for example in the setting of three-dimensional axisymmetric fluids without swirl, where it has been successfully applied [2].It does not seem however suited to studying several interacting vortices; in this case a combination with the Marchioro-Pulvirenti method has been used [4,5].Another technique effectively implemented to this context is the so-called gluing method [13,14].In these works, solutions with highly concentrated vorticities containing precise information on the vortex cores are constructed in the two-dimensional case, and helical vortex filaments are built in the three-dimensional setting.
In the context of viscous fluids, the connection with the Helmholtz-Kirchhoff point vortices has been shown for the Navier-Stokes equation [24,25,9], where the authors consider bounded vorticities initially sharply concentrated in separated regions and prove that they converge in the inviscid limit to the point vortices.A more singular initial datum was considered by Gallay [16], namely a collection of several point vortices.He proved that the corresponding Navier-Stokes vorticity concentrates, in terms of a weighted L 2 norm, near a combination of Lamb-Oseen vortices centered around a viscous regularization of the point vortex system (see also the nice review paper [17]).A Lamb-Oseen vortex (see expression (11)) is the solution corresponding to an initial vorticity concentrated on a single point, and it can in some ways be regarded as the fundamental solution to the two-dimensional Navier-Stokes vorticity equation.Gallay gives an estimate on the rate of concentration, proving that it is proportional to νt, where ν is the viscosity and t is the time.As consequence of this result, the convergence of the Navier-Stokes solution to the Helmholtz-Kirchhoff point vortex system for ν → 0 is also established.
It would be desirable to establish a relation between solutions to the Navier-Stokes equations and the (viscous) point vortex dynamics that has both features: First, it provides stability estimates in the sense of the Marchioro-Pulvirenti strategy showing that solutions that start close to point vortices remain close to vortices.Second, it describes the accurate shape of vortices that are actually spreading as a result of viscosity.
In the present work, we take a tiny step towards this goal by addressing the first property and by improving in several aspects the aforementioned works [24,25] by Marchioro on the Navier-Stokes equations: As a measure for vortex concentration, we work with Wasserstein distances.Weak notions of concentration are suitable in the viscous setting due to the fact that the viscosity spreads instantaneously over the full space.However, also in the inviscid case, weak notions of concentration are favorable as they allow for more general vortex configurations like elongated vortex regions or even long tentacles, as can be observed in chaotic or turbulent flows.We have considered concentration in the Wasserstein sense in our earlier works on the Euler equation [7,8].In the context of three-dimensional Euler filaments, geometric flat norms, which are related to Wasserstein distances in 2D, were considered [21].
Furthermore, we establish estimates that hold uniformly in the viscosity constant, independently from the scale of concentration.This way our results apply also to initial data that are given by a collection of point vortices as studied by Gallay [16], and our estimates describe the sharp rates of vortex spreading due to viscosity.
We finally consider configurations whose initial vorticity is allowed to be unbounded but in L p , for some p > 2, with almost no assumption on the magnitude of the L p norm.This same condition was studied [8] in the Euler setting and extends the corresponding requirement on the L ∞ norm that was considered in [25,9].
Organization of the paper.In Section 2 we introduce precisely the mathematical setting and state the main results.In Section 3 we present the proofs.

Mathematical setting
We study the dynamics of vortices in a viscous incompressible fluid in the twodimensional plane Ê 2 .This motion can be modelled via the Navier-Stokes equations in vorticity form.These read where u : (0, +∞) × Ê 2 → Ê 2 is the velocity field and ω : (0, +∞) × Ê 2 → Ê is the scalar vorticity field that measures the tendency of the fluid to rotate.The incompressibility condition is translated into the mathematical constraint that the velocity field has zero divergence, (2) The vorticity field can be computed from the velocity as the rotation ω = ∂ 1 u 2 −∂ 2 u 1 , and vice versa the velocity u can be recovered from ω via the Biot-Savart law where the convolution is understood in space.Here, the function K is the Biot-Savart kernel and z ⊥ denotes the counter-clockwise rotation of the vector z by 90 degrees.Finally, the constant ν > 0 is the kinematic viscosity of the fluid.We consider a compactly supported initial datum ω ∈ L p (Ê 2 ) for some p > 2. It is well-known (see [17] and references therein) that there exists a unique solution ω ∈ C (0, +∞); L 1 ∩ L ∞ (R 2 ) to equation (1).Such a solution has been shown to be smooth in space and time [1], and its Lebesgue norms are non-increasing in time, that is We study the situation where the initial vorticity is split into N components of definite sign where for every i ∈ {1, . . ., N} it holds either ωi ≥ 0 or ωi ≤ 0. The motion of the vorticity components obeys the advection-diffusion equation Hence their sign is preserved over time and, because of the uniqueness of the solutions to linear advection-diffusion equations, the vorticity ω remains expressed as the sum of its components, We observe that, thanks to the divergence-free condition (2) and to the absence of domain boundaries, the intensity, that is the space integral of each vortex component, is constant in time, Depending on the sign of the vorticity components, we will either have that a i > 0 or a i < 0, so that ω i /a i is a well-defined probability distribution.We suppose that the vortex components are sharply concentrated around N distinct points Ȳ1 , . . ., ȲN ∈ Ê 2 in the sense that for all i = 1, . . ., N. The parameter ε > 0 is the concentration scale, that will be assumed to be as small as needed.Here W 2 is the 2-Wasserstein distance, which, in the case that one of the measures considered is atomic, is just the square root of the variance, and measures the average "size" of the vortex region.We refer to Villani's book [32] for more information on Wasserstein distances.To understand the role of ε, notice that condition (4) is satisfied for example if ωi is supported in a ball of radius ε centred on Ȳi .An easy computation shows that the Wasserstein distance is minimized when the atomic measure is located on the center of vorticity Xi = 1 Since Wasserstein distances metrize the weak convergence in the sense of measures, see Theorem 7.12 in [32], assumption (4) implies that the rescaled vortex component ω i /a i converges weakly as ε goes to zero to an atomic measure.We suppose that the intensities a i are independent of ε and, accordingly, the L p norm of the components must diverge as ε → 0. Here, we assume that where γ is a fixed positive number.This is the same condition that we considered earlier in the inviscid setting [8], while similar requirements on the L ∞ norm were studied e.g. in [24,25,9] for the Navier-Stokes equations and in [26,6] for the Euler equations.Notice that the two assumptions (4) and (6) enforce that γ ≥ 2 − 2/p as a consequence of the interpolation 1 f p L p W 2 (f, δ X ) 2p−2 , which holds true for any probability distribution f .We furthermore observe that, since the L p norm of every component is non-increasing by equation ( 3), assumption (6) translates into a condition on each ω i for positive times, Finally, we make a last hypothesis on the initial configuration, namely we suppose that there is not much vorticity far away from their centers.More precisely, we assume that there exists a radius R ≫ ε and a constant β such that where m i (t, R) is the vorticity portion of the i-th component at time t lying at least at distance R from its center X i (t), that is, We will assume that for some δ ∈ [0, 1/2) and some constant R 0 .
Our goal is to estimate to what extent the vortex components remain close to the Helmholtz-Kirchhoff point-vortex system starting from Ȳ1 , . . ., ȲN , that is a collection of points that move according to the equations We will assume that the actual point vortices are not colliding in some time interval [0, T c ], and that their minimal distance Our first result gives a componentwise estimate.
Theorem 1.There exist ε 0 , ν 0 ∈ (0, 1), β 0 ∈ (1, ∞) and a time T ≤ T c such that for any concentration parameter ε ∈ (0, ε 0 ), viscosity ν ∈ (0, ν 0 ) and decay parameter β ∈ (β 0 , ∞) the following holds: Let Ȳ1 , . . ., ȲN ∈ R 2 be pairwise distinct points and consider their evolution Y 1 , . . ., Y N ∈ R 2 through the point-vortex system (10).Let ω i be an evolving vortex region (3) that is initially sharply concentrated around the point vortices (4), decaying away from the centers (8) and satisfying the bound (6).Then it holds for all i = 1, . . ., N and t ∈ [0, T ] that where C and Λ are constants that are independent of ε, ν and R.Moreover, T ≥ c for some constant c independent of ε, ν and R. If T c = ∞ and δ > 0 in (9), there is the stronger estimate The estimate in Theorem 1 features the sharp rate of viscous spreading √ νt: In the context of the Navier-Stokes equations, the support of each initial component is instantaneously distributed all over Ê 2 .This effect can be nicely observed when studying the evolution of the self-similar Lamb-Oseen vortex which solves the Navier-Stokes equation with initial datum δ 0 (i.e. the Dirac measure concentrated on the origin).Apparently, the solution is positive for all positive times, and the vortex core grows in time at the rate √ νt.This growth law can be expressed in terms of the Wasserstein distance, as can be verified via a straightforward scaling argument.Our result thus proves that the general "vortex components" described initially in (4) show a similar (optimal) growth.Moreover, in the limit ε → 0, our result describes the optimal spreading rates for the Navier-Stokes equations with atomic initial data.Theorem 1 can also be stated as an estimate on the full solution (rather than the single vortex components) in terms of the 1-Wasserstein distance W 1 .Thanks to the dual Kantorovich-Rubinstein representation see for example Theorem 1.14 in [32], W 1 is indeed well-defined merely if f and g have the same average or total mass, with no requirement on their sign.
Corollary 1.Under the assumptions of Theorem 1, it holds that for all t ≤ T .
This implies weak convergence in the sense of measures of the Navier-Stokes vorticity ω(t) to the "point-vortex measure" N i=1 a i δ Y i (t) if we consider the simultaneous limits ε → 0, ν → 0, with a convergence rate that, up to an exponential in time, remains of the order of ε, as was assumed for the initial configuration in (4), plus a correction term √ νt due to the expansion of the vortices caused by the viscosity.
We remark that the point-vortex system (10) can be regarded as a very weak solution to the Euler equation, as was observed by Schochet [29].Hence, this result may also be interpreted as a stability estimate between the Navier-Stokes and Euler equations.Such a bound could be rigorously obtained by a simple application of the triangle inequality when combined with our earlier result [8], if suitable estimates on the convergence of viscous to inviscid solutions were available.These would be optimal presumably if an order of √ νt on the W 1 distance between the Navier-Stokes and Euler vorticities could be achieved, but, even in the Yudovich setting of integrable and bounded vorticities, this has not yet been established in general, see e.g.[11,12,10,30].It is interesting to note that, in the setting considered in this paper, one sees the expected rate √ νt (at least for viscosities of the order of ε 2 /t) in the inviscid limit ν → 0.
We can also express the estimate in terms of the vorticity centers of the components and their velocities as follows.
Theorem 2. Under the assumptions of Theorem 1, it holds that for all t ≤ T and i = 1, . . ., N.
We finally comment on our estimate on time T in Theorem 1.In view of the exponential growth term in our concentration estimates, the statement of our theorems become meaningless if T ≫ log 1 ε+ √ ν .In this regard, the bound on T is satisfactory at least if the "outer vorticity" in (8) vanishes appropriately outside a decreasing ball around the vortex center, i.e., δ > 0 in (9).The logarithmic bound ceases to hold if we fix a radius independently of ε and ν, and it is not clear to us, how to overcome this restriction.Comparable logarithmic bounds on time were obtained earlier by Cetrone and Serafini [9] under stronger concentration assumptions.The exponential growth term results from a Gronwall argument.It is very likely that better estimates on time would require a completely different (and genuinely nonlinear) approach.

Proofs
We remark that we treat ε and ν as two independent parameters and, since we are interested in what happens when both of them are not too large, in the computations we will always assume that ε < 1 and ν < 1.Moreover, ε will always be (sufficiently) smaller than the other length parameters R and d, that is ε ≪ R, d, and we notice that R is chosen in (9) Before starting with the proofs, we introduce some notation.In the following, C will always denote a positive constant independent of t, ν and ε, possibly depending on other quantities, and whose precise value may change from line to line.When we write A B we mean that A ≤ CB for some C, and when we write A ∼ B we mean that A B and B A. Hence for example we will often neglect the vortex intensities a i .We denote the velocity generated by the i-th vorticity as , where the convolution is meant in space.Moreover, we define the "external field" F i acting on ω i as the velocity field generated by all other vorticities, that is, F i = j =i u j .In this way, the total velocity u = K * ω satisfies The 2-Wasserstein distance of the i-th component from its center will be indicated by In our first step, we derive a bound on the growth of the Wasserstein distance.
3.1 Proof of the estimate on W i .
We fix i ∈ {1, . . ., N}.It will be convenient to regularize the Biot-Savart kernel by cutting out the singularity at the origin.We do this with the help of a radial cut-off We then split the external field acting on ω i into two parts, being the regular part, and F B i the remainder In the next lemma, we show that F L i is indeed a Lipschitz function while and F B i is bounded with for any t ≥ 0 such that where θ = (p − 2)/(3p − 2).
Proof.The Lipschitz bound on F L i can be easily proven computing the gradient of the integrand, , where we used the scaling of the Biot-Savart kernel in the form and the properties of the cut-off function (13).
We concentrate now on proving the boundedness of (14).Therefore, for any θ ∈ (0, 1), it holds that where we have used the (generealized) Hölder inequality with exponent q chosen so that 1/q + θ + (1 − θ)/p = 1.Requiring q < 2 ensures that the first integral in ( 15) is of the order d (2−q)/q .Solving for θ gives and thus, the condition that θ is positive enforces q ∈ (p/(p − 1), 2).Choosing q to be, say, the mean value between p/(p − 1) and 2 yields θ = θ as in the statement of the lemma.Second, if x / ∈ B d/6 (X i ), we use the rougher (in fact, globally valid) bound This concludes the proof.
In principle, we could estimate the L p norms of the vorticity components in the previous lemma with the help of the a priori estimate (7) and the assumption on the initial data (6).For large times, however, it is more convenient to make use of the smoothing properties of the diffusive part in (3).More precisely, we have the following estimate.Lemma 2. For any i = 1, . . ., N, q ∈ [1, ∞] and any t > 0, it holds that This estimate is known to be true for solutions to the Navier-Stokes equation (see, e.g., Theorem 4.3 in [18]), but due to possible cancellation effects, we cannot deduce its validity for individual vortex components.Instead, we provide a short proof of this estimate, in which we follow an argumentation given in [15].
Proof.It is enough to consider the statement for q = 2 k .The statement for general q < ∞ follows by interpolation, while the statement for q = ∞ comes from taking the limit k → ∞.The case q = 1 = 2 0 is trivially true since then ω i (t) L 1 = a i .We suppose that k ≥ 1 or, equivalently, q ≥ 2 from here on.
We consider E q (t) := ω i (t) q L q and observe that it satisfies the identity under the evolution (3).Notice that we use here the fact that q ≥ 2, so that the dependency on q in this estimate can indeed be neglected.Making use of the 2D Nash inequality f 2 , the latter turns into the estimate It remains to apply an induction argument over k, and we suppose thus that the statement holds for q/2.Then (17) implies and an integration in time gives from which we easily deduce the statement of the lemma.
In the next lemma we use the properties of F L i and F B i that we just showed to deduce a bound on the evolution of the Wasserstein distance between ω i and X i .Lemma 3. Let θ ∈ (0, 1) be given as in Lemma 1.For any t ≥ 0 with property (14), it holds that Proof.Since time is fixed, we will often forget to write it in this proof.We observe first that, thanks to the assumption that the initial vorticity components are compactly supported, we know that W i (0) < ∞.Because of the smoothness of ω i , the Wasserstein distance remains finite and differentiable for positive times.We consider the evolution of the squared Wasserstein distance, and compute, using the evolution (3) of the vorticity components, multiple integrations by parts and the definition of the vorticity centers X i , The first integral on the right-hand side vanishes because K is odd and because z • K(z) = 0. Using the definition of X i again and decomposing The first integral is easy to estimate by W 2 i /d 2 because of the Lipschitz property of F L i .On the second integral we use the different bounds obtained for F B i .Splitting the integration domain into B d/6 (X i ) and its complement, we have where we used Jensen's and Hölder's inequalities.Observe that in the estimates above the assumption that the vortex components have definite sign was necessary to bound It remains only to apply Lemma 1 and the proof is complete.Now we come to the objective of this subsection, that is a Gronwall estimate on the Wasserstein distance.Our argument will be the following.We start from the differential inequality (18), and for small times we observe that the L p norm of the vorticity inherits the bound (6) from the initial datum.For larger times, the effect of the viscosity kicks in, and we use the decay estimate (16).In both cases, in order to balance (6) or (16), we need to make sure that the portion of vorticity far from the centers is small enough.For this reason, we have to ensure not only that the vorticity centers remain sufficiently far from each other as in (14), but also that most of each vorticity component remains concentrated around its center.We observe that, under the localization condition (4) and the defintion of d in ( 12), we know that min provided that ε ≤ d/8, which we will suppose from here on.We select then T = T (ε, ν, R) as the minimal time for which at least one of the following equalities holds true: where α = α(p, γ) < β is a positive parameter to be fixed later, cf. ( 21).Here we observe that, thanks to ( 4), ( 19) and ( 8), and because both the vortex centers X i and the outer vorticity m i are continuous in time, the time T is always strictly positive, T > 0, .We thus have for all t ≤ min{T, T c } and any i.Moreover, using the triangle inequality and ( 12), we also notice that the latter entails for all t ≤ min{T, T c } and any i = j.We derive the estimate on the growth of W i for times smaller than T , while in Subsection 3.3 we will show that the time T can be chosen independently of ε and ν.
Proof.Using conditions (20), (12) and assumption (6) together with (7) in estimate (18), we obtain, in the case It is sufficient to assume that α is large enough, so that all the exponents of ε are larger or equal to 1.This condition is guaranteed by ( 21) and hence, being that ε < 1, R 1 and θ < p−2 p , we have for these times.There exists thus a constant C > 0 such that and an integration in time and using d ∼ 1 yields for any t ≤ ε 2 /ν, where we applied assumption (4) and inequality (5).
For larger times, √ νt ≥ ε, together with conditions ( 20) and ( 12), we use estimate (16) in the differential inequality (18), and get Using condition ( 21) on α and the lower bound γ ≥ 2 − 2/p ensures that all the exponents of νt are larger than 1/2.Then thanks to the facts that νt ≤ 1 and d ∼ 1, we obtain the inequality for these times.Similarly as above, we rewrite this differential inequality as with a (possibly larger) constant C > 0. Integration in time over t ≥ ε 2 /ν and using d ∼ 1 once more then yields where we have used the previous bound (22).Combining both, we obtain our thesis.
From now on, we choose α such that (21) holds.Moreover, in what follows, we will write Λ = C/d 2 in the exponenential growth factor.

Derivation of the main estimates
So far, we have derived a bound on the spreading rate of the single components in terms of the second moment function W i (t), that is, the Wasserstein distance between the vorticity components and their centers.However, we would like to express it in relation to the point-vortices solving (10).Hence, in our next step we attempt to estimate the distance |X i − Y i | between the vorticity centers and the point-vortices.
Proof of Theorem 2. The proof is a simple adaptation of an argument in the proof of Theorem 2.2 from [9].Up to considering ω j /a j instead of ω j , we assume without loss of generality that ω j is nonnegative and that a j = 1 for all j.Moreover, we often forget about the t's.We compute where with ψ the same cutoff function that has been defined in (13).Let us consider A 1 first.Recalling Proposition 1, , splitting the integration domain in B d/6 (X i ) and its complement and using (12), by Lemma 1 we have We proceed similarly to the proof of Proposition 1.In the case that √ νt < ε, we use ( 7), ( 6), (20), as well as condition (21) on α to obtain |A 2 | ε for these times.If √ νt ≥ ε, we use estimates ( 16), (20) and condition (21) to obtain |A 2 | √ νt also in this case.We consider now A 3 .Because of the cutoff function 1 − ψ, we may restrict the integral to |y −Y on the domain of integration.Hence we can bound is bounded by 1, and thus |A 4 | j =i B d/6 (Y i ) ω j (y) dy.We observe that on the domain of integration there holds again thanks to (20).Then |A 4 | j =i B d/6 (X j ) c ω j (y) dy = j =i m j (t, d/6), and again we use condition (20) to bound m j with ε or √ νt respectively in the case √ νt < ε and √ νt ≥ ε, keeping in mind that α > 1 by (21).In conclusion, all the estimates obtained yield for all t ≤ T .Summing over i = 1, . . ., N and using a Gronwall argument, together with the initial condition that | Xi − Ȳi | ≤ ε, yields the result.
Proof of Theorem 1.This is an easy consequence of Proposition 1 and Theorem 2 thanks to the triangle inequality.
Proof of Corollary 1. From the properties of the Wasserstein distance W 1 there follows The thesis is then a consequence of Jensen's inequality and Theorem 1.
In the following proofs we assume without loss of generality that T ≤ T c because otherwise there is nothing left to prove.We distinguish two cases.
In the first case, we suppose that for some i.Proof.This is any easy consequence of the first estimate in Theorem 2. Indeed, by the virtue of (23), it holds that which yields the statement of the lemma upon taking the logarithm and choosing ε + √ ν sufficiently small.
In the second case, we assume that for some i ∈ {1, . . ., N} We endeavour to prove that T is bounded uniformly from below in ε and ν also in this case.In order to do this, we derive an estimate on the outer vorticity portion m i (T, d/6) by employing an iteration estimate similar to the one used in [8].This in turn is based on the iterative procedure developed by Marchioro and co-workers, see for example [25,6] but, differently from these works, we start from the assumption (8) that the initial vortex components are quickly decaying outside of balls of radius R, while [25,6] make use to a great extent of their hypothesis that this radius is of the order of ε.
To estimate m i (T, d/6), we introduce the "smoothened outer vorticity portion" µ i (t, ρ, r) defined in the following.We consider a smooth radially symmetric cut-off function η = η ρ,r such that The i-th smoothened outer vorticity portion is defined as It is easy to check that the functions m i and µ i satisfy the relation We show a differential inequality for µ i by making use of the fact that, thanks to Proposition 1, holds for all t ≤ T ≤ 1.
Lemma 5. Let i ∈ {1, . . ., N} be fixed.For any t ≤ T , r > 0 and ρ ∈ (R, d/6 − r) it holds that where Proof.Up to dividing by a i , we can suppose without loss of generality that ω i ≥ 0 and a i = 1.Moreover, we often neglect the time t and we write µ i (t, ρ) = µ i (t, ρ, r) for notational convenience.
Because ω i solves the advection-diffusion equation (3), and thanks to the fact that the second identity being true by the fact that the self-interaction term vanishes because of the property K(z) = −K(−z), we can compute Using once more K(z) = −K(−z) we may write We notice that the integrand in non-zero only if either ρ < |x − X i | < ρ + r or ρ < |y − X i | < ρ + r.We can therefore split the integration domain into the sets and we denote by I A 1 , I B 1 and I C 1 the contributions to I 1 due to the sets A, B and C respectively.The terms I A 1 and I B 1 can be estimated in almost the same way.Using the radial symmetry of the cut-off function and orthogolality, we notice that and thus, since on their domain of integration it holds that |x − y| ≥ ρ/2, using the scaling of the gradient of η from (25) and Proposition 1 , we see that and the estimate of I B 1 proceeds analogously.To estimate I C 1 we use the Lipschitz condition (25) and estimate (27) to obtain Now we consider I 2 .We split We use the Lipschitz-continuity of F L i given by Lemma 1, property (25), the triangle and Jensen's inequalities and Lemma 1 to write Now we consider I 22 .By our choice of r and ρ, it holds that ρ + r ≤ d/6, and thus the integrand vanishes when x / ∈ B d/6 (X i ).We can therefore split the integration domain of I 22 into the following sets, and investigate the two contributions I D 22 and I E 22 separately.On the set D, we use (25) so that Thanks to Lemma 1 and property (20), it holds that and then, by the non-increasing property (7) of the L p norm, the scaling assumption (6), condition (21) on α and the bound (16) on ω i , we find also because νt < 1. Hence We estimate the contribution I E 22 due to the set E as Now, thanks to the first inequality in (20) and Lemma 1, we have Invoking property (7), assumption ( 6), the heat-kernel-type bound ( 16), and condition ( 21) on α, we may hence write From this and (30), we see that The term I 3 can be easily estimated using the bound (25), In conclusion, thanks to relation (26) and r ≤ ρ 1, Integrating in time, we obtain the thesis.
Iterating estimate (28), we come to the lower bound of T .
Proposition 2. Suppose that (24) holds and β ≥ 3α.Then there exist ε 0 , ν 0 ∈ (0, 1) and a constant c such that for any ε ≤ ε 0 ad ν ≤ ν 0 it holds that T ≥ c.If T c = ∞ and δ > 0 in (9), there is the stronger estimate Before continuing, we want to make sure that such an M is consistent with all the hypotheses we made so far, i.e., (32), (34) and (35).Let us start by treating the case ε ≥ √ ν.
First, regarding (34), we notice that Hence, in order to guarantee (34), it is enough to require that α log 2 log 1 + σ log 2 α The latter can be strengthend by estimating the logarithm linearly, If now (9) holds true with δ = 0, that is R 1, we necessarily have to choose ξ = 0 and this estimate holds true for σ sufficiently small.Otherwise, if δ > 0, we are allowed to choose ξ = 1 and we have to ensure that for some universal constant c ≥ 0. This is possible whenever σ ≤ δ and ε sufficiently small.We now turn to condition (35), which can be rewritten as which can be strengthend if √ ν on the left-hand side is dropped and ξ is estimated by 1.In this case, choosing σ small does the job.We finally turn to the worst case of (32), which can be rewritten as Estimating ε ≤ ε + √ ν ≤ 2ε and arguing similarly as above, we notice that this estimate can be strengthend to In view of our defintion of M, this estimate can be rewritten as and δ > 0, so that ξ = 1.Since the case δ = 0 is weaker then the case δ > 0, cf. ( 8) and ( 9), the bound on T in (44) carries over to the case δ > 0. Estimate (37) thus becomes ν 8C

Appendix: An elmentary estimate
In this appendix, we prove the elementary auxiliary estimate for any m ∈ AE and s > 0. In fact, thanks to the binomial theorem, it is enough to