arxiv: v1 [math.ap] 23 Jul 2017

Transcription

1 GLOBAL REGULARIY FOR 1D EULERIAN DYNAMICS WIH SINGULAR INERACION FORCES ALEXANDER KISELEV AND CHANGHUI AN arxiv: v1 [math.ap] 23 Jul 2017 Abstract. he Euler-Poisson-Alignment EPA) system appears in mathematical biology and is used to model, in a hydrodynamic limit, a set agents interacting through mutual attraction/repulsion as well as alignment forces. We consider one-dimensional EPA system with a class of singular alignment terms as well as natural attraction/repulsion terms. he singularity of the alignment kernel produces an interesting effect regularizing the solutions of the equation and leading to global regularity for wide range of initial data. his was recently observedin [5]. Our goal in this paper is to generalize the result and toincorporate the attractive/repulsive potential. We prove that global regularity persists for these more general models. Contents 1. Introduction and statement of main results Self-organized dynamics with three-zone interactions Global regularity versus finite time blowup Euler-Poisson-Alignment system Euler dynamics with general three-zone interactions 5 2. Euler-Poisson-Alignment system A priori bounds Local wellposedness Global wellposedness Euler dynamics with general three-zone interactions General singular influence function General interaction potential 19 References Introduction and statement of main results We consider the following 1D system of pressureless Euler equations with nonlocal interaction forces t ρ+ x ρu) = 0, x R, t > 0, 1) t u+u x u = ψx y)uy,t) ux,t))ρy,t)dy x K ρ, 2) R Date: July 25,

2 2 ALEXANDER KISELEV AND CHANGHUI AN subject to initial density and velocity ρ,t),u,t)) t=0 = ρ 0,u 0 ). he term on the right hand side of 2) consists of two parts: an alignment interaction with communication weight ψ, and an attraction-repulsion interaction through a potential K Self-organized dynamics with three-zone interactions. System 1)-2) arises from many contexts in mathematical physics and biology. In particular, it serves as a macroscopic system in modeling collective behaviors of complex biological systems. he corresponding agent-based model has the form ẋ i = v i, m v i = 1 N N ψx i x j )v j v i ) 1 N j=1 N xi Kx i x j ), 3) j=1 where x i,v i ) N i=1 represent the position and velocity of agent i. he dynamics is governed by Newton s second law, with the interaction force modeled under a celebrated threezone framework proposed in [17], including long-range attraction, short-range repulsion, and mid-range alignment. he first part of the force describes the alignment interaction, where ψ characterizes the strength of the velocity alignment between two agents. Naturally, it is a decreasing function of the distance between agents. Such alignment force has been proposed by Cucker and Smale in [4]. he corresponding dynamics enjoys the flocking property [9], which is a common phenomenon observed in animal groups. he second part of the force represents the attraction-repulsion interaction. he sign of the force K determines whether the interaction is attractive or repulsive. his type of potential driven interaction force is widely considered in many physical and biological models, e.g. [6, 15]. Starting from the agent-based model 3), one can derive a kinetic representation of the system that describes the mean-field behavior as N, see [2, 10, 21]. hen, a variety of hydrodynamics limits can be obtained that capture the macroscopic behaviors in different regimes [8, 11, 16]. In particular, if we consider the mono-kinetic regime, the corresponding macroscopic system becomes 1)-2) Global regularity versus finite time blowup. We are interested in the global existence and regularity for the solution of the system 1)-2). Let us start with the case with no interaction forces, namely ψ = K 0. he system can be recognized as the pressureless Euler system. In particular, 2) becomes the classical inviscid Burgers equation, where smooth data forms shock discontinuity in finite time due to nonlinear convection u x u. ogether with 1), it is well-known that the solution generates singular shocks in finite time: ρx,t) at the position and time when shock occurs. With alignment force ψ 0 and K 0, the system is called the Euler-Alignment system. When ψ is Lipschitz, the system has been studied in [1, 20], where it is discovered that the alignment force tends to regularize the solution and prevent finite time blowup, but only for some initial data. his is so called critical threshold phenomenon: for subcritical initial data, the alignment force beats the nonlinear convection, and the solution is globally regular; while for supercritical initial data, the convection wins and the solution admits a finite time blowup.

3 1D EULERIAN DYNAMICS WIH SINGULAR INERACION FORCES 3 Another interesting and natural setting is when ψ is singular, taking the form ψx) = c α x 1+α, α > 0, 4) with c α beapositive constant. herange 0 < α 2 is most natural, and the case 0 < α < 1 is most interesting for the reasons explained later in this sub-section. he Euler-Alignment system corresponding to the choice 4) is studied in [5] for the periodic case. Without loss of generality, weset thescale and let = [ 1/2,1/2] betheperiodicdomain of size 1. he singular alignment force can be equivalently expressed as ψ α y)ux+y,t) ux,t))ρx+y,t)dy, 5) with the periodic influence function ψ α defined as ψ α x) = c α x+m 1+α, x \0}. 6) m Z Clearly, ψ α is singular at x = 0. Moreover, it has a positive lower bound ) 1 ψ m = ψ m α) := min ψ αx) = ψ α > 0. 7) x 2 his leads to the following fractional Euler-Alignment system t ρ+ x ρu) = 0, x, 8) t u+u x u = ψ α y)ux+y,t) ux,t))ρx+y,t)dy. 9) It is shown in [5] that system 8)-9) has a global smooth solution for all smooth initial data with ρ 0 > 0. his result is most interesting for the case 0 < α < 1 : if we set ρ 1 in 9), we get Burgers equation with fractional dissipation. It is well-known that in this case, there exist initial data leading to finite time blow up when 0 < α < 1; the 1 α 2 range leads to global regularity). However, it turns out that in the nonlinear disspation/allignment case described by 9), the singularity in the influence function and density modulation dominate the nonlinaer convection, for all initial data. his also contrasts with the case of Lischitz regular influence function ψ, where one has critical threshold in the phase space separating initial data leading to finite time blow up and to global regularity Euler-Poisson-Alignment system. Now, we take into account the attraction-repulsion force, namely K 0. We shall begin with a particular potential Nx) = k x 2. hepotential isthe1d Newtonian potential, andit isthekernel forthe1dpoisson equation, namely 2 xn ρ = kρ. Whenk > 0, thenewtonian force 2 xn ρisattractive, andwhenk < 0, thenewtonianforce is repulsive. We call the corresponding system Euler-Poisson-Alignment EPA) system. It has the form t ρ+ x ρu) = 0, 10) t u+u x u = x φ+ ψx y)uy,t) ux,t))ρy,t)dy, 11) R

4 4 ALEXANDER KISELEV AND CHANGHUI AN where the stream function φ = N ρ satisfies the Poisson equation 2 xφ = kρ. 12) When there is no alignment force ψ 0, the system coincides with the 1D pressureless Euler-Poisson equation, which has been extensively studied in [7]. he result is as follows: when k > 0, the attraction force together with convection drives the solution of Euler- Poisson equation to a finite time blowup for all smooth initial data; when k < 0, the repulsive force competes with the convection, and there exists a critical threshold on initial conditions which separates global regularity and finite time blowup. he EPA system 10)-12) is studied in [1], in the case when ψ is Lipschitz. When k < 0, a larger subcritical region of initial data is obtained that ensures global regularity. his implies that the alignment force helps repulsive potential to compete with the convection. However, it is also shown that when k > 0, the alignment force is too weak to compete with convection and attractive potential, so all smooth initial data lead to finite time blow up. Our first result concerns EPA system with singular alignment force, where the influence function has the form 4). he main goal is to understand whether the singular alignment can still regularize the solution when the Newtonian force is present. We shall study the system in the periodic setting. he 1D periodic Newtonian potential reads Nx) = k ) x, x. 13) It is the kernel of the Poisson equation with background, namely where ρ is the average density ρ = 1 x 2 N ρ) = kρ ρ), ρx,t)dx = ρ 0 x)dx. Note that ρ is conserved in time due to conservation of mass by evolution. he stream function φ in 11) satisfies the Poisson equation with constant background 2 xφ = kρ ρ). 14) he presence of the background ρ could change the behavior of the solution. For Euler- Poisson equation in periodic domain, namely 10)-11),14) with ψ 0, it is pointed out in [7] that the background has the tendency to balance both the convection and attractive forces. So for the attractive case k > 0, instead of finite time blowup for all initial data, a critical threshold is obtained. hough similar techniques in [1], one can derive critical thresholds for EPA system 10)- 11),14) with bounded Lipschitz influence function ψ. he EPA system with singular alignment force 5) and potential 13) reads t ρ+ x ρu) = 0, t u+u x u = x φ+ 15) ψ α y)ux+y,t) ux,t))ρx+y,t)dy, 2 xφ = kρ ρ).16) he following theorem shows that the singular alignment force dominates the Poisson force, and global regularity is obtained for all initial data. heorem 1.1. For α 0, 1), the fractional EPA system 15)-16) with smooth periodic initial data ρ 0,u 0 ) such that ρ 0 > 0 has a unique smooth solution.

5 1D EULERIAN DYNAMICS WIH SINGULAR INERACION FORCES 5 Remark 1.1. he proof can be easily extended to the range α 1 with more straightforward arguments for α > 1; see also [18] for a different approach. We focus on the 0 < α < 1 case in the rest of the paper. We note that the proof of global regularity in [5] is based, in particular, on rather precise algebraic structures that we will discuss below. Even though the interaction force we are adding is formally sub-critical, it is far from obvious that the fairly intricate arguments of [5] survive such perturbation Euler dynamics with general three-zone interactions. he results on EPA system can be extended to systems with more general interaction forces. In [1], critical thresholds are obtained for the system 1)-2), with Lipschitz influence function ψ, and regular potential K W 2,. In this paper, we will also consider the case of more general singular influence function ψ. More precisely, we assume that ψ ψ m > 0, and can be decomposed into two parts ψ = cψ α +ψ L, 17) where c > 0, ψ α is defined in 6), and ψ L is a bounded Lipschitz function. heorem 1.2. Consider system 1)-2) in the periodic setup t ρ+ x ρu) = 0, x, t > 0, t u+u x u = ψy)ux+y,t) ux,t))ρx+y,t)dy x K ρ, with smooth initial data ρ 0,u 0 ) such that ρ 0 > 0. Assume ψ is singular in the sense of 17), and K is a linear combination of Newtonian potential 13) and regular W 2, ) potential. hen, the system has a unique global smooth solution. We summarize the global behaviors of Euler equations with nonlocal interaction 1)-2) under different choices of interaction forces. Potential Alignment Name Domain Behaviors No No Euler R or Finite time blow up Lipshitz Euler-Alignment R or Critical threshold [1, 20] Singular Fractional EA Global regularity [5] Newtonian No Euler-Poisson R Finite time blow up [7] Critical threshold [7] Lipschitz EPA R Finite time blow up attractive) Critical threshold repulsive) [1] Critical threshold Singular Fractional EPA Global regularity heorem 1.1) General Lipshitz Euler-3Zone R or Critical thresholds [1] Singular Singular 3Zone Global regularity heorem 1.2)

6 6 ALEXANDER KISELEV AND CHANGHUI AN 2. Euler-Poisson-Alignment system In this section, we consider Euler-Poisson-Alignment system15)-16) with singular alignment force 5). Following the idea in [5], we let G = x u Λ α ρ, and calculate the dynamics of G using 15) and 16): t G = t x u t Λ α ρ = x u x u) kρ ρ)+ x Λ α ρu)+uλ α ρ)+λ α x ρu) = u x x u Λ α ρ) x u x u Λ α ρ) kρ ρ) = x Gu) kρ ρ). So, we can rewrite the dynamics in terms of ρ,g) as he velocity u can be recovered as t ρ+ x ρu) = 0, 18) t G+ x Gu) = kρ ρ), 19) x u = Λ α ρ+g. ux,t) = Λ α 1 x ρx,t) ρ)+ 1 x Gx,t)+I 0 t), 20) where I 0 can be determined by conservation of momentum. ρx,t)ux,t)dx = ρ 0 x)u 0 x)dx. 21) See [5] for detailed discussion A priori bounds. We first show an upper and lower bounds on density ρ for all finite times. For k = 0, a uniform in time bound is obtained in [5]. With the Newtonian potential, especially when k > 0, the attractive force definitely helps density concentration. Hence, the upper bound on ρ can be expected to grow in time. However, the bound we obtain in this section indicates that there is no finite time singular concentration on density, thanks to the singular alignment force. Let F = G/ρ. We can rewrite 18) as he first step is to obtain a bound on F. We calculate t F = ρ tg G t ρ ρ 2 his implies that t +u x )ρ = ρλ α ρ ρ 2 F. 22) = ρ xgu) kρ ρ)) G x ρu)) ρ 2 = u x F kρ ρ). ρ t +u x )F = k 1 ρ ). 23) ρ Denote Xx, t) the trajectory of the characterstic path starting at x, namely d Xx,t) = uxx,t),t), Xx,0) = x. 24) dt hen, we can solve for F along the characteristic path FXx,t),t) = F 0 x) kt+ t 0 k ρ ds. 25) ρxx,s),s)

7 1D EULERIAN DYNAMICS WIH SINGULAR INERACION FORCES 7 Define ρ m t) as the lower bound of ρ on time interval [0,t] hen, we get a bound on F from 25): ρ m t) = min s [0,] min x ρx,s). F,t) L F 0 L + k t+ k ρ t 0 1 ds. 26) ρ m s) herefore, in order to control F in L, we need a lower bound estimate on the density. heorem 2.1 Lower bound on density). Let ρ,u) be a strong solution to EPA system 15)16) with smooth periodic initial conditions ρ 0,u 0 ) such that ρ m 0) > 0. hen, there exist two positive constants A m and C m, depending only on the initial conditions, such that for any t 0, ρ m t) C m e Amt. 27) Proof. We depart from 22) and estimate Λ α ρ and F. For a fixed time t, denote x be a point where ρ attains its minimum. Note that x depends on t and it is not necessarily unique. he estimates below apply at any such point. We have Λ α ρx+y,t) ρx,t) ρx,t) = c α y 1+α dy = ψ α y) ρx+y,t) ρx,t) ) dy 28) ) ) ψ m ρx+y,t) ρx,t) dy = ψm ρ ρx,t). Here, we recall that ψ m is the positive lower bound of ψ α defined in 7). Combining 26) and 28), we obtain t ρx,t) ψ ) [ t ] 1 m ρ ρx,t) ψ m + F 0 L + k t+ k ρ ρ m s) ds ρx,t) 2. 29) We prove 27) by contradiction. For t = 0, the bound 27) holds if we let C m ρ m 0). Suppose 27) does not hold for all t 0. hen, there exists a finite time t 0 > 0 so that the inequality is violated for the first time at t = t 0 +. Pick any x = xt 0 ). Due to continuity of ρ, we know ρ m t 0 ) = ρx,t 0 ) = C m e Amt 0. 30) Plug in 30) to 29) and use the fact that 27) holds for all t [0,t 0 ]. We get t ρx,t 0 ) ρ m t 0 )[ ψm ρ ) t0 ) ] 1 ψ m + F 0 L + k t+ k ρ 0 ρ m s) ds ρ m t 0 ) ρ m t 0 )[ ψm ρ ) ψ m + F 0 L + k t 0 + k ρ ) ] e Amt 0 1) C m e Amt 0 A m C m [ ρ m t 0 ) ψ k ρ ) m ρ ψ m + F 0 L + k t 0 k ρ ) ] C m e Amt 0 A m A m C m [ ρ m t 0 ) ψ k ρ m ρ k C ) ) ] m k ρ + C m ψ m + F 0 L ) e Amt 0. A m ea m A m he right hand side is positive if we pick A m large enough and C m small enough. For instance, we can pick A m = k ψ m 1+ǫ), C m = minρ m 0),ǫe ρ}, 31) 0

8 8 ALEXANDER KISELEV AND CHANGHUI AN 1+ 4ψ m eψ m+ F 0 L ) 1 ). With this choice of A m and for any ǫ 0,ǫ ), where ǫ = 1 2 C m, we get t ρx,t 0 ) > 0. Now we obtain that ρx) < C m e Amt 0 < C m e Amt for some t < t 0. his contradicts our choice of t 0. Remark 2.1. he bound 27) with decay rate 31) is not necessarily sharp, but is enough for our purpose, as it eliminates the possibility of finite time creation of vacuum. One important observation is that for k = 0, we get A m = 0. In this case, the lower bound is uniform in time. Applying the lower bound 27) to 26), we immediately derive a bound on F F,t) L F 0 L + k t+ k ρ A m C m e Amt =: F M t). 32) Now, we are ready to obtain an upper bound on density ρ. heorem 2.2 Upper bound on density). Let ρ,u) be a strong solution to EPA system 15)16) with smooth periodic initial conditions ρ 0,u 0 ) such that ρ m 0) > 0. hen, there exist two positive constants A M and C M, depending only on the initial conditions, such that for any t 0 and x, ρx,t) ρ M t) := C M e A Mt. 33) Proof. We again departfrom22) andstart with alower boundestimate on Λ α ρ. For a fixed time t, denote x be a point where ρ attains its maximum. Applying nonlinear maximum principle by Constantin and Vicol [3], one can estimate Λ α ρ x,t) C 1 ρ x,t) 1+α, 34) if ρ x,t) 3 ρ. he constant C 1 only depends on initial conditions. One can consult [5] for more details of the estimate. Plugging the estimates 32) and 34) into 22), we obtain t ρ x,t) C 1 ρ x,t) 2+α +F M t)ρ x,t) 2. 35) It follows that t ρ x,t) < 0 if ρ x,t) > F M /C 1 ) 1/α. herefore, ) } FM t) 1/α ρx,t) ρ x,t) max ρ 0 L,3 ρ,, and 33) holds with A M = A m α, C M = max max x ρ 0x), 3 ρ, [ 1 C 1 C 1 F 0 L + k + k ρ )] } 1/α. ea m A m C m 2.2. Local wellposedness. With the apriori bounds, we state a local wellposedness result for the fractional EPA system 15)-16), as well as a Beale-Kato-Majda type necessary and sufficient condition to guarantee global wellposedness. he local wellposedness theory has been presented in detail in [5] for fractional Euler-Alignment system. We will show that presence of the Poisson force does not seriously affect the argument, no matter whether it is attractive or repulsive. We will only sketch the proof, indicating changes necessary.

9 1D EULERIAN DYNAMICS WIH SINGULAR INERACION FORCES 9 heorem 2.3 Local wellposedness). Consider EPA system 15)-16) with initial conditions ρ 0 and u 0 that satisfy ρ 0 H s ), min x ρ 0x) > 0, x u 0 Λ α ρ 0 H s α 2 ), with a sufficiently large even integer s > 0. hen, there exists 0 > 0 such that the EPA system has a unique strong solution ρx,t),ux,t) on [0, 0 ], with ρ C[0, 0 ],H s )) L 2 [0, 0 ],H s+α 2)), u C[0, 0 ],H s+1 α )). 36) Moreover, a necessary and sufficient condition for the solution to exist on a time interval [0,] is 0 x ρ,t) 2 L dt <. 37) Proof. We follow the proof in [5] and rewrite the equations 18) and 19) in terms of θ,g) where θ = ρ ρ. t θ+ x θu) = ρ x u, t G+ x Gu) = kθ, he velocity u is defined in 20). Given any > 0, we will obtain a differential inequality on Yt) := 1+ θ,t) 2 H s + G,t) 2 H s α 2, 38) for all t [0,]. hrough a commutator estimate [5, equation 3.23)], one can get 1 d 2dt θ 2 H s C 1+ 1 ) 1+ x θ 2 L ρ + G L )Yt) ρ m m 3 θ 2, 39) H s+α 2 where ρ m t) has a positive lower bound for t [0,] due to heorem 2.1. Also, G,t) L is bounded for t [0,] as G = Fρ and both F and ρ are bounded, see 32) and 33) respectively. We also compute 1 d 2dt G 2 = Λ s α 2G) Λ s α 2 Ḣ s α 2 x Gu))dx k Λ s α 2G) Λ s α 2θ)dx = I +II. he first term can be controlled by the following estimate [5, equation 3.25)] I ρ m 6 θ 2 +C 1+ 1 ) G 2 H s+α 2 L ρ + xθ L + G L ) G 2. 40) m H s α 2 he II term encodes the contribution of the attractive-repulsive potential. We have the following estimate II k G Ḣs α θ 2 Ḣ C k G s α 2 Ḣ θ s α 2 Hs C k Yt). 41) Combine 39), 40) and 41), we get d dt Yt) C1+ xθ,t) 2 L )Yt) ρ mt) θ 2 6, 42) H s+α 2 where C is a positive constant which might depend on.

10 10 ALEXANDER KISELEV AND CHANGHUI AN Applying Gronwall s inequality, we get Yt)+ 1 [ 6 min ρ mt) θ 2 Y0)exp C) t [0,] L 2 [0,];H s+α 2 )) 0 1+ x θ,s) 2 L )ds ], for all t [0,]. he right hand side is bounded as long as condition 37) is satisfied. herefore, θ C[0,],H s )) L 2 [0,],H s+α 2 )), G C[0,],H s α 2 )). his directly implies the regularity conditions on ρ in 36). he regularity conditions on u can also be easily obtained from 20) Global wellposedness. In this section, we prove that the Beale-Kato-Majda type condition 37) holds for any finite time. his will imply global wellposedness of the fractional EPA system and hence finish the proof of heorem 1.1. hroughout the section, we fix a time > 0 which is arbitrary). o derive a uniform L bound on x ρ, we argue that ρ,t) will obey certain modulus of continuity for t [0,]. Such method has been successfully used to obtain global regularity for 2D quasi-geostrophic equation with critical dissipation [14], fractal Burgers equation [13], as well as fractional Euler-Alignment system [5]. In all these examples, the solution has a certain scaling invariance property. Unfortunately, such property is not available for the fractional EPA system 15)-16). We note that the modulus method has been applied to subcritical perturbations destroying scaling before e.g. [12]). he argument in [12], however, relies on the specific structure of the perturbation, and cannot be readily ported to other settings. A novel feature compared to both [12] and [5] will be dependence of the modulus on time. his feature is linked to the possible decay of ρ m and growth of ρ L, and appears to be an intrinsic property of the problem. We use the same family of moduli of continuity as in [5], ξ ξ 1+α/2, 0 ξ < δ 1 ωξ) = γlogξ/δ)+δ δ 1+α/2 43), ξ δ, where γ,δ are small constants to be determined. Set ω B ξ) = ωbξ), where B is a large constant to be determined as well. Due to lack of scaling invariance, we will work directly on ω B. Bξ Bξ) 1+α/2, 0 ξ < B 1 δ ω B ξ) = γlog Bξ δ +δ δ1+α/2, ξ B 1 44) δ, We say that a function f obeys modulus of continuity ω if fx) fy) < ω x y ), x,y. Our plan is to find a ω B such that ρ,t) obeys ω B for all t [0,]. o construct ω B, we will first choose δ and γ which depend on initial conditions and, but not on B. hen, we will choose B that depend on,δ,γ as well as initial conditions. First, we would like to make sure that ρ 0 obeys ω B. Lemma 2.1. Let ρ 0 C 1 ). hen, ρ 0 obeys ω B if δ < 2 ρ 0 L, B > δ xρ 0 L exp x ρ 0 L 2 ρ 0 L 2 ρ0 L γ ). 45)

11 1D EULERIAN DYNAMICS WIH SINGULAR INERACION FORCES 11 Proof. We start with an elementary inequality ρ 0 x) ρ 0 y) min2 ρ 0 L, x ρ 0 L x y }. 46) As ω B is concave and monotone increasing, the right hand side of 46) is bounded by ω B x y ) if ) 2 ρ0 L ω B > 2 ρ 0 L. 47) x ρ 0 L Since ω B ξ) + as B +, 47) is satisfied by taking B large enough. Indeed, if δ and B satisfy 45), then ) ) 2 ρ0 L 2B ρ0 L ω B > γlog > 2 ρ 0 L. x ρ 0 L δ x ρ 0 L he following lemma describes the only possible breakthrough scenario for the modulus. Lemma 2.2. Suppose ρ 0 obeys a modulus of continuity ω B as in 44). If the solution ρx,t) violates ω B at some positive time, then there must exist t 1 > 0 and x 1 y 1 such that ρx 1,t 1 ) ρy 1,t 1 ) = ω B x 1 y 1 ), and ρ,t) obeys ω B for every 0 t < t 1. 48) he main point of the lemma is the existence of two distinct points where the solution touches the modulus as opposed to a single point x with ρx) = ω B 0) = B). his property is a consequence of ω B 0) = ; see [14] for more details. We will show that in the breakthrough scenario as above, t ρx 1,t 1 ) ρy 1,t 1 )) < 0, t 1 0,], 49) achieving a contradiction with the choice of time t 1 - and thus showing that the modulus ω B cannot be broken. ogether with Lemma 2.1 this implies that ρ,t) obeys ω B for all t [0,]. herefore, x ρ,t) L ω B0) = B, t [0,]. his proves global regularity of the fractional EPA system and ends the proof of heorem 1.1. he rest of the section is devoted to proof of 49). We fix t 1 and drop the time variable for simplicity. Let ξ = x 1 y 1. hen t ρx 1 ) ρy 1 )) = x ρx 1 )ux 1 ))+ x ρy 1 )uy 1 )) = ux 1 ) x ρx 1 ) uy 1 ) x ρy 1 ) ) ρx 1 ) ρy 1 ) ) x ux 1 ) ρy 1 ) x ux 1 ) x uy 1 ) ) = I +II +III. Decompose u into two parts u = u 1 +u 2 where hen, we can write 50) as u 1 x) = Λ α 1 x ρx) ρ), u 2x) = 1 x Gx)+I 0. t ρx 1 ) ρy 1 )) = I 1 +II 1 +III 1 +I 2 +II 2 +III 2, where I 1,II 1,III 1 represent the contributions from u 1, and I 2,II 2,III 2 represent the contribution from u 2. 50)

12 12 ALEXANDER KISELEV AND CHANGHUI AN For I 1,II 1,III 1, we proceed with an argument parallel to [5]. Let us recall the result. he following quantities play a role in the proof: ξ ωη) ) ωη) Ωξ) = c 1,α dη +ξ dη ; 0 ηα ξ η1+α ωξ) ω ξ η ) Aξ) = c 2,α dη; Dξ) = c 3,α ξ/2 0 R η 1+α 2ωξ) ωξ +2η) ωξ 2η) η 1+α dη + ξ/2 2ωξ) ωξ +2η)+ω2η ξ) η 1+α Lemma 2.3 [5, Lemma 4.4 and 4.5]). Let ρ,t) obey the modulus of continuity ω as in 43) for 0 t < t 1, and let x 1, y 1 be the breakthrough points at the first breakthrough time t 1. Suppose δ and γ are small constants such that δ < 1, γ δ δ1+α/2. 51) 2log2 hen, there exist positive constants C I, C II and C III, which may only depend on α, such that I 1 ω C I ξ, 0 < ξ < δ, ξ)ωξ), where Ωξ) C I ξ 1 α 52) ωξ), ξ δ. C II, 0 < ξ < δ, II 1 ωξ)aξ), where Aξ) C II γξ α 53), ξ δ. C III ξ 1 α/2, 0 < ξ < δ, III 1 ρ m Dξ), where Dξ) C III ωξ)ξ α 54), ξ δ. Applying the proof of Lemma 2.3 to the modulus of continuity ω B, we get the following estimates. Lemma 2.4. Let ρ,t) obey the modulus of continuity ω B as in 44) for 0 t < t 1, and let x 1, y 1 be the breakthrough points at the first breakthrough time t 1, as in 48). Suppose δ and γ are small constants satisfying 51). hen there exist positive constants C 2 and C 3, which may only depend on α, such that C 2 B 1+α ξ, 0 < ξ < B 1 δ, I 1, II 1 C 2 γω B ξ)ξ α, ξ B 1 55) δ, and III 1 ρ m D B ξ), D B ξ) := C 3 B 1+α/2 ξ 1 α/2, 0 < ξ < B 1 δ, C 3 ω B ξ)ξ α, ξ B 1 δ. Proof. hrough the same proof of Lemma 2.3 and replacing ω by ω B, one can obtain the following estimates similar to 52), 53) and 54): I 1 ω Bξ)Ω B ξ), II 1 ω B ξ)a B ξ), III 1 ρ m D B ξ). Here ω B is defined in 44), and Ω B ξ) = B α 1 ΩBξ), A B ξ) = B α ABξ), D B ξ) = B α DBξ). his directly implies 55) and 56) with C 2 = maxc I,C II } and C 3 = C III. dη ). 56)

13 then 1D EULERIAN DYNAMICS WIH SINGULAR INERACION FORCES 13 If we pick δ small enough so that δ < ) 2/α C3 ρ m ), 57) 4C 2 C 2 B 1+α ξ C 2 B 1+α/2 ξ 1 α/2 δ α/2 1 4 ρ md B ξ), ξ 0,B 1 δ). Also, pick γ small enough so that then herefore, we have γ < C 3 4C 2 ρ m ), 58) C 2 γω B ξ)ξ α 1 4 ρ md B ξ), ξ B 1 δ. I 1 +II 1 +III ρ md B ξ). 59) It remains to control I 2,II 2 and III 3. We start with the estimate on I 2. Lemma 2.5. Let ρ,t) obey the modulus of continuity ω B as in 44) for 0 t < t 1, and let x 1, y 1 be the breakthrough points at the first breakthrough time t 1, as in 48). Suppose δ and γ satisfy 51), and in addition δ < hen, ρm )C 3 6ρ M )F M ) ) 2/α,γ < αδ δ 1+α/2 ),and B > max 1,2δexp I ρ md B ξ). 6ρM )F M ) C 3 ρ m ) )}. 60) Proof. We start with the estimate I 2 x u 2 L ξω B ξ) = G L ξω B ξ) ρ M)F M )ξω B ξ). For ξ 0,B 1 δ), ω B ξ) < B. So, I 2 ρ M )F M )Bξ 1 6 ρ md B ξ), 61) provided that δ is small enough, satisfying 60), and B > 1. For ξ B 1 δ, since ρ is periodic and ω B is increasing, the breakthrough can not happen first at ξ > 1/2. So we only need to consider ξ B 1 δ,1/2]. As ω B ξ) = γ ξ in this range, we get I 2 ρ M )F M )γ. 62) On the other hand, compute d dξ D Bξ) = C 3 ξ α 1 αω B ξ)+γ) C 3 ξ α 1 αδ δ 1+α/2 )+γ) < 0, for all ξ B 1 δ, provided that γ is small enough, satisfying 60). herefore, ) B min D Bξ) = D B 1/2) C 3 γlog. 63) B 1 δ ξ 1/2 2δ

14 14 ALEXANDER KISELEV AND CHANGHUI AN Combining 62), 63) and the assumption on B in 60), we conclude I 2 ρ M )F M )γ C ) 3 B 6 γρ m)log 1 2δ 6 ρ md B ξ). 64) heestimates on II 2 and III 2 are more subtle. o proceed, it is convinient to decompose II 2 +III 2 in an alternative way II 2 +III 2 = ρx 1 ) x u 2 x 1 ) ρy 1 ) x u 2 y 1 ) ) = ρx 1 ) 2 Fx 1 ) ρy 1 ) 2 Fy 1 ) ) = ρx 1 ) 2 ρy 1 ) 2) Fx 1 ) ρy 1 ) 2 Fx 1 ) Fy 1 ) ) = IV +V. We first consider the case when ξ < B 1 δ. For IV, the estimate is similar to 61) IV = ω B ξ)ρx 1 )+ρy 1 )) Fx 1 ) 2ρ M F M Bξ 1 6 ρ md B ξ), 65) where the last inequality holds if δ is picked to be small enough, satisfying ) C 3 ρ m ) 2/α δ <. 66) 12ρ M )F M ) For V, we need the following lemma. Lemma 2.6. Let ρ,t) obey the modulus of continuity ω B with any B > 1 as in 44) for 0 t < t 1. hen, there exists a constant C F = C F ) such that Fx,t) Fy,t) C F )B x y, x,y, t [0,t 1 ]. 67) Proof. Recall the dynamics of F t F +u x F = k Let f = x F. Differentiate 68) with respect to x and get Let q = f/ρ. Using 15) and 69), we obtain 1 ρ ). 68) ρ t f + x uf) = k ρ xρ ρ 2. 69) t q +u x q = k ρ xρ ρ 3. It follows that t x ρxx,s),s) qxx,t),t) = q 0 x) k ρ 0 ρxx,s),s) 3 ds, where X is the trajectory of the characteristic path defined in 24). hen, since for t t 1, ρ,t) obeys ω B, we obtain the following estimate q,t) L q 0 L + k ρ t 0 B ρ m s) 3ds C )B, where the finite constant C depends on and initial data. his implies Fx) Fy) f L x y ρ M )C )Bξ =: C F )B x y.

15 1D EULERIAN DYNAMICS WIH SINGULAR INERACION FORCES 15 Applying the estimate 67) at the breakthrough points and using the upper bound on ρ 33), we get V ρ M ) 2 C F )Bξ < 1 6 ρ md B ξ), 70) where the second inequality holds by picking sufficiently small δ, satisifying ) C 3 ρ m ) 2/α δ < 6ρ M ) 2, 71) C F ) similar to the estimate in 61). Combining 59), 61), 65) and 70), we conclude that t ρx 1 ) ρx 2 )) < 0, ξ = x 1 x 2 < B 1 δ. Finally, we estimate II 2 +III 2 for ξ [B 1 δ,1/2]. As ρ and F are bounded, it is clear that II 2 +III 2 2ρ M ) 2 F M ) < 1 3 ρ md B ξ). 72) he second inequality holds by picking B large enough. his is due to the fact that D B ξ) is an increasing in B with lim B D B ξ) =. More precisely, using the bound 63), it suffices to pick 6ρM ) 2 ) F M ) B > 2δexp. 73) C 3 γρ m ) Combining 59), 64) and 72), we conclude that t ρx 1 ) ρx 2 )) < 0, ξ = x 1 x 2 [B 1 δ,1/2]. Let us summerize the procedure on the construction of the modulus of continuity ω B. First, we fix a time. hen, we pick a small parameter δ satisfying 45), 51), 57), 60), 66) and 71): δ < min 1, 2 ρ 0 L, x ρ 0 L C 3 ρ m ) max4c 2,12ρ M )F M ),6ρ M ) 2 C F )} ) 2/α }. 74) Next, we pick a small parameter γ satisfying 51), 58) and 60): ) } C3 1 γ < min ρ m ), min 4C 2 2log2,α δ δ 1+α/2 ). 75) Finally, we pick a large parameter B satisfying 45), 60) and 73): B > max 1, δ ) xρ 0 L 2 ρ0 L 6ρM ) 2 )} F M ) exp,2δexp. 76) 2 ρ 0 L γ C 3 γρ m ) Here we assume, without loss of generality, that γ 1 and ρ M ) 1 to simplify the expression. With this choice of ω B, we have shown that ρ,t) obeys ω B for all t [0,]. Hence, x ρ,t) L B, t [0,]. herefore, condition 37) is satisfied, and we obtain global regularity of the system. We end this section by the following remark.

16 16 ALEXANDER KISELEV AND CHANGHUI AN Remark 2.2. When k = 0, all the quantities ρ m,ρ M,F M and C F do not depend on. As a consequence, δ,γ and B do not depend on as well. herefore, x ρ,t) L B for any t 0. his estimate improves the result obtained in [5], where the bound on x ρ could grow in time. We note that stationary in time bound on x ρ for the Euler-Alignment model without Possion forcing) has been derived in [19] by a different argument. For k 0, with the singular attractive or repulsive force, our estimate on ρ m and ρ M is not uniform in time. We are able to obtain time-dependent bounds 27) and 33), where ρ m can decay exponentially in time, and ρ M and F M,C F ) can grow exponentially in time. From 74) and 75), we see that δ and γ decay exponentially in time. Finally, from 76), B grows double exponentially in time. herefore, we obtain a double exponential in time bound on x ρ,t) L. It is not clear whether such bound is optimal. We will leave it for future investigation. 3. Euler dynamics with general three-zone interactions In this section, we extend our global regularity result for EPA system to more general Euler dynamics with three-zone interactions. Recall the Euler-3Zone system under periodic setup t ρ+ x ρu) = 0, x, t > 0, 77) t u+u x u = ψy)ux+y,t) ux,t))ρy,t)dy x K ρ. 78) We will discuss the global wellposedness of the system with more general singular influence function ψ and interaction potential K General singular influence function. Consider a general influence function ψ which is positive ψ m := minψx) > 0, x and singular at origin. Recall the decomposition 17): we will consider the class of functions where one can decompose ψ into two parts ψ = cψ α +ψ L. Here ψ α is the singular power defined in 6), and ψ L is bounded and Lipschitz. In this case, let G = x u cλ α ρ+ψ L ρ. hen, the dynamics of G reads: t G = t x u c t Λ α ρ+ψ t ρ = x u x u)+c x Λ α ρu)+uλ α ρ ) x ψl ρu) uψ L ρ) ) 2 xxk ρ +cλ α x ρu) ψ L x ρu) = x u x u cλ α ρ+ψ L ρ) ) 2 xxk ρ = x Gu) 2 xxk ρ. herefore, ρ, G) still satisfy 18) and 19), with a different relation t ρ+ x ρu) = 0, t G+ x Gu) = 2 xxk ρ, 79) x u = Λ α ρ+g ψ L ρ.

17 1D EULERIAN DYNAMICS WIH SINGULAR INERACION FORCES 17 hen the velocity field u can be recovered as ux,t) = Λ α x 1 ρx,t) ρ)+ x 1 Gx,t) ψl ρx,t) ) +I 0 t), 80) where I 0 t) can be determined by conservation of momentum 21). he second term on the right hand side is well-defined since Gx,t) ψl ρx,t) ) dx = x ux,t) Λ α ρx,t) ) dx = 0, t 0. We can decompose u into two parts u = u S +u L, where u S is the singular part ) u S x,t) =Λ α x 1 ρx,t) ρ)+ x 1 Gx,t) Gx, 0)dx, x u S =Λ α ρ+g Gx, 0)dx, and u L is the Lipschitz part u L x,t) = x 1 ψ L ρx,t) x u L = ψ L ρ+ Gx, 0)dx. ) Gx, 0)dx +I 0 t), Now, we follow the same procedure as fractional EPA system to show global regularity of system79), 80). We first take the Newtonian potential 13). General interaction potentials will be discussed in the next section. he arguments below follow the same outline, so we focus on indicating changes. Step 1: Apriori lower bound on ρ. he statement and proof are identical to heorem 2.1, except that estimate 28) is replaced by cλ α ρx,t)+ψ L ρx,t) = cψα y)+ψ L y) ) ρx y,t) ρx,t) ) dy +ρx,t) ) ψ m ρx y,t) ρx,t) dy ρx,t) ψl L = ψ m ρ ψ m + ψ L L )ρx,t). Hence, estimate 29) becomes t ρx,t) ψ ) [ m ρ ρx,t) ψ m + ψ L L + F 0 L + k t+ k ρ t 0 ψ L y)dy 81) 82) ] 1 ρ m s) ds ρx,t) 2, where the only extra term ψ L L ρx,t) 2 is quadratic in ρ, and can be controlled by the linear term ψ ) m ρ ρx,t) if ρ m is small enough. Following the same proof, we obtain the lower bound 27) with coefficient A m,c m satisfying 31) for any ǫ 0,ǫ ), where ǫ = ψ m eψ m+ ψ L L + F 0 L ) 1 ). Step 2: Apriori upper bound on ρ. We follow the proof of heorem 2.2. he estimate 35) becomes d dt ρ x,t) C 1ρ x,t) 2+α +F M t)ρ x,t) 2 +ρ x,t) ψ L ρ x,t) C 1 ρ x,t) 2+α +F M t)ρ x,t) 2 + ψ L L ρρ x,t).

18 18 ALEXANDER KISELEV AND CHANGHUI AN Both second and third terms are dominated by the first term if ρ x,t) is big enough. In particular t ρ x,t) < 0 if ρ x,t) > max2f M /C 1 ) 1/α,2 ψ L L ρ/c 1 ) 1/1+α) }. herefore, ) 2FM t) 1/α ) } 2 ψl L ρ 1/1+α) ρx,t) ρ x,t) max ρ 0 L,3 ρ,,, 83) and 33) holds with A M = A m /α and [ 2 C M = max max ρ 0x), 3 ρ, F 0 L + k + k ρ )] 1 x C 1 ea m A m C m C 1 C 1 α, 2 ψl L ρ C 1 ) 1 } 1+α. Step 3: Local wellposedness. We write the system 79) 80) in terms of θ = ρ ρ and G as follows t θ + x θu S )+ x θu L ) = ρ x u S ρ x u L, 84) t G+ x Gu S )+ x Gu L ) = kθ, 85) where u S and u L are defined in 81) and 82) respectively. We proceed with a Gronwall estimate on the quantity Y in 38). he estimates in heorem 2.3 can be applied directly to the u S part. We will focus on the Lipschitz part u L. he procedure is similar to [1, heorem Appendix A.1]. We will summarize in below. For the term x θu L ), we have Λ s θ Λ s x θu L )dx = Λ s θ Λ s x θ u L dx+ Λ s θ [Λ s x,u L ]θ dx =: L 1 +L 2. We estimate the two terms one by one. For L 1, L 1 = Λ s θ) 2 ) x u L dx 2 1 Λ s θ) 2 x u L dx ψ L L ρ θ 2 Hs. 86) For L 2, we apply commutator estimate e.g. [1, Lemma Appendix A.1]) and get L 2 θ H s [Λ s x,u L ]θ L 2 θ H s x u L L θ H s + x u L H s θ L ) θ H s ψ L L ρ θ H s + ψ L L ρ H s θ L ) ψ L L 2 ρ+ θ L ) θ 2 H s ψ L L ρ. 87) Note that for the last inequality, we have used ρ H s θ H s + ρ H s = θ H s + ρ. For the term ρ x u L, ρ Λ s θ Λ s x u L dx ρ θ H s xψ L ) Λ s ρ) L 2 ρ x ψ L L θ H s θ H s+ ρ). 88) For the term x Gu L ), the estimate is similar to the term x θu L ). Λ s α 2G Λ s α 2 x Gu L )dx = Λ s α 2 G Λ s α 2 x G u L dx+ Λ s α 2 G [Λ s α 2 x,u L ]G dx =: L 4 +L 5. where L 4 = ) Λ s α 2G) 2 x u L dx 2 1 Λ s α 2G) 2 x u L dx ψ L L ρ G 2, 89) H s α 2

19 and L 5 G H s α 2 1D EULERIAN DYNAMICS WIH SINGULAR INERACION FORCES 19 [Λ s α 2 x,u L ]G L 2 G H s α 2 ) x u L L G H s α + 2 xu L H s G L ψ L L ρ G 2 + ψ L H s α 2 L ρ θ H s + ρ) G H s α 2 ] 12 ρ2 ψ L L ρ [2 G 2Hs α2 + θ 2Hs α ) 2 Combining 42), 86), 87), 88), 89), 90) and the fact that G,t) L is controlled from above by a finite growing in time) bound, we obtain that for all t [0,], d dt Yt) C)1+ xθ,t) 2 L )Yt) ρ mt) θ 2 6, H s+α 2 where the constant C depends on initial data and. he same Gronwall s inequality yields local wellposedness as well as BKM-type blowup condition 37). Step 4: Global wellposedness. o check the condition 37), we will use the procedure identical to that in section 2.3. Let us decompose u as in 80), u = u 1 +u 2 where u 1 x,t) = Λ α x 1 ρx,t) ρ), u 2 x,t) = x 1 Gx,t) ψl ρx,t) ) +I 0 t). 91) he only difference between our system 79) 80) and the EPA system is that there is an extra term in u 2. hroughout the proof in section 2.3, the only property of u 2 we have used is that x u 2 is bounded, namely x u 2,t) L ρ M )F M ) <, t [0,]. For our u 2 defined in 91), we also have a bound on x u 2 : x u 2,t) L = G,t) ψ L ρ,t) L ρ M )F M )+ ψ L L ρ <, t [0,]. Hence, global regularity follows from the same procedure by controlling the modulus of continuity General interaction potential. In this part, we consider system 77)-78) with a general interaction potential K W 2, ). his class of potentials is more regular than the Newtonian potential N defined in 13), as 2 xxn = kδ 0 1) L, where δ 0 is the Dirac delta at x = 0. We will show global wellposedness of Euler-3Zone system with W 2, potentials. he result automatically extends to systems with potentials that can be decomposed into a sum of a Newtonian potential and a W 2, potential. Now, let us assume K W 2, ). After the transformation, the dynamics for ρ,g) becomes 79), withvelocity field udefinedas 80). We shall runthroughthesameprocedure and point out the differences. Step 1: Apriori lower bound on ρ. Due to the change of the potential, the dynamics of F 23) becomes t +u x )F = 2 xxk ρ. 92) ρ herefore, we get t FXx,t),t) = F 0 x) 0 xx 2 K ρxx,s),s) ds, ρxx, s), s)

20 20 ALEXANDER KISELEV AND CHANGHUI AN where Xx,t) is the characteristic path defined in 24). hen, we obtain a bound F,t) L F 0 L + xx K L ρ t 0 1 ds, 93) ρ m s) which is similar as 26). In fact, it is a simpler boundas the right handside does not contain a linear term on t. he lower bound 27) follows then by the same argument, with A m = 2 xx K L, C m = min ψ m ρ m 0), ψ m ρ ψ m + ψ L L + F 0 L Step 2: Apriori upper bound on ρ. he upper bound estimate 83) can be obtained without any additional difficulties. Since we have F M t) = F 0 L + 2 xx K L ρ e Amt A m C m by 93) and the lower bound estimate on ρ, the upper bound 33) holds with A M = A m /α and C M = max max x ρ 0x), 3 ρ, [ 2 C 1 F 0 L + 2 xx K L ρ )] 1 A m C m }. α, 2 ψl L ρ C 1 ) 1 } 1+α. Step 3: Local wellposedness. Since the potential only enters the dynamics of G equation, so the system in terms of θ,g) is identical to 84)-85), except the right hand side of 85) is replaced by xxk 2 ρ. Hence, we only need to estimate this extra term. Λ s α 2G Λ s α 2 xx 2 K ρ) dx = Λ s α 2G xx 2 K 2ρ) Λs α dx 2 xx K L G H s α 2 θ H s xx K L Yt). he local wellposedness and BKM-type blowup condition 37) follow by applying the same Gronwall s inequality on Y. Step 4: Global wellposedness. he argument for EPA system can be directly applied to the general system as the ρ equations in both cases are the same. he different potential does change the estimates on ρ m,ρ M,F M,C F, which are needed to construct the modulus ω B. Since ρ m,ρ M and F M have been treated in the previous steps, we are left with estimating C F, namely proving Lemma 2.6 for the general system. Proof of Lemma 2.6. Let f = x F. Differentiate 92) with respect to x and get t f + x uf) = 3 xxxk ρ)ρ+ 2 xxk ρ) x ρ ρ 2. 94) Let q = f/ρ. Using 15) and 94), we obtain t q +u x q = 3 xxxk ρ)ρ+ 2 xxk ρ) x ρ ρ 3. 95) For t t 1, ρ,t) obeys ω B. hen x ρ,t) L ω B 0) = B. herefore, we can bound the right hand side of 95) as follows: 3 xxxk ρ)ρ+ 2 xxk ρ) x ρ ρ 3 2 xxk L x ρ L 1 ρ m t) xxk L ρ x ρ L ρ m t) 3

21 1D EULERIAN DYNAMICS WIH SINGULAR INERACION FORCES 21 B 2 xx K L 1 ρ m t) ρ ρ m t) 3 hen, we obtain the bound on q for all 0 t t 1 <, t ) 1 q,t) L q 0 L +B xxk 2 L ρ m t) 2 + ρ ρ m t) 3 ds C )B, where the finite constant C depends on 2 xxk L, and initial data. his implies Fx) Fy) f L x y ρ M )C )Bξ =: C F )B x y. Acknowledgment. his work has been partially supported by the NSF grants DMS and DMS References [1] José A Carrillo, Young-Pil Choi, Eitan admor, and Changhui an. Critical thresholds in 1d euler equations with nonlocal forces. Mathematical Models and Methods in Applied Sciences, 261): , [2] José A Carrillo, Massimo Fornasier, Jesús Rosado, and Giuseppe oscani. Asymptotic flocking dynamics for the kinetic Cucker-Smale model. SIAM Journal on Mathematical Analysis, 421): , [3] Peter Constantin and Vlad Vicol. Nonlinear maximum principles for dissipative linear nonlocal operators and applications. Geometric And Functional Analysis, 225): , [4] Felipe Cucker and Steve Smale. Emergent behavior in flocks. Automatic Control, IEEE ransactions on, 525): , [5] am Do, Alexander Kiselev, Lenya Ryzhik, and Changhui an. Global regularity for the fractional Euler alignment system. arxiv preprint arxiv: , [6] Maria R D Orsogna, Yao-Li Chuang, Andrea L Bertozzi, and Lincoln S Chayes. Self-propelled particles with soft-core interactions: patterns, stability, and collapse. Physical review letters, 9610):104302, [7] Shlomo Engelberg, Hailiang Liu, and Eitan admor. Critical thresholds in Euler-Poisson equations. Indiana University Mathematics Journal, 50: , [8] Razvan C Fetecau, Weiran Sun, and Changhui an. First-order aggregation models with alignment. Physica D: Nonlinear Phenomena, 325: , [9] Seung-Yeal Ha and Jian-Guo Liu. A simple proof of the Cucker-Smale flocking dynamics and mean-field limit. Communications in Mathematical Sciences, 72): , [10] Seung-Yeal Ha and Eitan admor. From particle to kinetic and hydrodynamic descriptions of flocking. Kinetic and Related Models, 13): , [11] rygve K Karper, Antoine Mellet, and Konstantina rivisa. Hydrodynamic limit of the kinetic Cucker- Smale flocking model. Mathematical Models and Methods in Applied Sciences, 2501): , [12] Alexander Kiselev and Fedor Nazarov. Global regularity for the critical dispersive dissipative surface quasi-geostrophic equation. Nonlinearity, 233):549, [13] Alexander Kiselev, Fedor Nazarov, and Roman Shterenberg. Blow up and regularity for fractal burgers equation. Dynamics of PDE, 53): , [14] Alexander Kiselev, Fedor Nazarov, and Alexander Volberg. Global well-posedness for the critical 2d dissipative quasi-geostrophic equation. Inventiones mathematicae, 1673): , [15] Alexander Mogilner and Leah Edelstein-Keshet. A non-local model for a swarm. Journal of Mathematical Biology, 386): , [16] David Poyato and Juan Soler. Euler-type equations and commutators in singular and hyperbolic limits of kinetic Cucker Smale models. Mathematical Models and Methods in Applied Sciences, 276): , [17] Craig W Reynolds. Flocks, herds and schools: A distributed behavioral model. ACM Siggraph Computer Graphics, 214):25 34, [18] Roman Shvydkoy and Eitan admor. Eulerian dynamics with a commutator forcing. ransactions of Mathematics and Its Applications, 11), [19] Roman Shvydkoy and Eitan admor. Eulerian dynamics with a commutator forcing III. fractional diffusion of order 0 < α < 1. arxiv preprint arxiv: , ).

22 22 ALEXANDER KISELEV AND CHANGHUI AN [20] Eitan admor and Changhui an. Critical thresholds in flocking hydrodynamics with non-local alignment. Philosophical ransactions of the Royal Society of London A: Mathematical, Physical and Engineering Sciences, ): , [21] Changhui an. A discontinuous Galerkin method on kinetic flocking models. Mathematical Models and Methods in Applied Sciences, 277): , Alexander Kiselev) Department of Mathematics, Duke University, 120 Science Drive, Durham, NC 27708, USA address: kiselev@math.duke.edu Changhui an) Department of Mathematics, Rice University, 6100 Main St., Houston, X 77005, USA address: ctan@rice.edu